Abstract
Food production is heavily dependent on soil phosphorus (P), a non-renewable mineral resource essential for plant growth and development. Alas, about 80% is unavailable for plant uptake. Arbuscular mycorrhizal fungi may promote soil P efficient use, although the mechanistic aspects are yet to be completely understood. In this study, plant and fungal variables involved in P acquisition were investigated in maize inbred lines, differing for mycorrhizal responsiveness and low-P tolerance, when inoculated with the symbiont Rhizoglomus irregulare (synonym Rhizophagus irregularis). The expression patterns of phosphate transporter (PT) genes in extraradical and intraradical mycelium (ERM/IRM) and in mycorrhizal and control maize roots were assessed, together with plant growth responses and ERM extent and structure. The diverse maize lines differed in plant and fungal accumulation patterns of PT transcripts, ERM phenotypic traits and plant performance. Mycorrhizal plants of the low-P tolerant maize line Mo17 displayed increased expression of roots and ERM PT genes, compared with the low-P susceptible line B73, which revealed larger ERM hyphal densities and interconnectedness. ERM structural traits showed significant correlations with plant/fungal expression levels of PT genes and mycorrhizal host benefit, suggesting that both structural and functional traits are differentially involved in the regulation of P foraging capacity in mycorrhizal networks.
Similar content being viewed by others
Introduction
World population is expected to reach 9.7 billion by 2050 and 10.9 billion in 21001 with a consequent increase in food demand, whose sustainable production will represent a pivotal challenge for agriculture. Current food crop production systems strongly depend on the use of chemical fertilizers, among which phosphorus (P), essential for plant growth and development, represents a non-renewable resource, whose reserves in mined rocks are estimated to decline within 100–200 years2.
P is a key macronutrient for plants, being involved in many biological processes, such as respiration, photosynthesis, phospholipid biosynthesis and in phosphorylation-based signaling mechanisms3. As a central component of nucleic acids, membrane phospholipids and energy metabolism in the form of ATP, P represents up to 0.2% of plant dry weight4. Plants absorb P as inorganic phosphate (Pi), but, although the total amount of P in the soil may be high, more than 80% is unavailable for plant uptake, being present in organic pool or immobilized by cations and clays or lost, eventually causing water eutrophication5. Accordingly, the development of strategies for enhancing P availability from the soil, reducing external P input in agricultural ecosystems, is essential for the sustainable production of food crops. Among plant root-associated fungal microbiota, arbuscular mycorrhizal (AM) fungi (AMF) can be exploited for this purpose6,7.
AMF are an ancient group of soil fungi belonging to the phylum Glomeromycota8, which establish beneficial associations with the roots of about 80% of plant species, including the most important food crops. AMF are key organisms in agricultural ecosystems, increasing plant performance and fitness and enhancing tolerance to different stress conditions in the most diversified environments9,10. AMF deliver also important ecosystem services, such as soil structure and carbon sequestration in the soil11,12, contributing to global warming mitigation13. In this mutualistic symbiosis AMF obtain carbon compounds by plants, in exchange of mineral nutrients—mainly P, N, S, K, Ca, Cu Zn—which are absorbed and translocated through a highly interconnected extraradical mycelium (ERM) that spreads in the soil from colonized roots14. ERM, extending up to 25 m g−1 soil15, improves the absorption of those soil nutrients which are out of reach or not accessible to the host plants. Such mycorrhizal P absorption pathway is often more effective than the direct uptake pathway through epidermal root cells16. So far, only a few phosphate transporter (PTs) genes have been characterized in AMF: one in Glomus versiforme (GvPT17), Funneliformis mosseae (formerly Glomus mosseae, GmosPT18) and Gigaspora margarita (GigmPT19) and seven in Rhizoglomus irregulare (formerly Glomus intraradices, synonym Rhizophagus irregularis RiPT1-RiPT720,21). Different studies reported that such transporters are expressed in ERM and in arbusculated cells18,19,20,22,23, demonstrating both AMF ability to acquire P from soil and to control P fluxes by intraradical mycelium (IRM). The final translocation to the host is carried out by specific plant PTs of the Phosphate transporter 1 (Pht1) family, identified in different plant species, including barrel medic, potato, tomato, rice and maize24. The transcription of Pht1 genes in the different species is modulated by external P availability and was found to be regulated or induced in mycorrhizal roots25,26,27,28,29,30,31.
Among the most important plant species cultivated worldwide, maize includes many varieties selected under sufficient input of P fertilization, potentially limiting plant susceptibility to mycorrhizal colonization. Several studies stressed the importance of screening a high number of genotypes, in order to assess the genetic variability of plant susceptibility to mycorrhizal colonization32,33,34,35,36,37. In maize only two studies reported a wide characterization of inbred lines able to grow in low-P soil conditions and to respond to mycorrhizal inoculation32,38. Studies aimed at evaluating the expression of fungal and plant P transporters in different lines of the same plant species represent a major challenge for the selection of plant and fungal genotypes to be used in sustainable crop production systems. Notwithstanding, scanty information is available on the role played by host plant genotype in the regulation of fungal P uptake and ERM growth and architecture. The expression of different fungal and plant PTs in ERM and mycorrhizal roots has been evaluated in a few species (Populus trichocarpa, Sorghum bicolor, Medicago truncatula and Zea mays23,38,39), while the differential expression patterns of plant PTs in different genotypes of the same species were assessed only in maize and wheat plants38,40.
In this work, with the aim of studying how plant genotypes modulate plant and fungal variables involved in P acquisition at low external P availability, four maize inbred lines, previously reported to differ for mycorrhizal responsiveness and low-P tolerance32,38, were grown in symbiosis with the AM fungal species Rhizoglomus irregulare (basionym Glomus irregulare, synonym Rhizophagus irregularis) in an in vivo experimental system. Such microcosm experimental system was specifically devised for the assessment of ERM molecular and structural traits, meanwhile allowing the detection of mycorrhizal host benefits and plant gene expression. To this aim, the maize lines B73 and Oh43 (characterized by high mycorrhizal responsiveness and low tolerance to P deficiency), Oh40B and Mo17 (characterized by low mycorrhizal responsiveness and high tolerance to P deficiency)32,38 were evaluated for (i) the expression patterns of fungal and plant PTs in ERM, IRM and mycorrhizal roots, and (ii) phenotypic traits related to host growth and P acquisition, and to ERM extent and structure.
Results
Expression of Z. mays Pht1 (ZmPht1) genes in different maize inbred lines
Using RT-qPCR the relative expressions of eight out of 13 ZmPht1 genes were reliably quantified in control and mycorrhizal maize roots (Figs. 1 and 2). ZmPht1s transcription levels of control plants varied among lines, showing the lowest values in Mo17 (see Supplementary Tables S1 and S2 online). Compared with Mo17, control plants of the other lines showed significantly higher expression levels of transcripts of ZmPht1;3 (in all lines), of ZmPht1;8, ZmPht1;11 and ZmPht1;12 (in B73) and of ZmPht1;9 (in Oh40B) (see Supplementary Table S1 online).
In mycorrhizal plants transcription levels were also different among lines. The analysis of relative expression data of mycorrhizal maize B73 indicated significantly lower expression for ZmPht1;3, ZmPht1;8, ZmPht1;11 and ZmPht1;12 compared with controls (Fig. 1). Moreover, in mycorrhizal B73 line, root expression was lower compared with most of the other lines, and significantly lower for the transporter genes ZmPht1;4 and ZmPht1;8 (P < 0.05) (Fig. 2). On the contrary, the inbred line Mo17 showed an opposite trend, significantly upregulating half of the assayed P transporter genes under AMF colonization, compared with controls (ZmPht1;4, ZmPht1;5, ZmPht1;8 and ZmPht1;11, Fig. 1). Interestingly, mycorrhizal Mo17 plants showed significantly higher expression in 5 out of 8 transporters, compared with those of B73 (Fig. 2). Mycorrhizal plants belonging to lines Oh40B and Oh43 showed variable expression levels, depending on the transporters (Figs. 1 and 2). Normalized datasets of ZmPht1s expression in control and mycorrhizal plants were separately submitted to PCA to highlight the differences in transcription patterns among the four lines: in the absence of mycorrhizal symbiosis, expression data of ZmPht1;5, ZmPht1;8, ZmPht1;11 and ZmPht1;12 were highly related to maize line B73 (Fig. 3a). In mycorrhizal plants, an opposite relationship between ZmPht1s expression and maize lines was shown, as most transporters were more related to Mo17, compared with the other lines (Fig. 3b).
Expression of R. irregulare PT (RiPT) genes in IRM and ERM in mycorrhizal maize inbred lines
The relative expression of five RiPTs (RiPT1, RiPT2, RiPT3, RiPT5 and RiPT7) was analyzed in ERM and IRM of mycorrhizal maize plants belonging to the different inbred lines (Fig. 4, Supplementary Table S3 online). RiPTs expression in ERM growing from B73 colonized roots mostly showed lower values compared with those of ERM associated to the other maize lines. On the contrary, the ERM growing from Mo17 roots showed the highest relative expression levels detected among maize lines. As to IRM, the expression of the five RiPTs did not show significant differences among the four maize lines (see Supplementary Table S3 online).
The expression levels of individual RiPTs, mostly lower in ERM compared with IRM, were not correlated in the two fungal structures, as indicated by a not significant Mantel’s correlation test (P = 0.62), obtained when comparing the two sets of expression data. Actually, a differential expression of the five RiPTs between ERM and IRM was detected within each inbred line, with significantly lower expression of RiPT1 and RiPT2 in ERM produced in symbiosis with three out of four maize lines (Fig. 4, Supplementary Table S4 online). PCA performed by combining ERM and IRM expression datasets indicated that the variance among maize lines was mostly explained for Oh40B and Mo17 by RiPTs expression in IRM and ERM, while B73, and two replicates out of three of Oh43, were located in the opposite quadrants of the biplot (see Supplementary Fig. S1 online).
Phenotypic traits of maize inbred lines and R. irregulare mycorrhizal networks
Mycorrhizal plants belonging to the four inbred lines did not show significant differences in colonization levels, percentages of arbuscules, vesicles, hyphae, root length and shoot dry biomass (Table 1, Supplementary Table S5 online). However, B73 and Oh43 displayed the largest increments in mycorrhizal plants shoot P content compared with the other two lines, which are less dependent on mycorrhizal symbiosis (Mo17 and Oh40B) (Table 1).
ERM growing from R. irregulare-colonized roots of the inbred line B73 showed the highest hyphal density, consistently with its branching density (Fig. 5), while ERM originating from roots of the line Mo17 was characterized by the lowest values, displaying significantly reduced hyphal (− 43.6%) and branching (− 42.9%) densities, compared with those of B73 (Table 2). ERM developing from Mo17 maize plants showed the highest values of hyphal exploration capacity.
Shoot dry weight and its increases in mycorrhizal plants, compared with controls, were significantly correlated with ERM hyphal density (r = 0.58 and 0.59, respectively, P = 0.04) and with ERM length (r = 0.63, P = 0.03, and r = 0.70, P = 0.01, respectively). Moreover, in mycorrhizal maize the increase in shoot P content, with respect to controls, was negatively correlated with the ERM explored area, while P percentages showed a negative relationship with ERM density (see Supplementary Fig. S2 online). Accordingly, Mantel’s test computed to assess the relationship among data of the whole set of ERM structural traits and mycorrhizal plant benefit traits (increments in shoot dry weight and P content, in comparison with control plants) showed a significant correlation (r = 0.33, P = 0.0079).
Correlation of RiPT and ZmPht1 transcripts with plant/fungal phenotypic traits
Results of the PCA performed on the averages of each relative expression dataset (ZmPht1s, IRM and ERM RiPTs) of mycorrhizal plants, on ERM phenotypic traits measurements (explored area, hyphal density, length, branching and anastomosis densities) and on mycorrhizal maize shoot dry weight and P increments, explained 70.6% of the total variance (PC1 and PC2), and showed differential responses of the four maize inbred lines (Fig. 6). The fraction of variation in ERM phenotypic variables explained by the two axes in the PCA was high, with 52% to 79% of variance explained by ERM anastomosis, branching, hyphal densities and length, on the first axis, and 86% of variance explained by ERM explored area on the second axis. Among plant phenotypic traits, only 11% of variation in shoot P increments was explained by the first, and 48% by the second axis, whereas 37% and 45% of variation in shoot dry weight increments and mycorrhizal colonization of roots, respectively, was explained by the first axis. Analysis of expression data indicated that variation in maize and ERM PTs explained by the first axis represented 74% and 75%, respectively, while IRM PTs showed a low contribution, 16% only, on the second axis. The separation of Mo17 samples was mostly explained by high levels of both plant and fungal P transporters expression (Fig. 6), as in the PCA performed on separate datasets of expression data (Fig. 3 and Supplementary Fig. S1 online). On the contrary, high values of some ERM and plant phenotypic traits (anastomosis/branch densities and mycorrhizal colonization/shoot growth) better explained separation of data from B73 line in the PCA. While the maize line Oh40B variation was not clearly explained by any of the traits studied, PCA suggested that P increases in shoots of mycorrhizal plants explained Oh43 samples separation, mostly along the second axis (Fig. 6).
The analyses also suggested that ERM RiPTs and root ZmPht1s expression were highly correlated, as confirmed by a significant Mantel’s test (r = 0.69, P = 0.001). Mantel’s tests did not show significant correlation among data of IRM RiPTs expression with those of ZmPht1 or ERM RiPTs (both P > 0.6). Mantel’s test computed to assess the relationship among data of ERM phenotypic traits and AMF/maize PTs expression levels showed significant correlations for both datasets (r = 0.30, P = 0.032 and r = 0.55, P = 0.0004, for ERM RiPTs and ZmPht1s, respectively).
Discussion
This study showed that, at low-P availability, the four assayed mycorrhizal maize inbred lines (B73, Mo17, Oh40B, Oh43) differed in their fungal/plant PTs expression patterns, in plant/ERM phenotypic traits, and in the relationships between molecular/phenotypic features. Interestingly, mycorrhizal plants of the low-P tolerant (and poorly-mycorrhizal responsive) maize line Mo17 displayed increased expression of roots and ERM phosphate transporters, compared with the low-P susceptible (and highly-mycorrhizal responsive) line B73 which revealed larger ERM hyphal and branch densities.
The inbred lines assayed here, when grown in the absence of mycorrhizal symbiosis, were characterized by differential responses to low-P status. Previous works reported that some ZmPht1 genes were induced by P starvation, although findings were not consistent among studies carried out on different plant species25,32 and experimental systems38,41,42,43,44,45,46. Here, high expression values of ZmPht1;1, ZmPht1;3, ZmPht1;9, ZmPht1;11 and ZmPht1;12 were found in control maize lines, suggesting that plants were responding to the low-P availability. Actually, the molecular mechanisms of P acquisition and distribution are still unclear in maize and the results of ZmPht1 genes expression in genotypes differing for their low-P tolerance are often inconsistent: in this work the low-P tolerant maize line Mo17 did not show higher ZmPht1 genes transcription, as induced by P deficiency, compared with the low-P susceptible line B73. Accordingly, no differential expression of ZmPht1;1, ZmPht1;2 and ZmPht1;4 was found among maize lines differing for their low-P tolerance42, while even downregulation of ZmPht1;1 and ZmPht1;2 was detected in response to P deprivation47. Besides the activity of this transporter family, further mechanisms adopted by plants to facilitate P uptake may play a major role in the response of low-P tolerant maize lines. Indeed, upregulation of genes, whose transcripts abundance modulates P allocation between above- and belowground tissues, as reported for ZmCCD10, may facilitate P flow to shoots, thus inducing roots to absorb and/or translocate P47. Moreover, some low-P tolerant maize lines showed variations in root morphology (increased topsoil foraging, root hairs and lateral roots) that may explain differential maize responses to low-P availability42.
Here, most ZmPht1 genes were downregulated in mycorrhizal B73 roots and upregulated in mycorrhizal Mo17 line. Upregulation of ZmPht1;4, ZmPht1;5 and ZmPht1;8 genes, which were previously suggested to be induced by either P starvation or mycorrhizal colonization27,38,43, was detected in mycorrhizal Mo17 line. In this work the expression of ZmPht1;1 and ZmPht1;3 genes was downregulated in all the mycorrhizal inbred lines tested, consistently with previous findings on other maize genotypes41,48,49.
Phosphorus content of Mo17 control plants was very high, compared with the other lines, showing its ability to acquire P even at low-P levels, and not to benefit from mycorrhizas in terms of plant growth and P increases. This finding confirms also previous studies on mycorrhizal response variability among genotypes, with little mycorrhizal benefits in host lines capable of low-P tolerance42,50. Surprisingly, mycorrhizal Mo17 plants showed larger expression of 5 out of 8 ZmPht1s, compared with lines B73 and Oh43, both displaying high mycorrhizal P increases. Such results are difficult to explain when considering only gene expression data, due to the wide variation in symbiotic upregulation of individual ZmPht1 genes, depending on the host and AMF genotype and on experimental systems, and to the lack of information on how external P levels affect the transcriptional regulation of ZmPht1 genes induced by both P-starvation and mycorrhizal symbiosis.
Maize is well known for its exceptional genome structural diversity, including copy number variants and presence/absence variants, and there is growing evidence for the role of structural variation in maize adaptation and phenotypic differences51. Whole-genome, array-based, comparative genomic hybridization revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. Several hundred copy number variation sequences have been observed among the two genotypes, as well as several thousand presence/absence variation sequences that are present in B73 but not Mo1752. Moreover, de novo assembly, annotation and comparative analysis of diverse maize genomes revealed previously unknown variation in gene content, genome structure, and methylation53. Such high structural and epigenetic differences among maize lines may explain differences in P transport activity observed here. Nevertheless, present data on ERM characteristics may provide further insights into the contrasting results obtained from the different maize lines.
Actually, ERM growing from mycorrhizal roots of the four inbred lines showed large differences in its structural traits, such as hyphal density, branching density and interconnectedness. An overall correlation between ERM structural traits and both ZmPht1s and ERM RiPTs expression was found in mycorrhizal plants, suggesting that, independently on host genotype specific modulation of plant and AMF transporters, ERM characteristics play a fundamental role in symbiotic nutritional fluxes. Here, B73 and Mo17 mycorrhizal plants showed statistically similar biomass, that was associated in B73 line to larger ERM hyphal density and interconnectedness, lower plant P levels and higher mycorrhizal P increases, and in Mo17 line, to higher ZmPht1s expression levels, larger host P percentages and lower mycorrhizal P increases. Such divergent behaviors, characterized by the combination of different traits, suggested the occurrence of complex and different response strategies to low-P availability. Data similar to those obtained for B73 on shoot P increment, ZmPht1s/RiPTs transcript levels and ERM hyphal densities were obtained in Oh43 mycorrhizal plants, confirming the identification of such lines as among the best mycorrhiza-responsive ones in low-P conditions32,54,55. Interestingly, previous findings showed that Oh43 was highly responsive to R. irregulare in terms of P uptake, and that this responsiveness was better explained by ERM abundance than mycorrhizal root colonization or ZmPht1 transcripts level38. In our experimental system, ERM produced by R. irregulare in symbiosis with Oh43 line did not show the largest values in the assessed phenotypic variables or in RiPT expression, compared with the other lines. Such findings may be ascribed to the variability of experimental systems, sampling plans and monitoring endpoints adopted to characterize ERM traits and investigate the role of ERM in host P acquisition56,57,58,59.
Here, ERM structural traits showed also a significant overall correlation with mycorrhizal host benefit traits (increases in growth and P content). Previous studies correlated the abundance of extraradical hyphae with the efficiency of fungal P uptake56,60 and a recent one confirmed the relationships between extraradical hyphal length and the contribution of mycorrhizal phosphate uptake pathway to maize P acquisition38,61. These relationships are worth investigating further, in order to disclose all ERM structural variables potentially relevant to the analysis of host/fungal P acquisition mechanisms. Nevertheless, they also confirm that fungal genotypes and multiple traits, both structural and functional, are differentially involved in the regulation of P foraging capacity and translocation in AMF networks. As plant P acquisition derives from the modulation of mycorrhizal and plant P acquisition pathways, the role of ERM abundance in plant P nutrition is strictly dependent on the assortment of plant/AMF genotypes and on the availability of soil P to plants. Differences detected in R. irregulare PT expression in symbiosis with different maize lines suggest that host genotypes may provide unequal lipid/sugars fluxes to the fungus, and in turn require diverse levels of mineral nutrients from the mycorrhiza-dependent pathway, probably affecting ERM growth and the foraging reward by the symbiotic fungus.
Our findings suggest that functional variables relevant to mycelial networks properties may be influenced by fungal foraging strategies, namely interconnectedness, nutrient transport, construction cost, and resilience, which have been studied in different fungal species62. Accordingly, the high interconnectedness, functional to multidirectional transport, of the fungal network formed by R. irregulare with B73 plants, which are less tolerant to P starvation, may have facilitated P exploitation of resources over exploration, allowing the host to acquire a very high proportion of its P from the fungus, downregulating most of its ZmPht1 genes. On the contrary, the less interconnected hyphal network growing from Mo17 roots suggests a stronger role in exploration than in resources exploitation, as the plant maintained high P levels in its tissues by upregulating most of its ZmPht1 genes.
Further studies should assess the genetic basis of mycorrhizal host responses, by coupling the quantification of phenotypic/functional traits in diverse maize genotypes with quantitative trait loci mapping, as previously reported for mycorrhizal colonization in wheat36,37. Findings from such studies would boost new targeted molecular breeding strategies and the improvement of our knowledge on the genetic mechanisms regulating P acquisition and ERM development and functioning.
Methods
Experimental set-up and growth conditions
Twenty seeds for each maize inbred line were germinated in sterile conditions in the dark at 25 °C after surface sterilization in 96% ethanol for 1 min, 3.5% sodium hypochlorite solution for 30 min, followed by several washing steps with sterile water. After 7–10 days, a half of maize plantlets for each line were inoculated with Rhizoglomus irregulare (Błaszk., Wubet, Renker & Buscot) Sieverd., G.A. Silva & Oehl (University of Pisa, isolate IMA6). The inoculum was obtained from Medicago sativa pot-culture soil after wet sieving through a 100-μm-mesh size sieve. Each root system was wrapped in a nylon net (41 μm mesh) in order to obtain a flat mesh pocket limiting root growth in the volume inside the nylon net (about 65 cm3). Wrapped inoculated plantlets and non-inoculated controls were transferred into perforated 14 cm Petri dishes, allowing the growth of maize shoots outside the system. The root systems were covered with sterile grit and moistened with 10 ml distilled water. After sealing with Parafilm, the Petri dishes were wrapped with aluminum foil, placed into sun-transparent bags and maintained in a growth chamber at 25 °C, with 16 h of light per day. The plants were watered weekly with 3 ml of modified Hoagland solution containing ½ strength of the standard concentration of KH2PO4.
After 7 weeks’ growth the nylon net was removed and roots were washed with distilled water, placed between two 13 cm diameter membranes (cellulose esters, 0.45 μm pore diameter size) to allow the growth of extraradical mycelial networks63, which were placed into new perforated Petri dishes containing sterile quartz grit and treated as described above. Mycorrhizal and control plants were maintained in the latter system and watered weekly with 3 ml of the modified Hoagland solution described above. On average, during the experiment 37.5 µM KH2PO4 were supplied to each plant in the growth solution.
The root sandwich systems were opened after 3 weeks’ growth (see Supplementary Fig. S3, S4 online), three biological replicates per each line were used for molecular analyses and three further biological replicates were used for ERM structural traits analyses, IRM assessment, total and mycorrhizal root length, shoot biomass production and shoot P content. Plant tissues P percentages were measured after sulphuric/perchloric acid digestion using the photometric method64.
Harvest and total RNA extraction
The ERM spreading from the roots and growing on the surfaces of the membranes was harvested on ice-cold sterile distilled water using two forceps under a dissecting microscope. Colonized roots were selected by visualizing autofluorescence of arbuscules under blue light, using an inverted fluorescence microscope, cut from the rest of the plant and collected on ice. Collected ERM and colonized roots, after removing excess water with a filter paper, were immediately weighed, and transferred in liquid nitrogen. Different aliquots of ERM (1.5 to 5 mg of fresh biomass) and of colonized and control roots (100 mg of fresh biomass) per each biological replicate were collected and stored at − 80 °C until RNA extraction. RNA from ERM was extracted using the MasterPure Complete DNA and RNA Purification Kit (Epicenter), according to the manufacturer’s instructions and eluted in TE with RiboGuard 40 U/µl. RNA from colonized and control roots was extracted using NucleoSpin RNA Plus kit (Macherey–Nagel) according to the manufacturer’s instructions and eluted in RNase-free water. RNA was quantified with the Qubit RNA HS Assay kit and stored at − 80 °C.
Quantitative gene expression analysis of Z. mays P transporter (ZmPht1) genes in maize roots
RNA extracted from colonized and non-colonized (control) maize roots was used for gene expression analysis of the Z. mays P transporters (ZmPht1). RNA extracts were DNase treated and cDNA from three biological replicates for each line was synthesized using High-Capacity cDNA Reverse Transcription Kits (Applied Biosystems, Italy), according to the manufacturer’s protocol. The synthesized cDNA was quantified with the Qubit RNA HS Assay kit and stored at − 20 °C. The PCR program was conducted in an Eppendorf Mastercycler thermal cycler according to these parameters: 95 °C for 3 min and then 35 cycles of 95 °C for 40 s, 57 °C for 40 s and 72 °C for 40 s. PCR reactions were performed in a final volume of 15 µl, containing 10 × GoTaq Reaction Buffer, 100 × BSA, 50 mM MgCl2, 10 µM of dNTPs, 10 μM of each primer, 5 unit of GoTaq DNA Polymerase, and 5 ng of cDNA. The PCR products were checked on 2.5% TBE agarose gel and visualized by EuroSafe staining.
The qRT-PCR was carried out in a final volume of 25 μl containing 12.5 μl of iQ SYBR Green Supermix 2 × (Bio-Rad, Italy), 2 μM of each primer and 5 ng cDNA template using a CFX Connect real-time PCR Detection System (Bio-Rad). The qRT-PCR program was used as follows: 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min and to calibrate gene expression data, a maize housekeeping gene ZmActin143, was used as an internal control. The relative transcript level was calculated as 2−ΔΔCT and the reactions were performed on three biological samples with two technical replicates. The specificity of qRT-PCR amplification was assessed using a melting curve analysis of 60 to 95 °C after the final PCR cycle65. The primers used were those reported by Liu et al.43.
Quantitative gene expression analysis of R. irregulare P transporter (RiPT) genes in ERM/IRM
To remove genomic DNA traces, the ERM and root RNA samples were treated with DNase, according to the manufacturer’s instructions. cDNAs (100 ng each sample) were synthesized using the iScript Select cDNA Synthesis kit (Bio-Rad) and stored at − 20 °C. For R. irregulare P transporters (RiPTs) gene expression analysis, the following primers were used: RiPT1, RiPT3, RiPT5 and RiPT739, while the pair RiPT2 F (5′-TTGGTGTCCTTACTGCACTC-3′) and RiPT2 R (5′- TGTTGTTGTGTTTGGACCGA-3′), was designed in Primer 3 (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) and tested in Amplify 3.1 (http://engels.genetics.wisc.edu/amplify) based on the partial RNA sequence of the R. irregulare phosphate transporter 2 (RiPT2)66.
Quantitative real time polymerase chain reaction (qRT-PCR) was performed using a CFX Connect real-time PCR Detection System. Reactions were carried out in a final volume of 20 μl containing 10 μl of iQ SYBR Green Supermix 2 × , 0.2 μM of each primer and 5 ng cDNA template. The PCR thermal program consisted in an initial incubation at 95 °C for 3 min, followed by 40 cycles of 95 °C for 30 s, 60 °C for 1 min and 72 °C for 30 s, where the fluorescence signal was measured21. At the end of each cycle, a melting curve was recorded (from 60 to 95 °C), to assess the specificity of the amplification product and therefore to exclude the formation of non-specific products by the primers65. The efficiency of the primers was tested by performing qRT-PCR on different cDNA dilutions. All reactions were performed at least two times on three biological and two technical replicates. All gene expression data were normalized to the expression of the elongation factor 1-α (ef1-α) using primers ef1-α F and ef1-α57. The relative levels of transporters transcription were calculated by using the 2−ΔCT method (ratio (reference/target) = 2CT(reference) − CT(target)), a variation of the 2−ΔΔCT method67.
Phenotypic traits of maize inbred lines and R. irregulare mycorrhizal networks
Three biological replicates for each maize line were used to measure shoot P percentages and shoot biomass of plants grown in the sandwich systems and to determine AMF root colonization and ERM structural traits. The shoot of each plant was cut and oven-dried at 50 °C until constant weight. Roots and ERM growing on the membranes were stained for the localization of succinate dehydrogenase (SDH) activity. After the SDH test, the same membranes were stained with Trypan blue in lactic acid (0.05%). These staining methods allowed the assessment of mean ERM density (expressed as mm of hyphal length per mm2 of membrane surface) and the area explored by ERM (expressed in mm2), under the dissecting microscope68. For each membrane, mean ERM density was measured with a grid eyepiece by evaluating the length of hyphae in 15 areas of 9.8 mm2, using the gridline intersect method, while the area explored by ERM was estimated with a transparent 5-mm-square grid. Total ERM length was calculated as the product of mean ERM density and area obtained from each replicate membrane. For each treatment and replicate membrane, five regions (32 × 20 mm each) covered by ERM were cut, mounted on microscope slides and examined under a Reichert-Jung Polyvar microscope to assess the number of hyphal contacts, anastomoses and branches68, which were referred to the hyphal length and to the observed area.
Mycorrhizal colonization was evaluated by clearing the roots with 10% KOH in a water bath at 80 °C for 15 min and subsequently staining them with Trypan blue in lactic acid (0.05%) after a treatment in 2% HCl for 10 min. Mycorrhizal and total root length, along with percentage of arbuscules, vesicles and hyphae, were estimated under a dissecting microscope using the grid line intersect method69.
Data analysis
Data of fungal and plant P transporter expression and data of mycelial growth and viability were logarithmically transformed, when needed, to fulfil ANOVA assumptions and subjected to one-way ANOVA with inbred line as the source of variation, followed by Tukey’s HSD post hoc test. Welch’s test for Unequal Variances, followed by Tamhane's post hoc test, was used when the variance homogeneity assumption was not satisfied. A probability level of P < 0.05 was considered for all tests. Percentage data were analyzed after arcsine transformation. Gene expression data analysed in this work are reported in Supplementary Tables S5 and S6 online. Multivariate (Mantel’s test, permutations 9999) and linear regression analyses were performed to determine the correlations among gene expression data and between them and phenotypic traits of both plants and ERM. Principal Component Analyses (PCA) were computed to assess the relationships between physiological and molecular data collected from four Z. mays inbred lines in symbiosis or not with R. irregulare. Statistical analyses and graphics were generated using the software tools R (package pheatmap), SPSS (version 23), PAST (version 4.03) and CANOCO (version 5.0).
Additional statements
The maize genotypes employed in this study are public lines. They were obtained from the Maize Genetics Cooperation Stock Center, U.S. National Plant Germplasm System (https://www.maizegdb.org/data_center/stock) and maintained by sibling at the Department of Sustainable Crop Production, Piacenza, Italy. All the experimental research studies on plants carried out in this work, including the collection of plant material, are complied with relevant institutional, national, and international guidelines and legislation.
Data availability
All data generated or analysed during this study are included in this published article (and its Supplementary Information files).
References
United Nations, Department of Economic and Social Affairs, Population Division World Population Prospects 2019: Highlights. ST/ESA/SER.A/423 (2019).
Cordell, D. & Whites, S. Life’s bottleneck: Sustaining the world’s phosphorus for a food secure future. Ann. Rev. Environ. Res. 39, 161–188 (2014).
Vance, C. P., Uhde-Stone, C. & Allan, D. L. Phosphorus acquisition and use: Critical adaptations by plants for securing a nonrenewable resource. New Phytol. 157, 423–447 (2003).
Schachtman, D. P., Reid, R. J. & Ayling, S. M. Phosphorus uptake by plants: From soil to cell. Plant Physiol. 116, 447–453 (1998).
Sbrana, C. et al. Mycorrhizal symbionts and associated bacteria: Potent allies to improve plant phosphorus availability and food security. Front. Microbiol. 12, 797381–797381. https://doi.org/10.3389/fmicb.2021.797381 (2021).
Schröder, J. J., Smit, A. L., Cordell, D. & Rosemarin, A. Improved phosphorus use efficiency in agriculture: A key requirement for its sustainable use. Chemosphere 84, 822–831 (2011).
Rouphael, Y. et al. Arbuscular mycorrhizal fungi act as biostimulants in horticultural crops. Sci. Hort. 196, 91–108 (2015).
Tedersoo, L., Tooming-Klunderud, A. & Anslan, S. PacBio metabarcoding of fungi and other eukaryotes: Biases and perspectives. New Phytol. 217, 1370–1385 (2018).
Begum, N. et al. Role of arbuscular mycorrhizal fungi in plant growth regulation: Implications in abiotic stress tolerance. Front. Plant Sci. 10, 1068. https://doi.org/10.3389/fpls.2019.01068 (2019).
Evelin, H., Devi, T. S., Gupta, S. & Kapoor, R. Mitigation of salinity stress in plants by arbuscular mycorrhizal symbiosis: Current understanding and new challenges. Front. Plant Sci. 10, 470. https://doi.org/10.3389/fpls.2019.00470 (2019).
Gianinazzi, S. et al. Agroecology: The key role of arbuscular mycorrhizas in ecosystem services. Mycorrhiza 20, 519–530 (2010).
Rillig, M. C. et al. Towards an integrated mycorrhizal technology: Harnessing mycorrhiza for sustainable intensification in agriculture. Front. Plant Sci. 7, 1625. https://doi.org/10.3389/fpls.2016.01625 (2016).
Verbruggen, E., Struyf, E. & Vicca, S. Can arbuscular mycorrhizal fungi speed up carbon sequestration by enhanced weathering?. Plants People Planet 3, 445–453. https://doi.org/10.1002/ppp3.10179 (2021).
Smith, S. E. & Read, D. J. Mycorrhizal Symbiosis (Academic Press, 2008).
Mikkelsen, B. L., Rosendahl, S. & Jakobsen, I. Underground resource allocation between individual networks of mycorrhizal fungi. New Phytol. 180, 890–898 (2008).
Wipf, D., Krajinski, F., van Tuinen, D., Recorbet, G. & Courty, P. E. Trading on the arbuscular mycorrhiza market: From arbuscule to common mycorrhizal networks. New Phytol. 223, 1127–1142 (2019).
Harrison, M. J. & van Buuren, M. L. A phosphate transporter from the mycorrhizal fungus Glomus versiforme. Nature 378, 626–629 (1995).
Benedetto, A., Magurno, F., Bonfante, P. & Lanfranco, L. Expression profiles of a phosphate transporter gene (GmosPT) from the endomycorrhizal fungus Glomus mosseae. Mycorrhiza 15, 620–627 (2005).
Xie, X. et al. Arbuscular mycorrhizal symbiosis requires a phosphate transceptor in the Gigaspora margarita fungal symbiont. Mol. Plant 9, 1583–1608 (2016).
Maldonado-Mendoza, I. E., Dewbre, G. R. & Harrison, M. J. A phosphate transporter gene from the extra-radical mycelium of an arbuscular mycorrhizal fungus Glomus intraradices is regulated in response to phosphate in the environment. Mol. Plant-Microbe Interact. 14, 1140–1148 (2001).
Walder, F., Boller, T., Wiemken, A. & Courty, P. E. Regulation of plants’ phosphate uptake in common mycorrhizal networks: Role of intraradical fungal phosphate transporters. Plant Signal Behav. 11, e1131372. https://doi.org/10.1080/15592324.2015.1131372 (2016).
Balestrini, R., Gómez-Ariza, J., Lanfranco, L. & Bonfante, P. Laser microdissection reveals that transcripts for five plant and one fungal phosphate transporter genes are contemporaneously present in arbusculated cells. Mol. Plant Microbe Interact. 20, 1055–1062 (2007).
Fiorilli, V., Lanfranco, L. & Bonfante, P. The expression of GintPT, the phosphate transporter of Rhizophagus irregularis, depends on the symbiotic status and phosphate availability. Planta 237, 1267–1277 (2013).
Bucher, M., Rausch, C. & Daram, P. Molecular and biochemical mechanisms of phosphorus uptake into plants. J. Plant Nutr. Soil Sci. 164, 209–217 (2001).
Rausch, C. et al. A phosphate transporter expressed in arbuscule-containing cells in potato. Nature 414, 462–465 (2001).
Harrison, M. J., Dewbre, G. R. & Liu, J. A phosphate transporter from Medicago truncatula involved in the acquisition of phosphate released by arbuscular mycorrhizal fungi. Plant Cell 14, 2413–2429 (2002).
Nagy, R. et al. The characterization of novel mycorrhiza-specific phosphate transporters from Lycopersicon esculentum and Solanum tuberosum uncovers functional redundancy in symbiotic phosphate transport in solanaceous species. Plant J. 42, 236–250 (2005).
Yang, S. Y. et al. Nonredundant regulation of rice arbuscular mycorrhizal symbiosis by two members of the PHOSPHATE TRANSPORTER1 gene family. Plant Cell 24, 4236–4251 (2012).
Xie, X. et al. Functional analysis of the novel mycorrhiza-specific phosphate transporter AsPT1 and PHT1 family from Astragalus sinicus during the arbuscular mycorrhizal symbiosis. New Phytol. 198, 836–852 (2013).
Smith, S. E., Jakobsen, I., Grønlund, M. & Smith, F. A. Roles of arbuscular mycorrhizas in plant phosphorus nutrition: Interactions between pathways of phosphorus uptake in arbuscular mycorrhizal roots have important implications for understanding and manipulating plant phosphorus acquisition. Plant Physiol. 156, 1050–1057 (2011).
Wang, Y., Wang, F., Lu, H., Liu, Y. & Mao, C. Phosphate uptake and transport in plants: An elaborate regulatory system. Plant Cell Physiol. 62, 564–572 (2021).
Kaeppler, S. M. et al. Variation among maize inbred lines and detection of quantitative trait loci for growth at low phosphorus and responsiveness to arbuscular mycorrhizal fungi. Crop Sci. 40, 358–364 (2000).
An, G. H. et al. How does arbuscular mycorrhizal colonization vary with host plant genotype? An example based on maize (Zea mays) germplasms. Plant Soil 327, 441–453 (2010).
Turrini, A. et al. Large variation in mycorrhizal colonization among wild accessions, cultivars, and inbreds of sunflower (Helianthus annuus L.). Euphytica 207, 331–342 (2016).
Wang, X., Zhao, S. & Bücking, H. Arbuscular mycorrhizal growth responses are fungal specific but do not differ between soybean genotypes with different phosphate efficiency. Ann. Bot. 118, 11–21 (2016).
De Vita, P. et al. Genetic markers associated to arbuscular mycorrhizal colonization in durum wheat. Sci. Rep. 8, 1–12. https://doi.org/10.1080/15592324.2015.1131372 (2018).
Ganugi, P. et al. Genetic variability assessment of 127 Triticum turgidum L. accessions for mycorrhizal susceptibility-related traits detection. Sci Rep. 11, 1. https://doi.org/10.1038/s41598-021-92837-1 (2021).
Sawers, R. J. et al. Phosphorus acquisition efficiency in arbuscular mycorrhizal maize is correlated with the abundance of root-external hyphae and the accumulation of transcripts encoding PHT1 phosphate transporters. New Phytol. 214, 632–643 (2017).
Calabrese, S. et al. Imbalanced regulation of fungal nutrient transports according to phosphate availability in a symbiocosm formed by poplar, sorghum, and Rhizophagus irregularis. Front. Plant Sci. 10, 1617. https://doi.org/10.3389/fpls.2019.01617 (2019).
Tian, H. et al. Influence of nutrient signals and carbon allocation on the expression of phosphate and nitrogen transporter genes in winter wheat (Triticum aestivum L.) roots colonized by arbuscular mycorrhizal fungi. PLoS ONE 12, e0172154. https://doi.org/10.1371/journal.pone.0172154 (2017).
Tian, H., Drijber, R. A., Zhang, J. L. & Li, X. L. Impact of long-term nitrogen fertilization and rotation with soybean on the diversity and phosphorus metabolism of indigenous arbuscular mycorrhizal fungi within the roots of maize (Zea mays L.). Agric. Ecosyst. Environ. 164, 3–61 (2013).
Gu, R. et al. Enhancing phosphorus uptake efficiency through QTL-based selection for root system architecture in maize. J. Genet. Genom. 43, 663–672 (2016).
Liu, F. et al. Systematic identification, evolution and expression analysis of the Zea mays PHT1 gene family reveals several new members involved in root colonization by arbuscular mycorrhizal fungi. Int. J. Mol. Sci. 17, 930 (2016).
Liu, F. et al. Identification and functional characterization of a maize phosphate transporter induced by mycorrhiza formation. Plant Cell Physiol. 59, 1683–1694 (2018).
Wang, F. et al. Maize ZmPT7 regulates Pi uptake and redistribution which is modulated by phosphorylation. Plant Biotech. J. 18, 2406–2419 (2020).
Sun, Y. & Zheng, H. Divergent molecular and physiological response of two maize breeding lines under phosphate deficiency. Plant Mol. Biol. Rep. 40, 197–207 (2021).
Zhong, Y. et al. ZmCCD10a encodes a distinct type of carotenoid cleavage dioxygenase and enhances plant tolerance to low phosphate. Plant Physiol. 184, 374–392 (2020).
Wright, D. P., Scholes, J. D., Read, D. J. & Rolfe, S. A. European and African maize cultivars differ in their physiological and molecular responses to mycorrhizal infection. New Phytol. 167, 881–896 (2005).
Nagy, R. et al. Differential regulation of five Pht1 phosphate transporters from maize (Zea mays L.). Plant Biol. 8, 186–197 (2006).
Zhu, Y. G., Smith, S. E. & Smith, F. A. Zinc (Zn)-phosphorus (P) interactions in two cultivars of spring wheat (Triticum aestivum L.) differing in P uptake efficiency. Ann. Bot. 88, 941–945 (2001).
Darracq, A. et al. Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants. BMC Genom. 19, 119 (2018).
Sun, S. et al. Extensive intraspecific gene order and gene structural variations between Mo17 and other maize genomes. Nat. Genet. 50, 1289–1295 (2018).
Hufford, M. B. et al. De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes. Science 37, 655–662 (2021).
Janos, D. P. Plant responsiveness to mycorrhizas differs from dependence upon mycorrhizas. Mycorrhiza 17, 75–91 (2007).
Sawers, R. J., Gebreselassie, M. N., Janos, D. P. & Paszkowski, U. Characterizing variation in mycorrhiza effect among diverse plant varieties. Theor. Appl. Genet. 120, 1029–1039 (2010).
Jakobsen, I., Gazey, C. & Abbott, L. K. Phosphate transport by communities of arbuscular mycorrhizal fungi in intact soil cores. New Phytol. 149, 95–103 (2001).
Lekberg, Y., Hammer, E. C. & Olsson, P. A. Plants as resource islands and storage units—Adopting the mycocentric view of arbuscular mycorrhizal networks. FEMS Microbiol. Ecol. 74, 336–345 (2010).
Pepe, A., Giovannetti, M. & Sbrana, C. Appressoria and phosphorus fluxes in mycorrhizal plants: Connections between soil-and plant-based hyphae. Mycorrhiza 30, 589–600 (2020).
van’t Padje, A. et al. Temporal tracking of quantum-dot apatite across in vitro mycorrhizal networks shows how host demand can influence fungal nutrient transfer strategies. ISME J. 15, 435–449 (2021).
Munkvold, L., Kjøller, R., Vestberg, M., Rosendahl, S. & Jakobsen, I. High functional diversity within species of arbuscular mycorrhizal fungi. New Phytol. 164, 357–364 (2004).
Chu, Q. et al. Soil plant-available phosphorus levels and maize genotypes determine the phosphorus acquisition efficiency and contribution of mycorrhizal pathway. Plant Soil 449, 357–371 (2020).
Aguilar-Trigueros, C. A., Boddy, L., Rillig, M. C. & Fricker, M. D. Network traits predict ecological strategies in fungi. ISME Commun. 2, 1–11 (2022).
Pepe, A., Sbrana, C., Ferrol, N. & Giovannetti, M. An in vivo whole-plant experimental system for the analysis of gene expression in extraradical mycorrhizal mycelium. Mycorrhiza 27, 659–668 (2017).
Jones, J. J. B., Benjamin, W. & Mills, H. A. Plant Analysis Handbook: A Practical Sampling, Preparation, Analysis, and Interpretation Guide (Micro-Macro Publishing, 1991).
Ririe, K. M., Rasmussen, R. P. & Wittwer, C. T. Product differentiation by analysis of DNA melting curves during the polymerase chain reaction. Anal. Biochem. 245, 154–160 (1997).
Vangelisti, A. et al. Gene expression in Rhizoglomus irregulare at two different time points of mycorrhiza establishment in Helianthus annuus roots, as revealed by RNA-seq analysis. Mycorrhiza 30, 373–387 (2020).
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25, 402–408 (2001).
Pepe, A., Giovannetti, M. & Sbrana, C. Different levels of hyphal self-incompatibility modulate interconnectedness of mycorrhizal networks in three arbuscular mycorrhizal fungi within the Glomeraceae. Mycorrhiza 26, 325–332 (2016).
Giovannetti, M. & Mosse, B. An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots. New Phytol. 84, 489–500 (1980).
Acknowledgements
This study received financial support from the University of Pisa (Fondi di Ateneo), Università Cattolica del Sacro Cuore, Piacenza, and CNR.
Author information
Authors and Affiliations
Contributions
A.L., A.M., A.T. and M.G. conceived the study. A.T., C.S. and L.G. planned and designed the experiments, A.L., A.T., C.S., L.A. and L.G. performed experiments and analysed data. A.T., C.S. and L.G. wrote the first draft and all authors contributed to the manuscript. This paper was part of L.G.’s doctoral thesis work at the University of Pisa.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Giovannini, L., Sbrana, C., Giovannetti, M. et al. Diverse mycorrhizal maize inbred lines differentially modulate mycelial traits and the expression of plant and fungal phosphate transporters. Sci Rep 12, 21279 (2022). https://doi.org/10.1038/s41598-022-25834-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-25834-7
This article is cited by
-
Bacterial community changes in the presence of AMF in the context of maize with low phosphorus content
Journal of Soils and Sediments (2024)
-
Mycorrhizal status and host genotype interact to shape plant nutrition in field grown maize (Zea mays ssp. mays)
Mycorrhiza (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.