# A fusarioid fungus forms mutualistic interactions with poplar trees that resemble ectomycorrhizal symbiosis

Ningning Yang, Xiaoliang Shan, Kexuan Wang, Junkun Lu, Ying Zhu, Redman S Regina, Russell J Rodriguez, Jiajia Yao, Francis M Martin, Zhilin Yuan

► To cite this version:

Ningning Yang, Xiaoliang Shan, Kexuan Wang, Junkun Lu, Ying Zhu, et al.. A fusarioid fungus forms mutualistic interactions with poplar trees that resemble ectomycorrhizal symbiosis. IMA Fungus, 2025, 16, pp.e143240. 10.3897/imafungus.16.143240. hal-05163188

HAL Id: hal-05163188

https://hal.inrae.fr/hal-05163188v1

Submitted on 15 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A fusarioid fungus forms mutualistic interactions with poplar trees that resemble ectomycorrhizal symbiosis

Ningning Yang<sup>1,2,3,4* ID</sup>, Xiaoliang Shan<sup>4,5* ID</sup>, Kexuan Wang<sup>4* ID</sup>, Junkun Lu<sup>1,6</sup>, Ying Zhu<sup>7</sup>, Redman S. Regina<sup>8</sup>, Russell J. Rodriguez<sup>8</sup>, Jiajia Yao<sup>4</sup>, Francis M. Martin<sup>4,9 ID</sup>, Zhilin Yuan<sup>1,3,4 ID</sup>

OPEN ACCESS

Academic editor: Paola Bonfante

Received: 2 December 2024

Accepted: 28 January 2025

Published: 7 March 2025

**Citation:** Yang N, Shan X, Wang K, Lu J, Zhu Y, Regina RS, Rodriguez RJ, Yao J, Martin FM, Yuan Z (2025) A fusarioid fungus forms mutualistic interactions with poplar trees that resemble ectomycorrhizal symbiosis. IMA Fungus 16: e143240. https://doi.org/10.3897/imafungus.16.143240

Copyright: © Ningning Yang et al.
This is an open access article distributed under
terms of the Creative Commons Attribution
License (Attribution 4.0 International – CC BY 4.0).

## Abstract

Fusarium species, recognised as global priority pathogens, frequently induce severe diseases in crops; however, certain species exhibit alternative symbiotic lifestyles and are either non-pathogenic or endophytic. In this study, we characterised the mutualistic relationship between the eFp isolate of *F. pseudograminearum* and five poplar species, resulting in formation root structures reminiscent of ectomycorrhizal (ECM) symbiosis. This functional symbiosis is evidenced by enhanced plant growth, reciprocal nutrient exchange, improved nitrogen and phosphorus uptake and upregulation of root sugar transporter gene expression (*PtSweet1*). Comparative and population genomics confirmed that eFp maintains a structurally similar genome, but exhibits significant divergence from ten conspecific pathogenic isolates. Notably, eFp enhanced the growth of diverse plant lineages (*Oryza*, *Arabidopsis*, *Pinus* and non-vascular liverworts), indicating a near-complete loss of virulence. Although this specialised symbiosis has only been established *in vitro*, it holds significant value in elucidating the evolutionary track from endophytic to mycorrhizal associations.

**Key words:** Convergent evolution, endophytic fungi, lifestyle switch, root-fungal symbioses, synthetic mutualism

## Introduction

Mycorrhizal fungi, the primary components of plant root mycobiomes, have co-evolved with plants for a minimum of 400 Mya and are essential for host nutrient uptake and enhanced fitness under adverse environmental conditions (Martin et al. 2017). An extension of this mutualism has been observed

<footnote n="1">1 State Key Laboratory of Tree Genetics and Breeding, Chinese Academy of Forestry, Beijing 100091, China</footnote>

<footnote n="2">2 College of Forestry, Nanjing Forestry University, Nanjing 210037, China</footnote>

<footnote n="3">3 Zhejiang Key Laboratory of Forest Genetics and Breeding, Hangzhou 311400, China</footnote>

<footnote n="4">4 Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Hangzhou 311400, China</footnote>

<footnote n="5">5 College of Plant Protection, Nanjing Agricultural University, Nanjing 21004, China</footnote>

<footnote n="6">6 Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou 510520, China</footnote>

<footnote n="7">7 Institute of Biology, Gansu Academy of Sciences, Lanzhou 730000, China</footnote>

<footnote n="8">8 Adaptive Symbiotic Technologies, University of Washington, Seattle, WA 98195, USA</footnote>

<footnote n="9">9 INRA, UMR 1136 INRA-Université de Lorraine 'Interactions Arbres/Microorganismes', Laboratoire d'Excellence ARBRE, Centre INRA-Lorraine, Champenoux, France
Corresponding authors: Zhilin Yuan (yuanzl@caf.ac.cn); Francis M. Martin (francis.martin@inrae.fr)</footnote>

<footnote>* These authors contributed equally to this work.</footnote>
in plant-endophytic fungal associations (Rodriguez et al. 2009), which are particularly significant for non-mycorrhizal plants or plants inhabiting extreme environments (Almario et al. 2017, 2022; Yuan et al. 2021). Analogous to what we observed in other types of mutualism, genotype-by-genotype by-environment (G × G × E) interactions, within the framework of geographic mosaic theory of co-evolution (Thompson 2019), could extensively impact the degree of root-fungus mutualism stability and effectiveness. Consequently, the outcome of these associations can fluctuate along the mutualism-parasitism continuum.

Several empirical studies have emphasised that individuals within mycorrhizal fungal populations exhibit substantial heritable genetic variations, contributing to physiological changes and indirectly influencing symbiotic traits. For instance, intraspecific genetic polymorphisms, specifically copy number variations (CNVs), have been shown to affect host growth in the arbuscular mycorrhizal fungus *Rhizophagus irregularis* (Wyss et al. 2016). Notably, the implications of within-fungus variation are complex in certain root-endophytic guilds. Accumulating evidence suggests that strains of endophytic fungi isolated from healthy plants are often genetically similar to their well-known pathogenic relatives (Fesel and Zuccaro 2016; Hacquard et al. 2016). Most endophytes belonging to the *Fusarium* and *Colletotrichum* genera are frequently observed in plants and experimental evidence has confirmed their positive effects on plant fitness, indicating symbiotic lifestyle switching (Redman et al. 2001; Hiruma et al. 2016; Lofgren et al. 2018). Transcriptomic signatures reveal the downregulation of pathogenicity-related genes or disruption of toxin biosynthesis (Hacquard et al. 2016). However, traits related to saprotrophic or pathogenic lifestyles may still be retained in the genome as they are potentially detrimental when transitioning to new host plants. This notion is supported by several recent studies (Hernandez-Escribano et al. 2018; Nelson et al. 2018; Tian et al. 2020).

Relationships spanning parasitic-endophytic interactions appear to be common and such endophytes have not yet diverged sufficiently from their pathogenic relatives (Hill et al. 2022). However, we observed that two endophytic *Fusarium* species, *F. culmorum* and *F. pseudograminearum* (eFp), isolated from the dunegrass *Leymus mollis* in a high-salinity beach habitat, exhibited plant growth-promoting traits and enhanced salt tolerance (Rodriguez et al. 2008; Redman et al. 2011). Notably, during the co-cultivation of *F. pseudograminearum* and poplar trees under gnotobiotic conditions, this fungus initiated the development of root structures resembling those developed during ectomycorrhizal (ECM) symbiosis and significantly promoted plant growth.

Thus, it appears that eFp acts as a specialised mutualist and we investigated the morphological and physiological attributes associated with this association and compared the similarities and differences with genuine ECM symbioses. We posited that *F. pseudograminearum* might interact with poplar trees in an isolate-specific manner. This hypothesis was examined by integrating data from axenic synthesis trials, anatomical and isotope tracing assays and experimental inoculations. Furthermore, it was hypothesised that *F. pseudograminearum* is genetically distinct from other conspecific pathogens.
# Materials and methods

## Poplar cuttings and fungal materials

To demonstrate the degree of potential compatible interactions between *F. pseudograminearum* isolate eFp and poplar trees, five poplar clones were used in this work including *P. tomentosa*, *P. trichocarpa*, *P. alba*, *P. alba* × *P. tremula* var. *glandulosa* 84k and *P. deltoides* × *P. euramericana* NL895. Inoculation tests were performed with clones of poplar cuttings that were micropropagated *in vitro* in the laboratory. The media composition for micropropagation of each clone is provided in Suppl. material 1: table S1. eFp was originally isolated from the healthy halophyte *Leymus mollis* (Rodriguez et al. 2008) and was routinely maintained on potato dextrose agar (PDA; BD Difco).

To further evaluate the plant growth-promoting activity of eFp in a wide range of plants, we selected the following plants for inoculation: *Arabidopsis thaliana* (Col-0 ecotype), rice (Oryza sativa), liverwort (*Marchantia polymorpha*, wild-type accession Takaragaike-1), *Pinus elliottii* and *Pinus massoniana*. The method for the propagation of *M. polymorpha* was provided by Redkar et al. (2022a) and seeds of other plants were surface-infected for germination and growth. Pathogenic *F. pseudograminearum* isolate Fp8 was used as a negative control for the plant growth assay.

## Artificial synthesis of eFp-poplar associations *in vitro*

It has been proposed that nutrient conditions can significantly influence mycorrhizal formation. Generally, high glucose, low phosphorus concentrations (50 µM KH₂PO₄) as well as a low concentration of 2-(N-morpholino) ethanesulphonic acid (MES) buffer (1.25 mM, pH 5.6) was applied to promote ectomycorrhization (Guerin-Laguette et al. 2000). Given the fast growth rate of eFp, the glucose concentration was kept at a 2 g l⁻¹. In addition, organic nitrogen, composed of a mixture of acidic, neutral and aromatic amino acids [glutamine (Glu), glycine (Gly), valine (Val), leucine (Leu) and phenylalanine (Phe)], was added to the base medium at a final concentration of 3.57 mM N.

## Microscopy and anatomical description of the eFp-root association

To characterise the infection pattern and growth of eFp within poplar roots, we collected the well-developed eFp-colonised root tips 40 days after inoculation and observed their potential functional structures. First, the gross morphology of the ECM-like root tips was examined using a 3D ultra-depth stereoscopic microscope VHX-5000 (Keyence, Osaka, Japan). Next, a subset of 80% (v/v) ethanol-fixed root samples was sectioned. Briefly, the roots were washed three times with PBS (pH 7.4) and embedded in Tissue OCT-Freeze Medium (Sakura Finetek USA, Inc., Torrance, CA, USA). Transverse sections (8–10 µm in thickness) were prepared using a Thermo Cryostar NX70 freezing microtome (Thermo Fisher Scientific, Walldorf, Germany) and dual staining of the fungal cell wall (wheat germ agglutinin-alexa fluor, WGA Alexa-Fluor 488, Thermo Fisher Scientific, MA, USA) and plant cell wall (Propidium Iodide, PI, Sigma-Aldrich, USA) was performed (Becker et al. 2018). Serial sectioning was performed to avoid the sectioning of artefacts. The samples were immersed
in staining solution (20 µg ml⁻¹ PI, 10 µg ml⁻¹ WGA Alexa-Fluor 488, and 0.1% Tween 20 dissolved in PBS) for 4 h. All sections were viewed at 400–1000 × magnification under a confocal laser scanning microscope (LSM 980 with Airy-scan 2; Carl Zeiss, Jena, Germany) equipped with ZEN2 software.

Measurement of poplar growth and nutrient uptake

Prior to fungal inoculation, uniform *P. tomentosa* cuttings were transferred to con-
tainers (10 cm in diameter and 18 cm in height) with nutrient agar and the poplar
roots were inoculated with eFp by spraying 200 µl of 1 × 10⁴ spores/ml into the
medium, while the fungus-free controls were mock-inoculated with sterile wa-
ter (five replicates, five plantlets per replicate per treatment). Nutrient agar was
composed of 1/10 Murashige and Skoog (MS) medium (Murashige and Skoog
1962) and 1/5 potato broth containing glucose and sucrose (2 g l⁻¹). Plants were
kept in a growth chamber under a 12 h light/12 h dark photoperiod at 25 °C with
a light intensity of approximately 5,000 lx. Forty days post-inoculation, the plant
height, stem diameter, shoot-root ratio and biomass of the treated plants were
measured and compared. The total root length, root surface area, root volume,
number of root tips and number of root forks of intact roots were rinsed with tap
water and scanned (Epson Expression 11000XL, Epson, CA, USA) using image
analysis software (WinRHIZO Pro version, Regent Instruments, Quebec, Can-
ada). Chlorophyll content was determined as described by Peng et al. (2021).
The nitrogen content of the shoot samples was determined using the Kjeldahl
nitrogen determination method with sulphuric acid-hydrogen peroxide digestion
(Bradstreet 1954) and the phosphorus content was determined using the mo-
lybdenum-antimony resistance colourimetric method (Tiessen and Moir 1993).

<sup>14</sup>C and <sup>15</sup>N isotopic tracer experiments

To further confirm the occurrence of mutualistic interactions between eFp
and *P. tomentosa* cuttings, stable isotope labelling experiments using ¹³C and
¹⁵N were performed to trace the bidirectional transfer of C and N between the
plants and the fungus (Kaiser et al. 2015). Detailed methods for establishing
microcosms can be found in Thoen et al. (2010). In brief, 90-mm split Petri dish-
es were used to create a plant and a fungal compartment in each microcosm,
with a barrier to prevent leakage of tracers between them. Nylon mesh (8.5 cm
long × 0.8 cm wide) in the middle of split Petri dishes acted as a barrier to avoid
the growth of plant roots into the fungal compartment. The plant compartment
was filled with 25 ml of nutrient agar and the roots cuttings were restricted to
the plant compartment. The shoot was outside the dish (a small hole made
by a hot scalpel) and the cutting was inoculated with a mycelial plug before
the ¹³C pulse, whereas a 5 mm diameter mycelial plug was only placed on the
¹⁵N-labelled agar in the fungal compartment. The fungal side of the ¹³C labelling
microcosm was loaded with 25 ml PDA. Two-day ¹³C pulse-labelling, produced
through a reaction between Ba¹³CO₃ (98 atom% ¹³C, Macklin, Shanghai, China)
and lactic acid (Wang et al. 2013), was supplied to a leaf chamber. For the ¹⁵N
pulse-chasing experiment, a small 35 mm Petri dish, preventing leakage of the
tracers to the surrounding agar, was equipped with 9.5 ml of malt extract agar
(MEA) containing the ¹⁵(NH₄)₂SO₄ (99 atom% ¹⁵N Macklin, Shanghai, China) (Lu
et al. 2020). To effectively hold the small dish in place, 16 ml of MEA was added to the fungal compartment and a 5 mm diameter mycelial plug was placed in the microcosm where the fungus grew (the control was mock-inoculated with a sterile PDA plug). Petri dishes were placed inside a square Petri dish with a 240 mm-side. A sterile moistened cotton-wool ball was placed inside the Petri dish to prevent wilting. The entire system was sealed with a double layer of Parafilm and the roots were covered with aluminium foil. Microcosms were maintained in a growth chamber for 2 weeks at 25 °C under a 14 h 12,000 lx light/10 h dark cycle. Given that shoot and leaf tissues vary in isotopic composition, the aboveground parts of the plant were divided into shoots and leaves to ensure comparability. Mycelia were also collected from fungal compartments. All the samples were dried at 65 °C for 48 h and ground into a fine powder.

Samples of ¹⁵N and ¹³C were measured using an IsoPrime100 isotope ratio mass spectrometer (Isoprime Ltd., Cheadle Hulme, UK) following the manufacturer's instructions. The δ¹⁵N and δ¹³C values were calculated using the following equation:

$$ \delta^{15}\text{N} \text{ or } \delta^{13}\text{C} (\%_{\infty}) = (R_{\text{sample}} / R_{\text{standard}} - 1) \times 1000 $$

where R is the ratio of ¹⁵N/¹⁴N or ¹³C/¹²C of the sample and the standard.

## Pangenomes of F. pseudograminearum and recombination analysis

Prior to our research, seven genome assemblies of F. pseudograminearum were publicly available in databases (Suppl. material 1: table S2; Gardiner et al. 2018). The genome of the Fp8 isolate stored in our laboratory was sequenced using an Illumina platform. The genome assemblies and annotations of the nine F. pseudograminearum isolates are shown in Suppl. material 1: table S2. The construction of nine F. pseudograminearum pangenomes was described by Yuan et al. (2021). Curves describing pangenome size and the mean number of core, dispensable and new genes were fitted in R using the 'nls' function (non-linear least squares). The curve was visualised using ggplot2 software. To further investigate the evolutionary history and potential recombination within F. pseudograminearum individuals, we adopted a network representation, allowing the integration of different conflicting phylogenies. We applied the NeighborNet method (Bryant and Moulton 2004) implemented in SplitsTree v. 4.13.1 (Huson 1998) on the matrix of pairwise genetic distances, calculated from the concatenated alignment of the 5,757 orthologs.

## Comparative phylogenomics

We used eight closely-related Fusarium species and ten isolates to construct the phylogeny. Species phylogeny was estimated using 5,974 single-copy orthologs by the Maximum Likelihood (ML) method using PhyML v.3.1, with the default model HKY85. Genes and gene clusters involved in secondary metabolism were predicted using antiSMASH v.4.0.2. The dbCAN2 (http://cys.bios.niu.edu/dbCAN2) tool was used to annotate the CAZyme repertoire. We further annotated several specific gene categories of nine F. pseudograminearum, including small secreted effectors, cytochrome P450s and GPCRs (Suppl. material 1: table S3), which are putatively involved in fungal virulence. Specifically,
SECRETOL was used to predict secreted proteins (http://genomics.cicbio-gune.es/SECRETOL/) and small cysteine-rich proteins (SSCPs) with fewer than 200 amino acids and more than 4% cysteine were considered small secreted effectors. The cytochrome P450s were identified using an online tool (http://drnelson.utmem.edu/CytochromeP450.html). G-protein-coupled receptors (GPCRs) were evaluated to verify the presence of seven transmembrane helices using TMHMM v.2.0, Phobius (http://www.cbs.dtu.dk/services/TMHMM/ and http://phobius.sbc.su.se/).

## Genomic structural variations between eFp and pathogenic *F. pseudograminearum* CS3096

Single polymorphisms (SNPs) and indels between the two genomes were extracted from whole-genome alignments using nucmer in the MUMmer v.3.9.4 package. We further calculated the SNP density in 20 kb non-overlapping windows across four chromosomes. Additionally, potentially large genomic rearrangements between the two assemblies were identified using SyRI v.1.6.

## Semi-quantitative measurements of H₂O₂ production in roots colonised by eFp and Fp8

For *in situ* detection of root H₂O₂ after exposure to the two *F. pseudograminearum* isolates, the colonised and non-colonised roots were immersed in 1 mg ml⁻¹ 3, 3'-diaminobenzidine (DAB, Sigma-Aldrich) at seven time points at 30 h intervals from 30 h to 210 h, then washed and bleached with acetic acid-glycerol-ethanol (1:1:3) (v/v/v) (Fester and Hause 2005). The intensity of DAB staining was quantified using ImageJ software.

## Sweet1 gene expression pattern in populus roots when interacting with eFp and Fp8

We conducted a time course transcriptomic experiment to measure the expression pattern of *PtSweet1*, a mycorrhiza-inducible gene. Total RNA was extracted from the *P. tomentosa*-colonised and non-colonised roots using the RNAprep Pure Plant Plus kit (TIANGEN, Beijing) according to the manufacturer's instructions. RNA integrity was confirmed by 1% agarose gel electrophoresis and RNA content was quantified using NanoDrop ND-2000 (NanoDrop, Wilmington, DE, USA). Aliquots containing 150 ng of total RNA were used for first-strand cDNA synthesis in a total volume of 20 µl, containing 1 µl oligo (dT) (Aidlab, Beijing), 1 µl gDNA Remover and 4 µl of 5×TRUE Reaction Mix, according to the manufacturer's instructions. After synthesis, aliquots were stored at -80 °C.

*PtSweet1*-specific primers (forward 5'-AACAAGTCTCTATTTCTTTGTAACA-3' and reverse 5'-CCATACCAAGCAGAAAGGA-3'), described by Zhang et al. (2016), were used in this study. The Actin (forward 5'-ACCCTCCAATCCAGACACTG-3' and reverse 5'-TTGCTGACCGTATGAGCAAG-3') (Zhang et al. 2016) and EF-1α (forward 5'-AGGTCCACCAACCTTGACTG-3' and reverse 5'-AGGTCCACCAACCTTGACTG-3') (Li et al. 2024) were used as reference genes to calibrate the expression of *PtSweet1*. Quantitative PCR was performed using 10 µl Q-PCR Master mix (Thermo Fisher Scientific), 8.5 µl cDNA, 0.5 µl ROXII and 10 pmol of
each primer in a Light Cycler 480 System (Roche, Basel, Switzerland). Real-time PCR experiments were performed with at least four independent RNA samples and the threshold cycle (CT) was determined in triplicate. The relative levels of transcription were calculated by using the 2⁻^ΔΔCt^ method (Livak and Schmittgen 2001). Negative controls without cDNA were used for all PCR reactions.

**In planta** – and in vitro fungal gene expressions related to toxin biosynthesis

To further determine whether the potential virulence of eFp was attenuated during symbiotic interactions with poplars, we investigated the transcriptional changes of gene clusters involved in the biosynthesis of deoxynivalenol (DON) and zearalenone (ZEN) using RNA-seq data. The identification of these two gene clusters has been previously described (Brown et al. 2004; Nahle et al. 2021). Two fungal growth conditions were established, including *in vitro* culturing on PDA plates and a symbiotic status *in planta* (40 days after inoculation, as described above). For the *in vitro* assay, fungal colonies were grown on PDA medium covered with cellophane membranes for 2 weeks at 24 °C. The mycelium grown on each plate was gently scraped from the cellophane, collected, snap-frozen in liquid nitrogen and stored at -80 °C until further analysis. For the *in planta* assay, fungal transcripts were measured from the ECM-like root tips. Three independent biological replicates were used for each growth condition. Raw sequencing reads were trimmed using fastp (Chen et al. 2018) and mapped to the eFp reference genome using Hisat2 (Kim et al. 2015). Uniquely mapped and multimapped reads were assigned and counted using a custom pipeline that integrated featureCounts (Liao et al. 2014), mmquant (Zytnicki 2017) and custom Python and R scripts. Raw read counts were normalised using the TMM normalisation approach to obtain Counts Per Million reads (CPMs) and further normalised by gene CDS lengths to obtain Fragments Per Kilobase of exon per million reads (FPKM) values using DESeq2 (Love et al. 2014) and edgeR (Robinson and Oshlack 2010).

**Data analysis and statistics**

Prior to statistical analyses, all phenotypic and physiological datasets were tested for normality and variance homogeneity using Shapiro-Wilk's and Levene's tests, respectively. All phenotypic and physiological datasets (plant phenotypes and nutrient uptake) were subjected to Student's t-test using IBM SPSS Statistics 20 software programme (SPSS Inc., http://www.spss.com.cn). A significant Student's t-test was performed at *P* = 0.05. All data are expressed as mean values with standard deviations (SD).

**Results**

**Morphological and anatomical structures of eFp-poplar associations**

Taking advantage of an *in vitro* system, *P. tomentosa* cuttings were inoculated by eFp under sterile conditions to determine the morphological and anatomical structures developed in this association. Three-weeks post-inoculation, short lateral roots were induced from the taproots and the mycelia often completely
covered the entire lateral roots, forming root structures resembling smooth ECM root tips, which were straight with monopodial pinnate branching and becoming greyish-brown over time (Fig. 1A, B and Suppl. material 1: fig. S1). The complete development phase of this unique poplar-eFp association is shown in the Suppl. material 1: Video. The control fine roots were devoid of any apparent fungal colonisation. WGA-PI staining of cross-sections revealed sparse and thin hyphae on the root surface, resembling a patchy fungal mantle, but formation of intercellular networks of hyphae around both epidermal and cortical cells (Langer et al. 2008) (Fig. 1C-E). Colonisation of the central cylinder did not occur. More specifically, the hyphae spread across the outer and inner cortex typically formed a highly branched labyrinthine structure (Fig. 1F), which maximises the available cell surface area and optimises cell-to-cell nutrient transfer (Vaario et al. 2000). Longitudinal sections of colonised roots also showed a characteristic fungal sheath consisting of several layers of tightly packed hyphae, which further invaded the root cap and epidermal cells intercellularly (Fig. 1G).

Subsequently, we evaluated the degree of symbiotic structure development in additional poplar species across three sections: *Leuce (P. alba* and *P. alba×P. glandulosa*), *Aigeiros (P. deltoides×P. euramericana)* and *Tacamahaca (P. trichocarpa)*. We observed similar structural features in these associations irrespective of poplar species (Fig. 2 and Suppl. material 1: fig. S1), although the colonisation density, progression and depth of penetration into the root cortex differed amongst poplar species (Fig. 2E). For example, empirical evidence suggests that eFp resulted in the development of a symbiotic structure more rapidly with the 84k clone than with other poplar species. These findings suggest that eFp may facilitate the formation of ECM-like root tips in a non-host-specific manner.

## eFp improves growth and nutrient uptake of *P. tomentosa*

It is evident that eFp inoculation indeed had a significant positive effect on poplar growth (Fig. 3A-C). Specifically, 40 days after inoculation, the stem diameter, biomass, total root surface, total root volume and number of root tips and forks of eFp-inoculated plants were significantly greater than those of non-inoculated cuttings, with increases of 12.53%, 18.32%, 29.28%, 52.91%, 41.95% and 31.71%, respectively (Student's t-test, *P* < 0.05) (Fig. 3D). Moreover, the chlorophyll content of the leaves of the inoculated plants was higher than that of the control plants (*P* < 0.05) (Fig. 3E), indicating enhanced photosynthetic efficiency upon eFp colonisation. Regarding nutrition acquisition status, the total N and P concentrations in shoots from inoculated plants were significantly higher than those from control plants (*P* < 0.01 and *P* < 0.001, respectively) (Fig. 3F). These data suggest that eFp is effective in enhancing plant N and P uptake.

## Two-way flow of nutrients inferred from ¹⁵N and ¹³C isotopic tracing

The bidirectional nature of nutrient exchange has been recognised as a proxy for mutualistic interactions. Therefore, we quantified the transfer of ¹⁵N and ¹³C labelled compounds and evaluated the marker content by calculating δ¹⁵N and δ¹³C for stems, leaves and mycelia (Akroume et al. 2019). In our experimental microcosm, the mycelium traversed the barrier, whereas the roots were unable to penetrate the fungal compartment. Consequently, the uptake
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Re-synthesis experiment between eFp and *P. tomentosa*: **A**, **B** the ectomycorrhizal-like structure between eFp and *P. tomentosa* was photographed using a camera and an ultra-depth of field microscope 40 days after inoculation. The ECM-like root tips were single to pinnately branched, clumped and white to dark brown. The short and swollen lateral roots were ensheathed with eFp mycelia. The photos were adjusted for clarity and contrast as needed; **C**, **D** sectioned ECM-like root tips showing extensive hyphal colonisation along the root epidermis and cortex; **E** longitudinal section of ECM-like root tips showing a dense fungal sheath and fungal invasion of the root cap and epidermis; **F**, **G** a magnified view of highly-branched labyrinthine structures formed by irregular hyphal ramification in the inner cortex. All scale bars are shown in each panel.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Induction of ECM-like root tips in five poplars with eFp: A phylogenetic relationships between five poplar trees from the three sections were constructed with 138 orthologous genes using RAxML v. 8.2.12; B growth status of cuttings of five poplars. The genomic data of the poplar species are presented in Suppl. material 1: table S5; C-E gross morphology of ECM-like root tips and a detail of cross-sections showing extensive hyphal colonisation across root epidermis and cortex in the five poplar-eFp associations. All scale bars are shown in each panel.</figcaption></figure>

of ¹⁵N in poplar roots was entirely contingent upon its transport via fungal mycelia. The results demonstrated that only cuttings inoculated with eFp incorporated ¹⁵N in both stems and leaves (δ¹⁵N decreased progressively from the stem to the leaf), indicating that ¹⁵N can be translocated from the fungus to the aboveground plant organs. As anticipated, isotopic signatures were not detected in the control plants, thereby excluding the possibility of ¹⁵N leakage from the small Petri dish containing the isotope.

As illustrated in Fig. 3G, concurrent experiments demonstrated ¹³C enrichment from ¹³C-labelled photo-assimilates in the mycelium. The detection of ¹³C in the fungal compartment suggests that plants transfer carbon compounds to their roots through the assimilation of carbon dioxide. Collectively, these findings indicate that eFp and poplar trees are capable of establishing bidirectional nutrient flow.

## Expression patterns of the poplar *PtSweet1* during interactions with eFp and Fp8

The plant sugar transporter gene *Sweet1* is used as a molecular marker for characterising ECM establishment because its expression is often upregulated during symbiotic interactions (Neb et al. 2017; Li et al. 2024). Similar to mycorrhizal symbioses, the *PtSweet1* transcript level exhibited significant upregulation at 10 and 15 d post-inoculation compared to non-colonised roots (Fig. 3H).
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

**Figure 3.** Plant growth promotion and bidirectional exchange of nutrients conferred by eFp: **A**, **B** colony appearance of eFp on PDA; **B** spore morphology is falcate and eFp produces only large conidia ranging from 0–6 septa; **C** comparison of the growth status between uninoculated and eFp-inoculated poplar cuttings; **D**, **E** statistical analyses of plant biomass, root development and chlorophyll content increased with eFp, which was analysed using Student's t-test; **F** comparison of total nitrogen and phosphorus concentrations in the shoots of the control and fungus-inoculated plants; **G** Isotope tracing elements (¹³C and ¹⁵N) to reveal the potential bidirectional flow of nutrients between the host poplar and eFp; **H** expression pattern of the *PtSweet1* gene during symbiotic interaction with eFp.

The expression pattern of *PtSweet1* in poplar 84k demonstrated a similar trend during symbiotic interactions with eFp (Suppl. material 1: fig. S2). In contrast, during the Fp8-poplar pathogenic interaction, *PtSweet1* expression increased rapidly at 5 d post-inoculation and subsequently decreased markedly at later time points (Suppl. material 1: fig. S2). This observation aligns with findings from plant-necrotrophic fungal interactions (Breia et al. 2021).
## eFp benefits a wide range of plant lineages

To evaluate the plant growth-promoting potential of eFp, we examined addition-
al plant lineages, including *Pinus massoniana* and *P. effiottii*, rice (*Oryza sativa*)
and *Arabidopsis thaliana*, as well as the non-vascular liverwort *M. polymorpha*
(Fig. 4A). These findings indicated that eFp enhanced the growth of all test-
ed plants without exhibiting any pathogenic symptoms, suggesting a possible
complete loss of virulence. In contrast, inoculation with strain Fp8 consistently
yielded negative effects on plant growth, with varying degrees of disease symp-
toms, which were more pronounced in non-woody plants.

## Difference in poplar root H₂O₂ production during interactions with eFp and Fp8

We subsequently measured *in situ* H₂O₂ production as an indicator of plant
response to fungal invasion during a 210 h post-inoculation time course
(Fig. 4B, C). The production of H₂O₂ is a characteristic feature of the successful
recognition of infection and activation of plant defence responses (Torres and
Dangl 2005). The roots generate H₂O₂ during colonisation by both symbiotic
and pathogenic fungi. Both colonised roots exhibited transient H₂O₂ elevation,
with the peak occurring at 120 h post-inoculation. H₂O₂ levels in the roots of
Fp8-infected plants were higher than those in the eFp-infected plants across
all time points, with significant differences observed at the four time points
(*P* < 0.05). This suggests that eFp infection appears to be considerably less
effective in triggering plant defence reactions.

## eFp retains a genome structurally similar, but is highly divergent to pathogenic *F. pseudograminearum*

To identify the potential genomic traits specific to eFp, we generated a chro-
mosome-level genome assembly. Phylogenetic analysis confirmed the place-
ment of eFp in the *Graminearum* clade and showed that eFp is conspecific to
the reference isolate *F. pseudograminearum* CS3096 (Fig. 5A; Gardiner et al.
2018). Three key findings emerged from comparative and population genomic
analyses. First, the eFp genome did not exhibit significant streamlining com-
pared with that of CS3096 and other genetically proximate taxa. Upon further
comparison of a set of specific gene categories related to lifestyle and viru-
lence, no apparent gene contraction or reduction in eFp expression was ob-
served (Fig. 5B, C and Suppl. material 1: table S3). One exception was the slight
decrease in the number of small secreted effectors. Despite these conserved
trait-associated gene families, transcriptomic responses may differ significant-
ly, as exemplified by the expression of two major toxins, deoxynivalenol (DON)
and zearalenone (ZEN) in eFp, which could maintain a minimal level during the
symbiotic interaction (Suppl. material 1: fig. S3). Second, the degree of synteny
between the eFp and CS3096 genomes was very high, with 11,485 syntenic
gene pairs identified, except for a region on chromosome 2 where a re-arrange-
ment (a large-scale inversion comprising approximately ~ 1.33 Mb covering >
14.5% of chromosome 2) was found and the number of genes in this inverted
synteny block contained a total of 500 genes (Fig. 5D).
<figure><img src="image_21.png" /><figcaption>**Figure 4.** Comparison of interaction outcomes of eFp and pathogenic Fp8 with vascular and non-vascular plants: **A** for plant-eFp interactions, a stable mutualism was recorded across a wide range of plant lineages, indicating an almost complete loss of virulence of eFp. In contrast, Fp8 is detrimental to plant growth and can even cause disease symptoms in some plants; **B** macroscopic and microscopic characteristics of eFp and Fp8 on different culture media (PDA, CYA and SNA), incubated at 25 °C for 7 days; **C** root H₂O₂ production at seven time points in poplar roots inoculated with eFp and Fp8, as determined by DAB staining. Hpi, hour post-inoculation. Error bars represent SDs of the replicates. Data were analysed using one-way analysis of variance (ANOVA). Statistical significance was set at *P* < 0.05.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>**Figure 5.** Comparison of genomic differences between eFp and pathogenic *Fusarium* spp.: A time-calibrated phylogeny of nine representative pathogenic *Fusarium* species, including the outgroup *Neonectria ditissima*; B, C comparison of the number of carbohydrate-activating enzyme genes and the number of clusters of secondary metabolite genes in the genomes of eFp and the eight pathogenic *Fusarium* species; D intra-species macro-synteny and re-arrangement plot depicting chromosome comparison between CS3096 and eFp using SyRI software. Syntenic blocks are indicated in grey. CS3096 chromosomes are shown in blue and eFp chromosomes are shown in orange. The four chromosomes were proportional to size; E whole-genome comparison of CS3096 and eFp. The black histogram bars represent the reference gene sequences of the four chromosomes of CS3096; the blue, orange, grey and red histogram bars represent the linear relationship between the gene sequences of the four chromosomes of eFp; the blue broken line (repeat) and the green folded line (SNP) represent the length proportion of the repeated sequences and the single-nucleotide polymorphisms in the reference genome CS3096 with a window of 20 kb, respectively; and the black dots (lost genes) are the missing genes relative to CS3096. SNPs located in the sub-telomeric regions of chromosomes were biased toward the ends of the chromosomes.</figcaption></figure>

Despite its high synteny, we identified an extensive number of single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) when mapping CS3096 (Fig. 5E and Suppl. material 1: table S4), indicating potential evolutionary innovations. A total of 97 genes were found to be absent and 150 genes were specific to eFp, based on pan-genome analysis of nine *F. pseudograminearum* individuals (Fig. 6A). However, no Gene Ontology (GO) terms were found to be significantly enriched. The *F. pseudograminearum* pan-genome was determined to be open, as evidenced by an estimated exponent $γ > 0$ (Fig. 6B–D). Furthermore, 68 novel genes were identified in each *F. pseudograminearum* genome. Additionally, reticulate (i.e. not tree-like) evolutionary relationships within the pathogenic group were demonstrated, suggesting a clear signal of recombination (Fig. 6E). This analysis also confirmed low divergence amongst pathogenic isolates, as evidenced by the short branch lengths within the network. In contrast, eFp exhibited a long interior branch that did not cluster with the other eight pathogenic individuals (Fig. 6E), confirming that the eFp followed a distinct evolutionary trajectory.
<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

**Figure 6.** Pan-genome analysis of *F. pseudograminearum*: A proportion of genes in the core, dispensable and private genomes; B estimation of *F. pseudograminearum* core genome. Parameter *κc* is the amplitude of the exponential decay, *τc* is the decay constant and *Ω* is the best-fit value of the core genome; C estimation of the pan-genome size of *F. pseudograminearum*. The curve is a least-squares fit of the power law to the median; D estimation of new genes: The number of individual-specific genes is plotted as a function of the number of genomes sequentially added. Parameters *κs* and *τs* are equivalent to *κc* and *τc*, respectively and *tg(θ)* measures the best-fit number of specific genes; E a SplitsTree network showing relationships of all *F. pseudograminearum* isolates with no outgroup, which was estimated using the NeighborNet method using pairwise genetic distances (F84 distances obtained from the concatenated 5,757 orthologs aligned in codons).
Discussion and conclusions

ECM-like associations can arise without prior co-evolution

In the present study, we characterised the mutualistic association established
between eFp and poplar trees. Comparative genomics of hundreds of ECM
fungi have shown that they share common genomic traits, such as reduced
number of plant cell wall-degrading enzymes and secondary metabolism (SM)
clusters and the proliferation of transposable elements, common landmarks of
the switch to symbiotic lifestyle in fungi (Kohler et al. 2015; Hess et al. 2018;
Miyauchi et al. 2020; Lebreton et al. 2021). However, the evolution of these
genomic traits along the saprotrophy-to-symbiotrophy continuum remains un-
clear. Our comparative genomic analyses suggest that, despite significant diver-
gence from pathogenic *F. pseudograminearum*, no single feature distinguishes
the eFp genome from the genomes of pathogenic strains. We did not observe
erosion of genome size or gene repertoire size between pathogenic *F. pseudo-
graminearum* and eFp individuals. Similarly, the repertoire of SM clusters and
CAZyme genes was similar amongst *F. pseudograminearum* strains and only a
slight decrease in the number of secreted effectors was observed in eFp. This
finding implies that eFp shares a limited set of genomic traits with typical ECM
fungi. These findings support previous studies demonstrating that CAZymes
and effector contents are not highly relevant for exploring genetic differences
in plant-associated lifestyles in fusarioid fungi (Hill et al. 2022). Moreover, at-
tenuation of virulence is a prerequisite for endophytic development. Although
pathogenic traits were retained in the present study, such as a complete set of
toxin biosynthetic clusters in the eFp genome, their expression was significant-
ly reduced during ECM development, rendering eFp less detrimental to its hosts
(Nunes-Alves 2015; Lofgren et al. 2018; Genre et al. 2020; Redkar et al. 2022b).

Our findings suggest that mycorrhizal plants can respond to non-adapt-
ed microbes, resulting in mutualistic endophytic associations without prior
co-evolution. A similar pattern was observed in the establishment of mutual-
ism between the alga *Chlamydomonas reinhardtii*-engineered filamentous fungi
(Aanen and Bisseling 2014). More specifically, the established symbiosis might
not require intimate co-evolution of the host and microbes, but rather might be
the result of the root perception of conserved microbial signals (Nod and Myc
factors), which are able to activate the symbiosis-specific host signalling path-
way, known as the common symbiotic signalling pathway (Skiada et al. 2020).
It is reasonable to assume that eFp could induce similar symbiotic responses
in ectomycorrhizal plants. Of note, at this time, the eFp-poplar associations
have only been re-synthesised under laboratory conditions. Unfortunately, the
ECM-like structure was absent in greenhouse and field experiments, indicating
that this specialised symbiosis requires specific nutritional requirements.

eFp serves as a beneficial mycobiont

Most *Fusarium* species are plant or human pathogens (Summerell 2019). However, emerging evidence suggests that endophytic and soil non-pathogenic *Fusarium* species are also ubiquitous (Maciá-Vicente et al. 2008), indicating their multifaceted lifestyles. The *F. oxysporum* (Fo) complex serves as a compelling model for understanding pathogenic and non-pathogenic relationships (Hom and Murray
2014; Brader et al. 2017) and it has become evident that not all Fo are pathogenic (Wang et al. 2020; Guo et al. 2021). For example, Fo47 exhibits both biocontrol and plant growth-promoting capabilities. This is not always the case, however, as such mutualistic interactions may fail. Redkar et al. (2022b) observed that Fo47 is detrimental to the non-vascular liverwort *M. polymorpha*, suggesting the retention of a set of core pathogenicity factors. In contrast, eFp confers benefits to a wide range of plants, reflecting near-complete loss of virulence. The principle of reciprocal nutrient exchange (Almario et al. 2022) is applicable to eFp-poplar association and eFp functions as a robust enhancer of plant growth. Overall, our data appear to support the notion that eFp has already entered a transitional state and could readily evolve into mutualistic endophytic associations and possibly into canonical mycorrhizal associations (Selosse et al. 2018, 2022; Wu and Cox 2021).

An extensive number of structural genomic variants, including deletions and point mutations, in combination with re-arrangements can greatly affect gene function and even create pseudogenes. These factors could have contributed to the uniqueness of the eFp genome structure. Likewise, some *Colletotrichum* species can also express both pathogenic and mutualistic interactions with plants resulting from subtle genetic variations (Hacquard et al. 2016), which facilitates a shift in lifestyle from harmful pathogens to mutualistic fungi (Hiruma et al. 2023). Therefore, it is suggested that the cryptic lifestyle of *Fusarium* remains largely unexplored.

## Endophyte-ECM plant interactions versus ECM fungus-non-ECM plant interactions

There has been an impressive amount of work on both endophytic and ECM associations. Model herbaceous plants and crops have been widely used as materials for endophytic fungal inoculation. By comparison, the interactions between ECM plants and endophytic fungi and those taking place between non-ECM plants and ECM fungi are rarely addressed. In the former system, there are increasing studies on poplar root-associated fungal endophytes and their symbiotic performance (Hacquard and Schadt 2015; Vélez et al. 2017). However, the anatomical structures of these associations have not been adequately investigated (Lacercat-Didier et al. 2016; Liao et al. 2019). To the best of our knowledge, this is the first report to show the establishment of fungal structures in root tips resembling those of ECM symbiosis by an endophytic *Fusarium*. In the latter system, it has been shown that some ECM fungi, such as *Tuber* species, colonise the roots of non-ECM herbaceous plants, form stable endophytic structures (Schneider-Maunoury et al. 2020) and improve the development of lateral roots (*Pisolithus tinctorius*-**Arabidopsis thaliana**) (Reboutier et al. 2002). These types of associations are relevant for investigating evolutionary transitions.

To create a broader catalogue of potential mycorrhizal fungi, recently, Kari-
man et al. (2024) advocated the necessity of expanding mycorrhizal boundaries
because some endophytes possess core features of mycorrhizal symbioses.
Some endophytes have isotopic patterns that resemble ECM fungi (Vane-
gas-León et al. 2019). One of the most compelling examples is the dark septate
endophytic fungus *Acephala macrosclerotiorum*, which establishes ECM associ-
ations (Lukešová et al. 2015). In parallel, the term ‘brassicoid mycorrhizas’ has
been coined by Almario et al. (2022) to describe the association formed between
*Arabidopsis thaliana* and the endophytic fungus *Colletotrichum tofieldiae* strain Ct61 (Hiruma et al. 2016). Moreover, a list of *Fusarium* species has been suggested as candidate mycorrhizal fungi in orchids (Jiang et al. 2019; Sisti et al. 2019). Thus, *C. tofieldiae* Ct61 and eFp may fit the category of mycorrhizal-like endophytes and the plant-fungus symbioses initiated by them support the 'waiting room hypothesis' (Selosse et al. 2022). Thus, a special focus on endophyte-ECM plant interactions is crucial to obtain structural and functional evidence of the mycorrhizal or endophytic status of root-associated fungi from ECM plants.

In conclusion, our work points out a novel type of specialised symbiosis developed in eFp-poplar associations that may mimic ectomycorrhizas, suggesting that eFp has undergone sufficient divergence and exhibits distinct symbiotic characteristics in its genome, although a plausible mechanism has not yet been fully elucidated. Further investigation into this unique form of mutualism will be facilitated through robust genetic manipulation of *Fusarium*. Further research carried out on this symbiotic model will yield more insights into the evolutionary transitions from endophytic to mycorrhizal associations, more broadly, understanding fungal guilds with potential dual ECM/endophyte niches.

## Acknowledgements

We are grateful to Prof. Zhang Feng (Lanzhou University, China) for helpful suggestions on *Sweet1 qRT-PCR* experiments.

## Additional information

### Conflict of interest

The authors have declared that no competing interests exist.

### Ethical statement

No ethical statement was reported.

### Adherence to national and international regulations

Not applicable.

### Funding

This work was funded by the "Ling Yan" Program for Tackling Key Problems in Zhejiang Province (No. 2022C02029), Biological Breeding-Major Projects (No. 2023ZD0405806) and Open Fund of the State Key Laboratory of Tree Genetics and Breeding (TGB2018001). This work was also supported by a grant from the French National Research Agency (ANR) as part of the "Investissements d'Avenir" programme (ANR-11-LABX-0002-01, Lab of Excellence ARBRE) (to FM).

### Author contributions

ZLY conceived and designed the study. NNY, JJY and YZ conducted the physiological and molecular studies and NNY also contributed to the figure preparation. KXW prepared the sections of the ECM roots. KXW and JKL conducted isotope-tracing experiments. XLS performed evolutionary and comparative genomic analyses. FMM, RR and RJR contributed to the data interpretation and provided useful suggestions. ZLY wrote the manuscript together with the FMM. All the authors approved the final version of the manuscript.
## Author ORCIDs

Ningning Yang https://orcid.org/0009-0006-3229-1748

Xiaoliang Shan https://orcid.org/0000-0002-7553-8046

Kexuan Wang https://orcid.org/0000-0002-8512-2175

Francis M. Martin https://orcid.org/0000-0002-4737-3715

Zhilin Yuan https://orcid.org/0000-0002-5799-8401

## Data availability

The RNA-seq reads from the six libraries were deposited in SRA with accession numbers SRR23815950-SRR23815955. The accession numbers for the assembled genome of eFp are CP064753-CP064759. The accession number for the assembled genome of Fp8 is SRR29263794. All data have been submitted under BioProject accession number PRJNA643593. The Source data (Suppl. material 2) and Suppl. material 1 are in the DRYAD repository (https://doi.org/10.5061/dryad.rjdfn2znz).

## References

Aanen DK, Bisseling T (2014) The birth of cooperation. Science 345: 29-30. https://doi.org/10.1126/science.1256542

Akroume E, Maillard F, Bach C, Hossann C, Brechet C, Angeli N, Zeller B, Saint-André L, Buée M (2019) First evidences that the ectomycorrhizal fungus *Paxillus involutus* mobilizes nitrogen and carbon from saprotrophic fungus necromass. Environmental Microbiology 21: 197-208. https://doi.org/10.1111/1462-2920.14440

Almario J, Jeena G, Wunder J, Langen G, Zuccaro A, Coupland G, Bucher M (2017) Root-associated fungal microbiota of nonmycorrhizal *Arabis alpina* and its contribution to plant phosphorus nutrition. Proceedings of the National Academy of Sciences of the United States of America 114: E9403-E9412. https://doi.org/10.1073/pnas.1710455114

Almario J, Fabiańska I, Saridis G, Bucher M (2022) Unearthing the plant-microbe quid pro quo in root associations with beneficial fungi. New Phytologist 234: 1967-1976. https://doi.org/10.1111/nph.18061

Becker Y, Green KA, Scott B, Becker AM (2018) Artificial inoculation of *Epichloë festucae* into *Lolium perenne*, and visualisation of endophytic and epiphyllous fungal growth. Bio-Protocol 8: e2990. https://doi.org/10.21769/BioProtoc.2990

Brader G, Compant S, Vescio K, Mitter B, Trognitz F, Ma LJ, Sessitsch A (2017) Ecology and genomic insights into plant-pathogenic and plant-nonpathogenic endophytes. Annual Review of Phytopathology 55: 61-83. https://doi.org/10.1146/annurev-phyto-080516-035641

Bradstreet RB (1954) Kjeldahl method for organic nitrogen. Analytical Chemistry 26: 185-187. https://doi.org/10.1021/ac60085a028

Breia R, Conde A, Badim H, Fortes AM, Gerós H, Granell A (2021) Plant SWEETs: from sugar transport to plant-pathogen interaction and more unexpected physiological roles. Plant Physiology 186: 836-852. https://doi.org/10.1093/plphys/kiab127

Brown DW, Dyer RB, McCormick SP, Kendra DF, Plattner RD (2004) Functional demarcation of the *Fusarium* core trichothecene gene cluster. Fungal Genetics and Biology 41: 454-462. https://doi.org/10.1016/j.fgb.2003.12.002

Bryant D, Moulton V (2004) NeighborNet: An agglomerative method for the construction of phylogenetic networks. Molecular Biology and Evolution 21: 255-265. https://doi.org/10.1093/molbev/msh018
Chen S, Zhou Y, Chen Y, Gu J (2018) fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34: i884–i890. https://doi.org/10.1093/bioinformatics/bty560

Fesel PH, Zuccaro A (2016) Dissecting endophytic lifestyle along the parasitism/mutualism continuum in *Arabidopsis*. Current Opinion in Microbiology 20: 103–112. https://doi.org/10.1016/j.mib.2016.05.008

Fester T, Hause G (2005) Accumulation of reactive oxygen species in arbuscular mycorrhizal roots. Mycorrhiza 15: 373–379. https://doi.org/10.1007/s00572-005-0363-4

Gardiner DM, Benfield AH, Stiller J, Stephen S, Aitken K, Liu C, Kazan K (2018) A high-resolution genetic map of the cereal crown rot pathogen *Fusarium pseudograminearum* provides a near-complete genome assembly. Molecular Plant Pathology 19: 217–226. https://doi.org/10.1111/mpp.12519

Genre A, Lanfranco L, Perotto S, Bonfante P (2020) Unique and common traits in mycorrhizal symbioses. Nature Reviews Microbiology 18: 649–660. https://doi.org/10.1038/s41579-020-0402-3

Guerin-Laguette A, Plassard C, Mousain D (2000) Effects of experimental conditions on mycorrhizal relationships between *Pinus sylvestris* and *Lactarius deliciosus* and unprecedented fruit-body formation of the Saffron milk cap under controlled soilless conditions. Canadian Journal of Microbiology 46: 790–799. https://doi.org/10.1139/w00-059

Guo L, Yu H, Wang B, Vescio K, Delulio GA, Yang H, Berg A, Zhang L, Steinberg C, Edel-Hermann V, Steinberg C, Kistler HC, Ma LJ (2021) Metatranscriptomic comparison of endophytic and pathogenic *Fusarium-Arabidopsis* interactions reveals plant transcriptional plasticity. Molecular Plant-Microbe Interactions 34: 1071–1083. https://doi.org/10.1094/MPMI-03-21-0063-R

Hacquard S, Schadt CW (2015) Towards a holistic understanding of the beneficial interactions across the *Populus* microbiome. New Phytologist 205: 1424–1430. https://doi.org/10.1111/nph.13133

Hacquard S, Kracher B, Hiruma K, Münch PC, Garrido-Oter R, Thon MR, Weimann A, Damm U, Dallery JF, Hainaut M, Henrissat B, Lespinet O, Sacristán S, Ver Loren van Themaat E, Kemen E, McHardy AC, Schulze-Lefert P, O'Connell RJ (2016) Survival trade-offs in plant roots during colonization by closely related beneficial and pathogenic fungi. Nature Communications 7: 11362. https://doi.org/10.1038/ncomms11362

Hernandez-Escribano L, Iturrixta E, Elvira-Recuenco M, Berbegal M, Campos JA, Renobales G, García I, Raposa R (2018) Herbaceous plants in the understory of a pitch canker-affected *Pinus radiata* plantation are endophytically infected with *Fusarium circinatum*. Fungal Ecology 32: 65–71. https://doi.org/10.1016/j.funeco.2017.12.001

Hess J, Skrede I, Chaib De, Mares M, Hainaut M, Henrissat B, Pringle A (2018) Rapid divergence of genome architectures following the origin of an ectomycorrhizal symbiosis in the genus *Amanita*. Molecular Biology and Evolution 35: 2786–2804. https://doi.org/10.1093/molbev/msy179

Hill R, Buggs RJA, Vu DT, Gaya E (2022) Lifestyle transitions in Fusarioid fungi are frequent and lack clear genomic signatures. Molecular Biology and Evolution 39: msac085. https://doi.org/10.1093/molbev/msac085

Hiruma K, Gerlach N, Sacristán S, Nakano RT, Hacquard S, Kracher B, Neumann U, Ramírez D, Bucher M, O'Connell RJ, Schulze-Lefert P (2016) Root endophyte *Colletotrichum tofieldiae* confers plant fitness benefits that are phosphate status dependent. Cell 165: 464–474. https://doi.org/10.1016/j.cell.2016.02.028

Hiruma K, Aoki S, Takino J, Higa T, Utami YD, Shiina A, Okamoto M, Nakamura M, Kawamura N, Ohmori Y, Sugita R, Tanoi K, Sato T, Oikawa H, Minami A, Iwasaki W, Saijo Y (2023) A fungal sesquiterpene biosynthesis gene cluster critical for mutualist-
pathogen transition in *Colletotrichum tofieldiae*. Nature Communications 14: 5288.
https://doi.org/10.1038/s41467-023-40867-w

Hom EF, Murray AW (2014) Niche engineering demonstrates a latent capacity for fungal-algal mutualism. Science 345: 94-98. https://doi.org/10.1126/science.1253320

Huson DH (1998) SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14: 68-73. https://doi.org/10.1093/bioinformatics/14.1.68

Jiang J, Zhang K, Cheng S, Nie Q, Zhou SX, Chen Q, Zhou J, Zhen X, Li XT, Zhen TW, Xu M, Hsiang T, Sun Z, Zhou Y (2019) *Fusarium oxysporum* KB-3 from *Bletilla striata*: an orchid mycorrhizal fungus. Mycorrhiza 29: 531-540. https://doi.org/10.1007/s00572-019-00904-3

Kariman K, Rengel Z, Pena R, Rahimlou S, Tibbett M (2024) The necessity to expand mycorrhizal boundaries: including the fungal endophytes that possess key mycorrhizal criteria. Pedosphere 34: 520-523. https://doi.org/10.1016/j.pedsph.2024.01.004

Kaiser C, Kilburn MR, Clode PL, Fuchslueger L, Koranda M, Cliff JB, Solaiman ZM, Murphy DV (2015) Exploring the transfer of recent plant photosynthates to soil microbes: mycorrhizal pathway vs direct root exudation. New Phytologist 205: 1537-1551.
https://doi.org/10.1111/nph.13138

Kim S, Lee K, Bae SJ, Hahn JS (2015) Promoters inducible by aromatic amino acids and γ-aminobutyrate (GABA) for metabolic engineering applications in *Saccharomyces cerevisiae*. Applied Microbiology and Biotechnology 99: 2705-2714. https://doi.org/10.1007/s00253-014-6303-5

Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, Canbäck B, Choi C, Cichocki N, Clum A, Colpaert J, Copeland A, Costa MD, Doré J, Floudas D, Gay G, Girlanda M, Henrissat B, Herrmann S, Hess J, Högberg N, Johansson T, Khouja HR, LaButti K, Lahrmann U, Levasseur A, Lindquist EA, Lipzen A, Marmeisse R, Martino E, Murat C, Ngan CY, Nehls U, Plett JM, Pringle A, Ohm RA, Perotto S, Peter M, Riley R, Rineau F, Ruytinx J, Salamov A, Shah F, Sun H, Tarkka M, Tritt A, Veneault-Fourrey C, Zuccaro A (2015) Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nature Genetics 47: 410-415. https://doi.org/10.1038/ng.3223

Lacercat-Didier L, Berthelot C, Foulon J, Errard A, Martino E, Chalot M, Blaudez D (2016) New mutualistic fungal endophytes isolated from poplar roots display high metal tolerance. Mycorrhiza 26: 657-671. https://doi.org/10.1007/s00572-016-0699-y

Langer I, Krpata D, Peintner U, Wenzel WW, Schweiger P (2008) Media formulation influences in vitro ectomycorrhizal synthesis on the European aspen *Populus tremula* L. Mycorrhiza 18: 297-307. https://doi.org/10.1007/s00572-008-0182-5

Lebreton A, Zeng Q, Miyauchi S, Kohler A, Dai YC, Martin FM (2021) Evolution of the mode of nutrition in symbiotic and saprotrophic fungi in forest ecosystems. Annual Review of Ecology, Evolution, and Systematics 52: 385-404. https://doi.org/10.1146/annurev-ecolsys-012021-114902

Li R, Shi WS, Zhang P, Ma JN, Zou R, Zhang XY, Kohler AN, Martin FM, Zhang F (2024) The poplar SWEET1c glucose transporter plays a key role in the ectomycorrhizal symbiosis. New Phytologist 244: 2518-2535. https://doi.org/10.1111/nph.20183

Liao HL, Bonito G, Rojas JA, Hameed K, Wu S, Schadt CW, Labbé J, Tuskan GA, Martin F, Grigoriev IV, Vilgalys R (2019) Fungal endophytes of *Populus trichocarpa* alter host phenotype, gene Expression, and rhizobiome composition. Molecular Plant-Microbe Interactions 32: 853-864. https://doi.org/10.1094/MPMI-05-18-0133-R

Liao Y, Smyth GK, Wei S (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30: 923-930. https://doi.org/10.1093/bioinformatics/btt656
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2^ΔΔCT^ method. Methods 25: 402-408. https://doi.org/10.1006/meth.2001.1262

Lofgren LA, LeBlanc NR, Certano AK, Nachtigall J, LaBine KM, Riddle J, Broz K, Dong Y, Bethan B, Kafer CW, Kistler HC (2018) *Fusarium graminearum*: pathogen or endophyte of North American grasses? New Phytologist 217: 1203-1212. https://doi.org/10.1111/nph.14894

Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology 15: 550. https://doi.org/10.1186/s13059-014-0550-8

Lu JK, Li ZS, Yang FC, Wang SK, Liang JF, He XH (2020) Concurrent carbon and nitrogen transfer between hemiparasite *Santalum album* and two N₂-fixing hosts in a sandalwood plantation. Forest Ecology and Management 464: 118060. https://doi.org/10.1016/j.foreco.2020.118060

Lukešová T, Kohout P, Větrovský T, Vohník M (2015) The potential of dark septate endophytes to form root symbioses with ectomycorrhizal and ericoid mycorrhizal middle European forest plants. PLOS ONE 10: e0124752. https://doi.org/10.1371/journal.pone.0124752

Maciá-Vicente JG, Jansson HB, Abdullah SK, Descals E, Salinas J, Lopez-Llorca LV (2008) Fungal root endophytes from natural vegetation in Mediterranean environments with special reference to *Fusarium* spp. FEMS Microbiology Ecology 64: 90-105. https://doi.org/10.1111/j.1574-6941.2007.00443.x

Martin FM, Uroz S, Barker DG (2017) Ancestral alliances: plant mutualistic symbioses with fungi and bacteria. Science 356: eaad4501. https://doi.org/10.1126/science.aad4501

Miyauchi S, Kiss E, Kuo A, Drula E, Kohler A, Sánchez-García M, Morin E, Andreopoulos B, Barry KW, Bonito G, Buée M, Carver A, Chen C, Cichocki N, Clum A, Culley D, Crous PW, Fauchery L, Girlanda M, Hayes RD, Kéri Z, LaButti K, Lipzen A, Lombard V, Magnuson J, Maillard F, Murat C, Nolan M, Ohm RA, Pangilinan J, Pereira MF, Perotto S, Peter M, Pfister S, Riley R, Sitrit Y, Stielow JB, Szöllősi G, Žifčáková L, Štursová M, Spatafora JW, Tedersoo L, Vaario LM, Yamada A, Yan M, Wang P, Xu J, Bruns T, Baldrian P, Vilgalys R, Dunand C, Henrissat B, Grigoriev IV, Hibbett D, Nagy LG, Martin FM (2020) Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits. Nature Communications 11: 5125. https://doi.org/10.1038/s41467-020-18795-w

Murashige T, Skoog F (1962) A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiologia Plantarum 15: 473-497. https://doi.org/10.1111/j.1399-3054.1962.tb08052.x

Nahle S, El Khoury A, Atoui A (2021) Current status on the molecular biology of zearalenone: its biosynthesis and molecular detection of zearalenone producing *Fusarium* species. European Journal of Plant Pathology 159: 247-258. https://doi.org/10.1007/s10658-020-02173-9

Neb D, Das A, Hintelmann A, Nehls U (2017) Composite poplars: a novel tool for ectomycorrhizal research. Plant Cell Reports 36: 1959-1970. https://doi.org/10.1007/s00299-017-2212-2

Nelson JM, Hauser DA, Hinson R, Shaw AJ (2018) A novel experimental system using the liv-erwort *Marchantia polymorpha* and its fungal endophytes reveals diverse and context-de-pendent effects. New Phytologist 218: 1217-1232. https://doi.org/10.1111/nph.15012

Nunes-Alves C (2015) A personalized symbiosis toolkit. Nature Reviews Microbiology 13: 187. https://doi.org/10.1038/nrmicro3464
Peng L, Shan XL, Yang YZ, Wang YC, Druzhinina IS, Pan XY, Jin W, He XH, Wang XY, Zhang XG, Martin FM, Yuan ZL (2021) Facultative symbiosis with a saprotrophic soil fungus promotes potassium uptake in American sweetgum trees. Plant, Cell & Environment 44: 2793-2809. https://doi.org/10.1111/pce.14053

Reboutier D, Bianchi M, Brault M, Roux C, Dauphin A, Rona JP, Legué V, Lapeyrie F, Bou-
teau F (2002) The indolic compound hypaphorine produced by ectomycorrhizal fun-
gus interferes with auxin action and evokes early responses in nonhost *Arabidopsis
thaliana*. Molecular Plant-Microbe Interactions 15: 932–938. https://doi.org/10.1094/
MPMI.2002.15.9.932

Redkar A, Gimenez Ibanez S, Sabale M, Zechmann B, Solano R, Di Pietro A (2022b)
*Marchantia polymorpha* model reveals conserved infection mechanisms in the vas-
cular wilt fungal pathogen *Fusarium oxysporum*. New Phytologist 234: 227–241.
https://doi.org/10.1111/nph.17909

Redkar A, Sabale M, Zuccaro A, Di Pietro A (2022a) Determinants of endophytic and pathogenic lifestyle in root colonizing fungi. Current Opinion in Plant Biology 67: 102226. https://doi.org/10.1016/j.pbi.2022.102226

Redman RS, Dunigan DD, Rodriguez RJ (2001) Fungal symbiosis: from mutualism to par-
asitism, who controls the outcome, host or invader? New Phytologist 151: 705–716.
https://doi.org/10.1046/j.0028-646x.2001.00210.x

Redman RS, Kim YO, Woodward CJ, Greer C, Espino L, Doty SL, Rodriguez RJ (2011)
Increased fitness of rice plants to abiotic stress via habitat adapted symbiosis: a
strategy for mitigating impacts of climate change. PLOS ONE 6: e14823. https://doi.
org/10.1371/journal.pone.0014823

Robinson MD, Oshlack A (2010) A scaling normalization method for differential expres-
sion analysis of RNA-seq data. Genome Biology 11: R25. https://doi.org/10.1186/
gb-2010-11-3-r25

Rodriguez RJ, Henson J, Van Volkenburgh E, Hoy M, Wright L, Beckwith F, Kim YO,
Redman RS (2008) Stress tolerance in plants via habitat-adapted symbiosis. ISME J
2: 404–416. https://doi.org/10.1038/ismej.2007.106

Rodriguez RJ, White Jr JF, Arnold AE, Redman RS (2009) Fungal endophytes: diversity and functional roles. New Phytologist 182: 314–330. https://doi.org/10.1111/j.1469-8137.2009.02773.x

Schneider-Maunoury L, Deveau A, Moreno M, Todesco F, Belmondo S, Murat C, Courty PE, Jąkalski M, Selosse MA (2020) Two ectomycorrhizal truffles, *Tuber melanosporum* and *T. aestivum*, endophytically colonise roots of non-ectomycorrhizal plants in natural environments. New Phytologist 225: 2542–2556. https://doi.org/10.1111/nph.16321

Selosse MA, Schneider-Maunoury L, Martos F (2018) Time to re-think fungal ecology? Fungal ecological niches are often prejudged. New Phytologist 217: 968–972. https://doi.org/10.1111/nph.14983

Selosse MA, Petrolli R, Mujica MI, Laurent L, Perez-Lamarque B, Figura T, Bourceret A, Jacquemyn H, Li T, Gao JY, Minasiewicz J, Martos F (2022) The waiting room hypothesis revisited by orchids: were orchid mycorrhizal fungi recruited among root endophytes? Annals of Botany 129: 259–270. https://doi.org/10.1093/aob/mcab134

Sisti LS, Flores-Borges DNA, Andrade SAL, Koehler S, Bonatelli ML, Mayer JLS (2019) The role of non-mycorrhizal fungi in germination of the mycoheterotrophic orchid *Pogo-niopsis schenckii* Cogn. Frontiers in Plant Science 10: 1589. https://doi.org/10.3389/fpls.2019.01589
Skiada V, Avramidou M, Bonfante P, Genre A, Papadopoulou KK (2020) An endophytic
*Fusarium*-legume association is partially dependent on the common symbiotic
signalling pathway. *New Phytologist* 226: 1429–1444. https://doi.org/10.1111/
nph.16457

Summerell BA (2019) Resolving *Fusarium*: current status of the genus. Annual Review of Phytopathology 57: 323–339. https://doi.org/10.1146/annurev-phyto-082718-100204

Thoen E, Harder CB, Kauserud H, Botnen SS, Vik U, Taylor A, Menkis A, Skrede I (2010) In vitro evidence of root colonization suggests ecological versatility in the genus *Mycena*. *New Phytologist* 227: 601–612. https://doi.org/10.1111/nph.16545

Thompson JN (2019) The geographic mosaic of coevolution. University of Chicago Press, Chicago.

Tian B, Xie J, FuY, Cheng J, Li BO, Chen T, Zhao Y, Gao Z, Yang P, Barbetti MJ, Tyler BM, Jiang DH (2020) A cosmopolitan fungal pathogen of dicots adopts an endophytic lifestyle on cereal crops and protects them from major fungal diseases. ISME J 14: 3120–3135. https://doi.org/10.1038/s41396-020-00744-6

Tiessen H, Moir JO (1993) Characterization of available P by sequential extraction. In: Carter MR (Ed.) Soil sampling and methods of analysis. Lewis Publishers, Michigan, 75–86. https://doi.org/10.1201/9781420005271.ch25

Torres MA, Dangl JL (2005) Functions of the respiratory burst oxidase in biotic interactions, abiotic stress and development. Current Opinion in Plant Biology 8: 397–403. https://doi.org/10.1016/j.pbi.2005.05.014

Vaario LM, Gill WM, Tanaka M, Ide Y, Suzuki K (2000) Aseptic ectomycorrhizal synthesis between *Abies firma* and *Cenococcum geophilum* in artificial culture. Mycoscience 41: 395–399. https://doi.org/10.1007/BF02463953

Vanegas-León ML, Sulzbacher MA, Rinaldi AC, Roy M, Selosse MA, Neves MA (2019) Are *Trechisporales* ectomycorrhizal or non-mycorrhizal root endophytes? Mycological Progress 18: 1231–1240. https://doi.org/10.1007/s11557-019-01519-w

Vélez JM, Tschaplinski TJ, Vilgalys R, Schadt CW, Bonito G, Hameed K, Hamilton CE (2017) Characterization of a novel, ubiquitous fungal endophyte from the rhizosphere and root endosphere of *Populus* trees. Fungal Ecology 27: 78–86. https://doi.org/10.1016/j.funeco.2017.03.001

Wang B, Yu H, Jia Y, Dong Q, Steinberg C, Alabouvette C, Edel-Hermann V, Kistler HC, Ye K, Ma LJ, Guo L (2020) Chromosome-scale genome assembly of *Fusarium oxysporum* strain Fo47, a fungal endophyte and biocontrol agent. Molecular Plant-Microbe Interactions 33: 1108–1111. https://doi.org/10.1094/MPMI-05-20-0116-A

Wang QX, Ning Y, Ji DG, Li SY, Hu H (2013) *In vitro* growth and carbon utilization of the green-leaved orchid *Dendrobium officinale* are promoted by mycorrhizal associations. Botanical Studies 54: 23. https://doi.org/10.1186/1999-3110-54-23

Wu B, Cox MP (2021) Comparative genomics reveals a core gene toolbox for lifestyle transitions in *Hypocreales* fungi. Environmental Microbiology 23: 3251–3264. https://doi.org/10.1111/1462-2920.15554

Wyss T, Masclaux FG, Rosikiewicz P, Pagni M, Sanders IR (2016) Population genomics reveals that within-fungus polymorphism is common and maintained in populations of the mycorrhizal fungus *Rhizophagus irregularis*. ISME J 10: 2514–2526. https://doi.org/10.1038/ismej.2016.29

Yuan ZL, Druzhinina IS, Gibbons JG, Zhong ZH, Van de Peer Y, Rodriguez RJ, Liu ZJ, Wang XY, Wei HS, Wu Q, Wang JY, Shi GH, Cai F, Peng L, Martin FM (2021) Divergence of a genomic island leads to the evolution of melanization in a halophyte root fungus. ISME J 15: 3468–3479. https://doi.org/10.1038/s41396-021-01023-8
Zhang L, Xu XD, Wang LJ, Lu MZ (2016) Identification and functional analysis of the sucrose transporter gene *PagSUT4*. *Scientia Silvae Sinicae* 52: 21–28.

Zytnicki M (2017) mmquant: How can multi-mapping reads be counted? BMC Bioinformatics 18: 411. https://doi.org/10.1186/s12859-017-1816-4

## Supplementary material 1

**Supplementary figures, tables and video**

Authors: Ningning Yang, Xiaoliang Shan, Kexuan Wang, Junkun Lu, Ying Zhu, Redman S. Regina, Russell J. Rodriguez, Jiajia Yao, Francis M. Martin, Zhilin Yuan

Data type: zip

Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbL/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.

Link: https://doi.org/10.3897/imafungus.16.143240.suppl1

## Supplementary material 2

**Source data**

Authors: Ningning Yang, Xiaoliang Shan, Kexuan Wang, Junkun Lu, Ying Zhu, Redman S. Regina, Russell J. Rodriguez, Jiajia Yao, Francis M. Martin, Zhilin Yuan

Data type: xlsx

Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbL/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.

Link: https://doi.org/10.3897/imafungus.16.143240.suppl2