INBREEDING INFLUENCES WITHIN‐BROOD HETEROZYGOSITY‐FITNESS CORRELATIONS (HFCS) IN AN ISOLATED PASSERINE POPULATION

Molecular estimates of inbreeding may be made using genetic markers such as microsatellites, however the interpretation of resulting heterozygosity‐fitness correlations (HFCs) with respect to inbreeding depression is not straightforward. We investigated the relationship between pedigree‐determined inbreeding coefficients (f) and HFCs in a closely monitored, reintroduced population of Stewart Island robins (Petroica australis rakiura) on Ulva Island, New Zealand. Using a full sibling design, we focused on differences in juvenile survival associated specifically with individual sibling variation in standardized multilocus heterozygosity (SH) when expected f was identical. We found that within broods, siblings with higher SH at microsatellite loci experienced a higher probability of juvenile survival. This effect, however, was detected primarily within broods that experienced inbreeding or when inbreeding had occurred in their pedigree histories (i.e., at the parents’ level). Thus we show, for the first time in a wild population, that the strength of an HFC is partially dependent on the presence of inbreeding events in the recent pedigree history. Our results illustrate the importance of realized effects of inbreeding on genetic variation and fitness and the value of full‐sibling designs for the study of HFCs in the context of small, inbred populations.     

The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species

In nonpedigreed wild populations, inbreeding depression is often quantified through the use of heterozygosity-fitness correlations (HFCs), based on molecular estimates of relatedness. Although such correlations are typically interpreted as evidence of inbreeding depression, by assuming that the marker heterozygosity is a proxy for genome-wide heterozygosity, theory predicts that these relationships should be difficult to detect. Until now, the vast majority of empirical research in this area has been performed on generally outbred, nonbottlenecked populations, but differences in population genetic processes may limit extrapolation of results to threatened populations. Here, we present an analysis of HFCs, and their implications for the interpretation of inbreeding, in a free-ranging pedigreed population of a bottlenecked species: the endangered takahe (Porphyrio hochstetteri). Pedigree-based inbreeding depression has already been detected in this species. Using 23 microsatellite loci, we observed only weak evidence of the expected relationship between multilocus heterozygosity and fitness at individual life-history stages (such as survival to hatching and fledging), and parameter estimates were imprecise (had high error). Furthermore, our molecular data set could not accurately predict the inbreeding status of individuals (as 'inbred' or 'outbred', determined from pedigrees), nor could we show that the observed HFCs were the result of genome-wide identity disequilibrium. These results may be attributed to high variance in heterozygosity within inbreeding classes. This study is an empirical example from a free-ranging endangered species, suggesting that even relatively large numbers (>20) of microsatellites may give poor precision for estimating individual genome-wide heterozygosity. We argue that pedigree methods remain the most effective method of quantifying inbreeding in wild populations, particularly those that have gone through severe bottlenecks.     

Heterozygosity-fitness correlations in a migratory bird: an analysis of inbreeding and single-locus effects

Studies in a multitude of taxa have described a correlation between heterozygosity and fitness and usually conclude that this is evidence for inbreeding depression. Here, we have used multilocus heterozygosity (MLH) estimates from 15 microsatellite markers to show evidence of heterozygosity-fitness correlations (HFCs) in a long-distance migratory bird, the light-bellied Brent goose. We found significant, positive heterozygosity-heterozygosity correlations between random subsets of the markers we employed, and no evidence that a model containing all loci as individual predictors in a multiple regression explained significantly more variation than a model with MLH as a single predictor. Collectively, these results lend support to the hypothesis that the HFCs we have observed are a function of inbreeding depression. However, we do find that fitness correlations are only detectable in years where population-level productivity is high enough for the reproductive asymmetry between high and low heterozygosity individuals to become apparent. We suggest that lack of evidence of heterozygosity-fitness correlations in animal systems may be because heterozygosity is a poor proxy measure of inbreeding, especially when employing low numbers of markers, but alternatively because the asymmetries between individuals of different heterozygosities may only be apparent when environmental effects on fitness are less pronounced.     

Multilocus heterozygosity and inbreeding in the Siberian jay

Because of its common negative association with fitness, inbreeding is a major concern in conservation biology. Traditionally it has been measured as individual inbreeding coefficient calculated from the pedigree, but recently multilocus heterozygosity estimates have become commonly used as proxies. However, theoretical and simulation studies have cast doubt on the validity of these surrogates especially when they are based on only a few molecular markers. Yet, empirical studies reporting the correlation between multilocus heterozygosity and inbreeding coefficient are rare. We studied this relationship in a wild Siberian jay (Perisoreus infaustus) population subject to a long-term field study over 30 years. The correlations between inbreeding coefficient and the employed heterozygosity measures—standardized heterozygosity and internal relatedness—based on 21 microsatellite loci were weak. These results together with results from theoretical and simulation studies caution against use of multilocus heterozygosity estimates to study inbreeding in natural populations.     

Local heterozygosity-fitness correlations with global positive effects on fitness in threespine stickleback

The complex interactions between genetic diversity and evolution have important implications in many biological areas including conservation, speciation, and mate choice. A common way to study these interactions is to look at heterozygosity-fitness correlations (HFCs). Until recently, HFCs based on noncoding markers were believed to result primarily from global inbreeding effects. However, accumulating theoretical and empirical evidence shows that HFCs may often result from genes being linked to the markers used (local effect). Moreover, local effect HFCs could differ from global inbreeding effects in their direction and occurrence. Consequently, the investigation of the structure and consequences of local HFCs is emerging as a new important goal in evolutionary biology. In this study of a wild threespine stickleback (Gasterosteus aculeatus) population, we first tested the presence of significant positive or negative local effects of heterozygosity at 30 microsatellites loci on five fitness components: survival, mating success, territoriality, length, and body condition. Then, we evaluated the direction and shape of total impact of local HFCs, and estimated the magnitude of the impacts on fitness using regression coefficients and selection differentials. We found that multilocus heterozygosity was not a reliable estimator of individual inbreeding coefficient, which supported the relevance of single-locus based analyses. Highly significant and temporally stable local HFCs were observed. These were mainly positive, but negative effects of heterozygosity were also found. Strong and opposite effects of heterozygosity are probably present in many populations, but may be blurred in HFC analyses looking for global effects only. In this population, both negative and positive HFCs are apparently driving mate preference by females, which is likely to contribute to the maintenance of both additive and nonadditive genetic variance.     

Heterozygosity at a single locus explains a large proportion of variation in two fitness‐related traits in great tits: a general or a local effect?

In natural populations, mating between relatives can have important fitness consequences due to the negative effects of reduced heterozygosity. Parental level of inbreeding or heterozygosity has been also found to influence the performance of offspring, via direct and indirect parental effects that are independent of the progeny own level of genetic diversity. In this study, we first analysed the effects of parental heterozygosity and relatedness (i.e. an estimate of offspring genetic diversity) on four traits related to offspring viability in great tits (Parus major) using 15 microsatellite markers. Second, we tested whether significant heterozygosity–fitness correlations (HFCs) were due to ‘local’ (i.e. linkage to genes influencing fitness) and/or ‘general’ (genome‐wide heterozygosity) effects. We found a significant negative relationship between parental genetic relatedness and hatching success, and maternal heterozygosity was positively associated with offspring body size. The characteristics of the studied populations (recent admixture, polygynous matings) together with the fact that we found evidence for identity disequilibrium across our set of neutral markers suggest that HFCs may have resulted from genome‐wide inbreeding depression. However, one locus (Ase18) had disproportionately large effects on the observed HFCs: heterozygosity at this locus had significant positive effects on hatching success and offspring size. It suggests that this marker may lie near to a functional locus under selection (i.e. a local effect) or, alternatively, heterozygosity at this locus might be correlated to heterozygosity across the genome due to the extensive ID found in our populations (i.e. a general effect). Collectively, our results lend support to both the general and local effect hypotheses and reinforce the view that HFCs lie on a continuum from inbreeding depression to those strictly due to linkage between marker loci and genes under selection.     

Heterozygosity-fitness correlations of conserved microsatellite markers in Kentish plovers Charadrius alexandrinus

Heterozygosity-fitness correlations (HFCs) are frequently used to examine the relationship between genetic diversity and fitness. Most studies have reported positive HFCs, although there is a strong bias towards investigating HFCs in genetically impoverished populations. We investigated HFCs in a large genetically diverse breeding population of Kentish plovers Charadrius alexandrinus in southern Turkey. This small shorebird exhibits highly variable mating and care systems, and it is becoming an ecological model species to understand breeding system evolution. Using 11 conserved and six anonymous microsatellite markers, we tested whether and how heterozygosity was associated with chick survival, tarsus and body mass growth controlling for nongenetic effects (chick sex, hatching date, length of biparental care and site quality) that influence survival and growth. There was no genome-wide effect of heterozygosity on fitness, and we did not find any significant effects of heterozygosity on growth rates. However, two of the 11 conserved markers displayed an association with offspring survival: one marker showed a positive HFC, whereas the other marker showed a negative HFC. Heterozygosity at three further conserved loci showed significant interaction with nongenetic variables. In contrast, heterozygosity based on anonymous microsatellite loci was not associated with fitness or growth. Markers that were correlated with chick survival were not more likely to be located in exons or introns than other markers that lacked this association.     

Adult survival and microsatellite diversity in possums: effects of major histocompatibility complex-linked microsatellite diversity but not multilocus inbreeding estimators

Adult survival is perhaps the fitness parameter most important to population growth in long-lived species. Intrinsic and extrinsic covariates of survival are therefore likely to be important drivers of population dynamics. We used long-term mark-recapture data to identify genetic, individual and environmental covariates of local survival in a natural population of mountain brushtail possums (Trichosurus cunninghami). Rainfall and intra-individual diversity at microsatellite DNA markers were associated with increased local survival of adults and juveniles. We contrasted the performance of several microsatellite heterozygosity measures, including internal relatedness (IR), homozygosity by loci (HL) and the mean multilocus estimate of the squared difference in microsatellite allele sizes within an individual (mean d ²). However, the strongest effect on survival was not associated with multilocus microsatellite diversity (which would indicate a genome-wide inbreeding effect), but a subset of two loci. This included a major histocompatibility complex (MHC)-linked marker and a putatively neutral microsatellite locus. For both loci, diversity measures incorporating allele size information had stronger associations with survival than measures based on heterozygosity, whether or not allele frequency information was included (such as IR). Increased survival was apparent among heterozygotes at the MHC-linked locus, but the benefits of heterozygosity to survival were reduced in heterozygotes with larger differences in allele size. The effect of heterozygosity on fitness-related traits was supported by data on endoparasites in a subset of the individuals studied in this population. There was no apparent density dependence in survival, nor an effect of sex, age or immigrant status. Our findings suggest that in the apparent absence of inbreeding, variation at specific loci can generate strong associations between fitness and diversity at linked markers.     

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

Heterozygosity–fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.     

Heterozygosity-fitness correlations within inbreeding classes: local or genome-wide effects?

Marker-based studies of inbreeding may lead to an enhanced understanding of inbreeding depression in natural populations, which is a major concern in conservation genetics. Correlations between marker heterozygosity and variation in fitness-associated traits—‘heterozygosity-fitness correlations’ (HFCs)—are of particular importance and have been widely applied in natural populations. In partially inbred populations, HFCs can be driven by selection against inbred individuals and thus reflect inbreeding depression. However, other explanations for HFCs also exist, such as functional effects of the markers per se or that the markers reveal selection on linked fitness genes due to extended linkage disequilibrium (LD) in the population. Accordingly, HFCs do not only arise in partially inbred populations, they may also occur within inbreeding classes such as families, i.e. in situations when there is no variation in the inbreeding coefficient. In this study we focus on the importance of LD for within-family HFCs, thereby aiming at enhancing our general understanding of HFCs. For non-coding markers, within-family HFCs have been proposed to be caused in two ways: either by ‘local effects’ at linked fitness genes in LD with the markers, or by ‘general effects’ due to a correlation between proportion of heterozygous markers (H _M ) and heterozygosity at genome-wide distributed fitness genes (H _GW ). To evaluate these contrasting hypotheses for within-family HFCs, we analysed simulated data sets of sexually reproducing populations with varying levels of LD. The results confirmed that segregation induces variation in both H _M and H _GW at a fixed level of inbreeding; as expected, the variation in H _M declined with increasing number of markers, whereas the variation in H _GW declined with decreasing LD. However, less intuitively, there was no positive correlation between the variation in H _M and H _GW within inbreeding classes when the local component of H _GW was accounted for (i.e. when the part of the chromosome in LD with the markers was excluded). This strongly suggests that within-family HFCs are not caused by general effects. Instead, our results support the idea that HFCs at a known level of inbreeding can be driven by local effects in populations with high to moderate LD. Note however that we define the local component of H _GW as the part of the chromosomes in LD with the markers. This implies that when LD is high, the local component will consist of a substantial part of the genome and thus provides a rather ‘genome-wide’ view. We caution against routinely interpreting positive HFCs as evidence of inbreeding depression and non-significant HFCs as lack thereof, especially when few markers are used.