The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic ...evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations.
Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females.
ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Sex-determining systems often undergo high rates of turnover but for reasons that remain largely obscure. Two recent evolutionary models assign key roles, respectively, to sex-antagonistic (SA) ...mutations occurring on autosomes and to deleterious mutations accumulating on sex chromosomes. These two models capture essential but distinct key features of sex-chromosome evolution; accordingly, they make different predictions and present distinct limitations. Here we show that a combination of features from the two models has the potential to generate endless cycles of sex-chromosome transitions: SA alleles accruing on a chromosome after it has been co-opted for sex induce an arrest of recombination; the ensuing accumulation of deleterious mutations will soon make a new transition ineluctable. The dynamics generated by these interactions share several important features with empirical data, namely, (i) that patterns of heterogamety tend to be conserved during transitions and (ii) that autosomes are not recruited randomly, with some chromosome pairs more likely than others to be co-opted for sex.
In sharp contrast with mammals and birds, many cold-blooded vertebrates present homomorphic sex chromosomes. Empirical evidence supports a role for frequent turnovers, which replace nonrecombining ...sex chromosomes before they have time to decay. Three main mechanisms have been proposed for such turnovers, relying either on neutral processes, sex-ratio selection, or intrinsic benefits of the new sex-determining genes (due, e.g., to linkage with sexually antagonistic mutations). Here, we suggest an additional mechanism, arising from the load of deleterious mutations that accumulate on nonrecombining sex chromosomes. In the absence of dosage compensation, this load should progressively lower survival rate in the heterogametic sex. Turnovers should occur when this cost outweighs the benefits gained from any sexually antagonistic genes carried by the nonrecombining sex chromosome. We use individual-based simulations of a Muller's ratchet process to test this prediction, and investigate how the relevant parameters (effective population size, strength and dominance of deleterious mutations, size of nonrecombining segment, and strength of sexually antagonistic selection) are expected to affect the rate of turnovers.
The recent advances of new genomic technologies have enabled the identification and characterization of sex chromosomes in an increasing number of nonmodel species, revealing that many plants and ...animals undergo frequent sex chromosome turnovers. What evolutionary forces drive these turnovers remains poorly understood, but it was recently proposed that drift might play a more important role than generally assumed. We analysed the dynamics of different types of turnovers using individual‐based simulations and show that when mediated by genetic drift, turnovers are usually easier to achieve than substitutions at neutral markers, but that their dynamics and relative likelihoods vary with the type of the resident and emergent sex chromosome system (XY and/or ZW) and the dominance relationships among the sex‐determining factors. Focusing on turnovers driven by epistatically dominant mutations, we find that drift‐mediated turnovers that preserve the heterogamety pattern are 2–4× more likely than those along which the heterogametic sex changes. This ratio nevertheless decreases along with effective population size and can even reverse in case of extreme polygyny. This can be attributed to a ‘drift‐induced’ selective force, known to influence transitions between male and female heterogamety, but which according to our study does not affect turnovers that preserve the heterogametic sex.
Sex determination is often seen as a dichotomous process: individual sex is assumed to be determined either by genetic (genotypic sex determination, GSD) or by environmental factors (environmental ...sex determination, ESD), most often temperature (temperature sex determination, TSD). We endorse an alternative view, which sees GSD and TSD as the ends of a continuum. Both effects interact a priori, because temperature can affect gene expression at any step along the sex-determination cascade. We propose to define sex-determination systems at the population- (rather than individual) level, via the proportion of variance in phenotypic sex stemming from genetic versus environmental factors, and we formalize this concept in a quantitative-genetics framework. Sex is seen as a threshold trait underlain by a liability factor, and reaction norms allow modeling interactions between genotypic and temperature effects (seen as the necessary consequences of thermodynamic constraints on the underlying physiological processes). As this formalization shows, temperature changes (due to e.g., climatic changes or range expansions) are expected to provoke turnovers in sex- determination mechanisms, by inducing large-scale sex reversal and thereby sex-ratio selection for alternative sex-determining genes. The frequency of turnovers and prevalence of homomorphic sex chromosomes in cold-blooded vertebrates might thus directly relate to the temperature dependence in sex-determination mechanisms.
Deleterious mutations accumulating on non-recombining Y chromosomes can drive XY to XY turnovers, as they allow to replace the old mutation-loaded Y by a new mutation-free one. The same process is ...thought to prevent XY to ZW turnovers, because the latter requires fixation of the ancestral Y, assuming dominance of the emergent feminizing mutation. Using individual-based simulations, we explored whether and how an epistatically dominant W allele can spread in a young XY system that gradually accumulates deleterious mutations. We also investigated how sexually antagonistic (SA) polymorphism on the ancestral sex chromosomes and the mechanism controlling X-Y recombination suppression affect these transitions. In contrast with XY to XY turnovers, XY to ZW turnovers cannot be favored by Y chromosome mutation load. If the arrest of X-Y recombination depends on genotypic sex, transitions are strongly hindered by deleterious mutations, and totally suppressed by very small SA cost, because deleterious mutations and female-detrimental SA alleles would have to fix with the Y. If, however, the arrest of X-Y recombination depends on phenotypic sex, X and Y recombine in XY ZW females, allowing for the purge of Y-linked deleterious mutations and loss of the SA polymorphism, causing XY to ZW turnovers to occur at the same rate as in the absence of deleterious and sex-antagonistic mutations. We generalize our results to other types of turnovers (e.g., triggered by non-dominant sex-determining mutations) and discuss their empirical relevance.
An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. ...Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species almost equal to141 thousand years ago (Kya), an exit out-of-Africa almost equal to51 Kya, and a recent colonization of the Americas almost equal to10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.
The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density ...SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Owing to technological advances in ancient DNA, it is now possible to sequence viruses from the past to track down their origin and evolution. However, ancient DNA data is considerably more degraded ...and contaminated than modern data making the identification of ancient viral genomes particularly challenging. Several methods to characterise the modern microbiome (and, within this, the virome) have been developed; in particular, tools that assign sequenced reads to specific taxa in order to characterise the organisms present in a sample of interest. While these existing tools are routinely used in modern data, their performance when applied to ancient microbiome data to screen for ancient viruses remains unknown. In this work, we conducted an extensive simulation study using public viral sequences to establish which tool is the most suitable to screen ancient samples for human DNA viruses. We compared the performance of four widely used classifiers, namely Centrifuge, Kraken2, DIAMOND and MetaPhlAn2, in correctly assigning sequencing reads to the corresponding viruses. To do so, we simulated reads by adding noise typical of ancient DNA to a set of publicly available human DNA viral sequences and to the human genome. We fragmented the DNA into different lengths, added sequencing error and C to T and G to A deamination substitutions at the read termini. Then we measured the resulting sensitivity and precision for all classifiers. Across most simulations, more than 228 out of the 233 simulated viruses were recovered by Centrifuge, Kraken2 and DIAMOND, in contrast to MetaPhlAn2 which recovered only around one third. Overall, Centrifuge and Kraken2 had the best performance with the highest values of sensitivity and precision. We found that deamination damage had little impact on the performance of the classifiers, less than the sequencing error and the length of the reads. Since Centrifuge can handle short reads (in contrast to DIAMOND and Kraken2 with default settings) and since it achieve the highest sensitivity and precision at the species level across all the simulations performed, it is our recommended tool. Regardless of the tool used, our simulations indicate that, for ancient human studies, users should use strict filters to remove all reads of potential human origin. Finally, we recommend that users verify which species are present in the database used, as it might happen that default databases lack sequences for viruses of interest.
Although it is generally accepted that geography is a major factor shaping human genetic differentiation, it is still disputed how much of this differentiation is a result of a simple process of ...isolation-by-distance, and if there are factors generating distinct clusters of genetic similarity. We address this question using a geographically explicit simulation framework coupled with an Approximate Bayesian Computation approach. Based on six simple summary statistics only, we estimated the most probable demographic parameters that shaped modern human evolution under an isolation by distance scenario, and found these were the following: an initial population in East Africa spread and grew from 4000 individuals to 5.7 million in about 132 000 years. Subsequent simulations with these estimates followed by cluster analyses produced results nearly identical to those obtained in real data. Thus, a simple diffusion model from East Africa explains a large portion of the genetic diversity patterns observed in modern humans. We argue that a model of isolation by distance along the continental landmasses might be the relevant null model to use when investigating selective effects in humans and probably many other species.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK