Interpreting the pervasive observation of U-shaped Site Frequency Spectra.

The standard neutral model of molecular evolution has traditionally been used as the null model for population genomics. We gathered a collection of 45 genome-wide site frequency spectra from a diverse set of species, most of which display an excess of low and high frequency variants compared to the expectation of the standard neutral model, resulting in U-shaped spectra. We show that multiple merger coalescent models often provide a better fit to these observations than the standard Kingman coalescent. Hence, in many circumstances these under-utilized models may serve as the more appropriate reference for genomic analyses. We further discuss the underlying evolutionary processes that may result in the widespread U-shape of frequency spectra.

CNCA aligns small annotated genomes.

BACKGROUND: To explore the evolutionary history of sequences, a sequence alignment is a first and necessary step, and its quality is crucial. In the context of the study of the proximal origins of SARS-CoV-2 coronavirus, we wanted to construct an alignment of genomes closely related to SARS-CoV-2 using both coding and non-coding sequences. To our knowledge, there is no tool that can be used to construct this type of alignment, which motivated the creation of CNCA. RESULTS: CNCA is a web tool that aligns annotated genomes from GenBank files. It generates a nucleotide alignment that is then updated based on the protein sequence alignment. The output final nucleotide alignment matches the protein alignment and guarantees no frameshift. CNCA was designed to align closely related small genome sequences up to 50 kb (typically viruses) for which the gene order is conserved. CONCLUSIONS: CNCA constructs multiple alignments of small genomes by integrating both coding and non-coding sequences. This preserves regions traditionally ignored in conventional back-translation methods, such as non-coding regions.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes.

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

A Minimal yet Flexible Likelihood Framework to Assess Correlated Evolution.

An evolutionary process is reflected in the sequence of changes of any trait (e.g., morphological or molecular) through time. Yet, a better understanding of evolution would be procured by characterizing correlated evolution, or when two or more evolutionary processes interact. Previously developed parametric methods often require significant computing time as they rely on the estimation of many parameters. Here, we propose a minimal likelihood framework modeling the joint evolution of two traits on a known phylogenetic tree. The type and strength of correlated evolution are characterized by a few parameters tuning mutation rates of each trait and interdependencies between these rates. The framework can be applied to study any discrete trait or character ranging from nucleotide substitution to gain or loss of a biological function. More specifically, it can be used to 1) test for independence between two evolutionary processes, 2) identify the type of interaction between them, and 3) estimate parameter values of the most likely model of interaction. In the current implementation, the method takes as input a phylogenetic tree with discrete evolutionary events mapped on its branches. The method then maximizes the likelihood for one or several chosen scenarios. The strengths and limits of the method, as well as its relative power compared to a few other methods, are assessed using both simulations and data from 16S rRNA sequences in a sample of 54 $\gamma$-enterobacteria. We show that, even with data sets of fewer than 100 species, the method performs well in parameter estimation and in evolutionary model selection. [Correlated evolution; maximum likelihood; model.].

Testing for population decline using maximal linkage disequilibrium blocks.

Only 6% of known species have a conservation status. Methods that assess conservation statuses are often based on individual counts and are thus too laborious to be generalized to all species. Population genomics methods that infer past variations in population size are easy to use but limited to the relatively distant past. Here we propose a population genomics approach that tests for recent population decline and may be used to assess species conservation statuses. More specifically, we study Maximal Recombination Free (MRF) blocks, that are segments of a sequence alignment inherited from a common ancestor without recombination. MRF blocks are relatively longer in small than in large populations. We use the distribution of MRF block lengths rescaled by their mean to test for recent population decline. However, because MRF blocks are difficult to detect, we also consider Maximal Linkage Disequilibrium (MLD) blocks, which are runs of single nucleotide polymorphisms compatible with a single tree. We develop a new method capable of inferring a very recent decline (e.g. with a detection power of 50% for populations whose size was halved to N, 0.05 ×N generations ago) from rescaled MLD block lengths. Our framework could serve as a basis for quantitative tools to assess conservation status in a wide range of species.

The genomic view of diversification.

The process of species diversification is traditionally summarized by a single tree, the species tree, whose reconstruction from molecular data is hindered by frequent conflicts between gene genealogies. Here, we argue that instead of seeing these conflicts as nuisances, we can exploit them to inform the diversification process itself. We adopt a gene-based view of diversification to model the ubiquitous presence of gene flow between diverging lineages, one of the most important processes explaining disagreements among gene trees. We propose a new framework for modelling the joint evolution of gene and species lineages relaxing the hierarchy between the species tree and gene trees inherent to the standard view, as embodied in a popular model known as the multispecies coalescent (MSC). We implement this framework in two alternative models called the gene-based diversification models (GBD): (a) GBD-forward following all evolving genomes through time and (b) GBD-backward based on coalescent theory. They feature four parameters tuning colonization, gene flow, genetic drift and genetic differentiation. We propose an inference method based on differences between gene trees. Applied to two empirical data sets prone to gene flow, we find better support for the GBD-backward model than for the MSC model. Along with the increasing awareness of the extent of gene flow, this work shows the importance of considering the richer signal contained in genomic histories, rather than in the mere species tree, to better apprehend the complex evolutionary history of species.

A Molecular Portrait of De Novo Genes in Yeasts.

New genes, with novel protein functions, can evolve "from scratch" out of intergenic sequences. These de novo genes can integrate the cell's genetic network and drive important phenotypic innovations. Therefore, identifying de novo genes and understanding how the transition from noncoding to coding occurs are key problems in evolutionary biology. However, identifying de novo genes is a difficult task, hampered by the presence of remote homologs, fast evolving sequences and erroneously annotated protein coding genes. To overcome these limitations, we developed a procedure that handles the usual pitfalls in de novo gene identification and predicted the emergence of 703 de novo gene candidates in 15 yeast species from 2 genera whose phylogeny spans at least 100 million years of evolution. We validated 85 candidates by proteomic data, providing new translation evidence for 25 of them through mass spectrometry experiments. We also unambiguously identified the mutations that enabled the transition from noncoding to coding for 30 Saccharomyces de novo genes. We established that de novo gene origination is a widespread phenomenon in yeasts, only a few being ultimately maintained by selection. We also found that de novo genes preferentially emerge next to divergent promoters in GC-rich intergenic regions where the probability of finding a fortuitous and transcribed ORF is the highest. Finally, we found a more than 3-fold enrichment of de novo genes at recombination hot spots, which are GC-rich and nucleosome-free regions, suggesting that meiotic recombination contributes to de novo gene emergence in yeasts.

Clinical and genetic study of Tunisian families with genetic generalized epilepsy: contribution of CACNA1H and MAST4 genes.

Genetic generalized epilepsies (GGE) (childhood absence epilepsy (CAE), juvenile myoclonic epilepsy (JME) and epilepsy with generalized tonic-clonic seizures (GTCS)) are mainly determined by genetic factors. Since few mutations were identified in rare families with autosomal dominant GGE, a polygenic inheritance was suspected in most patients. Recent studies on large American or European cohorts of sporadic cases showed that susceptibility genes were numerous although their variants were rare, making their identification difficult. Here, we reported clinical and genetic characteristics of 30 Tunisian GGE families, including 71 GGE patients. The phenotype was close to that in sporadic cases. Nineteen pedigrees had a homogeneous type of GGE (JME-CAE-CGTS), and 11 combined these epileptic syndromes. Rare non-synonymous variants were selected in probands using a targeted panel of 30 candidate genes and their segregation was determined in families. Molecular studies incriminated different genes, mainly CACNA1H and MAST4. The segregation of at least two variants in different genes in some pedigrees was compatible with the hypothesis of an oligogenic inheritance, which was in accordance with the relatively low frequency of consanguineous probands. Since at least 2 susceptibility genes were likely shared by different populations, genetic factors involved in the majority of Tunisian GGE families remain to be discovered. Their identification should be easier in families with a homogeneous type of GGE, in which an intra-familial genetic homogeneity could be suspected.

The neutral frequency spectrum of linked sites.

We introduce the conditional Site Frequency Spectrum (SFS) for a genomic region linked to a focal mutation of known frequency. An exact expression for its expected value is provided for the neutral model without recombination. Its relation with the expected SFS for two sites, 2-SFS, is discussed. These spectra derive from the coalescent approach of Fu (1995) for finite samples, which is reviewed. Remarkably simple expressions are obtained for the linked SFS of a large population, which are also solutions of the multi-allelic Kolmogorov equations. These formulae are the immediate extensions of the well known single site θ∕f neutral SFS. Besides the general interest in these spectra, they relate to relevant biological cases, such as structural variants and introgressions. As an application, a recipe to adapt Tajima's D and other SFS-based neutrality tests to a non-recombining region containing a neutral marker is presented.

Evolutionary constraints in fitness landscapes.

In the last years, several genotypic fitness landscapes-combinations of a small number of mutations-have been experimentally resolved. To learn about the general properties of "real" fitness landscapes, it is key to characterize these experimental landscapes via simple measures of their structure, related to evolutionary features. Some of the most relevant measures are based on the selectively acessible paths and their properties. In this paper, we present some measures of evolutionary constraints based on (i) the similarity between accessible paths and (ii) the abundance and characteristics of "chains" of obligatory mutations, that are paths going through genotypes with a single fitter neighbor. These measures have a clear evolutionary interpretation. Furthermore, we show that chains are only weakly correlated to classical measures of epistasis. In fact, some of these measures of constraint are non-monotonic in the amount of epistatic interactions, but have instead a maximum for intermediate values. Finally, we show how these measures shed light on evolutionary constraints and predictability in experimentally resolved landscapes.

Mass extinction in poorly known taxa.

Since the 1980s, many have suggested we are in the midst of a massive extinction crisis, yet only 799 (0.04%) of the 1.9 million known recent species are recorded as extinct, questioning the reality of the crisis. This low figure is due to the fact that the status of very few invertebrates, which represent the bulk of biodiversity, have been evaluated. Here we show, based on extrapolation from a random sample of land snail species via two independent approaches, that we may already have lost 7% (130,000 extinctions) of the species on Earth. However, this loss is masked by the emphasis on terrestrial vertebrates, the target of most conservation actions. Projections of species extinction rates are controversial because invertebrates are essentially excluded from these scenarios. Invertebrates can and must be assessed if we are to obtain a more realistic picture of the sixth extinction crisis.

The Impact of Selection, Gene Conversion, and Biased Sampling on the Assessment of Microbial Demography.

Recent studies have linked demographic changes and epidemiological patterns in bacterial populations using coalescent-based approaches. We identified 26 studies using skyline plots and found that 21 inferred overall population expansion. This surprising result led us to analyze the impact of natural selection, recombination (gene conversion), and sampling biases on demographic inference using skyline plots and site frequency spectra (SFS). Forward simulations based on biologically relevant parameters from Escherichia coli populations showed that theoretical arguments on the detrimental impact of recombination and especially natural selection on the reconstructed genealogies cannot be ignored in practice. In fact, both processes systematically lead to spurious interpretations of population expansion in skyline plots (and in SFS for selection). Weak purifying selection, and especially positive selection, had important effects on skyline plots, showing patterns akin to those of population expansions. State-of-the-art techniques to remove recombination further amplified these biases. We simulated three common sampling biases in microbiological research: uniform, clustered, and mixed sampling. Alone, or together with recombination and selection, they further mislead demographic inferences producing almost any possible skyline shape or SFS. Interestingly, sampling sub-populations also affected skyline plots and SFS, because the coalescent rates of populations and their sub-populations had different distributions. This study suggests that extreme caution is needed to infer demographic changes solely based on reconstructed genealogies. We suggest that the development of novel sampling strategies and the joint analyzes of diverse population genetic methods are strictly necessary to estimate demographic changes in populations where selection, recombination, and biased sampling are present.

Testing for Independence between Evolutionary Processes.

Evolutionary events co-occurring along phylogenetic trees usually point to complex adaptive phenomena, sometimes implicating epistasis. While a number of methods have been developed to account for co-occurrence of events on the same internal or external branch of an evolutionary tree, there is a need to account for the larger diversity of possible relative positions of events in a tree. Here we propose a method to quantify to what extent two or more evolutionary events are associated on a phylogenetic tree. The method is applicable to any discrete character, like substitutions within a coding sequence or gains/losses of a biological function. Our method uses a general approach to statistically test for significant associations between events along the tree, which encompasses both events inseparable on the same branch, and events genealogically ordered on different branches. It assumes that the phylogeny and themapping of branches is known without errors. We address this problem from the statistical viewpoint by a linear algebra representation of the localization of the evolutionary events on the tree.We compute the full probability distribution of the number of paired events occurring in the same branch or in different branches of the tree, under a null model of independence where each type of event occurs at a constant rate uniformly inthephylogenetic tree. The strengths andweaknesses of themethodare assessed via simulations;we then apply the method to explore the loss of cell motility in intracellular pathogens.

Measuring epistasis in fitness landscapes: The correlation of fitness effects of mutations.

Genotypic fitness landscapes are constructed by assessing the fitness of all possible combinations of a given number of mutations. In the last years, several experimental fitness landscapes have been completely resolved. As fitness landscapes are high-dimensional, simple measures of their structure are used as statistics in empirical applications. Epistasis is one of the most relevant features of fitness landscapes. Here we propose a new natural measure of the amount of epistasis based on the correlation of fitness effects of mutations. This measure has a natural interpretation, captures well the interaction between mutations and can be obtained analytically for most landscape models. We discuss how this measure is related to previous measures of epistasis (number of peaks, roughness/slope, fraction of sign epistasis, Fourier-Walsh spectrum) and how it can be easily extended to landscapes with missing data or with fitness ranks only. Furthermore, the dependence of the correlation of fitness effects on mutational distance contains interesting information about the patterns of epistasis. This dependence can be used to uncover the amount and nature of epistatic interactions in a landscape or to discriminate between different landscape models.

Mutational pattern of a sample from a critical branching population.

We study a universal object for the genealogy of a sample in populations with mutations: the critical birth-death process with Poissonian mutations, conditioned on its population size at a fixed time horizon. We show how this process arises as the law of the genealogy of a sample in a large class of nearly critical branching populations with rare mutations at birth, namely populations converging, in a large population asymptotic, towards the continuum random tree. We extend this model to populations with random foundation times, with (potentially improper) prior distributions [Formula: see text], [Formula: see text], including the so-called uniform ([Formula: see text]) and log-uniform ([Formula: see text]) priors. We first investigate the mutational patterns arising from these models, by studying the site frequency spectrum of a sample with fixed size, i.e. the number of mutations carried by k individuals in the sample. Explicit formulae for the expected frequency spectrum of a sample are provided, in the cases of a fixed foundation time, and of a uniform and log-uniform prior on the foundation time. Second, we establish the convergence in distribution, for large sample sizes, of the (suitably renormalized) tree spanned by the sample with prior [Formula: see text] on the time of origin. We finally prove that the limiting genealogies with different priors can all be embedded in the same realization of a given Poisson point measure.

The reproducibility of adaptation in the light of experimental evolution with whole genome sequencing.

A key question in evolutionary biology is the reproducibility of adaptation. This question can now be quantitatively analyzed using experimental evolution coupled to whole genome sequencing (WGS). With complete sequence data, one can assess convergence among replicate populations. In turn, convergence reflects the action of natural selection and also the breadth of the field of possible adaptive solutions. That is, it provides insight into how many genetic solutions or adaptive paths may lead to adaptation in a given environment. Convergence is both a property of an adaptive landscape and, reciprocally, a tool to study that landscape. In this chapter we present the links between convergence and the properties of adaptive landscapes with respect to two types of microbial experimental evolution. The first tries to reconstruct a full adaptive landscape using a handful of carefully identified mutations (the reductionist approach), while the second uses WGS of replicate experiments to infer properties of the adaptive landscape. Reductionist approaches have highlighted the importance of epistasis in shaping the adaptive landscape, but have also uncovered a wide diversity of landscape architectures. The WGS approach has uncovered a very high diversity of beneficial mutations that affect a limited set of genes or functions and also suggests some shortcomings of the reductionist approach. We conclude that convergence may be better defined at an integrated level, such as the genic level or even at a phenotypic level, and that integrated mechanistic models derived from systems biology may offer an interesting perspective for the analysis of convergence at all levels.

Chromosome-level genome assembly of the European green woodpecker Picus viridis.

The European green woodpecker, Picus viridis, is a widely distributed species found in the Western Palearctic region. Here, we assembled a highly contiguous genome assembly for this species using a combination of short- and long-read sequencing and scaffolded with chromatin conformation capture (Hi-C). The final genome assembly was 1.28 Gb and features a scaffold N50 of 37 Mb and a scaffold L50 of 39.165 Mb. The assembly incorporates 89.4% of the genes identified in birds in OrthoDB. Gene and repetitive content annotation on the assembly detected 15,805 genes and a ∼30.1% occurrence of repetitive elements, respectively. Analysis of synteny demonstrates the fragmented nature of the P. viridis genome when compared to the chicken (Gallus gallus). The assembly and annotations produced in this study will certainly help for further research into the genomics of P. viridis and the comparative evolution of woodpeckers. Five historical and seven contemporary samples have been resequenced and may give insights on the population history of this species.

A General Framework for Neutrality Tests Based on the Site Frequency Spectrum.

One of the main necessities for population geneticists is the availability of sensitive statistical tools that enable to accept or reject the standard Wright-Fisher model of neutral evolution. A number of statistical tests have been developed to detect specific deviations from the null frequency spectrum in different directions (e.g., Tajima's D, Fu and Li's F and D tests, Fay and Wu's H). A general framework exists to generate all neutrality tests that are linear functions of the frequency spectrum. In this framework, it is possible to develop a family of optimal tests with almost maximum power against a specific alternative evolutionary scenario. In this paper we provide a thorough discussion of the structure and properties of linear and nonlinear neutrality tests. First, we present the general framework for linear tests and emphasise the importance of the property of scalability with the sample size (that is, the interpretation of the tests should not depend on the sample size), which, if missing, can lead to errors in interpreting the data. After summarising the motivation and structure of linear optimal tests, we present a more general framework for the optimisation of linear tests, leading to a new family of tunable neutrality tests. In a further generalisation, we extend the framework to nonlinear neutrality tests and we derive nonlinear optimal tests for polynomials of any degree in the frequency spectrum.

Cultural transmission of reproductive success impacts genomic diversity, coalescent tree topologies, and demographic inferences.

Cultural transmission of reproductive success has been observed in many human populations as well as other animals. Cultural transmission of reproductive success consists of a positive correlation of nongenetic origin between the progeny size of parents and children. This correlation can result from various factors, such as the social influence of parents on their children, the increase of children's survival through allocare from uncles and aunts, or the transmission of resources. Here, we study the evolution of genomic diversity over time under cultural transmission of reproductive success. Cultural transmission of reproductive success has a threefold impact on population genetics: (1) the effective population size decreases when cultural transmission of reproductive success starts, mimicking a population contraction, and increases back to its original value when cultural transmission of reproductive success stops; (2) coalescent tree topologies are distorted under cultural transmission of reproductive success, with higher imbalance and a higher number of polytomies; and (3) branch lengths are reduced nonhomogenously, with a higher impact on older branches. Under long-lasting cultural transmission of reproductive success, the effective population size stabilizes but the distortion of tree topology and the nonhomogenous branch length reduction remain, yielding U-shaped site frequency spectra under a constant population size. We show that this yields a bias in site frequency spectra-based demographic inference. Considering that cultural transmission of reproductive success was detected in numerous human and animal populations worldwide, one should be cautious because inferring population past histories from genomic data can be biased by this cultural process.