Estimating genome-wide heterozygosity: effects of demographic history and marker type

The difficulty of detecting inbreeding depression and its effect on conservation decisions

Statistical inferences about inbreeding depression are often derived from analyses with low power and a high risk of failing to detect inbreeding depression. That risk is widely appreciated by scientists familiar with the relevant statistical and genetical theory, but may be overlooked and underappreciated by decision-makers. Consequently, there is value in demonstrating this risk using a real example. We use data from the wolf population on Isle Royale to demonstrate the difficulty of making reliable statistical inferences about inbreeding depression. This wolf population is known-by other methods-to have gone effectively extinct due to deleterious genetic processes associated with inbreeding. Beyond that demonstration, we use two case-studies-wolves on Isle Royale and vaquita (porpoises) from the Gulf of California, Mexico-to show how statistical inferences about inbreeding depression can affect conservation decisions. According to most decision theory, decisions depend importantly on: 1) probabilities that certain states exist (e.g. inbreeding depression is present) and 2) the utility assigned to various outcomes (e.g. the value of acting to mitigate inbreeding when it is present). The probabilities are provided by statistical inference; whereas utilities are almost entirely determined by normative values and judgements. Our analysis suggests that decisions to mitigate inbreeding depression are often driven more by utilities (normative values) than probabilities (statistical inferences). As such, advocates for mitigating inbreeding depression will benefit from better communicating to decision-makers the value of populations persisting and the extent to which decisions should depend on normative values.

Comparative Analysis of Microsatellite and SNP Markers for Genetic Management of Red Deer

The analysis of population genetic structure and individual multilocus heterozygosity are crucial for wildlife management and conservation. Microsatellite markers have traditionally been used to assess these genetic parameters. However, single-nucleotide polymorphisms (SNPs) are becoming increasingly popular. Our goal here was to determine to what extent SNPs can provide better insights than microsatellites into the overall genetic status and population genetic processes in the species. To this end, we genotyped 210 red deer (Cervus elaphus) in the Spanish wild population with both 11 microsatellites and 31,712 SNPs. We compared parameters related to population genetic structure and individual multilocus heterozygosity obtained with both types of markers. Our results showed correlations between parameters measured using both microsatellites and SNPs, particularly those related to the level of genetic diversity and genetic differentiation. However, we found notably lower precision of microsatellites in measuring the distribution of genetic diversity among individuals. We conclude that microsatellites can be used to monitor the overall genetic status and detect broad patterns in red deer populations. Nevertheless, the greater precision of SNPs in inferring genetic structure and multilocus heterozygosity leads us to encourage scientists and wildlife managers to prioritize their use whenever possible.

Population variability in X-chromosome inactivation across 9 mammalian species

One of the two X chromosomes in female mammals is epigenetically silenced in embryonic stem cells by X chromosome inactivation (XCI). This creates a mosaic of cells expressing either the maternal or the paternal X allele. The XCI ratio, the proportion of inactivated parental alleles, varies widely among individuals, representing the largest instance of epigenetic variability within mammalian populations. While various contributing factors to XCI variability are recognized, namely stochastic and/or genetic effects, their relative contributions are poorly understood. This is due in part to limited cross-species analysis, making it difficult to distinguish between generalizable or species-specific mechanisms for XCI ratio variability. To address this gap, we measured XCI ratios in nine mammalian species (9,143 individual samples), ranging from rodents to primates, and compared the strength of stochastic models or genetic factors for explaining XCI variability. Our results demonstrate the embryonic stochasticity of XCI is a general explanatory model for population XCI variability in mammals, while genetic factors play a minor role.

Multi-locus homozygosity promotes actuarial senescence in a wild mammal

Genome-wide homozygosity, caused for example by inbreeding, is expected to have deleterious effects on survival and/or reproduction. Evolutionary theory predicts that any fitness costs are likely to be detected in late life because natural selection will filter out negative impacts on younger individuals with greater reproductive value. Here we infer associations between multi-locus homozygosity (MLH), sex, disease and age-dependent mortality risks using Bayesian analysis of the life histories of wild European badgers Meles meles in a population naturally infected with Mycobacterium bovis (the causative agent of bovine tuberculosis [bTB]). We find important effects of MLH on all parameters of the Gompertz-Makeham mortality hazard function, but particularly in later life. Our findings confirm the predicted association between genomic homozygosity and actuarial senescence. Increased homozygosity is particularly associated with an earlier onset, and greater rates of actuarial senescence, regardless of sex. The association between homozygosity and actuarial senescence is further amplified among badgers putatively infected with bTB. These results recommend further investigation into the ecological and behavioural processes that result in genome-wide homozygosity, and focused work on whether homozygosity is harmful or beneficial during early life-stages.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

Evidence of an Absence of Inbreeding Depression in a Wild Population of Weddell Seals (Leptonychotes weddellii)

Inbreeding depression can reduce the viability of wild populations. Detecting inbreeding depression in the wild is difficult; developing accurate estimates of inbreeding can be time and labor intensive. In this study, we used a two-step modeling procedure to incorporate uncertainty inherent in estimating individual inbreeding coefficients from multilocus genotypes into estimates of inbreeding depression in a population of Weddell seals (Leptonychotes weddellii). The two-step modeling procedure presented in this paper provides a method for estimating the magnitude of a known source of error, which is assumed absent in classic regression models, and incorporating this error into inferences about inbreeding depression. The method is essentially an errors-in-variables regression with non-normal errors in both the dependent and independent variables. These models, therefore, allow for a better evaluation of the uncertainty surrounding the biological importance of inbreeding depression in non-pedigreed wild populations. For this study we genotyped 154 adult female seals from the population in Erebus Bay, Antarctica, at 29 microsatellite loci, 12 of which are novel. We used a statistical evidence approach to inference rather than hypothesis testing because the discovery of both low and high levels of inbreeding are of scientific interest. We found evidence for an absence of inbreeding depression in lifetime reproductive success, adult survival, age at maturity, and the reproductive interval of female seals in this population.

Polyandry and non-random fertilisation maintain long-term genetic diversity in an isolated island population of adders (Vipera berus)

Conservation genetic theory suggests that small and isolated populations should be subject to reduced genetic diversity i.e., heterozygosity and allelic diversity. Our 34 years study of an isolated island population of adders (Vipera berus) in southern Sweden challenges this notion. Despite a lack of gene flow and a yearly mean estimated reproductive adult population size of only 65 adult adders (range 12-171), the population maintains high levels of heterozygosity and allelic diversity similar to that observed in two mainland populations. Even a 14-year major "bottleneck" i.e., a reduction in adult adder numbers, encompassing at least four adder generations, did not result in any reduction in the island adders' heterozygosity and allelic diversity. Female adders are polyandrous, and fertilisation is non-random, which our empirical data and modelling suggest are underpinning the maintenance of the population's high level of heterozygosity. Our empirical results and subsequent modelling suggest that the positive genetic effects of polyandry in combination with non-random fertilisation, often overlooked in conservation genetic analyses, deserve greater consideration when predicting long-term survival of small and isolated populations.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

Relationship between genome-wide and MHC class I and II genetic diversity and complementarity in a nonhuman primate

Although mate choice is expected to favor partners with advantageous genetic properties, the relative importance of genome‐wide characteristics, such as overall heterozygosity or kinship, versus specific loci, is unknown. To disentangle genome‐wide and locus‐specific targets of mate choice, we must first understand congruence in global and local variation within the same individual. This study compares genetic diversity, both absolute and relative to other individuals (i.e., complementarity), assessed across the genome to that found at the major histocompatibility complex (MHC), a hyper‐variable gene family integral to immune system function and implicated in mate choice across species. Using DNA from 22 captive olive baboons (Papio anubis), we conducted double digest restriction site‐associated DNA sequencing to estimate genome‐wide heterozygosity and kinship, and sequenced two class I and two class II MHC loci. We found that genome‐wide diversity was not associated with MHC diversity, and that diversity at class I MHC loci was not correlated with diversity at class II loci. Additionally, kinship was a significant predictor of the number of MHC alleles shared between dyads at class II loci. Our results provide further evidence of the strong selective pressures maintaining genetic diversity at the MHC in comparison to other randomly selected sites throughout the genome. Furthermore, our results indicate that class II MHC disassortative mate choice may mediate inbreeding avoidance in this population. Our study suggests that mate choice favoring genome‐wide genetic diversity is not always synonymous with mate choice favoring MHC diversity, and highlights the importance of controlling for kinship when investigating MHC‐associated mate choice.

Ancient hybridization patterns between bighorn and thinhorn sheep

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

Genetic diversity and population structure of a global invader Mayweed chamomile (Anthemis cotula): management implications

Mayweed chamomile (Anthemis cotula) is a globally invasive, troublesome annual weed but knowledge of its genetic diversity, population structure in invaded regions and invasion patterns remains unstudied. Therefore, germplasm from 19 A. cotula populations (sites) from three geographically distinct invaded regions: the Walla Walla Basin (located in southern Washington) and the Palouse (located in both northern Idaho and eastern Washington), Pacific Northwest, USA and Kashmir Valley, India were grown in the greenhouse for DNA extraction and sequencing. A total of 18 829 single-nucleotide polymorphisms were called and filtered for each of 89 samples. Pairwise F ST, Nei's genetic distance, heterozygosity, Wright's inbreeding coefficient (F) and self-fertilization rates were estimated for populations within and among the three regions with a total of 19 populations comprised of 89 individuals. Overall measurements of genetic variation were low but significant among regions, populations and individuals. Despite the weak genetic structure, two main genetic clusters were evident, one comprised of populations from Palouse and Kashmir Valley, the other comprised of populations from the Walla Walla Basin. Significant selfing was observed in populations from the Walla Walla Basin and Palouse but not from Kashmir Valley, indicating that Mayweed chamomile in the Pacific Northwest, USA could persist with low pollinator or pollen donor densities. Although F ST values between the regions indicate Palouse populations are more closely related to Kashmir Valley than to Walla Walla Basin populations, based on Migrate-n analysis, panmixis was the most likely model, suggesting an unrestricted gene flow among all three regions. Our study indicated that Kashmir Valley populations either originated from or shared the origin with the Palouse populations, suggesting human-mediated migration of A. cotula between regions.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.

Genome-wide diversity and habitat underlie fine-scale phenotypic differentiation in the rainbow darter (Etheostoma caeruleum)

Adaptation to environmental change requires that populations harbor the necessary genetic variation to respond to selection. However, dispersal-limited species with fragmented populations and reduced genetic diversity may lack this variation and are at an increased risk of local extinction. In freshwater fish species, environmental change in the form of increased stream temperatures places many cold-water species at-risk. We present a study of rainbow darters (Etheostoma caeruleum) in which we evaluated the importance of genetic variation on adaptive potential and determined responses to extreme thermal stress. We compared fine-scale patterns of morphological and thermal tolerance differentiation across eight sites, including a unique lake habitat. We also inferred contemporary population structure using genomic data and characterized the relationship between individual genetic diversity and stress tolerance. We found site-specific variation in thermal tolerance that generally matched local conditions and morphological differences associated with lake-stream divergence. We detected patterns of population structure on a highly local spatial scale that could not be explained by isolation by distance or stream connectivity. Finally, we showed that individual thermal tolerance was positively correlated with genetic variation, suggesting that sites with increased genetic diversity may be better at tolerating novel stress. Our results highlight the importance of considering intraspecific variation in understanding population vulnerability and stress response.

An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern

BACKGROUND

Use of genomic tools to characterize wildlife populations has increased in recent years. In the past, genetic characterization has been accomplished with more traditional genetic tools (e.g., microsatellites). The explosion of genomic methods and the subsequent creation of large SNP datasets has led to the promise of increased precision in population genetic parameter estimates and identification of demographically and evolutionarily independent groups, as well as questions about the future usefulness of the more traditional genetic tools. At present, few empirical comparisons of population genetic parameters and clustering analyses performed with microsatellites and SNPs have been conducted.

RESULTS

Here we used microsatellite and SNP data generated from Gunnison sage-grouse (Centrocercus minimus) samples to evaluate concordance of the results obtained from each dataset for common metrics of genetic diversity (HO, HE, FIS, AR) and differentiation (FST, GST, DJost). Additionally, we evaluated clustering of individuals using putatively neutral (SNPs and microsatellites), putatively adaptive, and a combined dataset of putatively neutral and adaptive loci. We took particular interest in the conservation implications of any differences. Generally, we found high concordance between microsatellites and SNPs for HE, FIS, AR, and all differentiation estimates. Although there was strong correlation between metrics from SNPs and microsatellites, the magnitude of the diversity and differentiation metrics were quite different in some cases. Clustering analyses also showed similar patterns, though SNP data was able to cluster individuals into more distinct groups. Importantly, clustering analyses with SNP data suggest strong demographic independence among the six distinct populations of Gunnison sage-grouse with some indication of evolutionary independence in two or three populations; a finding that was not revealed by microsatellite data.

CONCLUSION

We demonstrate that SNPs have three main advantages over microsatellites: more precise estimates of population-level diversity, higher power to identify groups in clustering methods, and the ability to consider local adaptation. This study adds to a growing body of work comparing the use of SNPs and microsatellites to evaluate genetic diversity and differentiation for a species of conservation concern with relatively high population structure and using the most common method of obtaining SNP genotypes for non-model organisms.

Individual Heterozygosity Influences Arrival Times and Mating Success of Male Red-Breasted Flycatchers Ficedula parva

We examined the relationship between individual heterozygosity of male Red-breasted Flycatchers (Ficedula parva; a small long-distance migratory, socially monogamous bird species) and their mating success, arrival time and age. Using eight polymorphic microsatellite loci, we found that male heterozygosity is related to both mating success and arrival time, but not to age. Mated and earlier arriving males had higher heterozygosity than later arrivals and bachelors, but we did not find a relationship between age and individual heterozygosity of males. To our knowledge, this is the first study to provide evidence about the relationship between individual genetic diversity and arrival time, thus arrival time could be used as a signal of individual heterozygosity and quality.

Genome-wide association study identified candidate genes controlling continuous storage root formation and bulking in hexaploid sweetpotato

BACKGROUND

Continuous storage root formation and bulking (CSRFAB) in sweetpotato is an important trait from agronomic and biological perspectives. Information about the molecular mechanisms underlying CSRFAB traits is lacking.

RESULTS

Here, as a first step toward understanding the genetic basis of CSRFAB in sweetpotato, we performed a genome-wide association study (GWAS) using phenotypic data from four distinct developmental stages and 33,068 single nucleotide polymorphism (SNP) and insertion-deletion (indel) markers. Based on Bonferroni threshold (p-value < 5 × 10- 7), we identified 34 unique SNPs that were significantly associated with the complex trait of CSRFAB at 150 days after planting (DAP) and seven unique SNPs associated with discontinuous storage root formation and bulking (DCSRFAB) at 90 DAP. Importantly, most of the loci associated with these identified SNPs were located within genomic regions (using Ipomoea trifida reference genome) previously reported for quantitative trait loci (QTL) controlling similar traits. Based on these trait-associated SNPs, 12 and seven candidate genes were respectively annotated for CSRFAB and DCSRFAB traits. Congruent with the contrasting and inverse relationship between discontinuous and continuous storage root formation and bulking, a DCSRFAB-associated candidate gene regulates redox signaling, involved in auxin-mediated lateral root formation, while CSRFAB is enriched for genes controlling growth and senescence.

CONCLUSION

Candidate genes identified in this study have potential roles in cell wall remodeling, plant growth, senescence, stress, root development and redox signaling. These findings provide valuable insights into understanding the functional networks to develop strategies for sweetpotato yield improvement. The markers as well as candidate genes identified in this pioneering research for CSRFAB provide important genomic resources for sweetpotato and other root crops.

Heterozygosity-Fitness Correlations Reveal Inbreeding Depression in Neonatal Body Size in a Critically Endangered Rock Iguana

Inbreeding depression, though challenging to identify in nature, may play an important role in regulating the dynamics of small and isolated populations. Conversely, greater expression of genetic load can enhance opportunities for natural selection. Conditional expression concentrates these opportunities for selection and may lead to failure of detection. This study investigates the possibility for age-dependent expression of inbreeding depression in a critically endangered population of rock iguanas, Cyclura nubila caymanensis. We employ heterozygote-fitness correlations to examine the contributions of individual genetic factors to body size, a fitness-related trait. Nonsignificant reductions in homozygosity (up to 7%) were detected between neonates and individuals surviving past their first year, which may reflect natural absorption of inbreeding effects by this small, fecund population. The majority of variation in neonate body size was attributed to maternal or environmental effects (i.e., clutch identity and incubation length); however, heterozygosity across 22 microsatellite loci also contributed significantly and positively to model predictions. Conversely, effects of heterozygosity on fitness were not detectable when adults were examined, suggesting that inbreeding depression in body size may be age dependent in this taxon. Overall, these findings emphasize the importance of taking holistic, cross-generational approaches to genetic monitoring of endangered populations.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

No evidence of inbreeding depression in fast declining herds of migratory caribou

Identifying inbreeding depression early in small and declining populations is essential for management and conservation decisions. Correlations between heterozygosity and fitness (HFCs) provide a way to identify inbreeding depression without prior knowledge of kinship among individuals. In Northern Quebec and Labrador, the size of two herds of migratory caribou (Rivière-George, RG and Rivière-aux-Feuilles, RAF) has declined by one to two orders of magnitude in the last three decades. This raises the question of a possible increase in inbreeding depression originating from, and possibly contributing to, the demographic decline in those populations. Here, we tested for the association of genomic inbreeding indices (estimated with 22,073 SNPs) with body mass and survival in 400 caribou sampled in RG and RAF herds between 1996 and 2016. We found no association of individual heterozygosity or inbreeding coefficient with body mass or annual survival. Furthermore, those genomic inbreeding indices remained stable over the period monitored. These results suggest that the rapid and intense demographic decline of the herds did not cause inbreeding depression in those populations. Although we found no evidence for HFCs, if demographic decline continues, it is possible that such inbreeding depression would be triggered.

Abiotic genetic adaptation in the Fagaceae

Fagaceae can be found in tropical and temperate regions and contain species of major ecological and economic importance. In times of global climate change, tree populations need to adapt to rapidly changing environmental conditions. The predicted warmer and drier conditions will potentially result in locally maladapted populations. There is evidence that major genera of the Fagaceae are already negatively affected by climate change-related factors such as drought and associated biotic stressors. Therefore, knowledge of the mechanisms underlying adaptation is of great interest. In this review, we summarise current literature related to genetic adaptation to abiotic environmental conditions. We begin with an overview of genetic diversity in Fagaceae species and then summarise current knowledge related to drought stress tolerance, bud burst timing and frost tolerance in the Fagaceae. Finally, we discuss the role of hybridisation, epigenetics and phenotypic plasticity in adaptation.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Length polymorphisms at two candidate genes explain variation of migratory behaviors in blackpoll warblers (Setophaga striata)

Migratory behaviors such as the timing and duration of migration are genetically inherited and can be under strong natural selection, yet we still know very little about the specific genes or molecular pathways that control these behaviors. Studies in candidate genes Clock and Adcyap1 have revealed that both of these loci can be significantly correlated with migratory behaviors in birds, though observed relationships appear to vary across species. We investigated geographic genetic structure of Clock and Adcyap1 in four populations of blackpoll warblers (Setophaga striata), a Neotropical-Nearctic migrant that exhibits geographic variation in migratory timing and duration across its boreal breeding distribution. Further, we used data on migratory timing and duration, obtained from light-level geolocator trackers to investigate candidate genotype-phenotype relationships at the individual level. While we found no geographic structure in either candidate gene, we did find evidence that candidate gene lengths are correlated with five of the six migratory traits. Maximum Clock allele length was significantly and negatively associated with spring arrival date. Minimum Adcyap1 allele length was significantly and negatively associated with spring departure date and positively associated with fall arrival date at the wintering grounds. Additionally, we found a significant interaction between Clock and Adcyap1 allele lengths on both spring and fall migratory duration. Adcyap1 heterozygotes also had significantly shorter migration duration in both spring and fall compared to homozygotes. Our results support the growing body of evidence that Clock and Adcyap1 allele lengths are correlated with migratory behaviors in birds.

From microsatellites to single nucleotide polymorphisms for the genetic monitoring of a critically endangered sturgeon

The use of genetic information is crucial in conservation programs for the establishment of breeding plans and for the evaluation of restocking success. Short tandem repeats (STRs) have been the most widely used molecular markers in such programs, but next-generation sequencing approaches have prompted the transition to genome-wide markers such as single nucleotide polymorphisms (SNPs). Until now, most sturgeon species have been monitored using STRs. The low diversity found in the critically endangered European sturgeon (Acipenser sturio), however, makes its future genetic monitoring challenging, and the current resolution needs to be increased. Here, we describe the discovery of a highly informative set of 79 SNPs using double-digest restriction-associated DNA (ddRAD) sequencing and its validation by genotyping using the MassARRAY system. Comparing with STRs, the SNP panel proved to be highly efficient and reproducible, allowing for more accurate parentage and kinship assignments' on 192 juveniles of known pedigree and 40 wild-born adults. We explore the effectiveness of both markers to estimated relatedness and inbreeding, using simulated and empirical datasets. Interestingly, we found significant correlations between STRs and SNPs at individual heterozygosity and inbreeding that give support to a reasonable representation of whole genome diversity for both markers. These results are useful for the conservation program of A. sturio in building a comprehensive studbook, which will optimize conservation strategies. This approach also proves suitable for other case studies in which highly discriminatory genetic markers are needed to assess parentage and kinship.

Both candidate gene and neutral genetic diversity correlate with parasite resistance in female Mediterranean mouflon

BACKGROUND

Parasite infections can have substantial impacts on population dynamics and are accordingly a key challenge for wild population management. Here we studied genetic mechanisms driving parasite resistance in a large herbivore through a comprehensive approach combining measurements of neutral (16 microsatellites) and adaptive (MHC DRB1 exon 2) genetic diversity and two types of gastrointestinal parasites (nematodes and coccidia).

RESULTS

While accounting for other extrinsic and intrinsic predictors known to impact parasite load, we show that both neutral genetic diversity and DRB1 are associated with resistance to gastrointestinal nematodes. Intermediate levels of multi-locus heterozygosity maximized nematodes resistance, suggesting that both in- and outbreeding depression might occur in the population. DRB1 heterozygosity and specific alleles effects were detected, suggesting the occurrence of heterozygote advantage, rare-allele effects and/or fluctuating selection. On the contrary, no association was detected between genetic diversity and resistance to coccidia, indicating that different parasite classes are impacted by different genetic drivers.

CONCLUSIONS

This study provides important insights for large herbivores and wild sheep pathogen management, and in particular suggests that factors likely to impact genetic diversity and allelic frequencies, including global changes, are also expected to impact parasite resistance.

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness - Implications for brown trout conservation

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction-site-associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half- and full-siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual-level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual-level genotype information, such as quantifying relatedness and individual-level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.

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.

Disentangling the mechanisms of mate choice in a captive koala population

Successful captive breeding programs are crucial to the long-term survival of many threatened species. However, pair incompatibility (breeding failure) limits sustainability of many captive populations. Understanding whether the drivers of this incompatibility are behavioral, genetic, or a combination of both, is crucial to improving breeding programs. We used 28 years of pairing data from the San Diego Zoo koala colony, plus genetic analyses using both major histocompatibility complex (MHC)-linked and non-MHC-linked microsatellite markers, to show that both genetic and non-genetic factors can influence mating success. Male age was reconfirmed to be a contributing factor to the likelihood of a koala pair copulating. This trend could also be related to a pair's age difference, which was highly correlated with male age in our dataset. Familiarity was reconfirmed to increase the probability of a successful copulation. Our data provided evidence that females select mates based on MHC and genome-wide similarity. Male heterozygosity at MHC class II loci was associated with both pre- and post-copulatory female choice. Genome-wide similarity, and similarity at the MHC class II DAB locus, were also associated with female choice at the post-copulatory level. Finally, certain MHC-linked alleles were associated with either increased or decreased mating success. We predict that utilizing a variety of behavioral and MHC-dependent mate choice mechanisms improves female fitness through increased reproductive success. This study highlights the complexity of mate choice mechanisms in a species, and the importance of ascertaining mate choice mechanisms to improve the success of captive breeding programs.

RAD Sequencing and a Hybrid Antarctic Fur Seal Genome Assembly Reveal Rapidly Decaying Linkage Disequilibrium, Global Population Structure and Evidence for Inbreeding

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level by around 400 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organization.

A high-quality pedigree and genetic markers both reveal inbreeding depression for quality but not survival in a cooperative mammal

Inbreeding depression, the reduced fitness of offspring of closely related parents, is commonplace in both captive and wild populations and has important consequences for conservation and mating system evolution. However, because of the difficulty of collecting pedigree and life-history data from wild populations, relatively few studies have been able to compare inbreeding depression for traits at different points in the life cycle. Moreover, pedigrees give the expected proportion of the genome that is identical by descent (IBD_g ) whereas in theory with enough molecular markers realized IBD_g can be quantified directly. We therefore investigated inbreeding depression for multiple life-history traits in a wild population of banded mongooses using pedigree-based inbreeding coefficients (f_ped ) and standardized multilocus heterozygosity (sMLH) measured at 35-43 microsatellites. Within an information theoretic framework, we evaluated support for either f_ped or sMLH as inbreeding terms and used sequential regression to determine whether the residuals of sMLH on f_ped explain fitness variation above and beyond f_ped . We found no evidence of inbreeding depression for survival, either before or after nutritional independence. By contrast, inbreeding was negatively associated with two quality-related traits, yearling body mass and annual male reproductive success. Yearling body mass was associated with f_ped but not sMLH, while male annual reproductive success was best explained by both f_ped and residual sMLH. Thus, our study not only uncovers variation in the extent to which different traits show inbreeding depression, but also reveals trait-specific differences in the ability of pedigrees and molecular markers to explain fitness variation and suggests that for certain traits, genetic markers may capture variation in realized IBD_g above and beyond the pedigree expectation.

Age-dependent, negative heterozygosity-fitness correlations and local effects in an endangered Caribbean reptile, Iguana delicatissima

Inbreeding depression can have alarming impacts on threatened species with small population sizes. Assessing inbreeding has therefore become an important focus of conservation research. In this study, heterozygosity-fitness correlations (HFCs) were measured by genotyping 7 loci in 83 adult and 184 hatchling Lesser Antillean Iguanas, Iguana delicatissima, at a communal nesting site in Dominica to assess the role of inbreeding depression on hatchling fitness and recruitment to the adult population in this endangered species. We found insignificant correlations between multilocus heterozygosity and multiple fitness proxies in hatchlings and adults. Further, multilocus heterozygosity did not differ significantly between hatchlings and adults, which suggests that the survivorship of homozygous hatchlings does not differ markedly from that of their heterozygous counterparts. However, genotypes at two individual loci were correlated with hatching date, a finding consistent with the linkage between specific marker loci and segregating deleterious recessive alleles. These results provide only modest evidence that inbreeding depression influences the population dynamics of I. delicatissima on Dominica.

Genomic consequences of intensive inbreeding in an isolated wolf population

Inbreeding (mating between relatives) is a major concern for conservation as it decreases individual fitness and can increase the risk of population extinction. We used whole-genome resequencing of 97 grey wolves (Canis lupus) from the highly inbred Scandinavian wolf population to identify 'identical-by-descent' (IBD) chromosome segments as runs of homozygosity (ROH). This gave the high resolution required to precisely measure realized inbreeding as the IBD fraction of the genome in ROH (F _ROH). We found a striking pattern of complete or near-complete homozygosity of entire chromosomes in many individuals. The majority of individual inbreeding was due to long IBD segments (>5 cM) originating from ancestors ≤10 generations ago, with 10 genomic regions showing very few ROH and forming candidate regions for containing loci contributing strongly to inbreeding depression. Inbreeding estimated with an extensive pedigree (F _P) was strongly correlated with realized inbreeding measured with the entire genome (r ² = 0.86). However, inbreeding measured with the whole genome was more strongly correlated with multi-locus heterozygosity estimated with as few as 500 single nucleotide polymorphisms, and with F _ROH estimated with as few as 10,000 single nucleotide polymorphisms, than with F _P. These results document in fine detail the genomic consequences of intensive inbreeding in a population of conservation concern.

Population demography and heterozygosity-fitness correlations in natural guppy populations: An examination using sexually selected fitness traits

Heterozygosity-fitness correlations (HFCs) have been examined in a wide diversity of contexts, and the results are often used to infer the role of inbreeding in natural populations. Although population demography, reflected in population-level genetic parameters such as allelic diversity or identity disequilibrium, is expected to play a role in the emergence and detectability of HFCs, direct comparisons of variation in HFCs across many populations of the same species, with different genetic histories, are rare. Here, we examined the relationship between individual microsatellite heterozygosity and a range of sexually selected traits in 660 male guppies from 22 natural populations in Trinidad. Similar to previous studies, observed HFCs were weak overall. However, variation in HFCs among populations was high for some traits (although these variances were not statistically different from zero). Population-level genetic parameters, specifically genetic diversity levels (number of alleles, observed/expected heterozygosity) and measures of identity disequilibrium (g2 and heterozygosity-heterozygosity correlations), were not associated with variation in population-level HFCs. This latter result indicates that these metrics do not necessarily provide a reliable predictor of HFC effect sizes across populations. Importantly, diversity and identity disequilibrium statistics were not correlated, providing empirical evidence that these metrics capture different essential characteristics of populations. A complex genetic architecture likely underpins multiple fitness traits, including those associated with male fitness, which may have reduced our ability to detect HFCs in guppy populations. Further advances in this field would benefit from additional research to determine the demographic contexts in which HFCs are most likely to occur.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri

BACKGROUND

Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq).

RESULTS

Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.

CONCLUSIONS

We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

Genomics advances the study of inbreeding depression in the wild

Inbreeding depression (reduced fitness of individuals with related parents) has long been a major focus of ecology, evolution, and conservation biology. Despite decades of research, we still have a limited understanding of the strength, underlying genetic mechanisms, and demographic consequences of inbreeding depression in the wild. Studying inbreeding depression in natural populations has been hampered by the inability to precisely measure individual inbreeding. Fortunately, the rapidly increasing availability of high-throughput sequencing data means it is now feasible to measure the inbreeding of any individual with high precision. Here, we review how genomic data are advancing our understanding of inbreeding depression in the wild. Recent results show that individual inbreeding and inbreeding depression can be measured more precisely with genomic data than via traditional pedigree analysis. Additionally, the availability of genomic data has made it possible to pinpoint loci with large effects contributing to inbreeding depression in wild populations, although this will continue to be a challenging task in many study systems due to low statistical power. Now that reliably measuring individual inbreeding is no longer a limitation, a major focus of future studies should be to more accurately quantify effects of inbreeding depression on population growth and viability.

Meiotic recombination shapes precision of pedigree- and marker-based estimates of inbreeding

The proportion of an individual's genome that is identical by descent (GWIBD) can be estimated from pedigrees (inbreeding coefficient 'Pedigree F') or molecular markers ('Marker F'), but both estimators come with error. Assuming unrelated pedigree founders, Pedigree F is the expected proportion of GWIBD given a specific inbreeding constellation. Meiotic recombination introduces variation around that expectation (Mendelian noise) and related pedigree founders systematically bias Pedigree F downward. Marker F is an estimate of the actual proportion of GWIBD but it suffers from the sampling error of markers plus the error that occurs when a marker is homozygous without reflecting common ancestry (identical by state). We here show via simulation of a zebra finch and a human linkage map that three aspects of meiotic recombination (independent assortment of chromosomes, number of crossovers and their distribution along chromosomes) contribute to variation in GWIBD and thus the precision of Pedigree and Marker F. In zebra finches, where the genome contains large blocks that are rarely broken up by recombination, the Mendelian noise was large (nearly twofold larger s.d. values compared with humans) and Pedigree F thus less precise than in humans, where crossovers are distributed more uniformly along chromosomes. Effects of meiotic recombination on Marker F were reversed, such that the same number of molecular markers yielded more precise estimates of GWIBD in zebra finches than in humans. As a consequence, in species inheriting large blocks that rarely recombine, even small numbers of microsatellite markers will often be more informative about inbreeding and fitness than large pedigrees.

Identifying environmental correlates of intraspecific genetic variation

Genetic variation is critical to the persistence of populations and their capacity to adapt to environmental change. The distribution of genetic variation across a species' range can reveal critical information that is not necessarily represented in species occurrence or abundance patterns. We identified environmental factors associated with the amount of intraspecific, individual-based genetic variation across the range of a widespread freshwater fish species, the Murray cod Maccullochella peelii. We used two different approaches to statistically quantify the relative importance of predictor variables, allowing for nonlinear relationships: a random forest model and a Bayesian approach. The latter also accounted for population history. Both approaches identified associations between homozygosity by locus and both disturbance to the natural flow regime and mean annual flow. Homozygosity by locus was negatively associated with disturbance to the natural flow regime, suggesting that river reaches with more disturbed flow regimes may support larger, more genetically diverse populations. Our findings are consistent with the hypothesis that artificially induced perennial flows in regulated channels may provide greater and more consistent habitat and reduce the frequency of population bottlenecks that can occur frequently under the highly variable and unpredictable natural flow regime of the system. Although extensive river regulation across eastern Australia has not had an overall positive effect on Murray cod numbers over the past century, regulation may not represent the primary threat to Murray cod survival. Instead, pressures other than flow regulation may be more critical to the persistence of Murray cod (for example, reduced frequency of large floods, overfishing and chemical pollution).

Population genomic structure and linkage disequilibrium analysis of South African goat breeds using genome-wide SNP data

The sustainability of goat farming in marginal areas of southern Africa depends on local breeds that are adapted to specific agro-ecological conditions. Unimproved non-descript goats are the main genetic resources used for the development of commercial meat-type breeds of South Africa. Little is known about genetic diversity and the genetics of adaptation of these indigenous goat populations. This study investigated the genetic diversity, population structure and breed relations, linkage disequilibrium, effective population size and persistence of gametic phase in goat populations of South Africa. Three locally developed meat-type breeds of the Boer (n = 33), Savanna (n = 31), Kalahari Red (n = 40), a feral breed of Tankwa (n = 25) and unimproved non-descript village ecotypes (n = 110) from four goat-producing provinces of the Eastern Cape, KwaZulu-Natal, Limpopo and North West were assessed using the Illumina Goat 50K SNP Bead Chip assay. The proportion of SNPs with minor allele frequencies >0.05 ranged from 84.22% in the Tankwa to 97.58% in the Xhosa ecotype, with a mean of 0.32 ± 0.13 across populations. Principal components analysis, admixture and pairwise FST identified Tankwa as a genetically distinct population and supported clustering of the populations according to their historical origins. Genome-wide FST identified 101 markers potentially under positive selection in the Tankwa. Average linkage disequilibrium was highest in the Tankwa (r(2)  = 0.25 ± 0.26) and lowest in the village ecotypes (r(2) range = 0.09 ± 0.12 to 0.11 ± 0.14). We observed an effective population size of <150 for all populations 13 generations ago. The estimated correlations for all breed pairs were lower than 0.80 at marker distances >100 kb with the exception of those in Savanna and Tswana populations. This study highlights the high level of genetic diversity in South African indigenous goats as well as the utility of the genome-wide SNP marker panels in genetic studies of these populations.

Mate choice for neutral and MHC genetic characteristics in Alpine marmots: different targets in different contexts?

Sexual selection through female mate choice for genetic characteristics has been suggested to be an important evolutionary force maintaining genetic variation in animal populations. However, the genetic targets of female mate choice are not clearly identified and whether female mate choice is based on neutral genetic characteristics or on particular functional loci remains an open question. Here, we investigated the genetic targets of female mate choice in Alpine marmots (Marmota marmota), a socially monogamous mammal where extra‐pair paternity (EPP) occurs. We used 16 microsatellites to describe neutral genetic characteristics and two MHC loci belonging to MHC class I and II as functional genetic characteristics. Our results reveal that (1) neutral and MHC genetic characteristics convey different information in this species, (2) social pairs show a higher MHC class II dissimilarity than expected under random mate choice, and (3) the occurrence of EPP increases when social pairs present a high neutral genetic similarity or dissimilarity but also when they present low MHC class II dissimilarity. Thus, female mate choice is based on both neutral and MHC genetic characteristics, and the genetic characteristics targeted seem to be context dependent (i.e., the genes involved in social mate choice and genetic mate choice differ). We emphasize the need for empirical studies of mate choice in the wild using both neutral and MHC genetic characteristics because whether neutral and functional genetic characteristics convey similar information is not universal.