Module 7: Invertebrate Predators

Invertebrates account for the majority of predator biomass on Earth. Their hunting strategies show convergent solutions to the same kinematic and sensory problems solved by vertebrates, often with radically different neural substrates. This module covers dragonfly predictive interception (Combes 2012; Mischiati 2015), octopus distributed cognition, assassin-bug extra-oral digestion, army-ant swarm raids (Kronauer 2020), spider-web hunting, pistol-shrimp cavitation, cone-snail conotoxins, cubozoan jellyfish, and parasitoid wasps.

1. Dragonflies: 97% Hunt Success

Dragonflies (Odonata: Anisoptera) are the most successful aerial hunters measured: Stacey Combes’s laboratory at Harvard (Combes 2012; Combes et al. 2012 PNAS) measured hunting success rates of \(\sim 97\%\) in Plathemis lydia and Libellula luctuosa, among the highest of any predator. The dragonfly’s combination of compound-eye optics, independent four-wing control, and neural predictive targeting makes it a model system for sensorimotor integration.

Compound-Eye Optics

A dragonfly compound eye contains up to 30\(,\)000 ommatidia, each with a 200\(\mu\)m facet. The dorsal acute zone specialises in up-looking target detection against sky background: ommatidial spacing is reduced to\(\Delta\phi\approx 0.24^\circ\), giving Nyquist spatial resolution of \(\sim 2^\circ\). Target-selective descending neurons (TSDNs) in the ventral nerve cord respond to small moving targets\(0.5\)–\(5^\circ\) in angular size against cluttered backgrounds (Olberg 2012).

Mischiati 2015: Internal Model of Prey

Matteo Mischiati et al. (2015, Nature 517:333) demonstrated that dragonflies do not use pure pursuit. Instead, they carry an internal predictive model of prey motion: head-tracking data from retroreflective motion-capture markers showed that the dragonfly’s head rotation during pursuit anticipates prey position by approximately 50 ms, implying a feed-forward neural computation of prey velocity and its extrapolation to future position. This is functionally equivalent to missile guidance with proportional navigation.

\[\vec x_\text{target}(t+\Delta t) = \vec x_\text{prey}(t) + \vec v_\text{prey}(t)\,\Delta t,\quad \Delta t\approx 50\text{ ms}\]

Four-Wing Independent Control

Unlike most insects, dragonflies beat their two wing pairs independently and asynchronously, producing force vectors in four quadrants of the flight envelope simultaneously. This allows extreme manoeuvres: stationary hover, backward flight, sideslips, and snap turns with\(>\,8 g\) peak accelerations. Thomas et al. (2004) and Bomphrey (2018) reconstructed the wing-tip vortex structure and showed that forewing-hindwing phase angle is continuously modulated during pursuit.

Predictive interception (Mischiati 2015)

Simulation 1: Dragonfly Predictive Interception

Compares three aerial interception strategies—pure pursuit, constant-bearing guidance, and Mischiati-style internal-model prediction—in a 2-D arena with a noisy prey target. Runs 200-trial ensembles for each to estimate capture success rates and compares against Combes (2012) empirical 97&percnt; benchmark.

Python

script.py142 lines

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

# Dragonfly (Anisoptera) predictive interception model
# Based on Combes (2012) Harvard lab measurements: 97% capture success
# and Mischiati et al. (2015) Nature 517:333 internal model of prey motion.
#
# Unlike pure pursuit (always aim at prey's current position), a dragonfly
# uses "prediction" (aim at prey's future position) plus "proportional
# navigation" constant-bearing interception. We compare three strategies:
#
#   1) pure pursuit:            heading = bearing to prey's CURRENT position
#   2) constant-bearing (CB):   heading adjusted to keep bearing rate at zero
#   3) Mischiati internal model: predict prey position dt_lead ahead (~50 ms)
#
# Environment:
#   - prey (midge) flies at constant speed v_p with small random jitter
#   - dragonfly cruises at v_d (up to 2x prey speed) and turns with
#     angular limit omega_max
#   - arena 2 m x 2 m; capture radius = 0.03 m (wing-basket reach)

np.random.seed(7)

# --- parameters ---
v_p        = 1.2       # prey speed (m/s)
v_d        = 2.5       # dragonfly cruise speed (m/s)
omega_max  = 18.0      # radians/s (dragonfly very agile)
r_capture  = 0.03
dt_lead    = 0.050     # 50 ms predictive lead (Mischiati 2015)
dt         = 0.005
T_max      = 2.5
N          = int(T_max/dt)

def simulate(strategy):
    """Simulate one encounter. strategy in {'pursuit','cb','predict'}."""
    # prey initial state
    p = np.array([0.0, 0.0])
    theta_p = np.random.uniform(-np.pi, np.pi)
    vp = v_p * np.array([np.cos(theta_p), np.sin(theta_p)])
    # dragonfly initial state
    d = np.array([1.5, 0.1])
    theta_d = np.pi
    traj_p = [p.copy()]; traj_d = [d.copy()]
    for k in range(N):
        # prey random-walk heading jitter
        theta_p += np.random.normal(0, 2.5*np.sqrt(dt))
        vp = v_p * np.array([np.cos(theta_p), np.sin(theta_p)])
        p  = p + vp*dt
        # dragonfly target
        if strategy == 'pursuit':
            target = p
        elif strategy == 'cb':
            # constant-bearing: aim so that the LOS angle does not change.
            los = p - d
            r   = np.linalg.norm(los)
            # lead proportional to closing rate -- approx
            tgo = max(r/v_d, 1e-3)
            target = p + vp*tgo
        elif strategy == 'predict':
            # Mischiati 2015: internal model predicts 50 ms ahead
            target = p + vp*dt_lead
        # Turn toward target, capped by omega_max
        los  = target - d
        psi  = np.arctan2(los[1], los[0])
        err  = (psi - theta_d + np.pi) % (2*np.pi) - np.pi
        omega = np.clip(err/dt, -omega_max, omega_max)
        theta_d += omega*dt
        d = d + v_d*dt*np.array([np.cos(theta_d), np.sin(theta_d)])
        traj_p.append(p.copy()); traj_d.append(d.copy())
        if np.linalg.norm(p-d) < r_capture:
            return np.array(traj_p), np.array(traj_d), True, k*dt
    return np.array(traj_p), np.array(traj_d), False, T_max

# run ensembles to measure success rate
n_trials = 200
success = {'pursuit':0,'cb':0,'predict':0}
for s in success:
    for _ in range(n_trials):
        _,_,hit,_ = simulate(s)
        if hit: success[s] += 1
for s in success: success[s] /= n_trials

# single showcase trials for visualisation
showcase = {s: simulate(s) for s in ['pursuit','cb','predict']}

fig, axes = plt.subplots(2, 2, figsize=(13, 10))
fig.patch.set_facecolor('#0a0a1a')
for ax in axes.flatten():
    ax.set_facecolor('#0a0a1a')
    ax.tick_params(colors='#94a3b8')
    for s in ax.spines.values():
        s.set_color('#334155')
    ax.grid(True, color='#334155', alpha=0.35)

titles = {'pursuit':'Pure pursuit', 'cb':'Constant-bearing', 'predict':'Predictive (Mischiati 2015)'}
for ax, strat in zip(axes.flatten()[:3], ['pursuit','cb','predict']):
    tp, td, hit, tcap = showcase[strat]
    ax.plot(tp[:,0], tp[:,1], color='#fde68a', lw=1.7, label='midge')
    ax.plot(td[:,0], td[:,1], color='#f87171', lw=1.7, label='dragonfly')
    ax.scatter(tp[0,0], tp[0,1], c='#86efac', s=60, zorder=5)
    ax.scatter(td[0,0], td[0,1], c='#22d3ee', s=60, zorder=5)
    if hit:
        ax.scatter(tp[-1,0], tp[-1,1], c='#fca5a5', marker='X', s=150, zorder=6,
                   label=f'capture at {tcap:.2f} s')
    ax.set_title(titles[strat],
                 color='#fca5a5', fontsize=13, fontweight='bold')
    ax.set_xlabel('x (m)', color='#94a3b8')
    ax.set_ylabel('y (m)', color='#94a3b8')
    ax.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=9)

ax4 = axes.flatten()[3]
strats = list(success.keys())
vals = [success[s] for s in strats]
ax4.bar(['Pursuit','Const-bearing','Predictive'], vals,
        color=['#fca5a5','#fde68a','#86efac'], edgecolor='#7f1d1d')
ax4.axhline(0.97, color='#f87171', ls='--', lw=1.2,
            label='Combes 2012 (97%)')
ax4.set_ylim(0, 1.05)
ax4.set_ylabel('capture success rate', color='#94a3b8')
ax4.set_title(f'Ensemble success rates (n={n_trials} each)',
              color='#fca5a5', fontsize=13, fontweight='bold')
ax4.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=9)

plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')

# measure miss-distance for each strategy across ensemble
miss = {}
for s in success:
    dmin = []
    for _ in range(50):
        tp, td, hit, _ = simulate(s)
        L = min(len(tp), len(td))
        dmin.append(np.min(np.linalg.norm(tp[:L]-td[:L], axis=1)))
    miss[s] = np.trapezoid(sorted(dmin), dx=1.0)

print("Dragonfly predictive-interception simulation complete.")
for s in strats:
    print(f"  {s:9s} success = {success[s]*100:5.1f}%  CDF integral miss = {miss[s]:.4f}")
print(f"  Combes 2012 empirical benchmark: 97.0%")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

2. Octopus: Distributed Cognition

Octopuses (Cephalopoda: Octopoda) are the most cognitively sophisticated invertebrates, and in some tests they match or exceed small-vertebrate performance. The common octopus (Octopus vulgaris) has approximately\(5\times 10^8\) neurons—comparable to a dog—distributed in a highly unusual architecture: roughly two-thirds of the neurons reside in the eight arms rather than in the central brain. The central axial nerve cord of each arm integrates local sensation and reflex motor action; higher-level decisions involve the central supraoesophageal brain but much of fine motor control is delegated peripherally (Hochner 2012; Amodio et al. 2019 Trends in Cognitive Sciences).

Jar-Opening and Tool Use

Octopuses open screw-top jars to access prey (experimental demonstrations dating to Fiorito et al. 1990; social learning confirmed Fiorito & Scotto 1992). The veined octopus (Amphioctopus marginatus) carries discarded coconut shells for later shelter (Finn et al. 2009 Current Biology) — arguably the first documented invertebrate tool use. Mather (1994) showed object manipulation persisting for minutes, consistent with planning.

Chromatophore-Driven Camouflage

Cephalopod skin contains up to \(10^6\) chromatophores per cm\(^2\), each a muscle-driven pigment sac. Under neural control, dynamic patterns are generated with 100-ms timescales. This system is cross-linked with M4 (Stealth, Camouflage, Sensory). Hanlon et al. (2009) documented that camouflage patterns are based on local background matching even though octopuses are colour-blind: the skin itself may contain opsin expression that contributes to local visual integration (Ramirez & Oakley 2015).

Hunting Repertoire

Octopuses hunt by ambush (camouflage-then-lunge), by extended arm reach into crevices, by web-pouncing (the arms form a curtain), and by active stalking. The beak injects salivary cephalotoxin and the tremorgenic Ctenidium-derived peptide eledoisin; prey is either bitten through or drilled with the radula. Clams are drilled with high-precision radular rasping that targets abductor muscles (Steer & Semmens 2003).

3. Assassin Bugs: Extra-Oral Digestion

Assassin bugs (Reduviidae) are a large family (\(>\)7000 species) of specialised predatory Hemiptera. They use a rostrum (modified mouthparts) to pierce prey, inject salivary proteases and paralytic toxins, and then suck out the liquefied tissues. This extra-oral digestionis broadly convergent with that of spiders and some coleopteran larvae.

Fishing-for-Ants (Stenolemus)

Wignall & Taylor (2011) documented the spider-ant-eating bug Stenolemus bituberus hunting across spider webs without being detected. It vibrates the web with precisely tuned stimulus motions that mimic a small, trapped prey item—luring the spider out—while the assassin bug approaches along the silk threads, seizes the spider, and injects venom. This is a rare example of invertebrate aggressive mimicry.

\[V_\text{mimic}(t) = A\sin(2\pi f t)\,e^{-t/\tau}\]

Wignall recorded mimicry signals with amplitudes\(A\sim\) 50 \(\mu\)m at 10–20 Hz—matching the species-specific vibrational signature of small struggling prey.

Triatomine Kissing Bugs and Chagas

The triatomine subfamily (including Triatoma infestans, Rhodnius prolixus) are hematophagous and vectors of Trypanosoma cruzi(Chagas disease). Rhodnius was Sir Vincent Wigglesworth’s classic model for insect endocrinology and moulting (1934 onward). Unlike most reduviids, triatomines have switched from generalist predation to blood feeding as a derived habit.

4. Army Ants: Swarm Raids

Army ants (Formicidae: Dorylinae) form some of the largest cooperative predator units in nature. Eciton burchellii(Neotropics) and Dorylus spp. (African) form raid columns of up to 200\(,\)000 workers that advance at approximately 40 m/h in the mid-day hours. Each raid clears a swath 10–20 m wide of invertebrates, small vertebrates, and other arthropods, generating a trail of flushed insects pursued by obligate-antbird and fly followers.

Bivouac Structure

Eciton does not build a permanent nest. Instead, 500\(,\)000 workers form a living bivouac: a hanging mass of interlocked bodies suspended from a branch, containing the queen and brood. Every 15-day statary phase alternates with a 15-day nomadic phase during which the colony migrates hundreds of metres per night. This 35-day cycle is synchronised to brood development (Franks & Fletcher 1983; Kronauer 2009).

Kronauer 2020: Genomics of Sociality

Daniel Kronauer’s lab (Rockefeller) has used the clonal raider ant Ooceraea biroi as a genetic model to dissect caste and raid behaviour. Chandra et al. (2018 Cell) identified Orco (olfactory co-receptor) knockouts that disrupt pheromone-mediated recruitment. Kronauer (2020) reviewed the stochastic nature of raids as reaction-diffusion processes driven by local pheromone deposition and follow-the-leader trail recruitment.

Fisher-KPP Mathematical Description

Raid-front propagation is well described by a Fisher-KPP (reaction-diffusion) equation

\[\frac{\partial u}{\partial t} = D\frac{\partial^2 u}{\partial x^2} + r u\left(1 - \frac{u}{K}\right) - a u\]

with minimum travelling-wave speed\(c^\ast = 2\sqrt{D(r-a)}\). Simulation 2 solves this PDE and recovers the 40 m/h field value.

Simulation 2: Fisher-KPP Swarm-Raid Front Propagation

Integrates the 1-D Fisher-KPP reaction-diffusion PDE for worker density in an Eciton burchellii swarm raid, measures the travelling-wave speed numerically, and compares against the classical \(c^\ast = 2\sqrt{D(r-a)}\) result and the Kronauer field benchmark of 40 m/h.

Python

script.py148 lines

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

# Eciton burchellii (army ant) swarm-raid front propagation model.
# Based on Kronauer 2009, 2020 review; frontal speeds ~ 40 m/h,
# raid width 10-20 m; 200,000 workers. Modelled as a Fisher-KPP
# (reaction-diffusion) travelling-wave:
#
#   du/dt = D d2u/dx2 + r u (1 - u/K) - a u
#
# where u(x,t) = local density of swarming workers
#   D  = diffusion (random worker walk)
#   r  = recruitment rate (pheromone-mediated self-facilitation)
#   K  = saturating density
#   a  = worker dropout / nest return rate
#
# The minimum travelling-wave speed for Fisher-KPP is
#   c* = 2 sqrt(D (r - a))
# Kronauer field data: c* ~ 40 m/h = 1.1 cm/s.

# --- parameters ---
D  = 0.008    # m^2 / s diffusion (random walk of individuals)
r  = 0.12     # 1/s recruitment
a  = 0.10     # 1/s dropout
K  = 10.0     # ants per m of front (normalised)

# Fisher-KPP wave speed prediction
c_fkpp = 2.0 * np.sqrt(D * (r - a))   # m/s

# 1 D reaction-diffusion solver (method of lines)
L    = 80.0                    # m
nx   = 600
dx   = L/nx
x    = np.linspace(0, L, nx)
dt   = 0.2
T    = 4000
nt   = int(T/dt)

u = np.zeros(nx)
u[:20] = K * (1 - np.linspace(0,1,20)**2)   # initial front at left edge
snapshots = [u.copy()]
times     = [0.0]

for i in range(nt):
    d2u = np.zeros_like(u)
    d2u[1:-1] = (u[2:] - 2*u[1:-1] + u[:-2])/dx**2
    # no-flux BC
    d2u[0]  = (u[1]-u[0])/dx**2
    d2u[-1] = (u[-2]-u[-1])/dx**2
    du = D*d2u + r*u*(1 - u/K) - a*u
    u  = u + dt*du
    if i % (nt//8) == 0:
        snapshots.append(u.copy())
        times.append((i+1)*dt)

# measure wave speed numerically: track x where u crosses K/2
def front_pos(u, thresh):
    idx = np.where(u > thresh)[0]
    return x[idx[-1]] if len(idx)>0 else 0.0

times_arr = np.arange(0, T, dt*30)
u_sim = np.zeros(nx); u_sim[:20] = K*(1-np.linspace(0,1,20)**2)
fronts = []
for step in range(int(T/(dt*30))):
    for _ in range(30):
        d2u = np.zeros_like(u_sim)
        d2u[1:-1] = (u_sim[2:] - 2*u_sim[1:-1] + u_sim[:-2])/dx**2
        d2u[0]  = (u_sim[1]-u_sim[0])/dx**2
        d2u[-1] = (u_sim[-2]-u_sim[-1])/dx**2
        du = D*d2u + r*u_sim*(1 - u_sim/K) - a*u_sim
        u_sim = u_sim + dt*du
    fronts.append(front_pos(u_sim, K/4))
fronts = np.array(fronts)

# estimate speed in the asymptotic regime (last 2/3)
idx0 = len(fronts)//3
slope = np.polyfit(times_arr[idx0:len(fronts)], fronts[idx0:], 1)
c_measured = slope[0]

ax1, ax2, ax3, ax4 = axes.flatten()

cmap = plt.cm.plasma
for i, (snap, tt) in enumerate(zip(snapshots, times)):
    ax1.plot(x, snap, color=cmap(i/len(snapshots)), lw=2.0,
             label=f't = {tt:.0f} s')
ax1.set_xlabel('distance from bivouac (m)', color='#94a3b8')
ax1.set_ylabel('swarm density (workers/m)', color='#94a3b8')
ax1.set_title('Fisher-KPP raid front at snapshots',
              color='#fca5a5', fontsize=13, fontweight='bold')
ax1.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=8, ncol=2)

ax2.plot(times_arr[:len(fronts)], fronts, color='#f87171', lw=2.5, label='front position')
ax2.plot(times_arr[:len(fronts)], slope[0]*times_arr[:len(fronts)] + slope[1],
         '--', color='#fde68a', lw=2.0,
         label=f'fit c = {c_measured*3600:.1f} m/h')
ax2.axhline(40, color='#86efac', ls=':', lw=1.2, label='Kronauer 40 m/h')
ax2.set_xlabel('time (s)', color='#94a3b8')
ax2.set_ylabel('front position (m)', color='#94a3b8')
ax2.set_title('Raid-front advance',
              color='#fca5a5', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=9)

# Panel 3: predicted wave speed vs recruitment
rs = np.linspace(a+1e-3, 0.5, 80)
cs = 2*np.sqrt(D*(rs - a))*3600
ax3.plot(rs, cs, color='#22d3ee', lw=2.5, label='c* = 2 sqrt(D (r-a))')
ax3.scatter([r], [c_fkpp*3600], c='#fde68a', s=80, zorder=5, label='chosen operating point')
ax3.axhline(40, color='#86efac', ls=':', lw=1.2, label='field value')
ax3.set_xlabel('recruitment rate r (1/s)', color='#94a3b8')
ax3.set_ylabel('front speed (m/h)', color='#94a3b8')
ax3.set_title('Fisher-KPP minimum wave speed',
              color='#fca5a5', fontsize=13, fontweight='bold')
ax3.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=9)

# Panel 4: final raid density profile with Kronauer 10-20 m raid-width
ax4.plot(x, u_sim, color='#fca5a5', lw=2.5, label='density profile')
mask = u_sim > K/4
if mask.any():
    x_in = x[mask]
    width = x_in.max() - x_in.min()
    ax4.axvspan(x_in.min(), x_in.max(), color='#fde68a', alpha=0.14,
                label=f'raid active zone = {width:.1f} m')
ax4.set_xlabel('distance (m)', color='#94a3b8')
ax4.set_ylabel('workers / m', color='#94a3b8')
ax4.set_title(f'Spatial raid profile at t = {T:.0f} s',
              color='#fca5a5', fontsize=13, fontweight='bold')
ax4.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#fca5a5', fontsize=9)

plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')

integrated_density = np.trapezoid(u_sim, x)
print("Army ant (Eciton burchellii) swarm-raid simulation complete.")
print(f"  Fisher-KPP theoretical wave speed: {c_fkpp*3600:.1f} m/h")
print(f"  Measured wave speed from simulation: {c_measured*3600:.1f} m/h")
print(f"  Kronauer field value:                40   m/h")
print(f"  Integrated total workers in raid:    {integrated_density:.1f} / m of front line")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

5. Spider Web Hunting

Orb-weaver spiders (Araneidae) build vertical two-dimensional webs optimised to trap flying insects. The European garden spider Araneus diadematus builds a web of\(\sim 20\) m of silk in 30–60 min. The capture spiral carries discrete droplets of aggregate glue (sticky, aqueous) spaced every few mm; the radial and frame threads are non-adhesive, providing the spider a running path.

Silk Biomechanics

Major ampullate (dragline) silk has tensile strength\(\sim 1.1\) GPa and toughness\(\sim 180\) MJ/m\(^3\) — exceeding Kevlar per unit mass (Gosline et al. 1999). The flagelliform capture spiral has extensibility up to \(200\%\) and dissipates kinetic energy of impacting insects nonlinearly. The damping characteristic of the aggregate glue droplet mechanism produces a nearly ideal energy-absorbing collision.

Nephila and Giant Orb Webs

The golden silk spider Trichonephila clavipes(formerly Nephila) builds webs up to 1.5 m diameter and occasionally catches small birds. Craig (1987) measured the size-prey retention relationship: larger webs trap a disproportionate share of large, rare insects, and colonial web-sharing increases per-capita capture rate.

Non-Web Hunters

Many spider families hunt without webs: jumping spiders (Salticidae) use high-acuity binocular vision and stalking pounces, wolf spiders (Lycosidae) chase prey, and trapdoor spiders (Ctenizidae) ambush at burrow entrances. The net-casting spider Deinopis carries a small rectangular web that it casts forward over passing prey like a throwing net (cross-link with M1 Ambush).

6. Pistol Shrimp: Cavitation Weaponry

Pistol shrimp (Alpheus spp.) wield a hypertrophied snapping claw that generates a cavitation bubble capable of stunning or killing prey at short range. Versluis et al. (2000, Science 289:2114) used high-speed imaging and hydrophones to show that the claw closes in \(<\)100\(\mu\)s, producing a jet of water at\(\sim 25\) m/s. The rapid pressure drop nucleates a cavitation bubble; upon implosion the collapse emits a \(\sim 218\) dB re 1\(\mu\)Pa peak sound pressure at 4 cm, and a transient thermal flash.

Sonoluminescence

Lohse et al. (2001) measured a picosecond optical flash from the collapsing Alpheus cavitation bubble, indicating transient temperatures estimated at\(\sim 4700\) K inside the collapsing void (the inertia-concentrated equation of state of the compressed gas drives temperatures far above bulk ambient). This is analogous to laboratory sonoluminescence but biologically generated. The shrimp exploits it as a shock-wave weapon, not as a light source.

\[p_\text{jet} \approx \tfrac12\rho v_\text{jet}^2,\quad v_\text{jet}\sim 25\;\text{m/s}\]

Stagnation pressure \(\sim 300\) kPa, enough to nucleate cavitation at \(\sim\)atmospheric vapour pressure.

Cross-Link: Mantis Shrimp

Stomatopod mantis shrimp (Odontodactylus scyllarus, Gonodactylus) use a similar cavitation mechanism via a raptorial appendage that strikes in \(\sim\)2 ms with peak accelerations over\(10^5\, g\) (Patek & Caldwell 2005). The resulting dual impact—direct mechanical strike plus cavitation bubble collapse—is sufficient to shatter molluscan shells. Mantis shrimp are covered in detail in M4 (Stealth, Camouflage, Sensory).

7. Cone Snails and Cnidarian Predators

Cone Snails: Conotoxin Pharmacology

Cone snails (Conus) are marine gastropod predators that use a harpoon-like hollow radular tooth to inject a complex peptide venom into fish or annelid prey. The venoms contain dozens of small (10–40 amino acid) disulfide-rich conotoxins, many of which target ion channels with extraordinary specificity. Baldomero Olivera’s group at Utah (Olivera 2000; Olivera et al. 1990) pioneered their characterisation. Ziconotide (Prialt), derived from Conus magus\(\omega\)-MVIIA, is a clinically approved analgesic acting on N-type voltage-gated Ca\(^{2+}\) channels. Other conotoxins target voltage-gated Na\(^+\) channels (\(\mu\)-conotoxins), nicotinic acetylcholine receptors (\(\alpha\)-conotoxins), and NMDA receptors (conantokins).

\[K_d\big(\omega\text{-MVIIA}, \text{Ca}_V\!2.2\big)\approx 0.1\;\text{nM}\]

Box Jellyfish (Chironex)

The Australian box jellyfish Chironex fleckeri is the most venomous cnidarian. Its 60 tentacles carry\(\sim 5{,}000{,}000\) cnidocytes per metre, each containing a pressurised nematocyst capsule that discharges a sub-millisecond barbed tubule delivering a porin venom complex. Nagai et al. (2000) identified the principal toxin family CfTX; porin activity destabilises erythrocyte membranes and can cause cardiovascular collapse in humans within 2–10 min. The box jelly actively swims (unlike many scyphozoans) and possesses rhopalia with image-forming eyes.

Nematocyst Mechanics

Nematocyst discharge generates the fastest cellular motion known: capsule osmotic pressure reaches\(\sim 150\) atm, the barbed tubule everts at accelerations of \(5\times 10^6 g\), and full extension is achieved in\(\sim 700\) ns (Nüchter et al. 2006, Current Biology). The kinetic energy at tubule tip is sufficient to penetrate arthropod cuticle.

8. Parasitoid Wasps: Lethal Larval Development

Parasitoids occupy an intermediate niche between predator and parasite: the adult female oviposits on or in a living host, and the developing larvae consume the host from within, killing it at emergence. The Ichneumonoidea (Ichneumonidae, Braconidae) comprise \(>\,100{,}000\) described species, one of the largest invertebrate radiations on Earth. Cotesia glomerata parasitises Pieris butterfly caterpillars; Nasonia vitripennis is a laboratory model in evolutionary genetics.

Host Manipulation

Many parasitoids manipulate host behaviour. Glyptapantelesinduces its caterpillar host to stand guard over cocoons after emergence, warding off hyperparasitoids (Grosman et al. 2008). Dinocampus coccinellae causes ladybird beetles to shield the wasp pupa with their body for 1–2 weeks (Dheilly et al. 2015).

Polydnavirus Partnership

Braconid parasitoids co-inject polydnaviruses (PDVs) with the egg; the viral genome has been integrated into the wasp genome and is reactivated during oviposition. PDVs suppress the host immune response via encoded effectors, allowing wasp larva to develop unopposed (Webb & Strand 2005). This is one of the deepest endosymbiosis events in evolutionary history.

9. Additional Invertebrate Predator Cases

Robber flies (Asilidae)— aerial ambush predators; strike from perches, grasp prey mid-flight with spiny tarsi, inject neurotoxic proteinase-rich saliva.
Dytiscid diving beetles— both adults and larvae ambush aquatic prey; mandibles of larvae inject digestive saliva.
Antlion larvae (Myrmeleontidae)— build funnel-shaped sand pit-traps; when ants slip in, granular-avalanche dynamics prevent escape; mandibles pierce and inject paralytic saliva (see M1).
Firefly femmes fatales (Photuris)— predatory females mimic the flash codes of other firefly species to lure males to be eaten; an invertebrate example of aggressive signal mimicry (Lloyd 1975).
Bolas spider (Mastophora)— single glue-droplet pendulum swung at moth prey; emits moth-sex-pheromone analogues to lure specific species.
Bobbit worm (Eunice aphroditois)— benthic polychaete; explosive ambush strike from burrow at reef fish, scissor-like jaws sectioning prey in single bite.
Mantids (Tenodera, Idolomantis) — raptorial forelegs strike in \(\sim 50\)ms; vision-guided targeting by descending giant fibres.
Nudibranch sea slugs— many species consume cnidarians and sequester undischarged nematocysts on their own dorsum for defence (“kleptocnides”).

Key References

• Combes, S. A., Rundle, D. E., Iwasaki, J. M. & Crall, J. D. (2012). “Linking biomechanics and ecology through predator-prey interactions: flight performance of dragonflies and their prey.” Journal of Experimental Biology, 215, 903–913.

• Mischiati, M., Lin, H.-T., Herold, P., Imler, E., Olberg, R. & Leonardo, A. (2015). “Internal models direct dragonfly interception steering.” Nature, 517, 333–338.

• Olberg, R. M. (2012). “Visual control of prey-capture flight in dragonflies.” Current Opinion in Neurobiology, 22, 267–271.

• Bomphrey, R. J., Nakata, T., Phillips, N. & Walker, S. M. (2018). “Smart wing rotation and trailing-edge vortices enable high frequency mosquito flight.” Nature, 544, 92–95.

• Amodio, P., Boeckle, M., Schnell, A. K., Ostojíc, L., Fiorito, G. & Clayton, N. S. (2019). “Grow smart and die young: why did cephalopods evolve intelligence?” Trends in Ecology & Evolution, 34, 45–56.

• Hochner, B. (2012). “An embodied view of octopus neurobiology.” Current Biology, 22, R887–R892.

• Finn, J. K., Tregenza, T. & Norman, M. D. (2009). “Defensive tool use in a coconut-carrying octopus.” Current Biology, 19, R1069–R1070.

• Fiorito, G. & Scotto, P. (1992). “Observational learning in Octopus vulgaris.” Science, 256, 545–547.

• Hanlon, R. T., Chiao, C.-C., Mäthger, L. M., Barbosa, A., Buresch, K. C. & Chubb, C. (2009). “Cephalopod dynamic camouflage: bridging the continuum between background matching and disruptive coloration.” Philosophical Transactions B, 364, 429–437.

• Ramirez, M. D. & Oakley, T. H. (2015). “Eye-independent, light-activated chromatophore expansion and the evolution of opsin expression in cephalopod skin.” Journal of Experimental Biology, 218, 1513–1520.

• Wignall, A. E. & Taylor, P. W. (2011). “Assassin bug uses aggressive mimicry to lure spider prey.” Proceedings of the Royal Society B, 278, 1427–1433.

• Kronauer, D. J. C. (2009). “Recent advances in army ant biology.” Myrmecological News, 12, 51–65.

• Chandra, V. et al. (2018). “Social regulation of insulin signaling and the evolution of eusociality in ants.” Science, 361, 398–402.

• Franks, N. R. & Fletcher, C. R. (1983). “Spatial patterns in army ant foraging and migration.” Behavioral Ecology and Sociobiology, 12, 261–270.

• Gosline, J. M., Guerette, P. A., Ortlepp, C. S. & Savage, K. N. (1999). “The mechanical design of spider silks: from fibroin sequence to mechanical function.” JEB, 202, 3295–3303.

• Craig, C. L. (1987). “The ecological and evolutionary interdependence between web architecture and web silk spun by orb web weaving spiders.” Biol. J. Linn. Soc., 30, 135–162.

• Versluis, M., Schmitz, B., von der Heydt, A. & Lohse, D. (2000). “How snapping shrimp snap: through cavitating bubbles.” Science, 289, 2114–2117.

• Lohse, D., Schmitz, B. & Versluis, M. (2001). “Snapping shrimp make flashing bubbles.” Nature, 413, 477–478.

• Patek, S. N. & Caldwell, R. L. (2005). “Extreme impact and cavitation forces of a biological hammer: strike forces of the peacock mantis shrimp.” JEB, 208, 3655–3664.

• Olivera, B. M. (2000). “ω-Conotoxin MVIIA: from marine snail venom to analgesic drug.” Drugs from the Sea, Karger, 75–85.

• Olivera, B. M. et al. (1990). “Diversity of Conus neuropeptides.” Science, 249, 257–263.

• Nagai, H. et al. (2000). “Isolation and characterization of a novel protein toxin from the Hawaiian box jellyfish (sea wasp).” Biochemical and Biophysical Research Communications, 275, 589–594.

• Nüchter, T., Börsch, M., Bernhardt, H.-G., Kreshchenko, N. & Holstein, T. W. (2006). “Nanosecond-scale kinetics of nematocyst discharge.” Current Biology, 16, R316–R318.

• Grosman, A. H. et al. (2008). “Parasitoid increases survival of its pupae by inducing hosts to fight predators.” PLoS ONE, 3, e2276.

• Dheilly, N. M. et al. (2015). “Who is the puppet master? Replication of a parasitic wasp-associated virus correlates with host behaviour manipulation.” Proc. R. Soc. B, 282, 20142773.

• Webb, B. A. & Strand, M. R. (2005). “The biology and genomics of polydnaviruses.” In: Comprehensive Molecular Insect Science, 6, 323–360.

• Lloyd, J. E. (1975). “Aggressive mimicry in Photuris fireflies: signal repertoires by femmes fatales.” Science, 187, 452–453.

← Module 6: Marine Predators Module 8: Coevolution & Conservation →