Rates and fitness consequences of new mutations in humans

Polygenic selection to a changing optimum under self-fertilisation

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

Effects of urban-induced mutations on ecology, evolution and health

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

Habitat remediation followed by managed connectivity reduces unwanted changes in evolutionary trajectory of high extirpation risk populations

As we continue to convert green spaces into roadways and buildings, connectivity between populations and biodiversity will continue to decline. In threatened and endangered species, this trend is particularly concerning because the cessation of immigration can cause increased inbreeding and loss of genetic diversity, leading to lower adaptability and higher extirpation probabilities in these populations. Unfortunately, monitoring changes in genetic diversity from management actions such as assisted migration and predicting the extent of introduced genetic variation that is needed to prevent extirpation is difficult and costly in situ. Therefore, we designed an agent-based model to link population-wide genetic variability and the influx of unique alleles via immigration to population stability and extirpation outcomes. These models showed that management of connectivity can be critical in restoring at-risk populations and reducing the effects of inbreeding depression. However, the rescued populations were more similar to the migrant source population (average FST range 0.05-0.10) compared to the historical recipient population (average FST range 0.23-0.37). This means that these management actions not only recovered the populations from the effects of inbreeding depression, but they did so in a way that changed the evolutionary trajectory that was predicted and expected for these populations prior to the population crash. This change was most extreme in populations with the smallest population sizes, which are representative of critically endangered species that could reasonably be considered candidates for restored connectivity or translocation strategies. Understanding how these at-risk populations change in response to varying management interventions has broad implications for the long-term adaptability of these populations and can improve future efforts for protecting locally adapted allele complexes when connectivity is restored.

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are genetic novelties pivotal in mammalian evolution. Their phenotypic impacts and evolutionary pattern over time, however, remain elusive in humans due to the technical and ethical complexities in functional studies. By combining human gene age dating and Mendelian disease phenotyping, our research reveals a gradual increase in disease gene proportions with gene age. Logistic regression modeling indicates that this increase could be related to longer protein lengths and higher burdens of deleterious de novo germline variants (DNVs) for older genes. We also find a steady integration of new genes with biomedical phenotypes into the human genome over macroevolutionary timescales (~0.07% per million years). Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures across gene ages. Notably, young genes show significant enrichment in diseases related to the male reproductive system, indicating strong sexual selection. Young genes also exhibit disease-related functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, musculoskeletal phenotypes, and color vision. We further reveal a logistic growth pattern of pleiotropy over evolutionary time, indicating a diminishing marginal growth of new functions for older genes due to intensifying selective constraints over time. We propose a "pleiotropy-barrier" model that delineates higher potentials of phenotypic innovation for young genes than for older genes, a process subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

Constraining models of dominance for nonsynonymous mutations in the human genome

Dominance is a fundamental parameter in genetics, determining the dynamics of natural selection on deleterious and beneficial mutations, the patterns of genetic variation in natural populations, and the severity of inbreeding depression in a population. Despite this importance, dominance parameters remain poorly known, particularly in humans or other non-model organisms. A key reason for this lack of information about dominance is that it is extremely challenging to disentangle the selection coefficient (s) of a mutation from its dominance coefficient (h). Here, we explore dominance and selection parameters in humans by fitting models to the site frequency spectrum (SFS) for nonsynonymous mutations. When assuming a single dominance coefficient for all nonsynonymous mutations, we find that numerous h values can fit the data, so long as h is greater than ~0.15. Moreover, we also observe that theoretically-predicted models with a negative relationship between h and s can also fit the data well, including models with h=0.05 for strongly deleterious mutations. Finally, we use our estimated dominance and selection parameters to inform simulations revisiting the question of whether the out-of-Africa bottleneck has led to differences in genetic load between African and non-African human populations. These simulations suggest that the relative burden of genetic load in non-African populations depends on the dominance model assumed, with slight increases for more weakly recessive models and slight decreases shown for more strongly recessive models. Moreover, these results also demonstrate that models of partially recessive nonsynonymous mutations can explain the observed severity of inbreeding depression in humans, bridging the gap between molecular population genetics and direct measures of fitness in humans. Our work represents a comprehensive assessment of dominance and deleterious variation in humans, with implications for parameterizing models of deleterious variation in humans and other mammalian species.

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).

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are structural novelties pivotal in mammalian evolution. Their phenotypic impact on humans, however, remains elusive due to the technical and ethical complexities in functional studies. Through combining gene age dating with Mendelian disease phenotyping, our research reveals that new genes associated with disease phenotypes steadily integrate into the human genome at a rate of ~ 0.07% every million years over macroevolutionary timescales. Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures between young and old genes. Notably, young genes show significant enrichment in the male reproductive system, indicating strong sexual selection. Young genes also exhibit functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, bipedal locomotion, and color vision. Our findings further reveal increasing levels of pleiotropy over evolutionary time, which accompanies stronger selective constraints. We propose a "pleiotropy-barrier" model that delineates different potentials for phenotypic innovation between young and older genes subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

Effect of recombination on genetic diversity of Caenorhabditis elegans

Greater molecular divergence and genetic diversity are present in regions of high recombination in many species. Studies describing the correlation between variant abundance and recombination rate have long focused on recombination in the context of linked selection models, whereby interference between linked sites under positive or negative selection reduces genetic diversity in regions of low recombination. Here, we show that indels, especially those of intermediate sizes, are enriched relative to single nucleotide polymorphisms in regions of high recombination in C. elegans. To explain this phenomenon, we reintroduce an alternative model that emphasizes the mutagenic effect of recombination. To extend the analysis, we examine the variants with a phylogenetic context and discuss how different models could be examined together. The number of variants generated by recombination in natural populations could be substantial including possibly the majority of some indel subtypes. Our work highlights the potential importance of a mutagenic effect of recombination, which could have a significant role in the shaping of natural genetic diversity.

An empirical evaluation of the estimation of inbreeding depression from molecular markers under suboptimal conditions

Inbreeding depression (ID), the reduction in fitness due to inbreeding, is typically measured by the regression of the phenotypic values of individuals for a particular trait on their corresponding inbreeding coefficients (F). While genealogical records can provide these coefficients, they may be unavailable or incomplete, making molecular markers a useful alternative. The power to detect ID and its accuracy depend on the variation of F values of individuals, the sample sizes available, and the accuracy in the estimation of individual fitness traits and F values. In this study, we used Drosophila melanogaster to evaluate the effectiveness of molecular markers in estimating ID under suboptimal conditions. We generated two sets of 100 pairs of unrelated individuals from a large panmictic population and mated them for two generations to produce non-inbred and unrelated individuals (F = 0) and inbred individuals (full-sib progeny; F = 0.25). Using these expected genealogical F values, we calculated inbreeding depression for two fitness-related traits, pupae productivity and competitive fitness. We then sequenced the males from 17 non-inbred pairs and 17 inbred pairs to obtain their genomic inbreeding coefficients and estimate ID for the two traits. The scenario assumed was rather restrictive in terms of estimation of ID because: (1) the individuals belonged to the same generation of a large panmictic population, leading to low variation in individual F coefficients; (2) the sample sizes were small; and (3) the traits measured depended on both males and females while only males were sequenced. Despite the challenging conditions of our study, we found that molecular markers provided estimates of ID that were comparable to those obtained from simple pedigree estimations with larger sample sizes. The results therefore suggest that genomic measures of inbreeding are useful to provide estimates of inbreeding depression even under very challenging scenarios.

Quantifying the fraction of new mutations that are recessive lethal

The presence and impact of recessive lethal mutations have been widely documented in diploid outcrossing species. However, precise estimates of the proportion of new mutations that are recessive lethal remain limited. Here, we evaluate the performance of Fit∂a∂i, a commonly used method for inferring the distribution of fitness effects (DFE), in the presence of lethal mutations. Using simulations, we demonstrate that in both additive and recessive cases, inference of the deleterious nonlethal portion of the DFE is minimally affected by a small proportion (<10%) of lethal mutations. Additionally, we demonstrate that while Fit∂a∂i cannot estimate the fraction of recessive lethal mutations, Fit∂a∂i can accurately infer the fraction of additive lethal mutations. Finally, as an alternative approach to estimate the proportion of mutations that are recessive lethal, we employ models of mutation-selection-drift balance using existing genomic parameters and estimates of segregating recessive lethals for humans and Drosophila melanogaster. In both species, the segregating recessive lethal load can be explained by a very small fraction (<1%) of new nonsynonymous mutations being recessive lethal. Our results refute recent assertions of a much higher proportion of mutations being recessive lethal (4%-5%), while highlighting the need for additional information on the joint distribution of selection and dominance coefficients.

Mechanisms of robustness in gene regulatory networks involved in neural development

The functions of living organisms are affected by different kinds of perturbation, both internal and external, which in many cases have functional effects and phenotypic impact. The effects of these perturbations become particularly relevant for multicellular organisms with complex body patterns and cell type heterogeneity, where transcriptional programs controlled by gene regulatory networks determine, for example, the cell fate during embryonic development. Therefore, an essential aspect of development in these organisms is the ability to maintain the functionality of their genetic developmental programs even in the presence of genetic variation, changing environmental conditions and biochemical noise, a property commonly termed robustness. We discuss the implication of different molecular mechanisms of robustness involved in neurodevelopment, which is characterized by the interplay of many developmental programs at a molecular, cellular and systemic level. We specifically focus on processes affecting the function of gene regulatory networks, encompassing transcriptional regulatory elements and post-transcriptional processes such as miRNA-based regulation, but also higher order regulatory organization, such as gene network topology. We also present cases where impairment of robustness mechanisms can be associated with neurodevelopmental disorders, as well as reasons why understanding these mechanisms should represent an important part of the study of gene regulatory networks driving neural development.

Genomic Underpinnings of Population Persistence in Isle Royale Moose

Island ecosystems provide natural laboratories to assess the impacts of isolation on population persistence. However, most studies of persistence have focused on a single species, without comparisons to other organisms they interact with in the ecosystem. The case study of moose and gray wolves on Isle Royale allows for a direct contrast of genetic variation in isolated populations that have experienced dramatically differing population trajectories over the past decade. Whereas the Isle Royale wolf population recently declined nearly to extinction due to severe inbreeding depression, the moose population has thrived and continues to persist, despite having low genetic diversity and being isolated for ∼120 years. Here, we examine the patterns of genomic variation underlying the continued persistence of the Isle Royale moose population. We document high levels of inbreeding in the population, roughly as high as the wolf population at the time of its decline. However, inbreeding in the moose population manifests in the form of intermediate-length runs of homozygosity suggestive of historical inbreeding and purging, contrasting with the long runs of homozygosity observed in the smaller wolf population. Using simulations, we confirm that substantial purging has likely occurred in the moose population. However, we also document notable increases in genetic load, which could eventually threaten population viability over the long term. Overall, our results demonstrate a complex relationship between inbreeding, genetic diversity, and population viability that highlights the use of genomic datasets and computational simulation tools for understanding the factors enabling persistence in isolated populations.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

The evolution of personality disorders: A review of proposals

Personality disorders (PDs) are currently considered dysfunctions. However, personality differences are older than humanity and are ubiquitous in nature, from insects to higher primates. This suggests that a number of evolutionary mechanisms-other than dysfunctions-may be able to maintain stable behavioral variation in the gene pool. First of all, apparently maladaptive traits may actually improve fitness by enabling better survival or successful mating or reproduction, as exemplified by neuroticism, psychopathy, and narcissism. Furthermore, some PDs may harm important biological goals while facilitating others, or may be globally beneficial or detrimental depending on environmental circumstances or body condition. Alternatively, certain traits may form part of life history strategies: Coordinated suites of morphological, physiological and behavioral characters that optimize fitness through alternative routes and respond to selection as a whole. Still others may be vestigial adaptations that are no longer beneficial in present times. Finally, variation may be adaptative in and by itself, as it reduces competition for finite resources. These and other evolutionary mechanisms are reviewed and illustrated through human and non-human examples. Evolutionary theory is the best-substantiated explanatory framework across the life sciences, and may shed light on the question of why harmful personalities exist at all.

Sexual selection for males with beneficial mutations

Sexual selection is the process by which traits providing a mating advantage are favoured. Theoretical treatments of the evolution of sex by sexual selection propose that it operates by reducing the load of deleterious mutations. Here, we postulate instead that sexual selection primarily acts through females preferentially mating with males carrying beneficial mutations. We used simulation and analytical modelling to investigate the evolutionary dynamics of beneficial mutations in the presence of sexual selection. We found that female choice for males with beneficial mutations had a much greater impact on genetic quality than choice for males with low mutational load. We also relaxed the typical assumption of a fixed mutation rate. For deleterious mutations, mutation rate should always be minimized, but when rare beneficial mutations can occur, female choice for males with those rare beneficial mutations could overcome a decline in average fitness and allow an increase in mutation rate. We propose that sexual selection for beneficial mutations could overcome the ‘two-fold cost of sex’ much more readily than choice for males with low mutational load and may therefore be a more powerful explanation for the prevalence of sexual reproduction than the existing theory. If sexual selection results in higher fitness at higher mutation rates, and if the variability produced by mutation itself promotes sexual selection, then a feedback loop between these two factors could have had a decisive role in driving adaptation.

How much does the unguarded X contribute to sex differences in life span?

Females and males often have markedly different mortality rates and life spans, but it is unclear why these forms of sexual dimorphism evolve. The unguarded X hypothesis contends that dimorphic life spans arise from sex differences in X or Z chromosome copy number (i.e., one copy in the “heterogametic” sex; two copies in the “homogametic” sex), which leads to a disproportionate expression of deleterious mutations by the heterogametic sex (e.g., mammalian males; avian females). Although data on adult sex ratios and sex‐specific longevity are consistent with predictions of the unguarded X hypothesis, direct experimental evidence remains scant, and alternative explanations are difficult to rule out. Using a simple population genetic model, we show that the unguarded X effect on sex differential mortality is a function of several reasonably well‐studied evolutionary parameters, including the proportion of the genome that is sex linked, the genomic deleterious mutation rate, the mean dominance of deleterious mutations, the relative rates of mutation and strengths of selection in each sex, and the average effect of mutations on survival and longevity relative to their effects on fitness. We review published estimates of these parameters, parameterize our model with them, and show that unguarded X effects are too small to explain observed sex differences in life span across species. For example, sex differences in mean life span are known to often exceed 20% (e.g., in mammals), whereas our parameterized models predict unguarded X effects of a few percent (e.g., 1–3% in Drosophila and mammals). Indeed, these predicted unguarded X effects fall below statistical thresholds of detectability in most experiments, potentially explaining why direct tests of the hypothesis have generated little support for it. Our results suggest that evolution of sexually dimorphic life spans is predominantly attributable to other mechanisms, potentially including “toxic Y” effects and sexual dimorphism for optimal investment in survival versus reproduction.

Females and males are dimorphic for a wide range of traits, including the average lengths of their life spans. Sex differences in life span are both conspicuous and variable among species. For example, in mammals, females live ∼20% longer than males (on average), whereas in birds, males live ∼10% longer than females. One leading explanation for these patterns—the unguarded X hypothesis—argues that sex differences in life span emerge from the distinct sex chromosomes that females and males inherit. For many species, one sex (e.g., female mammals; male birds) carries two copies of each X‐linked gene, whereas the other carries one. Because harmful mutations are partially recessive, the sex with only one copy of the X is more prone to expressing them, and that sex should therefore have a shorter average life span. This prediction of the unguarded X hypothesis is qualitatively consistent with observations of sex‐ratio bias in adults and sexual dimorphism for longevity (e.g., mammalian males have one copy of the X and have shorter lives than females). However, there are other possible explanations for these patterns, making it unclear how much the unguarded X explains species diversity for sex‐specific longevity. We developed a mathematical model for the contribution of unguarded X effects to sex differences in survival and life span, and used data on mutation rates and their effects on survival and fitness to quantify the importance of the unguarded X across species. The model, when combined with current data, suggests that the unguarded X hypothesis cannot explain the conspicuous sex differences in life span that are commonly reported in animal species, particularly vertebrates. Our results suggest that the unguarded X is an unlikely general explanation for the evolution of sexually dimorphic life spans, which gives weight to alternative mechanisms, including “toxic Y” effects and sex differential selection via trade‐offs between survival and reproduction.

Cancer Risk and Mutational Patterns Following Organ Transplantation

The rapid development of medical technology and widespread application of immunosuppressive drugs have improved the success rate of organ transplantation significantly. However, the use of immunosuppressive agents increases the frequency of malignancy greatly. With the prospect of “precision medicine” for tumors and development of next-generation sequencing technology, more attention has been paid to the application of high-throughput sequencing technology in clinical oncology research, which is mainly applied to the early diagnosis of tumors and analysis of tumor-related genes. All generations of cancers carry somatic mutations, meanwhile, significant differences were observed in mutational signatures across tumors. Systematic sequencing of cancer genomes from patients after organ transplantation can reveal DNA damage and repair processes in exposed cancer cells and their precursors. In this review, we summarize the application of high-throughput sequencing and organoids in the field of organ transplantation, the mutational patterns of cancer genomes, and propose a new research strategy for understanding the mechanism of cancer following organ transplantation.

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.

Recommendations for improving statistical inference in population genomics

The field of population genomics has grown rapidly in response to the recent advent of affordable, large-scale sequencing technologies. As opposed to the situation during the majority of the 20th century, in which the development of theoretical and statistical population genetic insights outpaced the generation of data to which they could be applied, genomic data are now being produced at a far greater rate than they can be meaningfully analyzed and interpreted. With this wealth of data has come a tendency to focus on fitting specific (and often rather idiosyncratic) models to data, at the expense of a careful exploration of the range of possible underlying evolutionary processes. For example, the approach of directly investigating models of adaptive evolution in each newly sequenced population or species often neglects the fact that a thorough characterization of ubiquitous nonadaptive processes is a prerequisite for accurate inference. We here describe the perils of these tendencies, present our consensus views on current best practices in population genomic data analysis, and highlight areas of statistical inference and theory that are in need of further attention. Thereby, we argue for the importance of defining a biologically relevant baseline model tuned to the details of each new analysis, of skepticism and scrutiny in interpreting model fitting results, and of carefully defining addressable hypotheses and underlying uncertainties.

The relative fitness of the de novo variants in general Lithuanian population vs. in individuals with intellectual disability

The effect of a variant on an organism is always multifaceted and can be considered from multiple perspectives-biochemical, medical, or evolutionary. However, the relationship between the effects of amino acid substitution on protein activity, human health, and an individual's evolutionary fitness is not trivial. We uncover that the general Lithuanian population is characterized by a "mirror reflection" of the de novo variant fitness effect, confirming the theory of neutrality. Meanwhile, in the group of individuals with intellectual disability, compared with the reference exome de novo variants significantly changed the composition of the amino acid. Therefore, it predicts that, both in terms of the number of amino acids and changes in their relative fitness, the structure of the proteins encoded by the studied amino acids undergo significant changes following the de novo variant, leading to possible changes in protein function associated with phenotypic traits. These results suggest that the analysis of relative fitness of exome sequences with de novo variants can predict the future phenotype. Therefore even in those cases, then only a few of all functional prediction analysis tools predict a variant as damaging, the negative relative fitness or even adaptability of the genome variant should be carefully evaluated considering both its direct function and the global background of the possible disease-associated mechanism regardless of the phenotype being studied.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10-10, and is highly variable between MA lines, ranging from μ = 0.35 × 10-10 to μ = 131.7 × 10-10. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.

A simple expression for the strength of selection on recombination generated by interference among mutations

One of the most widely cited hypotheses to explain the evolutionary maintenance of genetic recombination states that the reshuffling of genotypes at meiosis increases the efficiency of natural selection by reducing interference among selected loci. However, and despite several decades of theoretical work, a quantitative estimation of the possible selective advantage of a mutant allele increasing chromosomal map length (the average number of cross-overs at meiosis) remains difficult. This article derives a simple expression for the strength of selection acting on a modifier gene affecting the genetic map length of a whole chromosome or genome undergoing recurrent mutation. In particular, it shows that indirect selection for recombination caused by interference among mutations is proportional to [Formula: see text], where [Formula: see text] is the effective population size, U is the deleterious mutation rate per chromosome, and R is the chromosome map length. Indirect selection is relatively insensitive to the fitness effects of deleterious alleles, epistasis, or the genetic architecture of recombination rate variation and may compensate for substantial costs associated with recombination when linkage is tight. However, its effect generally stays weak in large, highly recombining populations.

Population Modelling in Affective Disorders

Purpose of Review

The prevalence of affective disorders is on the rise. This upward trajectory leads to a substantial personal and societal cost. There is growing body of literature demonstrating decision-making impairments associated with affective disorders, and more studies are using computational modelling methods to infer underlying mechanisms of these impairments from participant choice behaviour. However, lack of population modelling suggests that data resources may still be underutilised.

Recent Findings

A number of recent studies associated major depression with abnormal risky decision-making as well as impairments in temporal discounting and social decision-making. These domains capture relevant aspects of real-life decision-making. Consequently, data from these studies can be used to define behavioural phenotypes for major depression.

Summary

The manuscript describes a detailed proposal for population modelling to capture changes in the prevalence rate of major depression. The population modelling approach can also identify which decision-making domains can account for a larger part of impairments in psychosocial functioning and how behavioural interventions built on computational principles can target these to improve real-life psychosocial functioning in patient groups.

Origin of Sex-Biased Mental Disorders: An Evolutionary Perspective

Sexual dimorphism or sex bias in diseases and mental disorders have two biological causes: sexual selection and sex hormones. We review the role of sexual selection theory and bring together decades of molecular studies on the variation and evolution of sex-biased genes and provide a theoretical basis for the causes of sex bias in disease and health. We present a Sexual Selection-Sex Hormone theory and show that male-driven evolution, including sexual selection, leads to: (1) increased male vulnerability due to negative pleiotropic effects associated with male-driven sexual selection and evolution; (2) increased rates of male-driven mutations and epimutations in response to early fitness gains and at the cost of late fitness; and (3) enhanced female immunity due to antagonistic responses to mutations that are beneficial to males but harmful to females, reducing female vulnerability to diseases and increasing the thresholds for disorders such as autism. Female-driven evolution, such as reproduction-related fluctuation in female sex hormones in association with stress and social condition, has been shown to be associated with increased risk of certain mental disorders such as major depression disorder in women. Bodies have history, cells have memories. An evolutionary framework, such as the Sexual Selection–Sex Hormone theory, provides a historical perspective for understanding how the differences in the sex-biased diseases and mental disorders have evolved over time. It has the potential to direct the development of novel preventive and treatment strategies.

Sites of chromosomal instability in the context of nuclear architecture and function

Chromosomal fragile sites are described as areas within the tightly packed mitotic chromatin that appear as breaks or gaps mostly tracing back to a loosened structure and not a real nicked break within the DNA molecule. Most facts about fragile sites result from studies in mitotic cells, mainly during metaphase and mainly in lymphocytes. Here, we synthesize facts about the genomic regions that are prone to form gaps and breaks on metaphase chromosomes in the context of interphase. We conclude that nuclear architecture shapes the activity profile of the cell, i.e. replication timing and transcriptional activity, thereby influencing genomic integrity during interphase with the potential to cause fragility in mitosis. We further propose fragile sites as examples of regions specifically positioned in the interphase nucleus with putative anchoring points at the nuclear lamina to enable a tightly regulated replication–transcription profile and diverse signalling functions in the cell. Consequently, fragility starts before the actual display as chromosomal breakage in metaphase to balance the initial contradiction of cellular overgrowth or malfunctioning and maintaining diversity in molecular evolution.

CpG-creating mutations are costly in many human viruses

Mutations can occur throughout the virus genome and may be beneficial, neutral or deleterious. We are interested in mutations that yield a C next to a G, producing CpG sites. CpG sites are rare in eukaryotic and viral genomes. For the eukaryotes, it is thought that CpG sites are rare because they are prone to mutation when methylated. In viruses, we know less about why CpG sites are rare. A previous study in HIV suggested that CpG-creating transition mutations are more costly than similar non-CpG-creating mutations. To determine if this is the case in other viruses, we analyzed the allele frequencies of CpG-creating and non-CpG-creating mutations across various strains, subtypes, and genes of viruses using existing data obtained from Genbank, HIV Databases, and Virus Pathogen Resource. Our results suggest that CpG sites are indeed costly for most viruses. By understanding the cost of CpG sites, we can obtain further insights into the evolution and adaptation of viruses.

Investigation of de novo variation in pediatric cardiomyopathy

Pediatric cardiomyopathies can be caused by variants in genes encoding the sarcomere and cytoskeleton in cardiomyocytes. Variants are typically inherited in an autosomal dominant manner with variable expressivity. De novo variants have been reported, however their overall frequency is largely unknown. We sought to determine the rate of de novo, pathogenic and likely pathogenic (P/LP) variants in children with a diagnosis of hypertrophic, dilated, or restrictive cardiomyopathy (HCM, DCM, or RCM), and to compare disease outcomes between individuals with and without a de novo variant. A retrospective record review identified 126 individuals with HCM (55%), DCM (37%), or RCM (8%) ≤18 years of age who had genetic testing. Overall, 50 (40%) had positive genetic testing and 18% of P/LP variants occurred de novo. The rate of de novo variation in those with RCM (80%) was higher than in those with HCM (9%) or DCM (20%). There was evidence of germline mosaicism in one family with RCM. Individuals with de novo variants were more likely than those without to have a history of arrhythmia (p = .049), sudden cardiac arrest (p = .024), hospitalization (p = .041), and cardiac transplantation (p = .030). The likelihood of de novo variation and impact on family risk and screening should be integrated into genetic counseling.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.

Base-excision repair deficiency alone or combined with increased oxidative stress does not increase mtDNA point mutations in mice

Mitochondrial DNA (mtDNA) mutations become more prevalent with age and are postulated to contribute to the ageing process. Point mutations of mtDNA have been suggested to originate from two main sources, i.e. replicative errors and oxidative damage, but the contribution of each of these processes is much discussed. To elucidate the origin of mtDNA mutations, we measured point mutation load in mice with deficient mitochondrial base-excision repair (BER) caused by knockout alleles preventing mitochondrial import of the DNA repair glycosylases OGG1 and MUTYH (Ogg1 dMTS, Mutyh dMTS). Surprisingly, we detected no increase in the mtDNA mutation load in old Ogg1 dMTS mice. As DNA repair is especially important in the germ line, we bred the BER deficient mice for five consecutive generations but found no increase in the mtDNA mutation load in these maternal lineages. To increase reactive oxygen species (ROS) levels and oxidative damage, we bred the Ogg1 dMTS mice with tissue specific Sod2 knockout mice. Although increased superoxide levels caused a plethora of changes in mitochondrial function, we did not detect any changes in the mutation load of mtDNA or mtRNA. Our results show that the importance of oxidative damage as a contributor of mtDNA mutations should be re-evaluated.

Somatic Mutagenesis in Mammals and Its Implications for Human Disease and Aging

DNA mutations as a consequence of errors during DNA damage repair, replication, or mitosis are the substrate for evolution. In multicellular organisms, mutations can occur in the germline and also in somatic tissues, where they are associated with cancer and other chronic diseases and possibly with aging. Recent advances in high-throughput sequencing have made it relatively easy to study germline de novo mutations, but in somatic cells, the vast majority of mutations are low-abundant and can be detected only in clonal lineages, such as tumors, or single cells. Here we review recent results on somatic mutations in normal human and animal tissues with a focus on their possible functional consequences.

Parental age in relation to the severity of cleft lip and/or palate

OBJECTIVE

To investigate parental age in relation to the severity of cleft diagnosis in a population-based consecutive sample of individuals with clefts.

SETTING/SAMPLE

Retrospective, consecutive.

MATERIAL/METHODS

The sample comprised 545 consecutive cases with nonsyndromic clefts (437 individuals with cleft lip with/without cleft palate [CL ± P] and 106 individuals with isolated cleft palate [CP]) and parental ages recorded. The groups were subdivided according to the extent/severity of cleft. Unilateral clefts were divided according to left/right sidedness. Logistic regression was used to analyse the association between parental age and risk of severe cleft separately for CL ± P and CP, as well as between parental age and risk of right-sided cleft.

RESULTS

In CL ± P, the risk of a complete cleft in the offspring increases with maternal age when the paternal age exceeds approximately 29 years. Moreover, the risk is higher when both parents are old than when both parents are young. In CP, no statistically significant results were identified. However, there were clear trends that indicated a similar pattern as that for CL ± P. No association was identified between increased parental age and the sidedness of clefts.

CONCLUSIONS

Parental age seems to contribute to cleft severity, as older parents showed a higher risk of having offspring with a complete cleft than younger parents.

Role of recombination and faithfulness to partner in sex chromosome degeneration

Sex determination in mammals is strongly linked to sex chromosomes. In most cases, females possess two copies of X chromosome while males have one X and one Y chromosome. It is assumed that these chromosomes originated from a pair of homologous autosomes, which diverged when recombination between them was suppressed. However, it is still debated why the sex chromosomes stopped recombining and how this process spread out over most part of the chromosomes. To study this problem, we developed a simulation model, in which the recombination rate between the sex chromosomes can freely evolve. We found that the suppression of recombination between the X and Y is spontaneous and proceeds very quickly during the evolution of population, which leads to the degeneration of the Y in males. Interestingly, the degeneration happens only when mating pairs are unfaithful. This evolutionary strategy purifies the X chromosome from defective alleles and leads to the larger number of females than males in the population. In consequence, the reproductive potential of the whole population increases. Our results imply that both the suppression of recombination and the degeneration of Y chromosome may be associated with reproductive strategy and favoured in polygamous populations with faithless mating partners.

An Upper Limit on the Functional Fraction of the Human Genome

For the human population to maintain a constant size from generation to generation, an increase in fertility must compensate for the reduction in the mean fitness of the population caused, among others, by deleterious mutations. The required increase in fertility due to this mutational load depends on the number of sites in the genome that are functional, the mutation rate, and the fraction of deleterious mutations among all mutations in functional regions. These dependencies and the fact that there exists a maximum tolerable replacement level fertility can be used to put an upper limit on the fraction of the human genome that can be functional. Mutational load considerations lead to the conclusion that the functional fraction within the human genome cannot exceed 15%.

Direct estimation of de novo mutation rates in a chimpanzee parent-offspring trio by ultra-deep whole genome sequencing

Mutations generate genetic variation and are a major driving force of evolution. Therefore, examining mutation rates and modes are essential for understanding the genetic basis of the physiology and evolution of organisms. Here, we aim to identify germline de novo mutations through the whole-genome surveyance of Mendelian inheritance error sites (MIEs), those not inherited through the Mendelian inheritance manner from either of the parents, using ultra-deep whole genome sequences (>150-fold) from a chimpanzee parent-offspring trio. We identified such 889 MIEs and classified them into four categories based on the pattern of inheritance and the sequence read depth: [i] de novo single nucleotide variants (SNVs), [ii] copy number neutral inherited variants, [iii] hemizygous deletion inherited variants, and [iv] de novo copy number variants (CNVs). From de novo SNV candidates, we estimated a germline de novo SNV mutation rate as 1.48 × 10^-8 per site per generation or 0.62 × 10^-9 per site per year. In summary, this study demonstrates the significance of ultra-deep whole genome sequencing not only for the direct estimation of mutation rates but also for discerning various mutation modes including de novo allelic conversion and de novo CNVs by identifying MIEs through the transmission of genomes from parents to offspring.

Older fathers' children have lower evolutionary fitness across four centuries and in four populations

Higher paternal age at offspring conception increases de novo genetic mutations. Based on evolutionary genetic theory we predicted older fathers' children, all else equal, would be less likely to survive and reproduce, i.e. have lower fitness. In sibling control studies, we find support for negative paternal age effects on offspring survival and reproductive success across four large populations with an aggregate N > 1.4 million. Three populations were pre-industrial (1670-1850) Western populations and showed negative paternal age effects on infant survival and offspring reproductive success. In twentieth-century Sweden, we found minuscule paternal age effects on survival, but found negative effects on reproductive success. Effects survived tests for key competing explanations, including maternal age and parental loss, but effects varied widely over different plausible model specifications and some competing explanations such as diminishing paternal investment and epigenetic mutations could not be tested. We can use our findings to aid in predicting the effect increasingly older parents in today's society will have on their children's survival and reproductive success. To the extent that we succeeded in isolating a mutation-driven effect of paternal age, our results can be understood to show that de novo mutations reduce offspring fitness across populations and time periods.

Older fathers' children have lower evolutionary fitness across four centuries and in four populations

Higher paternal age at offspring conception increases de novo genetic mutations. Based on evolutionary genetic theory we predicted older fathers' children, all else equal, would be less likely to survive and reproduce, i.e. have lower fitness. In sibling control studies, we find support for negative paternal age effects on offspring survival and reproductive success across four large populations with an aggregate N > 1.4 million. Three populations were pre-industrial (1670-1850) Western populations and showed negative paternal age effects on infant survival and offspring reproductive success. In twentieth-century Sweden, we found minuscule paternal age effects on survival, but found negative effects on reproductive success. Effects survived tests for key competing explanations, including maternal age and parental loss, but effects varied widely over different plausible model specifications and some competing explanations such as diminishing paternal investment and epigenetic mutations could not be tested. We can use our findings to aid in predicting the effect increasingly older parents in today's society will have on their children's survival and reproductive success. To the extent that we succeeded in isolating a mutation-driven effect of paternal age, our results can be understood to show that de novo mutations reduce offspring fitness across populations and time periods.

Genome instability: a conserved mechanism of ageing?

DNA is the carrier of genetic information and the primary template from which all cellular information is ultimately derived. Changes in the DNA information content through mutation generate diversity for evolution through natural selection but are also a source of deleterious effects. It has since long been hypothesized that mutation accumulation in somatic cells of multicellular organisms could causally contribute to age-related cellular degeneration and death. Assays to detect different types of mutations, from base substitutions to large chromosomal aberrations, have been developed and show unequivocally that mutations accumulate in different tissues and cell types of ageing humans and animals. More recently, next-generation sequencing-based methods have been developed to accurately determine the complete landscape of base substitution mutations in single cells. The first results show that the somatic mutation rate is much higher than the germline mutation rate and that base substitution loads in somatic cells are high enough to potentially affect cellular function.

Evaluation of chemical mutagenicity using next generation sequencing: A review

Mutations are heritable changes in the nucleotide sequence of DNA that can lead to many adverse effects. Genotoxicity assays have been used to identify chemical mutagenicity. Recently, next generation sequencing (NGS) has been used for this purpose. In this review, we present the progress in NGS application for assessing mutagenicity of chemicals, including the methods used for detecting the induced mutations, bioinformatics tools for analyzing the sequencing data, and chemicals whose mutagenicity has been evaluated using NGS. Available information suggests that NGS technology has unparalleled advantages for evaluating mutagenicity of chemicals can be applied for the next generation of mutagenicity tests.

A high-fidelity method for genomic sequencing of single somatic cells reveals a very high mutational burden

Postzygotic mutations in somatic cells lead to genome mosaicism and can be the cause of cancer, possibly other human diseases and aging. Somatic mutations are difficult to detect in bulk tissue samples. Here, we review the available assays for measuring somatic mutations, with a focus on recent single-cell, whole genome sequencing methods. Impact statement Somatic mutations cause cancer, possibly other diseases and aging. Yet, very little is known about the frequency of such mutations in vivo, their distribution across the genome, and their possible functional consequences other than cancer. Even in cancer, we do not know the heterogeneity of mutations within a tumor and if seemingly normal cells in its surroundings already have elevated mutation frequencies. Here, we review a new, whole genome amplification system that allows accurate quantification and characterization of single-cell mutational landscapes in human cells and tissues in relation to disease.

Genetic drift, selection and the evolution of the mutation rate

As one of the few cellular traits that can be quantified across the tree of life, DNA-replication fidelity provides an excellent platform for understanding fundamental evolutionary processes. Furthermore, because mutation is the ultimate source of all genetic variation, clarifying why mutation rates vary is crucial for understanding all areas of biology. A potentially revealing hypothesis for mutation-rate evolution is that natural selection primarily operates to improve replication fidelity, with the ultimate limits to what can be achieved set by the power of random genetic drift. This drift-barrier hypothesis is consistent with comparative measures of mutation rates, provides a simple explanation for the existence of error-prone polymerases and yields a formal counter-argument to the view that selection fine-tunes gene-specific mutation rates.

Hill-Robertson interference maintained by Red Queen dynamics favours the evolution of sex

Although it is well established theoretically that selective interference among mutations (Hill-Robertson interference) favours meiotic recombination, genomewide mean rates of mutation and strengths of selection appear too low to support this as the mechanism favouring recombination in nature. A possible solution to this discrepancy between theory and observation is that selection is at least intermittently very strong due to the antagonistic coevolution between a host and its parasites. The Red Queen theory posits that such coevolution generates fitness epistasis among loci, which generates negative linkage disequilibrium among beneficial mutations, which in turn favours recombination. This theory has received only limited support. However, Red Queen dynamics without epistasis may provide the ecological conditions that maintain strong and frequent selective interference in finite populations that indirectly selects for recombination. This hypothesis is developed here through the simulation of Red Queen dynamics. This approach required the development of a method to calculate the exact frequencies of multilocus haplotypes after recombination. Simulations show that recombination is favoured by the moderately weak selection of many loci involved in the interaction between a host and its parasites, which results in substitution rates that are compatible with empirical estimates. The model also reproduces the previously reported rapid increase in the rate of outcrossing in Caenorhabditis elegans coevolving with a bacterial pathogen.