A Population Genetic Model of Evolution of Host-Mate Attraction and Nonhost Repulsion in a Bark Beetle Pityogenes bidentatus

Algorithm 1

General code for algorithms to calculate the proportions P₁ to P₉ of the nine mated genotypes (AABB to aabb, see Table 1) for 200 generations based on initial proportions and nine constant survival factors S₁ to S₉ for individuals of these genotypes (an asterisk denotes multiplication).

AABB: P₁ = 0: S₁ = 1: AABb: P₂ = 0: S₂ = 0.8: AAbb: P₃ = 0: S₃ = 0.5

AaBB: P₄ = 0: S₄ = 0.9: AaBb: P₅ = 0: S₅ = 0.6: Aabb: P₆ = 0.0000001: S₆ = 0.3

aaBB: P₇ = 0: S₇ = 0.4: aaBb: P₈ = 0.0000001: S₈ = 0.3: aabb: P₉ = 0.9999998: S₉ = 0.2

For generation = 1 to 100

Row 1: AABB: G₁=16^∗P₁^∗P₁+16^∗P₁^∗P₂+16^∗P₁^∗P₄+8^∗P₁^∗P₅

AABb: G₂=16^∗P₁^∗P₂+32^∗P₁^∗P₃+8^∗P₁^∗P₅+ 16^∗P₁^∗P₆

AaBB: G₄=16^∗P₁^∗P₄+8^∗P₁^∗P₅+32^∗P₁^∗P₇+ 16^∗P₁^∗P₈

AaBb: G₅=8^∗P₁^∗P₅+ 16^∗P₁^∗P₆+ 16^∗P₁^∗P₈+32^∗P₁^∗P₉

Row 2: AABB: G₁=G₁+4^∗P₂^∗P₂+8^∗P₂^∗P₄+4^∗P₂^∗P₅

AABb: G₂=G₂+8^∗P₂^∗P₂+ 16^∗P₂*P₃+8^∗P₂^∗P₄ +8^∗P₂^∗P₅+8^∗P₂^∗P₆

AAbb: G₃=4^∗P₂^∗P₂+ 16^∗P₂^∗P₃+4^∗P₂^∗P₅+8^∗P₂^∗P₆

AaBB: G₄=G₄+8^∗P₂^∗P₄+4^∗P₂^∗P₅+16^∗P₂^∗P₇+8^∗P₂^∗P₈

AaBb: G₅=G₅+8^∗P₂^∗P₄+8^∗P₂^∗P₅+8^∗P₂^∗P₆+ 16^∗P₂^∗P₇+ 16^∗P₂^∗P₈+ 16^∗P₂^∗P₉

Aabb: G₆=4^∗P₂^∗P₅+8^∗P₂^∗P₆+8^∗P₂^∗P₈+ 16^∗P₂^∗P₉

Row 3: AAbb: G₃=G₃+16^∗P₃^∗P₃+8^∗P₃*P₅+16^∗P3^∗P₆

AABb: G₂=G₂+ 16^∗P₃^∗P₄+8^∗P₃^∗P₅

AaBb: G₅=G₅+ 16^∗P₃^∗P₄+8^∗P₃^∗P₅+32^∗P₃^∗P₇+ 16^∗P₃^∗P₈

Aabb: G₆=G₆+8^∗P₃^∗P₅+ 16^∗P₃^∗P₆+ 16^∗P₃^∗P₈+32^∗P₃^∗P₉

Row 4: AABB: G₁=G₁+4^∗P₄^∗P₄+4^∗P₄^∗P₅

AaBB: G₄=G₄+8^∗P₄^∗P₄+8^∗P₄^∗P₅+ 16^∗P₄^∗P₇+8^∗P₄^∗P₈

aaBB: G₇=4^∗P 4^∗P₄+4^∗P₄^∗P₅+ 16^∗P₄^∗P₇+8^∗P₄^∗P₈

AABb: G₂=G₂+4^∗P₄^∗P₅+8^∗P 4^∗P₆

AaBb: G₅=G₅+8^∗P₄^∗P₅+ 16^∗P₄^∗P₆+8^∗P₄^∗P₈+ 16^∗P 4^∗P₉

aaBb:G₈=4^∗P₄^∗P₅+8^∗P₄^∗P₆+8^∗P₄^∗P₈+16^∗P₄^∗P₉

Row 5: AABB: G₁=G₁+P₅^∗P₅: AABb: G₂=G₂+2^∗P₅^∗P₅+4^∗P₅^∗P₆

AAbb: G₃=G₃+P₅^∗P₅+4^∗P₅^∗P₆: AaBB: G₄=G₄+2^∗P₅^∗P₅+8^∗P₅^∗P₇+4^∗P₅^∗P₈

AaBb: G₅=G₅+4 ^∗P₅^∗P₅+8^∗P₅^∗P₆+8^∗P₅^∗P₇+8^∗P₅^∗P₈+8 ^∗P₅^∗P₉

Aabb: G₆=G₆+2^∗P₅^∗P₅+8^∗P₅^∗P₆+4^∗P₅^∗P₈+8^∗P₅^∗P₉

aaBB: G₇=G₇+P₅^∗P₅+8^∗P₅^∗P₇+4^∗P₅^∗P₈

aaBb: G₈=G₈+2^∗P₅^∗P₅+4 ^∗P₅^∗P₆+8^∗P₅^∗P₇+8^∗P₅^∗P₈+8^∗P₅^∗P₉

aabb: G₉=P₅^∗P₅+4^∗P₅^∗P₆+4^∗P₅^∗P₈+8^∗P₅^∗P₉

Row 6: AAbb: G₃=G₃+4^∗P₆^∗P₆: Aabb: G₆=G₆+8^∗P₆^∗P₆+8^∗P₆^∗P₈+16^∗P₆^∗P₉

aabb: G₉=G₉+4^∗P₆^∗P₆+8^∗P₆^∗P₈+ 16^∗P₆^∗P₉

AaBb: G₅=G₅+ 16^∗P₆^∗P₇+8^∗P₆^∗P₈ aaBb: G₈=G₈+ 16^∗P₆^∗P₇+8^∗P₆^∗P₈

Row 7: aaBB: G₇=G₇+16^∗P₇^∗P₇+16^∗P₇^∗P₈: aaBb: G₈=G₈+16^∗P₇^∗P₈+32^∗P₇^∗P₉

Row 8: aaBB: G₇=G₇+4^∗P₈^∗P₈: aaBb: G₈=G₈+8^∗P₈^∗P₈+16^∗P₈^∗P₉

aabb: G₉=G₉+4^∗P₈^∗P₈+ 16^∗P₈^∗P₉

Row 9: aabb: G₉=G₉+ 16^∗P₉^∗P₉

total = G₁^∗S₁ +G₂^∗S₂+G₃^∗S₃+G₄^∗S₄ +G₅^∗S₅+G₆^∗S₆+G₇^∗S₇+G₈^∗S₈+G₉^∗S₉

P₁ = G₁^∗S₁/total: P₂ = G₂^∗S₂/total: P₃ = G₃^∗S₃/total: P₄=G₄^∗S₄/total

P₅ = G₅^∗S₅/total: P₆ = G₆^∗S₆/total: P₇ = G₇^∗S₇/total: P₈ = G₈^∗S₇/total: P₉ = G₉^∗S₉/total

Next generation