Module 7

Reproduction & Development

Sexual cannibalism, sperm transfer mechanics, egg sac engineering, molting hydraulics, and sexual dimorphism

Spider reproduction involves extraordinary biomechanical challenges: males must transfer sperm via modified mouthparts (pedipalps) while potentially facing lethal aggression from larger females. Egg protection requires specialized silks that balance gas exchange with desiccation resistance. Growth demands repeated ecdysis (molting) driven by hydraulic pressure — a life-or-death process where mechanical failure is fatal. Sexual dimorphism in some species reaches extreme levels, with females outweighing males by 100-fold, creating unique evolutionary dynamics.

7.1 Sexual Cannibalism: Game Theory

In many spider species — notably Latrodectus (black widows), Argiope (garden spiders), and Nephila (golden orb weavers) — the female may cannibalize the male during or after copulation. From a biomechanical perspective, the male faces a critical optimization problem: maximizing lifetime reproductive success by balancing current mating investment against survival for future matings.

Evolutionary Game Theory Model

Consider a male who can choose mating duration \(t\). Longer copulation increases sperm transfer (and thus paternity share \(P(t)\)) but also increases cannibalism risk. The expected fitness of a male is:

Male Fitness Function

\[ W_{\text{male}}(t) = P(t) \cdot F_{\text{current}} + S(t) \cdot \sum_{k=1}^{n} P_k \cdot F_k \]

The first term is the fitness from the current mating: paternity share \(P(t)\) times the female's fecundity \(F_{\text{current}}\). The second term is expected future fitness: survival probability \(S(t)\) times the sum of expected payoffs from\(n\) future matings.

If we model paternity as a saturating function \(P(t) = P_{\max}(1 - e^{-\alpha t})\)and survival as a decreasing exponential \(S(t) = e^{-\beta t}\), the optimal mating duration \(t^*\) satisfies:

Optimal Mating Duration

\[ \frac{dW}{dt}\bigg|_{t^*} = 0 \implies \alpha P_{\max} e^{-\alpha t^*} F_{\text{current}} = \beta e^{-\beta t^*} W_{\text{future}} \]

This yields \(t^* = \frac{1}{\alpha - \beta}\ln\!\left(\frac{\alpha P_{\max} F_{\text{current}}}{\beta W_{\text{future}}}\right)\)when \(\alpha > \beta\). When future mating prospects are poor (\(W_{\text{future}} \to 0\)),\(t^* \to \infty\) — the male should mate as long as possible, even accepting cannibalism (the "adaptive suicide" strategy seen in Argiope aurantia).

Paternity Benefits of Cannibalism

In Argiope bruennichi, Welke & Schneider (2012) showed that cannibalized males achieve higher paternity than surviving males because: (1) the female is occupied eating, preventing remating, and (2) the male's pedipalp often breaks off, forming a mating plug. The payoff matrix:

Strategy

Current Paternity

Future Matings

Quick escape

~30%

0–2 additional

Extended mating

~60%

0–1 additional

Sacrifice (cannibalized)

~87%

7.2 Sperm Web & Pedipalpal Transfer

Spider males use a unique two-step sperm transfer process found nowhere else in the animal kingdom. First, the male spins a small sperm web (2–5 mm), deposits a droplet of semen from his genital pore (on the ventral abdomen), then dips his pedipalps into the droplet to load sperm into the pedipalpal bulbs — complex, inflatable structures that function as hydraulic syringes during copulation.

Capillary Uptake of Sperm

The pedipalpal bulb contains a coiled sperm duct (embolus) with an internal diameter of\(r \approx 5\text{--}20\,\mu\text{m}\). Sperm is drawn into this duct primarily by capillary action. The height to which a fluid rises in a capillary tube is given by the Jurin equation:

Capillary Rise (Jurin's Law)

\[ h = \frac{2\gamma \cos\theta}{\rho g r} \]

where \(\gamma \approx 0.05\,\text{N/m}\) is the surface tension of seminal fluid,\(\theta \approx 30°\) is the contact angle with chitin, \(\rho \approx 1050\,\text{kg/m}^3\)is the fluid density, \(g = 9.81\,\text{m/s}^2\), and \(r\) is the duct radius. For \(r = 10\,\mu\text{m}\):\(h = \frac{2 \times 0.05 \times \cos 30°}{1050 \times 9.81 \times 10^{-5}} \approx 0.84\,\text{m}\)— far more than needed, confirming capillary forces are sufficient.

7.3 Egg Sac Silk Engineering

Tubuliform (egg sac) silk is produced by the tubuliform glands and is mechanically distinct from dragline or capture silk. It must simultaneously: (1) resist physical damage from predators and parasitoid wasps, (2) block UV radiation that would damage developing embryos, (3) allow gas exchange (O\(_2\) in, CO\(_2\) out), and (4) minimize water loss while preventing flooding.

Water Vapor Transmission

The egg sac membrane acts as a selective barrier. Its water vapor transmission rate (WVTR) can be modeled using Fick's first law for diffusion through a porous membrane:

Water Vapor Transmission Rate

\[ \text{WVTR} = \frac{D \cdot S \cdot \Delta p}{L} = \frac{P \cdot \Delta p}{L} \]

where \(D\) is the diffusion coefficient, \(S\) is the solubility coefficient,\(P = DS\) is the permeability, \(\Delta p\) is the water vapor pressure difference across the membrane, and \(L\) is the membrane thickness. Measured values for tubuliform silk: \(P \approx 2 \times 10^{-12}\,\text{mol}\cdot\text{m}/(\text{m}^2\cdot\text{s}\cdot\text{Pa})\), giving WVTR \(\approx 200\,\text{g}/\text{m}^2/\text{day}\) at 25°C, 50% RH.

Mechanical Properties

Tubuliform silk fibers have unique mechanical properties optimized for impact protection:

Tensile strength: ~500 MPa (lower than dragline's ~1.1 GPa but still strong)
Extensibility: ~30% strain at break (similar to Kevlar)
Toughness: ~100 MJ/m\(^3\) — highest toughness-to-weight ratio of any silk
Fiber diameter: 5–8 \(\mu\)m (thicker than dragline's 1–4 \(\mu\)m)

UV Protection

The outer silk layer of the egg sac acts as a UV filter, absorbing wavelengths below 350 nm that could damage developing embryonic DNA. The UV attenuation follows the Beer–Lambert law through the silk matrix:

UV Attenuation through Egg Sac Silk

\[ I(\lambda) = I_0(\lambda) \cdot e^{-\alpha(\lambda) \cdot L_{\text{eff}}} \]

where \(\alpha(\lambda)\) is the wavelength-dependent absorption coefficient of the silk protein matrix. At \(\lambda = 300\,\text{nm}\),\(\alpha \approx 50\,\text{mm}^{-1}\), giving >99% UV blocking through a 0.1 mm silk layer. The aromatic amino acids (tyrosine, tryptophan, phenylalanine) in tubuliform spidroin provide intrinsic UV absorption centered at 280 nm.

Thermal Regulation

The multi-layered egg sac also functions as a thermal buffer. The air gaps between silk layers create insulation, and the overall thermal resistance can be modeled as resistors in series:

\[ R_{\text{thermal}} = \sum_{i} \frac{L_i}{k_i A} = \frac{L_{\text{silk}}}{k_{\text{silk}} A} + \frac{L_{\text{air}}}{k_{\text{air}} A} \]

Silk thermal conductivity \(k_{\text{silk}} \approx 0.25\,\text{W/(m\cdot K)}\) and air \(k_{\text{air}} \approx 0.025\,\text{W/(m\cdot K)}\). The air gaps dominate thermal resistance, maintaining egg temperature 2–5°C above ambient during cool nights through retained metabolic heat from developing embryos.

7.4 Molting Mechanics: Hydraulic Ecdysis

Spiders grow by molting (ecdysis) — shedding the old exoskeleton (exuvium) and expanding the new, soft cuticle before it hardens. Young spiders molt approximately 5–10 times before reaching adulthood. The process is critically dependent on hemolymph pressure: the spider must generate sufficient internal pressure to split the old exoskeleton along predetermined fracture lines (ecdysial sutures) and then inflate the new, larger exoskeleton.

Pressure Vessel Analysis

The old exoskeleton can be modeled as a thin-walled pressure vessel. The internal pressure needed to initiate fracture along the ecdysial suture is determined by the hoop stress criterion:

Thin-Walled Pressure Vessel

\[ P_{\text{burst}} = \frac{\sigma_{\text{cuticle}} \cdot t}{r} \]

where \(\sigma_{\text{cuticle}}\) is the tensile strength of the cuticle at the ecdysial suture (~10–50 MPa, deliberately weakened by enzymatic digestion),\(t\) is the cuticle thickness, and \(r\) is the local radius of curvature. For a spider with \(r = 2\,\text{mm}\), \(t = 20\,\mu\text{m}\),\(\sigma = 20\,\text{MPa}\):\(P = 20 \times 10^6 \times 20 \times 10^{-6} / 2 \times 10^{-3} = 200\,\text{kPa}\).

Scaling of Molting Pressure

As spiders grow, the required burst pressure changes with body size. If cuticle thickness scales as \(t \propto r^{0.8}\) (negative allometry — cuticle gets relatively thinner) and suture strength remains roughly constant, then:

\[ P_{\text{burst}} \propto \frac{t}{r} \propto \frac{r^{0.8}}{r} = r^{-0.2} \]

This means larger spiders require proportionally less pressure to molt, which is fortunate because maximum hemolymph pressure also has physiological limits. However, the absolute volume of hemolymph required increases as \(V \propto r^3\), making each successive molt more metabolically costly.

Molt Failure

If the spider cannot generate sufficient hemolymph pressure — due to dehydration, injury, or metabolic insufficiency — it becomes trapped in the old exoskeleton and dies. In captive tarantulas, molt failure rates of 5–15% have been observed, often correlated with low humidity (desiccation of the ecdysial membrane) or nutritional deficiency.

Sexual Cannibalism Game Theory & Molting Pressure Analysis

Python

script.py186 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 10))

# ========== Panel 1: Cannibalism Game Theory Payoff ==========
ax1 = axes[0, 0]

t = np.linspace(0.1, 20, 200)  # mating duration (minutes)
alpha = 0.3   # paternity saturation rate
P_max = 0.9   # max paternity share
beta_vals = [0.05, 0.1, 0.2, 0.4]  # cannibalism hazard rates

F_current = 100  # current female fecundity (eggs)
W_future_vals = [0, 50, 100, 200]  # expected future fitness

# Paternity function
P = P_max * (1 - np.exp(-alpha * t))

# Plot fitness for different future prospects
colors = ['#ef4444', '#f97316', '#a78bfa', '#22c55e']
for i, W_f in enumerate(W_future_vals):
    beta = 0.15  # fixed hazard rate
    S = np.exp(-beta * t)  # survival probability
    W = P * F_current + S * W_f
    ax1.plot(t, W, color=colors[i], linewidth=2, label=f'W_future = {W_f}')
    opt_idx = np.argmax(W)
    ax1.plot(t[opt_idx], W[opt_idx], 'o', color=colors[i], markersize=10,
             markeredgecolor='white', markeredgewidth=1.5, zorder=5)

ax1.set_xlabel('Mating Duration (min)', fontsize=12, color='white')
ax1.set_ylabel('Expected Fitness', fontsize=12, color='white')
ax1.set_title('Optimal Mating Duration vs Future Prospects', fontsize=14, fontweight='bold', color='#c4b5fd')
ax1.legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#7c3aed', labelcolor='white')
ax1.set_facecolor('#0a0a1a')
ax1.tick_params(colors='white')
for spine in ax1.spines.values():
    spine.set_color('#7c3aed')

# ========== Panel 2: Survival vs Paternity Trade-off ==========
ax2 = axes[0, 1]

t_range = np.linspace(0.1, 15, 200)
beta = 0.15
P_t = P_max * (1 - np.exp(-alpha * t_range))
S_t = np.exp(-beta * t_range)

ax2_twin = ax2.twinx()

l1, = ax2.plot(t_range, P_t * 100, color='#a78bfa', linewidth=2.5, label='Paternity share')
l2, = ax2_twin.plot(t_range, S_t * 100, color='#f472b6', linewidth=2.5, label='Survival prob.')

# Shade the danger zone
ax2.fill_between(t_range, 0, P_t * 100, where=S_t < 0.5, alpha=0.1, color='#ef4444')

ax2.set_xlabel('Mating Duration (min)', fontsize=12, color='white')
ax2.set_ylabel('Paternity Share (%)', fontsize=12, color='#a78bfa')
ax2_twin.set_ylabel('Survival Probability (%)', fontsize=12, color='#f472b6')
ax2.set_title('The Cannibalism Trade-off', fontsize=14, fontweight='bold', color='#c4b5fd')
ax2.set_facecolor('#0a0a1a')
ax2.tick_params(axis='y', colors='#a78bfa')
ax2.tick_params(axis='x', colors='white')
ax2_twin.tick_params(colors='#f472b6')
for spine in ax2.spines.values():
    spine.set_color('#7c3aed')
for spine in ax2_twin.spines.values():
    spine.set_color('#7c3aed')

ax2.legend(handles=[l1, l2], fontsize=10, facecolor='#1a1a2e', edgecolor='#7c3aed',
           labelcolor='white', loc='center right')

# ========== Panel 3: Molting Pressure vs Body Size ==========
ax3 = axes[1, 0]

radii = np.linspace(0.2, 15, 100)  # body radius in mm
sigma_suture = 20e6  # suture strength (Pa) - weakened cuticle

# Cuticle thickness scaling: t = k * r^0.8
k_t = 10e-6 / (1.0**0.8)  # calibrated: 10 um at r=1mm
t_cuticle = k_t * (radii**0.8)  # meters (radii in mm -> convert)
t_cuticle_m = k_t * ((radii * 1e-3)**0.8)
r_m = radii * 1e-3

P_burst = sigma_suture * t_cuticle_m / r_m  # Pa
P_burst_kPa = P_burst / 1000

# Maximum hemolymph pressure (physiological limit)
P_max_hemo = np.ones_like(radii) * 60  # kPa (typical max for spiders)

# Different suture strengths
for sigma, ls, label in [(10e6, '--', '\u03C3 = 10 MPa'),
                          (20e6, '-', '\u03C3 = 20 MPa'),
                          (40e6, ':', '\u03C3 = 40 MPa')]:
    P_b = sigma * t_cuticle_m / r_m / 1000
    ax3.plot(radii, P_b, linewidth=2, linestyle=ls, color='#a78bfa', label=label)

ax3.axhline(y=60, color='#ef4444', linewidth=2, linestyle='--', label='Max hemolymph pressure')
ax3.fill_between(radii, 60, 300, alpha=0.1, color='#ef4444')
ax3.text(10, 70, 'MOLT FAILURE ZONE', color='#ef4444', fontsize=11, fontweight='bold')

# Annotate instars
instar_radii = [0.3, 0.5, 0.8, 1.2, 1.8, 2.5, 3.5, 5.0, 7.0, 10.0]
for j, r_inst in enumerate(instar_radii[:8]):
    P_inst = 20e6 * k_t * ((r_inst * 1e-3)**0.8) / (r_inst * 1e-3) / 1000
    if P_inst < 250:
        ax3.plot(r_inst, P_inst, 's', color='#fbbf24', markersize=7, zorder=5)
        if j % 2 == 0:
            ax3.annotate(f'Instar {j+1}', (r_inst, P_inst), textcoords="offset points",
                        xytext=(10, 5), fontsize=9, color='#fbbf24')

ax3.set_xlabel('Body Radius (mm)', fontsize=12, color='white')
ax3.set_ylabel('Required Burst Pressure (kPa)', fontsize=12, color='white')
ax3.set_title('Molting Pressure vs Body Size', fontsize=14, fontweight='bold', color='#c4b5fd')
ax3.legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#7c3aed', labelcolor='white')
ax3.set_facecolor('#0a0a1a')
ax3.tick_params(colors='white')
ax3.set_ylim(0, 250)
for spine in ax3.spines.values():
    spine.set_color('#7c3aed')

# ========== Panel 4: Sexual Dimorphism - Rensch's Rule ==========
ax4 = axes[1, 1]

# Spider sexual dimorphism data (log mass female vs male)
np.random.seed(123)

# Simulate data: in spiders, Rensch's rule is REVERSED
# Larger species show MORE female-biased dimorphism
n_species = 60
log_female_mass = np.random.uniform(-1, 3, n_species)

# Reversed Rensch: slope < 1 on log-log (male size doesn't keep up)
slope = 0.75  # reversed Rensch
intercept = -0.3
noise = np.random.normal(0, 0.2, n_species)
log_male_mass = slope * log_female_mass + intercept + noise

# Some extreme dimorphism species (Nephila-like)
extreme_idx = log_female_mass > 2.0
log_male_mass[extreme_idx] -= 0.5  # males even smaller relative

# Color by dimorphism ratio
dimorphism = log_female_mass - log_male_mass
colors_scat = plt.cm.viridis((dimorphism - dimorphism.min()) / (dimorphism.max() - dimorphism.min()))

ax4.scatter(log_female_mass, log_male_mass, c=dimorphism, cmap='viridis',
           s=50, alpha=0.8, edgecolors='white', linewidth=0.3)

# Regression line
fit = np.polyfit(log_female_mass, log_male_mass, 1)
x_fit = np.linspace(-1, 3.5, 100)
ax4.plot(x_fit, np.polyval(fit, x_fit), color='#f472b6', linewidth=2,
         label=f'Slope = {fit[0]:.2f} (reversed Rensch)')

# Isometry line
ax4.plot(x_fit, x_fit, color='#4c1d95', linewidth=1.5, linestyle='--', label='Isometry (slope = 1)')

# Annotate Nephila
ax4.annotate('Nephila\n(100x dimorphism)', xy=(2.8, 0.5), fontsize=10,
            color='#fbbf24', fontweight='bold',
            arrowprops=dict(arrowstyle='->', color='#fbbf24'))

cb = plt.colorbar(ax4.collections[0], ax=ax4, pad=0.02)
cb.set_label('Female/Male size ratio (log)', color='white', fontsize=10)
cb.ax.tick_params(colors='white')

ax4.set_xlabel('log(Female Mass, mg)', fontsize=12, color='white')
ax4.set_ylabel('log(Male Mass, mg)', fontsize=12, color='white')
ax4.set_title("Sexual Dimorphism: Reversed Rensch's Rule", fontsize=14, fontweight='bold', color='#c4b5fd')
ax4.legend(fontsize=9, facecolor='#1a1a2e', edgecolor='#7c3aed', labelcolor='white')
ax4.set_facecolor('#0a0a1a')
ax4.tick_params(colors='white')
for spine in ax4.spines.values():
    spine.set_color('#7c3aed')

fig.patch.set_facecolor('#0a0a1a')
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a1a')
plt.close()
print("Reproduction & Development Analysis Complete")
print(f"Optimal mating duration (W_future=0): {t[np.argmax(P_max*(1-np.exp(-alpha*t))*F_current)]:.1f} min (max = sacrifice)")
print(f"Optimal mating duration (W_future=200): {t[np.argmax(P_max*(1-np.exp(-alpha*t))*F_current + np.exp(-0.15*t)*200)]:.1f} min")
print(f"Burst pressure at r=2mm: {20e6 * k_t * ((2e-3)**0.8) / 2e-3 / 1000:.1f} kPa")
print(f"Rensch slope: {fit[0]:.2f} (< 1 = reversed)")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

7.5 Extreme Sexual Dimorphism

Spiders exhibit some of the most extreme sexual size dimorphism (SSD) in the animal kingdom. In Nephila pilipes, females can reach 40 mm body length and weigh over 1 g, while males are only 5–6 mm and weigh ~10 mg — a 100-fold mass difference. This is female-biasedSSD, the reverse of the typical vertebrate pattern.

Rensch's Rule — and Its Reversal

Rensch's rule states that in male-biased SSD taxa, larger species show greater dimorphism (the allometric slope of male size on female size is > 1 on log-log scale). In spiders, this is reversed: larger species show greaterfemale-biased dimorphism, with allometric slopes of 0.6–0.8.

Allometric Scaling of SSD

\[ \log(M_{\text{male}}) = a + b \cdot \log(M_{\text{female}}), \qquad b < 1 \text{ (reversed Rensch)} \]

Empirical values: \(b \approx 0.73\) across orb-weaving spiders (Hormiga et al. 2000). Two hypotheses: (1) fecundity selection drives female gigantism (egg production scales as\(M^{0.7\text{--}1.0}\)), and (2) selection favors small male size for faster maturation and greater mate-searching agility (the "gravity hypothesis" — small males climb faster on vertical webs).

The Gravity Hypothesis

Moya-Laraño et al. (2002) proposed that small males are selectively favored because they can climb faster on vertical surfaces. For a spider of mass \(m\) climbing a vertical web, the metabolic cost of transport scales as:

\[ v_{\text{climb}} \propto m^{-1/3}, \qquad \text{COT} = \frac{mg}{v \cdot \eta} \]

Smaller males climb faster (\(v \propto m^{-1/3}\)) and arrive at females sooner, gaining a first-male sperm precedence advantage. In orb weavers where scramble competition determines mating success, this creates strong directional selection for male miniaturization.

References

Welke, K.W. & Schneider, J.M. (2012). Sexual cannibalism benefits offspring survival. Animal Behaviour, 83(1), 201–207.
Andrade, M.C.B. (1996). Sexual selection for male sacrifice in the Australian redback spider. Science, 271(5245), 70–72.
Foelix, R.F. (2011). Biology of Spiders (3rd ed.). Oxford University Press.
Hormiga, G., Scharff, N. & Coddington, J.A. (2000). The phylogenetic basis of sexual size dimorphism in orb-weaving spiders. Systematic Biology, 49(3), 435–462.
Moya-Laraño, J., Halaj, J. & Wise, D.H. (2002). Climbing to reach females: Romeo should be small. Evolution, 56(2), 420–425.
Blackledge, T.A. & Hayashi, C.Y. (2006). Silken toolkits: biomechanics of silk fibers spun by the orb web spider Argiope argentata. Journal of Experimental Biology, 209(13), 2452–2461.
Kuntner, M., Coddington, J.A. & Hormiga, G. (2008). Phylogeny of extant nephilid orb-weaving spiders: testing morphological and ethological homologies. Cladistics, 24(2), 147–217.
Schwager, E.E. et al. (2017). The house spider genome reveals an ancient whole-genome duplication during arachnid evolution. BMC Biology, 15, 62.
Herberstein, M.E. (ed.) (2011). Spider Behaviour: Flexibility and Versatility. Cambridge University Press.
Schneider, J.M. & Fromhage, L. (2010). Monogynous mating strategies in spiders. Animal Behaviour and Cognition, 71, 441–456.

M6. Vision & Color M8. Evolution & Diversity