Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.     

Nucleotide diversity in gorillas

Comparison of the levels of nucleotide diversity in humans and apes may provide valuable information for inferring the demographic history of these species, the effect of social structure on genetic diversity, patterns of past migration, and signatures of past selection events. Previous DNA sequence data from both the mitochondrial and the nuclear genomes suggested a much higher level of nucleotide diversity in the African apes than in humans. Noting that the nuclear DNA data from the apes were very limited, we previously conducted a DNA polymorphism study in humans and another in chimpanzees and bonobos, using 50 DNA segments randomly chosen from the noncoding, nonrepetitive parts of the human genome. The data revealed that the nucleotide diversity (pi) in bonobos (0.077%) is actually lower than that in humans (0.087%) and that pi in chimpanzees (0.134%) is only 50% higher than that in humans. In the present study we sequenced the same 50 segments in 15 western lowland gorillas and estimated pi to be 0.158%. This is the highest value among the African apes but is only about two times higher than that in humans. Interestingly, available mtDNA sequence data also suggest a twofold higher nucleotide diversity in gorillas than in humans, but suggest a threefold higher nucleotide diversity in chimpanzees than in humans. The higher mtDNA diversity in chimpanzees might be due to the unique pattern in the evolution of chimpanzee mtDNA. From the nuclear DNA pi values, we estimated that the long-term effective population sizes of humans, bonobos, chimpanzees, and gorillas are, respectively, 10,400, 12,300, 21,300, and 25,200.     

Evidence for a complex demographic history of chimpanzees

To characterize patterns of genomic variation in central chimpanzees(Pan troglodytes troglodytes) and gain insight into their evolution,we sequenced nine unlinked, intergenic regions, representinga total of 19,000 base pairs, in 14 individuals. When theseDNA sequences are compared with homologous sequences previouslycollected in humans and in western chimpanzees (Pan troglodytesverus), nucleotide diversity is higher in central chimpanzeesthan in western chimpanzees or in humans. Consistent with alarger effective population size of central chimpanzees, levelsof linkage disequilibrium are lower than in humans. Patternsof linkage disequilibrium further suggest that homologous geneconversion may be an important contributor to genetic exchangeat short distances, in agreement with a previous study of thesame DNA sequences in humans. In central chimpanzees, but notin western chimpanzees, the allele frequency spectrum is significantlyskewed towards rare alleles, pointing to population size changesor fine-scale population structure. Strikingly, the extent ofgenetic differentiation between western and central chimpanzeesis much stronger than what is seen between human populations.This suggests that careful attention should be paid to geographicsampling in studies of chimpanzee genetic variation.     

Diversification of bitter taste receptor gene family in western chimpanzees

In mammals, bitter taste is mediated by T2R genes, which belong to the large family of seven transmembrane G protein-coupled receptors. Because T2Rs are directly involved in the interaction between mammals and their dietary sources, it is likely that these genes evolved to reflect species' specific diets during mammalian evolution. Here, we investigated the sequences of all 28 putative functional chimpanzee T2R genes (cT2Rs) in 46 western chimpanzees to compare the intraspecies variations in chimpanzees to those already known for all 25 human functional T2R genes (hT2Rs). The numbers of functional genes varied among individuals in western chimpanzees, and most chimpanzees had two or three more functional genes than humans. Similarly to hT2Rs, cT2Rs showed high nucleotide diversity along with a large number of amino acid substitutions. Comparison of the nucleotide substitution patterns in cT2Rs with those in five cT2R pseudogenes and 14 autosomal intergenic noncoding regions among the same individuals revealed that the evolution of cT2R genes was almost identical to that of putative neutral regions with slight but significantly positive Tajima's D values, suggesting that selective constraint on these genes was relaxed with weak balancing selection. These trends have resulted in the occurrence of various divergent alleles of T2Rs within the western chimpanzee populations and in heterozygous individuals who might have the ability to taste a broader range of substances.     

Comparative nuclear and mitochondrial genome diversity in humans and chimpanzees

Restriction mapping and sequencing have shown that humans have substantially lower levels of mitochondrial genome diversity (d) than chimpanzees. In contrast, humans have substantially higher levels of heterozygosity (H) at protein-coding loci, suggesting a higher level of diversity in the nuclear genome. To investigate the discrepancy further, we sequenced a segment of the mitochondrial genome control region (CR) from 49 chimpanzees. The majority of these were from the Pan troglodytes versus subspecies, which was underrepresented in previous studies. We also estimated the average heterozygosity at 60 short tandem repeat (STR) loci in both species. For a total sample of 115 chimpanzees, d = 0.075 +/0 0.037, compared to 0.020 +/- 0.011 for a sample of 1,554 humans. The heterozygosity of human STR loci is significantly higher than that of chimpanzees. Thus, the higher level of nuclear genome diversity relative to mitochondrial genome diversity in humans is not restricted to protein-coding loci. It seems that humans, not chimpanzees, have an unusual d/H ratio, since the ratio in chimpanzees is similar to that in other catarrhines. This discrepancy in the relative levels of nuclear and mitochondrial genome diversity in the two species cannot be explained by differences in mutation rate. However, it may result from a combination of factors such as a difference in the extent of sex ratio disparity, the greater effect of population subdivision on mitochondrial than on nuclear genome diversity, a difference in the relative levels of male and female migration among subpopulations, diversifying selection acting to increase variation in the nuclear genome, and/or directional selection acting to reduce variation in the mitochondrial genome.      

Common chimpanzees have greater diversity than humans at two of the three highly polymorphic MHC class I genes

 MHC class I polymorphism improves the defense of vertebrate species against viruses and other intracellular pathogens. To see how polymorphism at the same class I genes can evolve in different species we compared the MHC-A, MHC-B, and MHC-C loci of common chimpanzees and humans. Diversity in 23 Patr-A, 32 Patr-B, and 18 Patr-C alleles obtained from study of 48 chimpanzees was compared to diversity in 66 HLA-A, 149 HLA-B, and 41 HLA-C alleles obtained from a study of over 1 million humans. At each locus, alleles group hierarchically into families and then lineages. No alleles or families are shared by the two species, commonality being seen only at the lineage level. The overall nucleotide sequence diversity of MHC class I is estimated to be greater for modern chimpanzees than humans. Considering the numbers of lineages, families, and alleles, Patr-B and Patr-C have greater diversity than the HLA-B and HLA-C, respectively. In contrast, Patr-A has less polymorphism than HLA-A, due to the absence of A2 lineage alleles. The results are consistent with ancestral humans having passed through a narrower population bottleneck than chimpanzees, and with pathogen-mediated selection having favored either preservation of A2 lineage alleles on the human line and/or their extinction on the chimpanzee line. Received: 8 December 1999 / Accepted: 30 December 1999     

Human-chimpanzee DNA sequence variation in the four major genes of the renin angiotensin system

The renin angiotensin system (RAS) is involved in blood pressure control and water/sodium metabolism. The genes encoding the proteins of this system are candidate genes for essential hypertension. The RAS involves four main molecules: angiotensinogen, renin, angiotensin I-converting enzyme, and the angiotensin II type 1 receptor (encoded by the genes AGT, REN, DCP1, and AGTR1, respectively). We performed a molecular screening over 17,037 bp of the coding and 5' and 3' untranslated regions of these genes, from three to six common chimpanzees. We identified 44 single-nucleotide polymorphisms (SNPs) in chimpanzee samples, including 18 coding-region SNPs, 5 of which led to an amino acid replacement. We observed common and different features at various sites (synonymous, nonsynonymous, and noncoding) within and between the four chimpanzee genes: (1) the nucleotide diversity at noncoding sites was similar; (2) the nucleotide diversity at nonsynonymous sites was low, probably reflecting purifying selection, except for the AGT gene; (3) the nucleotide diversity at synonymous sites, which was dependent on the G+C content at the third position of the codon, was high, except for the AGTR1 gene. Comparison of the chimpanzee SNPs with those previously reported for humans identified 119 sites with fixed differences (including 62 coding sites, 17 of which resulted in amino acid differences between the species). Analysis of polymorphism within species and divergence between species shed light on the evolutionary constraints on these genes. In particular, comparison of the pattern of mutation at polymorphic and fixed sites between humans and chimpanzees suggested that the high G+C content of the DCP1 gene was maintained by positive selection at its silent sites. Finally, we propose 68 ancestral alleles for the human RAS genes and discuss the implications for their use in future hypertension-susceptibility association studies.     

Modern African ape populations as genetic and demographic models of the last common ancestor of humans, chimpanzees, and gorillas

In order to fully understand human evolutionary history through the use of molecular data, it is essential to include our closest relatives as a comparison. We provide here estimates of nucleotide diversity and effective population size of modern African ape species using data from several independent noncoding nuclear loci, and use these estimates to make predictions about the nature of the ancestral population that eventually gave rise to the living species of African apes, including humans. Chimpanzees, bonobos, and gorillas possess two to three times more nucleotide diversity than modern humans. We hypothesize that the last common ancestor (LCA) of these species had an effective population size more similar to modern apes than modern humans. In addition, estimated dates for the divergence of the Homo, Pan, and Gorilla lineages suggest that the LCA may have had stronger geographic structuring to its mtDNA than its nuclear DNA, perhaps indicative of strong female philopatry or a dispersal system analogous to gorillas, where females disperse only short distances from their natal group. Synthesizing different classes of data, and the inferences drawn from them, allows us to predict some of the genetic and demographic properties of the LCA of humans, chimpanzees, and gorillas.     

Similar numbers but different repertoires of olfactory receptor genes in humans and chimpanzees

Animals recognize their external world through the detection of tens of thousands of chemical odorants. Olfactory receptor (OR) genes encode proteins for detecting odorant molecules and form the largest multigene family in mammals. It is known that humans have fewer OR genes and a higher fraction of OR pseudogenes than mice or dogs. To investigate whether these features are human specific or common to all higher primates, we identified nearly complete sets of OR genes from the chimpanzee and macaque genomes and compared them with the human OR genes. In contrast to previous studies, here we show that the number of OR genes ( approximately 810) and the fraction of pseudogenes (51%) in chimpanzees are very similar to those in humans, though macaques have considerably fewer OR genes. The pseudogenization rates and the numbers of genes affected by positive selection are also similar between humans and chimpanzees. Moreover, the most recent common ancestor between humans and chimpanzees had a larger number of functional OR genes (>500) and a lower fraction of pseudogenes (41%) than its descendents, suggesting that the OR gene repertoires are in a phase of deterioration in both lineages. Interestingly, despite the close evolutionary relationship between the 2 species, approximately 25% of their functional gene repertoires are species specific due to massive gene losses. These findings suggest that the tempo of evolution of OR genes is similar between humans and chimpanzees, but the OR gene repertoires are quite different between them. This difference might be responsible for the species-specific ability of odor perception.     

Characterization of Enteroviruses from Non-Human Primates in Cameroon Revealed Virus Types Widespread in Humans along with Candidate New Types and Species

Enteroviruses (EVs) infecting African Non-Human Primates (NHP) are still poorly documented. This study was designed to characterize the genetic diversity of EVs among captive and wild NHP in Cameroon and to compare this diversity with that found in humans. Stool specimens were collected in April 2008 in NHP housed in sanctuaries in Yaounde and neighborhoods. Moreover, stool specimens collected from wild NHP from June 2006 to October 2008 in the southern rain forest of Cameroon were considered. RNAs purified directly from stool samples were screened for EVs using a sensitive RT-nested PCR targeting the VP1 capsid coding gene whose nucleotide sequence was used for molecular typing. Captive chimpanzees (Pan troglodytes) and gorillas (Gorilla gorilla) were primarily infected by EV types already reported in humans in Cameroon and elsewhere: Coxsackievirus A13 and A24, Echovirus 15 and 29, and EV-B82. Moreover EV-A119, a novel virus type recently described in humans in central and west Africa, was also found in a captive Chimpanzee. EV-A76, which is a widespread virus in humans, was identified in wild chimpanzees, thus suggesting its adaptation and parallel circulation in human and NHP populations in Cameroon. Interestingly, some EVs harbored by wild NHP were genetically distinct from all existing types and were thus assigned as new types. One chimpanzee-derived virus was tentatively assigned as EV-J121 in the EV-J species. In addition, two EVs from wild monkeys provisionally registered as EV-122 and EV-123 were found to belong to a candidate new species. Overall, this study indicates that the genetic diversity of EVs among NHP is more important than previously known and could be the source of future new emerging human viral diseases.     

Nonneutral Mitochondrial DNA Variation in Humans and Chimpanzees

We sequenced the NADH dehydrogenase subunit 3 (ND3) gene from a sample of 61 humans, five common chimpanzees, and one gorilla to test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution. Within humans and within chimpanzees, the ratio of replacement to silent nucleotide substitutions was higher than observed in comparisons between species, contrary to neutral expectations. To test the generality of this result, we reanalyzed published human RFLP data from the entire mitochondrial genome. Gains of restriction sites relative to a known human mtDNA sequence were used to infer unambiguous nucleotide substitutions. We also compared the complete mtDNA sequences of three humans. Both the RFLP data and the sequence data reveal a higher ratio of replacement to silent nucleotide substitutions within humans than is seen between species. This pattern is observed at most or all human mitochondrial genes and is inconsistent with a strictly neutral model. These data suggest that many mitochondrial protein polymorphisms are slightly deleterious, consistent with studies of human mitochondrial diseases.     

Using reporter gene assays to identify cis regulatory differences between humans and chimpanzees

Most phenotypic differences between human and chimpanzee are likely to result from differences in gene regulation, rather than changes to protein-coding regions. To date, however, only a handful of human-chimpanzee nucleotide differences leading to changes in gene regulation have been identified. To hone in on differences in regulatory elements between human and chimpanzee, we focused on 10 genes that were previously found to be differentially expressed between the two species. We then designed reporter gene assays for the putative human and chimpanzee promoters of the 10 genes. Of seven promoters that we found to be active in human liver cell lines, human and chimpanzee promoters had significantly different activity in four cases, three of which recapitulated the gene expression difference seen in the microarray experiment. For these three genes, we were therefore able to demonstrate that a change in cis influences expression differences between humans and chimpanzees. Moreover, using site-directed mutagenesis on one construct, the promoter for the DDA3 gene, we were able to identify three nucleotides that together lead to a cis regulatory difference between the species. High-throughput application of this approach can provide a map of regulatory element differences between humans and our close evolutionary relatives.     

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.     

How chimpanzees look at pictures: a comparative eye-tracking study

Surprisingly little is known about the eye movements of chimpanzees, despite the potential contribution of such knowledge to comparative cognition studies. Here, we present the first examination of eye tracking in chimpanzees. We recorded the eye movements of chimpanzees as they viewed naturalistic pictures containing a full-body image of a chimpanzee, a human or another mammal; results were compared with those from humans. We found a striking similarity in viewing patterns between the two species. Both chimpanzees and humans looked at the animal figures for longer than at the background and at the face region for longer than at other parts of the body. The face region was detected at first sight by both species when they were shown pictures of chimpanzees and of humans. However, the eye movements of chimpanzees also exhibited distinct differences from those of humans; the former shifted the fixation location more quickly and more broadly than the latter. In addition, the average duration of fixation on the face region was shorter in chimpanzees than in humans. Overall, our results clearly demonstrate the eye-movement strategies common to the two primate species and also suggest several notable differences manifested during the observation of pictures of scenes and body forms.     

Genomic analysis of common chimpanzee major histocompatibility complex class I genes

To investigate how MHC class I genes have changed in the approximately 5 million years since chimpanzees and humans diverged, we characterized six genomic fragments ranging in size from 5.1 to 6.1 kb, each containing the complete coding region, introns, and flanking regions of one of the following chimpanzee class I genes: Patr-A, Patr-E, Patr-F, Patr-G, Patr-H, and Patr-J. In humans, these genes are closely linked within the class I region and are representatives of three distinct functional categories of class I genes: the highly polymorphic Ia genes (HLA-A), the conserved Ib genes (HLA-E, HLA-F, and HLA-G), and the class I pseudogenes (HLA-H and HLA-J). Southern blot analysis of chimpanzee and human class I genes produced nearly identical patterns, suggesting that the organization and linkage of these genes differs little in the two species. Comparison of the chimpanzee fragment sequences with their human orthologues revealed structural conservation of these genes yet differences in their degree of functional constraint. This is apparent in the location and nature of the amino acid changes between species and the substantial differences in levels of divergence at functional and nonfunctional sites. Additionally, there is no correlation between patterns of divergence at these sites and intraspecific variation, an observation explained by either appreciable gene conversion or high levels of recombination, the latter unlikely given the observed strong linkage disequilibrium of these loci.     

Natural selection on the olfactory receptor gene family in humans and chimpanzees

The olfactory receptor (OR) genes constitute the largest gene family in mammalian genomes. Humans have >1,000 OR genes, of which only ~40% have an intact coding region and are therefore putatively functional. In contrast, the fraction of intact OR genes in the genomes of the great apes is significantly greater (68%–72%), suggesting that selective pressures on the OR repertoire vary among these species. We have examined the evolutionary forces that shaped the OR gene family in humans and chimpanzees by resequencing 20 OR genes in 16 humans, 16 chimpanzees, and one orangutan. We compared the variation at the OR genes with that at intergenic regions. In both humans and chimpanzees, OR pseudogenes seem to evolve neutrally. In chimpanzees, patterns of variability are consistent with purifying selection acting on intact OR genes, whereas, in humans, there is suggestive evidence for positive selection acting on intact OR genes. These observations are likely due to differences in lifestyle, between humans and great apes, that have led to distinct sensory needs.     

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.     

Positive and negative selection on noncoding DNA close to protein-coding genes in wild house mice

During the past two decades, evidence has accumulated of adaptive evolution within protein-coding genes in a variety of species. However, with the exception of Drosophila and humans, little is known about the extent of adaptive evolution in noncoding DNA. Here, we study regions upstream and downstream of protein-coding genes in the house mouse Mus musculus castaneus, a species that has a much larger effective population size (N(e)) than humans. We analyze polymorphism data for 78 genes from 15 wild-caught M. m. castaneus individuals and divergence to a closely related species, Mus famulus. We find high levels of nucleotide diversity and moderate levels of selective constraint in upstream and downstream regions compared with nonsynonymous sites of protein-coding genes. From the polymorphism data, we estimate the distribution of fitness effects (DFE) of new mutations and infer that most new mutations in upstream and downstream regions behave as effectively neutral and that only a small fraction is strongly negatively selected. We also estimate the fraction of substitutions that have been driven to fixation by positive selection (α) and the ratio of adaptive to neutral divergence (ω(α)). We find that α for upstream and downstream regions (～ 10%) is much lower than α for nonsynonymous sites (～ 50%). However, ω(α) estimates are very similar for nonsynonymous sites (～ 10%) and upstream and downstream regions (～ 5%). We conclude that negative selection operating in upstream and downstream regions of M. m. castaneus is weak and that the low values of α for upstream and downstream regions relative to nonsynonymous sites are most likely due to the presence of a higher proportion of neutrally evolving sites and not due to lower absolute rates of adaptive substitution.     

Human and chimpanzee gene expression differences replicated in mice fed different diets

Although the human diet is markedly different from the diets of closely related primate species, the influence of diet on phenotypic and genetic differences between humans and other primates is unknown. In this study, we analyzed gene expression in laboratory mice fed diets typical of humans and of chimpanzees. The effects of human diets were found to be significantly different from that of a chimpanzee diet in the mouse liver, but not in the brain. Importantly, 10% of the genes that differ in their expression between humans and chimpanzee livers differed also between the livers of mice fed the human and chimpanzee diets. Furthermore, both the promoter sequences and the amino acid sequences of these diet-related genes carry more differences between humans and chimpanzees than random genes. Our results suggest that the mouse can be used to study at least some aspects of human-specific traits.     

Lack of genetic diversity across diverse immune genes in an endangered mammal,the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll‐like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome‐level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole‐genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29–220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long‐term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad‐scale immunogenetic diversity analysis in threatened species.