Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Integrated analysis of zone-specific protein and metabolite profiles within nitrogen-fixing Medicago truncatula-Sinorhizobium medicae nodules

  • Aaron J. Ogden,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Molecular Plant Science Program, Washington State University, Pullman, Washington, United States of America, Institute of Biological Chemistry, Washington State University, Pullman, Washington, United States of America

  • Mahmoud Gargouri ,

    Contributed equally to this work with: Mahmoud Gargouri, JeongJin Park

    Roles Formal analysis, Methodology, Validation

    Current address: Centre de Biotechnologie de Borj Cédria, BP 901 Hammam-Lif, Ben Arous, Tunisia

    Affiliation Institute of Biological Chemistry, Washington State University, Pullman, Washington, United States of America

  • JeongJin Park ,

    Contributed equally to this work with: Mahmoud Gargouri, JeongJin Park

    Roles Formal analysis, Methodology, Validation

    Current address: W.M. Keck Biomedical Mass Spectrometry Lab, University of Virginia, Charlottesville, VA, United States of America

    Affiliation Institute of Biological Chemistry, Washington State University, Pullman, Washington, United States of America

  • David R. Gang,

    Roles Funding acquisition, Resources, Writing – review & editing

    Affiliations Molecular Plant Science Program, Washington State University, Pullman, Washington, United States of America, Institute of Biological Chemistry, Washington State University, Pullman, Washington, United States of America

  • Michael L. Kahn

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

    kahn@wsu.edu

    Affiliations Molecular Plant Science Program, Washington State University, Pullman, Washington, United States of America, Institute of Biological Chemistry, Washington State University, Pullman, Washington, United States of America

Abstract

Symbiotic nitrogen fixation (SNF) between rhizobia and legumes requires metabolic coordination within specialized root organs called nodules. Nodules formed in the symbiosis between S. medicae and barrel medic (M. truncatula) are indeterminate, cylindrical, and contain spatially distinct developmental zones. Bacteria in the infection zone II (ZII), interzone II-III (IZ), and nitrogen fixation zone III (ZIII) represent different stages in the metabolic progression from free-living bacteria into nitrogen fixing bacteroids. To better understand the coordination of plant and bacterial metabolism within the nodule, we used liquid and gas chromatography coupled to tandem mass spectrometry (MS) to observe protein and metabolite profiles representative of ZII, IZ, ZIII, whole-nodule, and primary root. Our MS-based approach confidently identified 361 S. medicae proteins and 888 M. truncatula proteins, as well as 160 metabolites from each tissue. The data are consistent with several organ- and zone-specific protein and metabolite localization patterns characterized previously. We used our comprehensive dataset to demonstrate how multiple branches of primary metabolism are coordinated between symbionts and zones, including central carbon, fatty acid, and amino acid metabolism. For example, M. truncatula glycolysis enzymes accumulate from zone I to zone III within the nodule, while equivalent S. medicae enzymes decrease in abundance. We also show the localization of S. medicae's transition to dicarboxylic acid-dependent carbon metabolism within the IZ. The spatial abundance patterns of S. medicae fatty acid (FA) biosynthesis enzymes indicate an increased demand for FA production in the IZ and ZIII as compared to ZI. These observations provide a resource for those seeking to understand coordinated physiological changes during the development of SNF.

Introduction

The largest contribution of nitrogen into the biosphere results from symbiotic interactions between legumes and soil bacteria of the family Rhizobiaceae [13]. This process, symbiotic nitrogen fixation (SNF), requires both plants and bacteria to drastically alter their metabolisms in order to establish the unusual conditions needed for SNF [4]. The developmental transition leading to SNF is initiated by nitrogen deficient legumes, which release flavonoids and other diffusible metabolites into the rhizosphere [5]. Rhizobia perceive these signals and respond by excreting Nod factors, lipochitooligosaccharides (LCOs) that interact with plant lysin-motif (LysM) receptors to stimulate the common symbiotic pathway (CSP) [6,7].

The CSP allows rhizobia to invade plant root hairs through tubular plant membrane channels called infection threads (ITs) [8]. Throughout symbiosis, ITs continue to traffic rhizobia into adjacent plant cells where they are released via endocytosis [9]. Endocytotic rhizobia are maintained within a plant-derived membrane, together referred to as the symbiosome, which prevents direct contact between symbionts. Within the infected plant cells, rhizobia begin a metabolic transition from bacteria to bacteroids, a differentiated state that is capable of fixing nitrogen [10].

As infection threads form, the CSP also initiates nodule organogenesis. For galegoid legumes, such as Medicago truncatula, nodule organogenesis involves the formation of a persistent plant meristem [11]. Persistence of a meristem has substantial consequences for both nodule structure and function. Most notably, the meristem produces a cylindrical mature nodule that contains five developmentally and physiologically distinguishable zones [1215]. The first and smallest zone (ZI) is a thin layer of 2–5 cells at the tip of the meristem in which there are few rhizobia, typically contained in ITs [15]. Inward from ZI is the larger infection zone II (ZII), containing ITs and recently endocytized rhizobia. Immediately proximal to ZII is the interzone II-III (IZ), in which the bacteria have partially transitioned into terminally differentiated bacteroids, and the plant cells prepare to support nitrogen fixation by inducing proteins like leghemoglobin and glutamine synthetase, and establishing the needed microaerobic conditions in which O2 levels are of the order of 5–10 nM [10,16,17]. SNF is considered mature and most active in fixation zone III (ZIII). ZIII has a rhizobial population composed primarily of bacteroids and is the largest zone of the nodules studied in this work. ZIII is characterized by proteins involved in fixating and assimilating nitrogen, as well as oxygen binding and amino acid metabolism [1821]. After five or more weeks of growth, a senescent zone, ZIV, appears in the region of the nodule proximal to the root [15,16].

Despite being considered metabolically distinct, the extent to which these zones differ metabolically is undefined. Previous attempts to understand how metabolism is coordinated as SNF develops have probed the physiology of whole-nodules, nodules differing in age, and nodules arrested at particular growth stages and have typically followed the presence of one or a few mRNAs or proteins [22]. These studies generally do not profile the many physiological changes that occur as the nodule tissue differentiates. Recently, Roux et al. (2014) generated transcriptomes of each zone by combining laser capture microdissection (LCM) and Illumina platform-based RNA-seq [23]. While Roux et al. (2014) provided a comprehensive resource of zone-specific transcriptional changes, our understanding of SNFs metabolic development requires a complementary analysis of proteins and metabolites during development.

We report here the use of nano-ultra performance liquid chromatography (UPLC) coupled to label-free mass spectrometry (MS) to catalogue the protein profiles of whole nodules, roots, and three developmental fractions of nodules formed in the M. truncatula-S. medicae plant-bacterial system. Furthermore, we use gas chromatography (GC) in line with MS to create primary metabolite profiles for the three sub-nodule fractions. Unlike RNA, metabolites and proteins cannot be amplified during analysis and thus more tissue is needed than can be easily obtained using LCM. This limitation was resolved by pooling the sections collected from multiple nodules prior to metabolite or protein extraction. Accurate and reproducible capture of the target zones during sectioning was made possible by using the drastic accumulation of pigmented leghemoglobin as a marker. Leghemoglobin is characteristically induced in the IZ [10,16,17] and the transition zone separates ZI and ZII from the nitrogen fixation ZIII. Our methodology was validated by comparing our data with previously described zone-specific biomarkers in closely related indeterminate nodule systems. Furthermore, the high reproducibility and clear statistical distinction between each fraction was shown using multivariate statistical analyses. After validation, we aim to elaborate the coordinated changes in central carbon, fatty acid, and amino acid metabolism between symbionts.

Materials & methods

Plant and bacterial growth conditions

M. truncatula A17 Jemalong seeds were scarified and sterilized prior to germination using a protocol adapted from [24]. Seeds were first soaked in concentrated H2SO4 for five min followed by five washes with sterile distilled H2O (sdH2O). Seeds were then incubated overnight at 4°C in a solution of 100 μg/ml ampicillin, 100 μg/ml kanamycin, and 10 μg/ml rifampicin. Antibiotics were removed the following morning with five sdH20 washes, and the seeds were spread across agar filled petri plates containing yeast mannitol broth [25,26]. Allowing the seeds to vernalize and germinate on rich media validates the sterilization steps, and seeds showing contamination were discarded. To synchronize germination, seeds were vernalized in the dark at 4°C for 48 hr, after which they were relocated to a germination chamber for 24 hr at 20°C. Seedlings were then transferred to sterilized trays containing perlite at a density of 50 seedlings/tray. Perlite trays were then inoculated with either plant nutrient solution (PNS) +3mM urea, PNS, or PNS containing rhizobia for positive, negative, and treatment conditions, respectively [27]. S. medicae used for plant inoculation was cultured at 30°C overnight in minimal mannitol broth supplemented with 0.05% NH4Cl, and added to PNS to a final optical density of 0.1. After inoculation, trays were transferred into insect exclusion tents (Bugdorm 2120F) and grown in a greenhouse at 25°C under a light/dark cycle of 16/8. Plants were watered as needed and fertilized three times a week with PNS or PNS + 3 mM urea for negative and positive controls, respectively.

Tissue harvesting

Whole nodules, primary root, and each sub-nodule fraction was collected 30 days post inoculation (dpi) of M. truncatula with S. medicae. To isolate the three zones, each nodule was hand sectioned under binocular magnification into three parts using the distinct red color caused by leghemoglobin accumulation as a marker of the metabolic shift from non-fixing zones ZI and ZII (Fraction I), the early fixing/ interzone (Fraction II), and the mature fixation zone ZIII (Fraction III). Fraction I was isolated from the distal end of the nodule, approximately one millimeter before the tissue became red. To collect Fraction II, a second excision was made approximately 0.5 millimeters proximal to the initial site of leghemoglobin accumulation. Finally, Fraction III was collected by excising approximately five millimeters from the previous cut. The three Fractions therefore consisted of a white zone Fraction I, a pink zone Fraction II, and a dark red zone Fraction III. Samples were sectioned in this manner, immediately flash frozen in liquid nitrogen, and stored at -80°C prior to protein or metabolite extraction. Approximately 50 plants were used for each replicate resulting in ~50 mg of tissue in all cases prior to homogenization and extraction.

Metabolite extraction and derivatization

All samples were homogenized by mechanical disruption via bead beating (Tissuelyser II Retsch MM400, Qiagen). Samples were returned to liquid nitrogen every 60 seconds to prevent degradation. Following homogenization, samples were extracted with 500 μL of a mixture of methanol, isopropyl alcohol, and water (5:2:2). Samples were then vortexed briefly, incubated for five min in a sonicating water bath (5510 R-MT, Branson), and centrifuged at 16,000 x g for five min. The supernatant was then collected and dried under vacuum without heat. The resulting sample was extracted using 500 μL of a 1:1 acetonitrile and H2O, briefly vortexed, and sonicated for five min. As before, the extracts were centrifuged, and the supernatant dried under vacuum. Derivatization was carried out by first incubating the samples at 30°C for 1.5 hr in five μL of 20 mg methylhydroxylamine hydrochloride (Sigma-M50400) in 500 μL pyridine. Each sample was then derivatized by adding 45 μL of N-methyl-N-trimethylsilylfluoroacetamide with 1% trimethylchlorosilane (Sigma-69478) and incubating at 37°C for 30 min. Samples were then transferred to GC vials and analyzed the same day.

Protein extraction and digestion

Protein extraction was adapted from protocols provided by [28]. First, tissue was homogenized via bead beating (Tissuelyser II Retsch MM400, Qiagen). Each sample was then precipitated with 1.0 ml 10% TCA/acetone and centrifuged at 16,000 x g for 5 min at 4°C. The resulting pellet was then washed with 80% methanol with 0.1 M ammonium acetate and centrifuged at 16,000 x g for 5 min at 4°C. Finally, the pellet was washed with 80% acetone and pelleted at 16,000 x g for 5 min at 4°C. After air-drying the samples to remove residual acetone, each pellet was mixed thoroughly with 0.5 ml/0.1 g starting material of 1:1 phenol (Sigma-P4557)/sdH20. After centrifugation at 16,000 x g for 5 min, the lower phenol phase was transferred to a new tube. Protein was precipitated overnight at 4°C by adding 0.1 M ammonium acetate in methanol at a ratio of 5:1 methanol/phenol. Protein was then pelleted at 16,000 x g for 10 min, air-dried, and resuspended in 50 mM ammonium bicarbonate pH 7.5. Protein was quantified with a Qubit Protein Assay Kit (Life Technologies, Carlsbad, CA) in compliance with the manufacturer’s protocol. Disulfide bonds were reduced using 100 mM dithiothreitol (DTT) at a ratio of 1:10 DTT/sample volume and incubated at 60°C for 30 min. Cysteine bonds were then alkylated with 300 mM iodoacetamide at the same volume ratio for 30 min at room temperature. Finally, protein was digested using trypsin (Promega V511A) at a 1:50 ratio of trypsin/protein, and incubated at 37°C for 8 hr.

GC-MS parameters and processing

The relative abundance of primary metabolites was determined by gas chromatography-time of flight mass spectrometry (GC-TOFMS) using a Pegasus IV (LECO, St. Joseph, MI). Separation was performed using an RXI-5Sil MS with Integra-Guard column (30m x 0.25 mm ID x 0.25 μm film thickness (Restek, GmbH, Bad Homburg, Germany). Column conditions were held initially at 50°C for 1 min following a ramp to 330°C at 20°C/min. All injections were 1 μL in splitless mode with helium as a carrier gas flowing at 1 ml/min. Finally, ChromaTOF software version 4.41 with LECO/Fiehn Metabolomics database was used for processing. Peak intensity values were normalized by the starting fresh weight of each sample. For statistical analysis, features with >50% missing values were omitted and the remaining missing values were imputed using the K-nearest neighbor method [29]. Values were then gLOG transformed according to Metaboanalyst [30]. These normalized data were also used in two-tailed t-tests assuming equal variance and ANOVA.

LC-MS/MS parameters and processing

Nano-ultra performance liquid chromatography-ion mobility time-of-flight tandem mass spectrometry (nano-UPLC-IMS-TOFMS/MS) was used to evaluate the relative abundance of protein within each sample. Chromatographic separation was achieved by injecting 100 ng of protein into a Nanoacquity (Waters, Milford MA) reverse phase C18 column at a flow rate of 0.45 μL/min with a gradient of 5% ACN to 95% ACN in H2O with 0.1% formic acid over 85 min. The Nanoacquity was in line with a Synapt G2S (Waters, Milford MA) equipped with an electrospray ionization (ESI) source operating in positive ion mode. High energy MS/MS function included a transfer collision energy ramp from 15–50 eV. The lock mass using Leu-enkephaline was acquired every 30 seconds. Comprehensive databases for peptide identification were constructed using UniProtKB, Swissprot-Prot, and TrEMBL for both M. truncatula and S. medicae [31,32]. Data processing was performed using TransOmics V1.1 (Waters, Milford MA). Data filtering required that peaks appear in at least two of three biological replicates, i.e., 66.6% consistency and was subsequently normalized by total ion current (TIC). It is known that the population/proportion of rhizobia within the nodule increases from FI to FII and FIII. This became apparent as the contribution of bacterial peptide intensities to the TIC increased from FI to FII and FIII. We accounted for population differences by scaling the intensity of each feature such that their contribution to the TIC was equivalent across samples (S1 Table). For statistical analysis, features with >50% missing values were omitted from analysis, and the remaining missing values were imputed using K-nearest neighbor method [29]. Values were then gLOG transformed according to Metaboanalyst [30]. These normalized data were also used in two-tailed t-tests assuming equal variance and ANOVA.

Results

Highly distinct protein profiles of five plant and bacterial SNF tissues

To evaluate the physiology of tissues associated with SNF, we developed a method for generating comprehensive protein profiles using label free nano-ultra performance liquid chromatography (nano-UPLC) coupled to mass spectrometry (MS). Tissues chosen for analysis included primary root, whole nodule, and three sub-nodule fractions of M. truncatula inoculated with S. medicae (Fig 1). Fraction one (FI) consisted of the colorless nodule tip, distal to the leghemoglobin accumulation indicative of IZ [10,16,17]. It is important to note that our FI contains both infection ZII and meristematic ZI. However, because ZII is much larger than ZI, our FI is likely representative of ZII, and we confirm this expectation later. The second fraction (FII) consists of an approximately 0.5 mm thick section taken at the site of leghemoglobin accumulation and represents the IZ. Fraction three (FIII) was taken from the center of the nodule proximal to the FII section, and represents the fixation ZIII. Finally, we also analyzed S. medicae inoculated primary roots collected 30 dpi.

thumbnail
Fig 1. Zone capture method.

A longitudinal section of a representative M. truncatula-S. medicae root nodule, with red lines indicate the location of incision when separating each zone. Also included in the analysis are whole nodules and primary root.

https://doi.org/10.1371/journal.pone.0180894.g001

Proteomic analysis of these five tissues identified 1,249 proteins with ≥2 unique peptides, consisting of 361 and 888 proteins mapping to S. medicae and M. truncatula, respectively (S2 Table). Table 1 displays examples of correlation between our fraction-specific MS data and previous localizations using mRNA hybridization, immunolocalization, proteomics, reporter gene fusion expression patterns, or transcriptomics [3337]. For example, MtN1 (P93331_MEDTR) and ENOD 20 (MTR_4g130800) were significantly more abundant in FI compared to FII and FIII, a finding that is consistent with observations in the literature [38,39]. Furthermore, our FII protein levels agree with those previously described as IZ specific, including M. truncatula's nodule specific cysteine-rich peptide 333 (A7KHG3_MEDTR) and nodulin 22 (MTR_3g055450), as well as the nodule-specific glycine rich proteins 2A and 2B (MTR_2g042480 and MTR_2g042510, respectively). Finally, proteins known to be most abundant in nitrogen fixing zone III, including nitrogenase components (NifH), oxygen transporters (Lb1, MTR_5g066070), and calcium signaling machinery (MtCaML1, MTR_3g055570) are maximally accumulated in FIII. These results validate our FI, FII, and FIII as accurately representing ZII, IZ, and ZIII, respectively.

thumbnail
Table 1. Comparison of fraction-specific protein abundances with previously characterized zone-specific protein or transcript localizations.

https://doi.org/10.1371/journal.pone.0180894.t001

We hypothesize that previous experiments that collected protein profiles of whole-nodule lysates were in fact sampling the largest zone in the nodule, ZIII, which corresponds to our FIII. We also sought to evaluate whether the protein profiles of our FIII samples would be statistically distinguishable from the whole-nodule sample by principal component analysis (PCA). Fig 2A shows the PCA scoreplot generated from the protein profiles of each sample. Component 1 comprises 44.7% of data variation and is caused by differences between FI, FIII, and root samples, while component 2 shows 20.2% of data variation is predominately caused by differences between the root, FII, and FIII. The only tissues with protein profiles that did not differ substantially were FIII and the whole nodule, which is consistent with the fact that most of the nodule was composed of FIII. These observations also confirm our ability to reproducibly capture tissues with distinct protein profiles using our sectioning method. Analysis of each the three sub-nodule fractions alone demonstrates the majority of data variation was caused by differences between FI and FIII with intermediate values in FII at both metabolite and protein level (Fig 2B and 2C).

thumbnail
Fig 2. Multivariate analysis.

(a) Scoreplot results from PCA using protein profiles of all tissues sampled. (b) PLS-DA scoreplot of fraction-specific protein profiles. (c) PLS-DA scoreplot of fraction-specific metabolite profiles. In all cases dashed circles indicate a 95% confidence region. FI, Fraction I; FII, Fraction II; FIII, Fraction III.

https://doi.org/10.1371/journal.pone.0180894.g002

After confirming that our FI, FII, and FIII fractions accurately and reproducibly represent ZII, IZ, and ZIII, we compared changes in protein profiles of each tissue. Fig 3A and 3B shows the number of proteins with significant intensity differences between consecutive fractions (t-test, p-value < 0.05). As expected, S. medicae proteins were highly up-regulated in the whole nodule as compared to the S. medicae inoculated root. At the level of sub-nodule fractions, 75% of S. medicae proteins with significant differences between FI and FII were upregulated in FII, suggesting a large metabolic transition occurring in FII (IZ according to previous annotation). Similarly, more M. truncatula proteins changed significantly in FII than between any other sub-nodule fractions (Fig 3B). Our analysis identified 70 S. medicae and 189 M. truncatula proteins with significant differences between FII and the surrounding fractions FI and FIII (Fig 3A and 3B). These data provide further evidence of a coordinated and drastic metabolic change occurring in FII. Finally, the number of proteins identified in different ontological categories is illustrated in Fig 3C and displays the depth of our approach.

thumbnail
Fig 3. Summary of protein and metabolite profile data.

(a)-(c) Number of proteins or metabolites showing significant intensity differences in pair wise tissue comparisons (p-value < 0.05). Light grey bars indicate the number of features with greater abundance in distal tissues, while dark grey bars indicate the number of proteins with greater abundance in proximal tissues. (d) Total proteins identified by ontology. Light grey and dark grey bars correspond to S. medicae and M. truncatula proteins, respectively. FI, Fraction I; FII, Fraction II; FIII, Fraction III; TCA, Tricarboxylic Acid Cycle; ETC, Electron Transport Chain; EFs, Elongation Factors; TFs, Transcription Factors; NSPs, Nodule specific proteins.

https://doi.org/10.1371/journal.pone.0180894.g003

Proteome analysis reveals coordinated changes during nodule development

Glycolysis is reciprocally regulated between symbionts as development proceeds.

Proteome data from each tissue was used to examine coordination of carbohydrate metabolism between symbionts during development. The S. medicae phosphofructokinase (PFK) Smed_2106, an enzyme that regulates glycolysis, was 3-fold less abundant in FIII compared to FI, with an intermediate abundance in FII (p-value 0.04). Conversely, the M. truncatula 6-PFK MTR_070s0026 was 47-fold more abundant in FIII compared to FI. Other M. truncatula PFK isozymes showed a similar trend, with MTR_6g090130 increasing 8-fold from FI to FII and remaining abundant in FIII. Interestingly, another rate limiting glycolytic enzyme, pyruvate kinase (PK) (MTR_1g061630), was significantly 10 fold less abundant in FIII compared to FI (p-value 0.02). These two changes are likely to divert the flow of phosphoenolpyruvate away from away from NAD regeneration via pyruvate dehydrogenase and the plant TCA cycle, diverting carbon flow through phosphoenolpyruvate carboxylase (PEPC) (MTR_2g092930 & MTR_2g076670) which were 17-fold and 7-fold more abundant in FIII compared to FI (p-value 0.005 and 0.001, respectively). The pattern of increased M. truncatula PFKs and decreased S. medicae PFK from FI to FIII suggests a coordinated change in carbohydrate metabolism between symbionts. M. truncatula carbohydrate catabolism likely increased in FIII to produce C-4 dicarboxylic acids via PEPC in support of bacteroid metabolism. Decreased S. medicae PFK in FIII suggests a decreased use of carbohydrate carbon sources by bacteroids and further supports the role of C-4 dicarboxylic acids as their major carbon source [40,41].

Spatial localization of S. medicae's transition to TCA dependent carbon metabolism.

Fig 4 displays the relative abundance of S. medicae and M. truncatula TCA cycle enzymes in each nodule fraction. Many S. medicae TCA enzymes, including aconitate hydratase (Smed_3083), succinyl-CoA Ligase (Smed_2942), and malate synthase (Smed_3270), and malate dehydrogenase (Smed_2944), were ≥2 fold more abundant in FII and FIII as compared to FI (p-value <0.05). These observations suggest S. medicae transitions to greater reliance on the TCA cycle in the interzone FII. With the noted exception of malate dehydrogenase, all other M. truncatula TCA cycle enzymes were unchanged or decreased from FI to FIII. These findings indicate S. medicae transitions from carbohydrate to dicarboxylic acid dependent metabolism in FII.

thumbnail
Fig 4.

Relative abundance of TCA cycle enzymes from (a) M. truncatula and (b) S. medicae. Each heat map box corresponds to FI-FIII. The color for each box corresponds to the abundance of a protein expressed as a percentage of the total amount of that protein found in all three fractions.

https://doi.org/10.1371/journal.pone.0180894.g004

Coordinated changes in branched chain amino acid metabolism.

MS identified 26 proteins involved in amino acid metabolism, 21 of which show significant differences in abundance between fractions (ANOVA, p-value <0.05). Of particular interest are proteins necessary for synthesis of branched chain amino acids (BCAAs), which are critical to SNF in closely related systems [4244]. S. medicae's 3-isopropylmalate dehydratase LeuC (Smed_3059) and ketol-acid reductoisomerase IlvC (Smed_1962) were 17- and 6-fold less abundant in FIII compared to FI, respectively (p-value <0.0001). These findings indicate that S. medicae BCAA synthesis decreases during bacteroid development, similar to patterns observed in other SNF systems.

Interestingly, a similar trend is observed in abundance patterns of the M. truncatula BCAA synthesis enzymes 2-isopropylmalate synthase (MTR_1g116500), 3-isopropylmalate dehydratase (MTR_1g045750), and ketol-acid reductoisomerase (MTR_4g079780), that show a 2 fold decrease in FIII compared to FI [45]. Despite these observations, however, our GC-MS metabolite profiling analysis showed that levels of leucine, isoleucine, and valine do not significantly change between fractions, suggesting that the 2-fold decrease in plant BCAA biosynthetic enzymes did not lead to decreased BCAA pools, or that alternative mechanisms exist to maintain BCAA pools between zones (S3 Table).

Dynamic zone-specific accumulation of fatty acid biosynthesis and modification enzymes

Proteomic analysis identified 20 proteins involved in fatty acid (FA) metabolism whose abundance significantly changed between the three fractions (ANOVA, p-value < 0.05). The early steps of S. medicae’s FA biosynthesis are catalyzed by the enzymes FabB (Smed_3449) and Acetyl-CoA carboxylase subunits (Smed_0952 and Smed_0953), which were 15-, 3-, and 500-fold more abundant in FIII compared to FI (p-value<0.05). These observations suggest S. medicae requires increased FA production as SNF develops. The M. truncatula FA elongation enzyme 3-oxoacyl-ACP synthase (MTR_122s0008) was also more abundant in FIII compared to FI (p-value = 0.0001).

In addition to FA biosynthesis, we also observed fraction-specific differences in FA modifying enzymes of the lipoxygenase protein family. Lipoxygenases are known to alter gene expression, exhibit anti-microbial activity, and produce the plant hormone jasmonate [46]. M. truncatula lipoxygenases, annotated Lip1-5, differed in their patterns of accumulation across fractions. Lip1 and Lip2 were highly abundant in FI and decreased significantly across subsequent fractions. Conversely, Lip4 and Lip5 were most abundant in FIII and decreased 9- and 35-fold in FI, respectively. Lip3, however, was maximally abundant in FII, and decreased in abundance in the adjacent FI and FIII fractions (p-value = 0.01). These localizations suggest a distinct function for different lipoxygenases and modified FAs during nodule development. These changes in localization of FA biosynthesis and modifying proteins were consistent with the distribution of the corresponding mRNAs [23].

Highly distinct metabolite profiles of SNF tissues

To provide additional insight into the physiology of nodule development, we generated primary metabolite profiles for each sub-nodule fraction using gas chromatography coupled to mass spectrometry (GC-MS). Our method successfully identified 160 metabolites, many of which showed a significant difference in abundance between different fractions (Fig 3C). These 68 compounds include 8 fatty acids, 18 sugars, 6 dicarboxylic acids, 11 amino acids, and 3 sterols. The transition zone FII contained 82 metabolites with an average intensity that was intermediate between FI and FIII. Fig 2C displays the PLS-DA scoreplot generated from these profiles and identifies metabolites whose abundance correlates with the various fractions [47]. Component one of the PLS-DA scoreplot accounts for 50.0% of the data variation, and was due primarily to differences between FI and FIII.

Glycolysis and TCA metabolite abundances support proteome findings

The glycolytic intermediates fructose-6-phosphate and glucose-6-phosphate were 2- and 6-fold less abundant, respectively, in FIII compared to FI, indicating a greater need for glycolysis in the nitrogen fixing regions of the nodule (p-value <0.02). This was consistent with the observed increase of M. truncatula PFKs in FIII. Additionally, pyruvate was 2-fold less abundant in FIII than FI (p-value = 0.001). Furthermore, the intensity of malate significantly decreased in FIII and FII as compared to FI. The combination of decreased malate with increased levels of S. medicae dicarboxylic acid consuming enzymes, beginning specifically in FII, was consistent with a transition from glycolysis to TCA dependent metabolism. Altogether, the levels of sugar phosphates and dicarboxylic acids were consistent with the proteome changes observed, suggesting a net use of plant glycolysis up to PEP in FIII as compared to FI, accompanied by malate consumption via the TCA cycle in bacteroids.

Discussion

Sub-nodule protein profiles reveal coordinated changes in plant and bacterial metabolism during nodule development

Indeterminate nodules contain conspicuous developmental zones along their longitudinal axis, each of which represents a particular stage in the metabolic development of SNF. While many previous experiments have outlined aspects of SNF metabolism using protein and metabolite profiling methods, they have not focused on the coordinated changes that occur in each zone as the nodule develops. Transcriptomes have recently been generated for each zone using a closely related system, M. truncatula-S. meliloti [23]. Here, we describe a complementary study that applies a label free LC-MS technique to create protein profiles for whole-nodules, three sub-nodule fractions, and roots of the M. truncatula-S. medicae symbiosis. To yield additional insight about the consequences of the changes in protein abundance, we also collected GC-MS primary metabolite profiles from each sub-nodule fraction.

The accuracy of our zone isolation method was validated using previously described biomarkers for each zone. Furthermore, the ability of our methods to capture distinct protein and metabolite profiles of the different fractions was shown using multivariate analysis including PCA and PLS-DA. We then demonstrated that we were able to follow important changes in relative abundance of plant and bacterial enzymes across the three sub-nodule zones with an emphasis on carbohydrate, central carbon, fatty acid, and amino acid metabolism. For each of these metabolic processes, the spatial patterns of sub-nodule protein abundance were largely in agreement with zone-specific transcript levels observed by Roux et al. (2014). Furthermore, the integration of metabolite profiles from each fraction corroborated many of our proteome observations. Fig 5 shows a metabolic pathway model that captures the dynamic and coordinated nature of multiple metabolic processes including glycolysis, TCA, fatty acid biosynthesis, and others. The data set in S2 Table would allow other pathways to be examined.

thumbnail
Fig 5. Comprehensive model of coordinated metabolism of carbohydrates, central carbon, fatty acids, amino acids, and others.

6PFK, 6 phosphofructokinase; ACC1, Acetyl-CoA carboxylase 1; ACC2, Acetyl-CoA carboxylase 2; AcpP, Acyl carrier protein; AKGDH, Alpha ketoglutarate dehydrogenase; Asn, Asparagine; ASNS, Asparagine synthase; Asp, Aspartate; AST, Aspartate transaminase; Cit, Citrate; ENO, Enolase; F1,6BP, Fructose 1,6-bisphosphate; F6P, Fructose 6-phosphate; FabB, 3-oxoacyl-(acyl carrier protein) synthase I; FbpA, Fructose-bisphosphate aldolase; FI, Fraction I; FII, Fraction II; FIII, Fraction III; G2P, Glyceraldehyde 2-phosphate; G3P, Glyceraldehyde 3-phosphate; GAPDH, Glyceraldehyde phosphate dehydrogenase; Glc, Glucose; Glc1P, Glucose 1-phosphate; Glc6P, Glucose 6-phosphate; Gln, Glutamine; Glu, Glutamate; GS, Glutamine synthetase; GSI, Glutamine synthetase I; ICDH, Isocitrate dehydrogenase; Mal, Malate; MDH, Malate dehydrogenase; NifD, Nitrogenase molybdenum-iron protein alpha chain; NifH, Nitrogenase reductase; NifK, Nitrogenase molybdenum-iron protein beta chain; NifN, Nitrogenase molybdenum-cofactor biosynthetic protein; OAA, Oxaloacetate; OGD, Oxoglutarate dehydrognease; PC, Pyruvate carboxylase; PDH, Pyruvate dehydrogenase; PEP, Phosphoenolpyruvate; PEPC, Phosphoenolpyruvate carboxylase; PFK, Phosphofructokinase; PGM, Phosphoglucomutase; PK, Pyruvate kinase; SDH, Succinate dehydrogenase; Suc, Sucrose; TCA, Tricarboxylic acid cycle; TIM, Triose phosphate isomerase; α-KG, alpha ketoglutarate.

https://doi.org/10.1371/journal.pone.0180894.g005

As SNF develops, multiple changes in fraction-specific protein profiles confirmed current paradigms of symbiosis. For example, the role of C-4 dicarboxylic acids as the primary carbon source for terminally differentiated bacteroids was supported by observed decreases in S. medicae glycolysis enzymes and increases in TCA enzymes from FI to FIII. These changes appear to be coordinated with increased accumulation of M. truncatula glycolysis enzymes from FI to FIII, likely to accommodate growing bacteroid use of dicarboxylic acids as their primary carbon source. The pattern of proteins involved in the transition of both plant and bacterial pathways to support the dicarboxylic acid dependency of bacteroids was localized to our equivalent of the IZ, FII, and therefore identified this metabolic shift as a novel characteristic of the IZ. It is interesting to note that M. truncatula pyruvate kinase (MTR_1g061630) was significantly less abundant in FIII compared to FI. This change was different than that of the M. truncatula PEPC isozymes, including MTR_2g092930 and MTR_2g076670, which both showed a significant 17- and 7-fold increase in FIII as compared to FI. This is consistent with the idea that regenerating the NAD+ needed to complete a glycolytic fermentation in the M. truncatula compartments was carried out by diverting phosphoenolpyruvate to oxaloacetate via PEPC and then to malate using malate dehydrogenase. Malate would then be fed to the bacteroids as a dicarboxylic acid source.

We were interested to see fraction-specific (i.e. zone-specific) abundance patterns of fatty acid biosynthetic (FAB) and modifying enzymes. The marked accumulation of S. medicae FAB proteins in FIII as compared to FI strongly suggests that the bacteria are producing fatty acids. Increased FA synthesis was likely caused by requirements for membrane growth characteristic of bacteroids or may be used as a sink for surplus reductant. This idea was supported by significant decreases in free FAs in FIII as compared to FI, including linoleic, lauric, stearic, and lignoceric acid, which we propose were being esterified to glycerol for incorporation into membranes. The fraction-specific patterns of multiple M. truncatula lipoxygenases further support the involvement of modified lipids in nodule development.

Understanding the changes in symbiont metabolism that occur during indeterminate nodule development has been the goal of many projects in the rhizobia research community. Early attempts at achieving this were limited to monitoring the location of one or a small set of gene products within the nodule via histological techniques. More recently, rhizobia mutants that cause nodules to become arrested at particular stages of development have been examined at the protein and transcript level. Advances in protein, transcript, and metabolite extraction methods, coupled with ever increasing instrument sensitivity is only now allowing us to observe the metabolic changes occurring within a wildtype indeterminate nodule. This work represents a substantial contribution toward understanding wild-type symbiont development at the protein and metabolite levels.

Supporting information

S1 Table. Proteome and metabolome dataset normalizations.

https://doi.org/10.1371/journal.pone.0180894.s001

(XLSX)

S2 Table. Proteome dataset of root, FI-FIII, and nodule tissues.

https://doi.org/10.1371/journal.pone.0180894.s002

(XLSX)

S3 Table. Metabolome dataset of root, FI-FIII, and nodule tissues.

https://doi.org/10.1371/journal.pone.0180894.s003

(XLSX)

Acknowledgments

We would like to thank Dr. Svetlana Yurgel for her contributions toward the conceptualization and execution of this work. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [48] partner repository with the dataset identifier PXD006683.

References

  1. 1. Herridge DF, Peoples MB, Boddey RM. Global inputs of biological nitrogen fixation in agricultural systems. Plant Soil. 2008;311(1–2):1–18.
  2. 2. Smil V. Nitrogen and food production: proteins for human diets. AMBIO J Hum Environ. 2002;31(2):126–131.
  3. 3. Long S, Kahn M, Seefeldt L, Tsay Y-F, Kopriva S. Chapter 16: Nitrogen and Sulfur. In: Biochemistry and Molecular Biology of Plants. John Wiley & Sons; p. 711–68.
  4. 4. Oldroyd GE, Downie JA. Coordinating nodule morphogenesis with rhizobial infection in legumes. Annu Rev Plant Biol. 2008;59:519–546. pmid:18444906
  5. 5. Haag AF, Arnold MFF, Myka KK, Kerscher B, Dall’Angelo S, Zanda M, et al. Molecular insights into bacteroid development during Rhizobium–legume symbiosis. FEMS Microbiol Rev. 2012;364–383.
  6. 6. Rolfe BG, Gresshoff PM. Genetic analysis of legume nodule initiation. Annu Rev Plant Physiol Plant Mol Biol. 1988;39(1):297–319.
  7. 7. Gough C, Cullimore J. Lipo-chitooligosaccharide signaling in endosymbiotic plant-microbe interactions. Mol Plant-Microbe Interact MPMI. 2011 Aug;24(8):867–78. pmid:21469937
  8. 8. Ferguson BJ, Indrasumunar A, Hayashi S, Lin Meng-Han, Lin Yu-Hsiang, Reid DE, et al. Molecular Analysis of Legume Nodule Development and Autoregulation. J Integr Plant Biol. 2010 Jan;52(1):61–76. pmid:20074141
  9. 9. Gage DJ. Infection and invasion of roots by symbiotic, nitrogen-fixing rhizobia during nodulation of temperate legumes. Microbiol Mol Biol Rev. 2004;68(2):280–300. pmid:15187185
  10. 10. Mylona P, Pawlowski K, Bisseling T. Symbiotic Nitrogen Fixation. Plant Cell. 1995 Jul;7(7):869–85. pmid:12242391
  11. 11. Patriarca EJ, Tate R, Ferraioli S, Iaccarino M. Organogenesis of legume root nodules. Int Rev Cytol. 2004;234:201–262. pmid:15066376
  12. 12. Foucher F, Kondorosi E. Cell cycle regulation in the course of nodule organogenesis in Medicago. Plant Mol Biol. 2000 Aug;43(5–6):773–86. pmid:11089876
  13. 13. Van de Velde W, Guerra JCP, De Keyser A, De Rycke R, Rombauts S, Maunoury N, et al. Aging in legume symbiosis. A molecular view on nodule senescence in Medicago truncatula. Plant Physiol. 2006;141(2):711–720. pmid:16648219
  14. 14. Monahan-Giovanelli H, Pinedo CA, Gage DJ. Architecture of infection thread networks in developing root nodules induced by the symbiotic bacterium Sinorhizobium meliloti on Medicago truncatula. Plant Physiol. 2006 Feb;140(2):661–70. pmid:16384905
  15. 15. Vasse J, De Billy F, Camut S, Truchet G. Correlation between ultrastructural differentiation of bacteroids and nitrogen fixation in alfalfa nodules. J Bacteriol. 1990;172(8):4295–4306. pmid:2376562
  16. 16. Soupène E, Foussard M, Boistard P, Truchet G, Batut J. Oxygen as a key developmental regulator of Rhizobium meliloti N2-fixation gene expression within the alfalfa root nodule. Proc Natl Acad Sci. 1995;92(9):3759–3763. pmid:7731979
  17. 17. de Billy F, Barker DG, Gallusci P, Truchet G. Leghaemoglobin gene transcription is triggered in a single cell layer in the indeterminate nitrogen-fixing root nodule of alfalfa. Plant J. 1991 Jul 1;1(1):27–35.
  18. 18. Temple SJ, Vance CP, Stephen Gantt J. Glutamate synthase and nitrogen assimilation. Trends Plant Sci. 1998 Feb;3(2):51–6.
  19. 19. Lam H-M, Coschigano KT, Oliveira IC, Melo-Oliveira R, Coruzzi GM. The Molecular-Genetics of Nitrogen Assimilation into Amino Acids in Higher Plants. Annu Rev Plant Physiol Plant Mol Biol. 1996;47(1):569–93.
  20. 20. Udvardi M, Poole PS. Transport and Metabolism in Legume-Rhizobia Symbioses. Annu Rev Plant Biol. 2013;64(1):781–805.
  21. 21. Sulieman S, Tran L-SP. Asparagine: an amide of particular distinction in the regulation of symbiotic nitrogen fixation of legumes. Crit Rev Biotechnol. 2013;33(3):309–327. pmid:22793647
  22. 22. Marx H, Minogue CE, Jayaraman D, Richards AL, Kwiecien NW, Siahpirani AF, et al. A proteomic atlas of the legume Medicago truncatula and its nitrogen-fixing endosymbiont Sinorhizobium meliloti. Nat Biotechnol. 2016 Nov;34(11):1198–205. pmid:27748755
  23. 23. Roux B, Rodde N, Jardinaud M-F, Timmers T, Sauviac L, Cottret L, et al. An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing. Plant J. 2014 Mar 1;77(6):817–37. pmid:24483147
  24. 24. Barker DG, Pfaff T, Moreau D, Groves E, Ruffel S, Lepetit M, et al. Growing M. truncatula: choice of substrates and growth conditions. Medicago Truncatula Handb Samuel Roberts Noble Found Ardmore OK Httpwww Noble OrgMedicagoHandbook [Internet]. 2006 [cited 2014 Jan 21]; Available from: https://www.noble.org/Global/medicagohandbook/pdf/GrowingMedicagotruncatula.pdf
  25. 25. Vincent JM. A manual for the practical study of the root-nodule bacteria. Man Pract Study Root-Nodule Bact [Internet]. 1970 [cited 2014 Feb 27]; Available from: http://www.cabdirect.org/abstracts/19710700726.html
  26. 26. Hoben HJ, Somasegaran P. Comparison of the pour, spread, and drop plate methods for enumeration of Rhizobium spp. in inoculants made from presterilized peat. Appl Environ Microbiol. 1982;44(5):1246–1247. pmid:16346141
  27. 27. McDermott TR, Kahn ML. Cloning and mutagenesis of the Rhizobium meliloti isocitrate dehydrogenase gene. J Bacteriol. 1992 Jul 1;174(14):4790–7. pmid:1320616
  28. 28. Wang W, Vignani R, Scali M, Cresti M. A universal and rapid protocol for protein extraction from recalcitrant plant tissues for proteomic analysis. Electrophoresis. 2006;27(13):2782–2786. pmid:16732618
  29. 29. Stacklies W, Redestig H, Scholz M, Walther D, Selbig J. pcaMethods—a bioconductor package providing PCA methods for incomplete data. Bioinformatics. 2007;23(9):1164–1167. pmid:17344241
  30. 30. Xia J, Psychogios N, Young N, Wishart DS. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009;37(suppl 2):W652–W660.
  31. 31. The UniProt Consortium. Activities at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2013 Nov 18;42(D1):D191–8.
  32. 32. Apweiler R, Martin MJ, O’Donovan C, Magrane M, Alam-Faruque Y, Antunes R, et al. Ongoing and future developments at the Universal Protein Resource. Nucleic Acids Res. 2011;39(Database issue):D214–9. pmid:21051339
  33. 33. Catalano CM, Lane WS, Sherrier DJ. Biochemical characterization of symbiosome membrane proteins from Medicago truncatula root nodules. ELECTROPHORESIS. 2004 Feb 1;25(3):519–31. pmid:14760646
  34. 34. Liu J, Miller SS, Graham M, Bucciarelli B, Catalano CM, Sherrier DJ, et al. Recruitment of Novel Calcium-Binding Proteins for Root Nodule Symbiosis in Medicago truncatula. Plant Physiol. 2006 May 1;141(1):167–77. pmid:16543412
  35. 35. Lefebvre B, Timmers T, Mbengue M, Moreau S, Hervé C, Tóth K, et al. A remorin protein interacts with symbiotic receptors and regulates bacterial infection. Proc Natl Acad Sci. 2010 Feb 2;107(5):2343–8. pmid:20133878
  36. 36. Yoshioka H, Gregerson RG, Samac DA, Hoevens KCM, Trepp G, Gantt JS, et al. Aspartate Aminotransferase in Alfalfa Nodules: Localization of mRNA During Effective and Ineffective Nodule Development and Promoter Analysis. Mol Plant Microbe Interact. 1999 Apr 1;12(4):263–74.
  37. 37. Queiroux C, Washburn BK, Davis OM, Stewart J, Brewer TE, Lyons MR, et al. A comparative genomics screen identifies a Sinorhizobium meliloti 1021 sodM-like gene strongly expressed within host plant nodules. BMC Microbiol. 2012;12:74. pmid:22587634
  38. 38. Gamas P, de Billy F, Truchet G. Symbiosis-Specific Expression of Two Medicago truncatula Nodulin Genes, MtN1 and MtN13, Encoding Products Homologous to Plant Defense Proteins. Mol Plant Microbe Interact. 1998 May 1;11(5):393–403. pmid:9574507
  39. 39. Vernoud V, Journet E-P, Barker DG. MtENOD20, a Nod Factor-Inducible Molecular Marker for Root Cortical Cell Activation. Mol Plant Microbe Interact. 1999 Jul 1;12(7):604–14.
  40. 40. Yurgel SN, Kahn ML. Dicarboxylate transport by rhizobia. FEMS Microbiol Rev. 2004;28(4):489–501. pmid:15374663
  41. 41. Yurgel SN, Kahn ML. Sinorhizobium meliloti dctA mutants with partial ability to transport dicarboxylic acids. J Bacteriol. 2005;187(3):1161–1172. pmid:15659691
  42. 42. Lodwig EM, Hosie AHF, Bourdès A, Findlay K, Allaway D, Karunakaran R, et al. Amino-acid cycling drives nitrogen fixation in the legume–Rhizobium symbiosis. Nature. 2003 Apr 17;422(6933):722–6. pmid:12700763
  43. 43. de las Nieves Peltzer M, Roques N, Poinsot V, Aguilar OM, Batut J, Capela D. Auxotrophy accounts for nodulation defect of most Sinorhizobium meliloti mutants in the branched-chain amino acid biosynthesis pathway. Mol Plant Microbe Interact. 2008;21(9):1232–1241. pmid:18700827
  44. 44. Prell J, White JP, Bourdes A, Bunnewell S, Bongaerts RJ, Poole PS. Legumes regulate Rhizobium bacteroid development and persistence by the supply of branched-chain amino acids. Proc Natl Acad Sci. 2009 Jul 28;106(30):12477–82. pmid:19597156
  45. 45. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. pmid:10592173
  46. 46. López MA, Vicente J, Kulasekaran S, Vellosillo T, Martínez M, Irigoyen ML, et al. Antagonistic role of 9-lipoxygenase-derived oxylipins and ethylene in the control of oxidative stress, lipid peroxidation and plant defence. Plant J. 2011 Aug 1;67(3):447–58. pmid:21481031
  47. 47. Ramadan Z, Jacobs D, Grigorov M, Kochhar S. Metabolic profiling using principal component analysis, discriminant partial least squares, and genetic algorithms. Talanta. 2006 Feb 28;68(5):1683–91. pmid:18970515
  48. 48. Vizcaíno JA, Csordas A, del-Toro N, Dianes JA, Griss J, Lavidas I, et al. 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 2016;44(D1):D447–D456. pmid:26527722