Module 7: Stress Biophysics โ€” Drought, Cold & Pathogens

Trees are sessile organisms that must endure extreme and unpredictable abiotic stresses โ€” months of drought, temperatures from โˆ’40ยฐC to +45ยฐC โ€” and constant biotic threats from pathogens and herbivores. This module examines the molecular machinery that perceives stress, orchestrates adaptive responses, and maintains cellular integrity under conditions that would rapidly kill most organisms.

7.1 ABA Biosynthesis & Stomatal Signaling Cascade

Abscisic acid (ABA) is the central hormone of abiotic stress response in plants, governing stomatal closure, seed dormancy, and the induction of stress-tolerance genes. Its biosynthesis begins in plastids via cleavage of C40 carotenoids (violaxanthin/neoxanthin) to form xanthoxin, which is then converted to ABA in the cytosol.

NCED: The Committed Step

9-cis-epoxycarotenoid dioxygenase (NCED) is the committed and rate-limiting step of ABA biosynthesis. It is a non-heme iron dioxygenase localized to the chloroplast stroma that cleaves 9โ€ฒ-cis-neoxanthin or 9-cis-violaxanthin to yield xanthoxin (C15) and a C25 byproduct. In trees under drought stress, NCEDtranscript abundance increases dramatically within 30โ€“60 minutes and can rise 10โ€“100-fold. NCED activity is also activated post-translationally by reactive oxygen species (ROS) via oxidation of critical cysteine residues, creating a rapid ABA burst in response to cellular redox signals.

ABA Signaling Pathway: PYR/PYL โ†’ PP2C โ†’ SnRK2 โ†’ SLAC1ABAPYR/PYLReceptor (soluble)bindsPP2C(ABI1, ABI2, HAB1)ABA-PYR/PYLinhibits PP2CSnRK2.2/3/6(OST1/SRK2E)dephosphorylates(inactivates)โœ— (ABA blocks)SLAC1Anion channel (S-type)phosphorylates(activates)Stomatal closureCl-/malate effluxK+ efflux (GORK)ABF/AREB TFsABA-responsive genesactivatesWithout ABA:PP2C active โ†’ dephosphorylates/inactivates SnRK2With ABA:PYR/PYL-ABA complex inhibits PP2C โ†’ SnRK2 activeโ†’ SLAC1 phosphorylated โ†’ Cl-/mal2- effluxโ†’ guard cell turgor loss โ†’ stomatal closure

Core ABA signaling module. Without ABA: PP2C phosphatases inactivate SnRK2 kinases. With ABA: PYR/PYL receptors bind ABA and simultaneously inhibit PP2C (co-receptor mechanism), releasing SnRK2 autophosphorylation โ†’ activation โ†’ SLAC1 anion channel phosphorylation โ†’ stomatal closure.

Osmotic Adjustment: van't Hoff Equation

Under drought, cells accumulate compatible solutes to lower their osmotic potential, maintaining cell turgor as external water potential declines. The van't Hoff equation relates osmotic potential to solute concentration:

\[ \pi = -RT \sum_i c_i \]

where \( \pi \) is osmotic potential (Pa), \( R = 8.314 \) J molโปยน Kโปยน,\( T \) is temperature (K), and \( c_i \) is molar concentration of solute \( i \) (mol mโปยณ). For real solutions, an osmotic coefficient\( \phi \) is introduced: \( \pi = -\phi RT \sum c_i \).

Compatible solutes are small, highly soluble molecules that do not interfere with normal metabolism even at high concentrations:

Proline

Universal drought/salt stress osmolyte; also ROS scavenger; synthesized from glutamate via P5CS enzyme; can reach 100 mM in stressed cells

Glycine betaine

Major osmolyte in mangroves, halophytes, some grasses; synthesized from choline via CMO/BADH; stabilizes PSII and membrane proteins

Sorbitol / Mannitol

Sugar alcohols; common in woody Rosaceae (sorbitol in apple/cherry phloem sap); mannitol in celery, olive; also function as antioxidants

7.2 Freeze-Thaw Embolism: Ice Nucleation & Bubble Dissolution

Freeze-thaw embolism is a distinct vulnerability mechanism from drought embolism. When xylem sap freezes, dissolved gases are expelled from the ice lattice (ice has very low gas solubility). Upon thawing, these gas bubbles must either dissolve or remain as emboli that block water flow.

Ice Nucleation in Xylem

Water in xylem conduits can remain supercooled to โˆ’5ยฐC to โˆ’10ยฐC in narrow vessels before heterogeneous ice nucleation (catalyzed by ice-nucleating proteins on cell walls, bacteria, or mineral surfaces). In wide-diameter vessels typical of ring-porous species (Quercus, Fraxinus), nucleation occurs at higher temperatures (โˆ’1 to โˆ’3ยฐC) due to a lower surface-area-to-volume ratio affecting nucleation kinetics. Once nucleation begins in one vessel, ice propagates through pit membranes to adjacent conduits via freeze-fracture of the pit membrane (a catastrophic process in spring wood).

Critical Bubble Dissolution: Young-Laplace Derivation

On thawing, a spherical gas bubble of radius \( r \) (diameter \( d \)) experiences a Laplace pressure difference between the gas inside and the liquid outside. For a spherical interface with surface tension \( \gamma \):

Consider a virtual displacement: expand the bubble by \( dr \). The work done against surface tension is:

\[ dW_{surface} = \gamma \cdot d(4\pi r^2) = 8\pi\gamma r\,dr \]

The work done by the excess pressure \( \Delta P = P_{in} - P_{out} \) expanding the volume \( dV = 4\pi r^2 dr \) is:

\[ dW_{pressure} = \Delta P \cdot 4\pi r^2 dr \]

At mechanical equilibrium, \( dW_{surface} = dW_{pressure} \):

\[ \Delta P = P_{gas} - P_{liquid} = \frac{4\gamma}{d} \]

Young-Laplace equation for a spherical bubble (factor of 2 for each interface: \( 2 \times 2\gamma/r = 4\gamma/r \))

For a bubble to dissolve, the xylem must be under sufficient tension to overcome this Laplace pressure. The critical dissolution condition:

\[ |\Psi_{xylem}| > \Delta P_{crit} = \frac{4\gamma}{d} \]

For \( \gamma = 0.0728 \) N mโปยน and a typical bubble of \( d = 0.3 \) ฮผm:\( \Delta P_{crit} = 4 \times 0.0728 / (0.3 \times 10^{-6}) \approx 0.97 \) MPa. Typical xylem tensions in temperate trees during spring refilling are just 0.1โ€“0.3 MPa โ€” insufficient to dissolve large bubbles. Bubble dissolution relies instead on active osmotic mechanisms: starch-to-sugar conversion in ray parenchyma creates a high osmotic pressure that drives water into vessels, generating a local positive pressure (root pressure or stem pressure) that compresses and dissolves bubbles.

Ring-Porous vs Diffuse-Porous Vulnerability

Ring-Porous (Quercus, Fraxinus, Ulmus)

Large earlywood vessels (50โ€“300 ฮผm diameter) are extremely vulnerable to freeze-thaw embolism: nearly 100% are embolized after a single freeze event. These species produce new earlywood each spring to replace embolized vessels โ€” ring-porosity is an adaptation to predictable annual embolism, not a vulnerability. Functional xylem is only 1โ€“2 rings wide.

Diffuse-Porous (Betula, Populus, Acer)

Smaller vessels (20โ€“80 ฮผm) throughout growth ring. More resistant to freeze-thaw embolism due to smaller bubble size; can dissolve more easily under spring tensions. Multiple growth rings remain functional. Birch and maple generate high root/stem pressures in spring (maple syrup harvest) to dissolve winter embolism. The pressure arises from starch hydrolysis + COโ‚‚ dissolution in stem parenchyma.

Freeze-Thaw Embolism Cycle in Xylem Vessel1. Liquid (winter)Water under tensionpsi = -1 to -2 MPaFreeze2. Ice (frozen)ICE CRYSTALGas expelled from latticeThaw3. Post-thawGASBUBBLEembolism!Root/stempressure4. RefilledBubble dissolved+P from ray parenchymaCritical dissolution: |psi_xylem| must exceed 4*gamma/dRing-porous: large bubbles, usually not dissolved by tension aloneDiffuse-porous: small bubbles, dissolved by moderate tension or root/stem pressureConifer tracheids: small diameter = high resistance to freeze-thaw embolism

7.3 Cold Acclimation: Membrane Adaptation & Cryoprotection

Cold acclimation is the process by which plants exposed to low (but above-freezing) temperatures progressively develop tolerance to subsequent freezing. In trees, acclimation begins in early autumn triggered by short photoperiod and low temperature, and involves coordinated changes at molecular, membrane, and cellular levels.

Membrane Lipid Desaturation

Cold temperatures dramatically increase membrane lipid phase transition temperatures, risking rigidification of membranes at low temperatures. To compensate, FAD (FATTY ACID DESATURASE) enzymes introduce additional double bonds into fatty acyl chains:

FAD2 (ER)

18:1โ†’18:2 (ฮ”12 desaturase), adds second double bond at C12

FAD3 (ER)

18:2โ†’18:3 (ฯ‰-3 desaturase), adds third double bond at C15

FAD7/FAD8 (plastid)

18:2โ†’18:3 in plastid phospholipids; FAD8 cold-induced specifically

In cold-acclimated trees, the ratio of unsaturated:saturated fatty acids in plasma membrane phosphatidylglycerol can shift from ~1:3 to ~3:1, maintaining membrane fluidity at โˆ’10ยฐC that would otherwise be solid at those temperatures.

LEA Proteins & Cryoprotection

LATE EMBRYOGENESIS ABUNDANT (LEA) proteins are intrinsically disordered proteins (IDPs) induced by ABA and low temperature. Their protective mechanisms include: (1) vitrification โ€” they form a highly viscous glass-like matrix in dehydrated cells, preventing membrane fusion and protein aggregation; (2) molecular shielding โ€” they coat membrane surfaces, replacing water molecules and preventing lipid phase separation during freeze-induced dehydration; (3) ion buffering โ€” their multiple charged residues bind ions that would otherwise reach damaging concentrations as intracellular water freezes.

7.4 Pathogen Defense: PTI, ETI & the Zigzag Model

Plant immunity evolved as a two-tiered system (Jones & Dangl, 2006): Pattern-Triggered Immunity (PTI) provides a baseline defense against all microbes, while Effector-Triggered Immunity (ETI) provides a more powerful response against adapted pathogens recognized via their secreted virulence effectors.

The Zigzag Model (Jones & Dangl 2006)

โ‘ 

PAMP Recognition (PTI): PRRs (Pattern Recognition Receptors) detect PAMPs/MAMPs (flagellin, chitin, EF-Tu, LPS) โ†’ MAPK cascade โ†’ PR gene expression, callose deposition, ROS burst, stomatal closure

โ‘ก

Effector Delivery (ETS): Adapted pathogens deliver type III secretion effectors (bacteria) or RXLR effectors (oomycetes) into plant cells to suppress PTI

โ‘ข

Effector Recognition (ETI): NBS-LRR resistance (R) proteins detect effectors directly or via the "guard hypothesis" (guardee modification) โ†’ HR (hypersensitive response), strong defense amplification

โ‘ฃ

Effector Evolution (ETSยฒ): New pathogen effectors suppress ETI; new R proteins evolve in an evolutionary arms race

SAR: Systemic Acquired Resistance

SAR is a whole-plant defense state, established after local pathogen attack, that primes distal tissues for faster and stronger defense activation. The signal molecule salicylic acid (SA) is central:

Pathogen attack โ†’ SA biosynthesis (ICS1/PBS3 pathway: chorismate โ†’ SA) โ†’ SA โ†’ methyl SA (MeSA, volatile) โ†’ systemic transport via phloem โ†’ distal leaves: MeSA โ†’ SA (SABP2 methylesterase) โ†’ NPR1 activation โ†’ TGA TF complex โ†’ PR1, PR2, PR5 gene expression โ†’ broad-spectrum resistance

NPR1 (NON-EXPRESSOR OF PR GENES 1) is the SA receptor and master regulator of SAR. Without SA, NPR1 exists as an oligomer in the cytoplasm, held together by intermolecular disulfide bonds. SA binding to NPR1 (and the related NPR3/NPR4) triggers monomerization (via reduction of Cys82 and Cys216 by TRX) and nuclear import, where NPR1 binds TGA transcription factors to activate PR genes.

Jasmonate Signaling: LOX-AOS-AOC Pathway & COI1-JAZ

Jasmonate (JA and its bioactive conjugate JA-Ile) is the primary defense hormone against necrotrophic pathogens and herbivores. Biosynthesis begins in the chloroplast from ฮฑ-linolenic acid (18:3):

\[ \text{18:3-PE} \xrightarrow{LOX2} \text{13-HPOT} \xrightarrow{AOS} \text{12,13-EOT} \xrightarrow{AOC} \text{OPDA} \xrightarrow{\text{OPR3 (perox.)}} \text{JA} \xrightarrow{JAR1} \text{JA-Ile} \]

The bioactive JA-Ile is perceived by the SCFCOI1 E3 ubiquitin ligase: JA-Ile acts as a "molecular glue" bridging COI1 (F-box subunit) and JAZ (JASMONATE ZIM-domain) repressor proteins. This triggers JAZ ubiquitination and 26S proteasome degradation, releasing MYC2/MYC3/MYC4 transcription factors to activate jasmonate-responsive genes including:

Terpene synthases (TPS) โ€” direct defense volatiles

VSP2 (Vegetative Storage Protein 2) โ€” anti-nutritional to insects

WRKY TFs โ€” amplification cascade

JAZ genes themselves โ€” negative feedback loop

LOX2 โ€” amplification of JA burst

PI1/PI2 (Proteinase Inhibitors) โ€” anti-herbivore

Python: Vulnerability Curves, Drought Dynamics & Bubble Dissolution

Three interconnected simulations: (1) species-specific Weibull vulnerability curves showing P50 and P88 for ring-porous vs diffuse-porous vs tracheid-based xylem, (2) dynamic simulation of drought-induced conductivity loss comparing isohydric (Picea-like) vs anisohydric (Quercus-like) stomatal strategies, and (3) the Young-Laplace critical bubble dissolution pressure as a function of bubble diameter, explaining why ring-porous species cannot dissolve freeze-thaw bubbles by tension alone.

Python
script.py117 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server