The Hebeloma cylindrosporum HcPT2 Pi transporter plays a key role in ectomycorrhizal symbiosis.

Through a mutualistic relationship with woody plant roots, ectomycorrhizal fungi provide growth-limiting nutrients, including inorganic phosphate (Pi), to their host. Reciprocal trades occur at the Hartig net, which is the symbiotic interface of ectomycorrhizas where the two partners are symplasmically isolated. Fungal Pi must be exported to the symbiotic interface, but the proteins facilitating this transfer are unknown. In the present study, we combined transcriptomic, microscopy, whole plant physiology, X-ray fluorescence mapping, 32 P labeling and fungal genetic approaches to unravel the role of HcPT2, a fungal Pi transporter, during the Hebeloma cylindrosporum-Pinus pinaster ectomycorrhizal association. We localized HcPT2 in the extra-radical hyphae and the Hartig net and demonstrated its determinant role for both the establishment of ectomycorrhizas and Pi allocation towards P. pinaster. We showed that the host plant induces HcPT2 expression and that the artificial overexpression of HcPT2 is sufficient to significantly enhance Pi export towards the central cylinder. Together, our results reveal that HcPT2 plays an important role in ectomycorrhizal symbiosis, affecting both Pi influx in the mycelium and efflux towards roots under the control of P. pinaster.

The ectomycorrhizal basidiomycete Hebeloma cylindrosporum undergoes early waves of transcriptional reprogramming prior to symbiotic structures differentiation.

To clarify the early molecular interaction between ectomycorrhizal partners, we performed a RNA-Seq study of transcriptome reprogramming of the basidiomycete Hebeloma cylindrosporum before symbiotic structure differentiation with Pinus pinaster. Mycorrhiza transcriptome was studied for comparison. By reference to asymbiotic mycelium, 47 and 46 genes were specifically upregulated over fivefold (p ≤ 0.05) upon rhizosphere colonization and root adhesion respectively. Other 45 were upregulated throughout the symbiotic interaction, from rhizosphere colonization to differentiated mycorrhizas, whereas 274 were specifically upregulated in mycorrhizas. Although exoproteome represents 5.6% of H. cylindrosporum proteome, 38.5% of the genes upregulated upon pre-infectious root colonization encoded extracellular proteins. The proportion decreased to 23.5% in mycorrhizas. At all studied time points, mycorrhiza-induced small secreted proteins (MiSSPs), representing potential effectors, were over-represented among upregulated genes. This was also the case for carbohydrate-active enzymes (CAZymes). Several CAZymes were upregulated at all studied stages of the interaction. Consistent with a role in fungal morphogenesis and symbiotic interface differentiation, CAZymes over-expressed before and upon root attachment targeted fungal and both fungal and plant polysaccharides respectively. Different hydrophobins were upregulated upon early root adhesion, in mycorrhizas or throughout interaction. The functional classification of genes upregulated only in mycorrhizas pointed to intense metabolic activity and nutritional exchanges.

Comparative genomics, proteomics and transcriptomics give new insight into the exoproteome of the basidiomycete Hebeloma cylindrosporum and its involvement in ectomycorrhizal symbiosis.

Extracellular proteins play crucial roles in the interaction between mycorrhizal fungi and their environment. Computational prediction and experimental detection allowed identification of 869 proteins constituting the exoproteome of Hebeloma cylindrosporum. Small secreted proteins (SSPs) and carbohydrate-active enzymes (CAZymes) were the two major classes of extracellular proteins. Twenty-eight per cent of the SSPs were secreted by free-living mycelia and five of the 10 most abundant extracellular proteins were SSPs. By contrast, 63-75% of enzymes involved in nutrient acquisition were secreted. A total of 150 extracellular protein-coding genes were differentially expressed between mycorrhizas and free-living mycelia. SSPs were the most affected. External environmental conditions also affected expression of 199 exoproteome genes in mycorrhizas. SSPs displayed different patterns of regulation in response to presence of a host plant or other environmental signals. Several of the genes most overexpressed in the presence of organic matter encoded oxidoreductases. Hebeloma cylindrosporum has not fully lost its ancestral saprotrophic capacities but rather adapted them not to harm its hosts and to use soil organic nitrogen. The complex and divergent patterns of regulation of SSPs in response to a symbiotic partner and/or organic matter suggest various roles in the biology of mycorrhizal fungi.

Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists.

To elucidate the genetic bases of mycorrhizal lifestyle evolution, we sequenced new fungal genomes, including 13 ectomycorrhizal (ECM), orchid (ORM) and ericoid (ERM) species, and five saprotrophs, which we analyzed along with other fungal genomes. Ectomycorrhizal fungi have a reduced complement of genes encoding plant cell wall-degrading enzymes (PCWDEs), as compared to their ancestral wood decayers. Nevertheless, they have retained a unique array of PCWDEs, thus suggesting that they possess diverse abilities to decompose lignocellulose. Similar functional categories of nonorthologous genes are induced in symbiosis. Of induced genes, 7-38% are orphan genes, including genes that encode secreted effector-like proteins. Convergent evolution of the mycorrhizal habit in fungi occurred via the repeated evolution of a 'symbiosis toolkit', with reduced numbers of PCWDEs and lineage-specific suites of mycorrhiza-induced genes.

A fungal conserved gene from the basidiomycete Hebeloma cylindrosporum is essential for efficient ectomycorrhiza formation.

We used Agrobacterium-mediated insertional mutagenesis to identify genes in the ectomycorrhizal fungus Hebeloma cylindrosporum that are essential for efficient mycorrhiza formation. One of the mutants presented a dramatically reduced ability to form ectomycorrhizas when grown in the presence of Pinus pinaster. It failed to form mycorrhizas in the presence of glucose at 0.5 g liter(-1), a condition favorable for mycorrhiza formation by the wild-type strain. However, it formed few mycorrhizas when glucose was replaced by fructose or when glucose concentration was increased to 1 g liter(-1). Scanning electron microscopy examination of these mycorrhizas revealed that this mutant was unable to differentiate true fungal sheath and Hartig net. Molecular analyses showed that the single-copy disrupting T-DNA was integrated 6,884 bp downstream from the start codon, of an open reading frame potentially encoding a 3,096-amino-acid-long protein. This gene, which we named HcMycE1, has orthologs in numerous fungi as well as different other eukaryotic microorganisms. RNAi inactivation of HcMycE1 in the wild-type strain also led to a mycorrhizal defect, demonstrating that the nonmycorrhizal phenotype of the mutant was due to mutagenic T-DNA integration in HcMycE1. In the wild-type strain colonizing P. pinaster roots, HcMycE1 was transiently upregulated before symbiotic structure differentiation. Together with the inability of the mutant to differentiate these structures, this suggests that HcMycE1 plays a crucial role upstream of the fungal sheath and Hartig net differentiation. This study provides the first characterization of a fungal mutant altered in mycorrhizal ability.

The impact of population structure on genomic prediction in stratified populations.

Impacts of population structure on the evaluation of genomic heritability and prediction were investigated and quantified using high-density markers in diverse panels in rice and maize. Population structure is an important factor affecting estimation of genomic heritability and assessment of genomic prediction in stratified populations. In this study, our first objective was to assess effects of population structure on estimations of genomic heritability using the diversity panels in rice and maize. Results indicate population structure explained 33 and 7.5% of genomic heritability for rice and maize, respectively, depending on traits, with the remaining heritability explained by within-subpopulation variation. Estimates of within-subpopulation heritability were higher than that derived from quantitative trait loci identified in genome-wide association studies, suggesting 65% improvement in genetic gains. The second objective was to evaluate effects of population structure on genomic prediction using cross-validation experiments. When population structure exists in both training and validation sets, correcting for population structure led to a significant decrease in accuracy with genomic prediction. In contrast, when prediction was limited to a specific subpopulation, population structure showed little effect on accuracy and within-subpopulation genetic variance dominated predictions. Finally, effects of genomic heritability on genomic prediction were investigated. Accuracies with genomic prediction increased with genomic heritability in both training and validation sets, with the former showing a slightly greater impact. In summary, our results suggest that the population structure contribution to genomic prediction varies based on prediction strategies, and is also affected by the genetic architectures of traits and populations. In practical breeding, these conclusions may be helpful to better understand and utilize the different genetic resources in genomic prediction.

Accuracy of across-environment genome-wide prediction in maize nested association mapping populations.

Most of previous empirical studies with genome-wide prediction were focused on within-environment prediction based on a single-environment (SE) model. In this study, we evaluated accuracy improvements of across-environment prediction by using genetic and residual covariance across correlated environments. Predictions with a multienvironment (ME) model were evaluated for two corn polygenic leaf structure traits, leaf length and leaf width, based on within-population (WP) and across-population (AP) experiments using a large maize nested association mapping data set consisting of 25 populations of recombinant inbred-lines. To make our study more applicable to plant breeding, two cross-validation schemes were used by evaluating accuracies of (CV1) predicting unobserved phenotypes of untested lines and (CV2) predicting unobserved phenotypes of lines that have been evaluated in some environments but not others. We concluded that (1) genome-wide prediction provided greater prediction accuracies than traditional quantitative trait loci-based prediction in both WP and AP and provided more advantages over quantitative trait loci -based prediction for WP than for AP. (2) Prediction accuracy with ME was significantly greater than that attained by SE in CV1 and CV2, and gains with ME over SE were greater in CV2 than in CV1. These gains were also greater in WP than in AP in both CV1 and CV2. (3) Gains with ME over SE attributed to genetic correlation between environments, with little effect from residual correlation. Impacts of marker density on predictions also were investigated in this study.

Evaluation of genome-wide selection efficiency in maize nested association mapping populations.

In comparison to conventional marker-assisted selection (MAS), which utilizes only a subset of genetic markers associated with a trait to predict breeding values (BVs), genome-wide selection (GWS) improves prediction accuracies by incorporating all markers into a model simultaneously. This strategy avoids risks of missing quantitative trait loci (QTL) with small effects. Here, we evaluated the accuracy of prediction for three corn flowering traits days to silking, days to anthesis, and anthesis-silking interval with GWS based on cross-validation experiments using a large data set of 25 nested association mapping populations in maize (Zea mays). We found that GWS via ridge regression-best linear unbiased prediction (RR-BLUP) gave significantly higher predictions compared to MAS utilizing composite interval mapping (CIM). The CIM method may be selected over multiple linear regression to decrease over-estimations of the efficiency of GWS over a MAS strategy. The RR-BLUP method was the preferred method for estimating marker effects in GWS with prediction accuracies comparable to or greater than BayesA and BayesB. The accuracy with RR-BLUP increased with training sample proportion, marker density, and heritability until it reached a plateau. In general, gains in accuracy with RR-BLUP over CIM increased with decreases of these factors. Compared to training sample proportion, the accuracy of prediction with RR-BLUP was relatively insensitive to marker density.

The dominant Hc.Sdh (R) carboxin-resistance gene of the ectomycorrhizal fungus Hebeloma cylindrosporum as a selectable marker for transformation.

In an attempt to get a marker gene suitable for genetical transformation of the ectomycorrhizal fungus Hebeloma cylindrosporum, the gene Hc.Sdh (R) that confers carboxin-resistance was isolated from a UV mutant of this fungus. It encodes a mutant allele of the Fe-S subunit of the succinate dehydrogenase gene that carries a single amino acid substitution known to confer carboxin-resistance. This gene was successfully used as the selective marker to transform, via Agrobacterium tumefaciens, monokaryotic and dikaryotic strains of H. cylindrosporum. We also successfully transformed hygromycin-resistant insertional mutants. Transformation yielded mitotically stable carboxin-resistant mycelia. This procedure produced transformants, the growth of which was not affected by 2 microg l(-1) carboxin, whereas wild-type strains were unable to grow in the presence of 0.1 microg l(-1) of this fungicide. This makes the carboxin-resistance cassette much more discriminating than the hygromycin-resistance one. PCR amplification and Southern blot hybridisation indicated that more than 90% of the tested carboxin-resistant mycelia contained the Hc.Sdh (R) cassette, usually as a single copy. The AGL-1 strain of A. tumefaciens was a much less efficient donor than LBA 1126; the former yielded ca. 0-30% transformation frequency, depending on fungal strain and resistance cassette used, whereas the latter yielded ca. 60-95%.

Nonmycorrhizal (myc-) mutants of Hebeloma cylindrosporum obtained through insertional mutagenesis.

Polyethylene glycol-mediated transformation of protoplasts was used as a method for insertional mutagenesis to obtain mutants of the ectomycorrhizal fungus Hebeloma cylindrosporum impaired in symbiotic ability. Following restriction enzyme-mediated integration or conventional plasmid insertion, a library of 1,725 hygromycin-resistant monokaryotic transformants was generated and screened for the symbiotic defect, using Pinus pinaster seedlings as host plants. A total of 51 transformants displaying a dramatically reduced mycorrhizal ability were identified. Among them, 29 were nonmycorrhizal (myc-), but only 10 of them had integrated one or several copies of the transforming plasmid in their genome. Light and scanning electron microscopy observations of pine roots inoculated with myc- mutants suggested that we selected mutants blocked at early stages of interaction between partners or at the stage of Hartig net formation. Myc- mutants with plasmid insertions were crossed with a compatible wild-type monokaryon and allowed to fruit. Monokaryotic progenies were obtained in three independent crosses and were analyzed for symbiotic activity and plasmid insertion. In all three progenies, a 1:1 myc-:myc+ segregation ratio was observed, suggesting that each myc- phenotype resulted from a single gene mutation. However, for none of the three mutants, the myc- phenotype segregated with any of the plasmid insertions. Our results support the idea that master genes, the products of which are essential for symbiosis establishment, do exist in ectomycorrhizal fungi.

Characterisation and expression analysis of a nitrate transporter and nitrite reductase genes, two members of a gene cluster for nitrate assimilation from the symbiotic basidiomycete Hebeloma cylindrosporum.

Symbiotic ectomycorrhizal fungi contribute to the nitrogen nutrition of their host-plants but little information is available on the molecular control of their nitrogen metabolism. We cloned and characterised genes encoding a nitrite reductase and a nitrate transporter in the ectomycorrhizal basidiomycete Hebeloma cylindrosporum. These two genes are divergently transcribed and linked to a previously cloned nitrate reductase gene, thus demonstrating that nitrate assimilation gene clusters occur in homobasidiomycetes. The nitrate transporter polypeptide (NRT2) is characterised by 12 transmembrane domains and presents both a long putative intracellular loop and a short C-terminal tail, two structural features which distinguish fungal high-affinity transporters from their plant homologues. In different wild-type genetic backgrounds, transcription of the two genes was repressed by ammonium and was strongly stimulated not only in the presence of nitrate but also in the presence of organic nitrogen sources or under nitrogen deficiency.

Agrobacterium tumefaciens-mediated transformation as a tool for insertional mutagenesis in the symbiotic ectomycorrhizal fungus Hebeloma cylindrosporum.

We transformed haploid mycelium of Hebeloma cylindrosporum via Agrobacterium tumefaciens and optimised the procedure to develop a new tool for insertional mutagenesis in this fungus. Southern blot analysis of 83 randomly selected transformants showed that they all contained plasmid inserts. Each of them showed a unique hybridisation pattern, suggesting that integration was random in the fungal genome. Sixty percent of transformants obtained in the presence of bacteria pre-treated with acetosyringone integrated a single transferred DNA copy. Thermal asymmetric interlaced polymerase chain reaction allowed us to recover the left border and the right border junctions in 85% and 15% of transformants analysed, respectively. Results show that A. tumefaciens-mediated transformation may be a powerful tool for insertional mutagenesis in H. cylindrosporum.

Isolation and characterization of nitrate reductase deficient mutants of the ectomycorrhizal fungus Hebeloma cylindrosporum.

To clarify the role of the fungal nitrate assimilation pathway in nitrate reduction by mycorrhizal plants, nitrate reductase (NR)-deficient (NR- ) mutants of the ectomycorrhizal basidiomycete Hebeloma cylindrosporum Romagnesi have been selected. These mutants were produced by u.v. mutagenesis on protoplasts originating from homokaryotic mycelia belonging to complementary mating types of this heterothallic tetrapolar species. Chlorate-resistant mutants were first selected in the presence of different nitrogen (N) sources in the culture medium. Among 1495 chlorate resistant mycelia, 30 failed to grow on nitrate and lacked a detectable NR activity. Growth tests on different N sources suggested that the NR activity of all the different mutants is specifically impaired as a result of mutations in either the gene coding for NR apoprotein or genes controlling the synthesis of the molybdenum cofactor. Furthermore, restoration of NR activity in some of the dikaryons obtained after crosses between the different mutant mycelia suggested that not all the selected mutations mapped in the same gene. Utilization of N on a NH4 15 NO3 medium was studied for two mutant strains and their corresponding wild-type homokaryons. None of the mutants could use nitrate whereas 15 N enrichment values indicated that 13-27% of N present in 13-d-old wild-type mycelia originated from nitrate. Apparently, the mutant mycelia do not compensate their inability to use nitrate by a more efficient use of ammonium. These different NR mutants still form mycorrhizas with the habitual host plant, Pinus pinaster (Ait.), making them suitable for study of the contribution of the fungal nitrate assimilation pathway to nitrate assimilation by mycorrhizal plants.