Analyzing the soybean transcriptome during autoregulation of mycorrhization identifies the transcription factors GmNF-YA1a/b as positive regulators of arbuscular mycorrhization
© Schaarschmidt et al.; licensee BioMed Central Ltd. 2013
Received: 26 February 2013
Accepted: 18 June 2013
Published: 18 June 2013
Similarly to the legume-rhizobia symbiosis, the arbuscular mycorrhiza interaction is controlled by autoregulation representing a feedback inhibition involving the CLAVATA1-like receptor kinase NARK in shoots. However, little is known about signals and targets down-stream of NARK. To find NARK-related transcriptional changes in mycorrhizal soybean (Glycine max) plants, we analyzed wild-type and two nark mutant lines interacting with the arbuscular mycorrhiza fungus Rhizophagus irregularis.
Affymetrix GeneChip analysis of non-inoculated and partially inoculated plants in a split-root system identified genes with potential regulation by arbuscular mycorrhiza or NARK. Most transcriptional changes occur locally during arbuscular mycorrhiza symbiosis and independently of NARK. RT-qPCR analysis verified nine genes as NARK-dependently regulated. Most of them have lower expression in roots or shoots of wild type compared to nark mutants, including genes encoding the receptor kinase GmSIK1, proteins with putative function as ornithine acetyl transferase, and a DEAD box RNA helicase. A predicted annexin named GmAnnx1a is differentially regulated by NARK and arbuscular mycorrhiza in distinct plant organs. Two putative CCAAT-binding transcription factor genes named GmNF-YA1a and GmNF-YA1b are down-regulated NARK-dependently in non-infected roots of mycorrhizal wild-type plants and functional gene analysis confirmed a positive role for these genes in the development of an arbuscular mycorrhiza symbiosis.
Our results indicate GmNF-YA1a/b as positive regulators in arbuscular mycorrhiza establishment, whose expression is down-regulated by NARK in the autoregulated root tissue thereby diminishing subsequent infections. Genes regulated independently of arbuscular mycorrhization by NARK support an additional function of NARK in symbioses-independent mechanisms.
KeywordsAffymetrix annexin GeneChip autoregulation arbuscular mycorrhiza CCAAT-binding transcription factor NF-Y Rhizophagus irregularis Glycine max (soybean) quantitative RT-PCR split-root system
Plants have a long success story in hosting microsymbionts in their roots to improve their supply with mineral nutrients, particularly the two important macronutrients phosphate and nitrogen. The arbuscular mycorrhiza (AM) symbiosis, an interaction of plants with fungi of the phylum Glomeromycota , probably co-evolved with the early land plants around 450 million years ago (for overview see ). Nowadays, the majority of land plants can form an AM symbiosis that is characterized by the exchange of phosphate against monosaccharides . Moreover, mycorrhizal plants can also benefit from improved availability of other minerals and water, and from induced abiotic and biotic stress tolerance, all contributing to higher plant biodiversity and productivity of ecosystems . AM fungi are obligate biotrophs colonizing the root cortex of plants where they can form, depending on the plant and fungal species, inter- and intracellular hyphae, highly branched intracellular hyphae called arbuscules and/or intracellular hyphal coils, and vesicles that serve as storage organs (for overview see ).
In addition to the wide-spread AM interaction, a few plant families including leguminous and actinorhizal plants have evolved the capability to interact with nitrogen-fixing bacteria like rhizobia and Frankia, respectively. The bacterial symbionts are hosted intracellularly in specialized organs called root nodules. Studies on plant genes involved in the establishment of these intracellular root-microbe symbioses indicated a common evolutionary origin. This led to the assumption that genes involved in the AM symbiosis might have been recruited to allow the interaction with nitrogen-fixing bacteria (for review see [6–9]).
In leguminous plants, a whole set of genes are known to be essential to successfully establish the AM symbiosis and the legume-rhizobia symbiosis, referred to here as nodulation (for review see [7–10]). Establishment of both endosymbioses is initiated by an intense signal exchange between the partners. Nod and Myc factors contain lipochitooligosaccharide signals that are perceived by plant receptor kinases (RKs): in the case of nodulation by LysM RK(s) and in the case of an AM symbiosis by a still unknown receptor, which somehow interacts with the Nod factor receptor complex . The signals are further processed by an early common signaling pathway including a plasma membrane-bound leucine rich repeat (LRR) RK that was found to be also essential for actinorhiza formation (for review see ). Activation of nuclear cation channels and induction of nucleoporin(-related) proteins are involved in generating specific oscillations of the Ca2+ concentrations in the nucleoplasm and perinuclear cytoplasm in both, AM and nodulation. In both symbioses, the Ca2+ signal is translated by a Ca2+ calmodulin-dependent protein kinase (CCaMK) that activates in interaction with other proteins specific transcription factors (TFs) (for overview see [8–10]).
To reduce carbon losses and to maintain the mutualistic character of the symbioses, the plant tightly controls the infection with both heterotrophic endosymbionts. One common regulatory mechanism that limits the number of successful infection events is called autoregulation (autoregulation of nodulation (AON); autoregulation of mycorrhization (AOM)) (for review see [12–16]). The general autoregulation system of legumes comprises systemic and long-distance feedback inhibition initiated by early signals of the plant-microbe interaction suppressing subsequent infections. The key signal mediator of autoregulation is a CLAVATA1 (CLV1)-like RK, hereafter referred to as NARK (Nodulation Autoregulation Receptor Kinase) . Grafting and split-root experiments revealed that NARK acts in the shoot, limiting infections systemically in the entire root system [18, 19]. Mutant plants with defective NARK are characterized by a supernodulating phenotype and a nitrate- and acid-tolerant nodulation [20, 21] and can also exhibit increased mycorrhizal colonization and higher arbuscule abundance [22–25]. However, in the AM symbiosis, the effect of autoregulation in the wild-type is often less obvious, because the fungus can spread longitudinally within the root cortex without forming new infections once it has entered the root. Nevertheless, the autoregulation effect becomes apparent in split-root plants in which colonization caused by previous and subsequent infections can be distinguished [26, 27].
Nod factor application and cross-infections with rhizobia and AM fungi demonstrated initiation of the general autoregulation system by common early signals . Corresponding to the Arabidopsis thaliana CLV system, CLE peptides are suggested as root-derived signals activating NARK in the shoot (for review see [14, 15]). However, AM-induced CLEs are so far not described. The shoot-derived inhibitor (SDI), acting downstream of NARK in AON, has been characterized biochemically as a heat-stable, ethanol-soluble, low-molecular weight molecule which is unlikely an RNA or protein [29, 30]. A putative receptor of SDI in the root might be encoded by TML . NARK is described to affect phytohormone balances including reduction of the shoot-to-root auxin transport and of jasmonic acid biosynthesis in the shoot [32–34]. Recently, it was shown that NARK also affects the ubiquitin-dependent protein degradation pathway via regulation of GmUFD1a .
Although AON and AOM share common elements, not all nark mutants characterized by supernodulation exhibit an equivalent lack of AOM . This indicates the existence of additional factors that are specific in controlling the establishment of both endosymbioses. In the present study, transcriptional changes occurring during mycorrhization in soybean (Glycine max) plants were investigated. To find NARK-regulated genes involved in AOM, we performed Affymetrix GeneChip analyses with split-root plants of Bragg wild-type and nark mutants partially inoculated with the AM fungus Rhizophagus irregularis. Selected NARK-regulated, AM-dependent and AM-independent, genes were verified using RT-qPCR. Among these genes, we identified two putative CCAAT-binding TF subunits, named GmNF-YA1a and GmNF-YA1b, that were downregulated by NARK in an AM-dependent manner. The CCAAT sequence is one of the most common cis-acting elements and was found in approximately 30% of 502 eukaryotic promoters  (for overview see [38–40]). In higher eukaryotes, CCAAT-boxes are found in all types of promoters, including constitutively expressed and inducible promoters. Proteins that bind to the CCAAT-box or to the inverted box have been characterized inter alia for mammals, yeast, filamentous fungi, and plants, and include proteins of the nuclear factor Y (NF-Y) family. DNA binding activities of NF-Ys, which are also called CBFs (CCAAT-binding factors) or HAPs (heme activating proteins), require a high conservation of the five nucleotides. NF-Ys usually act as heterotrimers. Two of the subunits (NF-YB, NF-YC/CBF-A, CBF-C/HAP3, HAP5) form a stable heterodimer, which can then interact with the third subunit (NF-YA/CBF-B/HAP2) allowing binding to the promoter region. In case of yeast, an additional acid subunit, HAP4, activates the HAP2,3,5 complex. In plants, NF-Ys can function in diverse processes (for overview see [41, 42]) and a few NF-Y subunits have previously found to be upregulated during root endosymbioses [43–53]. Though, no NF-YA subunit that is involved in AM symbiosis or suppressed during root symbioses has been described so far.
Results and discussion
Experimental setup and mycorrhization phenotype of wild-type and nark mutant plants
To ensure that the autoregulation system was activated at the time-point of first harvest, some of the plants were then inoculated at the second root-part and mycorrhization was analyzed after another 3 weeks. At this time-point, the mRNA level of RiBTub1 in the initially inoculated root-parts was significantly higher in both nark mutants compared to the wild-type (Figure 2c; 'previously inoculated'). No RiBTub1 transcript accumulation was detected in corresponding root-parts of control plants. In the later inoculated root-parts, the RiBTub1 mRNA level indicates that subsequent infections led to approximately 40% lower colonization rates in wild-type plants previously inoculated at the other root-part compared to wild-type plants without previous infection (Figure 2c; 'subsequently inoculated'). In contrast, both nark mutants, nts382 and nts1007, showed a defective autoregulation. Here, the RiBTub1 transcript level was not affected by previous inoculation. For validation of NARK-regulated gene expression, the experiment was repeated twice (split-root experiment II and III) reflecting the mycorrhization phenotypes and the AOM effect found in the split-root experiment I (Figure S2 in Additional file 1).
For the nark mutant nts1007 a missing autoregulation has been previously described for nodulation as well as for mycorrhization [17, 22]. This mutant has a nonsense mutation that truncates the RK protein at glutamine residue 106 (Q106*) eliminating most of the LRRs and the entire kinase domain . In contrast, the nts382 mutant (Q920*) has a nonsense mutation in the C-terminal part of the kinase domain  and might still exhibit some kinase activity. Previous studies on the En6500 mutant (K606*) and some other uncharacterized nark mutants of the soybean cv. Enrei indicated a lacking AOM in these supernodulating mutants, suggesting that a lower NARK activity is sufficient to trigger AOM than AON . However, in the present study, we found the Bragg nark mutant nts382 to be, in addition to the defective AON , also severely affected in AOM (Figure 2, Figure S2 in Additional file 1). Because in nts382 the truncated NARK contains even more of the kinase domain than in En6500, differences in AOM between nark mutants of Bragg and of Enrei might somehow be rather cultivar-related.
Because nark mutants are characterized by a nitrogen tolerant nodulation, we additionally analyzed the mycorrhization phenotype under different phosphate and nitrogen supply conditions (Figure S3 in Additional file 1). In wild-type plants, AM fungal colonization was suppressed by increasing fertilization with phosphate. In addition, both nark mutants showed a decrease in mycorrhizal colonization in response to increasing levels of phosphate (Figure S3a in Additional file 1). Under the tested conditions, however, a suppression of the mycorrhizal colonization by nitrogen was not found - neither in nark mutants nor in the wild-type (Figure S3b in Additional file 1). The lower mycorrhization under strong nitrate limitation might be caused by reduced carbon supply of the fungus as plants showed chlorosis and growth reduction (Figure S3c in Additional file 1).
Analyzing the soybean transcriptome during AOM by the Affymetrix GeneChip
Transcript profiling using the Affymetrix GeneChip for soybean was conducted with material of Bragg wild-type and of the nark mutant nts1007, both genotypes either non-inoculated or partially inoculated with R. irregularis. The nark mutant nts1007 was previously shown to be defective in AOM  and used for Affymetrix analysis of nodulated seedlings . To find genes regulated by AM and NARK, the Affymetrix data were screened by (multiple) pairwise comparisons using the dChip software  as described in the Materials and methods section.
In contrast to the strong changes in gene expression of mycorrhizal root tissue and to a previous analysis on systemic changes in mycorrhizal M. truncatula , only some soybean genes were found to show a systemic AM-response in shoots and/or non-colonized root-parts of mycorrhizal plants under the given criteria (Figures 3b and 3c). The overall expression patterns of most of these genes, however, hardly indicate a putative regulation by NARK (data not shown). Thus, to find additional candidates potentially regulated by NARK, the criteria of pairwise comparisons were weakened and putative NARK-regulated genes controlled in a rather AM-independent manner were included. The Affymetrix data were screened again by performing multiple pairwise comparisons with dChip and by using the clustering software CLANS . A distinct gene cluster reflecting a putative regulation by NARK was not found by CLANS. However, using CLANS to screen the genes for a specific expression pattern, we found mostly the same genes as found by the multiple pairwise comparisons with dChip. After screening, candidates for further analysis were selected manually due to their overall expression pattern within all samples that indicate a putative AM-dependent or AM-independent regulation by NARK (further information is provided in Table S2 in Additional file 2). Among these candidates, several genes were identified with a putative function as early nodulin, TF, zinc finger protein, cytochrome P450, calmodulin, pathogenesis-related protein, and with transport function including a putative upregulated amino acid transporter (Glyma19g22590).
NARK-regulated candidate genes validated by RT-qPCR
Validated NARK regulation of soybean genes in R. irregularis-inoculated plants.
Affymetrix probe set(s)
Predicted gene function/Name
Ornithine acetyl transferase
Stress-induced receptor-like kinase GmSIK1
DEAD box RNA helicase
Symbioses-dependent control of gene expression by NARK might to some extent occur in specific cell-types and/or transiently only. Transient regulation could particularly take place in cases of induction or de-novo synthesis of gene products. A specific temporal- and/or spatial-restricted NARK regulation might also explain the low number of genes identified in our study. In contrast to a previous analysis of nodulated seedlings , no suppression of JA-biosynthesis or -response genes was found by AM or by NARK. This might be due to the different developmental stages of the plants analyzed. Even in the work by Kinkema and Gresshoff , the differences between non-inoculated wild-type and nts1007 were less pronounced with increasing age of plants, showing mostly no changes in 17-day-old plants. Such transiently and developmentally controlled changes might be less reflected in our 27-day-old soybean plants. However, whether the changes previously found in nodulated soybeans were not detectable in mycorrhizal plants because they are cell-type specific, transiently and/or developmentally regulated or because they are specific for nodulation remains unknown.
NARK-regulated genes lacking obvious regulation by AM
Transcript accumulation of the stress-induced receptor-like kinase GmSIK1 [UniGen:EU669879] in mycorrhizal plants was on average around four times lower in roots and shoots of wild-type compared to both nark mutants (Figure 4b). However, due to a strong variation in non-mycorrhizal wild-type, no significant AM dependency of NARK regulation was detected. Except from shoot tips that contain relatively high transcript levels, transcript accumulation of GmSIK1 increased with increasing age of above- and below-ground tissue (Figure S6h in Additional file 1), indicating a role of this abiotic stress-induced RK in processes linked to plant development or senescence.
According to our RT-qPCR data, transcript levels of the other putative NARK-suppressed genes were overall around 30% to 50% reduced in wild-type compared to nark mutant tissue. One of the genes that were rather slightly downregulated by NARK in roots and shoots is a predicted DEAD box RNA helicase (Glyma17g09270) (Figure 4c). Gene transcripts were highest in tips of shoots and roots (Figure S6j in Additional file 1). Glyma17g09270 is targeted by the two Affymetrix probe sets GmaAffx.68580.1.S1_at and GmaAffx.46141.1.S1_at. The latter, which as well reflects NARK regulation of target gene(s), might also detect the homologous gene Glyma05g02590. RT-qPCR analysis does not indicate gene regulation by NARK for this homolog (data not shown). DEAD box RNA helicases are widely occurring proteins with diverse activities including ATP binding and ATPase activity as well as RNA binding and unwinding or transport activity and are associated with nearly all processes that involve RNA, including development and stress responses [70, 71]. Moreover, two genes of unknown function, Glyma10g35000 and Glyma07g36986, were downregulated by NARK in roots (Figures S5b, c in Additional file 1). Glyma07g36986 mRNA levels were also slightly lower in shoots of mycorrhizal wild-type compared to both nark mutants.
In general, AM-independent gene regulation by NARK indicates an additional role of NARK in other, symbioses-independent processes as previously been proposed . Differences in polar auxin transport or JA biosynthesis are also found between non-inoculated nark mutant and wild-type plants [22, 32–34]. This indicates that even non-inoculated plants exhibit a certain NARK activity. Indeed, symbioses-independent phenotypes of nark mutants of L. japonicus and M. truncatula [72, 73] support a function of NARK in root growth and/or lateral root formation. Part of the NARK-regulated genes identified in the present study was characterized by a development-dependent expression pattern. Hence, such genes might be involved in NARK-controlled signaling in plant development and/or in plant stress responses not linked to AM formation. However, even genes lacking obvious AM dependency in NARK regulation could play a role in AM symbiosis and in AOM as indicated by gene expression data of the putative annexin Glyma15g38010.
GmAnn1a as a differentially NARK- and AM-regulated annexin
GmNF-YA1a/b as AM-dependently NARK-regulated TF subunits
BLASTing the soybean genome  revealed that the target sequence of GmaAffx.40657.1.S1_at strongly aligns with transcript sequences of two homologous genes encoding putative CCAAT-binding TFs: Glyma03g36140 and Glyma19g38800. Both have an overall sequence identity of 90% to 94% and 86% to 91% on predicted coding sequence level and amino acid level, respectively (Figure S8 in Additional file 1). Ranges are due to alternative transcripts of Glyma19g38800. The putative amino acid sequences of both genes contain conserved protein domains supporting a function as CCAAT-binding TF of the NF-YA family (Figure S8b in Additional file 1). Thus, the homologous TF genes corresponding to GmaAffx.40657.1.S1_at were designated as GmNF-YA1 with GmNF-YA1a and GmNF-YA1b for Glyma03g36140 and Glyma19g38800, respectively. RT-qPCR analysis with gene-specific primers confirmed that both homologs, GmNF-YA1a and GmNF-YA1b, were AM-dependently suppressed by NARK in roots to a similar extent (Figure 6d, e).
The soybean transcription factor database v2.0 of the Center for Bioinformatics (CBI) at Peking University [80, 81] contains in total 83 NF-Y proteins including 21 of the NF-YA family, 40 of the NF-YB family, and 22 TFs of the NF-YC family. According to the soybean Affymetrix GeneChip annotation of the NSF-funded project 'Gene Networks in Seed Development' , the array contains 36 probe sets targeting at least 34 different putative NF-Y TF genes (11x NF-YA, 10x NF-YB, 13x NF-YC). Our Affymetrix data did not indicate that one of the other putative NF-Y genes is NARK-regulated in autoregulated roots like GmNF-YA1a/b (Table S3 in Additional file 2). One putative NF-YB (Glyma20g00240, targeted by Gma.8502.1.S1_at) exhibited a higher mRNA level in non-mycorrhizal wild-type roots compared to mycorrhizal root-parts, however, root signals of non-mycorrhizal nts1007 were similar to that of the mycorrhizal wild-type contradicting a role in AOM (Figure S9a in Additional file 1). Another putative NF-YB (Glyma02g46970, targeted by Gma.12719.1.S1_at) showed slightly lower signals in shoots of mycorrhizal wild-type compared to shoots of non-mycorrhizal wild-type and of nts1007 (Figure S9b in Additional file 1). In contrast to the AM-dependently downregulated GmNF-YA1a/b genes, two putative NF-Y genes were found to be locally inducible by AM: Glyma09g07960 (Gma.3174.1.S1_at; NF-YA family) had slightly higher signals and Glyma12g34510 (GmaAffx.66953.1.S1_at; NF-YC family) was strongly induced in mycorrhizal root-parts (Figures S9c,d in Additional file 1). The closest homologs to GmNF-YA1a/b, Glyma02g35190 and Glyma10g10240 , appeared in our experimental setup not to be AM-induced, but developmentally regulated (Figure S10 in Additional file 1).
In addition to the here identified NF-Y genes regulated in roots of mycorrhizal plants, few NF-Ys were previously described to be locally induced by root endosymbionts. Most of them have been shown to be upregulated in nodulation, including the NF-YA genes Glyma12g34510, Glyma02g35190, and Glyma10g10240 of soybean , MtHAP2-1 and MtHAP2-2 of M. truncatula [44–46, 51], and LjCBF-A01/LjNF-YA1 and LjCBF-A22 of L. japonicus [43, 52], the NF-YB1 gene from L. japonicus , as well as the NF-YC genes PvNF-YC1 and PvNF-YC2 of Phaseolus vulgaris [50, 53] and MT007765 of M. truncatula . AM-inducible NF-YB (Mtr.4282.1.S1_at) and NF-YC genes (MtCBf1/Mtr.51511.1.S1_at and MtCbf2/Mtr.16863.1.S1_at) were previously found in M. truncatula . Promoter studies indicate a rather general role for these genes in the coordination of AM fungal colonization during all stages of AM . However, AM-induced NF-YA genes or other NF-Y genes downregulated during mycorrhization have so far not been described.
RNAi-mediated downregulation of GmNF-YA1a/b decreases AM
Transcript analysis of NF-YA and phosphate transporter genes in A. rhizogenes-transformed roots
Glyma10g10240 (putative NF-YA)
Glyma02g35190 (putative NF-YA)
Glyma13g08720 (phosphate transporter)
Glyma14g28780 (phosphate transporter)
Glyma14g36650 (phosphate transporter)
Simultaneously to reduced GmNF-YA1a/b mRNA accumulation, GmNF-YA1-RNAi expressing roots had on average a by 55% reduced colonization rate than the empty vector controls as determined by transcript analysis of the fungal marker gene RiBTub1 and by microscopic analysis. Such a reduced AM fungal colonization upon suppression of GmNF-YA1 mRNA accumulation occurred in both nark mutants as well as in wild-type plants (Figures 7a, b). Microscopy revealed that the mycorrhization was in a relatively early stage of establishment and that the rate of root colonization was mainly a consequence of individual infection events instead of excessive longitudinal spreading of the fungus after penetration into the root (Figure 7c). Thus, the degree of colonization indicates a lower number of infection events in RNAi plants than in the empty vector controls. To additionally assess the function of the AM symbiosis in the transformed roots, we measured transcript levels of AM-inducible phosphate transporters. In M. truncatula, MtPT4 represents a phosphate transporter gene that is specifically induced by AM and indispensable for the symbiosis [83, 84]. Soybean phosphate transporters were previously analyzed identifying Glyma13g08720, Glyma14g28780, and Glyma14g36650 as closest homologs to MtPT4  and as AM-inducible genes . On average, the three AM-inducible phosphate transporter genes were suppressed by approximately 60% in roots expressing the GmNF-YA1-RNAi construct compared to the empty vector control (Table 2). These data clearly indicate a function of GmNF-YA1a/b as positive regulators in AM formation.
The decreased formation of AM in GmNF-YA1-RNAi roots is in line with the gene regulation of GmNF-YA1a/b observed in this study: expression of GmNF-YA1 was highest in non-inoculated roots (Figure 6c) and was systemically suppressed in wild-type roots during an AM interaction (Figure 6) thereby contributing to AOM leading to a reduced mycorrhizal colonization. Moreover, RNAi-mediated suppression of GmNF-YA1a/b in roots reduced AM formation not only in the wild-type but also in the nark mutant background pointing to the fact that GmNF-YA1a/b acts downstream of NARK.
Overall, our gene expression analysis in split-root experiments provides for the first time insights into transcriptional regulation by NARK during the autoregulation of arbuscular mycorrhization - a systemic regulatory mechanism controlling the most widespread and ancient interaction between vascular terrestrial plants and microbes. By RT-qPCR analysis of plant material of independent experiments, we verified nine genes to be regulated by NARK. Gene expression of annexins indicate putative involvement in AM establishment and possibly also in regulation of AM. Furthermore, we identified NF-YA genes regulated during AM, including two NARK-regulated ones named GmNF-YA1a and GmNF-YA1b. So far, no information on a role of NF-YAs during AM symbiosis or autoregulation was available. By our transcript and functional gene analysis, we designate the newly identified TF subunits GmNF-YA1a/b as positive regulators of mycorrhization, which might act as targets for NARK-mediated AOM to restrict new infections events.
Materials and methods
Plant material, germination of seeds, and growth conditions
Seeds of soybean cv. Bragg (wild-type and two allelic nark mutant lines, nts382 and nts1007) were surface-sterilized with 1.5% sodium hypochlorite for 5 min. After washing with distilled water, seeds were germinated in wet expanded clay of 2 mm to 5 mm particle size (Original LamstedtTon; Fibo ExClay, Lamstedt, Germany). The same substrate was used for inoculum production and all experiments described below. If not mentioned otherwise, plants were grown in a plant growth chamber under 16 h light (200 µmol m-2 s-1) at 26°C and 8 h dark at 22°C with a relative humidity of 60% to 65%. All plants were watered with distilled water and fertilized with Long Ashton nutrient solution containing 20% of the regular phosphate content .
Fungal material and inoculation of plants
The AM fungus R. irregularis (formerly Glomus intraradices ) Schenk & Smith isolate 49  was enriched by previous co-cultivation with leek (Allium porrum cv. Elefant) in a greenhouse. As inoculum, the freshly harvested substrate was used. Roots or root-parts of soybean plants were inoculated with R. irregularis by careful removal of the previous substrate and transfer into expanded clay containing 10% to 20% (v/v) R. irregularis-inoculum. Non-mycorrhizal roots or root-parts were transferred in the same way into pure expanded clay.
In total, three independent split-root experiments (I-III) were performed. In each experiment, Bragg wild-type, nts382, and nts1007 plants were cultivated in parallel as follows: after 6 days of germination, the main root of the seedlings was cut off around 2 cm below the shoot basis, and seedlings were continued to be cultivated under the conditions described. After 2 days, plants were transferred to the split-root system by dividing the arising lateral roots onto two individual pots connected to each other (each 9 × 9 × 9.5 cm) (Figure S1 in Additional file 1). In doing so, mycorrhizal plants were inoculated on one root-part with 20% (v/v) inoculum. The other root-part and root-parts of control plants were not inoculated (Figure 1). Nineteen days after inoculation, 50% of the plants were harvested.
To check plants for AOM activation at the first harvest time-point, the non-harvested plants were transferred to bigger pots (each 11 × 11 × 11.5 cm). In doing so, the previously non-inoculated root-part of mycorrhizal plants and one root-part of control plants were inoculated with 10% (v/v) inoculum (Figure 1 and Figure S1 in Additional file 1). Plants were harvested 21 days after subsequent inoculation.
In experiment II and III some additional control and mycorrhizal plants were cultivated in bigger pots without subsequent inoculation and harvested 40 days after inoculation, and above- and below-ground plant tissues were harvested as described below.
Root transformation experiments
For RNAi suppression of GmNF-YA1a/b, a fragment of GmNF-YA1 of 428 bp size was amplified using PCR Super Mix High Fidelity (Invitrogen, Carlsbad, CA, USA) (for primer binding sites see Figure S8a in Additional file 1). To allow directional cloning of the fragment in sense and antisense direction, the restriction sites BamHI, SpeI and SwaI, AscI, respectively, were attached during PCR (for primer sequences see Table S4 in Additional file 1). Gel-purified fragments were ligated into the pGEM-T Easy vector (Promega, Madison, WI, USA). Complementary sense and antisense arms, showing 100% nucleotide identity to Glyma03g36140, were cloned into the pRNAi vector  under the CaMV 35S promoter. The resulting construct was cloned into the pRedRootII vector that contains the gene for DsRed1 in the T-DNA region (kindly provided by R. Geurts, Wageningen, The Netherlands) using KpnI and PacI, and finally transformed into A. rhizogenes K599 . As control, A. rhizogenes K599 cells were transformed with the empty pRedRootII vector.
Root transformation was performed according to the protocol published by Kereszt et al.  using a bacterial suspension that was injected into the hypocotyl close to the cotyledonary node of 5-day-old-seedlings of soybean wild-type, nts382, and nts1007. All seedlings were kept in a humid chamber under 100% humidity at 12 h light (200 µmol m-2 s-1, 28°C) and 12 h dark (25°C). After 3 and 4 weeks, newly developed hairy roots were screened twice for fluorescence of DsRed1 using a fluorescence stereomicroscope equipped with a DsRed filter (Leica 409 MZ FLIII; Leica Microsystems, Wetzlar, Germany). Two to four hairy roots that showed the strongest fluorescence were kept, and the other hairy or wild-type roots were removed. Hairy roots were covered with expanded clay and plants were further cultivated in a greenhouse (16 h sun light supplemented with artificial light, 22°C; 8 h dark, 20°C) acclimating them slowly to lower humidity. One day after second screening, all plants were inoculated with 15% (v/v) inoculum, further cultivated in the greenhouse, and harvested 3 weeks later.
Harvest of plants and plant tissues, isolation of total RNA
Split-root plants were rapidly taken out of the substrate, which was carefully removed from the root-parts. Root-parts and shoots were separated and immediately flash frozen in liquid nitrogen. Samples were homogenized in liquid nitrogen and stored at -80°C until extraction. Roots of A. rhizogenes-transformed plants were harvested in the same way. For staining of fungal structures (see below), a median cross-section of each root sample of approximately 1 cm length was taken before freezing.
Different tissues of some non-inoculated plants of split-root experiment II and III were harvested 7 weeks after sowing. For that purpose, shoots were rapidly dissected into shoot tips, young (sink) leaves, fully developed (source) leaves, yellowing (senescing) leaves, and stalks. For each plant, roots were rapidly removed from the expanded clay and six to ten root-tips, six to ten root-hair zones, and the remaining root material were selected. All tissues were immediately flash frozen, homogenized in liquid nitrogen and stored at -80°C.
Total RNA was isolated using the Qiagen RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) including DNA digestion (RNase-free DNase Set; Qiagen) according to the supplier's protocol.
Staining of fungal structures and analyzing the degree of mycorrhization
Colonization with R. irregularis was analyzed in a representative cross-section of each root sample. Fungal structures in the root pieces of around 1 cm length were stained according to Vierheilig et al.  using 5% (v/v) ink (Sheaffer Skrip jet black; Sheaffer Manufacturing, Shelton, CT, USA) in 2% (v/v) acetic acid. The degree of mycorrhization was determined for at least 50 root-pieces per root sample using a stereomicroscope (Stemi 2000-C; Carl Zeiss, Jena, Germany). In root-transformation experiments, all root-pieces per root sample were analyzed. Pictures were taken with a digital microscope system (VHX-1000 equipped with a VH-Z250R zoom lens, Keyence, Osaka, Japan).
RNA samples of three plants per treatment of the split-root experiment I were analyzed using Affymetrix GeneChips. Synthesis and purification of cDNA; synthesis, labeling, purification, quality control, and fragmentation of cRNA; and hybridization, washing, and scanning of the chips were done by the Affymetrix service partner Atlas Biolabs (Berlin, Germany) according to the supplier's protocols. Data of the binary files were processed and analyzed using the DNA-chip analyzer software package dChip . Arrays were normalized on perfect match and mismatch probe level separately for roots and shoots with a root/shoot array with median overall intensity as baseline array. Afterwards, the model-based expression values were calculated for all arrays together. To find genes that are regulated by AM and NARK, pairwise comparisons and combined pairwise comparisons were performed with dChip. If not mentioned otherwise, the recommended presettings of the software were used. At first, we searched for R. irregularis-regulated genes in wild-type and in the nts1007 mutant. To find genes locally induced or suppressed after R. irregularis-inoculation, we performed the following comparisons, independently for wild-type and nts1007: mycorrhizal versus non-mycorrhizal root-parts of mycorrhizal plants and versus non-mycorrhizal roots of control plants (Figure 3a). To find genes systemically induced or suppressed after R. irregularis-inoculation, comparisons were performed as follows, again each for wild-type and nts1007: (I) shoots of mycorrhizal plants versus shoots of non-mycorrhizal plants (Figure 3b), (II) non-inoculated (autoregulated) root-parts of mycorrhizal plants versus roots of non-inoculated control plants (Figure 3c). The criteria for all comparisons performed with dChip were larger than two-fold changes with P ≤0.1 (unpaired and paired t-tests). For nts1007, only those genes were counted, which showed changes compared to control nts1007 as well as to control wild-type tissue. To avoid under-estimation of the standard error of group means, the measurement error was considered leading to less significant P values than the standard analysis. The identification of additional candidate genes was performed using the clustering software CLANS .
Putative gene functions and categorizations are in accordance with the soybean Affymetrix GeneChip annotation of the NSF-funded project 'Gene Networks in Seed Development' . Candidate genes without annotation were analyzed by sequence homology via searching the NCBI nucleotide collection with the BLAST tool . Selected candidate genes (Table S1 in Additional file 2) were further analyzed by RT-qPCR. Target genes of Affymetrix probe sets were checked by BLAST analysis of the Affymetrix target sequence against the soybean genome  and updated with the latest version (Phytozome v9.0) .
The entire Affymetrix GeneChip data from this study are publicly available at the plant expression database PLEXdb  under the accession number GM53  and at the Gene Expression Omnibus repository  under the GEO Series accession number GSE44685 .
Determination of R. irregularis and soybean transcripts using RT-qPCR
Transcript levels of putative NARK-regulated soybean genes and of R. irregularis β-Tubulin1 (RiBTub1, GenBank:AY326320.1; for primer information see Table S4 in Additional file 1) were quantified by RT-qPCR. First strand cDNA synthesis of 1 µg RNA was performed in a final volume of 20 µL with M-MLV Reverse Transcriptase, RNase H Minus, Point Mutant (Promega) according to the supplier's protocol using oligo(dT)19 primer. RT- samples were prepared in the same way using water instead of enzyme. QPCR primers for soybean candidate genes were designed with the Primer Express software (Applied Biosystems, Foster City, CA, USA) using the corresponding target sequences of the soybean Affymetrix GeneChip (for primer sequences see Table S1 in Additional file 2; primers of validated NARK-regulated genes are also listed in Table S4 in Additional file 1). As reference, soybean ubiquitin GmSUBI-1 transcripts were measured [GenBank:NM_001248971.1] (for primer information see Table S4 in Additional file 1).
For qPCR, 3 μL of 1:10 diluted cDNA were mixed with Brilliant II SYBR Green QPCR Master Mix (Agilent Technologies, Santa Clara, CA, USA), supplemented with the reference dye ROX (final concentration: 30 nM), 2 pmol of forward primer and 2 pmol of reverse primer in a final volume of 10 μL in two to three independent technical replicates. As negative controls, water and RT- samples (1:10 diluted) were used instead of cDNA. Fluorescence of SYBR Green I and of ROX were measured using the Mx3000P and Mx3005P QPCR systems (Agilent Technologies) and the following PCR program: enzyme activation (95°C for 10 min), 40 cycles of amplification (95°C for 30 s, 60°C for 1 min with measurement of fluorescence at the end of this step), followed by dissociation (95°C for 1 min, 60°C for 30 s, heating up to 95°C with a heating rate of 0.1°C s-1 and continuous measurement). Data were evaluated with the MxPro software (Agilent Technologies). For each sample, mRNA levels of target genes were normalized to GmSUBI-1 mRNA by the comparative Ct (2-ΔCt) method .
All data are derived from at least three biological replicates per experiment. The exact number of replicates used for each analysis is given in the corresponding figure legend. If not mentioned otherwise, data are tested for significant differences with P ≤0.05 by the Student's t-tests. In case of multiple testing, Bonferroni correction was applied.
Supplemental material is available for this article. The Affymetrix GeneChip data from this study have been submitted to the plant expression database PLEXdb  and are available under the accession number GM53  and have in addition been deposited in the Gene Expression Omnibus  accessible through GEO Series accession number GSE44685 .
autoregulation of mycorrhization
nuclear factor Y
reverse transcriptase quantitative polymerase chain reaction
We thank the German Academic Exchange Service (DAAD) for a post-doctoral fellowship for SS at the lab of PMG initiating this cooperation. We greatly acknowledge Brett Ferguson and Alain Tissier for critically reading the manuscript.
- Schüßler A, Schwarzott D, Walker C: A new fungal phylum, the Glomeromycota: phylogeny and evolution. Mycol Res. 2001, 105: 1413-1421. 10.1017/S0953756201005196.View ArticleGoogle Scholar
- Brundrett MC: Coevolution of roots and mycorrhizas of land plants. New Phytol. 2002, 154: 275-304. 10.1046/j.1469-8137.2002.00397.x.View ArticleGoogle Scholar
- Kiers ET, Duhamel M, Beesetty Y, Mensah JA, Franken O, Verbruggen E, Fellbaum CR, Kowalchuk GA, Hart MM, Bago A, Palmer TM, West SA, Vandenkoornhuyse P, Jansa J, Bücking H: Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis. Science. 2011, 333: 880-882. 10.1126/science.1208473.View ArticlePubMedGoogle Scholar
- van der Heijden MGA, Klironomos JN, Ursic M, Moutoglis P, Streitwolf-Engel R, Boller T, Wiemken A, Sanders IR: Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature. 1998, 396: 69-72. 10.1038/23932.View ArticleGoogle Scholar
- Smith FA, Smith SE: Structural diversity in (vesicular)-arbuscular mycorrhizal symbioses. New Phytol. 1997, 137: 373-388. 10.1046/j.1469-8137.1997.00848.x.View ArticleGoogle Scholar
- Kereszt A, Kondorosi E: Unlocking the door to invasion. Science. 2011, 331: 865-866. 10.1126/science.1202342.View ArticlePubMedGoogle Scholar
- Markmann K, Parniske M: Evolution of root endosymbiosis with bacteria: how novel are nodules?. Trends Plant Sci. 2009, 14: 77-86.View ArticlePubMedGoogle Scholar
- Oldroyd GED, Harrison MJ, Paszkowski U: Reprogramming plant cells for endosymbiosis. Science. 2009, 324: 753-754. 10.1126/science.1171644.View ArticlePubMedGoogle Scholar
- Streng A, op den Camp R, Bisseling T, Geurts R: Evolutionary origin of rhizobium Nod factor signaling. Plant Signal Behav. 2011, 6: 1510-1514. 10.4161/psb.6.10.17444.View ArticlePubMedPubMed CentralGoogle Scholar
- Gough C, Cullimore J: Lipo-chitooligosaccharide signaling in endosymbiotic plant-microbe interactions. Mol Plant Microbe Interact. 2011, 24: 867-878. 10.1094/MPMI-01-11-0019.View ArticlePubMedGoogle Scholar
- Maillet F, Poinsot V, Andre O, Puech-Pages V, Haouy A, Gueunier M, Cromer L, Giraudet D, Formey D, Niebel A, Martinez EA, Driguez H, Becard G, Denarie J: Fungal lipochitooligosaccharide symbiotic signals in arbuscular mycorrhiza. Nature. 2011, 469: 58-63. 10.1038/nature09622.View ArticlePubMedGoogle Scholar
- Ferguson BJ, Indrasumunar A, Hayashi S, Lin M-H, Lin Y-H, Reid DE, Gresshoff PM: Molecular analysis of legume nodule development and autoregulation. J Integr Plant Biol. 2010, 52: 61-76. 10.1111/j.1744-7909.2010.00899.x.View ArticlePubMedGoogle Scholar
- Hause B, Schaarschmidt S: The role of jasmonates in mutualistic symbioses between plants and soil-born microorganisms. Phytochemistry. 2009, 70: 1589-1599. 10.1016/j.phytochem.2009.07.003.View ArticlePubMedGoogle Scholar
- Mortier V, Holsters M, Goormachtig S: Never too many? How legumes control nodule numbers. Plant Cell Environ. 2012, 35: 245-258. 10.1111/j.1365-3040.2011.02406.x.View ArticlePubMedGoogle Scholar
- Reid DE, Ferguson BJ, Hayashi S, Lin Y-H, Gresshoff PM: Molecular mechanisms controlling legume autoregulation of nodulation. Ann Bot. 2011, 108: 789-795. 10.1093/aob/mcr205.View ArticlePubMedPubMed CentralGoogle Scholar
- Staehelin C, Xie Z-P, Illana A, Vierheilig H: Long-distance transport of signals during symbiosis. Are nodule formation and mycorrhization autoregulated in a similar way?. Plant Signal Behav. 2011, 6: 372-377. 10.4161/psb.6.3.13881.View ArticlePubMedPubMed CentralGoogle Scholar
- Searle IR, Men AE, Laniya TS, Buzas DM, Iturbe-Ormaetxe I, Carroll BJ, Gresshoff PM: Long-distance signaling in nodulation directed by a CLAVATA1-like receptor kinase. Science. 2003, 299: 109-112. 10.1126/science.1077937.View ArticlePubMedGoogle Scholar
- Caetano-Anolles G, Gresshoff PM: Early induction of feedback regulatory responses governing nodulation in soybean. Plant Sci. 1990, 71: 69-10.1016/0168-9452(90)90069-Z.View ArticleGoogle Scholar
- Delves AC, Mathews A, Day DA, Carter AS, Carroll BJ, Gresshoff PM: Regulation of the soybean-rhizobium nodule symbiosis by shoot and root factors. Plant Physiol. 1986, 82: 588-590. 10.1104/pp.82.2.588.View ArticlePubMedPubMed CentralGoogle Scholar
- Carroll BJ, McNeil DL, Gresshoff PM: A supernodulation and nitrate-tolerant symbiotic (nts) soybean mutant. Plant Physiol. 1985, 78: 34-40. 10.1104/pp.78.1.34.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin M-H, Gresshoff PM, Ferguson BJ: Systemic regulation of soybean nodulation by acidic growth conditions. Plant Physiol. 2012, 160: 2028-2039. 10.1104/pp.112.204149.View ArticlePubMedPubMed CentralGoogle Scholar
- Meixner C, Ludwig-Mueller J, Miersch O, Gresshoff P, Staehelin C, Vierheilig H: Lack of mycorrhizal autoregulation and phytohormonal changes in the supernodulating soybean mutant nts1007. Planta. 2005, 222: 709-715. 10.1007/s00425-005-0003-4.View ArticlePubMedGoogle Scholar
- Sakamoto K, Nohara Y: Soybean (Glycine max [L.] Merr.) shoots systemically control arbuscule formation in mycorrhizal symbiosis. Soil Sci Plant Nutr. 2009, 55: 252-257. 10.1111/j.1747-0765.2009.00358.x.View ArticleGoogle Scholar
- Shrihari PC, Sakamoto K, Inubushi K, Akao S: Interaction between supernodulating or non-nodulating mutants of soybean and two arbuscular mycorrhizal fungi. Mycorrhiza. 2000, 10: 101-106. 10.1007/s005720000064.View ArticleGoogle Scholar
- Solaiman MZ, Senoo K, Kawaguchi M, Imaizumi-Anraku H, Akao S, Tanaka A, Obata H: Characterization of mycorrhizas formed by Glomus sp. on roots of hypernodulating mutants of Lotus japonicus. J Plant Res. 2000, 113: 443-448. 10.1007/PL00013953.View ArticleGoogle Scholar
- Pearson JN, Abbott LK, Jasper DA: Mediation of competition between two colonizing VA mycorrhizal fungi by the host plant. New Phytol. 1993, 123: 93-98.View ArticleGoogle Scholar
- Vierheilig H, Garcia-Garrido JM, Wyss U, Piché Y: Systemic suppression of mycorrhizal colonization of barley roots already colonized by AM fungi. Soil Biol Biochem. 2000, 32: 589-595. 10.1016/S0038-0717(99)00155-8.View ArticleGoogle Scholar
- Catford J-G, Staehelin C, Lerat S, Piche Y, Vierheilig H: Suppression of arbuscular mycorrhizal colonization and nodulation in split-root systems of alfalfa after pre-inoculation and treatment with Nod factors. J Exp Bot. 2003, 54: 1481-1487. 10.1093/jxb/erg156.View ArticlePubMedGoogle Scholar
- Kenjo T, Yamaya H, Arima Y: Shoot-synthesized nodulation-restricting substances of wild-type soybean present in two different high-performance liquid chromatography peaks of the ethanol-soluble medium-polarity fraction. Soil Sci Plant Nutr. 2010, 56: 399-406. 10.1111/j.1747-0765.2010.00467.x.View ArticleGoogle Scholar
- Lin Y-H, Ferguson BJ, Kereszt A, Gresshoff PM: Suppression of hypernodulation in soybean by a leaf-extracted, NARK- and Nod factor-dependent, low molecular mass fraction. New Phytol. 2010, 185: 1074-1086. 10.1111/j.1469-8137.2009.03163.x.View ArticlePubMedGoogle Scholar
- Magori S, Oka-Kira E, Shibata S, Umehara Y, Kouchi H, Hase Y, Tanaka A, Sato S, Tabata S, Kawaguchi M: TOO MUCH LOVE, a root regulator associated with the long-distance control of nodulation in Lotus japonicus. Mol Plant Microbe Interact. 2009, 22: 259-268. 10.1094/MPMI-22-3-0259.View ArticlePubMedGoogle Scholar
- Kinkema M, Gresshoff PM: Investigation of downstream signals of the soybean autoregulation of nodulation receptor kinase GmNARK. Mol Plant Microbe Interact. 2008, 21: 1337-1348. 10.1094/MPMI-21-10-1337.View ArticlePubMedGoogle Scholar
- Seo HS, Li J, Lee S-Y, Yu J-W, Kim K-H, Lee S-H, Lee I-J, Paek N-C: The hypernodulating nts mutation induces jasmonate synthetic pathway in soybean leaves. Mol Cells. 2007, 24: 185-193.PubMedGoogle Scholar
- van Noorden GE, Ross JJ, Reid JB, Rolfe BG, Mathesius U: Defective long-distance auxin transport regulation in the Medicago truncatula super numeric nodules mutant. Plant Physiol. 2006, 140: 1494-1506. 10.1104/pp.105.075879.View ArticlePubMedPubMed CentralGoogle Scholar
- Reid DE, Hayashi S, Lorenc M, Stiller J, Edwards D, Gresshoff PM, Ferguson BJ: Identification of systemic responses in soybean nodulation by xylem sap feeding and complete transcriptome sequencing reveal a novel component of the autoregulation pathway. Plant Biotechnol J. 2012, 10: 680-689. 10.1111/j.1467-7652.2012.00706.x.View ArticlePubMedGoogle Scholar
- Meixner C, Vegvarib G, Ludwig-Müller J, Gagnond H, Steinkellner S, Staehelin C, Gresshoff P, Vierheilig H: Two defined alleles of the LRR receptor kinase GmNARK in supernodulating soybean govern differing autoregulation of mycorrhization. Physiol Plant. 2007, 130: 261-270. 10.1111/j.1399-3054.2007.00903.x.View ArticleGoogle Scholar
- Bucher P: Weight matrix descriptions of four eukaryotic RNA polymerase II promoter elements derived from 502 unrelated promoter sequences. J Mol Biol. 1990, 212: 563-578. 10.1016/0022-2836(90)90223-9.View ArticlePubMedGoogle Scholar
- Dolfini D, Gatta R, Mantovani R: NF-Y and the transcriptional activation of CCAAT promoters. Crit Rev Biochem Mol Biol. 2012, 47: 29-49. 10.3109/10409238.2011.628970.View ArticlePubMedGoogle Scholar
- Edwards D, Murray JAH, Smith AG: Multiple genes encoding the conserved CCAAT-box transcription factor complex are expressed in Arabidopsis. Plant Physiol. 1998, 117: 1015-1022. 10.1104/pp.117.3.1015.View ArticlePubMedPubMed CentralGoogle Scholar
- Mantovani R: The molecular biology of the CCAAT-binding factor NF-Y. Gene. 1999, 239: 15-27. 10.1016/S0378-1119(99)00368-6.View ArticlePubMedGoogle Scholar
- Laloum T, De Mita S, Gamas P, Baudin M, Niebel A: CCAAT-box binding transcription factors in plants: Y so many?. Trends Plant Sci. 2012, 1001: NO of Pages 10Google Scholar
- Petroni K, Kumimoto RW, Gnesutta N, Calvenzani V, Fornari M, Tonelli C, Holt BF, Mantovani R: The promiscuous life of plant NUCLEAR FACTOR Y transcription factors. Plant Cell Online. 2012, 24: 4777-4792. 10.1105/tpc.112.105734.View ArticleGoogle Scholar
- Asamizu E, Shimoda Y, Kouchi H, Tabata S, Sato S: A positive regulatory role for LjERF1 in the nodulation process is revealed by systematic analysis of nodule-associated transcription factors of Lotus japonicus. Plant Physiol. 2008, 147: 2030-2040. 10.1104/pp.108.118141.View ArticlePubMedPubMed CentralGoogle Scholar
- Combier JP, de Billy F, Gamas P, Niebel A, Rivas S: Trans-regulation of the expression of the transcription factor MtHAP2-1 by a uORF controls root nodule development. Genes Dev. 2008, 22: 1549-1559. 10.1101/gad.461808.View ArticlePubMedPubMed CentralGoogle Scholar
- Combier J-P, Frugier F, de Billy F, Boualem A, El-Yahyaoui F, Moreau S, Vernié T, Ott T, Gamas P, Crespi M, Niebel A: MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev. 2006, 20: 3084-3088. 10.1101/gad.402806.View ArticlePubMedPubMed CentralGoogle Scholar
- El Yahyaoui F, Küster H, Ben Amor B, Hohnjec N, Pühler A, Becker A, Gouzy J, Vernie T, Gough C, Niebel A, Godiard L, Gamas P: Expression profiling in Medicago truncatula identifies more than 750 genes differentially expressed during nodulation, including many potential regulators of the symbiotic program. Plant Physiol. 2004, 136: 3159-3176. 10.1104/pp.104.043612.View ArticlePubMedPubMed CentralGoogle Scholar
- Hayashi S, Reid DE, Lorenc MT, Stiller J, Edwards D, Gresshoff PM, Ferguson BJ: Transient Nod factor-dependent gene expression in the nodulation-competent zone of soybean (Glycine max [L.] Merr.) roots. Plant Biotechnol J. 2012, 10: 995-1010. 10.1111/j.1467-7652.2012.00729.x.View ArticlePubMedGoogle Scholar
- Hogekamp C, Arndt D, Pereira PA, Becker JD, Hohnjec N, Küster H: Laser microdissection unravels cell-type-specific transcription in arbuscular mycorrhizal roots, including CAAT-box transcription factor gene expression correlating with fungal contact and spread. Plant Physiol. 2011, 157: 2023-2043. 10.1104/pp.111.186635.View ArticlePubMedPubMed CentralGoogle Scholar
- Libault M, Joshi T, Takahashi K, Hurley-Sommer A, Puricelli K, Blake S, Finger RE, Taylor CG, Xu D, Nguyen HT, Stacey G: Large-scale analysis of putative soybean regulatory gene expression identifies a Myb gene involved in soybean nodule development. Plant Physiol. 2009, 151: 1207-1220. 10.1104/pp.109.144030.View ArticlePubMedPubMed CentralGoogle Scholar
- Meschini EP, Blanco FA, Zanetti ME, Beker MP, Küster H, Pühler A, Aguilar OM: Host genes involved in nodulation preference in common bean (Phaseolus vulgaris)-Rhizobium etli symbiosis revealed by suppressive subtractive hybridization. Mol Plant Microbe Interact. 2008, 21: 459-468. 10.1094/MPMI-21-4-0459.View ArticlePubMedGoogle Scholar
- Moreau S, Verdenaud M, Ott T, Letort S, de Billy F, Niebel A, Gouzy J, de Carvalho-Niebel F, Gamas P: Transcription reprogramming during root nodule development in Medicago truncatula. PLoS ONE. 2011, 6: e16463-10.1371/journal.pone.0016463.View ArticlePubMedPubMed CentralGoogle Scholar
- Soyano T, Kouchi H, Hirota A, Hayashi M: NODULE INCEPTION directly targets NF-Y subunit genes to regulate essential processes of root nodule development in Lotus japonicus. PLoS Genet. 2013, 9: e1003352-10.1371/journal.pgen.1003352.View ArticlePubMedPubMed CentralGoogle Scholar
- Zanetti ME, Blanco FA, Beker MP, Battaglia M, Aguilar OM: A C subunit of the plant nuclear factor NF-Y required for rhizobial infection and nodule development affects partner selection in the common bean-Rhizobium etli symbiosis. Plant Cell Online. 2010, 22: 4142-4157. 10.1105/tpc.110.079137.View ArticleGoogle Scholar
- Vierheilig H: Further root colonization by arbuscular mycorrhizal fungi in already mycorrhizal plants is suppressed after a critical level of root colonization. J Plant Physiol. 2004, 161: 339-341. 10.1078/0176-1617-01097.View ArticlePubMedGoogle Scholar
- Corradi N, Kuhn G, Sanders IR: Monophyly of β-tubulin and H + -ATPase gene variants in Glomus intraradices: consequences for molecular evolutionary studies of AM fungal genes. Fungal Genet Biol. 2004, 41: 262-273. 10.1016/j.fgb.2003.11.001.View ArticlePubMedGoogle Scholar
- Floß DS, Hause B, Lange PR, Küster H, Strack D, Walter MH: Knock-down of the MEP pathway isogene 1-deoxy-d-xylulose 5-phosphate synthase 2 inhibits formation of arbuscular mycorrhiza-induced apocarotenoids, and abolishes normal expression of mycorrhiza-specific plant marker genes. Plant J. 2008, 56: 86-100. 10.1111/j.1365-313X.2008.03575.x.View ArticlePubMedGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: Expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98: 31-36. 10.1073/pnas.98.1.31.View ArticlePubMedGoogle Scholar
- Hohnjec N, Vieweg MF, Puhler A, Becker A, Kuster H: Overlaps in the transcriptional profiles of Medicago truncatula roots inoculated with two different Glomus fungi provide insights into the genetic program activated during arbuscular mycorrhiza. Plant Physiol. 2005, 137: 1283-1301. 10.1104/pp.104.056572.View ArticlePubMedPubMed CentralGoogle Scholar
- Manthey K, Krajinski F, Hohnjec N, Firnhaber C, Pühler A, Perlick AM, Küster H: Transcriptome profiling in root nodules and arbuscular mycorrhiza identifies a collection of novel genes induced during Medicago truncatula root endosymbioses. Mol Plant Microbe Interact. 2004, 17: 1063-1077. 10.1094/MPMI.2004.17.10.1063.View ArticlePubMedGoogle Scholar
- Wulf A, Manthey K, Doll J, Perlick AM, Linke B, Bekel T, Meyer F, Franken P, Kuster H, Krajinski F: Transcriptional changes in response to arbuscular mycorrhiza development in the model plant Medicago truncatula. Mol Plant Microbe Interact. 2003, 16: 306-314. 10.1094/MPMI.2003.16.4.306.View ArticlePubMedGoogle Scholar
- Guether M, Balestrini R, Hannah M, He J, Udvardi MK, Bonfante P: Genome-wide reprogramming of regulatory networks, transport, cell wall and membrane biogenesis during arbuscular mycorrhizal symbiosis in Lotus japonicus. New Phytol. 2009, 182: 200-212. 10.1111/j.1469-8137.2008.02725.x.View ArticlePubMedGoogle Scholar
- Liu J, Blaylock LA, Endre G, Cho J, Town CD, VandenBosch KA, Harrison MJ: Transcript profiling coupled with spatial expression analyses reveals genes involved in distinct developmental stages of an arbuscular mycorrhizal symbiosis. Plant Cell Online. 2003, 15: 2106-2123. 10.1105/tpc.014183.View ArticleGoogle Scholar
- Liu J, Maldonado-Mendoza I, Lopez-Meyer M, Cheung F, Town CD, Harrison MJ: Arbuscular mycorrhizal symbiosis is accompanied by local and systemic alterations in gene expression and an increase in disease resistance in the shoots. Plant J. 2007, 50: 529-544. 10.1111/j.1365-313X.2007.03069.x.View ArticlePubMedGoogle Scholar
- Güimil S, Chang H-S, Zhu T, Sesma A, Osbourn A, Roux C, Ioannidis V, Oakeley EJ, Docquier M, Descombes P, Briggs SP, Paszkowski U: Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization. Proc Natl Acad Sci USA. 2005, 102: 8066-8070. 10.1073/pnas.0502999102.View ArticlePubMedPubMed CentralGoogle Scholar
- Frickey T, Weiller G: Analyzing microarray data using CLANS. Bioinformatics. 2007, 23: 1170-1171. 10.1093/bioinformatics/btm079.View ArticlePubMedGoogle Scholar
- Terakado J, Yoneyama T, Fujihara S: Shoot-applied polyamines suppress nodule formation in soybean (Glycine max). J Plant Physiol. 2006, 163: 497-505. 10.1016/j.jplph.2005.05.007.View ArticlePubMedGoogle Scholar
- Adio AM, Casteel CL, De Vos M, Kim JH, Joshi V, Li B, Juéry C, Daron J, Kliebenstein DJ, Jander G: Biosynthesis and defensive function of N δ -acetylornithine, a jasmonate-induced Arabidopsis metabolite. Plant Cell Online. 2011, 23: 3303-3318. 10.1105/tpc.111.088989.View ArticleGoogle Scholar
- Kalamaki MS, Georgios Merkouropoulos , Kanellis AK: Can ornithine accumulation modulate abiotic stress tolerance in Arabidopsis?. Plant Signal Behav. 2009, 4: 1099-1101. 10.4161/psb.4.11.9873.View ArticlePubMedPubMed CentralGoogle Scholar
- Takahashi T, Kakehi J-I: Polyamines: ubiquitous polycations with unique roles in growth and stress responses. Ann Bot. 2010, 105: 1-6. 10.1093/aob/mcp259.View ArticlePubMedGoogle Scholar
- Cordin O, Banroques J, Tanner NK, Linder P: The DEAD-box protein family of RNA helicases. Gene. 2006, 367: 17-37.View ArticlePubMedGoogle Scholar
- Gong Z, Dong C-H, Lee H, Zhu J, Xiong L, Gong D, Stevenson B, Zhub J-K: A DEAD box RNA helicase is essential for mRNA export and important for development and stress responses in Arabidopsis. Plant Cell. 2005, 17: 256-267. 10.1105/tpc.104.027557.View ArticlePubMedPubMed CentralGoogle Scholar
- Schnabel E, Journet E-P, Carvalho-Niebel F, Duc G, Frugoli J: The Medicago truncatula SUNN gene encodes a CLV1-like leucine-rich repeat receptor kinase that regulates nodule number and root length. Plant Mol Biol. 2005, 58: 809-822. 10.1007/s11103-005-8102-y.View ArticlePubMedGoogle Scholar
- Wopereis J, Pajuelo E, Dazzo FB, Jiang Q, Gresshoff PM, Bruijn FJd, Stougaard J, Szczyglowski K: Short root mutant of Lotus japonicus with a dramatically altered symbiotic phenotype. Plant J. 2000, 23: 97-114. 10.1046/j.1365-313x.2000.00799.x.View ArticlePubMedGoogle Scholar
- Konopka-Postupolska D, Clark G, Hofmann A: Structure, function and membrane interactions of plant annexins: An update. Plant Sci. 2011, 181: 230-241. 10.1016/j.plantsci.2011.05.013.View ArticlePubMedGoogle Scholar
- Talukdar T, Gorecka K, de Carvalho-Niebel F, Downie J, Cullimore J, Pikula S: Annexins-calcium- and membrane-binding proteins in the plant kingdom: potential role in nodulation and mycorrhization in Medicago truncatula. Acta Biochim Pol. 2009, 56: 199-210.PubMedGoogle Scholar
- Amiour N, Recorbet G, Robert F, Gianinazzi S, Dumas-Gaudot E: Mutations in DMI3 and SUNN modify the appressorium-responsive root proteome in arbuscular mycorrhiza. Mol Plant Microbe Interact. 2006, 19: 988-997. 10.1094/MPMI-19-0988.View ArticlePubMedGoogle Scholar
- de Carvalho-Niebel F, Timmers ACJ, Chabaud M, Defaux-Petras A, Barker DG: The Nod factor-elicited annexin MtAnn1 is preferentially localised at the nuclear periphery in symbiotically activated root tissues of Medicago truncatula. Plant J. 2002, 32: 343-352. 10.1046/j.1365-313X.2002.01429.x.View ArticlePubMedGoogle Scholar
- Benning UF, Tamot B, Guelette BS, Hoffmann-Benning S: New aspects of phloem-mediated long-distance lipid signaling in plants. Front Plant Sci. 2012, 3: 53-View ArticlePubMedPubMed CentralGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, et al: Genome sequence of the palaeopolyploid soybean. Nature. 2010, 463: 178-183. 10.1038/nature08670.View ArticlePubMedGoogle Scholar
- SoyDB. [http://planttfdb.cbi.edu.cn/index.php?sp=Gma]
- Wang Z, Libault M, Joshi T, Valliyodan B, Nguyen H, Xu D, Stacey G, Cheng J: SoyDB: a knowledge database of soybean transcription factors. BMC Plant Biol. 2010, 10: 14-10.1186/1471-2229-10-14.View ArticlePubMedPubMed CentralGoogle Scholar
- 'Gene Networks in Seed Development' Soybean IVT Array Annotation. [http://www.seedgenenetwork.net/annotate]
- Harrison M, Dewbre G, Liu J: A phosphate transporter from Medicago truncatula involved in the acquisition of phosphate released by arbuscular mycorrhizal fungi. Plant Cell. 2002, 14: 2413-2429. 10.1105/tpc.004861.View ArticlePubMedPubMed CentralGoogle Scholar
- Javot H, Penmetsa R, Terzaghi N, Cook D, Harrison M: A Medicago truncatula phosphate transporter indispensable for the arbuscular mycorrhizal symbiosis. Proc Natl Acad Sci USA. 2007, 30: 1720-1725.View ArticleGoogle Scholar
- Qin L, Zhao J, Tian J, Chen L, Sun Z, Guo Y, Lu X, Gu M, Xu G, Liao H: The high-affinity phosphate transporter GmPT5 regulates phosphate transport to nodules and nodulation in soybean. Plant Physiol. 2012, 159: 1634-1643. 10.1104/pp.112.199786.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura Y, Kobae Y, Mizuno T, Hata S: Identification and expression analysis of arbuscular mycorrhiza-inducible phosphate transporter genes of soybean. Biosci Biotechnol Biochem. 2012, 76: 309-313. 10.1271/bbb.110684.View ArticlePubMedGoogle Scholar
- Catford J, Staehelin C, Larose G, Piché Y, Vierheilig H: Systemically suppressed isoflavonoids and their stimulating effects on nodulation and mycorrhization in alfalfa split-root systems. Plant Soil. 2006, 285: 257-266. 10.1007/s11104-006-9012-8.View ArticleGoogle Scholar
- Hewitt EJ: Sand and water culture methods used in the study of plant nutrition. 1966, Farnhan Royal: Commonwealth Agricultural BureauxGoogle Scholar
- Schüßler A, Walker C: The Glomeromycota. A species list with new families and new genera. 2010, Arthur Schüßler & Christopher Walker, Gloucester. Published in December 2010 in libraries at The Royal Botanic Garden Edinburgh, The Royal Botanic Garden Kew, Botanische Staatssammlung Munich, and Oregon State UniversityGoogle Scholar
- Maier W, Peipp H, Schmidt J, Wray V, Strack D: Levels of a terpenoid glycoside (blumenin) and cell wall-bound phenolics in some cereal mycorrhizas. Plant Physiol. 1995, 109: 465-470. 10.1104/pp.109.2.465.View ArticlePubMedPubMed CentralGoogle Scholar
- Limpens E, Ramos J, Franken C, Raz V, Compaan B, Franssen H, Bisseling T, Geurts R: RNA interference in Agrobacterium rhizogenes-transformed roots of Arabidopsis and Medicago truncatula. J Exp Bot. 2004, 55: 983-992. 10.1093/jxb/erh122.View ArticlePubMedGoogle Scholar
- Kereszt A, Li D, Indrasumunar A, Nguyen CDT, Nontachaiyapoom S, Kinkema M, Gresshoff PM: Agrobacterium rhizogenes-mediated transformation of soybean to study root biology. Nat Protoc. 2007, 2: 948-952. 10.1038/nprot.2007.141.View ArticlePubMedGoogle Scholar
- Vierheilig H, Coughlan AP, Wyss U, Piché Y: Ink and vinegar, a simple staining technique for arbuscular-mycorrhizal fungi. Appl Environ Microbiol. 1998, 64: 5004-5007.PubMedPubMed CentralGoogle Scholar
- NCBI BLAST. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
- Phytozome. [http://www.phytozome.net]
- Dash S, Van Hemert J, Hong L, Wise RP, Dickerson JA: PLEXdb: gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2012, 40: D1194-D1201. 10.1093/nar/gkr938.View ArticlePubMedGoogle Scholar
- PLEXdb Experiment GM53. [http://www.plexdb.org/modules/PD_browse/experiment_browser.php?experiment=GM53]
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.View ArticlePubMedPubMed CentralGoogle Scholar
- GEO Series GSE44685. [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE44685]
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.