Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus

Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g₂. We calculated g₂ together with F_ST and pairwise relatedness (F_ij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise F_ij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g₂. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g₂ is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g₂, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity.