The present study is aimed at disclosing metabolic profile alterations in the leaves of the Medicago lupulina MlS-1 line that result from high-efficiency arbuscular mycorrhiza (AM) symbiosis formed with Rhizophagus irregularis under condition of a low phosphorus level in the substrate. A highly effective AM symbiosis was established in the period from the stooling to the shoot branching initiation stage (the efficiency in stem height exceeded 200%). Mycorrhization led to a more intensive accumulation of phosphates (glycerophosphoglycerol and inorganic phosphate) in M. lupulina leaves. Metabolic spectra were detected with GS-MS analysis. The application of complex mathematical analyses made it possible to identify the clustering of various groups of 320 metabolites and thus demonstrate the central importance of the carbohydrate and carboxylate-amino acid clusters. The results obtained indicate a delay in the metabolic development of mycorrhized plants. Thus, AM not only accelerates the transition between plant developmental stages but delays biochemical “maturation” mainly in the form of a lag of sugar accumulation in comparison with non-mycorrhized plants. Several methods of statistical modeling proved that, at least with respect to determining the metabolic status of host-plant leaves, stages of phenological development have priority over calendar age.

The nature of plant–fungi interaction at early stages of arbuscular mycorrhiza (AM) development is still a puzzling problem. To investigate the processes behind this interaction, we used the Medicago lupulina MlS-1 line that forms high-efficient AM symbiosis with Rhizophagus irregularis. AM fungus actively colonizes the root system of the host plant and contributes to the formation of effective AM as characterized by a high mycorrhizal growth response (MGR) in the host plant. The present study is aimed at distinguishing the alterations in the M. lupulina root metabolic profile as an indicative marker of effective symbiosis. We examined the root metabolome at the 14th and 24th day after sowing and inoculation (DAS) with low substrate phosphorus levels. A GS-MS analysis detected 316 metabolites. Results indicated that profiles of M. lupulina root metabolites differed from those in leaves previously detected. The roots contained fewer sugars and organic acids. Hence, compounds supporting the growth of mycorrhizal fungus (especially amino acids, specific lipids, and carbohydrates) accumulated, and their presence coincided with intensive development of AM structures. Mycorrhization determined the root metabolite profile to a greater extent than host plant development. The obtained data highlight the importance of active plant–fungi metabolic interaction at early stages of host plant development for the determination of symbiotic efficiency.

The study is aimed at revealing the effects of Rhizophagus irregularis inoculation on the transcriptome of Medicago lupulina leaves at the early (second leaf formation) and later (flowering) stages of plant development. A pot experiment was conducted under conditions of low phosphorus (P) level in the substrate. M. lupulina plants were characterized by high mycorrhizal growth response and mycorrhization parameters. Library sequencing was performed on the Illumina HiseqXTen platform. Significant changes in the expression of 4863 (p_adj < 0.01) genes from 34049 functionally annotated genes were shown by Massive Analysis of cDNA Ends (MACE-Seq). GO enrichment analysis using the Kolmogorov–Smirnov test was performed, and 244 functional GO groups were identified, including genes contributing to the development of effective AM symbiosis. The Mercator online tool was used to assign functional classes of differentially expressed genes (DEGs). The early stage was characterized by the presence of six functional classes that included only upregulated GO groups, such as genes of carbohydrate metabolism, cellular respiration, nutrient uptake, photosynthesis, protein biosynthesis, and solute transport. At the later stage (flowering), the number of stimulated GO groups was reduced to photosynthesis and protein biosynthesis. All DEGs of the GO:0016036 group were downregulated because AM plants had higher resistance to phosphate starvation. For the first time, the upregulation of genes encoding thioredoxin in AM plant leaves was shown. It was supposed to reduce ROS level and thus, consequently, enhance the mechanisms of antioxidant protection in M. lupulina plants under conditions of low phosphorus level. Taken together, the obtained results indicate genes that are the most important for the effective symbiosis with M. lupulina and might be engaged in other plant species.