(A) Joint and marginal distributions of ages (TMRCA in generations on a log10 scale) and lengths of TMRCA-segments (in bp on a log10 scale) in a selfing population (σ=0.95) with a stepwise change …
Probability of a break in the genealogy to occur in a coalescent tree of sample size two in the absence of selfing (red); in presence of 90% selfing (orange); in the case of a transition from …
(A) Explicit selfing implemented in a forward-in-time Wright-Fisher model (slim3). Population size is constant and selfing rate shifts from outcrossing (σ=0, green) to predominant selfing (σ=0.95, …
(A–I) Joint and marginal distributions of ages in generations and lengths of TMRCA-segments in a population with a constant population size and a shift from outcrossing (green) to predominant …
Inference of times of transition from outcrossing (σ=0.1) to predominant selfing (σ=0.99) using neutral simulations. The x-axes represent the true value of tσ in units of log10(generations) and the y…
Best-case convergence of teSMC using 10 sequences (i.e. haploid genomes) of 100 Mb (green) when population undergoes a bottleneck (true sizes are indicated in gray) with either variation of selfing …
Best-case convergence of teSMC using different combinations of sample sizes (, , or sequences; i.e. haploid genomes) and sequence lengths (L=10 Mb or L=100 Mb, see legend). Population size is …
Comparisons between true (black lines) and estimated selfing rates and population sizes estimated by teSMC for 10 replicates. Here, simulations were done using a constant population size (N=40,000) …
(A–P) Comparisons between true (black lines) and estimated selfing rates and population sizes estimated by teSMC for 10 replicates. Here, simulations were done as in Figure 2—figure supplement 3 …
Estimated population sizes by teSMC using 20 sequences (i.e. haploid genomes) of length 5 Mb (red) when population size is constant and set to 100,000 (black line) with a constant selfing value of …
Estimated selfing rates through time estimated by teSMC using 20 sequences (i.e. haploid genomes) of length 5 Mb (red) when population size is constant and set to 100,000 (black line) with a …
Best-case convergence of teSMC under four different scenarios of selfing transition. (A) Slow transition from outcrossing to selfing. (B) Transition from selfing to outcrossing. (C) Transition from …
Performance of teSMC on simulated sequence data under four different scenarios of selfing variations. (A) Slow transition from outcrossing to selfing. (B) Transition from selfing to outcrossing. (C) …
(A) Demographic model 1 in the model choice analysis: one population with a single transition from predominant selfing to predominant outcrossing. (B) Demographic model 2 in the model choice …
(A) The transition matrix of pairwise diversity (TMwin) between adjacent non-overlapping 10 kb windows measured for 1 Mb of data simulated with the population-size change model (Figure 3B) in an …
Parameter re-estimation of the three additional parameters of the model described in Figure 3. (A, B, C) Re-estimation of the population size on 100 datasets simulated under a model with constant …
Inference of times of transition from outcrossing (σ=0.1) to predominant selfing (σ=0.99) using (A) teSMC and (B) tsABC. Simulations were done under constant population size and negative selection …
Parameter re-estimation of the three remaining parameters of the model (see Figure 3A). (A, B, C) Re-estimation of the population size on 100 datasets simulated under a model with constant …
(A) Inferred transitions from outcrossing to selfing for three independent genetic clusters of A. thaliana from the 1001 genomes project (CEU, IBnr, Relict) using tsABC and teSMC. The 95% CI (CEU) …
Panels show, for the five chromosomes of A. thaliana, the nucleotide diversity (pi) calculated with tskit using a sliding window of 100 kb. Results are shown only for the CEU sample used (n=99) for …
(A, B) Inferred transitions from outcrossing to selfing for three ancestry groups of A. thaliana from the 1001 genomes project (CEU: central European, IBnr: Iberia non-relict, Relict) using teSMC. (C…
Estimations were conducted for three different ancestry groups: central Europe (CEU), Iberian non-relicts (IBnr), and Relicts using both teSMC and tsABC. The 95% CI of all jointly inferred …
Method | Population | Sample size | Mode | 95% credibility interval | |
---|---|---|---|---|---|
teSMC | CEU | 20 | 697,490 | NA | NA |
teSMC | IBnr | 20 | 713,421 | NA | NA |
teSMC | Relicts | 17 | 749,668 | NA | NA |
tsABC | CEU | 99 | 707,995 | 443,486 | 973,841 |
tsABC | IBnr | 66 | 756,976 | 397,049 | 988,708 |
tsABC | Relicts | 17 | 592,321 | 386,406 | 934,499 |
Supplementary tables.
Parameters for simulated datasets to investigate the performance of tsABC. Table S2: Parameter priors used for the performance analysis of tsABC. Table S3: Parameter priors used to estimate the transition from outcrossing to selfing in A. thaliana. Table S4: A. thaliana samples with >95% cluster membership (https://1001genomes.github.io/admixture-map/) used to estimate the transition from outcrossing to selfing obtained from the 1001 genomes website (https://1001genomes.org/data/GMI-MPI/releases/v3.1/SNP_matrix_imputed_hdf5/1001_SNP_MATRIX.tar.gz). Table S5: Genomic regions of A. thaliana in TAIR10 used for the tsABC and teSMC analyses. Table S6: Jointly estimated parameters in the context of transitions from predominant outcrossing to predominant selfing in A. thaliana.