Module 1: Plant Chemical Ecology

1.1 Allelopathy: Chemical Warfare Between Plants

Allelopathy — the chemical inhibition of one plant by another — was first rigorously documented by Molisch (1937) and has since been identified in hundreds of species. Allelopathic compounds are released via root exudation, leaf leachate, volatile emission, and decomposition of residues. They can persist in soil for weeks to months, creating zones of inhibition around the producing plant.

Juglone: A Case Study in Allelopathic Biochemistry

Juglone (5-hydroxy-1,4-naphthoquinone) is produced by black walnut (Juglans nigra) and related species. It is synthesized via the shikimate pathway and released primarily as the non-toxic glycoside hydrojuglone, which is oxidized in soil to the toxic quinone form. Juglone's mechanism of action involves:

• Inhibition of NADH dehydrogenase (Complex I) in mitochondria, blocking electron transport and ATP synthesis
• Redox cycling: juglone accepts electrons to form a semiquinone radical, which transfers electrons to O₂, generating superoxide (\(\text{O}_2^{-\cdot}\)) and causing oxidative stress
• Root growth inhibition through disruption of auxin-mediated cell elongation

Sorgoleone, produced by Sorghum bicolor, is a p-benzoquinone that inhibits photosystem II (PSII) electron transport by competing with plastoquinone for the Q_B binding site on the D1 protein. Its inhibition constant:

making sorgoleone more potent than the commercial herbicide atrazine at the same target.

Diffusion-Reaction Model for Allelopathic Zones

The spatial extent of allelopathic inhibition depends on the balance between chemical diffusion through soil, degradation by soil microbes, and production rate. The governing partial differential equation for concentration \(c(\mathbf{x}, t)\) is:

where \(D\) is the effective diffusion coefficient in soil (typically\(10^{-6}\text{-}10^{-5}\) cm²/s for small molecules in saturated soil),\(\lambda\) is the first-order degradation rate constant, and\(S(\mathbf{x})\) is the source term (production rate at position\(\mathbf{x}\)).

For a point source in radial symmetry at steady state (\(\partial c/\partial t = 0\)), the equation becomes:

This is a modified Bessel equation with solution:

where \(K_0\) is the modified Bessel function of the second kind. The characteristic decay length — the effective allelopathic radius — is:

For juglone in typical forest soil (\(D \approx 5 \times 10^{-6}\) cm²/s,\(\lambda \approx 10^{-7}\) s⁻¹), the allelopathic radius is:

This matches field observations that many plant species (tomato, alfalfa, blueberry) fail to thrive within approximately 7-15 m of mature black walnut trees. The actual inhibition zone depends on root architecture, soil hydrology, and the threshold concentration for biological effect.

1.2 Volatile Organic Compounds (VOCs) as Ecological Signals

Plants emit enormous quantities of volatile organic compounds (VOCs) — an estimated 1000 Tg C/year globally, equivalent to ~1.5% of total terrestrial GPP. The major biogenic VOCs include isoprene (C5), monoterpenes (C10), sesquiterpenes (C15), and various oxygenated compounds (methanol, acetone, acetaldehyde).

Isoprene: Thermotolerance and Beyond

Isoprene (2-methyl-1,3-butadiene, C₅H₈) is the most abundantly emitted biogenic VOC, produced by isoprene synthase in chloroplasts from dimethylallyl diphosphate (DMAPP) via the methylerythritol phosphate (MEP) pathway. Its primary function appears to be thermotolerance: isoprene partitions into thylakoid membranes, stabilizing them against heat-induced phase transitions.

The energy cost is substantial: 6 carbon equivalents (3 CO₂ + 20 ATP + 14 NADPH) per isoprene molecule, representing 2-10% of recently fixed carbon. The benefit is measurable: isoprene-emitting leaves survive temperatures 5-8 C higher than non-emitting leaves from the same species.

Tritrophic Signaling: Plants Recruit Bodyguards

When herbivores damage plant tissues, the plant releases a distinct blend of herbivore-induced plant volatiles (HIPVs) — typically a mixture of green leaf volatiles (C6 aldehydes and alcohols from lipoxygenase pathway), monoterpenes, sesquiterpenes, and the "universal" herbivore signal methyl salicylate (MeSA). These HIPVs attract the natural enemies of the herbivore (predators, parasitoid wasps), creating a tritrophic signaling cascade:

Remarkably, the HIPV blend is herbivore-specific: maize damaged by Spodoptera littoralisreleases a different terpenoid profile than maize damaged by Mythimna separata, and parasitoid wasps can discriminate between them (Turlings et al. 1990).

Henry's Law: VOC Partitioning

The emission of VOCs from leaf interior to atmosphere is governed by Henry's law, which describes the equilibrium partitioning of a volatile compound between liquid (mesophyll cell water) and gas (substomatal cavity) phases:

where \(p_i\) is the partial pressure of compound \(i\) in the gas phase, \(H_i\) is the Henry's law constant, and \(x_i\) is the mole fraction in the liquid phase. The dimensionless Henry's law constant (air-water partition coefficient) is:

The temperature dependence of Henry's constant follows the van't Hoff equation:

Integration gives:

For isoprene, \(\Delta H_\text{sol} \approx -30\) kJ/mol, meaning Henry's constant increases sharply with temperature. At 25 C, \(K_H \approx 8\) for isoprene — strongly favoring gas phase partitioning. This thermodynamic basis explains the exponential increase in isoprene emission with temperature (independently of the enzymatic temperature response of isoprene synthase).

1.3 Mycorrhizal Networks ("Wood Wide Web")

Approximately 90% of land plants form mycorrhizal associations — mutualistic symbioses where fungal hyphae colonize plant roots and extend into the surrounding soil. The two dominant types are arbuscular mycorrhizae (AM, Glomeromycota, ~80% of plant species) and ectomycorrhizae(ECM, Basidiomycota and Ascomycota, dominant in temperate/boreal forests).

Common Mycorrhizal Networks (CMNs)

When a single fungal mycelium colonizes the roots of multiple trees, it creates a common mycorrhizal network (CMN) — an underground conduit for resource transfer. Suzanne Simard's groundbreaking experiments (1997) used isotopic tracers (\(^{13}\text{C}\) and \(^{14}\text{C}\)) to demonstrate bidirectional carbon transfer between paper birch (Betula papyrifera) and Douglas fir (Pseudotsuga menziesii) through shared ectomycorrhizal fungi (Rhizopogon spp.).

Key findings from Simard's research program include:

• Net carbon transfer from birch (in full sun) to fir (shaded), averaging 3-10% of net photosynthesis
• "Mother trees" (large, old, highly connected hub trees) provision seedlings through CMNs, especially kin
• Defense signaling: trees attacked by herbivores or pathogens release jasmonic acid and other signals through CMNs, priming defenses in connected neighbors
• Network topology: CMNs follow scale-free network architecture, with a few highly connected hubs and many peripheral nodes

Fick's Law for Carbon Flux Through Hyphae

Carbon transfer through a hyphal strand can be modeled using Fick's first law of diffusion. The flux \(J\) (mol m⁻² s⁻¹) through a hyphal conduit is:

where \(D_\text{eff}\) is the effective diffusion coefficient through hyphal cytoplasm and \(dC/dx\) is the concentration gradient of sugar (as trehalose or mannitol — the major carbon forms in fungal cytoplasm) along the hypha. For a cylindrical hypha of radius \(a\) and length \(L\), the total flux rate is:

For typical ectomycorrhizal hyphae (\(a \approx 3\) \(\mu\)m,\(L \approx 1\) m, \(D_\text{eff} \approx 5 \times 10^{-10}\) m²/s):

However, actual carbon flux through CMNs is orders of magnitude higher than passive diffusion predicts. This indicates active transport mechanisms: cytoplasmic streaming (motor protein-driven bulk flow at ~1-5 \(\mu\)m/s), vesicular transport, and peristaltic contractions of hyphal walls. The effective "diffusion coefficient" including active transport can be 10³-10⁴ times higher than molecular diffusion alone.

Source-Sink Dynamics

Carbon flows from source trees (high photosynthetic rate, excess carbohydrate) to sink trees (shaded, stressed, or young seedlings) along the concentration gradient. The steady-state carbon balance for a tree connected to \(n\) neighbors via CMN is:

where \(P_i\) is photosynthetic input, \(R_i\) is respiratory consumption, and the sum represents net carbon exchange with all connected neighbors. This creates a network diffusion equation that can be analyzed using graph Laplacian methods.

1.4 Phytochemical Diversity as Defense Strategy

A single plant species can produce hundreds to thousands of secondary metabolites. Nicotiana tabacum produces over 2,500 known compounds; the genus Solanumcontains species producing complex cocktails of steroidal glycoalkaloids, phenolics, and protease inhibitors. Why such extraordinary diversity?

The Screening Hypothesis

Jones & Firn (1991) proposed the screening hypothesis: because the probability that any single compound is bioactive against a given threat is low, plants benefit from producing a large diversity of compounds. If the probability that compound \(k\) is effective against herbivore \(j\) is\(p_k\), then the probability that at least one compound in a cocktail of\(n\) compounds is effective is:

For small \(\bar{p}\) and large \(n\), this approaches:

With \(\bar{p} = 0.01\) and \(n = 500\) compounds, the probability of having at least one effective defense is \(P = 1 - e^{-5} = 99.3\%\). The metabolic cost scales linearly with \(n\) while the benefit saturates, predicting an optimal diversity level.

Optimal Defense Theory (ODT)

McKey (1974) and Rhoades (1979) formalized optimal defense theory: plants should allocate defenses to maximize fitness, concentrating chemical defenses in tissues with the highest value (reproductive organs, young leaves) and the highest vulnerability. The fitness function is:

where \(h(D)\) is the fraction of tissue lost to herbivory (a decreasing function of defense \(D\)) and \(c(D)\) is the metabolic cost of producing defenses (an increasing function of \(D\)). Differentiating and setting \(dW/dD = 0\):

The optimality condition becomes:

The left side is the marginal proportional benefit of defense (proportional reduction in herbivory per unit defense), and the right side is the marginal proportional cost. At the optimum, these are equal. This predicts that:

• Young leaves (high value, high vulnerability) should be more defended than old leaves — confirmed by higher alkaloid concentrations in young leaves of Nicotiana
• Seeds and flowers should be heavily defended — confirmed by high concentrations of cardiac glycosides in milkweed flowers
• Roots in nutrient-poor soils (higher replacement cost) should be better defended — confirmed by higher tannin levels in heath plants
• Plants under high herbivore pressure should invest more in defense — confirmed by induction experiments with jasmonic acid

Constitutive vs Induced Defenses

Plants deploy both constitutive defenses (always present: lignin, trichomes, basal alkaloids) and induced defenses (produced only upon attack: protease inhibitors, pathogenesis-related proteins, phytoalexins). The jasmonate signaling pathway mediates most herbivore-induced defenses:

Induced defenses save metabolic resources when herbivores are absent, but incur a time lag (typically 12-48 hours). The optimal balance between constitutive and induced defenses depends on the frequency and predictability of herbivore attack: high-frequency, predictable attack favors constitutive defense; rare, unpredictable attack favors induced defense.

Allelopathic Zone Diagram

The black walnut tree releases juglone from roots, bark, and fallen leaves. The compound diffuses through soil, creating a concentration gradient that inhibits sensitive plants within the characteristic diffusion-degradation radius \(r^* = \sqrt{D/\lambda}\).

Python: Allelopathy, VOCs, Mycorrhizae & Defense

This four-panel simulation explores: (1) the 2D steady-state juglone concentration field from a point-source black walnut, using the Bessel-function solution of the diffusion-degradation equation; (2) VOC emission rates vs temperature for isoprene, monoterpenes, and sesquiterpenes using the Guenther algorithm; (3) a simulated mycorrhizal carbon transfer network showing Fick's law flux between trees of different carbon status; and (4) optimal defense theory showing the fitness-maximizing level of chemical defense investment.

References

1. Simard, S.W., Perry, D.A., Jones, M.D., Myrold, D.D., Durall, D.M. & Molina, R. (1997). Net transfer of carbon between ectomycorrhizal tree species in the field. Nature, 388(6642), 579-582.
2. Turlings, T.C., Tumlinson, J.H. & Lewis, W.J. (1990). Exploitation of herbivore-induced plant odors by host-seeking parasitic wasps. Science, 250(4985), 1251-1253.
3. Guenther, A., Hewitt, C.N., Erickson, D., et al. (1995). A global model of natural volatile organic compound emissions. Journal of Geophysical Research, 100(D5), 8873-8892.
4. Rice, E.L. (1984). Allelopathy. 2nd ed. Academic Press, Orlando.
5. Jones, C.G. & Firn, R.D. (1991). On the evolution of plant secondary chemical diversity. Philosophical Transactions of the Royal Society B, 333(1267), 273-280.
6. McKey, D. (1974). Adaptive patterns of alkaloid physiology. American Naturalist, 108(961), 305-320.
7. Rhoades, D.F. (1979). Evolution of plant chemical defense against herbivores. In: Rosenthal, G.A. & Janzen, D.H. (eds), Herbivores: Their Interaction with Secondary Plant Metabolites, pp. 3-54. Academic Press.
8. Molisch, H. (1937). Der Einfluss einer Pflanze auf die andere — Allelopathie. Gustav Fischer, Jena.
9. Sharkey, T.D. & Yeh, S. (2001). Isoprene emission from plants. Annual Review of Plant Physiology and Plant Molecular Biology, 52, 407-436.