Module 5: Root Biochemistry & Rhizosphere

5.1 Root Architecture & Lateral Root Formation

Root system architecture (RSA) determines the volume of soil explored and thus the efficiency of water and nutrient capture. RSA is plastic — roots proliferate in nutrient-rich patches (foraging response) and avoid toxic zones. Lateral roots (LRs) originate from pericycle cells opposite to xylem poles in a process exquisitely controlled by auxin.

Auxin Maxima and LBD Transcription Factors

The PIN-FORMED (PIN) auxin efflux carriers, particularly PIN1 and PIN2, create an auxin concentration gradient along the root axis: maximum at the root tip (quiescent centre, columella initials), declining toward the elongation zone. Lateral root initiation is triggered by local auxin maxima in pericycle founder cells. The auxin receptor TIR1 (an F-box protein within the SCF^TIR1 E3 ubiquitin ligase complex) binds auxin:

LATERAL ORGAN BOUNDARIES DOMAIN (LBD) transcription factors LBD16, LBD18, and LBD29 are direct targets of AUXIN RESPONSE FACTOR 7 (ARF7) and ARF19. LBD16 activates cell-cycle genes (CYCD3;1, CDKA;1) in pericycle founder cells, driving anticlinal divisions that produce the lateral root primordium (LRP). The LRP passes through overlying cortex and epidermis via cell-wall remodeling by polygalacturonases and expansins, emerging as a new root axis 2–4 days after initiation in Arabidopsis.

Nitrogen-Responsive Root Foraging

Local NO₃⁻ supply stimulates lateral root elongation independent of shoot-derived signals. The dual-affinity nitrate transporter NRT1.1 (CHL1/NPF6.3) acts as a nitrate sensor: at low [NO₃⁻] it is phosphorylated at Thr101 by the kinase CIPK23, switching it to high-affinity mode and simultaneously triggering auxin import into lateral root tips, promoting LR elongation. This represents a single protein acting as both transporter and transceptor (transporter + receptor).

5.2 Ion Channel Biophysics & Membrane Potential

Goldman-Hodgkin-Katz Equation: Full Derivation

The resting membrane potential of root hair and cortical cells arises from the selective permeability of the plasma membrane to different ions. The Goldman-Hodgkin-Katz (GHK) voltage equation provides a quantitative description under the constant-field assumption (electric field uniform across membrane).

Starting from the Nernst-Planck equation for ionic flux density \( J_i \) of ion species \( i \) with charge \( z_i \):

Under the constant-field assumption \( d\phi/dx = -V_m/L \) (where \( L \) is membrane thickness, \( V_m = \phi_{in} - \phi_{out} \)), the ODE for each ion becomes:

This linear first-order ODE has the integrating factor \( \mu = e^{-z_i F V_m x/(RTL)} \). Applying boundary conditions \( C_i(0) = C_i^{out} \), \( C_i(L) = C_i^{in} \) and solving:

where \( P_i = D_i \beta_i / L \) is the permeability (\( \beta_i \) = partition coefficient). At electrochemical equilibrium, the net current is zero. For the biologically relevant ions K⁺, Na⁺ (cations, \( z=+1 \)) and Cl⁻ (anion, \( z=-1 \)):

In plant root cells, typical values are: \( [K^+]_o = 1 \) mM, \( [K^+]_i = 100 \) mM,\( [Na^+]_o = 1 \) mM, \( [Na^+]_i = 10 \) mM, with \( P_K \gg P_{Na} \). This typically yields \( V_m \approx -120 \) to \( -180 \) mV, which is the thermodynamic driving force for K⁺ uptake via inward-rectifying channels.

Note that plant cells maintain a far more negative membrane potential than animal cells (~−120 to −180 mV vs ~−70 mV) because the H⁺-ATPase (proton pump) actively extrudes H⁺, generating a large electrochemical gradient that powers nutrient uptake. The H⁺-ATPase consumes ~40% of root ATP production.

K⁺ Channels: KAT1 and SKOR

Two major K⁺ channel families mediate K⁺ fluxes across the root plasma membrane: KAT1 (K⁺ channel in Arabidopsis thaliana 1), an inward-rectifying channel (Shaker family), and SKOR (Stelar K⁺ Outward Rectifier), which mediates K⁺ release into the xylem. Both exhibit strong voltage dependence described by a two-state Boltzmann model:

where \( V_{1/2} \) is the half-activation voltage and \( k_s \) is the slope factor (mV per e-fold change). For KAT1: \( V_{1/2} \approx -100 \) mV, \( k_s \approx -15 \) mV (activates upon hyperpolarization). For SKOR: \( V_{1/2} \approx -40 \) mV,\( k_s \approx +20 \) mV (activates upon depolarization, driven by K⁺ concentration in the stele).

KAT1 has a characteristic domain structure: six transmembrane helices (S1–S6), a voltage-sensing domain (S1–S4) with positively charged Arg/Lys residues in S4, a pore-forming loop between S5 and S6 containing the selectivity filter (GYGV motif), and a cytoplasmic CNBD (cyclic nucleotide binding domain) whose occupation by cAMP slightly modulates gating. The gating charge \( z_g \approx 1.9 \)elementary charges per subunit, with four subunits per functional channel, giving a macroscopic gating valence of ~7.6.

5.3 Root Cross-Section: Anatomy & Transport Pathways

Water and solutes move from soil to xylem via three pathways: apoplastic (cell walls), symplastic (plasmodesmata), and transmembrane. The Casparian strip forces all solutes through the symplast at the endodermis, enabling selectivity. A second barrier — the suberin lamellae — develops in older roots, creating two checkpoints.

Root cross-section: epidermis (outermost, with root hairs), cortex, endodermis with Casparian strip (dashed yellow, forces symplastic transport), pericycle (site of lateral root initiation), and the vascular stele with alternating xylem (blue) and phloem (orange) poles.

5.4 Nitrate Transport: Dual-Affinity Systems

Nitrogen is typically the most limiting macronutrient for tree growth. Uptake of inorganic nitrogen (NO₃⁻, NH₄⁺) from soil solution involves a suite of transporter families with complementary concentration ranges and regulatory modes.

NRT1.1 (NPF6.3): The Transceptor

NRT1.1 is the archetypal plant transceptor: a proton-coupled nitrate transporter that also functions as a nitrate sensor. It belongs to the NPF (Nitrate Peptide transporter Family, also called NRT1/PTR). NRT1.1 exhibits dual-affinity kinetics with a transition governed by phosphorylation at Thr101 by CIPK23 (CBL-interacting protein kinase 23):

NRT2 Family: High-Affinity NO₃⁻ Transport

NRT2.1 is the primary high-affinity transporter (HATS) for NO₃⁻ uptake at concentrations below ~1 mM (typical for forest soils). It is an H⁺/NO₃⁻ co-transporter requiring the accessory protein NAR2.1 (NRT3.1) for membrane trafficking and function. The stoichiometry is 2 H⁺ per NO₃⁻, driven by the plasma membrane H⁺ electrochemical gradient. NRT2.1 expression is induced by nitrogen starvation and suppressed by high internal amino acid concentrations (feedback by glutamine via the NLP transcription factor pathway).

Overall Michaelis-Menten kinetics for the combined HATS + LATS system:

At soil [NO₃⁻] around 0.01–0.05 mM (common in forest floors), the HATS dominates and operates well below saturation, making root surface area and transporter expression rate-limiting. At agricultural NO₃⁻ levels (1–10 mM), the LATS operates near saturation.

5.5 Phosphate Acquisition: PHT1 Transporters & Mycorrhizae

Phosphorus is often the primary limiting nutrient in forest soils because inorganic phosphate (Pi) binds tightly to iron/aluminium oxides (acid soils) or calcium minerals (calcareous soils), limiting its concentration in soil solution to just 1–10 μM. Root hairs partially compensate through the depletion zone effect, but the primary adaptation in trees is mycorrhizal symbiosis.

PHT1 Transporters

PHT1 (Phosphate Transporter 1) family members are H⁺/Pi co-transporters belonging to the MFS (Major Facilitator Superfamily). They have 12 transmembrane domains and operate with a stoichiometry of 2–4 H⁺ per H₂PO₄⁻ (at soil pH). The electrochemical driving force is provided entirely by the H⁺-ATPase. In trees, multiple PHT1 genes are expressed with organ-specific patterns: PHT1;1 and PHT1;4 are expressed in root hairs (high affinity, Km ~3 μM), while PHT1;2 is expressed in inner cortex. Under Pi starvation, PHT1 expression is upregulated 10–100-fold via PHR1 (PHOSPHATE STARVATION RESPONSE 1) — a MYB-CC transcription factor that also induces root architecture remodeling.

Arbuscular Mycorrhizal (AM) Symbiosis

Over 80% of land plant species form arbuscular mycorrhizal (AM) symbioses with Glomeromycota fungi. The fungal hyphae extend far beyond the root depletion zone, effectively increasing the soil volume explored by 100–1000-fold. Pi acquired by the extraradical mycelium is converted to polyphosphate (polyP) granules (10–20 Pi units) and transported via cytoplasmic streaming to the periarbuscular membrane interface, where it is hydrolyzed by vacuolar phosphatases to free Pi for delivery to the plant.

Ectomycorrhizal Networks in Forests

Ectomycorrhizal (ECM) fungi (Basidiomycota and Ascomycota) form a Hartig net — a hyphal network between cortical cells — and a fungal mantle surrounding the root tip. Unlike AM fungi, ECM hyphae do not penetrate cell walls. ECM fungi include Amanita, Suillus, Paxillus, and Cenococcum.

Common Mycorrhizal Networks (CMN) arise when a single fungal individual simultaneously colonizes multiple trees (even different species). CMNs have been demonstrated to transfer ¹³C-labelled carbon from large "mother trees" to seedlings in the understory. The mechanism involves both mass flow along water potential gradients and active cytoplasmic streaming. The ecological implications are profound: CMNs may buffer seedling carbon deficits in deep shade, and may mediate below-ground signaling of herbivore attack.

5.6 Root Exudates: Chemical Ecology of the Rhizosphere

Trees release 5–21% of their net photosynthetically fixed carbon as root exudates: low-molecular-weight organic acids, sugars, amino acids, and specialized signaling compounds. These exudates profoundly shape the rhizosphere microbiome and determine nutrient availability.

Frankia Nitrogen Fixation in Alder

Actinorhizal plants (including Alnus spp., alder) form nitrogen-fixing nodules with Frankia actinobacteria. Alder trees can fix 40–300 kg N ha⁻¹ yr⁻¹, making them ecosystem engineers of riparian and early-successional forests.

The nitrogenase complex in Frankia consists of dinitrogenase (MoFe protein, α₂β₂ tetramer with P-clusters and FeMo-cofactor) and dinitrogenase reductase (Fe protein, homodimer). The overall reaction:

The extreme O₂ sensitivity of nitrogenase (irreversible oxidative inactivation of Fe-S clusters) is managed in Frankia by specialized vesicles — thick-walled, lipid-bilayer-surrounded spherical structures at hyphal tips — that create an O₂ diffusion barrier. The vesicle wall contains multiple hopanoid lipid layers, analogous to the cyanobacterial heterocyst envelope. Fixed nitrogen is exported as NH₄⁺ or citrulline to the plant cytoplasm, where glutamine synthetase (GS) and glutamate synthase (GOGAT) assimilate it.

Rhizosphere Microbiome & PGPR

The rhizosphere harbors a microbial biomass 10–100× higher than bulk soil. Plant Growth Promoting Rhizobacteria (PGPR) enhance growth through multiple mechanisms: (1) IAA production (Pseudomonas, Azospirillum) that stimulates root hair initiation; (2) phosphate solubilization via gluconic acid production; (3) siderophore production (pyoverdines, enterobactin) that chelate Fe³⁺, increasing iron availability; (4) ACC deaminase activity (cleaves ethylene precursor 1-aminocyclopropane-1-carboxylic acid), reducing stress ethylene; (5) induction of Induced Systemic Resistance (ISR) — a primed defense state mediated by jasmonate/ethylene signaling that confers broad-spectrum resistance to foliar pathogens without direct PR gene expression.

Python: Nutrient Uptake Kinetics & Depletion Zone Modeling

The following simulation models: (1) dual-affinity NO₃⁻ uptake kinetics (NRT1.1 vs NRT2.1), (2) steady-state radial concentration profile in the rhizosphere depletion zone using the cylindrical diffusion equation, and (3) comparison of Pi uptake via root hairs (PHT1) vs arbuscular mycorrhizal fungi, illustrating the kinetic advantage of the mycorrhizal pathway at the vanishingly low Pi concentrations typical of forest soils.