A robust sequencing assay of a thousand amplicons for the high-throughput population monitoring of Alpine ibex immunogenetics

Genetic variation is a major factor determining susceptibility to diseases. Polymorphism at the major histocompatibility complex (MHC) and other immune function loci can underlie health and reproductive success of individuals. Endangered species of low population size could be severely compromised to evolve disease resistance due to reduced adaptive variation. A major impediment to screen adaptive genetic variation in wild species is the difficulty to comprehensively genotype immune-related loci based on low input material. Here, we design and validate a targeted amplicon sequencing assay to parallelize the analysis of a thousand loci of the MHC, other immunity-related genes, and genome-wide markers for the assessment of population structure. We apply the approach to Alpine ibex, one of the most successful examples of restoration of a large mammal in Europe. We used 51 whole genome sequenced individuals to select representative target SNPs. We integrated SNP call data from four related species for amplification robustness and genotyped 158 Alpine ibex individuals for validation. We show that the genome-wide markers perform equally well at resolving population structure as RAD-seq or low-coverage genome sequencing datasets with orders of magnitude more markers. The targeted amplicon sequencing assay is robust to >100-fold variation in input DNA quantity and generates useful genotype information from fecal samples. The amplicon marker set also identified recent species hybridization events with domestic goats. The immune loci show unexpectedly high degrees of differentiation within the species. Our assay strategy can realistically be implemented into population genetic surveys of a large range of species.