Stoichiometry of carbon, nitrogen and phosphorus is closely linked to trophic modes in orchids
BMC Plant Biology volume 23, Article number: 422 (2023)
Mycorrhiza is a ubiquitous form of symbiosis based on the mutual, beneficial exchange of resources between roots of autotrophic (AT) plants and heterotrophic soil fungi throughout a complex network of fungal mycelium. Mycoheterotrophic (MH) and mixotrophic (MX) plants can parasitise this system, gaining all or some (respectively) required nutrients without known reciprocity to the fungus. We applied, for the first time, an ecological stoichiometry framework to test whether trophic mode of plants influences their elemental carbon (C), nitrogen (N), and phosphorus (P) composition and may provide clues about their biology and evolution within the framework of mycorrhizal network functioning.
We analysed C:N:P stoichiometry of 24 temperate orchid species and P concentration of 135 species from 45 plant families sampled throughout temperate and intertropical zones representing the three trophic modes (AT, MX and MH). Welch’s one-way ANOVA and PERMANOVA were used to compare mean nutrient values and their proportions among trophic modes, phylogeny, and climate zones. Nutrient concentration and stoichiometry significantly differentiate trophic modes in orchids. Mean foliar C:N:P stoichiometry showed a gradual increase of N and P concentration and a decrease of C: nutrients ratio along the trophic gradient AT < MX < MH, with surprisingly high P requirements of MH orchids. Although P concentration in orchids showed the trophy-dependent pattern regardless of climatic zone, P concentration was not a universal indicator of trophic modes, as shown by ericaceous MH and MX plants.
The results imply that there are different evolutionary pathways of adaptation to mycoheterotrophic nutrient acquisition, and that the high nutrient requirements of MH orchids compared to MH plants from other families may represent a higher cost to the fungal partner and consequently lead to the high fungal specificity observed in MH orchids.
Mycorrhiza is a ubiquitous form of symbiosis, since over 90% of plant species are estimated to form mutualistic relationships with soil fungi to access essential nutrients. In exchange for carbon (C) produced by plants in photosynthesis, soil fungi supply water and micro- and macroelements . It is estimated that approximately 80% of nitrogen (N) and phosphorus (P), which are considered most limiting to plant growth [2,3,4,5], is provided by mycorrhizal fungi [6,7,8]. It has been frequently reported that hyphae of mycorrhizal fungi interconnect roots of various plants resulting in a complex web called Common Mycorrhizal Network (CMN), where nutrients and signalling compounds can be exchanged between connected plants [9, 10]. Functional evidence of such transfer is provided by mycoheterotrophic plants . They are devoid of chlorophyll and derive all their nutrients, including organic C, from autotrophic plants via mycelium of shared mycorrhizal fungi and are therefore also called epiparasites [12,13,14]. In addition to MH plants, there are also green, photosynthetically active plants that may derive fungal C with little or no reciprocity to the mycorrhizal partner . In this way they cover deficiencies of C resulting from low photosynthetic efficiency in the low-light habitats they occupy. Because they tap on two C sources, they are called mixotrophic or partially mycoheterotrophic [16,17,18,19,20].
Evaluation and quantification of resource translocation between plants sharing CMN are an ongoing challenge . Although many experiments rely on stable isotopes to detect C and N nutrient flow within CMN, ecological stoichiometry (ES) also provides a useful framework linking differences in the composition of an organism’s multiple chemical elements and ecological interactions in an ecosystem to offer predictions about mycorrhizal functioning [22,23,24,25]. Indeed, the stoichiometry of the exchange of plant C for fungal P and N has been found to play a vital role in the rate and efficiency of resource exchange between mycorrhizal partners shaping not only the strengths of interactions, but also their mutualistic versus parasitic nature [22, 26]. Going further, N:K stoichiometry of soil as well as soil N and P concentrations have been identified as key determinants of the local occurrence and density of MH plants in tropical rainforest [27, 28]. A growing body of evidence also suggest that the rate and direction of resource transfer in CMNs reflect plant nutrient requirements [25, 29, 30]. This, within the ES framework, implies that both biochemical and physiological mechanisms of specialisation of MH plants to their fungal partner(s) could be reflected in the elemental composition of a plant.
Determining the elemental composition of MH plants may therefore help to understand the requirements of parasites of CMN for the main nutrients and ecology of this trophic mode, though such a determination has yet to be made.
All orchid species are MH during germination of their almost reserveless, dust-like seeds and subsequent early ontogenesis [12, 13, 31]. Only a minority of known orchid species stay fully MH in adulthood. Most of them develop green leaves when adult and become photosynthetic and autotrophic (AT) or mixotrophic (MX) if they continue to derive C from associated fungi while remaining photosynthetically active. Sometimes they produce variegated individuals , but rarely do MX orchids form achlorophyllous, non-photosynthetic variants called albinos, which survive fully on MH nutrition [17, 29].
Although most orchids are green, this family encompasses ca. 250 MH species and outnumbers other plant families in which this trophic mode has been identified [13, 33, 34]. In all, orchids are thus a good model to study the elemental composition of the organism in relation to trophic modes.
Interestingly, MH and MX orchids have been found to exhibit unusually high nitrogen concentration compared to autotrophic reference plants from the same habitats [14, 35,36,37]. Given the interdependence of N and P in biochemical functions and cellular components, N and P ratios are tightly linked , and we expect that MX and especially MH orchids are also enriched in P. Many terrestrial orchids, especially MH, have a simplified root system , which indicates a heavy reliance on mycorrhizal pathways for P and other nutrient uptake. Given that P plays a key role in the establishment and sustainability of plant-fungal symbioses [40,41,42] while also being a limiting factor of plant growth globally [5, 43, 44], high P and N requirement revealed by MH orchids may be a crucial factor shaping orchids’ ecological interactions. Functional diversity in plant–fungus combinations may lead to selection for fungi with required nutritional capability, as a driver of evolution of partner specificity and niche diversification [45, 46]. However, in contrast to N, the P concentrations in orchids as well as other plants showing various trophic modes remain largely unexplored [28, 47].
Here we introduce the first data on C:N:P stoichiometry for temperate orchids displaying three trophic modes, namely AT, MX and MH. We also present P concentration analyses extended to other plant groups with different trophic modes in both temperate and intertropical zones. Our questions are as follows:  Does elemental stoichiometry vary among orchids representing different trophic modes ? Do MX and MH stoichiometric C:N:P signatures reflect C shortage ? Is the tissue concentration of P influenced by trophic modes or is it family-dependent ? Can P concentration serve as an additional explanatory variable to characterise trophic modes?
Materials and methods
Sample collection and element measurements
Plant samples were collected and analysed in three sets. The first set (#1) comprised 200 samples collected between 2004 and 2020 from 10 localities covering both temperate and neotropical regions. The second set (#2) encompassed 602 samples collected from 13 populations of terrestrial plants growing in the European temperate zone harvested from 2017 to 2020. Taxonomic identification of the samples was carried out by the specialists listed in Additional file 1; Table S1A. Samples for which it was possible to make a voucher were deposited in Herbarium Universitetas Gedanensis – UGDA (Gdańsk, Poland), their numbers are given in the Additional file 1; Table S2. We used both sample sets to examine P content in a total of 802 samples representing 135 species and 45 families. Among them, orchids constituted 343 biological samples from 39 species, spanning the three trophic modes (AT, MX and MH). The third set (#3) mostly comprised orchid samples and a selection (26 samples) of reference autotrophic, non-orchid species (REF) from set #2. These samples were reanalysed for C and N concentrations to examine the C:N:P stoichiometry of three trophic modes in temperate orchids. Due to insufficient material for this analysis, each sample was a mix of two biological samples for a given species per given population (Additional file 1; Table S1B). In this set of 138 samples, AT orchids constituted 66% of samples, while MX and MH orchids accounted for 25% and 9%, respectively, which roughly reflects the proportion of species numbers in the orchid flora of temperate regions .
At each sampling site, leaves (or stem in the case of MH) were collected from 1 to 5 different individuals representing all trophic modes present locally, with autotrophic species serving as a reference. All plants were collected at the flowering stage, since P and N concentrations vary during a plant’s lifetime . They were immediately placed in silica gel and, after transfer to the laboratory, oven-dried at 70 °C until constant weight and ground to a fine powder using a Quiagen TissueLyzer II.
Samples from set #1 were digested with nitric acid (HNO3, 65%), hydrogen peroxide (H2O2, 30%), and subsequently phosphate was analysed with a NIVO microplate reader (Perkin-Elmer France SAS, Villebon sur Yvette, France). The samples collected in Brazil and Japan were analysed in the country of origin. Perchloric acid was used for digestion of samples from set #2 and phosphate was analysed colorimetrically using flow-injection analysis (FIA Lachat QC8500, Lachat Instruments, USA). The P concentration data (set #1 and #2) were analysed together because the impact of the method on the variance was one hundred times smaller than the differences between the study groups. C and N percentages of plant tissue from set #3 were measured on a NC Elemental analyser (ThermoQuest, Germany). P concentration values for a given sample mix were obtained by averaging the results from the relevant samples from #2 P analysis. Units for leaf C, N and P were expressed as weight per gram of dry mass (g DM) and elemental ratio as molar units.
Species were assigned to AT, MX or MH trophic modes based on literature data on 13 C isotopic signature (enrichment in 13 C; Additional file 1; Table S3). There are some claims, based on 18O and 2 H isotopic signature , that some AT orchids with 13 C isotopic signature typical for autotrophs could be MX. Yet, these statements rely on evidence for a nutrient flow from the fungus to the orchid without data on the reverse flow: thus, there is no data on the direction of the net flow and thus no direct evidence for MX nutrition. We therefore restricted the MX category to the orchids and plants displaying characteristic 13 C isotopic signature, for which ability to extract fungal carbon is further supported by the existence of achlorophyllous variants . Moreover, green plant species for which no stable isotope data were found in the literature were assumed AT, likely by large the most frequent status in orchids.
For P concentration analysis (sets #1 and #2 pooled), samples were assigned to the following categories: trophy (AT, MX, MH), phylogeny (Orchidaceae, Ericaceae, B&G - Burmaniaceae and Gentianaceae) and climate (TROP - intertropical zone and TEMP - temperate latitudes; see, raw data in Additional file 3). The “climate” category acknowledges the well-known latitudinal change of N and P tissue concentration in plants . In this study, the relation between P concentrations and latitudinal distribution of plant-sampling sites were described and tested by linear regression and Pearson correlation. Welch’s one-way ANOVA was used to compare mean nutrient values and their proportions among the categories studied due to unequal variances between groups, followed by the Games-Howell post-hoc test. All statistical analyses were performed in R  and based on logarithmically transformed, x’= log10(x) P data.
For C:N:P stoichiometry (set #3) samples were analysed as one factor: trophy with three categories represented by orchids (AT, MX and MH); autotrophic plants from various families were a reference (REF). To illustrate the level of multivariate similarity between the studied categories, non-metric multidimensional scaling (n-MDS) was performed based on the matrix of Euclidean distances between the samples. Prior to n-MDS, all values of C:N:P were standardised by variable and logarithmically transformed, x’= log(x + 1). Based on the same distance matrix, the differences between groups were tested by one-way PERMANOVA with post-hoc tests (999 permutations).
The growth rate hypothesis (GRH) posits that organisms with higher growth rates require disproportionately higher P than N concentrations, resulting in a scaling exponent of N to P concentrations below unity (i.e. α < 1.0) . We conducted reduced major axis (SMA) regression to study the pattern of the leaf N versus P scaling relationship in relation to trophic modes AT and MH in orchids and reference AT plants (REF) for comparison. To perform the analyses, N and P concentration data from set #3 were logarithmically transformed, x’= log10(x). Analyses were performed in R package SMAtr 3 .
Patterns of leaf C, N, P, and their stoichiometric ratios according to trophic modes
Based on the results of the stoichiometric data pool, the average C:N:P molar ratio for AT, MX and MH orchids was 552:26:1, 503:31:1, and 220:14:1, respectively (see the raw data in Additional file 3). The trophy factor explained 44.8% of variation in the dataset (Fig. 1, Additional file 2; Table S4) and significantly (p = 0.001; Table 1) differentiated trophic groups based on their mean nutrient stoichiometric proportion: C:N (AT = 22.1, MX = 16.5, MH = 12.9 ), C:P (AT = 558, MX = 422, MH = 174), N:P (AT = 26.8, MX = 25.5, MH = 13.4) (Fig. 2, Additional file 2; Table S5). MH orchids significantly differentiated from AT, MX and REF with the lowest ratios of C:N and C:P and N:P (Fig. 2, Additional file 2; Table S5). The values of nutrient concentration and their stoichiometric ratios for MX trophic mode formed a continuum between AT and MH, except for N:P, which was similar in all green plants (AT, MX orchids and REF). In orchids, the highest variation coefficient (calculated as standard deviation / mean) of most parameters in MX further confirmed their intermediate position (Figs. 1 and 2, Additional file 2; Table S5).
The mean elemental concentrations of orchid biomass significantly differed between trophic modes: N (AT = 23.7, MX = 31.6; MH = 39.4 mg N per g DM) and P (AT = 2.01, MX = 3.23; MH = 6.70 mg P per g DM), except for C (Fig. 2, Additional file 2; Table S5, S6). Hence, mean tissue concentration values for N and P increased markedly along the trophic gradient AT < MX < MH in orchids. The tendency also applied to a narrowed dataset including only the tribe Neottieae (Additional file 2; Table S5) and was also visible in the differences in average value of nutrient concentration in green MX (N = 27.4 mg/g DM, P = 3.20 mg/g DM) versus albino MH (N = 48.9 mg/g DM, P = 7.63 mg/g DM) variants of Cephalanthera damasonium (Mill.) Druce. It is noteworthy that temperate AT orchids have slightly but significantly higher mean N concentration (N = 23.6 mg per g DM; p = 0.019) compared to reference autotrophic non-orchid species (REF) (N = 19.9 mg/g DM). Thus, we may claim that terrestrial, temperate AT orchids tend to accumulate more N than AT non-orchids.
The scaling exponent of N to P concentrations dropped gradually from AT reference plants (α = 0.814) through AT orchids (α = 0.761) to reach the lowest value (α = 0.668) in MH orchids (Additional file 2; Fig. S1, Table S9). Although slopes of regression lines did not differ significantly between groups (likelihood ratio = 0.335; df = 2; p = 0.845), regression line shifts were significant (Wald = 116.1; df = 2; p < 0.001).
Phosphorus concentration differences through trophy and phylogeny
Analyses of pooled data from sets #1 and #2 demonstrated that P concentration in orchids increases along the trophic gradient AT < MX < MH and this trend stays salient irrespective of geographical range, as shown by significant differences between trophic modes in terrestrial orchids both in intertropical and temperate zones (Fig. 3, Additional file 2; Fig S2, Table S7). A significant relation between P concentrations and latitudinal distribution of plant-sampling sites (p < 0.001; Additional file 2; Fig. S2) revealed that P concentration increased with latitude. This tendency was probably the reason why the concentration of P in tropical MH orchids was significantly lower than in temperate ones (p < 0.001; Fig. 3, Additional file 2; Table S7).
AT terrestrial orchids from both temperate and tropical zones showed slightly but significantly higher P concentrations (p = 0.001) compared to reference AT non-orchid species (REF, Fig. 3, Additional file 2; Fig. S2, Table S7). A similar trend was not detected in set #3 (p = 0.28; Additional file 2; Table S6), which may be the result of a smaller data set than sets #1 and #2 pooled (Fig. 3, Additional file 2; Table S7). Thus, we may claim that terrestrial AT orchids regardless of climatic zone tended to accumulate more P than non-orchid autotrophic plants. Significantly higher P content compared to the AT reference (p < 0.001) is also a feature of tropical, non-orchid MH represented by Gentianaceae and Burmaniaceae (Fig. 3; B&G), which, unlike orchids, are associated with arbuscular mycorrhizal fungi. In contrast, temperate Ericaceae, which like orchids are associated with asco- and basidiomycetes, showed no differences, either in comparison with REF plants (p > 0.999) or between MX and MH within the family (p = 0.999) (Fig. 3, Additional file 2; Table S7).
Nutrient concentration and stoichiometry vary between trophic modes in orchids
In this study, we provide for the first time the C:N:P stoichiometry of temperate orchids representing different trophic modes (AT, MX and MH). Various studies have confirmed that an increasing reliance on MH nutrition entails an increase in isotopic 15 N enrichment  and leaf N concentration in orchids [17, 29, 35, 36, 54]. Not only did we confirm the trend of growing N concentration along the trophic gradient (AT < MX < MH), but we also found the same trend for P. The analyses indicated that the mean N (39.4 mg/g DM) and P (6.70 mg/g DM) concentrations in leaves of MH orchids were nearly two (for N) and three times (for P) higher than in AT orchids and REF. Comparison of the respective values with the global average for temperate zone herbs (N = 22.4 mg g− 1 and P = 1.74 mg g− 1; ) also confirms the tendency of temperate MH orchids to have high nutrient accumulation. Remarkably, even AT orchids showed higher N and P concentrations than AT reference plants. One could wonder whether this may result from an admixture of truly AT orchids with hidden MX species. This is especially relevant, since some orchids associated with ‘rhizoctonias’ and usually considered AT are now claimed to receive organic matter from the ‘rhizoctonias’ fungi, based on their 18O and 2H isotopic signature : yet the reverse flow from orchid to fungus remains unevaluated, so that the MX status of such orchids remains speculative. To exclude the possibility of hidden MX species into account, we analysed the frequency distribution of N (and P) concentration values in the AT orchid group: bimodality is expected if this group is not homogenous, yet a single mode was observed (Additional file 2; Fig. S3). In addition, the low (lowest of all tested groups) standard deviation of N and P concentrations among AT orchids further supports the claim of homogeneity (Additional file 2; Table S5). We rather believe that the extreme nutrient contents reflect a physiological particularity of this family, very possibly connected with their mycorrhizal fungi - the ‘rhizoctonias’, which are not known to be mycorrhizal in other plants . The relatively large genome size typical of Orchidaceae  may also contribute to the higher average N and P levels in AT orchids compared to AT reference plants, as nucleic acids are rich in these elements. Although this hypothesis needs to be thoroughly tested, it is unlikely to affect calculations between trophic groups within the Orchidaceae as the average genome size is subfamily-specific [56,57,58] rather than trophic mode-specific (the largest genome in this study is represented by AT orchid - Cypripedium calceolus L.; see Additional file 2; Table S8).
The MX group of species formed an almost continuous gradient of nutrient concentration between AT and MH orchids, confirming its flexible mechanism of nutrient acquisition and use [18, 19, 59,60,61]. Surprisingly, Hynson et al.  showed no significant differences between MX and MH with respect to N concentration: we believe that this was due to the combination of temperate MH orchid species and the (sub)tropical ones, since the latter may have lowered the mean N concentration value due to the common feature of decreasing plant N concentrations towards the equator .
The phylogenetic signal connected with genetically and physiologically determined growth form, size and structure may have a major impact on plant stoichiometry [4, 62, 63]. However, we planned our studies to minimise bias by sampling within closely related species, applying a uniform sampling protocol and collecting plants at the same phenological stage to control for seasonal variation in foliar nutrient concentration . We therefore expect the influence of the phylogenetic signal on the data to be negligible. The fact that all MX and MH samples (and 16% of AT samples) belong to the subfamily Epidendroideae indicates that, at least at the subfamily level, trophic modes, rather than phylogeny, predicted stoichiometry and N and P concentrations (Additional file 2; Table S5). This, in addition to the significant differences in P concentration between AT and MH orchids in the tropics, suggests a general trend in orchids, irrespective of climate zone and phylogenetic traits.
A possible explanation of the high N concentration in MX and MH orchids is acquisition of C in the form of amino acids [54, 65,66,67]. However, the reason for high P levels in orchids is unclear. At least the inorganic form of P can be transferred by fungi to AT orchids [47, 68]. Lysed fungal hyphae, which are on average richer in N and P than green plant tissue [69, 70], may also be the source of these nutrients [71, 72]. Indeed, differences in N concentration of fungal hyphae were found to influence N gain of MX Epipactis species . However, digestion of the fungal matter inside root cells is commonly observed in many non-orchid MH  and AT plants , leaving an open question about the significance of this process in gaining exceptionally high P concentrations in orchids. However, it is a fact that the value of P concentration per unit of N concentration increases faster in MH orchids than in AT plants, including orchids, indicating differences in nutrient acquisition and/or recycling.
It is presumed that in the face of C deficiency, MH orchids reprogramme their metabolism to use N for energy [54, 75]. This raises the question of whether high P concentration shapes metabolism in a different way compared to AT plants or simply enlarges the P storage pool. Studies determining the contribution of the storage (inorganic P), structural (e.g., phospholipids) and biochemically functional (RNA, ATP) forms of P in MH orchids would clarify the functional significance of their high P concentrations. However, the lowest value of both scaling exponents of N to P and N:P coefficient of variation observed in MH orchids may point to high nutrient demand connected with high growth rate , which was also suggested based on transcriptome analyses .
Higher concentrations of P and N in orchids than in non-orchid MH plants  implying high nutrient demand for fungal community and CMN also raise the question of their implications for the evolution of orchid fungal specificity, demography and distribution. More studies showing the patterns of MH orchid demography in relation to the abundance and distribution of mycorrhizal fungi, taking into account environmental variables, as in , but also considering the stoichiometry of the soil and organisms (plants and fungi) supporting MH growth, could help to better understand the functioning MH orchids within CMN.
Mycoheterotrophic lifestyle affects orchid stoichiometry
The low stoichiometric ratios (C:N, C:P and N:P) characteristic of MH orchids may result from passive reflection of those of their mycorrhizal fungi, as has been observed in fungi parasitising algae . However, the higher mean (molar) C:N (14.4) and C:P (348) values in ectomycorrhizal fungi , from which partners of temperate MH orchids are recruited, may suggest the presence of a stoichiometric imbalance in the transfer pathway. To compensate for imbalance, consumers recycle C and nutrients differently using multiple mechanisms, including excretion of elements that are in excess . In MH orchids, the demand for C met by uptake of N-rich sources of C  may result in an excess of N. Ammonium ions released at the orchid-fungus symbiotic interface as described in [66, 79] may simply be the result of excretion of excess N. Whether this is a sign of mutualism, as suggested by  or a way to mask “cheating” is debatable. Nevertheless, the actual existence of a stoichiometric imbalance in relation to C:N:P would have to be checked directly on the orchid-fungus pair.
While the concentration of C does not differ significantly between AT and MH orchids, various MH features suggest parsimonious C management: the typically small size of MH plants, their leaflessness and the short life-span of the stem , which limit the need for C-costly structural investment to build mechanically resistant cell walls . Carbon limitation may be the main reason for extremely slow growth (up to several years) in the early ontogenesis phase when all orchids show MH nutrition . Therefore, the nutritional composition of the adult MH plant, i.e. high P concentration and low (< 10) N:P ratio, as well as the low scaling component that usually is the hallmark of a fast-growing organism [23, 52, 81] may seem paradoxical. At least two possibilities can be considered here:  C, rather than N and/or P, is the main growth limiting factor ;  C is not the limiting factor during the flowering stage, e.g. due to the release of C stored in the rhizome or roots, and the flowering shoot of MH plants is indeed a fast-growing organ. In Epipogium aphllum Sw. the flowering shoot grows 4 cm per day and the plant completes its life cycle in a fortnight . Presumed fast growth rate during flowering shoot expansion may imply a high nutrient demand, which may make the MH plant a strong physiological sink in mycorrhizal networks [83,84,85] further leading to high nutrient concentration. Although other physiological mechanisms associated with nutrient mass flow, such as high respiration demand of heterotrophic stem  or high transpiration rate  may contribute to the effect. More analyses of the vegetative parts of MH plants, which are slowly growing and may have less pronounced N and P demands, may test the second hypothesis. Either way, the requirements of the plant rather than fungal species have a stronger influence on resource flow, as exemplified by the different N concentrations of albino (MH) and green (MX) variants of C. damasonium sharing mycorrhizal fungi . Higher N and P concentrations in the albino than in the green variant of C. damasonium found in our data also confirm this pattern.
Phosphorus level is driven by trophic mode and phylogeny
Extended P concentration data across phylogeny, trophic modes, and climatic zones from sets #1 and #2 confirmed that P concentration in leaves of terrestrial orchids increases along the AT < MX < MH trophic gradient. This trend applies to all climatic zones, as despite generally lower P concentrations in tropical MH orchids than temperate ones, they show the same 3-fold higher P concentration relative to non-orchid reference AT plants. The lower P concentration in tropical MH orchids compared to temperate ones may simply be due to differences in P availability, which is lower in old, rain-washed tropical soils [5, 86]. This is supported by the same decreasing trend of N:P ratio with latitude observed in fungi .
However, P concentration is not a universal indicator of MH plant dependency on fungal resources, as exemplified by the lack of differences between ericaceous MH, MX and reference AT plants as well as high differences between Burmaniaceae, Gentianaceae and Orchidaceae. These findings may indicate different ecophysiological adaptation for nutrient acquisition from the fungi, developed during evolution on the path to mycoheterotrophy. Interestingly, it has already been shown that orchid and ericaceous MH species exhibit different enrichments in the stable isotopes of C and N as well as N concentration .
Trophic modes in orchids i.e. autotrophy (AT) mixotrophy (MX) and mycoheterotrophy (MH), shape both the demand for C and nutrients (N and P) and their stoichiometry. A markedly increasing demand for N along the trophic gradient in orchids agrees with the data on the acquisition of N-rich C substrates. High P concentration relative to AT plants is an inherent feature of orchids regardless of climatic zone but does not characterise MH species from other families, indicating different, family-dependent factors shaping adaptation to mycoheterotrophic nutrient acquisition. Future studies should address the relationship between MH plant partner specificity and plant nutrient demand and the stoichiometric balance of resource flow between fungi and MH plants.
The datasets supporting the conclusions of this article are included within the article and its additional files.
Smith SE, Read D. In: Smith SE, Read D, editors. Mycorrhizal symbiosis [Internet]. Third edit. London: Academic Press; 2008. pp. 1–9.
Ågren GI, Weih M. Multi-dimensional plant element Stoichiometry—Looking beyond Carbon, Nitrogen, and Phosphorus. Front Plant Sci. 2020;11:23.
Elser JJ, Bracken MES, Cleland EE, Gruner DS, Harpole WS, Hillebrand H, et al. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems. Ecol Lett. 2007;10(12):1135–42.
Reich PB, Oleksyn J. Global patterns of plant leaf N and P in relation to temperature and latitude. Proc Natl Acad Sci U S A. 2004;101(30):11001–6.
Wright SJ. Low phosphorus levels limit carbon capture by amazonian forests. Nature Publishing Group; 2022.
Behie SW, Bidochka MJ. Nutrient transfer in plant-fungal symbioses. Trends Plant Sci [Internet]. 2014;19(11):734–40.
George E, Marschner H, Jakobsen I. Role of arbuscular mycorrhizal fungi in uptake of phosphorus and nitrogen from soil. Crit Rev Biotechnol. 1995;15(3–4):257–70.
van der Heijden MGA, Martin FM, Selosse MA, Sanders IR. Mycorrhizal ecology and evolution: the past, the present, and the future. New Phytol. 2015;205(4):1406–23.
Gorzelak MA, Asay AK, Pickles BJ, Simard SW. Inter-plant communication through mycorrhizal networks mediates complex adaptive behaviour in plant communities. AoB Plants [Internet]. 2015;7.
Cahanovitc R, Livne-Luzon S, Angel R, Klein T. Ectomycorrhizal fungi mediate belowground carbon transfer between pines and oaks. ISME J. 2022;16(5):1420–9.
McKendrick SL, Leake JR, Taylor DL, Read DJ. Symbiotic germination and development of myco-heterotrophic plants in nature: Ontogeny of Corallorhiza trifida and characterization of its mycorrhizal fungi. New Phytol. 2000;145(3):523–37.
Leake JR. The biology of myco-heterotrophic (‘saprophytic’) plants. New Phytol. 1994;127(2):171–216.
Merckx VSFT, Mycoheterotrophy. The biology of plants living on fungi. Mycoheterotrophy Biol Plants Living Fungi. 2013;1–356.
Gebauer G, Meyer M. 15 N and 13 C natural abundance of autotrophic and myco-heterotrophic orchids provides insight into nitrogen and carbon gain from fungal association. New Phytol. 2003;160(1):209–23.
Lallemand F, Figura T, Damesin C, Fresneau C, Griveau C, Fontaine N, et al. Mixotrophic orchids do not use photosynthates for perennial underground organs. New Phytol. 2019;221(1):12–7.
Hynson N, Madsen T, Selosse M, Adam I, Ogura-Tsujita Y, Roy M et al. The physiological ecology of mycoheterotrophy. In: Mycoheterotrophy: The Biology of Plants Living on Fungi. 2013. p. 297–342.
Julou T, Burghardt B, Gebauer G, Berveiller D, Damesin C, Selosse MA. Mixotrophy in orchids: insights from a comparative study of green individuals and nonphotosynthetic individuals of Cephalanthera damasonium. New Phytol. 2005;166(2):639–53.
Preiss K, Adam IKU, Gebauer G. Irradiance governs exploitation of fungi: fine-tuning of carbon gain by two partially myco-heterotrophic orchids. Proc R Soc B Biol Sci. 2010;277(1686):1333–6.
Selosse M-A, Roy M. Green plants that feed on fungi: facts and questions about mixotrophy. Trends Plant Sci. 2009;14(2):64–70.
Stöckel M, Těšitelová T, Jersáková J, Bidartondo MI, Gebauer G. Carbon and nitrogen gain during the growth of orchid seedlings in nature. New Phytol. 2014;202(2):606–15.
Figueiredo AF, Boy J, Guggenberger G. Common Mycorrhizae Network: a review of the Theories and Mechanisms behind Underground interactions. Front Fungal Biol. 2021;2:1–13.
Johnson NC. Resource Stoichiometry elucidates the structure and function of Arbuscular Mycorrhizas across Scales. Source New Phytol. 2010;185(3):631–47.
Sterner R, Elser JJ. Ecological stoichiometry: the Biology of Elements from Molecules to the Biosphere. In 2002. p. 439.
Van de Waal DB, Elser JJ, Martiny AC, Sterner RW, Cotner JB. Editorial: progress in ecological stoichiometry. Front Microbiol. 2018;9:1–5.
Yang X, Ma Y, Zhang J, Bai H, Shen Y. How arbuscular mycorrhizal fungi drives herbaceous plants’ C: N: P stoichiometry? A meta-analysis. Sci Total Environ. 2023;862:160807.
Schwartz MW, Hoeksema JD. Specialization and resource trade: biological markets as a model of mutualisms. Ecology. 1998;79(3):1029–38.
Gomes SIF, van Bodegom PM, Merckx VSFT, Soudzilovskaia NA. Environmental drivers for cheaters of arbuscular mycorrhizal symbiosis in tropical rainforests. New Phytol. 2019;223(3):1575–83.
Sheldrake M, Rosenstock NP, Revillini D, Olsson PA, Joseph Wright S, Turner BL. A phosphorus threshold for mycoheterotrophic plants in tropical forests. Proc R Soc B Biol Sci. 2017;284(1848).
Roy M, Gonneau C, Rocheteau A, Berveiller D, Thomas JC, Damesin C, et al. Why do mixotrophic plants stay green? A comparison between green and achlorophyllous orchid individuals in situ. Ecol Monogr. 2013;83(1):95–117.
Dawson HR, Shek KL, Maxwell TM, Reed PB, Bomfim B, Bridgham SD et al. Agnostic fungi: plant functional types and tissue stoichiometry explain nutrient transfer in common arbuscular mycorrhizal networks of temperate grasslands. bioRxiv [Internet]. 2023; Available from: https://www.biorxiv.org/content/early/2023/06/01/2022.10.05.511035.
Rasmussen HN. Terrestrial orchids: from seed to mycotrophic plant. Cambridge, UK: Cambridge University Press; 1995. p. 433.
Stöckel M, Meyer C, Gebauer G. The degree of mycoheterotrophic carbon gain in green, variegated and vegetative albino individuals of Cephalanthera damasonium is related to leaf chlorophyll concentrations. New Phytol. 2011;189(3):790–6.
Jacquemyn H, Merckx VSFT. Mycorrhizal symbioses and the evolution of trophic modes in plants. J Ecol. 2019;107(4):1567–81.
Perez-Lamarque B, Selosse MA, Öpik M, Morlon H, Martos F. Cheating in arbuscular mycorrhizal mutualism: a network and phylogenetic analysis of mycoheterotrophy. New Phytol. 2020;226(6):1822–35.
Abadie JC, Püttsepp Ü, Gebauer G, Faccio A, Bonfante P, Selosse MA. Cephalanthera longifolia (Neottieae, Orchidaceae) is mixotrophic: a comparative study between green and nonphotosynthetic individuals. Can J Bot. 2006;84(9):1462–77.
Gonneau C, Jersáková J, de Tredern E, Till-Bottraud I, Saarinen K, Sauve M, et al. Photosynthesis in perennial mixotrophic Epipactis spp. (Orchidaceae) contributes more to shoot and fruit biomass than to hypogeous survival. J Ecol. 2014;102(5):1183–94.
Hynson NA, Schiebold JMI, Gebauer G. Plant family identity distinguishes patterns of carbon and nitrogen stable isotope abundance and nitrogen concentration in mycoheterotrophic plants associated with ectomycorrhizal fungi. Ann Bot. 2016;118(3):467–79.
Reich PB, Oleksyn J, Wright IJ, Niklas KJ, Hedin L, Elser JJ. Evidence of a general 2/3-power law of scaling leaf nitrogen to phosphorus among major plant groups and biomes. Proc R Soc B Biol Sci [Internet]. 2010;277(1683):877–83.
Imhof S, Massicotte H, Melville L, Peterson R. Subterranean Morphology and Mycorrhizal Structures. In: Mycoheterotrophy: The Biology of Plants Living on Fungi. 2013. p. 157–214.
Balzergue C, Chabaud M, Barker D, Bécard G, Rochange S. High phosphate reduces host ability to develop arbuscular mycorrhizal symbiosis without affecting root calcium spiking responses to the fungus. Front Plant Sci. 2013;4.426.
Bhalla K, Qu X, Kretschmer M, Kronstad JW. The phosphate language of fungi. Trends Microbiol [Internet]. 2022;30(4):338–49.
Das D, Paries M, Hobecker K, Gigl M, Dawid C, Lam HM et al. PHOSPHATE STARVATION RESPONSE transcription factors enable arbuscular mycorrhiza symbiosis. Nat Commun. 2022;13(1).
Allen K, Fisher JB, Phillips RP, Powers JS, Brzostek ER. Modeling the Carbon cost of Plant Nitrogen and Phosphorus Uptake Across Temperate and Tropical forests the fixation and uptake of. Nitrogen Model. 2020;3(May):1–12.
Hou E, Luo Y, Kuang Y, Chen C, Lu X, Jiang L, et al. Global meta-analysis shows pervasive phosphorus limitation of aboveground plant production in natural terrestrial ecosystems. Nat Commun [Internet]. 2020;11(1):1–9.
Nurfadilah S, Swarts ND, Dixon KW, Lambers H, Merritt DJ. Variation in nutrient-acquisition patterns by mycorrhizal fungi of rare and common orchids explains diversification in a global biodiversity hotspot. Ann Bot. 2013;111(6):1233–41.
Figura T, Tylová E, Jersáková J, Vohník M, Ponert J. Fungal symbionts may modulate nitrate inhibitory effect on orchid seed germination. Mycorrhiza. 2021;31(2):231–41.
Davis B, Lim WH, Lambers H, Dixon KW, Read DJ. Inorganic phosphorus nutrition in green-leaved terrestrial orchid seedlings. Ann Bot. 2022;129(6):669–78.
Kühn R, Pedersen H, Cribb P. Field guide to the orchids of Europe and the Mediterranean. Kew Publishing; 2019.
Zhang H, Wu H, Yu Q, Wang Z, Wei C, Long M, et al. Sampling date, Leaf Age and Root size: implications for the study of Plant C:N:P Stoichiometry. PLoS One [Internet]. 2013;8(4):1–8.
Gebauer G, Preiss K, Gebauer AC. Partial mycoheterotrophy is more widespread among orchids than previously assumed. New Phytol. 2016;211(1):11–5.
R Core Team. R: a Language and Environment for Statistical Computing. Vienna, Austria; 2013.
Tian D, Yan Z, Niklas KJ, Han W, Kattge J, Reich PB, et al. Global leaf nitrogen and phosphorus stoichiometry and their scaling exponent. Natl Sci Rev. 2018;5(5):728–39.
Warton DI, Duursma RA, Falster DS, Taskinen S. Smatr 3– an R package for estimation and inference about allometric lines. Methods Ecol Evol. 2012;3(2):257–9.
Lallemand F, Martin-Magniette ML, Gilard F, Gakière B, Launay-Avon A, Delannoy É, et al. In situ transcriptomic and metabolomic study of the loss of photosynthesis in the leaves of mixotrophic plants exploiting fungi. Plant J. 2019;98(5):826–41.
Selosse MA, Petrolli R, Mujica MI, Laurent L, Perez-Lamarque B, Figura T, et al. The Waiting Room Hypothesis revisited by orchids: were orchid mycorrhizal fungi recruited among root endophytes? Ann Bot. 2022;129(3):259–70.
Leitch IJ, Kahandawala I, Suda J, Hanson L, Ingrouille MJ, Chase MW, et al. Genome size diversity in orchids: consequences and evolution. Ann Bot. 2009;104(3):469–81.
Jersáková J, Trávníček P, Kubátová B, Krejčíková J, Urfus T, Liu Z-J, et al. Genome size variation in Orchidaceae subfamily Apostasioideae: filling the phylogenetic gap. Bot J Linn Soc. 2013;172(1):95–105.
Chumová Z, Záveská E, Hloušková P, Ponert J, Schmidt P-A, Čertner M, et al. Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids. Plant J. 2021;107(2):511–24.
Jacquemyn H, Brys R, Waud M, Evans A, Figura T, Selosse MA. Mycorrhizal Communities and Isotope Signatures in two partially mycoheterotrophic Orchids. Front Plant Sci. 2021;12:1–9.
May M, Jąkalski M, Novotná A, Dietel J, Ayasse M, Lallemand F, et al. Three-year pot culture of Epipactis helleborine reveals autotrophic survival, without mycorrhizal networks, in a mixotrophic species. Mycorrhiza. 2020;30(1):51–61.
Rasmussen HN, Rasmussen FN. Orchid mycorrhiza: implications of a mycophagous life style. Oikos. 2009;118(3):334–45.
Elser JJ, Fagan WF, Kerkhoff AJ, Swenson NG, Enquist BJ. Biological stoichiometry of plant production: metabolism, scaling and ecological response to global change. New Phytol. 2010;186(3):593–608.
He J-S, Wang X, Schmid B, Flynn D, Li X, Reich P, et al. Taxonomic identity, phylogeny, climate and soil fertility as drivers of leaf traits across chinese grassland biomes. J Plant Res. 2010;123:551–61.
Milla R, Palacio S, Maestro-Martínez M, Montserrat-Martí G. Leaf exchange in a Mediterranean shrub: water, nutrient, non-structural carbohydrate and osmolyte dynamics. Tree Physiol. 2007;27(7):951–60.
Dearnaley JDW, Cameron DD. Nitrogen transport in the orchid mycorrhizal symbiosis – further evidence for a mutualistic association. Volume 213. New Phytologist: Blackwell Publishing Ltd; 2017. pp. 10–2.
Fochi V, Chitarra W, Kohler A, Voyron S, Singan VR, Lindquist EA, et al. Fungal and plant gene expression in the Tulasnella calospora–serapias vomeracea symbiosis provides clues about nitrogen pathways in orchid mycorrhizas. New Phytol. 2017;213(1):365–79.
Suetsugu K, Yamato M, Miura C, Yamaguchi K, Takahashi K, Ida Y, et al. Comparison of green and albino individuals of the partially mycoheterotrophic orchid Epipactis helleborine on molecular identities of mycorrhizal fungi, nutritional modes and gene expression in mycorrhizal roots. Mol Ecol. 2017;26(6):1652–69.
Cameron DD, Johnson I, Leake JR, Read DJ. Mycorrhizal acquisition of inorganic phosphorus by the green-leaved terrestrial orchid Goodyera repens. Ann Bot. 2007;99(5):831–4.
Zhang J, Elser JJ, Carbon. Nitrogen: Phosphorus stoichiometry in fungi: a meta-analysis. Front Microbiol. 2017;8:1–9.
Elser JJ, Fagan WF, Denno RF, Dobberfuhl DR, Folarin A, Huberty A, et al. Nutritional constraints in terrestrial and freshwater food webs. Nature. 2000;408(6812):578–80.
Bougoure J, Ludwig M, Brundrett M, Cliff J, Clode P, Kilburn M, et al. High-resolution secondary ion mass spectrometry analysis of carbon dynamics in mycorrhizas formed by an obligately myco-heterotrophic orchid. Plant Cell Environ. 2014;37(5):1223–30.
Kuga Y, Sakamoto N, Yurimoto H. Stable isotope cellular imaging reveals that both live and degenerating fungal pelotons transfer carbon and nitrogen to orchid protocorms. New Phytol. 2014;202(2):594–605.
Schiebold JMI, Bidartondo MI, Karasch P, Gravendeel B, Gebauer G. You are what you get from your fungi: Nitrogen stable isotope patterns in Epipactis species. Ann Bot. 2017;119(7):1085–95.
Pumplin N, Harrison MJ. Live-cell imaging reveals Periarbuscular membrane domains and Organelle Location in Medicago truncatula roots during Arbuscular Mycorrhizal Symbiosis. Plant Physiol. 2009;151(2):809–19.
Jąkalski M, Minasiewicz J, Caius J, May M, Selosse M-A, Delannoy E. The genomic impact of Mycoheterotrophy in Orchids. Front Plant Sci. 2021;12:1–16.
McCormick MK, Whigham DF, O’Neill JP, Becker JJ, Werner S, Rasmussen HN, et al. Abundance and distribution of Corallorhiza odontorhiza reflect variations in climate and ectomycorrhizae. Ecol Monogr [Internet]. 2009;79(4):619–35.
Frenken T, Wierenga J, Gsell AS, van Donk E, Rohrlack T, Van de Waal DB. Changes in N:P supply ratios affect the ecological stoichiometry of a toxic cyanobacterium and its fungal parasite. Front Microbiol. 2017;8:1–11.
Zechmeister-Boltenstern S, Keiblinger KM, Mooshammer M, Peñuelas J, Richter A, Sardans J, et al. The application of ecological stoichiometry to plant–microbial–soil organic matter transformations. Ecol Monogr. 2015;85(2):133–55.
Jiang Y, Hu X, Yuan Y, Guo X, Chase MW, Ge S, et al. The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions. BMC Plant Biol. 2022;22(1):1–14.
Xing K, Zhao M, Niinemets Ü, Niu S, Tian J, Jiang Y. Relationships between Leaf Carbon and Macronutrients across Woody Species and Forest Ecosystems Highlight how Carbon is allocated to Leaf Structural function. 2021;12:1–13.
Rees TAV, Raven JA. The maximum growth rate hypothesis is correct for eukaryotic photosynthetic organisms, but not cyanobacteria. New Phytol [Internet]. 2021;230(2):601–11.
Krawczyk E. Wybrane aspekty biologii storzana bezlistnego (Epipogium aphyllum, Orchidaceae) w Polsce „Some aspects of ghost orchid’s biology (Epipogium aphyllum, Orchidaceae) in Poland. University of Gdańsk; 2016.
Selosse MA, Richard F, He X, Simard SW. Mycorrhizal networks: des liaisons dangereuses? Trends Ecol Evol. 2006;21(11):621–8.
Simard SW, Beiler KJ, Bingham MA, Deslippe JR, Philip LJ, Teste FP. Mycorrhizal networks: mechanisms, ecology and modelling. Fungal Biol Rev. 2012;26(1):39–60.
Van Der Heijden MGA, Horton TR. Socialism in soil? The importance of mycorrhizal fungal networks for facilitation in natural ecosystems. J Ecol. 2009;97(6):1139–50.
Vitousek PM, Farrington H. Nutrient limitation and Soil Development: experimental test of a Biogeochemical Theory. Biogeohem. 1997;37:63–75.
We want to thank Raphael Buatois, Joana Pimentel, Kateřina Kučerová and Ondřej Žampach for assistance in the laboratory work, LabEx CeMEB, an ANR “Investissements d’avenir” program (ANR-10-LABX-04- 01) for providing platform for phosphorus analysis in France. Yosuke Matsuda for help in organising plant collection, species identification and phosphorus analysis in Japan. Florent Martos and Allain Rousteau for plant species identification from Reunion and Guadelupe (respectively), Michał Filipiak for discussion and comments on the manuscript and David March for English corrections.
The work was supported by budget subsidy for the research capacity of the Faculty of Biology; University of Gdańsk (531-D110-D585-21); TF was supported by Foundation of Faculty of Science; Charles University; JJ was supported by Czech Science Foundation (18-11378 S); MFB and M-AS were supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-PROTAX); The costs of sample collection in Poland were financed by National Science Cantre, Poland (2015/18/A/NZ8/00149) granted to M-AS.
All methods were performed in accordance with the relevant guidelines and regulations.
Ethics approval of research involving plants
The authors declared that experimental research works on the plants described in this paper comply with institutional, national, and international guidelines. Field studies were conducted in accordance with local legislation and permission were granted for sampling from relevant authorities.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Minasiewicz, J., Zwolicki, A., Figura, T. et al. Stoichiometry of carbon, nitrogen and phosphorus is closely linked to trophic modes in orchids. BMC Plant Biol 23, 422 (2023). https://doi.org/10.1186/s12870-023-04436-z
- Trophic modes