X chromosomes in SLiM 5

9 views
Skip to first unread message

Filip Ruzicka

unread,
Mar 26, 2026, 5:09:34 AM (11 days ago) Mar 26
to slim-discuss
Hi,

I just wanted to clarify the default behaviour with recombination and sex chromosomes in SLiM 5, when using options initializeChromosome(..., type="X") and initializeRecombinationRate(...). 

On the one hand, initializeChromosome() with "type="X" assumes haploidy in males ("an X chromosome that is diploid (XX) in females, haploid (X–) in males." per the manual), which I assume means no recombination in males. On the other hand, initializeRecombinationRate() without using the sex option assumes recombination in both sexes ("the supplied recombination rate map will be used for both sexes", per the manual). 

My guess is that the default initializeChromosome("type=X") behaviour "overrides" the default initializeRecombinationRate() behaviour? And there is in fact no recombination on an X chromosome in males? I just wanted to double-check if this is true. Below I paste my initialisation set-up. 

Thanks!
Filip


initialize() {
if (exists("slimgui")){
defineConstant("LO",999999); // number of bp
defineConstant("N",500); // population size of individuals
defineConstant("ST",N*10); // start time at which beneficial mutation is introduced
defineConstant("END",N*110); // end simulation 
defineConstant("MU",1e-8); // mutation rate
defineConstant("R",1e-8); // recombination rate
defineConstant("HB",0.5); // dominance coefficient for the beneficial mutation
defineConstant("SB",0.01); // selection coefficient for the beneficial mutation
defineConstant("SD",-0.02); // selection coefficient for the deleterious mutation
defineConstant("C",0.05); // fraction of constrained sites
defineConstant("HD",0.5); // dominance coefficient for the deleterious mutation
defineConstant("R_RATIO",1); // X/A ratio of recombination rates
defineConstant("MU_RATIO",1); // X/A ratio of mutation rates
}
initializeSex();

// m1 mutation type: neutral
initializeMutationType("m1", 0.5, "f", 0.0);
// m2 mutation type: deleterious
initializeMutationType("m2", HD, "f", SD);

// g1 genomic element type containing neutral and deleterious alleles
initializeGenomicElementType("g1", c(m1,m2), c(1-C,C));
// initialise an X chromosome
initializeChromosome(1, LO, type="X", symbol="X");
// initialise genomic element
initializeGenomicElement(g1);
// initialise mutation rate
initializeMutationRate(MU*MU_RATIO);
// initialise recombination rate
initializeRecombinationRate(R*R_RATIO);

// m3 mutation type: beneficial
initializeMutationType("m3", HB, "f", SB);

}


Peter Ralph

unread,
Mar 26, 2026, 11:02:37 AM (11 days ago) Mar 26
to slim-d...@googlegroups.com
Hello, Filip - yes, this is correct. (There's no other chromosome for it to recombine with!)

Happy slimulating,
  Peter

From: slim-d...@googlegroups.com <slim-d...@googlegroups.com> on behalf of Filip Ruzicka <ruzick...@gmail.com>
Sent: Thursday, March 26, 2026 2:09 AM
To: slim-discuss <slim-d...@googlegroups.com>
Subject: X chromosomes in SLiM 5
 
Hi, I just wanted to clarify the default behaviour with recombination and sex chromosomes in SLiM 5, when using options initializeChromosome(. . . , type="X") and initializeRecombinationRate(. . . ). On the one hand, initializeChromosome() with "type="X"
ZjQcmQRYFpfptBannerStart
You haven't previously corresponded with this sender.
Use caution with links and attachments. Learn more about this email warning tag.
 
ZjQcmQRYFpfptBannerEnd
--
SLiM forward genetic simulation: http://messerlab.org/slim/
---
You received this message because you are subscribed to the Google Groups "slim-discuss" group.
To unsubscribe from this group and stop receiving emails from it, send an email to slim-discuss...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/slim-discuss/0e5b9ac6-f424-4296-89ae-403cb1078e34n%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages