Cooperative social behaviors, such as parental care, have long been hypothesized to relax selection leading to the accumulation of genetic variation in populations. Although the idea has been discussed for decades, there has been relatively little experimental work to investigate how social behavior contributes to genetic variation in populations. Here, we investigate how parental care can shape molecular genetic variation in the subsocial insect, Nicrophorus vespilloides. Using whole-genome sequencing of populations that had evolved in the presence or absence of parental care for 30 generations, we show that parental care maintains levels of standing genetic variation. In contrast, under a harsh environment without care, strong directional selection caused a reduction in genetic variation. Furthermore, we show that adaptation to the loss of care is associated with genetic divergence between populations at loci related to stress, morphological development, and transcriptional regulation. These data reveal how social behavior is linked to the genetic processes that shape and maintain genetic diversity within populations, and provides rare empirical evidence for an old hypothesis.

Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of \"inbreeding,\" calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Dingoes arrived in Australia during the mid-Holocene and are the top-order terrestrial predator on the continent. Although dingoes subsequently spread across the continent, the initial founding population(s) could have been small. We investigated this hypothesis by sequencing the whole genomes of three dingoes and also obtaining the genome data from nine additional dingoes and 56 canines, including wolves, village dogs and breed dogs, and examined the signatures of bottlenecks and founder effects. We found that the nucleotide diversity of dingoes was low, 36% less than highly inbred breed dogs and 3.3 times lower than wolves. The number of runs of homozygosity (RoH) segments in dingoes was 1.6-4.7 times higher than in other canines. While examining deleterious mutational load, we observed that dingoes carried elevated ratios of nonsynonymous-to-synonymous diversities, significantly higher numbers of homozygous deleterious Single Nucleotide Variants (SNVs), and increased numbers of loss of function SNVs, compared to breed dogs, village dogs, and wolves. Our findings can be explained by bottlenecks and founder effects during the establishment of dingoes in mainland Australia. These findings highlight the need for conservation-based management of dingoes and the need for wildlife managers to be cognisant of these findings when considering the use of lethal control measures across the landscape.

Each new human has an expected Ud = 2 - 10 new deleterious mutations. This deluge of deleterious mutations cannot all be purged, and therefore accumulate in a declining fitness ratchet. Using a novel simulation framework designed to efficiently handle genome-wide linkage disequilibria across many segregating sites, we find that rarer, beneficial mutations of larger effect are sufficient to compensate fitness declines due to the fixation of many slightly deleterious mutations. Drift barrier theory posits a similar asymmetric pattern of fixations to explain ratcheting genome size and complexity, but in our theory, the cause is Ud > 1 rather than small population size. In our simulations, Ud ~2 - 10 generates high within-population variance in relative fitness; two individuals will typically differ in fitness by 15-40%. Ud ~2 - 10 also slows net adaptation by ~13%-39%. Surprisingly, fixation rates are more sensitive to changes in the beneficial than the deleterious mutation rate, e.g. a 10% increase in overall mutation rate leads to faster adaptation; this puts to rest dysgenic fears about increasing mutation rates due to rising paternal age.

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Deleterious variants are selected against but can linger in populations at low frequencies for long periods of time, decreasing fitness and contributing to disease burden in humans and other species. Deleterious variants occur at low frequency but distinguishing deleterious variants from low-frequency neutral variation is challenging based on population genomics data alone. As a result, we have little sense of the number and identity of deleterious variants in wild populations. For haplodiploid species, it has been hypothesised that deleterious alleles will be directly exposed to selection in haploid males, but selection can be masked in diploid females when deleterious variants are recessive, resulting in more efficient purging of deleterious mutations in males. Therefore, comparisons of the differences between haploid and diploid genomes from the same population may be a useful method for inferring rare deleterious variants. This study provides the first formal test of this hypothesis. Using wild populations of Northern paper wasps (Polistes fuscatus), we find that males have fewer missense and nonsense variants per generation than females from the same population. Allele frequency differences are especially pronounced for rare missense and nonsense variants and these differences lead to a lower mutational load in males than females. Based on these data we infer that many highly deleterious mutations are segregating in the paper wasp population. Stronger selection against deleterious alleles in haploid males may have implications for adaptation in other haplodiploid insects and provides evidence that wild populations harbour abundant deleterious variants.

Benevolent social behaviours, such as parental care, are thought to enable mildly deleterious mutations to persist. We tested this prediction experimentally using the burying beetle Nicrophorus vespilloides, an insect with biparental care. For 20 generations, we allowed replicate experimental burying beetle populations to evolve either with post-hatching care (\'Full Care\' populations) or without it (\'No Care\' populations). We then established new lineages, seeded from these experimental populations, which we inbred to assess their mutation load. Outbred lineages served as controls. We also tested whether the deleterious effects of a greater mutation load could be concealed by parental care by allowing half the lineages to receive post-hatching care, while half did not. We found that inbred lineages from the Full Care populations went extinct more quickly than inbred lineages from the No Care populations-but only when offspring received no post-hatching care. We infer that Full Care lineages carried a greater mutation load, but that the associated deleterious effects on fitness could be overcome if larvae received parental care. We suggest that the increased mutation load caused by parental care increases a population\'s dependence upon care. This could explain why care is seldom lost once it has evolved.

Mitochondrial DNA (mtDNA) deletions underpin mitochondrial dysfunction in human tissues in aging and disease. The multicopy nature of the mitochondrial genome means these mtDNA deletions can occur in varying mutation loads. At low levels, these deletions have no impact, but once the proportion of molecules harbouring a deletion exceeds a threshold level, then dysfunction occurs. The location of the breakpoints and the size of the deletion impact upon the mutation threshold required to cause deficiency of an oxidative phosphorylation complex, and this varies for each of the different complexes. Furthermore, mutation load and deletion species can vary between adjacent cells in a tissue, with a mosaic pattern of mitochondrial dysfunction observed. As such, it is often important for understanding human aging and disease to be able to characterise the mutation load, breakpoints and size of deletion(s) from a single human cell. Here, we detail protocols for laser micro-dissection and single cell lysis from tissues, and the subsequent analysis of deletion size, breakpoints and mutation load using long-range PCR, mtDNA sequencing and real-time PCR, respectively.

Small populations with limited range are often threatened by inbreeding and reduced genetic diversity, which can reduce fitness and exacerbate population decline. One of the most extreme natural examples is the Devils Hole pupfish (Cyprinodon diabolis), an iconic and critically endangered species with the smallest known range of any vertebrate. This species has experienced severe declines in population size over the last 30 years and suffered major bottlenecks in 2007 and 2013, when the population shrunk to 38 and 35 individuals, respectively. Here, we analysed 30 resequenced genomes of desert pupfishes from Death Valley, Ash Meadows and surrounding areas to examine the genomic consequences of small population size. We found extremely high levels of inbreeding (FROH = 0.34-0.81) and an increased amount of potentially deleterious genetic variation in the Devils Hole pupfish as compared to other species, including unique, fixed loss-of-function alleles and deletions in genes associated with sperm motility and hypoxia. Additionally, we successfully resequenced a formalin-fixed museum specimen from 1980 and found that the population was already highly inbred prior to recent known bottlenecks. We thus document severe inbreeding and increased mutation load in the Devils Hole pupfish and identify candidate deleterious variants to inform management of this conservation icon.