Limited MHC polymorphism in the southern elephant seal: implications for MHC evolution and marine mammal population biology

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Novel insights on aquatic mammal MHC evolution: Evidence from manatee DQB diversity

The low diversity in marine mammal major histocompatibility complex (MHC) appears to support the hypothesis of reduced pathogen selective pressure in aquatic systems compared to terrestrial environments. However, the lack of characterization of the aquatic and evolutionarily distant Sirenia precludes drawing more generalized conclusions. Therefore, we aimed to characterize the MHC DQB diversity of two manatee species and compare it with those reported for marine mammals. Our results identified 12 and 6 alleles in T. inunguis and T. manatus, respectively. Alleles show high rates of nonsynonymous substitutions, suggesting loci are evolving under positive selection. Among aquatic mammals, Pinnipeda DQB had smaller numbers of alleles, higher synonymous substitution rate, and a dN/dS ratio closer to 1, suggesting it may be evolving under more relaxed selection compared to fully aquatic mammals. This contradicts one of the predictions of the hypothesis that aquatic environments impose reduced pathogen pressure to mammalian immune system. These results suggest that the unique evolutionary trajectories of mammalian MHC may impose challenges in drawing ecoevolutionary conclusions from comparisons across distant vertebrate lineages.

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.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

Divergent Evolution of TRC Genes in Mammalian Niche Adaptation

Mammals inhabit a wide variety of ecological niches, which in turn can be affected by various ecological factors, especially in relation to immunity. The canonical TRC repertoire (TRAC, TRBC, TRGC, and TRDC) codes C regions of T cell receptor chains that form the primary antigen receptors involved in the activation of cellular immunity. At present, little is known about the correlation between the evolution of mammalian TRC genes and ecological factors. In this study, four types canonical of TRC genes were identified from 37 mammalian species. Phylogenetic comparative methods (phyANOVA and PGLS) and selective pressure analyses among different groups of ecological factors (habitat, diet, and sociality) were carried out. The results showed that habitat was the major ecological factor shaping mammalian TRC repertoires. Specifically, trade-off between TRGC numbers and positive selection of TRAC and the balanced evolutionary rates between TRAC and TRDC genes were speculated as two main mechanisms in adaption to habitat and sociality. Overall, our study suggested divergent mechanisms for the evolution of TRCs, prompting mammalian immunity adaptions within diverse niches.

Genomic analysis of MHC-based mate choice in the monogamous California mouse

Variation at Major Histocompatibility Complex (MHC) genes is thought to be an important mechanism underlying mate choice in vertebrates, with individuals typically predicted to prefer MHC-dissimilar reproductive partners. However, analyses based on individual MHC loci have generated contradictory results regarding the role of these genes in mate-choice decisions. To provide a more comprehensive assessment of relationships between MHC variation and mating behavior, we used an exome capture strategy to characterize variability at 13 MHC loci, 312 innate immune system genes, and 1044 nonimmune genes in 25 obligate monogamous pairs of California mice (Peromyscus californicus) from 2 free-living populations of this species in Monterey County, California. Pairwise genotypic comparisons and analyses of SNP-based allelic differences failed to detect disassortative mating based on MHC variability; reproductive partners were not more dissimilar than randomly generated male-female pairs at MHC, innate or nonimmune loci. Within populations, individuals tended to be more closely related at MHC genes than at innate or nonimmune genes. Consistent with the functional role of immunogenes, the 2 study populations were highly differentiated at MHC and innate genes but not at nonimmune loci. Collectively, our results suggest that MHC genetic variation in California mice reflects local differences in pathogen exposure rather than disassortative mating based on variability at MHC Class I and II genes.

Genomic organization and adaptive evolution of IGHC genes in marine mammals

•

The genomic organization of IGHC genes in marine mammal is similar to that of terrestrial relatives.

•

The number of the immunoglobulin heavy chain constant region genes vary among different mammals.

•

Different levels of selective pressures were detected between marine and terrestrial mammalian lineages.

Immunoglobulins are important elements of the adaptive immune system that bind to an immense variety of microbial antigens to neutralize infectivity and specify effector functions. In the present study, the immunoglobulin heavy chain constant region (IGHC) genes from marine mammals were identified and compared with those of their terrestrial relatives to explore their genomic organization and evolutionary characteristics. The genomic organization of marine mammal IGHC genes was shown to be conservative with other eutherian mammals. Stronger signals of positive selection on IGHC were revealed in terrestrial mammals than that in marine mammals with the branch-site model, displaying different selective pressure, which might suggest their divergent adaptations to contrasted environments.

Effects of Population Bottleneck and Balancing Selection on the Chinese Alligator Are Revealed by Locus-Specific Characterization of MHC Genes

Chinese alligator (Alligator sinensis) is an endangered freshwater crocodilian endemic to China, which experienced a severe bottleneck about 30 years ago. In this study, we developed locus-specific primers to investigate the polymorphism of 3 major histocompatibility complex (MHC) loci in 3 Chinese alligator populations, in combination with 6 neutral microsatellite markers as a contrast. We found the genetic trace for the bottleneck effect on the endangered Chinese alligator: the low allelic diversity (2 alleles at each locus), the low nucleotide substitution rate (no more than 0.009) at all sites, the deviation from Hardy-Weinberg Equilibrium/heterozygote deficiency, and the significant Tajima’s D values, indicating the MHC class I and class II loci being at different stages of bottleneck. We also obtained 3 pieces of evidence for balancing selection on this severely bottlenecked reptile: an obvious excess of nonsynonymous substitutions over synonymous at the antigen-binding positions, the mean synonymous substitution rate of MHC exons significantly higher than mean nucleotide substitution rate of introns, and the differentiation coefficient F_ST of MHC loci significantly lower than that of microsatellite loci. Consequently, we emphasize that the Chinese alligator holds a pretty low adaptive ability and requires scientific conservation strategies to ensure the long-term population development.

Major histocompatibility complex variation and the evolution of resistance to amphibian chytridiomycosis

Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has been implicated in population declines and species extinctions of amphibians around the world. Susceptibility to the disease varies both within and among species, most likely attributable to heritable immunogenetic variation. Analyses of transcriptional expression in hosts following their infection by Bd reveal complex responses. Species resistant to Bd generally show evidence of stronger innate and adaptive immune system responses. Major histocompatibility complex (MHC) class I and class II genes of some susceptible species are up-regulated following host infection by Bd, but resistant species show no comparable changes in transcriptional expression. Bd-resistant species share similar pocket conformations within the MHC-II antigen-binding groove. Among susceptible species, survivors of epizootics bear alleles encoding these conformations. Individuals with homozygous resistance alleles appear to benefit by enhanced resistance, especially in environmental conditions that promote pathogen virulence. Subjects that are repeatedly infected and subsequently cleared of Bd can develop an acquired immune response to the pathogen. Strong directional selection for MHC alleles that encode resistance to Bd may deplete genetic variation necessary to respond to other pathogens. Resistance to chytridiomycosis incurs life-history costs that require further study.

Organization and characteristics of the major histocompatibility complex class II region in the Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis)

Little is known about the major histocompatibility complex (MHC) in the genome of Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) (YFP) or other cetaceans. In this study, a high-quality YFP bacterial artificial chromosome (BAC) library was constructed. We then determined the organization and characterization of YFP MHC class II region by screening the BAC library, followed by sequencing and assembly of positive BAC clones. The YFP MHC class II region consists of two segregated contigs (218,725 bp and 328,435 bp respectively) that include only eight expressed MHC class II genes, three pseudo MHC genes and twelve non-MHC genes. The YFP has fewer MHC class II genes than ruminants, showing locus reduction in DRB, DQA, DQB, and loss of DY. In addition, phylogenic and evolutionary analyses indicated that the DRB, DQA and DQB genes might have undergone birth-and-death evolution, whereas the DQB gene might have evolved under positive selection in cetaceans. These findings provide an essential foundation for future work, such as estimating MHC genetic variation in the YFP or other cetaceans. This work is the first report on the MHC class II region in cetaceans and offers valuable information for understanding the evolution of MHC genome in cetaceans.

Assembly and characterization of the MHC class I region of the Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis)

The Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis; YFP) is the sole freshwater subspecies of N. asiaeorientalis and is now critically endangered. Major histocompatibility complex (MHC) is a family of highly polymorphic genes that play an important immunological role in antigen presentation in the vertebrates. Currently, however, little is known about MHC region in the genome of the YFP, which hampers conservation genetics and evolutionary ecology study using MHC genes. In this work, a nucleotide sequence of 774,811 bp covering the YFP MHC class I region was obtained by screening a YFP bacterial artificial chromosome (BAC) library, followed by sequencing and assembly of positive BAC clones. A total of 45 genes were successfully annotated, of which four were MHC class I genes. There are high similarities among the four YFP MHC class I genes (>94%). Divergence in the coding region of the four YFP MHC class I genes is mainly localized to exons 2 and 3, which encode the antigen-binding sites of MHC class I genes. Additionally, comparison of the MHC structure in YFP to those of cattle, sheep, and pig showed that MHC class I genes are located in genome regions with regard to the conserved genes, and the YFP contains the fewest MHC class I genes among these species. This is the first report characterizing a cetacean MHC class I region and describing its organization, which would be valuable for further investigation of adaptation in natural populations of the YFP and other cetaceans.

Morphological and genomic comparisons of Hawaiian and Japanese Black-footed Albatrosses (Phoebastria nigripes) using double digest RADseq: implications for conservation

Evaluating the genetic and demographic independence of populations of threatened species is important for determining appropriate conservation measures, but different technologies can yield different conclusions. Despite multiple studies, the taxonomic status and extent of gene flow between the main breeding populations of Black-footed Albatross (Phoebastria nigripes), a Near-Threatened philopatric seabird, are still controversial. Here, we employ double digest RADseq to quantify the extent of genomewide divergence and gene flow in this species. Our genomewide data set of 9760 loci containing 3455 single nucleotide polymorphisms yielded estimates of genetic diversity and gene flow that were generally robust across seven different filtering and sampling protocols and suggest a low level of genomic variation (θ per site = ∼0.00002-0.00028), with estimates of effective population size (N e = ∼500-15 881) falling far below current census size. Genetic differentiation was small but detectable between Japan and Hawaii (F ST ≈ 0.038-0.049), with no F ST outliers. Additionally, using museum specimens, we found that effect sizes of morphological differences by sex or population rarely exceeded 4%. These patterns suggest that the Hawaiian and Japanese populations exhibit small but significant differences and should be considered separate management units, although the evolutionary and adaptive consequences of this differentiation remain to be identified.

Stress physiology in marine mammals: how well do they fit the terrestrial model?

Stressors are commonly accepted as the causal factors, either internal or external, that evoke physiological responses to mediate the impact of the stressor. The majority of research on the physiological stress response, and costs incurred to an animal, has focused on terrestrial species. This review presents current knowledge on the physiology of the stress response in a lesser studied group of mammals, the marine mammals. Marine mammals are an artificial or pseudo grouping from a taxonomical perspective, as this group represents several distinct and diverse orders of mammals. However, they all are fully or semi-aquatic animals and have experienced selective pressures that have shaped their physiology in a manner that differs from terrestrial relatives. What these differences are and how they relate to the stress response is an efflorescent topic of study. The identification of the many facets of the stress response is critical to marine mammal management and conservation efforts. Anthropogenic stressors in marine ecosystems, including ocean noise, pollution, and fisheries interactions, are increasing and the dramatic responses of some marine mammals to these stressors have elevated concerns over the impact of human-related activities on a diverse group of animals that are difficult to monitor. This review covers the physiology of the stress response in marine mammals and places it in context of what is known from research on terrestrial mammals, particularly with respect to mediator activity that diverges from generalized terrestrial models. Challenges in conducting research on stress physiology in marine mammals are discussed and ways to overcome these challenges in the future are suggested.

The impact of translocations on neutral and functional genetic diversity within and among populations of the Seychelles warbler

Translocations are an increasingly common tool in conservation. The maintenance of genetic diversity through translocation is critical for both the short- and long-term persistence of populations and species. However, the relative spatio-temporal impacts of translocations on neutral and functional genetic diversity, and how this affects genetic structure among the conserved populations overall, have received little investigation. We compared the impact of translocating different numbers of founders on both microsatellite and major histocompatibility complex (MHC) class I diversity over a 23-year period in the Seychelles warbler (Acrocephalus sechellensis). We found low and stable microsatellite and MHC diversity in the source population and evidence for only a limited loss of either type of diversity in the four new populations. However, we found evidence of significant, but low to moderate, genetic differentiation between populations, with those populations established with fewer founders clustering separately. Stochastic genetic capture (as opposed to subsequent drift) was the main determinant of translocated population diversity. Furthermore, a strong correlation between microsatellite and MHC differentiation suggested that neutral processes outweighed selection in shaping MHC diversity in the new populations. These data provide important insights into how to optimize the use of translocation as a conservation tool.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

De novo assembly and characterization of spotted seal Phoca largha transcriptome using Illumina paired-end sequencing

Spotted seal (Phoca largha) is categorized as a critically endangered species in China. The aim of this study was to investigate spotted seal transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from the spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In the transcriptome, 193 unigenes were assigned to defense mechanisms. Three unigenes encoded MHC class I and 17 unigenes encoded MHC class II. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified the expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results would contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Evolution of the MHC-DQB exon 2 in marine and terrestrial mammals

On the basis of a general low polymorphism, several studies suggest that balancing selection in the class II major histocompatibility complex (MHC) is weaker in marine mammals as compared with terrestrial mammals. We investigated such differential selection among Cetacea, Artiodactyla, and Primates at exon 2 of MHC-DQB gene by contrasting indicators of molecular evolution such as occurrence of transpecific polymorphisms, patterns of phylogenetic branch lengths by codon position, rates of nonsynonymous and synonymous substitutions as well as accumulation of variable sites on the sampling of alleles. These indicators were compared between the DQB and the mitochondrial cytochrome b gene (cytb) as a reference of neutral expectations and differences between molecular clocks resulting from life history and historical demography. All indicators showed that the influence of balancing selection on the DQB is more variable and overall weaker for cetaceans. In our sampling, ziphiids, the sperm whale, monodontids and the finless porpoise formed a group with lower DQB polymorphism, while mysticetes exhibited a higher DQB variation similar to that of terrestrial mammals as well as higher occurrence of transpecific polymorphisms. Different dolphins appeared in the two groups. Larger variation of selection on the cetacean DQB could be related to greater stochasticity in their historical demography and thus, to a greater complexity of the general ecology and disease processes of these animals.

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.

MHC genotype and near-deterministic mortality in grey seals

The Major Histocompatability Complex (MHC) is one of the best known and best characterised components of the immune system, yet its functions remain somewhat enigmatic, including both anti-pathogen activity and kin recognition. To explore the importance of the MHC relative to literally hundreds of other components of the immune system, we compared MHC genotype frequencies between pups and adults in the grey seal (Halichoerus grypus), one of many marine mammals that exhibit low allelic diversity. We find that one allele is strongly associated with pup survival, pups being more likely to be found dead if they lack it, while total allele number is a remarkably strong predictor of survivorship to adulthood. We estimate that approximately 70% of mortality can be attributed to the MHC. Our study therefore shows that low MHC allele diversity belies its critical role in determining whether a weaned pup negotiates disease to become a breeding adult.

Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers

Next-generation sequencing provides a powerful new approach for developing functional genomic tools for non-model species. This study aims to investigate the spotted seal (Phoca largha) transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

The functional gene diversity in natural populations over postglacial areas: the shaping mechanisms behind genetic composition of longnose dace (Rhinichthys cataractae) in northeastern North America

The diversity of functional genes and the related processes are important issues for conservation biology. This is especially relevant for populations that have suffered from demographic reduction as a consequence of the processes of postglacial colonization. In this perspective, the aims of the present study are (1) to quantify the genetic diversity of functional genes and (2) to disentangle the long- and short-term effects of natural selection that shapes genetic diversity from those of drift, mutation, and allopatric fragmentation. This research was conducted using an extensive genetic polymorphism analysis of populations of longnose dace (Rhinichthys cataractae) living over an area once covered by Pleistocene glaciations. The sequence and diversity of one exon of three genes (MHC IIβ, growth hormone, and trypsin) were jointly analyzed with non-coding nuclear loci from 27 populations; these populations were sampled over four major basins of northeastern North America. The survey revealed a surprisingly low allelic richness, especially for the MHC gene, considering the number of individuals and populations sampled. The results suggest that there is a complex mixture of different evolutionary processes shaping the level of polymorphism among longnose dace. While our study underlines the importance of the short-term effects of neutral processes and the major impact of post-glacial colonization on gene diversity, locally dependent balancing selection was detected on MHC. From this perspective, our results support an understanding of the importance of drift on functional gene diversity but also highlight the transient effects of natural selection on allelic composition, even in populations that show drastic reduction of genetic diversity.

Reduced MHC and neutral variation in the Galápagos hawk, an island endemic

BACKGROUND

Genes at the major histocompatibility complex (MHC) are known for high levels of polymorphism maintained by balancing selection. In small or bottlenecked populations, however, genetic drift may be strong enough to overwhelm the effect of balancing selection, resulting in reduced MHC variability. In this study we investigated MHC evolution in two recently diverged bird species: the endemic Galápagos hawk (Buteo galapagoensis), which occurs in small, isolated island populations, and its widespread mainland relative, the Swainson's hawk (B. swainsoni).

RESULTS

We amplified at least two MHC class II B gene copies in each species. We recovered only three different sequences from 32 Galápagos hawks, while we amplified 20 unique sequences in 20 Swainson's hawks. Most of the sequences clustered into two groups in a phylogenetic network, with one group likely representing pseudogenes or nonclassical loci. Neutral genetic diversity at 17 microsatellite loci was also reduced in the Galápagos hawk compared to the Swainson's hawk.

CONCLUSIONS

The corresponding loss in neutral diversity suggests that the reduced variability present at Galápagos hawk MHC class II B genes compared to the Swainson's hawk is primarily due to a founder event followed by ongoing genetic drift in small populations. However, purifying selection could also explain the low number of MHC alleles present. This lack of variation at genes involved in the adaptive immune response could be cause for concern should novel diseases reach the archipelago.

Geographic variation of the major histocompatibility complex in Eastern Atlantic grey seals (Halichoerus grypus)

Pathogen-driven balancing selection maintains high genetic diversity in many vertebrates, particularly in the major histocompatibility complex (MHC) immune system gene family, which is often associated with disease susceptibility. In large natural populations where subpopulations face different pathogen pressures, the MHC should show greater genetic differentiation within a species than neutral markers. We examined genetic diversity at the MHC-DQB locus and nine putatively neutral microsatellite markers in grey seals (Halichoerus grypus) from eight United Kingdom (UK) colonies, the Faeroe Islands and Sable Island, Canada. Five DQB alleles were identified in grey seals, which varied in prevalence across the grey seal range. Among the seal colonies, significant differences in DQB allele and haplotype frequencies and in average DQB heterozygosity were observed. Additionally, the DQB gene exhibited greater differentiation among colonies compared with neutral markers, yet a weaker pattern of isolation by distance (IBD). After correcting for the underlying IBD pattern, subpopulations breeding in similar habitats were more similar to one another in DQB allele frequencies than populations breeding in different habitats, but the same did not hold true for microsatellites, suggesting that habitat-specific pathogen pressure influences MHC evolution. Overall, the data are consistent with selection at MHC-DQB loci in grey seals with both varying selective pressures and geographic population structure appearing to influence the DQB genetic composition of breeding colonies.

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.

MHC gene copy number variation in Tasmanian devils: implications for the spread of a contagious cancer

Tasmanian devils face extinction owing to the emergence of a contagious cancer. Devil facial tumour disease (DFTD) is a clonal cancer spread owing to a lack of major histocompatibility complex (MHC) barriers in Tasmanian devil populations. We present a comprehensive screen of MHC diversity in devils and identify 25 MHC types and 53 novel sequences, but conclude that overall levels of MHC diversity at the sequence level are low. The majority of MHC Class I variation can be explained by allelic copy number variation with two to seven sequence variants identified per individual. MHC sequences are divided into two distinct groups based on sequence similarity. DFTD cells and most devils have sequences from both groups. Twenty per cent of individuals have a restricted MHC repertoire and contain only group I or only group II sequences. Counterintuitively, we postulate that the immune system of individuals with a restricted MHC repertoire may recognize foreign MHC antigens on the surface of the DFTD cell. The implication of these results for management of DFTD and this endangered species are discussed.

Sequence polymorphism and geographical variation at a positively selected MHC-DRB gene in the finless porpoise (Neophocaena phocaenoides): implication for recent differentiation of the Yangtze finless porpoise?

Sequence polymorphism at the MHC class II DRB locus was investigated in three finless porpoise (Neophocaena phocaenoides) populations in Chinese waters. Intragenic recombination and strong positive selection were the main forces in generating sequence diversity in the DRB gene. MHC sequence diversity changed significantly along the study period. Significant decrease in heterozygosity and lost alleles have been detected in the Yangtze River population and South China Sea population since 1990. Furthermore, there is a trend of increasing population differentiation over time. Especially, the genetic differentiation between the Yangtze River population and the Yellow Sea population was very low prior to 1990 (F (ST) = 0.036, P = 0.009), but became very significant after 1990 (F (ST) = 0.134, P < 0.001), suggesting a recent augmentation of genetic differentiation between both populations probably in a relatively short-term period. Porpoises from the Yangtze River displayed divergent frequencies of shared and private alleles from those displayed by two marine populations, which suggest that the former riverine population has been under a different selection regime (characteristic of a fresh water environment) than that of its marine counterparts.

High MHC DQB variation and asymmetric allelic distribution in the endangered Yangtze finless porpoise, Neophocaena phocaenoides asiaeorientalis

The endangered Yangtze finless porpoise is found in the middle and lower reaches of the Yangtze River and its adjoining big lakes. To explore the major histocompatibility complex (MHC) genetic diversity and allelic distribution patterns across its range, we investigated variation at DQB exon 2. From 76 porpoises, we identified 18 DQB sequences. The freshwater Yangtze populations had much higher allelic diversity than marine populations. Among these freshwater populations, the middle-reach population had higher allelic diversity than the lower-reach population. The high DQB diversity level, relative to that of a neutral mtDNA locus, suggests that balancing selection is acting at the DQB gene and that rapid evolution and local positive selection play critical roles in generating and retaining high MHC diversity in the freshwater population. As the balancing selection might be driven by environmental pathogens, we suggest that maintaining MHC variation should be a high priority in the conservation and management of this endangered population, especially as an ex situ conservation strategy.

Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations

Patterns of spatio-temporal genetic variation at a class II major histocompatibility complex (MHC) locus and multiple microsatellite loci were analysed within and between three water vole metapopulations in Scotland, UK. Comparisons of MHC and microsatellite spatial genetic differentiation, based on standardised tests between two demographically asynchronous zones within a metapopulation, suggested that spatial MHC variation was affected by balancing selection, directional selection and random genetic drift, but that the relative effects of these microevolutionary forces vary temporally. At the metapopulation level, between-year differentiation for MHC loci was significantly correlated with that of microsatellites, signifying that neutral factors such as migration and drift were primarily responsible for overall temporal genetic change at the metapopulation scale. Between metapopulations, patterns of genetic differentiation implied that, at large spatial scales, MHC variation was primarily affected by directional selection and drift. Levels of MHC heterozygosity in excess of Hardy-Weinberg expectations were consistent with overdominant balancing selection operating on MHC variation within metapopulations. However, this effect was not constant among all samples, indicating temporal variation in the strength of selection relative to other factors. The results highlight the benefit of contrasting variation at MHC with neutral markers to separate the effects of stochastic and deterministic microevolutionary forces, and add to a growing body of evidence showing that the mode and relative strength of selection acting on MHC diversity varies both spatially and temporally.

Selection, diversity and evolutionary patterns of the MHC class II DAB in free-ranging Neotropical marsupials

BACKGROUND

Research on the genetic architecture and diversity of the MHC has focused mainly on eutherian mammals, birds and fish. So far, studies on model marsupials used in laboratory investigations indicated very little or even no variation in MHC class II genes. However, natural levels of diversity and selection are unknown in marsupials as studies on wild populations are virtually absent. We used two endemic South American mouse opossums, Gracilinanus microtarsus and Marmosops incanus, to investigate characteristic features of MHC selection. This study is the first investigation of MHC selection in free-ranging Neotropical marsupials. In addition, the evolutionary history of MHC lineages within the group of marsupials was examined.

RESULTS

G. microtarsus showed extensive levels of MHC diversity within and among individuals as 47 MHC-DAB alleles and high levels of sequence divergence were detected at a minimum of four loci. Positively selected codon sites were identified, of which most were congruent with human antigen binding sites. The diversity in M. incanus was rather low with only eight observed alleles at presumably two loci. However, these alleles also revealed high sequence divergence. Again, positive selection was identified on specific codon sites, all congruent with human ABS and with positively selected sites observed in G. microtarsus. In a phylogenetic comparison alleles of M. incanus interspersed widely within alleles of G. microtarsus with four alleles being present in both species.

CONCLUSION

Our investigations revealed extensive MHC class II polymorphism in a natural marsupial population, contrary to previous assumptions. Furthermore, our study confirms for the first time in marsupials the presence of three characteristic features common at MHC loci of eutherian mammals, birds and fish: large allelic sequence divergence, positive selection on specific sites and trans-specific polymorphism.

Clinal variation in MHC diversity with temperature: evidence for the role of host-pathogen interaction on local adaptation in Atlantic salmon

In vertebrates, variability at genes of the Major Histocompatibility Complex (MHC) represents an important adaptation for pathogen resistance, whereby high allelic diversity confers resistance to a greater number of pathogens. Pathogens can maintain diversifying selection pressure on their host's immune system that can vary in intensity based on pathogen richness, pathogen virulence, and length of the cohabitation period, which tend to increase with temperature. In this study, we tested the hypothesis that genetic diversity of MHC increases with temperature along a latitudinal gradient in response to pathogen selective pressure in the wild. A total of 1549 Atlantic salmon from 34 rivers were sampled between 46 degrees N and 58 degrees N in Eastern Canada. The results supported our working hypothesis. In contrast to the overall pattern observed at microsatellites, MHC class II allelic diversity increased with temperature, thus creating a latitudinal gradient. The observed temperature gradient was more pronounced for MHC amino acids of the peptide-binding region (PBR), a region that specifically binds to pathogens, than for the non-PBR. For the subset of rivers analyzed for bacterial diversity, MHC amino acid diversity of the PBR also increased significantly with bacterial diversity in each river. A comparison of the relative influence of temperature and bacterial diversity revealed that the latter could have a predominant role on MHC PBR variability. However, temperature was also identified as an important selective agent maintaining MHC diversity in the wild. Based on the bacteria results and given the putative role of temperature in shaping large-scale patterns of pathogen diversity and virulence, bacterial diversity is a plausible selection mechanism explaining the observed association between temperature and MHC variability. Therefore, we propose that genetic diversity at MHC class II represents local adaptation to cope with pathogen diversity in rivers associated with different thermal regimes. This study illuminates the link between selection pressure from the environment, host immune adaptation, and the large-scale genetic population structure for a nonmodel vertebrate in the wild.

Low MHC variation in the endangered Galápagos penguin (Spheniscus mendiculus)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.

Genetic Variation of the MHCDQBLocus in the Finless Porpoise (Neophocaena phocaenoides)

The Major Histocompatibility Complex (MHC) is a large multigene coding for glycoproteins that play a key role in the initiation of immune responses in vertebrates. The exon 2 region of the MHC DQB locus was analyzed using 160 finless porpoises from 5 populations in Japanese waters. The 5 populations were based on a previous mitochondrial DNA control region analysis, which showed distinct geographical separation. Eight DQB alleles were detected, and the geographical distribution of the alleles indicated that most of them are shared among the populations. Heterozygosity of the DQB alleles in each population ranged from 0.55 to 0.78, and for all 5 populations was 0.78. Low MHC variability is not a common feature in marine mammals, but the finless porpoise populations inhabiting coastal waters had a relatively high MHC heterozygosity. Balancing selection in the MHC DQB alleles of the finless porpoise was indicated by the higher rate of nonsynonymous than synonymous substitutions for PBR; however, an excess of hetrozygotes compared to expectation was not observed. This suggests that the MHC DQB locus in the finless porpoise may have been under balancing selection for a long evolutionary time period, and is influenced by genetic drift beyond the effect of balancing selection for short time periods in small local populations.

Major histocompatibility complex locus DRA polymorphism in the endangered Sorraia horse and related breeds

The major histocompatibility complex (MHC) genes play well-defined roles in eliciting immune responses and combating infectious diseases. This genetic system is among the most polymorphic. The extent of genetic variation within a population has been directly correlated with fitness for many traits. The MHC class II locus DRA polymorphism was analysed in the endangered Sorraia horse, two other Portuguese and four New World horse breeds considered to be historically close to the Sorraia. Comparison of the Sorraia with other breeds demonstrated less MHC variation among Sorraia horses. If DRA polymorphism provides greater disease resistance, selective breeding to increase MHC polymorphism may increase fitness of this population.

Understanding the phylogeographic patterns of European hedgehogs, Erinaceus concolor and E. europaeus using the MHC

The genome of the European hedgehog, Erinaceus concolor and E. europaeus, shows a strong signal of cycles of restriction to glacial refugia and postglacial expansion. Patterns of expansion, however, differ for mitochondrial DNA (mtDNA) and preliminary analysis of nuclear markers. In this study, we determine phylogeographic patterns in the hedgehog using two loci of the major histocompatibility complex (MHC), isolated for the first time in hedgehogs. These genes show long persistence times and high polymorphism in many species because of the actions of balancing selection. Among 84 individuals screened for variation, only two DQA alleles were identified in each species, but 10 DQB alleles were found in E. concolor and six in E. europaeus. A strong effect of demography on patterns of DQB variability is observed, with only weak evidence of balancing selection. While data from mtDNA clearly subdivide both species into monophyletic subgroups, the MHC data delineate only E. concolor into distinct subgroups, supporting the preliminary findings of other nuclear markers. Together with differences in variability, this suggests that the refugia history and/or expansion patterns of E. concolor and E. europaeus differ.

Sequence variation and gene duplication at MHC DQB loci of baiji (Lipotes vexillifer), a Chinese river dolphin

The major histocompatibility complex (MHC) is a fundamental part of the vertebrate immune system, and the high variability in many MHC genes is thought to play an important role in the recognition of parasites. Baiji (Lipotes vexillifer) is one of the most endangered species in the world. Its wild population has declined to fewer than 100 individuals and has a very high risk of becoming extinct in the near future. In this study we present a first step in the molecular characterization of a DQB-like locus of baiji by nucleotide sequence analysis of the polymorphic exon 2 segments. In the examined 172 bp sequences from a group of 18 incidentally captured or stranded individuals, 48 variable sites were determined and 43 alleles were identified, many of which were represented by only one clone. Three to seven alleles were found in each individual, suggesting gene duplications. No deletion, insertion, or exceptional stop codon was detected, suggesting these alleles function in vivo. Phylogenetic reconstruction using neighbor joining grouped the 43 alleles into two distinct lineages, differing by seven nucleotides and four amino acids. Substitutions of amino acids tend to be clustered around sites postulated to be responsible for selective peptide recognition. In the peptide-binding region (PBR) of the DQB locus, the average number of nonsynonymous substitutions per site is greater than that of synonymous substitutions per site (0.1962 versus 0.0256, respectively). Nucleotide and amino acid sequences both showed a relatively high level of similarity (nucleotides 90.6%; amino acids 80.6%) to those of beluga whale (Delphinapterus leucas) and narwhal (Monodon monoceros). The high level of baiji MHC polymorphism revealed in the present study has not been reported in other cetaceans and could be a consequence of the small baiji population adapting to freshwater with a relatively high level of pathogens.

Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae)

The Chatham Island black robin, Petroica traversi, is a highly inbred, endangered passerine with extremely low levels of variation at hypervariable neutral DNA markers. In this study we investigated variation in major histocompatibility complex (MHC) class II genes in both the black robin and its nonendangered relative, the South Island robin Petroica australis australis. Previous studies have shown that Petroica have at least four expressed class II B MHC genes. In this study, the sequences of introns flanking exon 2 of these loci were characterized to design primers for peptide-binding region (PBR) sequence analysis. Intron sequences were comprised of varying numbers of repeated units, with highly conserved regions immediately flanking exon 2. Polymerase chain reaction primers designed to this region amplified three or four sequences per black robin individual, and eight to 14 sequences per South Island robin individual. MHC genes are fitness-related genes thought to be under balancing selection, so they may be more likely to retain variation in bottlenecked populations. To test this, we compared MHC variation in the black robin with artificially bottlenecked populations of South Island robin, and with their respective source populations, using restriction fragment length polymorphism analyses and DNA sequencing of the PBR. Our results indicate that the black robin is monomorphic at class II B MHC loci, while both source and bottlenecked populations of South Island robin have retained moderate levels of variation. Comparison of MHC variation with minisatellite DNA variation indicates that genetic drift outweighs balancing selection in determining MHC diversity in the bottlenecked populations. However, balancing selection appears to influence MHC diversity over evolutionary timescales, and the effects of gene conversion are evident.