Tropospheric ozone is a major air pollutant that significantly damages crop production around the world. Crop metabolic responses to rising chronic ozone stress have not been well-studied in the field, especially in C4 crops. In this study, we investigated the metabolomic profile of leaves from two diverse maize (Zea mays) inbred lines and the hybrid cross during exposure to season-long elevated ozone (~100 nL L−1) in the field using free air concentration enrichment (FACE) to identify key biochemical responses of maize to elevated ozone. Senescence, measured by loss of chlorophyll content, was accelerated in the hybrid line, B73 x Mo17, but not in either inbred line (B73 or Mo17). Untargeted metabolomic profiling further revealed that inbred and hybrid lines of maize differed in metabolic responses to ozone. A significant difference in the metabolite profile of hybrid leaves exposed to elevated ozone occurred as leaves aged, but no age-dependent difference in leaf metabolite profiles between ozone conditions was measured in the inbred lines. Phytosterols and α-tocopherol levels increased in B73 x Mo17 leaves as they aged, and to a significantly greater degree in elevated ozone stress. These metabolites are involved in membrane stabilization and chloroplast reactive oxygen species (ROS) quenching. The hybrid line also showed significant yield loss at elevated ozone, which the inbred lines did not. This suggests that the hybrid maize line was more sensitive to ozone exposure than the inbred lines, and up-regulated metabolic pathways to stabilize membranes and quench ROS in response to chronic ozone stress.

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the NAM pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

The maize pangenome has demonstrate large amounts of presence/absence variation and it has been hypothesized that presence/absence variation contributes to stress response. To uncover whether the observed genetic variation in physiological response to elevated ozone (a secondary air pollutant that causes significant crop yield losses) concentration is due to variation in genic content, and/or variation in gene expression, we examine the impact of sustained elevated ozone concentration on the leaf tissue from 5 diverse maize inbred genotypes (B73, Mo17, Hp301, C123, NC338). Analysis of long reads from the transcriptomes of the 10 conditions found expressed genes in the leaf are part of the shared genome, with 94.5% of expressed genes from syntenic loci. Quantitative analysis of short reads from 120 plants (twelve from each condition) found limited transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes), while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. The genes underpinning the divergence of B73 from the other 4 genotypes implicates ethylene signaling consistent with some findings in Arabidopsis. For the 82 of the 83 genes differentially expressed among all 5 genotypes and the 788 of 789 genes differentially expressed in 4 genotypes (excluding B73) in sensitivity to ozone is associated with oxidative stress tolerance being associated with a weaker response to a reactive oxygen species (ROS) signal and suggests that genetic variation in downstream processes is key to ozone tolerance.