Abstract
Air pollution caused by tropospheric ozone contributes to the decline of forest ecosystems; for instance, sacred fir, Abies religiosa (Kunth) Schltdl. & Cham. forests in the peri-urban region of Mexico City. Individual trees within these forests exhibit variation in their response to ozone exposure, including the severity of visible symptoms in needles. Using RNA-Seq metatranscriptomic data and ITS2 metabarcoding, we investigated whether symptom variation correlates with the taxonomic and functional composition of fungal mycobiomes from needles collected in this highly polluted area in the surroundings of Mexico City. Our findings indicate that ozone-related symptoms do not significantly correlate with changes in the taxonomic composition of fungal mycobiomes. However, genes coding for 30 putative proteins were differentially expressed in the mycobiome of asymptomatic needles, including eight genes previously associated with resistance to oxidative stress. These results suggest that fungal communities likely play a role in mitigating the oxidative burst caused by tropospheric ozone in sacred fir. Our study illustrates the feasibility of using RNA-Seq data, accessible from global sequence repositories, for the characterization of fungal communities associated with plant tissues, including their gene expression.