Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes1

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Although difficult to test empirically, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering" hypothesis which suggests that expansion of the sex-linked region (SLR) reduces the homozygous expression of deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of partially recessive deleterious mutational variation for the evolution of otherwise neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. Both autosomal and SLR-expanding inversions face a race against time: lightly-loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, after which their chances of fixing become negligible. In contrast, initially unloaded inversions eventually become neutral as their deleterious load reaches the same equilibrium as non-inverted haplotypes. Despite the differences in inheritance and indirect selection, SLR-expanding inversions exhibit similar evolutionary dynamics to autosomal inversions over many biologically plausible parameter conditions. Differences emerge when the population average mutation load is quite high; in this case large autosomal inversions that are lucky enough to be mutation-free can rise to intermediate to high frequencies where selection in homozygotes becomes important (Y-linked inversions never appear as homozygous karyotypes); conditions requiring either high mutation rates, highly recessive deleterious mutations, weak selection, or a combination thereof.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counterintuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here, we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify seven deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild nonhuman primates.

Coevolution of longevity and female germline maintenance

An often-overlooked aspect of life-history optimization is the allocation of resources to protect the germline and secure safe transmission of genetic information. While failure to do so renders significant fitness consequences in future generations, germline maintenance comes with substantial costs. Thus, germline allocation should trade off with other life-history decisions and be optimized in accordance with an organism's reproductive schedule. Here, we tested this hypothesis by studying germline maintenance in lines of seed beetle, selected for early (E) or late (L) reproduction for 350 and 240 generations, respectively. Female animals provide maintenance and screening of male gametes in their reproductive tract and oocytes. Here, we reveal the ability of young and aged E- and L-females to provide this form of germline maintenance by mating them to males with ejaculates with artificially elevated levels of protein and DNA damage. We find that germline maintenance in E-females peaks at young age and then declines, while the opposite is true for L-females, in accordance with the age of reproduction in the respective regime. These findings identify the central role of allocation to secure germline integrity in life-history evolution and highlight how females can play a crucial role in mitigating the effects of male germline decisions on mutation rate and offspring quality.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

High genetic load without purging in caribou, a diverse species at risk

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

Background Selection From Unlinked Sites Causes Nonindependent Evolution of Deleterious Mutations

Background selection describes the reduction in neutral diversity caused by selection against deleterious alleles at other loci. It is typically assumed that the purging of deleterious alleles affects linked neutral variants, and indeed simulations typically only treat a genomic window. However, background selection at unlinked loci also depresses neutral diversity. In agreement with previous analytical approximations, in our simulations of a human-like genome with a realistically high genome-wide deleterious mutation rate, the effects of unlinked background selection exceed those of linked background selection. Background selection reduces neutral genetic diversity by a factor that is independent of census population size. Outside of genic regions, the strength of background selection increases with the mean selection coefficient, contradicting the linked theory but in agreement with the unlinked theory. Neutral diversity within genic regions is fairly independent of the strength of selection. Deleterious genetic load among haploid individuals is underdispersed, indicating nonindependent evolution of deleterious mutations. Empirical evidence for underdispersion was previously interpreted as evidence for global epistasis, but we recover it from a non-epistatic model.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the species' designation as critically endangered. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS is unknown. To address these uncertainties, we sequenced and analyzed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Genetic divergence between sites is minimal, and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Despite this, however, demographic inference confirms a severe population bottleneck beginning around the time of WS onset, highlighting the temporal offset that may occur between a population collapse and its potential impact on genetic diversity. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of genetic structure, uniform diversity, and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

Using Computational Simulations to Model Deleterious Variation and Genetic Load in Natural Populations

AbstractDeleterious genetic variation is abundant in wild populations, and understanding the ecological and conservation implications of such variation is an area of active research. Genomic methods are increasingly used to quantify the impacts of deleterious variation in natural populations; however, these approaches remain limited by an inability to accurately predict the selective and dominance effects of mutations. Computational simulations of deleterious variation offer a complementary tool that can help overcome these limitations, although such approaches have yet to be widely employed. In this perspective article, we aim to encourage ecological and conservation genomics researchers to adopt greater use of computational simulations to aid in deepening our understanding of deleterious variation in natural populations. We first provide an overview of the components of a simulation of deleterious variation, describing the key parameters involved in such models. Next, we discuss several approaches for validating simulation models. Finally, we compare and validate several recently proposed deleterious mutation models, demonstrating that models based on estimates of selection parameters from experimental systems are biased toward highly deleterious mutations. We describe a new model that is supported by multiple orthogonal lines of evidence and provide example scripts for implementing this model (https://github.com/ckyriazis/simulations_review).

Genetic Load and Adaptive Potential of a Recovered Avian Species that Narrowly Avoided Extinction

High genetic diversity is a good predictor of long-term population viability, yet some species persevere despite having low genetic diversity. Here we study the genomic erosion of the Seychelles paradise flycatcher (Terpsiphone corvina), a species that narrowly avoided extinction after having declined to 28 individuals in the 1960s. The species recovered unassisted to over 250 individuals in the 1990s and was downlisted from Critically Endangered to Vulnerable in the International Union for the Conservation of Nature Red List in 2020. By comparing historical, prebottleneck (130+ years old) and modern genomes, we uncovered a 10-fold loss of genetic diversity. Highly deleterious mutations were partly purged during the bottleneck, but mildly deleterious mutations accumulated. The genome shows signs of historical inbreeding during the bottleneck in the 1960s, but low levels of recent inbreeding after demographic recovery. Computer simulations suggest that the species long-term small Ne reduced the masked genetic load and made the species more resilient to inbreeding and extinction. However, the reduction in genetic diversity due to the chronically small Ne and the severe bottleneck is likely to have reduced the species adaptive potential to face environmental change, which together with a higher load, compromises its long-term population viability. Thus, small ancestral Ne offers short-term bottleneck resilience but hampers long-term adaptability to environmental shifts. In light of rapid global rates of population decline, our work shows that species can continue to suffer the effect of their decline even after recovery, highlighting the importance of considering genomic erosion and computer modeling in conservation assessments.

Human deleterious mutation rate implies high fitness variance, with declining mean fitness compensated by rarer beneficial mutations of larger effect

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.

Model worms: knowledge gains and risks associated with the use of model species in parasitological research

Model parasite species, whose entire life cycle can be completed in the laboratory and maintained for multiple generations, have played a fundamental role in our understanding of host–parasite interactions. Yet, keeping parasites in laboratory conditions may expose them to unnatural evolutionary pressures, and using laboratory cultures for research is therefore not without limitations. Using 2 widely-used model helminth species, the cestode Hymenolepis diminuta and the nematode Heligmosomoides polygyrus, I illustrate the caution needed when interpreting experimental results on model species. I first review more than 1200 experimental studies published on these species in the past 4 decades, to determine which research areas they have contributed to. This is followed by an examination of the institutional laboratory cultures that have provided the parasites used in these studies. Some of these have persisted for decades and accounted for a substantial proportion of published studies, whereas others have been short-lived. Using information provided by the curators of active cultures, I summarize data on their origins and maintenance conditions. Finally, I discuss how laboratory cultures may have been subject to the influence of evolutionary genetic processes, such as founder effects, genetic drift and inbreeding. I also address the possibility that serial passage through laboratory hosts across multiple generations has exerted artificial selection on several parasite traits, resulting in genetic and phenotypic divergence among laboratory cultures, and between these cultures and natural parasite populations. I conclude with recommendations for the continued usage of laboratory helminth cultures aimed at maximizing their important contribution to parasitological research.

Relative rates of evolution of male-beneficial nuclear compensatory mutations and male-harming Mother's Curse mitochondrial alleles

Mother's Curse alleles represent a significant source of potential male fitness defects. The maternal inheritance of mutations with the pattern of sex-specific fitness effects, s♀>0>s♂, allows Mother's Curse alleles to spread through a population even though they reduce male fitness. Although the mitochondrial genomes of animals contain only a handful of protein-coding genes, mutations in many of these genes have been shown to have a direct effect on male fertility. The evolutionary process of nuclear compensation is hypothesized to counteract the male-limited mitochondrial defects that spread via Mother's Curse. Here we use population genetic models to investigate the evolution of compensatory autosomal nuclear mutations that act to restore the loss of fitness caused by mitochondrial mutation pressures. We derive the rate of male fitness deterioration by Mother's Curse and the rate of restoration by nuclear compensatory evolution. We find that the rate of nuclear gene compensation is many times slower than that of its deterioration by cytoplasmic mutation pressure, resulting in a significant lag in the recovery of male fitness. Thus, the numbers of nuclear genes capable of restoring male mitochondrial fitness defects must be large in order to sustain male fitness in the face of mutation pressures.

Population Rescue through an Increase in the Selfing Rate under Pollen Limitation: Plasticity versus Evolution

AbstractIncreased rates of self-fertilization offer reproductive assurance when plant populations experience pollen limitation, but self-fertilization may reduce fitness by exposing deleterious mutations. If an environmental change responsible for pollen limitation also induces plastic mating system shifts toward self-pollination, the reproductive assurance benefit and inbreeding depression cost of increased self-fertilization occur immediately, while the benefit and cost happen more gradually when increased self-fertilization occur through evolution. I built eco-evolutionary models to explore the demographic and genetic conditions in which higher self-fertilization by plasticity and/or evolution rescues populations, following deficits due to a sudden onset of pollen limitation. Rescue is most likely under an intermediate level of selfing rate increase, either through plasticity or evolution, and this critical level of selfing rate increase is higher under stronger pollen limitation. Generally, rescue is more likely through plasticity than through evolution. Under weak pollen limitation, rescue by enhanced self-fertilization may mainly occur through purging of deleterious mutations rather than reproductive assurance. The selfing rate increase conferring the highest rescue probability is lower when the initial population size is smaller. This article shows the importance of plasticity during plant population rescue and offers insights for future studies of the evolution of mating system plasticity.

Inferring the distributions of fitness effects and proportions of strongly deleterious mutations

The distribution of fitness effects is a key property in evolutionary genetics as it has implications for several evolutionary phenomena including the evolution of sex and mating systems, the rate of adaptive evolution, and the prevalence of deleterious mutations. Despite the distribution of fitness effects being extensively studied, the effects of strongly deleterious mutations are difficult to infer since such mutations are unlikely to be present in a sample of haplotypes, so genetic data may contain very little information about them. Recent work has attempted to correct for this issue by expanding the classic gamma-distributed model to explicitly account for strongly deleterious mutations. Here, we use simulations to investigate one such method, adding a parameter (plth) to capture the proportion of strongly deleterious mutations. We show that plth can improve the model fit when applied to individual species but underestimates the true proportion of strongly deleterious mutations. The parameter can also artificially maximize the likelihood when used to jointly infer a distribution of fitness effects from multiple species. As plth and related parameters are used in current inference algorithms, our results are relevant with respect to avoiding model artifacts and improving future tools for inferring the distribution of fitness effects.

Demographic history and genomic consequences of 10,000 generations of isolation in a wild mammal

Increased human activities caused the isolation of populations in many species-often associated with genetic depletion and negative fitness effects. The effects of isolation are predicted by theory, but long-term data from natural populations are scarce. We show, with full genome sequences, that common voles (Microtus arvalis) in the Orkney archipelago have remained genetically isolated from conspecifics in continental Europe since their introduction by humans over 5,000 years ago. Modern Orkney vole populations are genetically highly differentiated from continental conspecifics as a result of genetic drift processes. Colonization likely started on the biggest Orkney island and vole populations on smaller islands were gradually split off, without signs of secondary admixture. Despite having large modern population sizes, Orkney voles are genetically depauperate and successive introductions to smaller islands resulted in further reduction of genetic diversity. We detected high levels of fixation of predicted deleterious variation compared with continental populations, particularly on smaller islands, yet the fitness effects realized in nature are unknown. Simulations showed that predominantly mildly deleterious mutations were fixed in populations, while highly deleterious mutations were purged early in the history of the Orkney population. Relaxation of selection overall due to benign environmental conditions on the islands and the effects of soft selection may have contributed to the repeated, successful establishment of Orkney voles despite potential fitness loss. Furthermore, the specific life history of these small mammals, resulting in relatively large population sizes, has probably been important for their long-term persistence in full isolation.

Fluctuating selection and the determinants of genetic variation

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

Methodological challenges in the genomic analysis of an endangered mammal population with low genetic diversity

Recently, populations of various species with very low genetic diversity have been discovered. Some of these persist in the long term, but others could face extinction due to accelerated loss of fitness. In this work, we characterize 45 individuals of one of these populations, belonging to the Iberian desman (Galemys pyrenaicus). For this, we used the ddRADseq technique, which generated 1421 SNPs. The heterozygosity values of the analyzed individuals were among the lowest recorded for mammals, ranging from 26 to 91 SNPs/Mb. Furthermore, the individuals from one of the localities, highly isolated due to strong barriers, presented extremely high inbreeding coefficients, with values above 0.7. Under this scenario of low genetic diversity and elevated inbreeding levels, some individuals appeared to be almost genetically identical. We used different methods and simulations to determine if genetic identification and parentage analysis were possible in this population. Only one of the methods, which does not assume population homogeneity, was able to identify all individuals correctly. Therefore, genetically impoverished populations pose a great methodological challenge for their genetic study. However, these populations are of primary scientific and conservation interest, so it is essential to characterize them genetically and improve genomic methodologies for their research.

Severe inbreeding, increased mutation load and gene loss-of-function in the critically endangered Devils Hole pupfish

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 (F_ROH = 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.

Eco-evolutionary extinction and recolonization dynamics reduce genetic load and increase time to extinction in highly inbred populations

Understanding how genetic and ecological effects can interact to shape genetic loads within and across local populations is key to understanding ongoing persistence of systems that should otherwise be susceptible to extinction through mutational meltdown. Classic theory predicts short persistence times for metapopulations comprising small local populations with low connectivity, due to accumulation of deleterious mutations. Yet, some such systems have persisted over evolutionary time, implying the existence of mechanisms that allow metapopulations to avoid mutational meltdown. We first hypothesize a mechanism by which the combination of stochasticity in the numbers and types of mutations arising locally (genetic stochasticity), resulting local extinction, and recolonization through evolving dispersal facilitates metapopulation persistence. We then test this mechanism using a spatially and genetically explicit individual-based model. We show that genetic stochasticity in highly structured metapopulations can result in local extinctions, which can favor increased dispersal, thus allowing recolonization of empty habitat patches. This causes fluctuations in metapopulation size and transient gene flow, which reduces genetic load and increases metapopulation persistence over evolutionary time. Our suggested mechanism and simulation results provide an explanation for the conundrum presented by the continued persistence of highly structured populations with inbreeding mating systems that occur in diverse taxa.

Modeling minimum viable population size with multiple genetic problems of small populations

An important goal for conservation is to define minimum viable population (MVP) sizes for long-term persistence of a species. There is increasing evidence of the role of genetics in population extinction; thus, conservation practitioners are starting to consider the effects of deleterious mutations (DM), in particular the effects of inbreeding depression on fitness. We sought to develop methods to account for genetic problems other than inbreeding depression in MVP estimates, quantify the effect of the interaction of multiple genetic problems on MVP sizes, and find ways to reduce the arbitrariness of time and persistence probability thresholds in MVP analyses. To do so, we developed ecoevolutionary quantitative models to track population size and levels of genetic diversity. We assumed a biallelic multilocus genome with loci under single or multiple, interacting genetic forces. We included mutation-selection-drift balance (for loci with DM) and 3 forms of balancing selection for loci for which variation is lost through genetic drift. We defined MVP size as the lowest population size that avoids an ecoevolutionary extinction vortex. For populations affected by only balancing selection, MVP size decreased rapidly as mutation rates increased. For populations affected by mutation-selection-drift balance, the MVP size increased rapidly. In addition, MVP sizes increased rapidly as the number of loci increased under the same or different selection mechanisms until even arbitrarily large populations could not survive. In the case of fixed number of loci under selection, interaction of genetic problems did not always increase MVP sizes. To further enhance understanding about interaction of genetic problems, there is need for more empirical studies to reveal how different genetic processes interact in the genome.

Challenges in quantifying genome erosion for conservation

Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.

Heterogeneity in viral populations increases the rate of deleterious mutation accumulation

RNA viruses have high mutation rates, with the majority of mutations being deleterious. We examine patterns of deleterious mutation accumulation over multiple rounds of viral replication, with a focus on how cellular coinfection and heterogeneity in viral output affect these patterns. Specifically, using agent-based intercellular simulations we find, in agreement with previous studies, that coinfection of cells by viruses relaxes the strength of purifying selection, and thereby increases the rate of deleterious mutation accumulation. We further find that cellular heterogeneity in viral output exacerbates the rate of deleterious mutation accumulation, regardless of whether this heterogeneity in viral output is stochastic or is due to variation in cellular multiplicity of infection. These results highlight the need to consider the unique life histories of viruses and their population structure to better understand observed patterns of viral evolution.

Genomic Assessment of Cancer Susceptibility in the Threatened Catalina Island Fox (Urocyon littoralis catalinae)

Small effective population sizes raise the probability of extinction by increasing the frequency of potentially deleterious alleles and reducing fitness. However, the extent to which cancers play a role in the fitness reduction of genetically depauperate wildlife populations is unknown. Santa Catalina island foxes (Urocyon littoralis catalinae) sampled in 2007-2008 have a high prevalence of ceruminous gland tumors, which was not detected in the population prior to a recent bottleneck caused by a canine distemper epidemic. The disease appears to be associated with inflammation from chronic ear mite (Otodectes) infections and secondary elevated levels of Staphyloccus pseudointermedius bacterial infections. However, no other environmental factors to date have been found to be associated with elevated cancer risk in this population. Here, we used whole genome sequencing of the case and control individuals from two islands to identify candidate loci associated with cancer based on genetic divergence, nucleotide diversity, allele frequency spectrum, and runs of homozygosity. We identified several candidate loci based on genomic signatures and putative gene functions, suggesting that cancer susceptibility in this population may be polygenic. Due to the efforts of a recovery program and weak fitness effects of late-onset disease, the population size has increased, which may allow selection to be more effective in removing these presumably slightly deleterious alleles. Long-term monitoring of the disease alleles, as well as overall genetic diversity, will provide crucial information for the long-term persistence of this threatened population.

Mutation accumulation opposes polymorphism: supergenes and the curious case of balanced lethals

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes

Many organisms have sex chromosomes with large nonrecombining regions that have expanded stepwise, generating "evolutionary strata" of differentiation. The reasons for this remain poorly understood, but the principal hypotheses proposed to date are based on antagonistic selection due to differences between sexes. However, it has proved difficult to obtain empirical evidence of a role for sexually antagonistic selection in extending recombination suppression, and antagonistic selection has been shown to be unlikely to account for the evolutionary strata observed on fungal mating-type chromosomes. We show here, by mathematical modeling and stochastic simulation, that recombination suppression on sex chromosomes and around supergenes can expand under a wide range of parameter values simply because it shelters recessive deleterious mutations, which are ubiquitous in genomes. Permanently heterozygous alleles, such as the male-determining allele in XY systems, protect linked chromosomal inversions against the expression of their recessive mutation load, leading to the successive accumulation of inversions around these alleles without antagonistic selection. Similar results were obtained with models assuming recombination-suppressing mechanisms other than chromosomal inversions and for supergenes other than sex chromosomes, including those without XY-like asymmetry, such as fungal mating-type chromosomes. However, inversions capturing a permanently heterozygous allele were found to be less likely to spread when the mutation load segregating in populations was lower (e.g., under large effective population sizes or low mutation rates). This may explain why sex chromosomes remain homomorphic in some organisms but are highly divergent in others. Here, we model a simple and testable hypothesis explaining the stepwise extensions of recombination suppression on sex chromosomes, mating-type chromosomes, and supergenes in general.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Despite limited empirical evidence, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering hypothesis" that suggests that expansion of the sex-linked region (SLR) reduces homozygous expression of partially recessive deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of deleterious mutational variation for the evolution of neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. We find that SLR-expanding inversions face a race against time: lightly loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, and must fix before this window of opportunity closes. The outcome of this race is strongly influenced by inversion size, the mutation rate, and the dominance coefficient of deleterious mutations. Yet, small inversions have elevated fixation probabilities relative to neutral expectations for biologically plausible parameter values. Our results demonstrate that deleterious genetic variation can plausibly drive recombination suppression in small steps and would be most consistent with empirical patterns of small evolutionary strata or gradual recombination arrest.

Competitive ability depends on mating system and ploidy level across Capsella species

BACKGROUND AND AIMS

Self-fertilization is often associated with ecological traits corresponding to the ruderal strategy, and selfers are expected to be less competitive than outcrossers, either because of a colonization/competition trade-off or because of the deleterious genetic effects of selfing. Range expansion could reduce further competitive ability while polyploidy could mitigate the effects of selfing. If pollinators are not limited, individual fitness is thus expected to be higher in outcrossers than in selfers and, within selfers, in polyploids than in diploids. Although often proposed in the botanical literature and also suggested by meta-analyses, these predictions have not been directly tested yet.

METHODS

In order to compare fitness and the competitive ability of four Capsella species with a different mating system and ploidy level, we combined two complementary experiments. First, we carried out an experiment outdoors in north-west Greece, i.e. within the range of the obligate outcrossing species, C. grandiflora, where several life history traits were measured under two different disturbance treatments, weeded plots vs. unweeded plots. To better control competition and to remove potential effects of local adaptation of the outcrosser, we also performed a similar competition experiment but under growth chamber conditions.

KEY RESULTS

In the outdoor experiment, disturbance of the environment did not affect the phenotype in any of the four species. For most traits, the obligate outcrossing species performed better than all selfing species. In contrast, polyploids did not survive or reproduce better than diploids. Under controlled conditions, as in the field experiment, the outcrosser had a higher fitness than selfing species and was less affected by competition. Finally, contrary to the outdoor experiment where the two behaved identically, polyploid selfers were less affected by competition than diploid selfes.

CONCLUSIONS

In the Capsella genus, selfing induces lower fitness than outcrossing and can also reduce competitive ability. The effect of polyploidy is, however, unclear. These results highlight the possible roles of ecological context in the evolution of selfing species.

Accounting for the genetic load in assisted reproductive technology

The genetic load in the human genome has important ramifications for assisted reproductive technology (ART), human reproduction and fertility more generally. Here, we discuss these topics in the light of evolutionary genetic theory, the technological revolution in ART and the advances in the fields of genomics and bioinformatics.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Obligate sexual reproduction of a homothallic fungus closely related to the Cryptococcus pathogenic species complex

<title>eLife digest</title>.

Fungi are enigmatic organisms that flourish in soil, on decaying plants, or during infection of animals or plants. Growing in myriad forms, from single-celled yeast to multicellular molds and mushrooms, fungi have also evolved a variety of strategies to reproduce. Normally, fungi reproduce in one of two ways: either they reproduce asexually, with one individual producing a new individual identical to itself, or they reproduce sexually, with two individuals of different 'mating types' contributing to produce a new individual. However, individuals of some species exhibit 'homothallism' or self-fertility: these individuals can produce reproductive cells that are universally compatible, and therefore can reproduce sexually with themselves or with any other cell in the population.

Homothallism has evolved multiple times throughout the fungal kingdom, suggesting it confers advantage when population numbers are low or mates are hard to find. Yet some homothallic fungi been overlooked compared to heterothallic species, whose mating types have been well characterised. Understanding the genetic basis of homothallism and how it evolved in different species can provide insights into pathogenic species that cause fungal disease.

With that in mind, Passer, Clancey et al. explored the genetic basis of homothallism in Cryptococcus depauperatus, a close relative of C. neoformans, a species that causes fungal infections in humans. A combination of genetic sequencing techniques and experiments were applied to analyse, compare, and manipulate C. depauperatus' genome to see how this species evolved self-fertility.

Passer, Clancey et al. showed that C. depauperatus evolved the ability to reproduce sexually by itself via a unique evolutionary pathway. The result is a form of homothallism never reported in fungi before. C. depauperatus lost some of the genes that control mating in other species of fungi, and acquired genes from the opposing mating types of a heterothallic ancestor to become self-fertile.

Passer, Clancey et al. also found that, unlike other Cryptococcus species that switch between asexual and sexual reproduction, C. depauperatus grows only as long, branching filaments called hyphae, a sexual form. The species reproduces sexually with itself throughout its life cycle and is unable to produce a yeast (asexual) form, in contrast to other closely related species.

This work offers new insights into how different modes of sexual reproduction have evolved in fungi. It also provides another interesting case of how genome plasticity and evolutionary pressures can produce similar outcomes, homothallism, via different evolutionary paths. Lastly, assembling the complete genome of C. depauperatus will foster comparative studies between pathogenic and non-pathogenic Cryptococcus species.

Utilizing evolutionary conservation to detect deleterious mutations and improve genomic prediction in cassava

Introduction

Cassava (Manihot esculenta) is an annual root crop which provides the major source of calories for over half a billion people around the world. Since its domestication ~10,000 years ago, cassava has been largely clonally propagated through stem cuttings. Minimal sexual recombination has led to an accumulation of deleterious mutations made evident by heavy inbreeding depression.

Methods

To locate and characterize these deleterious mutations, and to measure selection pressure across the cassava genome, we aligned 52 related Euphorbiaceae and other related species representing millions of years of evolution. With single base-pair resolution of genetic conservation, we used protein structure models, amino acid impact, and evolutionary conservation across the Euphorbiaceae to estimate evolutionary constraint. With known deleterious mutations, we aimed to improve genomic evaluations of plant performance through genomic prediction. We first tested this hypothesis through simulation utilizing multi-kernel GBLUP to predict simulated phenotypes across separate populations of cassava.

Results

Simulations showed a sizable increase of prediction accuracy when incorporating functional variants in the model when the trait was determined by<100 quantitative trait loci (QTL). Utilizing deleterious mutations and functional weights informed through evolutionary conservation, we saw improvements in genomic prediction accuracy that were dependent on trait and prediction.

Conclusion

We showed the potential for using evolutionary information to track functional variation across the genome, in order to improve whole genome trait prediction. We anticipate that continued work to improve genotype accuracy and deleterious mutation assessment will lead to improved genomic assessments of cassava clones.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo's range were evaluated in this study. Using available microsatellite data (2-17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Populations, Traits, and Their Spatial Structure in Humans

The spatial distribution of genetic variants is jointly determined by geography, past demographic processes, natural selection, and its interplay with environmental variation. A fraction of these genetic variants are "causal alleles" that affect the manifestation of a complex trait. The effect exerted by these causal alleles on complex traits can be independent or dependent on the environment. Understanding the evolutionary processes that shape the spatial structure of causal alleles is key to comprehend the spatial distribution of complex traits. Natural selection, past population size changes, range expansions, consanguinity, assortative mating, archaic introgression, admixture, and the environment can alter the frequencies, effect sizes, and heterozygosities of causal alleles. This provides a genetic axis along which complex traits can vary. However, complex traits also vary along biogeographical and sociocultural axes which are often correlated with genetic axes in complex ways. The purpose of this review is to consider these genetic and environmental axes in concert and examine the ways they can help us decipher the variation in complex traits that is visible in humans today. This initiative necessarily implies a discussion of populations, traits, the ability to infer and interpret "genetic" components of complex traits, and how these have been impacted by adaptive events. In this review, we provide a history-aware discussion on these topics using both the recent and more distant past of our academic discipline and its relevant contexts.

The Alignment of Natural and Sexual Selection

Sexual selection has the potential to decrease mean fitness in a population through an array of costs to nonsexual fitness. These costs may be offset when sexual selection favors individuals with high nonsexual fitness, causing the alignment of sexual and natural selection. We review the many laboratory experiments that have manipulated mating systems aimed at quantifying the net effects of sexual selection on mean fitness. These must be interpreted in light of population history and the diversity of ways manipulations have altered sexual interactions, sexual conflict, and sexual and natural selection. Theory and data suggest a net benefit is more likely when sexually concordant genetic variation is enhanced and that ecological context can mediate the relative importance of these different effects. Comparative studies have independently examined the consequences of sexual selection for population/species persistence. These provide little indication of a benefit, and interpreting these higher-level responses is challenging.

Strong sexual selection fails to protect against inbreeding-driven extinction in a moth

Sexual selection is predicted to influence population persistence because skew in male reproductive success may facilitate the purging of mutation load. We manipulated the strength of sexual selection in populations of Indian meal moths, Plodia interpunctella, by adjusting adult sex ratios to be either male- or female-biased, leading to strong and weak sexual selection in males, respectively. After between 19 and 22 generations of experimental evolution, we examined whether mutation load differed between these populations by enforcing successive generations of inbreeding, tracking extinction events, offspring viability and assaying the effect of inbreeding on male mating success and female choice. We found no effect of the strength of sexual selection on the rate of extinction or offspring viability. We did, however, find changes in both male mating success and female choice, with both being influenced by the sex ratio treatment and the number of generations of inbreeding. Males from male-biased populations were more successful at mating with stock females, and mating success declined rapidly with inbreeding regardless of sex ratio treatment. Females from male-biased populations were less likely to mate with stock males at the onset of the experiment, but tended to mate more frequently with increasing inbreeding compared to females from female-biased populations. Our results demonstrate that while mating behaviors have diverged between male-biased and female-biased lines mutation loads remained similar. This suggests that the benefits of sexual selection to population fitness may be low or slow to accumulate under the benign environmental conditions in which these populations evolved.

Overt and concealed genetic loads revealed by QTL mapping of genotype-dependent viability in the Pacific oyster Crassostrea gigas

Understanding the genetic bases of inbreeding depression, heterosis, and genetic load is integral to understanding how genetic diversity is maintained in natural populations. The Pacific oyster Crassostrea gigas, like many long-lived plants, has high fecundity and high early mortality (type-III survivorship), manifesting a large, overt, genetic load; the oyster harbors an even greater concealed genetic load revealed by inbreeding. Here, we map viability QTL (vQTL) in six interrelated F2 oyster families, using high-density linkage maps of single nucleotide polymorphisms generated by genotyping-by-sequencing (GBS) methods. Altogether, we detect 70 vQTL and provisionally infer 89 causal mutations, 11 to 20 per family. Genetic mortality caused by independent (unlinked) vQTL ranges from 94.2% to 97.8% across families, consistent with previous reports. High-density maps provide better resolution of genetic mechanisms, however. Models of one causal mutation present in both identical-by-descent (IBD) homozygotes and heterozygotes fit genotype frequencies at 37 vQTL; consistent with the mutation-selection balance theory of genetic load, 20 are highly deleterious, completely recessive mutations and 17 are less deleterious, partially dominant mutations. Another 22 vQTL require pairs of recessive or partially dominant causal mutations, half showing selection against recessive mutations linked in repulsion, producing pseudo-overdominance. Only eight vQTL appear to support the overdominance theory of genetic load, with deficiencies of both IBD homozygotes, but at least four of these are likely caused by pseudo-overdominance. Evidence for epistasis is absent. A high mutation rate, random genetic drift, and pseudo-overdominance may explain both the oyster's extremely high genetic diversity and a high genetic load maintained primarily by mutation-selection balance.

Secular changes in human reproduction and assisted reproductive technologies

Since the middle to late 20th century the majority of children born in the developing world have been likely to enter into post-reproductive age. Currently, child mortality is at its lowest level in human history. While more children are living to post reproductive age, approximately 15% of couples are experiencing infecundity. This is either a result of one or both members of the couple being infecund, or, despite both being fecund, the interaction between them prevents fertility for some reason. Assisted reproductive technologies have provided many infertile couples an opportunity to have children. Assisted reproductive technologies operate by intervening and manipulating gametic and intrauterine natural selection. This paper discusses the possible influence of assisted reproductive technologies on child development. This paper outlines some of the reported changes in children resulting from assisted reproductive technologies. Although, few people are either aware or care about possible long term consequences of relaxed natural selection contributed by medical intervention (i.e. assisted reproductive technologies) we have little understanding to what extent such medical interference may affect long term fitness in humans.

Multilevel selection favors fragmentation modes that maintain cooperative interactions in multispecies communities

Reproduction is one of the requirements for evolution and a defining feature of life. Yet, across the tree of life, organisms reproduce in many different ways. Groups of cells (e.g., multicellular organisms, colonial microbes, or multispecies biofilms) divide by releasing propagules that can be single-celled or multicellular. What conditions determine the number and size of reproductive propagules? In multicellular organisms, existing theory suggests that single-cell propagules prevent the accumulation of deleterious mutations (e.g., cheaters). However, groups of cells, such as biofilms, sometimes contain multiple metabolically interdependent species. This creates a reproductive dilemma: small daughter groups, which prevent the accumulation of cheaters, are also unlikely to contain the species diversity that is required for ecological success. Here, we developed an individual-based, multilevel selection model to investigate how such multi-species groups can resolve this dilemma. By tracking the dynamics of groups of cells that reproduce by fragmenting into smaller groups, we identified fragmentation modes that can maintain cooperative interactions. We systematically varied the fragmentation mode and calculated the maximum mutation rate that communities can withstand before being driven to extinction by the accumulation of cheaters. We find that for groups consisting of a single species, the optimal fragmentation mode consists of releasing single-cell propagules. For multi-species groups we find various optimal strategies. With migration between groups, single-cell propagules are favored. Without migration, larger propagules sizes are optimal; in this case, group-size dependent fissioning rates can prevent the accumulation of cheaters. Our work shows that multi-species groups can evolve reproductive strategies that allow them to maintain cooperative interactions.

Larval development in the Pacific oyster and the impacts of ocean acidification: Differential genetic effects in wild and domesticated stocks

The adaptive capacity of marine calcifiers to ocean acidification (OA) is a topic of great interest to evolutionary biologists and ecologists. Previous studies have provided evidence to suggest that larval resilience to high pCO2 seawater for these species is a trait with a genetic basis and variability in natural populations. To date, however, it remains unclear how the selective effects of OA occur within the context of complex genetic interactions underpinning larval development in many of the most vulnerable taxa. Here we evaluated phenotypic and genetic changes during larval development of Pacific oysters (Crassostrea gigas) reared in ambient (~400 µatm) and high (~1600 µatm) pCO2 conditions, both in domesticated and naturalized "wild" oysters from the Pacific Northwest, USA. Using pooled DNA samples, we determined changes in allele frequencies across larval development, from early "D-stage" larvae to metamorphosed juveniles (spat), in both groups and environments. Domesticated larvae had ~26% fewer loci with changing allele frequencies across developmental stages and <50% as many loci affected by acidified culture conditions, compared to larvae from wild broodstock. Functional enrichment analyses of genetic markers with significant changes in allele frequency revealed that the structure and function of cellular membranes were disproportionately affected by high pCO2 conditions in both groups. These results indicate the potential for a rapid adaptive response of oyster populations to OA conditions; however, underlying genetic changes associated with larval development differ between these wild and domesticated oyster stocks and influence their adaptive responses to OA conditions.

Selection in males purges the mutation load on female fitness

Theory predicts that the ability of selection and recombination to purge mutation load is enhanced if selection against deleterious genetic variants operates more strongly in males than females. However, direct empirical support for this tenet is limited, in part because traditional quantitative genetic approaches allow dominance and intermediate-frequency polymorphisms to obscure the effects of the many rare and partially recessive deleterious alleles that make up the main part of a population's mutation load. Here, we exposed the partially recessive genetic load of a population of Callosobruchus maculatus seed beetles via successive generations of inbreeding, and quantified its effects by measuring heterosis-the increase in fitness experienced when masking the effects of deleterious alleles by heterozygosity-in a fully factorial sex-specific diallel cross among 16 inbred strains. Competitive lifetime reproductive success (i.e., fitness) was measured in male and female outcrossed F1s as well as inbred parental "selfs," and we estimated the 4 × 4 male-female inbred-outbred genetic covariance matrix for fitness using Bayesian Markov chain Monte Carlo simulations of a custom-made general linear mixed effects model. We found that heterosis estimated independently in males and females was highly genetically correlated among strains, and that heterosis was strongly negatively genetically correlated to outbred male, but not female, fitness. This suggests that genetic variation for fitness in males, but not in females, reflects the amount of (partially) recessive deleterious alleles segregating at mutation-selection balance in this population. The population's mutation load therefore has greater potential to be purged via selection in males. These findings contribute to our understanding of the prevalence of sexual reproduction in nature and the maintenance of genetic variation in fitness-related traits.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

Mistranslation reduces mutation load in evolving proteins through negative epistasis with DNA mutations

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population’s mutation load, mistranslation can affect an important determinant of evolvability.

Mutational Load and the Functional Fraction of the Human Genome

The fraction of the human genome that is functional is a question of both evolutionary and practical importance. Studies of sequence divergence have suggested that the functional fraction of the human genome is likely to be no more than ∼15%. In contrast, the ENCODE project, a systematic effort to map regions of transcription, transcription factor association, chromatin structure, and histone modification, assigned function to 80% of the human genome. In this article, we examine whether and how an analysis based on mutational load might set a limit on the functional fraction. In order to do so, we characterize the distribution of fitness of a large, finite, diploid population at mutation-selection equilibrium. In particular, if mean fitness is ∼1, the fitness of the fittest individual likely to occur cannot be unreasonably high. We find that at equilibrium, the distribution of log fitness has variance nus, where u is the per-base deleterious mutation rate, n is the number of functional sites (and hence incorporates the functional fraction f), and s is the selection coefficient of deleterious mutations. In a large (N=109) reproducing population, the fitness of the fittest individual likely to exist is ∼e5nus. These results apply to both additive and recessive fitness schemes. Our approach is different from previous work that compared mean fitness at mutation-selection equilibrium with the fitness of an individual who has no deleterious mutations; we show that such an individual is exceedingly unlikely to exist. We find that the functional fraction is not very likely to be limited substantially by mutational load, and that any such limit, if it exists, depends strongly on the selection coefficients of new deleterious mutations.

Measurements of lethal and nonlethal inbreeding depression inform the de novo domestication of Silphium integrifolium

PREMISE

Inbreeding depression, or the reduction in fitness of progeny with related parents, has the potential to adversely affect the long-term viability of both wild and captive plant populations. Silphium integrifolium, a prairie plant native to the central United States, has been identified as a potential candidate for domestication as a perennial oilseed crop. Little is known about the potential for inbreeding depression in this species, but it is expected to be nonnegligible because S. integrifolium is both perennial and self-incompatible. Here, we measure lethal inbreeding depression expressed through embryo deaths, and nonlethal inbreeding depression expressed through changes in vigor and fitness phenotypes of progeny.

METHODS

First, we made controlled crosses among related and unrelated individuals to determine the effect of two different levels of inbreeding on seed production. Then, we grew inbred and outbred progeny from this population to reproductive maturity and measured 11 key traits.

RESULTS

We found that within an improved S. integrifolium population, individuals carried an average of slightly less than one lethal allele per gamete. In progeny, significant inbreeding depression was observed in at least one family for eight of the 11 measured traits.

CONCLUSIONS

Inbreeding depression is likely to be an important challenge to S. integrifolium domestication, reducing overall population fecundity and values for important phenotypes. These effects may grow worse as selection reduces effective population size. We recommend several strategies for S. integrifolium breeding to help mitigate these problems.

Does geographic variation in thermal tolerance in Daphnia represent trade-offs or conditional neutrality?

Geographic variation in thermal tolerance in Daphnia seems to represent genetic load at the loci specifically responsible for heat tolerance resulting from conditional neutrality. We see no evidence of trade-offs between fitness-related traits at 25 °C vs. 10 °C or between two algal diets across Daphnia magna clones from a variety of locations representing the opposite ends of the distribution of long-term heat tolerance. Likewise, we found no evidence of within-environment trade-offs between heat tolerance and fitness-related traits in any of the environments. Neither short-term and long-term heat tolerance shows any consistent relationship with lipid fluorescence polarization and lipid peroxidation across clones or environments. Pervasive positive correlations between fitness-related traits indicate differences in genetic load rather than trade-off based local adaptation or thermal specialization. For heat tolerance such differences may be caused by either relaxation of stabilizing selection due to lower exposure to high temperature extremes, i.e., conditional neutrality, or by small effective population size followed by the recent range expansion.