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.

1.5 Terpenoid Diversity: The Largest Class of Natural Products

Terpenoids (isoprenoids) constitute the largest and most structurally diverse class of natural products, with over 80,000 known structures. Every terpenoid is built from the same simple 5-carbon precursor — the isoprene unit — yet the combinatorial explosion of chain length, cyclization patterns, and oxidative modifications produces an astronomical diversity of final structures.

The Universal C5 Building Block: IPP and DMAPP

All terpenoids are assembled from two isomeric C5 units: isopentenyl diphosphate (IPP) and its isomer dimethylallyl diphosphate (DMAPP). These are produced by two independent biosynthetic pathways:

• MVA (mevalonate) pathway — operates in the cytosol; starting from acetyl-CoA via HMG-CoA reductase (the target of statin drugs). Primarily supplies IPP/DMAPP for sesquiterpenes, triterpenes, and sterols
• MEP (methylerythritol phosphate) pathway — operates in plastids; starting from pyruvate and glyceraldehyde-3-phosphate via DXS (1-deoxy-D-xylulose 5-phosphate synthase). Supplies IPP/DMAPP for monoterpenes, diterpenes, carotenoids, and the phytol tail of chlorophyll

The two pathways are compartmentalized but not completely independent — limited exchange of IPP occurs between cytosol and plastid, allowing metabolic cross-talk.

Terpenoid Classification by Chain Length

Sequential head-to-tail condensation of IPP units by prenyltransferases builds longer prenyl chains. Each additional IPP doubles the potential for structural diversification:

• Hemiterpenes (C5): isoprene itself — the simplest terpenoid, emitted by plants for thermotolerance
• Monoterpenes (C10): from geranyl diphosphate (GPP = IPP + DMAPP). Examples: menthol, limonene, pinene, camphor. Major components of essential oils and floral scent
• Sesquiterpenes (C15): from farnesyl diphosphate (FPP = GPP + IPP). Examples: farnesol, artemisinin (antimalarial), \(\beta\)-caryophyllene. The largest subclass with ~7,000 known structures
• Diterpenes (C20): from geranylgeranyl diphosphate (GGPP = FPP + IPP). Examples: gibberellins (plant hormones), taxol (anticancer), phytol
• Triterpenes (C30): from squalene (2 x FPP). Examples: sterols, saponins, limonoids. ~20,000 known structures
• Tetraterpenes (C40): from phytoene (2 x GGPP). Examples: carotenoids (\(\beta\)-carotene, lycopene, lutein)

Combinatorial Diversity

The staggering structural diversity of terpenoids arises from the combinatorial explosion of post-assembly modifications. For a terpenoid skeleton of \(n\) carbon atoms, the potential diversity can be estimated from:

For sesquiterpenes (C15): ~20 basic skeletal types (from different cyclization patterns of FPP) \(\times\) ~10 oxidation states per skeleton \(\times\) ~5-20 decorations (glycosylation, acylation, methylation) yields a theoretical space of ~1,000-4,000 structures from a single C15 chain — matching the ~7,000 observed sesquiterpenes when considering stereochemistry.

The key enzymes generating this diversity are the terpene synthases (TPSs) — a remarkably evolvable enzyme family where single amino acid changes can redirect cyclization cascades to produce entirely different carbon skeletons. This evolvability makes terpenoid diversity a particularly fertile substrate for natural selection in chemical defense.

1.6 Phenolic Compounds & the Shikimate Pathway

Plant phenolics — encompassing phenylpropanoids, flavonoids, stilbenes, lignans, and lignin — constitute a vast class of secondary metabolites all derived from the shikimate pathway. This pathway connects primary carbohydrate metabolism to aromatic amino acid biosynthesis and is the gateway to the entire phenylpropanoid arsenal.

From Shikimate to Defense: The Biosynthetic Route

The shikimate pathway converts phosphoenolpyruvate (PEP, from glycolysis) and erythrose-4-phosphate (E4P, from the pentose phosphate pathway) into chorismate through seven enzymatic steps. Chorismate is then channeled into three aromatic amino acids:

The gateway enzyme to phenylpropanoid metabolism is phenylalanine ammonia-lyase (PAL), which catalyzes the non-oxidative deamination of phenylalanine to trans-cinnamic acid:

PAL is one of the most important regulatory enzymes in plant secondary metabolism. Its expression is induced by UV light, pathogen attack, wounding, and nutrient stress. PAL activity diverts up to 20% of photosynthetically fixed carbon into the phenylpropanoid pathway in stressed plants, representing a massive metabolic commitment.

UV-B Protection by Flavonoids

Flavonoids are among the most important UV-screening compounds in plants. They accumulate in the upper epidermis of leaves where they absorb damaging UV-B radiation (280-315 nm) before it reaches the photosynthetic mesophyll. The physics of UV absorption follows the Beer-Lambert law:

where \(A\) is absorbance (dimensionless), \(\varepsilon\) is the molar extinction coefficient (L mol⁻¹ cm⁻¹), \(c\) is the concentration of the absorbing species, and \(l\) is the path length through the absorbing medium. The fraction of UV radiation transmitted through the epidermal layer is:

Typical molar extinction coefficients for flavonoids at their absorption maxima:

• Quercetin (flavonol): \(\varepsilon_{370} \approx 27{,}000\) L mol⁻¹ cm⁻¹ (Band I); \(\varepsilon_{256} \approx 21{,}000\) (Band II)
• Luteolin (flavone): \(\varepsilon_{348} \approx 18{,}000\) L mol⁻¹ cm⁻¹
• Sinapic acid (hydroxycinnamic): \(\varepsilon_{325} \approx 15{,}000\) L mol⁻¹ cm⁻¹

With typical epidermal flavonoid concentrations of \(\sim 0.5\) mM and epidermal thickness of \(\sim 20\) \(\mu\)m (0.002 cm), the absorbance is:

This seems low, but epidermal cells concentrate flavonoids in vacuoles where local concentrations can reach 10-100 mM, yielding \(A > 1\) (blocking \(>90\%\)of UV). Furthermore, multiple layers of epidermis and the cuticular wax layer contribute additional absorption.

Tannin-Protein Binding and Precipitation

Tannins are complex polyphenols that function as quantitative defenses by binding and precipitating dietary proteins in herbivore guts. The tannin-protein precipitation depends on the tannin/protein ratio. At low tannin concentrations, each tannin molecule binds to multiple sites on different protein molecules, forming soluble complexes. At higher concentrations, cross-linking between proteins causes precipitation. The precipitation curve follows:

where \(r = [\text{tannin}]/[\text{protein}]\) is the tannin/protein ratio,\(r_{50}\) is the ratio at which 50% of maximum precipitation occurs, and\(n\) is the Hill coefficient reflecting cooperativity. At very high tannin concentrations, the precipitate can re-dissolve (the "redissolution effect") as excess tannin saturates binding sites without cross-linking.

1.7 Induced vs Constitutive Defenses: The Jasmonate Switch

Plants face a fundamental strategic choice in chemical defense: maintain defenses continuously (constitutive) or produce them only when attacked (induced). Each strategy has distinct costs and benefits, and the optimal choice depends on the ecology of the plant-herbivore interaction.

Jasmonic Acid (JA) Signaling Pathway

The jasmonic acid (JA) signaling pathway is the master regulatory circuit for induced defenses against herbivores. The signal transduction cascade proceeds as follows:

1. Wounding: herbivore feeding damages cell membranes, releasing linolenic acid (18:3) from membrane phospholipids via phospholipase A
2. JA biosynthesis: linolenic acid is converted to 12-oxo-phytodienoic acid (OPDA) via the octadecanoid pathway (LOX → AOS → AOC), then reduced and \(\beta\)-oxidized to jasmonic acid
3. JA-Ile conjugation: JA is conjugated to isoleucine by JAR1 to form the bioactive jasmonoyl-isoleucine (JA-Ile)
4. JAZ repressor degradation: JA-Ile promotes binding of JAZ (Jasmonate ZIM-domain) repressor proteins to the F-box protein COI1 (CORONATINE INSENSITIVE 1), targeting JAZ for ubiquitin-mediated proteasomal degradation
5. Defense gene activation: removal of JAZ repressors liberates MYC2 and other transcription factors, activating defense genes including protease inhibitors, polyphenol oxidase, and terpenoid synthases

Optimal Defense Timing: Cost-Benefit Analysis

The total cost of defense over the plant's lifetime depends on the strategy chosen. For a constitutive defense maintained over total time \(t_\text{total}\):

where \(c_1\) is the cost per unit time. For an induced defense activated only during attacks (with frequency \(f\) and duration \(t_\text{attack}\)):

where \(c_2\) is the per-unit-time cost of induced defense (typically higher than \(c_1\) due to rapid mobilization) and \(c_\text{delay}\) is the cost incurred during the lag period before defenses are fully deployed (tissue lost to herbivory during the 12-48 hour induction window).

Induced defense is cheaper when \(C_\text{induced} < C_\text{constitutive}\):

This predicts that induced defense is favored when attack frequency is low, consistent with the observation that long-lived trees in species-rich forests (where any individual herbivore species is rare) rely more heavily on inducible defenses than do fast-growing herbs in herbivore-dominated grasslands.

Systemic Acquired Resistance (SAR)

In addition to the JA pathway for herbivore defense, plants employ systemic acquired resistance (SAR) mediated by salicylic acid (SA) for pathogen defense. Local infection triggers SA biosynthesis, which is converted to methyl salicylate (MeSA) for systemic transport through the phloem. In distal tissues, MeSA is hydrolyzed back to SA, activating NPR1 (NONEXPRESSOR OF PR GENES 1), which induces pathogenesis-related (PR) protein expression.

The JA and SA pathways are often antagonistic: activating JA-mediated herbivore defenses suppresses SA-mediated pathogen defenses and vice versa. This creates a defense dilemma — plants under simultaneous attack by herbivores and pathogens must prioritize one defense over the other, a vulnerability exploited by some pathogens that manipulate JA signaling to suppress SA-dependent immunity.

Python: Terpenoid Diversity, UV Protection & Defense Optimization

This three-panel simulation explores: (1) the distribution of known terpenoid compounds across structural classes from hemiterpenes (C5) to tetraterpenes (C40), showing how IPP/DMAPP units combine to generate extraordinary chemical diversity; (2) UV absorption spectra of major flavonoid classes showing their role as UV sunscreens, with the UV-B and UV-A spectral regions highlighted; and (3) the cost-benefit comparison of constitutive vs induced defense strategies as a function of herbivore attack frequency, revealing the crossover point where the optimal strategy switches.

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.
10. Gershenzon, J. & Dudareva, N. (2007). The function of terpene natural products in the natural world. Nature Chemical Biology, 3(7), 408-414.
11. Vogt, T. (2010). Phenylpropanoid biosynthesis. Molecular Plant, 3(1), 2-20.
12. Agati, G. & Tattini, M. (2010). Multiple functional roles of flavonoids in photoprotection. New Phytologist, 186(4), 786-793.
13. Wasternack, C. & Hause, B. (2013). Jasmonates: biosynthesis, perception, signal transduction and action in plant stress responses. Annals of Botany, 111(6), 1021-1058.
14. Kessler, A. & Baldwin, I.T. (2002). Plant responses to insect herbivory: the emerging molecular analysis. Annual Review of Plant Biology, 53, 299-328.
15. Hagerman, A.E. & Butler, L.G. (1981). The specificity of proanthocyanidin-protein interactions. Journal of Biological Chemistry, 256(9), 4494-4497.