Uncoupling the Effects of Highland Maize Chromosomal Inversion Inv4m from Leaf Phosphorus Deficiency Responses

Technical Summary: Uncoupling the Effects of Highland Maize Chromosomal Inversion Inv4m from Leaf Phosphorus Deficiency Responses

Problem Statement
Local adaptation in maize involves the selection of adaptive alleles, often preserved by chromosomal inversions that act as "supergenes" by suppressing recombination. The Inv4m inversion (13 Mb) is prevalent in highland maize and its wild relative Zea mays ssp. mexicana, and is nearly fixed in populations above 2,500 meters. While Inv4m is strongly associated with local adaptation to high-altitude environments (influencing flowering time and cold response), the specific physiological mechanisms driving its adaptive value remain unclear. The highland environments of the Trans-Mexican volcanic belt are characterized by acidic volcanic soils (Andosols) with low phosphorus (P) availability. Previous studies suggested Inv4m might carry alleles enhancing P acquisition or utilization, given the presence of genes like ZmPho1;2a within the inversion. However, it was unknown whether Inv4m confers a specific advantage under P starvation or if its adaptive value is driven by constitutive developmental effects independent of nutrient stress. Furthermore, the interaction between P starvation and the natural gradient of sequential leaf senescence in maize remains mechanistically poorly understood.

Methodology
The study employed a multi-omics approach to dissect the effects of Inv4m and P deficiency.

• Genetic Material: Near-Introgression Lines (NILs) were developed by backcrossing a highland maize variety (MICH21, carrying the Inv4m highland haplotype) into the temperate B73 genetic background for eight generations. This created pairs of lines differing primarily at the Inv4m locus (Inv4m vs. Control/CTRL).
• Experimental Design: NILs were grown in a field trial at Rocksprings, PA, under two phosphorus treatments: sufficient (+P) and deficient (-P).
• Sampling Strategy: To capture the interplay between nutrient stress and developmental stage, specific leaf positions were sampled at 63 Days After Planting (approx. V13 stage). Four leaves per plant were collected, ranging from the youngest fully expanded "collar" leaf (Node 0) to older basal leaves (Nodes 2, 4, and 6), representing a gradient of sequential senescence.
• Omics Analyses:
  • Phenotyping: Measured anthesis, silking, biomass accumulation (stover dry weight), kernel weight, and cob diameter.
  • Ionomics: Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES) was used to quantify mineral concentrations (P, Zn, Ca, S, Mg, Mn, K, Fe) in stover and grain.
  • Transcriptomics: RNA-Seq was performed on the four leaf positions to analyze gene expression.
  • Lipidomics: HPLC-MS was used to profile membrane lipids and storage lipids.
• Statistical Modeling: Linear models were fitted to test for main effects of Phosphorus, Genotype (Inv4m), Leaf Stage, and their interactions (e.g., Leaf × P, Genotype × P).

Key Results

1. Phenotypic Independence of Inv4m from P Stress: Phosphorus deficiency caused severe phenotypic penalties, including delayed flowering, reduced biomass, and lower kernel weight, in both Inv4m and CTRL lines. Crucially, there was no significant Genotype × Phosphorus interaction for primary agronomic traits. The Inv4m inversion did not alter the magnitude or direction of the P starvation response. The only significant interaction was observed for cob diameter, where Inv4m lines produced thinner cobs under P deficiency, but this was a secondary trait.
2. Conserved Ionomic Responses: P deficiency induced strong shifts in mineral accumulation (e.g., decreased P in stover and grain, increased Zn in stover) that were largely shared between Inv4m and CTRL lines. While Inv4m showed additive effects on certain mineral contents (e.g., lower stover Ca, Zn, and K), these were not specific to the P stress response.
3. Transcriptomic and Lipidomic Gradients Driven by Leaf Age: The molecular response to P starvation was heavily modulated by leaf developmental stage.
   • Senescence Gradient: Older leaves exhibited a transcriptional profile characterized by downregulation of photosynthesis genes (e.g., lhcb10, rca3) and upregulation of senescence-associated genes (SAGs).
   • Amplified Stress Response: The magnitude of the P starvation response increased with leaf age. Older leaves showed a "coordinated decommissioning" of light harvesting and carbon fixation machinery, alongside an acceleration of nutrient remobilization pathways.
   • Lipid Remodeling: P starvation triggered a canonical shift from phospholipids to galactolipids and the accumulation of triacylglycerols (TAGs). This remodeling was age-dependent; for instance, the depletion of phosphatidylcholine PC34:2 and the accumulation of specific TAGs (e.g., TG56:6) were exacerbated in older leaves.
   • LPE Divergence: Unlike in 16:3 plants (e.g., Arabidopsis) where lysophosphatidylethanolamines (LPE) act as senescence suppressors, LPE accumulation in this 18:3 plant (maize) correlated with leaf age and senescence, suggesting a species-specific metabolic role.
4. Absence of Inv4m Modulation: No significant interaction between Inv4m genotype and P deficiency was detected in the transcriptome, lipidome, or ionome. Genes within the inversion, including the candidate ZmPho1;2a, did not show differential regulation specific to the Inv4m genotype under P stress.

Significance and Claims
The paper concludes that the adaptive value of the Inv4m inversion in highland maize is likely not driven by an enhanced capacity to cope with phosphorus deficiency. Instead, the study suggests that Inv4m operates through constitutive developmental effects (such as phenology and flowering time) rather than by modifying specific nutrient stress response pathways.

The study provides a detailed characterization of how phosphorus starvation interacts with natural leaf senescence in maize. It identifies two coordinately regulated molecular programs triggered by P limitation:

1. Light harvesting shutdown: A progressive downregulation of photosynthetic capacity in older leaves.
2. Accelerated senescence: An upregulation of nutrient remobilization and catabolic processes in older tissues.

The authors argue that older leaves integrate nutrient limitation into the natural progression toward senescence, prioritizing the remobilization of resources to younger, photosynthetically active tissues and developing reproductive structures. This multi-omics framework clarifies that the "P starvation response" in maize is not a uniform whole-plant reaction but a spatially and temporally dynamic process heavily dependent on leaf developmental stage. The findings suggest that crop improvement strategies might benefit from optimizing the timing and tissue specificity of senescence responses rather than assuming uniform stress tolerance mechanisms.