Immunogenetic Diversity Research Articles

Classical HLA alleles are associated with prevalent and persistent cervical high-risk HPV infection in African women

BackgroundPersistent cervical high-risk human papillomavirus (hrHPV) infection is a necessary cause of cervical cancer. However, the host genetic factors underlying its risk are not well understood. We hypothesized that immunogenetic variation plays a role in hrHPV infection and persistence. Therefore, we conducted a study of classical HLA alleles and their association with hrHPV infection and persistence among women. MethodsWe characterized HPV infection using SPF10/LiPA25in Nigerian women at baseline and at 6 months follow-up visits in 2014. hrHPV infection was prevalent if at least one carcinogenic HPV genotype was detected at the baseline visit and persistent if at least one carcinogenic HPV genotype was detected at the baseline and follow-up visits. Classical HLA alleles were imputed from genotypes in the MHC region using the HLA genotype imputation with attribute bagging (HIBAG) algorithm. HLA association tests were conducted under additive genetic models. ResultsThe mean (±SD) age of the 517 study participants was 38 (±8) years, 48% were HIV negative, 24% were hrHPV positive at baseline and 10% had persistent hrHPV infections. In multivariate regression models adjusted for age, HIV status and the first principal component, DQA1*01:02 and DQA1*02:01 were positively associated with prevalent but not persistent hrHPV infections, while DQA1*05:01 was negatively associated with prevalent hrHPV but positively associated with persistent cervical hrHPV infections. Four haplotypes (A*30:01-DQA1*05:01, B*07:02-C*07:02, B*07:02-DQA1*05:01 and C*07:02-DQA1*05:01) were significantly associated with prevalent cervical hrHPV infections and several haplotypes that included the DQA1*05:01 allelic variant were significantly associated with persistent cervical hrHPV infections. Six amino acid positions on DQα1 were associated with prevalent but not persistent cervical hrHPV infections. ConclusionsIn this first study to investigate the association between HLA alleles and persistent hrHPV in African women, we identified important risk alleles that merit further investigation. Our findings provide new insights into risk factors for hrHPV infection in African ancestry women.

PC-10 Classical HLA alleles are associated with prevalent and persistent cervical high-risk HPV infection in African women

Background: Persistent cervical high-risk human papillomavirus (hrHPV) infection is a necessary cause of cervical cancer. However, the genetic factors underlying its risk are not well understood. We hypothesized that immunogenetic variation plays a role in hrHPV infection and persistence. Therefore, we conducted a study of classical HLA alleles and their association with hrHPV infection and persistence. Methods: We characterized HPV infection using SPF25/LiPA10 in Nigerian women at baseline and at 6 months follow-up visits in 2014. hrHPV infection was prevalent if at least one carcinogenic HPV genotype was detected in a sample provided at the baseline visit and persistent if at least one carcinogenic HPV genotype was detected in each of the samples provided at the baseline and follow-up visits. Genotyping was performed with the Illumina Multi-Ethnic Genotyping Array. Classical HLA alleles were imputed from genotypes in the MHC region using the HLA genotype imputation with attribute bagging (HIBAG) algorithm. HLA association tests were conducted under additive genetic models. P values were adjusted for multiple comparisons with the false discovery rate (FDR) method. Results: The mean (±SD) age of the study participants was 38 (±8) years, 48% (247/517) were HIV negative, 24% (125/517) were hrHPV positive at baseline and 10% (51/517) had persistent hrHPV infections. In multivariate regression models adjusted for age, HIV status and the first principal component, DQA1*01:02 (OR 1.46; 95% CI: 1.07 to 1.97; P = 0.03). and DQA1*02:01 (OR 1.89; 95% CI: 1.12 to 3.15; P = 0.03) were positively associated with prevalent but not persistent hrHPV infections, while DQA1*05:01 was negatively associated with prevalent hrHPV (OR 0.45; 95% CI: 0.24 to 0.85; P = 0.03) but positively associated with persistent cervical hrHPV infections (OR 2.46; 95% CI: 1.27 to 4.77; P = 0.04). Four haplotypes (A*30:01 - DQA1*05:01, B*07:02 - C*07:02, B*07:02 - DQA1*05:01 and C*07:02 - DQA1*05:01) were significantly associated with prevalent cervical hrHPV infections and several haplotypes that included the DQA1*05:01 allelic variant (including B*15:10 - DQA1*05:01, DQA1*05:01 - DQB1*02:01 and DQA1*05:01 - DRB1*03:01) were significantly associated with persistent cervical hrHPV infections. Six amino acids located on the DQA1 gene were associated with prevalent but not persistent cervical hrHPV infections. Conclusions: This study, the first to investigate the association between HLA alleles and persistent hrHPV risk in African women, has identified important risk alleles that merit further study. Our findings provide new insights into risk factors for hrHPV infection in women of African ancestry.

Implications of management on immunogenetic variation in the European grayling (Thymallus thymallus)

Implications of management on immunogenetic variation in the European grayling (Thymallus thymallus)

Next-generation sequencing reveals new insights about gene usage and CDR-H3 composition in the horse antibody repertoire

Horse serum antibodies have been used for greater than a century for the treatment and prophylaxis of infectious diseases and envenomations. Little is known, however, about the immunogenetic diversity that produces horse serum antibodies. Here, we employed next-generation sequencing for a first-in-kind comprehensive analysis of the equine B-cell repertoire. Nearly 45,000 and 30,000 clonotypes were obtained for the heavy-chain (IGH) and lambda light-chain (IGL) loci, respectively. We observed skewed use of the common subgroups IGHV2 (92.49%) and IGLV8 (82.50%), consistent with previous reports, but also novel use of the rare genes IGHV6S1 and IGLV4S2. CDR-H3 amino acid composition revealed different amino acid patterns at positions 106 and 116 compared to human, rabbit, and mouse, suggesting that an extended conformation predominates among horse CDR-H3 loops. Our analysis provides new insights regarding the mechanisms employed to generate antibody diversity in the horse, and could be applicable to the optimized design of synthetic antibodies intended for future therapeutic use.

Avian Toll-like receptor allelic diversity far exceeds human polymorphism: an insight from domestic chicken breeds

Immune genes show remarkable levels of adaptive variation shaped by pathogen-mediated selection. Compared to humans, however, population polymorphism in animals has been understudied. To provide an insight into immunogenetic diversity in birds, we sequenced complete protein-coding regions of all Toll-like receptor (TLR) genes with direct orthology between mammals and birds (TLR3, TLR4, TLR5 and TLR7) in 110 domestic chickens from 25 breeds and compared their variability with a corresponding human dataset. Chicken TLRs (chTLRs) exhibit on average nine-times higher nucleotide diversity than human TLRs (hTLRs). Increased potentially functional non-synonymous variability is found in chTLR ligand-binding ectodomains. While we identified seven sites in chTLRs under positive selection and found evidence for convergence between alleles, no selection or convergence was detected in hTLRs. Up to six-times more alleles were identified in fowl (70 chTLR4 alleles vs. 11 hTLR4 alleles). In chTLRs, high numbers of alleles are shared between the breeds and the allelic frequencies are more equal than in hTLRs. These differences may have an important impact on infectious disease resistance and host-parasite co-evolution. Though adaptation through high genetic variation is typical for acquired immunity (e.g. MHC), our results show striking levels of intraspecific polymorphism also in poultry innate immune receptors.

Female assortative mate choice functionally validates synthesized male odours of evolving stickleback river-lake ecotypes.

During mate choice decisions, females of many vertebrates use male olfactory cues to achieve immunogenetic optimality of their offspring. Three-spined sticklebacks (Gasterosteus aculeatus) populating habitats that differ in their parasite communities evolve locally adapted combinations of genetic variants encoded at the major histocompatibility complex (MHC). Such adaptation confers optimal resistance to the local parasite fauna. Immunogenetic signatures co-evolved with local parasites favour population-specific assortative mate choice behaviour. Previous studies have shown that female sticklebacks evaluate male MHC-associated olfactory cues during the process of mate choice, but how habitat-specific information is exchanged between males and females has remained elusive. Here, we directly demonstrate the molecular nature of the olfactory cue providing habitat-specific information. Under controlled laboratory conditions, females that are ready to mate prefer mixtures of synthetic MHC peptide ligands mimicking the optimal allele number of their original population. These results imply that female sticklebacks can determine the number of MHC alleles of their prospective mates, compare it to their own immunogenetic status, and, if optimal with respect to the immunogenetic complementarity, accept the male as mate. Our results suggest a potentially common mechanism of ecological speciation in vertebrates that is based on the olfactory assessment of habitat-specific immunogenetic diversity.

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.

Neutral variation does not predict immunogenetic variation in the European grayling (Thymallus thymallus)-implications for management.

Preservation of genetic diversity is critical to successful conservation, and there is increasing demand for the inclusion of ecologically meaningful genetic information in management decisions. Supportive breeding programmes are increasingly implemented to combat declines in many species, yet their effect on adaptive genetic variation is understudied. This is despite the fact that supportive breeding may interfere with natural evolutionary processes. Here, we assessed the performance of neutral and adaptive markers (major histocompatibility complex; MHC) to inform management of European grayling (Thymallus thymallus), which routinely involves supplementation of natural populations with hatchery-reared fish (stocking). This study is the first to characterize MH II DAA and DAB loci in grayling and to investigate immune genetic variation in relation to management practice in this species. High-throughput Illumina sequencing of "introduced," "stocked native" and "non-stocked native" populations revealed significantly higher levels of allelic richness and heterozygosity for MH markers than microsatellites exclusively in non-stocked native populations. Likewise, significantly lower differentiation at the MH II than for microsatellites was apparent when considering non-stocked native populations, but not stocked populations. We developed a simulation model to test the effects of relaxation of selection during the early life stage within captivity. Dependent on the census population size and stocking intensity, there may be long-term effects of stocking on MH II, but not neutral genetic diversity. This is consistent with our empirical results. This study highlights the necessity for considering adaptive genetic variation in conservation decisions and raises concerns about the efficiency of stocking as a management practice.

Development of a genotype-by-sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure.

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.

Profiling the immunome of little brown myotis provides a yardstick for measuring the genetic response to white-nose syndrome.

White‐nose syndrome (WNS) has devastated populations of hibernating bats in eastern North America, leading to emergency conservation listings for several species including the previously ubiquitous little brown myotis (Myotis lucifugus). However, some bat populations near the epicenter of the WNS panzootic appear to be stabilizing after initial precipitous declines, which could reflect a selective immunogenetic sweep. To investigate the hypothesis that WNS exerts significant selection on the immunome of affected bat populations, we developed a novel, high‐throughput sequence capture assay targeting 138 adaptive, intrinsic, and innate immunity genes of putative adaptive significance, as well as their respective regulatory regions (~370 kbp of genomic sequence/individual). We used the assay to explore baseline immunogenetic variation in M. lucifugus and to investigate whether particular immune genes/variants are associated with WNS susceptibility. We also used our assay to detect 1,038 putatively neutral single nucleotide polymorphisms and characterize contemporary population structure, providing context for the identification of local immunogenetic adaptation. Sequence capture provided a cost‐effective, “all‐in‐one” assay to test for neutral genetic and immunogenetic structure and revealed fine‐scale, baseline immunogenetic differentiation between sampling sites <600 km apart. We identified functional immunogenetic variants in M. lucifugus associated with WNS susceptibility. This study lays the foundations for future investigations of rangewide immunogenetic adaptation to WNS in M. lucifugus and provides a blueprint for studies of evolutionary rescue in other host–pathogen systems.

Prelude to a panzootic: Gene flow and immunogenetic variation in northern little brown myotis vulnerable to bat white-nose syndrome

The fungus that causes bat white-nose syndrome (WNS) recently leaped from eastern North America to the Pacific Coast. The pathogen’s spread is associated with the genetic population structure of a host ( Myotis lucifugus). To understand the fine-scale neutral and immunogenetic variation among northern populations of M. lucifugus, we sampled 1142 individuals across the species’ northern range. We used genotypes at 11 microsatellite loci to reveal the genetic structure of, and directional gene flow among, populations to predict the likely future spread of the pathogen in the northwest and to estimate effective population size ( Ne). We also pyrosequenced the DRB1-like exon 2 of the class II major histocompatibility complex (MHC) in 160 individuals to explore immunogenetic selection by WNS. We identified three major neutral genetic clusters: Eastern, Montane Cordillera (and adjacent sampling areas), and Haida Gwaii, with admixture at intermediate areas and significant substructure west of the prairies. Estimates of Newere unexpectedly low (289–16 000). Haida Gwaii may provide temporary refuge from WNS, but the western mountain ranges are not barriers to its dispersal in M. lucifugus and are unlikely to slow its spread. Our major histocompatibility complex (MHC) data suggest potential selection by WNS on the MHC, but gene duplication limited the immunogenetic analyses.

Correction: Prelude to a panzootic: Gene flow and immunogenetic variation in northern little brown myotis vulnerable to bat white-nose syndrome

Correction: Prelude to a panzootic: Gene flow and immunogenetic variation in northern little brown myotis vulnerable to bat white-nose syndrome

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.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout.

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

Toll-like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the 'ghost of selection past' and quantify the 'drift debt'.

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll-like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent-based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the 'ghost of selection past'. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such 'drift debt' occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.

Linking innate immunogenetic variation with phenotypic traits in a wild population of tree swallows, Tachycineta bicolor

Linking innate immunogenetic variation with phenotypic traits in a wild population of tree swallows, Tachycineta bicolor

Reciprocal translocation of small numbers of inbred individuals rescues immunogenetic diversity.

Genetic rescue can reduce inbreeding depression and increase fitness of small populations, even when the donor populations are highly inbred. In a recent experiment involving two inbred island populations of the New Zealand South Island robin, Petroica australis, reciprocal translocations improved microsatellite diversity and individual fitness. While microsatellite loci may reflect patterns of genome-wide diversity, they generally do not indicate the specific genetic regions responsible for increased fitness. We tested the effectiveness of this reciprocal translocation for rescuing diversity of two immunogenetic regions: Toll-like receptor (TLR) and major histocompatibility complex (MHC) genes. We found that the relatively small number of migrants (seven and ten per island) effectively brought the characteristic TLR gene diversity of each source population into the recipient population. However, when migrants transmitted TLR alleles that were already present at high frequency in the recipient population, it was possible for offspring of mixed heritage to have decreased gene diversity compared to recipient population diversity prior to translocation. In contrast to TLRs, we did not observe substantial changes in MHC allelic diversity following translocation, with limited evidence of a decrease in differentiation, perhaps because most MHC alleles were observed at both sites prior to the translocation. Overall, we conclude that small numbers of migrants may successfully restore the diversity of immunogenetic loci with few alleles, but that translocating larger numbers of animals would provide additional opportunity for the genetic rescue of highly polymorphic immunity regions, such as the MHC, even when the source population is inbred.

Genome-wide standing variation facilitates long-term response to bidirectional selection for antibody response in chickens

BackgroundLong-term selection experiments provide a powerful approach to gain empirical insights into adaptation, allowing researchers to uncover the targets of selection and infer their contributions to the mode and tempo of adaptation. Here we implement a pooled genome re-sequencing approach to investigate the consequences of 39 generations of bidirectional selection in White Leghorn chickens on a humoral immune trait: antibody response to sheep red blood cells.ResultsWe observed wide genome involvement in response to this selection regime. Many genomic regions were highly differentiated resulting from this experimental selection regime, an involvement of up to 20% of the chicken genome (208.8 Mb). While genetic drift has certainly contributed to this, we implement gene ontology, association analysis and population simulations to increase our confidence in candidate selective sweeps. Three strong candidate genes, MHC, SEMA5A and TGFBR2, are also presented.ConclusionsThe extensive genomic changes highlight the polygenic genetic architecture of antibody response in these chicken populations, which are derived from a common founder population, demonstrating the extent of standing immunogenetic variation available at the onset of selection.

Lack of diversity at innate immunity Toll-like receptor genes in the Critically Endangered White-winged Flufftail (Sarothrura ayresi).

The White-winged Flufftail (Sarothrura ayresi) population is listed as globally Critically Endangered. White-winged Flufftails are only known to occur, with any regularity, in the high-altitude wetlands of South Africa and Ethiopia. Threats to the species include the limited number of suitable breeding sites in Ethiopia and severe habitat degradation and loss both in Ethiopia and South Africa. Toll-like receptors (TLRs) are increasingly being studied in a variety of taxa as a broader approach to determine functional genetic diversity. In this study, we confirm low genetic diversity in the innate immune regions of the White-winged Flufftail similar to that observed in other bird species that have undergone population bottlenecks. Low TLR diversity in White-winged Flufftail indicates that this species is more likely to be threatened by changes to the environment that would potentially expose the species to new diseases. Thus, conservation efforts should be directed towards maintaining pristine habitat for White-winged Flufftail in its current distribution range. To date, no studies on immunogenetic variation in White-winged Flufftail have been conducted and to our knowledge, this is the first study of TLR genetic diversity in a critically endangered species.

Diversity in the Toll-Like Receptor Genes of the African Penguin (Spheniscus demersus).

The African penguin, Spheniscus demersus, is listed as Endangered by the IUCN Red List of Threatened Species due to the drastic reduction in population numbers over the last 20 years. To date, the only studies on immunogenetic variation in penguins have been conducted on the major histocompatibility complex (MHC) genes. It was shown in humans that up to half of the genetic variability in immune responses to pathogens are located in non-MHC genes. Toll-like receptors (TLRs) are now increasingly being studied in a variety of taxa as a broader approach to determine functional genetic diversity. In this study, we confirm low genetic diversity in the innate immune region of African penguins similar to that observed in New Zealand robin that has undergone several severe population bottlenecks. Single nucleotide polymorphism (SNP) diversity across TLRs varied between ex situ and in situ penguins with the number of non-synonymous alterations in ex situ populations (n = 14) being reduced in comparison to in situ populations (n = 16). Maintaining adaptive diversity is of vital importance in the assurance populations as these animals may potentially be used in the future for re-introductions. Therefore, this study provides essential data on immune gene diversity in penguins and will assist in providing an additional monitoring tool for African penguin in the wild, as well as to monitor diversity in ex situ populations and to ensure that diversity found in the in situ populations are captured in the assurance populations.