Skip to content

Releases: popsim-consortium/stdpopsim

Major feature release

19 Mar 02:35
d7e5acc
Compare
Choose a tag to compare

Major feature release adding documented, stable support for simulating selection (via SLiM),
as well as new species and models.

Bug fixes:

  • Incorrect scaling of the DroMel/LognormalPlusPositive_R16 DFE led to negative selection
    coefficients that were too large (many on the order of -1e3); scaling was fixed for the
    log scale (@clararehmann, #1699)

  • Factor of two error fixed in HomSap/Gamma_K17 DFE (@RyanGutenkunst, #1478)

  • Discretization for exponential growth using the SLiM engine led to the final size
    in growing/declining populations not matching (by a small amount) what is listed
    in the model; improved discretization scheme (@petrelharp, #1622)

  • Updates to SLiM support: updated the active flag in the SLiM code to be integer.`

Breaking changes:

  • The 2018 Browning et al HomSap demographic model previously named
    "AmericanAdmixture_4B11" is now named "AmericanAdmixture_4B18"; the old name
    works but throws a deprecation warning. (@petrelharp, #1603)

  • The time_units attribute of tree sequence metadata is now set to "generations"
    for SLiM output, avoiding the "time units mismatch warning".
    (@nspope, #1567)

  • Previously, Contigs specified with left and right arguments produced
    simulations having coordinates shifted so they were relative to left. This
    is no longer the case: coordinates in resulting simulations retain the
    original chromosome's coordinates. As a result, regions outside [left, right) will now have missing data. To produce simulation output with the
    previous behavior, use the .trim() method of the resulting TreeSequence.
    So, the Contig's .original_coordinates attribute is now called
    .coordinates, and several methods (e.g., Contig.basic_contig()) now take
    left and right arguments as well as length. Finally, arguments such as
    relative_coordinates=True/False to Contig.dfe_breakpoints() are no longer
    necessary and are deprecated. (@petrelharp, #1570)

  • To add the possibility for a relationship between dominance and selection
    coefficient, now each stdpopsim MutationType might have more than one
    underlying SLiM mutation type; so, where this is recorded in top-level
    metadata (under ts.metadata["stdpopsim"]["DFEs"]) is now a list
    instead of a single value. This will not affect anyone who is not
    parsing the metadata related to DFEs.

  • SLiM extended events and selective sweep infrastructure have been
    moved from the stdpopsim.ext namespace into stdpopsim proper.

  • The length_multiplier option to Species.get_contig is deprecated and
    prints a warning. The options left and right should be used to truncate a
    contig, instead. (@nspope, #1605)

  • Added assembly_source and assembly_build_version properties to
    Genome objects, so that different genomes could be based on
    different sources (previously, all were with reference to the Ensembl
    release 103). (@andrewkern, #646

New features:

  • Relationship between dominance and selection coefficient:
    Added the dominance_coeff_list argument to MutationType, allowing
    for DFEs with a discretized relationship between h and s.
    (@petrelharp, #1498)

  • Model specific recombination rates:
    Added the recombination_rate attribute to demographic model,
    allowing for model specific recombination rates in addition to
    the species default rates.
    (@gregorgorjanc, #1591)

  • Added support for the SLiM engine on Windows. (@petrelharp, #1571)

New species:

New DFEs:

New demographic models:

New annotations:

New QCs:

Alpha release for 0.2.1: additional DFEs, species, and maintenance

06 Jul 16:44
Compare
Choose a tag to compare

This is an alpha release in preparation for an 0.2.1 release, incorporating a number of additional DFEs and species, and some maintenance.

Major feature release

02 Nov 17:13
d2e4371
Compare
Choose a tag to compare

Major feature release adding many new species and models, as well as support
for simulating selection via SLiM.

Bug fixes:

  • Parameters in the HomSap/Zigzag_1S14 model now match those in Schiffels &
    Durbin (2014) (@grahamgower, #750).

  • Recombination rate for DroMel chr4 changed to 0
    (@izabelcavassim, #1092).

  • Per-chromosome mean recombination rates for HomSap were incorrectly
    calculated (@nspope, #1345).

Breaking changes:

  • Removed GeneticMap class from public API (@jeromekelleher, #713).

  • Samples are now specified via population/individual pairs, using
    species/chromosome ploidy. The old API for specifying haploid samples via
    population index has been retained, but is deprecated and will be
    removed at some point (@nspope, #1361).

New species:

New models:

New genetic maps:

New features:

Additions to CLI:

  • Sample specification has switched from positional and haploid (e.g. stdpopsim HomSap -d OutOfAfrica_3G09 6 0 10) to named with species-specific ploidy (equivalent to stdpopsim HomSap -d OutOfAfrica_3G09 YRI:3 CEU:0 CHB:5). Positional sample specification is still supported but will raise a deprecation warning.

  • Arguments --dfe, --dfe-interval, --dfe-bed-file, --help-dfe
    for specifying DFEs (@izabelcavassim, #1052).

  • Arguments --help-annotations, --dfe-annotation for associating annotation
    tracks with DFEs (@andrewkern, #1117).

  • Argument --length for simulating from a generic contig
    (@apragsdale, #664).

  • Arguments --inclusion-mask, --exclusion-mask for masking simulated sequences
    (@apragsdale, #664).

  • Arguments --left and --right for simulating an interval on a named chromosome
    (@nspope, #1348)

  • Argument --keep-mutation-ids-as-alleles retains SLiM mutation IDs for
    allele codes instead of converting these to nucleotides (@nspope, #1356).

Catalog maintenance infrastructure:

New annotations:

New DFEs:

Important bugfix and feature release

29 May 08:19
f6c583d
Compare
Choose a tag to compare

Important bugfix and feature release, recommended for all users.

Significant errors in the HomSap/OutOfAfrica_3G09 and HomSap/OutOfAfrica_2T12 models have been fixed. We recommend repeating any analyses performed using these models. See here for more details on the error in the three population Out of Africa model and analysis of the differences from the correct model.

The recombination rate for AraTha was also off by a factor of 10. We recommend repeating any analyses performed using this species.

Bug fixes:

  • Fix error in HomSap/OutOfAfrica_3G09 model, in which migration between ancestral African and European populations was allowed to continue in the most ancient time period (@apragsdale, #496, #516).

  • Fix similar error in HomSap/OutOfAfrica_2T12 model (@ndukler, #520, #516).

  • Fix recombination rate estimate for AraTha (@grahamgower, #537, #527), which was off by a factor of 10.

  • Require attrs >=19.10 (@grahamgower, #399, #394)

New species:

New models:

New features:

Breaking changes:

  • Move the –quiet/-q command line option to the top-level. Previously we would write stdpopsim HomSap -q 10 whereas we now write stdpopsim -q HomSap. (@jeromekelleher, #515, #547)

  • The long form --verbosity argument has been changed to --verbose (#547).

  • Removed DroMel chrM (@grahamgower, #528, #405).

Bugfix release

02 Jan 18:18
ec99cb2
Compare
Choose a tag to compare

Bugfix release. Fixes some distribution issues and temporarily removes the PonPyg species.

Bug fixes:

New features:

Breaking changes:

  • Temporarily remove the PonPyg species from the catalog to provide time to fix issues with genomes and multi-species models (#365).

Initial release

18 Dec 18:26
41db8d5
Compare
Choose a tag to compare

Initial release of the stdpopsim resource.