# Integrated Gene Regulatory Network Analysis Reveals Coordinated Transcriptional Reprogramming in the Arabidopsis thaliana–Trichoderma atroviride Interaction

**Authors:** Evelyn Sánchez, Lorena Melet, José D. Fernández, Tomás C. Moyano, Jonathan Canan, Diego Pérez-Stuardo, Sebastián Reyes-Cerpa, Nathan R. Johnson, Consuelo Olivares-Yáñez, Elena A. Vidal

PMC · DOI: 10.3390/plants15040533 · 2026-02-08

## TL;DR

This study explores how plants and fungi coordinate their gene activity during a mutualistic interaction, revealing key regulatory networks and small RNA dynamics.

## Contribution

The paper introduces integrated gene regulatory networks and sRNA analysis to uncover coordinated transcriptional reprogramming in a plant-fungus mutualism.

## Key findings

- Arabidopsis activates hypoxia, stress, and root development pathways during interaction with Trichoderma atroviride.
- Fungal transcriptional reprogramming focuses on membrane transport and primary metabolism, coordinated by key transcription factors.
- T. atroviride shows strong induction of 20–22 nt sRNAs that target host genes related to root function and immunity.

## Abstract

Mutualistic interactions between plants and beneficial fungi rely on extensive transcriptional reprogramming in both partners, yet the underlying regulatory mechanisms coordinating these responses remain incompletely understood. Here, we combined a transcriptomics analysis with a gene regulatory network (GRN) inference to dissect the interaction between Arabidopsis thaliana and the growth-promoting fungus Trichoderma atroviride. At an early but established stage of colonization (72 h post-inoculation), we identified widespread transcriptional changes in both of the organisms, including host activation of hypoxia, stress and root development-related pathways in Arabidopsis, and fungal reprogramming of membrane transport and primary metabolism. Using DNA-binding motifs and GENIE3-based regulatory inference, we reconstructed interaction-specific GRNs for each species. The subnetworks focused on the main differentially expressed biological processes and uncovered ERF-, WRKY-, NAC- and DOF-centered modules linking hypoxia responses with developmental remodeling in the plant, whereas the putative orthologs of TFs involved in developmental and stress-related TFs in fungi, such as CrzA, RME1, NsdC, PacC and RPN4, formed a regulatory core coordinating fungal transport and metabolic adjustment. In parallel, we uncovered contrasting sRNA dynamics between the partners. While the Arabidopsis sRNA changes were limited, T. atroviride exhibited a strong induction of 20–22 nt sRNAs, including a small set of high-confidence sRNA–mRNA interactions targeting host genes involved in root function and immunity. Together, our results extend previous pathway-based descriptions of the Arabidopsis–Trichoderma mutualism and provide a systems-level, testable framework for how coordinated regulatory programs in both of the partners support the interaction.

## Linked entities

- **Genes:** ERF (ETS2 repressor factor) [NCBI Gene 2077], WRKY (probable WRKY transcription factor 33) [NCBI Gene 103865671], XK (X-linked Kx blood group antigen, Kell and VPS13A binding protein) [NCBI Gene 7504], DOF (Days on feed) [NCBI Gene 100531217], crzA (protein crzA) [NCBI Gene 2872018], rme-1 (Dynamin-type G domain-containing protein) [NCBI Gene 179010], nsdC (nsdC) [NCBI Gene 64847334], pacC (protein pacC) [NCBI Gene 2873843], PSMD9 (proteasome 26S subunit, non-ATPase 9) [NCBI Gene 5715]
- **Species:** Arabidopsis thaliana (taxon 3702), Trichoderma atroviride (taxon 63577)

## Full-text entities

- **Genes:** GBF1 (G-box binding factor 1) [NCBI Gene 829826] {aka AP22.21, AP22_21, AtGBF1, G-box binding factor 1}, BZIP25 (basic leucine zipper 25) [NCBI Gene 824627] {aka ARABIDOPSIS THALIANA BASIC LEUCINE ZIPPER 25, ATBZIP25, BASIC LEUCINE ZIPPER O2 HOMOLOG 4, BZO2H4, basic leucine zipper 25}, GATA12 (GATA transcription factor 12) [NCBI Gene 832652] {aka F18A17.80, F18A17_80, GATA transcription factor 12}, MYB51 (myb domain protein 51) [NCBI Gene 838438] {aka AtMYB51, BW51A, BW51B, F25I16.9, F25I16_9, HIG1}, FLOT1 (SPFH/Band 7/PHB domain-containing membrane-associated protein family) [NCBI Gene 832596] {aka F21J6.108, F21J6_108, flotillin 1}, AT4G19880 (Glutathione S-transferase family protein) [NCBI Gene 827734] {aka T16H5.240, T16H5_240}, ERF6 (ethylene responsive element binding factor 6) [NCBI Gene 827463] {aka ATERF6, DL4780C, ERF-6-6, FCAALL.120, ethylene responsive element binding factor 6}, ERF7 (ethylene response factor 7) [NCBI Gene 821575] {aka ATERF-7, ATERF7, ethylene response factor 7}, OPR1 (12-oxophytodienoate reductase 1) [NCBI Gene 844001] {aka 12-OXOPHYTODIENOATE REDUCTASE 1, 12-oxophytodienoate reductase 1, ARABIDOPSIS 12-OXOPHYTODIENOATE REDUCTASE 1, ATOPR1, F28O16.5, F28O16_5}, ERF71 (Integrase-type DNA-binding superfamily protein) [NCBI Gene 819365] {aka Arabidopsis thaliana ethylene response factor 71, AtERF71, HRE2, HYPOXIA RESPONSIVE ERF (ETHYLENE RESPONSE FACTOR) 2, T30B22.18, ethylene response factor 71}, WRKY8 (WRKY DNA-binding protein 8) [NCBI Gene 834678] {aka ARABIDOPSIS THALIANA WRKY DNA-BINDING PROTEIN 8, ATWRKY8, MPL12.15, MPL12_15, WRKY DNA-binding protein 8}, GATA16 (GATA transcription factor 16) [NCBI Gene 834990] {aka GATA transcription factor 16, K21P3.18, K21P3_18}, WRKY75 (WRKY DNA-binding protein 75) [NCBI Gene 831147] {aka ARABIDOPSIS THALIANA WRKY DNA-BINDING PROTEIN 75, ATWRKY75, T19L5.40, T19L5_40, WRKY DNA-binding protein 75}, TAA1 (uncharacterized protein) [NCBI Gene 843393] {aka CKRC1, F24J13.13, F24J13_13, SAV3, SHADE AVOIDANCE 3, WEAK ETHYLENE INSENSITIVE 8}, TCH3 (Calcium-binding EF hand family protein) [NCBI Gene 818709] {aka ARABIDOPSIS THALIANA CALMODULIN LIKE 4, ATCAL4, CALMODULIN-RELATED PROTEIN, CML12, T3K9.13, T3K9_13}, MIR168a (ncRNA) [NCBI Gene 28720148] {aka MICRORNA 168, MIR168, microRNA168A, p_MI0000210}, AT1G57630 (Toll-Interleukin-Resistance (TIR) domain family protein) [NCBI Gene 842139] {aka T8L23.10, T8L23_10}, NPR1 (regulatory protein (NPR1)) [NCBI Gene 842733] {aka ARABIDOPSIS NONEXPRESSER OF PR GENES 1, ATNPR1, F15H21.6, F15H21_6, NIM1, NON-INDUCIBLE IMMUNITY 1}, AGL44 (AGAMOUS-like 44) [NCBI Gene 815907] {aka AGAMOUS-like 44, ANR1, ARABIDOPSIS NITRATE REGULATED 1, F15N24.5, F15N24_5}, HRS1 (myb-like transcription factor family protein) [NCBI Gene 837890] {aka HYPERSENSITIVITY TO LOW PI-ELICITED PRIMARY ROOT  SHORTENING 1, HYPERSENSITIVITY TO LOW PI-ELICITED PRIMARY ROOT SHORTENING 1, T6J4.6, T6J4_6}, NTM1 (NAC with transmembrane motif1) [NCBI Gene 828149] {aka ANAC068, Arabidopsis NAC domain containing protein 68, F11O4.4, F11O4_4, NAC with transmembrane motif1}, MYB108 (myb domain protein 108) [NCBI Gene 819827] {aka AtMYB108, BOS1, BOTRYTIS-SUSCEPTIBLE1, F5E6.18, F5E6_18, myb domain protein 108}, AT1G02360 (Chitinase family protein) [NCBI Gene 837724] {aka T6A9.5, T6A9_5}, ANAC087 (NAC domain containing protein 87) [NCBI Gene 831945] {aka Arabidopsis NAC domain containing protein 87, MRG7.23, MRG7_23}, MIR403 (ncRNA) [NCBI Gene 5007968] {aka microRNA403, p_MI0001072}, S6K2 (serine/threonine protein kinase 2) [NCBI Gene 820019] {aka ARABIDOPSIS THALIANA PROTEIN KINASE 2, ARABIDOPSIS THALIANA SERINE/THREONINE PROTEIN KINASE 2, ATPK19, ATPK2, ATS6K2, Arabidopsis thaliana protein kinase 19}, YSL8 (YELLOW STRIPE like 8) [NCBI Gene 841257] {aka F11A17.8, F11A17_8, YELLOW STRIPE like 8}, CML38 (calmodulin-like 38) [NCBI Gene 843998] {aka F28O16.2, F28O16_2, calmodulin-like 38}, TIR (toll/interleukin-1 receptor-like protein) [NCBI Gene 843624] {aka AtTN10, F3N23.13, F3N23_13, TIR-nucleotide binding site family 10, TN10, toll/interleukin-1 receptor-like}, GSTF8 (glutathione S-transferase phi 8) [NCBI Gene 819386] {aka ATGSTF5, ATGSTF8, Arabidopsis thaliana glutathione S-transferase phi 8, F17A22.12, GLUTATHIONE S-TRANSFERASE (CLASS PHI) 5, GLUTATHIONE S-TRANSFERASE 16}, GATA3 (GATA transcription factor 3) [NCBI Gene 829620] {aka GATA transcription factor 3, T4L20.260, T4L20_260}, RVE8 (Homeodomain-like superfamily protein) [NCBI Gene 820117] {aka AT3G09610, LCL5, LHY-CCA1-LIKE5, REVEILLE 8}, AT4G30380 (Barwin-related endoglucanase) [NCBI Gene 829161] {aka EXLB2, EXPANSIN-LIKE B2 PRECURSOR, F17I23.280, F17I23_280}, TDR1 (Integrase-type DNA-binding superfamily protein) [NCBI Gene 821901] {aka AtERF98, AtTDR1, ERF98, Ethylene Response Factor 98, Transcriptional Regulator of Defense Response 1}, ERF8 (ethylene response factor 8) [NCBI Gene 841751] {aka ATERF-8, ATERF8, ERF TRANSCRIPTION FACTOR8, ETHYLENE RESPONSE ELEMENT BINDING FACTOR 4, F8L10.19, ethylene response factor 8}, RAP2.9 (uncharacterized protein) [NCBI Gene 826148] {aka DEAR5, DREB AND EAR MOTIF PROTEIN 5, related to AP2 9}, DOF1 (DOF zinc finger protein 1) [NCBI Gene 841595] {aka ADOF1, DOF zinc finger protein 1, F19C24.9, F19C24_9}, GNC (GATA type zinc finger transcription factor family protein) [NCBI Gene 835788] {aka GATA, GATA TRANSCRIPTION FACTOR 21, GATA21, MPI10.2, MPI10_2, carbon metabolism-involved}, DREB2C (Integrase-type DNA-binding superfamily protein) [NCBI Gene 818627] {aka AtERF48, T7M7.18}, SPL (sporocyteless (SPL)) [NCBI Gene 828841] {aka F27G19.11, NOZZLE, NZZ, SPOROCYTELESS}, AT1G14550 (Peroxidase superfamily protein) [NCBI Gene 838017] {aka F14L17.33, F14L17_33}, WRKY18 (WRKY DNA-binding protein 18) [NCBI Gene 829308] {aka ARABIDOPSIS THALIANA WRKY DNA-BINDING PROTEIN 18, ATWRKY18, F11C18.16, WRKY DNA-binding protein 18}, MYB15 (myb domain protein 15) [NCBI Gene 821904] {aka ATMYB15, ATY19, MYB DOMAIN PROTEIN 15, myb domain protein 15}, WRKY22 (WRKY family transcription factor) [NCBI Gene 827896] {aka AtWRKY22, F2N1.6, F2N1_6}, ERF73 (Integrase-type DNA-binding superfamily protein) [NCBI Gene 843568] {aka AtERF73, HRE1, HYPOXIA RESPONSIVE ERF (ETHYLENE RESPONSE FACTOR) 1, T10D10.17, T10D10_17, ethylene response factor 73}, GATA19 (GATA transcription factor 19) [NCBI Gene 829814] {aka AP22.1, AP22_1, GATA transcription factor 19, HANL2, hanaba taranu like 2}, NAC038 (NAC domain containing protein 38) [NCBI Gene 816979] {aka ANAC038, ANAC039, Arabidopsis NAC domain containing protein 39, NAC domain containing protein 38, T28I24.16, T28I24_16}, DREB2B (DRE/CRT-binding protein 2B) [NCBI Gene 820273] {aka DEHYDRATION-RESPONSIVE ELEMENT BINDING PROTEIN 2, DRE/CRT-binding protein 2B, DREB2}, AUX1 (Transmembrane amino acid transporter family protein) [NCBI Gene 818390] {aka AUXIN RESISTANT 1, AtAUX1, F16M14.5, F16M14_5, MAP1, MODIFIER OF ARF7/NPH4 PHENOTYPES 1}, HSFB2A (heat shock transcription factor B2A) [NCBI Gene 836322] {aka ARABIDOPSIS THALIANA HEAT SHOCK TRANSCRIPTION FACTOR  B2A, ARABIDOPSIS THALIANA HEAT SHOCK TRANSCRIPTION FACTOR B2A, AT-HSFB2A, MTG10.7, MTG10_7, heat shock transcription factor  B2A}, WRKY6 (WRKY family transcription factor) [NCBI Gene 842527] {aka ATWRKY6, F19K23.22, F19K23_22}, HAM3 (GRAS family transcription factor) [NCBI Gene 828208] {aka ARABIDOPSIS THALIANA HAIRY MERISTEM 3, ATHAM3, F6N15.20, F6N15_20, HAIRY MERISTEM 3, LOM3}, CRF4 (cytokinin response factor 4) [NCBI Gene 828908] {aka T13J8.60, T13J8_60, cytokinin response factor 4}, DREB26 (Integrase-type DNA-binding superfamily protein) [NCBI Gene 838793] {aka T26F17.14, T26F17_14, dehydration response element-binding protein 26}, BES1 (Brassinosteroid signaling positive regulator (BZR1) family protein) [NCBI Gene 838518] {aka 107 PROTEIN, BRASSINAZOLE-RESISTANT 2, BRI1-EMS-SUPPRESSOR 1, BZR2, F18O14.7, F18O14_7}, MYB3R-5 (myb domain protein 3r-5) [NCBI Gene 831802] {aka ARABIDOPSIS THALIANA MYB DOMAIN PROTEIN 3R5, ATMYB3R5, T1E22.80, T1E22_80, myb domain protein 3r-5}, MAPKKK13 (mitogen-activated protein kinase kinase kinase 13) [NCBI Gene 837226] {aka F10K1.14, F10K1_14, mitogen-activated protein kinase kinase kinase 13}, MYB68 (myb domain protein 68) [NCBI Gene 836708] {aka ATMYB68, MPA24.14, MPA24_14, MYB DOMAIN PROTEIN 68, myb domain protein 68}, HSFA7A (heat shock transcription factor A7A) [NCBI Gene 824354] {aka ARABIDOPSIS THALIANA HEAT SHOCK TRANSCRIPTION FACTOR  A7A, ARABIDOPSIS THALIANA HEAT SHOCK TRANSCRIPTION FACTOR A7A, AT-HSFA7A, heat shock transcription factor  A7A, heat shock transcription factor A7A}, YUC1 (Flavin-binding monooxygenase family protein) [NCBI Gene 829389] {aka L23H3.20, L23H3_20, YUC, YUCCA, YUCCA 1}
- **Diseases:** fungal (MESH:D009181), Hypoxia (MESH:D000860), injury to (MESH:D014947)
- **Chemicals:** MES (MESH:C004550), BR (MESH:D060406), GA (MESH:D005708), JA (MESH:C011006), glyoxylate (MESH:C031150), GA4 (MESH:C532593), auxin (MESH:D007210), ROS (MESH:D017382), Ca (MESH:D002118), GA1 (MESH:C011258), CK (MESH:D003583), sucrose (MESH:D013395), Fe (MESH:D007501), cysteine (MESH:D003545), lipid (MESH:D008055), gibberellin (MESH:D005875), sterol (MESH:D013261), ET (MESH:C036216), glycan (MESH:D011134), urea (MESH:D014508), N (MESH:D009584), amino acid (MESH:D000596), starch (MESH:D013213), carbon (MESH:D002244), carbohydrate (MESH:D002241), SA (MESH:D012459), pyruvate (MESH:D019289), methionine (MESH:D008715), GA53 (-), nylon (MESH:D009757), ABA (MESH:D000040), DAP (MESH:C041756), sugar (MESH:D000073893), inositol phosphate (MESH:D007295), nitrate (MESH:D009566), phosphate (MESH:D010710), salt (MESH:D012492), P (MESH:D010758), 6-pentyl-2H-pyran-2-one (MESH:C008798), oxygen (MESH:D010100), Zn (MESH:D015032)
- **Species:** Botrytis cinerea (gray fruit mold, species) [taxon 40559], Neurospora crassa (species) [taxon 5141], Trichoderma reesei (species) [taxon 51453], Alternaria sect. Alternaria (section) [taxon 2499237], Solanum tuberosum (potatoes, species) [taxon 4113], Oryza sativa (Asian cultivated rice, species) [taxon 4530], Trichoderma atroviride (species) [taxon 63577], Trichoderma (genus) [taxon 5543], Pyricularia oryzae (rice blast fungus, species) [taxon 318829], Homo sapiens (human, species) [taxon 9606], Trichoderma virens (species) [taxon 29875], Solanum lycopersicum (tomato, species) [taxon 4081], Olea europaea (common olive, species) [taxon 4146], Arabidopsis thaliana (mouse-ear cress, species) [taxon 3702]

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12944188/full.md

---
Source: https://tomesphere.com/paper/PMC12944188