- Research article
- Open Access
Role of the GRAS transcription factor ATA/RAM1 in the transcriptional reprogramming of arbuscular mycorrhiza in Petunia hybrida
BMC Genomics volume 18, Article number: 589 (2017)
Development of arbuscular mycorrhiza (AM) requires a fundamental reprogramming of root cells for symbiosis. This involves the induction of hundreds of genes in the host. A recently identified GRAS-type transcription factor in Petunia hybrida, ATA/RAM1, is required for the induction of host genes during AM, and for morphogenesis of the fungal endosymbiont. To better understand the role of RAM1 in symbiosis, we set out to identify all genes that depend on activation by RAM1 in mycorrhizal roots.
We have carried out a transcript profiling experiment by RNAseq of mycorrhizal plants vs. non-mycorrhizal controls in wild type and ram1 mutants. The results show that the expression of early genes required for AM, such as the strigolactone biosynthetic genes and the common symbiosis signalling genes, is independent of RAM1. In contrast, genes that are involved at later stages of symbiosis, for example for nutrient exchange in cortex cells, require RAM1 for induction. RAM1 itself is highly induced in mycorrhizal roots together with many other transcription factors, in particular GRAS proteins.
Since RAM1 has previously been shown to be directly activated by the common symbiosis signalling pathway through CYCLOPS, we conclude that it acts as an early transcriptional switch that induces many AM-related genes, among them genes that are essential for the development of arbuscules, such as STR, STR2, RAM2, and PT4, besides hundreds of additional RAM1-dependent genes the role of which in symbiosis remains to be explored. Taken together, these results indicate that the defect in the morphogenesis of the fungal arbuscules in ram1 mutants may be an indirect consequence of functional defects in the host, which interfere with nutrient exchange and possibly other functions on which the fungus depends.
Arbuscular mycorrhiza (AM) represents a mutualistic interaction between most land plants and fungi (Glomeromycota) that culminates in the exchange of nutrients between the two partners. For the plant, the main benefit of AM is the increased supply with essential nutrient elements such as phosphorus (P), nitrogen (N), sulfur (S), copper (Cu), and zinc (Zn) [1, 2]. The functional changes associated with AM development require fundamental reprogramming of root cells, to allow them to be colonized by AM fungi and to generate the symbiotic machinery required for long-term compatibility and nutrient exchange. The establishment of AM proceeds through several stages, i.) the exchange of diffusible signals, ii.) mutual recognition, iii.) symbiotic signaling to the nucleus in cells of the host iv.) induction of AM-related genes in the host, v.) cellular rearrangement that allows accommodation of the AM fungus, and, vi.) construction of the symbiotic interface (reviewed in [3,4,5]).
The development of AM has been characterized in considerable detail in several AM host model species, in particular Medicago truncatula, Lotus japonicus, rice (Oryza sativa), and petunia (Petunia hybrida). The first presymbiotic signal is the constitutively released root-born phytohormone strigolactone (SL), which triggers hyphal branching and induces mitochondrial metabolism in AM fungi [6,7,8,9]. SL biosynthesis has been well characterized in petunia, pea (Pisum sativum), Arabidopsis thaliana and rice , and an SL transporter has been identified in petunia, which promotes AM development . AM fungi, in turn, release chitin-derived signal molecules (myc factors) [12, 13] that are thought to be recognized in the plant host by lysin motif (LysM) receptor like kinases, in analogy to nod factor receptors in the root nodule symbiosis (RNS) of the legumes .
Upon recognition of these symbiotic signals, a signaling cascade is triggered which is required for both AM and RNS, and therefore is referred to as the common symbiosis signaling pathway (CSSP). The CSSP involves the symbiosis receptor-like kinase SYMRK at the plasma membrane , and a characteristic repetitive calcium transient in and around the nucleus (calcium spiking) . This calcium signal is perceived and integrated by a dedicated nuclear calcium- and calmodulin-dependent protein kinase (CCaMK) . CCaMK interacts with, and activates, the transcription factor (TF) CYCLOPS (CYC) , which triggers a transcriptional cascade resulting in the induction of hundreds of genes that are required for the establishment of the symbiotic interaction.
Induction of symbiotic gene expression during AM has been studied in several host model species such as M. truncatula, rice, petunia, L. japonicus, and tomato (Solanum lycopersicum) [18,19,20,21,22,23,24]. These studies revealed that a set of conserved genes is induced in mycorrhizal roots, indicative of a conserved AM-related transcriptional program in different host species . The best-studied example of such a conserved gene encodes an AM-specific phosphate transporter (PT), represented by StPT3 in potato (Solanum tuberosum) , PhPT4 in petunia , MtPT4 in M. truncatula , LjPT3 in L. japonicus , PtPT11 in Populus trichocarpa , OsPT11 in rice [30, 31], and SbPT11 in Sorghum bicolor . A large group of AM-induced or AM-specific genes encode proteases , the function of which remains poorly understood.
Since hundreds of genes have to be coordinately activated in mycorrhizal roots, TFs can be expected to play a critical role in AM-related gene induction. Two TFs have been characterized in considerable detail, CYCLOPS (CYC) , and RAM1 [34,35,36,37]. CYC directly interacts with CCaMK, and therefore may act as one of the earliest AM-related TFs . One of the functions of CYC is to activate RAM1 by direct binding to the RAM1 promoter at a conserved palindromic promoter element known as the GGCGCC box . Since AM development in cyc mutants can be restored by overexpression of RAM1 , it appears that RAM1 is downstream of CYC, and that CYC acts mainly or exclusively through RAM1.
While an earlier study proposed that RAM1 may act very early in presymbiotic signaling , recent reports suggest that its role may rather be late in intracellular accommodation of the fungal partner in host cells, and in the establishment of the symbiotic interface during the formation of the highly branched arbuscules [35,36,37]. Consistent with this idea, RAM1 expression is highly induced in mycorrhizal roots [34,35,36,37].
To reveal at which level of the AM interaction RAM1 acts, an RNAseq experiment was performed on mycorrhizal roots of wild type P. hybrida, and in ram1 mutants , using a nurse plant inoculation system to force fungal infection in both genotypes. We then explored by comparative transcriptome profiling the effect of the ram1 mutation on the expression of groups of genes that are known to act at distinct stages early or late during AM development. In particular, we evaluated the expression of SL-related genes, CSSP components (common SYM genes), infection-related genes, and genes related to arbuscule functioning .
Our results show that the genes involved at early stages (SL-related genes and components of the CSSP) are not affected by the ram1 mutation. In contrast, most AM-induced genes exhibited strongly reduced levels of expression in ram1 mutants. In particular genes involved in arbuscule functioning, for example PT4 and STR, are strongly dependent on RAM1. Genes involved in infection take an intermediate position, since they are partially affected by ram1 mutations.
An interesting outcome of our study is the finding that a large number of GRAS genes are induced in mycorrhizal roots. In particular, a dicot-specific clade that was recently described as the MIG clade  is highly induced in mycorrhizal roots.
Mutations in the transcription factor gene RAM1 of P. hybrida result in abortion of AM , in particular of arbuscule formation (Fig. 1a, b). Since a transcription factor in the host is unlikely to act directly on the AM fungus, we reasoned that the mutant phenotype could be explained by defects in the induction of essential symbiosis-related genes in the plant. In order to explore the global transcriptional changes associated with AM symbiosis in P. hybrida, and to identify genes of the host the expression of which depends on RAM1, an RNAseq analysis was performed with petunia wild type and ram1-2 mutants (further referred to as ram1). To achieve comparable mycorrhizal colonization levels in wild type and the AM-defective ram1 mutant , nurse plant inoculation was used, which results in wild type levels of overall colonization in mutants, although the qualitative aspects of the mutant phenotype in symbiotic development are retained [37, 39]. Mycorrhizal roots of three biological replicate plants with an average colonization level of 20% ± 5% and 15% ± 8% in wild type and ram1, respectively, were used along with three non-colonized controls. In each treatment, ca. 19 million RNAseq reads were obtained (Additional file 1).
Global gene expression pattern in wild type and ram1
The petunia genome is predicted to contain at least 32′928 genes , of which 22′954 were found to be expressed in the roots (Additional file 2), based on the criterion that they were detected by at least 10 reads in at least one of the four conditions. Based on Baggerly’s test , 7′050 genes exhibited a significant difference in gene expression between the mycorrhizal wild type roots and non-mycorrhizal controls (Additional file 2). Among these, 1321 genes were 2-fold induced and 797 were 2-fold repressed. Applying a 5-fold threshold for the expression ratio, 650 genes were induced and 91 genes were repressed, respectively.
In ram1 mutants, 3992 genes were significantly regulated based on Baggerly’s test, of which 475 were 2-fold induced and 717 2-fold repressed in mycorrhizal roots vs. controls (Additional file 2). Applying a 5-fold threshold, only 201 genes were induced and 51 repressed, respectively (Additional file 2). Of the 5-fold induced genes in wild type and ram1, 430 and 151, respectively, exhibited background levels of expression in controls roots, indicating that they are AM-specific. In order to focus on robustly induced AM-dependent genes, the further analysis was restricted to genes with a change in expression ratio of at least 5-fold. The complete set of gene expression data is available in Additional file 3, and the complete list of the 650 genes with the induction ratio >5 is provided in Additional file 4.
ram1 shows residual expression of AM-related genes
In order to address commonalities between the expression patterns in wild type and ram1, we tested how many genes were induced in both genotypes (Fig. 1c-e). Among the 650 genes induced in the wild type, 182 were also induced in ram1 (Fig. 1c). Among the 201 genes induced in ram1, only 19 were not induced in the wild type. These results suggest that the genes induced in ram1 represent essentially a subgroup of the genes induced in the wild type. In agreement with this assumption, no gene showed an inverse gene induction ratio (induced in wild type and repressed in ram1, or vice versa) (Fig. 1c). No gene was repressed both in wild type and ram1 (Fig. 1c), consistent with previous findings that gene repression is not a consistent feature in AM-related reprogramming of gene expression . Not only more genes were induced in wild type compared to ram1, but also considerably more were specifically expressed in mycorrhizal roots (straight line of points in Fig. 1 d, e; Additional file 2).
Functional categories induced in mycorrhizal roots
Among the 650 genes induced >5-fold in the wild type, the two highest GO functional categories were ‘metabolic process’ and ‘transcription’, followed by ‘protein metabolic process’, ‘transport’, ‘secondary metabolism’, and ‘signaling’ (Fig. 2). This pattern of induction is compatible with previous transcriptomic studies in various AM hosts [18,19,20,21,22,23,24]. Based on the large number of AM-induced genes (>600), it is not surprising that many AM-induced genes belong to the category ‘transcription’. Considering the residual gene induction in ram1, a variable fraction of genes were still induced in the mutant to a certain degree (Fig. 2), although in the categories “metabolic process”, “protein metabolism”, lipid metabolism’, and ‘cell wall’, very low numbers of genes were induced in ram1 (Fig. 2).
In order to appreciate the global changes in gene expression in ram1 in a more quantitative way, the induction of all AM-induced genes (>5×) was plotted for the wild type and ram1 ordered by categories as in Fig. 2 (Additional file 5). While induction ratios from 5-fold to over 6000-fold are identically counted as “induced” in Figs. 1 and 2, this representation allows better to compare the effect of ram1 mutation on AM-related gene induction on a global scale and in a quantitative way. Additional file 5 reveals the strong general reduction in gene expression of most AM-inducible genes in ram1. All categories, and most individual genes, exhibited strongly reduced expression levels in ram1 compared to the wild type, indicating that the majority of all AM-inducible genes depends, directly or indirectly, on RAM1 for induction.
GO term enrichment analysis among differentially expressed genes in wild type and ram1
In order to test which pathways are differentially regulated by AM in the wild type and ram1, we determined the fraction of differentially expressed genes in groups that belong to defined GO categories, relative to the total number of genes in the respective categories based on the complete genome sequence. This analysis of GO term enrichment revealed several pathways to be significantly overrepresented in the AM-induced genes (Additional file 6). Based on the statistical significance of over-expressed genes in mycorrhizal pathways, it became apparent that the main metabolic pathways in the two genotypes strongly differ (Additional file 6). Only few pathways that were overrepresented in the wild type were also overrepresented in the mutant. In particular, the highly induced ammonium uptake pathway was not induced in ram1, whereas GO categories in protein ubiquitination and in defense and response to biotic stress remained induced in the mutants like in the wild type (Additional file 6). Similarly, the repression in the wild type of sulfate and phosphate transporter pathways, which may reflect the inhibition of the constitutive uptake pathways in mycorrhizal roots, was not observed in ram1 (Additional file 6). These results indicate that many symbiotic functions are particularly dependent on RAM1, either directly or, indirectly, due to the abnormal development of arbuscules. Thus, these results shed new light on major arbuscule-related cellular pathways, and on their dependence on RAM1.
Regulation of genes encoding transcription factors
Since transcriptional activation is a central feature of AM, and since the gene category with the second most AM-induced genes was “transcription”, we explored the entire set of genes encoding transcription factors (TFs) in the petunia genome , and we determined the fraction of the individual TF classes that are induced in mycorrhizal roots. Most of the AM-induced TFs belonged to the GRAS family of TFs (Additional file 7), where 26% of all members were AM-inducible (Fig. 3). After the GRAS TFs, the class with the second-largest fraction of AM-inducible members was the AP2/ERF-type TFs (Fig. 3; Additional file 7), consistent with previous reports . Among a total of 53 AM-induced TFs, 19 belong to the GRAS family (Additional file 7), 10 are classified as AP/ERF-type TFs, while the remaining 16 TFs are distributed over seven other classes of TFs (Fig. 3; Additional file 7). While the largest TF class in petunia is represented by MYB TFs (Additional file 7), only few of them were induced in AM (Fig. 3; Additional file 7).
GRAS transcription factors induced in AM
GRAS proteins represent a large family of TFs with a wide range of functions in development and plant-microbe interactions . In order to classify the GRAS members regulated during AM symbiosis, the phylogeny of all petunia GRAS-type TFs was established according to . Petunia GRAS proteins can be classified in 13 subfamilies (Fig. 4) . The 19 AM-inducible GRAS proteins are distributed over 6 subfamilies (Fig. 4).
An outstanding group is the Pt20/MIG subfamily, in which nine of the 11 members turned out to be induced in mycorrhizal roots (Fig. 4, Additional file 8). The MIG family was recently found in M. truncatula as a group of AM-induced genes with at least one member (MIG1) being required for arbuscule development . The 11 MIG genes of petunia are distributed over 5 contigs in the genome (indicated by similar contig names in Fig. 4), i.e. some of them are arranged in close proximity in the genome, similar to the M. truncatula genes MIG1-MIG3 that are organized in a cluster as tandem repeats . The MIG subfamily is specific to dicotyledonous plants species, but is absent from Arabidopsis (Additional file 8).
The RAM1/RAD1 subfamily is also absent from Arabidopsis but it exists both in monocot and dicot species that are competent to engage in AM (Additional file 9). This phylogenetic pattern is indicative of evolutionary selection for an AM-related function [44,45,46]. RAM1 and RAD1 belong to two closely related sister clades that each has at least one AM-induced member in P. hybrida, M. truncatula, and L. japonicus (Fig. 4; Additional file 9). A third clade with an AM-related phylogeny is the Os19 subfamily that has members in monocot and dicot species, but not Arabidopsis (Additional file 10). Additional AM-induced members where identified in the SCARECORW (SCR) clade (Additional file 11), in the SHORT ROOT (SHR) clade (Additional file 12), as well as in the SCL3 clade (Additional file 13).
Since the GRAS genes were most prominent among the AM-inducible TF genes (Fig. 3), we further explored the expression of a subset of 11 GRAS genes. Quantitative real-time PCR (qPCR) analysis revealed that 8 of them were highly induced, or exclusively expressed, in mycorrhizal roots (Additional file 14). Interestingly, the three remaining genes were significantly repressed in mycorrhizal roots, relative to the controls (Additional file 14). To test whether these GRAS genes may also have a role in aerial organs, we also determined their expression in leaves, buds (shoot tips), stems, petals stamens and carpels. These experiments showed that the AM-induced GRAS genes (RAM1, GRAS1-GRAS6) were primarily or exclusively expressed in mycorrhizal roots (Additional file 15). The AM-repressed GRAS genes (GRAS11-GRAS13), on the other hand, exhibited an ubiquitous expression patterns, including the floral organs. These results show that several GRAS genes are dedicated to AM symbiosis, as it was shown in L. japonicus .
Role of RAM1 during early stages of AM
One of the goals of this study was to identify the AM-inducible genes that depend on activation of RAM1, and from this information to interprete the defect of ram1 mutants in AM symbiosis. In the first report on RAM1, its function has been proposed to be related to early presymbiotic signaling through the action on its target RAM2 [34, 48]. Recent studies, however, suggested that RAM1 might act later by controlling fungal morphogenesis during the formation of arbuscules in cells of the root cortex [35,36,37]. Since the biological functions of TFs are mediated by their target genes, we assessed well-defined groups of AM-related genes that can serve as markers for distinct stages of the interaction .
The earliest genes with a role in AM are involved in the biosynthesis and in the release of the phytohormone strigolactone (SL), which stimulates AM fungi [10, 11, 49]. A second set of genes is required for signaling upon mutual recognition between the two partners in AM, as well as in root nodule symbiosis, and therefore they are known as CSSP genes [3,4,5]. Both these groups of genes were only marginally affected at the transcriptional level in ram1 mutants, indicating that their expression does not depend on RAM1 (Figs. 5 and 6).
Role of RAM1 during arbuscule development
Many of the AM-induced genes are thought to have a specific symbiosis-related function. Some of them appear to function in cellular infection (penetration) and in intracellular accommodation of AM fungi . Among these, VAPYRIN and EXO70i are the best characterized examples [50,51,52,53].
At later stages, the symbiotic function of arbuscules requires the activity of specific nutrient transporters such as the phosphate transporter PT4, which is localized to the periarbuscular membrane, where it mediates the transfer of inorganic phosphate delivered by the fungus into the symbiotic interface . Other proteins, such as STR and STR2 have been shown to be expressed in cells with arbuscules and to be localized to the periarbuscular membrane as well, but their function in symbiosis is unknown . Finally, for genes such as RAM2, it is clear that they are expressed in cells with arbuscules, while the function and the localization of the protein is elusive. Nevertheless, all these genes are established markers for cells with functional arbuscules, and they are essential for AM.
In order to test for the function of RAM1 at later stages of AM, we assessed the expression of later marker genes in ram1 mutants. The markers VAPYRIN and EXO70i, which are involved in arbuscule morphogenesis and branching, show a basal level of expression and are induced in mycorrhizal roots [50, 52, 53]. In ram1, they showed a similar behavior as in the wild type, but with considerably lower levels of induction (Fig. 7). A second component of the exocyst, EXO84, as well as an AM-responsive gene encoding a microtubule-associated protein (MAP gene), which might be involved in arbuscule morphogenesis [55, 56], showed very similar expression patterns as EXO70i (Fig. 7).
In contrast to this expression pattern, the expression of marker genes with late functions in arbusculated cells is essentially restricted to cells with functional arbuscules. In particular, PT4, AMT2, RAM2, STR, and the duplicated STR2A and STR2B, which are not expressed in non-mycorrhizal control roots (Fig. 8), were not expressed in mycorrhizal ram1 mutant roots (Fig. 8). Hence, the induction of these late AM-related marker genes depends on RAM1 and/or fully developed arbuscules. Taken together, these results indicate that the induction of late AM-specific genes strictly depends on RAM1, whereas earlier genes, including genes required for signaling, cellular penetration, and intracellular accommodation show little or intermediate dependence on RAM1. Based on the essential function during AM of RAM1-dependent genes such as PT4, RAM2, STR and STR2, the defect of ram1 mutants in arbuscule formation and in the establishment of a functional symbiosis can be explained by the lack of induction of these genes in ram1.
In order to obtain indepent evidence for the gene expression patterns described above, we employed qPCR. Selected representatives were chosen from the four gene groups based on their potential importance for AM symbiosis. D27 and CCD8 are involved in SL biosynthesis, SYMRK and CYCLOPS are components of the CSSP, VAPYRIN is essential for infection, whereas PT4 and RAM2 are essential for arbuscule functioning. Comparison of their expression showed that the early genes (D27, CCD8, SYMRK, CYCLOPS) were only marginally affected by the ram1 mutation, whereas PT4 and RAM2 induction in mycorrhizal roots was undetectable in ram1 mutants, as in non-colonized roots. These results are in line with our RNAseq data (Figs. 5, 6, 7 and 8), as well as with the literature on RAM1 [34,35,36,37, 47].
Establishment of AM
Endosymbioses of plants such as AM and RNS involve intracellular accommodation of fungal and bacterial mutualists, respectively, in root cells. In the case of RNS, a new organ, the nodule, is formed which is optimized for efficient symbiotic function. In order to allow root cells to be infected and to establish the cellular machinery required for new symbiosis-related functions, such as suppression of defense, intracellular accommodation of the microbe, and nutrient exchange, the host cells need to be fundamentally reprogrammed. The reprogramming of root cells during the establishment of AM has been described in several host species at the cellular level [55, 57,58,59,60], as well as at the transcriptomic level [18,19,20,21,22,23,24]. Both approaches revealed considerable commonalities of AM in different hosts, indicating a high degree of conservation. Consistent with these observations, a conserved signalling pathway, the CSSP, is required for the development of AM as well as of RNS [3,4,5].
Transcriptional networks in AM
The conserved patterns of gene induction in AM  suggests that a conserved machinery reprograms root cells for symbiosis, and indeed, the conserved function of the GRAS protein RAM1 in gene regulation and arbuscule morphogenesis in various AM-competent angiosperms indicates that RAM1 may represent a central component in this machinery [34,35,36,37, 47]. Little is known about the regulatory cascade that induces AM-related genes downstream of the CSSP. The best-characterized TFs with a role in AM symbiosis are RAM1 and CYCLOPS (CYC). Recent evidence suggests that CYC is activated by direct interaction with CCaMK, which responds to calcium spiking in the nucleus . Activated CYC targets the promoter motif GCCGGC that occurs in many AM-related genes of M. truncatula . In particular, CYC induces the transcription of RAM1, which in turn can activate downstream AM-related target genes [34,35,36,37, 47].
Role of GRAS transcription factors in AM symbiosis
Transcriptional reprogramming often relies on cascades of TFs with upstream master regulators, and downstream transcriptional networks involving multiple TFs for functionally relevant groups of target genes. CYC is likely to function as a master switch, since it is directly activated by the CSSP, and it is itself not transcriptionally induced during AM, just as the other components of the CSSP (Fig. 6). Hence, it does not require transcriptional activation for its function in AM. In contrast, RAM1 is strongly induced during AM. It is a direct target of CYC , and the fact that cyc mutants transformed with a RAM1 overexpression construct exhibit restored AM colonization  suggests that CYC acts primarily through RAM1. Therefore, RAM1 may represent a subordinate master regulator that activates the donwstream AM-related genes required for AM functioning. Indeed, previous results have shown that induction of several AM marker genes depends on RAM1 [34,35,36,37, 47].
Here, we show that GRAS proteins represent an important part of the transcriptional network required for AM, as it is the case for L. japonicus . An interesting case is the MIG family of TFs, which represent a dicot-specific GRAS subfamily with AM-induction in M. truncatula, L. japonicus, as well as in petunia (Additional file 8). The fact that MIG genes are arranged in tandem repeats both, in M. truncatula as well as in petunia, suggests that they have been subject to repeated gene duplications. Interestingly, the phylogenetic patterns between the two species suggests that the MIG duplications were species-specific and therefore more recent then the common ancestor of the two species (Additional file 8). Although DELLA proteins have been shown to be important for AM in monocots and dicots [36, 61, 62], no DELLA gene was induced in petunia (Fig. 4), indicating that their expression may be constitutively active independent of the mycorrhizal status.
RAM1 regulates late genes in AM with a role in symbiotic functioning
We have used several sets of marker genes for specific stages of AM to map the point of action of RAM1 in AM. Based on our comparative transcriptomic analysis, RAM1 has no role in the regulation of early genes such as the presymbiotic SL-related genes, and the CSSP genes (Figs. 5 and 6). This conclusion is consistent with a role of RAM1 downstream of CYC and downstream of the CSSP. A partial dependence of the infection markers VAPYRIN and EXO70i (Fig. 7) suggests that RAM1 may contribute, directly or indirectly to the induction of these genes.
Considering later genes that encode proteins with essential roles in nutrient exchange and in arbuscule functioning, and that are expressed only in cells with developing arbuscules, it appears that RAM1 directly or indirectly controls a wide array of these target genes (Additional file 3; Fig. 8). Interestingly, recent evidence suggests that some of these genes (in particular RAM2, STR, and STR2) may be related to the supply of the fungal partner with lipid-related substrates [63, 64]. Taken together, these findings and results from previous studies [34,35,36,37, 47], suggest that RAM1 acts as a central element in the transcriptional activation of arbuscule-related genes.
RAM1 encodes a GRAS transcription factor that is essential for the induction of the majority of AM-inducible genes. In particular, late arbuscule-related genes that encode essential functions in nutrient exchange are strictly dependent on RAM1. Since these genes are required for the establishment of AM, the defect of ram1 mutants in AM development can be explained by the lack of induction of these essential target genes. An interesting outcome of our study was the observation that many TFs, in particular GRAS TFs are induced in mycorrhizal roots, of which the majority belongs to phylogenetic clades that are absent from the non-mycorrhizal model plant Arabidopsis, thus indicating that they may have been under selection for a function in symbiosis.
Plant material, growth conditions, and inoculation
Petunia hybrida wild type W115 and petunia ram1-2  seeds were surface-sterilized with sodium hypochlorite (1.3%, 10 min), rinsed five times in sterile water and germinated on seedling substrate (Klasmann, http://www.klasmann-deilmann.com). After 4 weeks, plantlets were transferred to 10x10cm pots containing a sterilized mixture of 75% sand with 25% unfertilized soil and inoculated with around 30 g of soil inoculum of Rhizophagus irregularis (MUCL 43204). Plants were grown in growth chambers with a day/night cycle of 12 h (25 °C)/12 h (20 °C). Plants were fertilized weekly with a solution containing 3 mM MgSO4, 0.75 mM KNO3, 0.87 mM KCl, 0.2 mM KH2PO4, 1.52 mM Ca(NO3)2, 0.02 mM NaFeEDTA, 11 μM MnSO4, 1 μM ZnSO4, 30 μM H3BO3, 0.96 μM CuSO4, 0.03 μM (NH4)6Mo7O24, and 0.01 μM Na2MoO4.
Sample preparation and analysis
Plant roots were harvested 15 days post infection (dpi). Roots were removed from substrate under tap water and cut into ~1 cm small pieces. Two subsamples of about 100 mg were immediately frozen in liquid nitrogen, stored at −80 °C and lyophilised. One subsample of about 100 mg was taken for root colonization measurements.
Total RNA was extracted from lyophilized roots using the guanidinium thiocyanate-phenol-chloroform method  and DNAse treated on column (DNA, RNA and protein purification, Macherey-Nagel, Switzerland). Total RNA was quantified with the Qbit RNA BR Assay kit and purity was estimated using the Nanodrop (ND-1000, Witec, Switzerland). Preparation of 12 libraries and 2 × 100 bp pair-end Illumina HiSeq mRNA sequencing (RNA-Seq) was performed by Beckman Coulter Genomics (Danvers, MA, USA). RNAseq raw data were submitted to GEO (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=gzofymqunhctjcd&acc=GSE96896).
Raw reads were trimmed for quality and aligned to the Petunia axillaris v4 reference transcripts from solgenomics (https://solgenomics.net) using the suite CLC Genomics Workbench v9. For mapping, the minimum length fraction was 0.9, the minimum similarity fraction 0.9 and the maximum number of hits for a read was set to 10. The unique and total mapped number of reads for each transcript was determined, and then normalized to RPKM (Reads Per Kilobase of exon model per Million mapped reads). Intact pairs were counted as two, broken pairs as one. Fold change values were calculated by proportion-based test statistics  with a False Discovery Rate (FDR) correction for multiple testing . Genes were considered differentially expressed if the induction ratio was either 2-fold induced or repressed (−2 ≤ or ≥2) with an FDR-corrected p-value of p < 0.05, and if the difference between mycorrhizal and control tretments was at least 10 reads.
Gene set enrichment was performed with CLC Genomic Workbench by using the unconditional GOstats test of  based on the Petunia axillaris v4 annotation. This test measures the extent to which the annotation categories of features in up (≥2) or down (≤ − 2) differentially expressed gene lists are over- or under-represented relative to those of the features in the total Petunia gene repertoire, The “p-value” corresponds to the tail probability of the hyper geometric distribution. This value indicates whether functionally related gene groups defined by GO terms are significantly over- or underrepresented among genes that are differentially regulated compared to the gene set in the complete genome.
Real time PCR
Reverse transcription was performed on 1 μg of RNA with the Omniscript reverse transcription (RT) kit (Qiagen) using a mix of oligo(dT) and random primers (Promega) according to the manufacturer’s instructions. Quantitative Real-Time PCR was performed on a mic qPCR cycler (Bio Molecular Systems) with the SensiMix SYBR Hi-ROX Kit (Bioline). Relative expression was calculated using the GAPDH gene as reference gene . Results shown in Additional files 14 and 15 were obtained with RNA isolated from shoot tips (bud), leaves, pistils, stems, stamens, petals, control roots (roots-NM), and mycorrhizal roots harvested 5 weeks after inoculation with R. irregularis MUCL43204 (roots-M). Results shown in Additional file 16 were obtained using the same RNA samples as the RNAseq experiment.
For phylogenetic analysis, the GRAS amino acid sequences were aligned with ClustalW (http://www.ebi.ac.uk/Tools/msa/clustalw2/) using the following multiple alignment parameters: gap opening penalty 15, gap extension penalty 0.3, and delay divergent sequences set to 25%; and the Gonnet series was selected as the protein weight matrix. Neighbor joining trees were constructed using Poisson correction model for distance computation in MEGA4 . Bootstrap analysis was carried out with 1000 replicates. Branch lengths (drawn in the horizontal dimension only) are proportional to phylogenetic distances.
Calcium- and calmodulin-dependent protein kinase
Common symbiosis signalling pathway
GAI, RGA, SCR
Required for arbuscular Mycorrhiza1
George E. Nutrient uptake - contributions of arbuscular mycorrhizal fungi to plant mineral nutrition. In: Kapulnik Y, Douds DD, editors. Arbuscular mycorrhizas: physiology and function. Dordrecht: Kluwer Academic Publishers; 2000. p. 307–43.
Smith SE, Read DJ. Mycorrhizal Symbiosis. 3rd ed. New York: Academic Press; 2008.
Gutjahr C, Parniske M. Cell and developmental biology of arbuscular mycorrhiza symbiosis. Annu Rev Cell Dev Biol. 2013;29:593–617.
Harrison MJ. Cellular programs for arbuscular mycorrhizal symbiosis. Curr Opin Plant Biol. 2012;15(6):691–8.
Oldroyd GED. Speak, friend, and enter: signalling systems that promote beneficial symbiotic associations in plants. Nat Rev Microbiol. 2013;11(4):252–63.
Akiyama K, Matsuzaki K, Hayashi H. Plant sesquiterpenes induce hyphal branching in arbuscular mycorrhizal fungi. Nature. 2005;435(7043):824–7.
Besserer A, Puech-Pagès V, Kiefer P, Gomez-Roldan V, Jauneau A, Roy S, Portais JC, Roux C, Bécard G, Séjalon-Delmas N. Strigolactones stimulate arbuscular mycorrhizal fungi by activating mitochondria. PLoS Biol. 2006;4(7):1239–47.
Besserer A, Becard G, Jauneau A, Roux C, Sejalon-Delmas N. GR24, a synthetic analog of strigolactones, stimulates the mitosis and growth of the arbuscular mycorrhizal fungus Gigaspora rosea by boosting its energy metabolism. Plant Physiol. 2008;148(1):402–13.
Besserer A, Bécard G, Roux C, Séjalon-Delmas N. Role of mitochondria in the response of arbuscular mycorrhizal fungi to strigolactones. Plant Signal Behav. 2009;4(1):75–7.
Al-Babili S, Bouwmeester HJ. Strigolactones, a novel carotenoid-derived plant hormone. Annu Rev Plant Biol. 2015;66:161–86.
Kretzschmar T, Kohlen W, Sasse J, Borghi L, Schlegel M, Bachelier JB, Reinhardt D, Bours R, Bouwmeester HJ, Martinoia E. A petunia ABC protein controls strigolactone-dependent symbiotic signalling and branching. Nature. 2012;483(7389):341–6.
Genre A, Chabaud M, Balzergue C, Puech-Pagès V, Novero M, Rey T, Fournier J, Rochange S, Bécard G, Bonfante P, et al. Short-chain chitin oligomers from arbuscular mycorrhizal fungi trigger nuclear Ca2+ spiking in Medicago truncatula roots and their production is enhanced by strigolactone. New Phytol. 2013;198(1):179–89.
Maillet F, Poinsot V, André O, Puech-Pagès V, Haouy A, Gueunier M, Cromer L, Giraudet D, Formey D, Niebel A, et al. Fungal lipochitooligosaccharide symbiotic signals in arbuscular mycorrhiza. Nature. 2011;469:58–64.
Antolin-Llovera M, Petutsching EK, Ried MK, Lipka V, Nurnberger T, Robatzek S, Parniske M. Knowing your friends and foes - plant receptor-like kinases as initiators of symbiosis or defence. New Phytol. 2014;204(4):791–802.
Oldroyd GED, Downie JA. Nuclear calcium changes at the core of symbiosis signalling. Curr Opin Plant Biol. 2006;9(4):351–7.
Singh S, Parniske M. Activation of calcium- and calmodulin-dependent protein kinase (CCaMK), the central regulator of plant root endosymbiosis. Curr Opin Plant Biol. 2012;15(4):444–53.
Singh S, Katzer K, Lambert J, Cerri M, Parniske M. CYCLOPS, a DNA-binding transcriptional activator, orchestrates symbiotic root nodule development. Cell Host Microbe. 2014;15(2):139–52.
Breuillin F, Schramm J, Hajirezaei M, Ahkami A, Favre P, Druege U, Hause B, Bucher M, Kretzschmar T, Bossolini E, et al. Phosphate systemically inhibits development of arbuscular mycorrhiza in Petunia hybrida and represses genes involved in mycorrhizal functioning. Plant J. 2010;64:1002–17.
Guether M, Balestrini R, Hannah M, He J, Udvardi M, Bonfante P. Genome-wide reprogramming of regulatory networks, cell wall and membrane biogenesis during arbuscular-mycorrhizal symbiosis in Lotus japonicus. New Phytol. 2009;182(1):200–12.
Güimil S, Chang HS, Zhu T, Sesma A, Osbourn A, Roux C, Ionnidis V, Oakeley EJ, Docquier M, Descombes P, et al. Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization. Proc Natl Acad Sci U S A. 2005;102(22):8066–70.
Hohnjec N, Vieweg ME, Puhler A, Becker A, Küster 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(4):1283–301.
Liu JY, 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. 2003;15(9):2106–23.
Manthey K, Krajinski F, Hohnjec N, Firnhaber C, Puhler 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(10):1063–77.
Siciliano V, Genre A, Balestrini R, Cappellazzo G, de Wit P, Bonfante P. Transcriptome analysis of arbuscular mycorrhizal roots during development of the prepenetration apparatus. Plant Physiol. 2007;144(3):1455–66.
Rausch C, Daram P, Brunner S, Jansa J, Laloi M, Leggewie G, Amrhein N, Bucher M. A phosphate transporter expressed in arbuscule-containing cells in potato. Nature. 2001;414(6862):462–6.
Wegmüller S, Svistoonoff S, Reinhardt D, Stuurman J, Amrhein N, Bucher M. A transgenic dTph1 insertional mutagenesis system for forward genetics in mycorrhizal phosphate transport of petunia. Plant J. 2008;54(6):1115–27.
Harrison MJ, Dewbre GR, Liu JY. A phosphate transporter from Medicago truncatula involved in the acquisiton of phosphate released by arbuscular mycorrhizal fungi. Plant Cell. 2002;14(10):2413–29.
Maeda D, Ashida K, Iguchi K, Chechetka SA, Hijikata A, Okusako Y, Deguchi Y, Izui K, Hata S. Knockdown of an arbuscular mycorrhiza-inducible phosphate transporter gene of Lotus Japonicus suppresses mutualistic symbiosis. Plant Cell Physiol. 2006;47(7):807–17.
Loth-Pereda V, Orsini E, Courty PE, Lota F, Kohler A, Diss L, Blaudez D, Chalot M, Nehls U, Bucher M, et al. Structure and expression profile of the phosphate Pht1 transporter gene family in mycorrhizal Populus trichocarpa. Plant Physiol. 2011;156(4):2141–54.
Paszkowski U, Kroken S, Roux C, Briggs SP. Rice phosphate transporters include an evolutionarily divergent gene specifically activated in arbuscular mycorrhizal symbiosis. Proc Natl Acad Sci U S A. 2002;99(20):13324–9.
Yang SY, Gronlund M, Jakobsen I, Grotemeyer MS, Rentsch D, Miyao A, Hirochika H, Kumar CS, Sundaresan V, Salamin N, et al. Nonredundant regulation of rice arbuscular mycorrhizal symbiosis by two members of the PHOSPHATE TRANSPORTER1 gene family. Plant Cell. 2012;24(10):4236–51.
Walder F, Brule D, Koegel S, Wiemken A, Boller T, Courty PE. Plant phosphorus acquisition in a common mycorrhizal network: regulation of phosphate transporter genes of the Pht1 family in sorghum and flax. New Phytol. 2015;205(4):1632–45.
Takeda N, Kistner C, Kosuta S, Winzer T, Pitzschke A, Groth M, Sato S, Kaneko T, Tabata S, Parniske M. Proteases in plant root symbiosis. Phytochemistry. 2007;68(1):111–21.
Gobbato E, Marsh JF, Vernie T, Wang E, Maillet F, Kim J, Miller JB, Sun J, Bano SA, Ratet P, et al. A GRAS-type transcription factor with a specific function in mycorrhizal signaling. Curr Biol. 2012;22(23):2236–41.
Park H-J, Floss DS, Levesque-Tremblay V, Bravo A, Harrison MJ. Hyphal branching during arbuscule development requires Reduced Arbuscular Mycorrhiza1. Plant Physiol. 2015;169(4):2774–88.
Pimprikar P, Carbonnel S, Paries M, Katzer K, Klingl V, Bohmer MJ, Karl L, Floss DS, Harrison MJ, Parniske M, et al. A CCaMK-CYCLOPS-DELLA complex activates transcription of RAM1 to regulate arbuscule branching. Curr Biol. 2016;26:987–98.
Rich MK, Schorderet M, Bapaume L, Falquet L, Morel P, Vandenbussche M, Reinhardt D. A petunia GRAS transcription factor controls symbiotic gene expression and fungal morphogenesis in arbuscular mycorrhiza. Plant Physiol. 2015;168:788–97.
Heck C, Kuhn H, Heidt S, Walter S, Rieger N, Requena N. Symbiotic fungi control plant root cortex development through the novel GRAS transcription factor MIG1. Curr Biol. 2016;26(20):2770–8.
Sekhara Reddy DMR, Schorderet M, Feller U, Reinhardt D. A petunia mutant affected in intracellular accommodation and morphogenesis of arbuscular mycorrhizal fungi. Plant J. 2007;51:739–50.
Bombarely A, Moser M, Consortium TP. Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat Plants. 2016;2:16074–82.
Baggerly KA, Deng L, Morris JS, Aldaz CM. Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics. 2003;19(12):1477–83.
Hirsch S, Oldroyd GED. GRAS-domain transcription factors that regulate plant development. Plant Signal Behav. 2009;4(8):698–700.
Liu XY, Widmer A. Genome-wide comparative analysis of the GRAS gene family in Populus, Arabidopsis and rice. Plant Mol Biol Rep. 2014;32(6):1129–45.
Bravo A, York T, Pumplin N, Mueller LA, Harrison MJ. Genes conserved for arbuscular mycorrhizal symbiosis identified through phylogenomics. Nat Plants. 2016;2(2):15208.
Favre P, Bapaume L, Bossolini E, Delorenzi L, Falquet L, Reinhardt D. A novel bioinformatics pipeline to discover genes related to arbuscular mycorrhizal symbiosis based on their evolutionary conservation pattern among higher plants. BMC Plant Biol. 2014;14:333.
Delaux P-M, Varala K, Edger PP, Coruzzi GM, Pires JC, Ane J-M. Comparative phylogenomics uncovers the impact of symbiotic associations on host genome evolution. PLoS Genet. 2014;10(7):e1004487.
Xue L, Cui H, Buer B, Vijayakumar V, Delaux P-M, Junkermann S, Bucher M. Network of GRAS transcription factors involved in the control of arbuscule development in Lotus japonicus. Plant Physiol. 2015;167(3):854–71.
Wang ET, Schornack S, Marsh JF, Gobbato E, Schwessinger B, Eastmond P, Schultze M, Kamoun S, Oldroyd GED. A common signaling process that promotes mycorrhizal and oomycete colonization of plants. Curr Biol. 2012;22(23):2242–6.
Akiyama K, Hayashi H. Strigolactones: chemical signals for fungal symbionts and parasitic weeds in plant roots. Ann Bot. 2006;97(6):925–31.
Feddermann N, Duvvuru Muni RR, Zeier T, Stuurman J, Ercolin F, Schorderet M, Reinhardt D. The PAM1 gene of petunia, required for intracellular accommodation and morphogenesis of arbuscular mycorrhizal fungi, encodes a homologue of VAPYRIN. Plant J. 2010;64(3):470–81.
Murray JD, Duvvuru Muni R, Torres-Jerez I, Tang Y, Allen S, Andriankaja M, Li G, Laxmi A, Cheng X, Wen J, et al. Vapyrin, a gene essential for intracellular progression of arbuscular mycorrhizal symbiosis, is also essential for infection by rhizobia in the nodule symbiosis of Medicago truncatula. Plant J. 2011;65(2):244–52.
Pumplin N, Mondo SJ, Topp S, Starker CG, Gantt JS, Harrison MJ. Medicago truncatula Vapyrin is a novel protein required for arbuscular mycorrhizal symbiosis. Plant J. 2010;61(3):482–94.
Zhang XC, Pumplin N, Ivanov S, Harrison MJ. EXO70I is required for development of a sub-domain of the periarbuscular membrane during arbuscular mycorrhizal symbiosis. Curr Biol. 2015;25(16):2189–95.
Zhang Q, Blaylock LA, Harrison MJ. Two Medicago truncatula half-ABC transporters are essential for arbuscule development in arbuscular mycorrhizal symbiosis. Plant Cell. 2010;22(5):1483–97.
Genre A, Bonfante P. Actin versus tubulin configuration in arbuscule-containing cells from mycorrhizal tobacco roots. New Phytol. 1998;140(4):745–52.
Genre A, Bonfante P. Cytoskeleton-related proteins in tobacco mycorrhizal cells: gamma-tubulin and clathrin localisation. Eur J Histochem. 1999;43(2):105–11.
Bonfante-Fasolo P. Anatomy and morphology of VA mycorrhizae. In: Powell CL, Bagyaraj DJ, editors. VA Mycorrhizae, vol. 57. Boca Raton: CRC Press; 1984. p. 5–33.
Genre A, Chabaud M, Faccio A, Barker DG, Bonfante P. Prepenetration apparatus assembly precedes and predicts the colonization patterns of arbuscular mycorrhizal fungi within the root cortex of both Medicago truncatula and Daucus carota. Plant Cell. 2008;20(5):1407–20.
Genre A, Chabaud M, Timmers T, Bonfante P, Barker DG. Arbuscular mycorrhizal fungi elicit a novel intracellular apparatus in Medicago truncatula root epidermal cells before infection. Plant Cell. 2005;17(12):3489–99.
Pumplin N, Harrison MJ. Live-cell imaging reveals periarbuscular membrane domains and organelle location in Medicago truncatula roots during arbuscular mycorrhizal symbiosis. Plant Physiol. 2009;151(2):809–19.
Floss DS, Levy JG, Levesque-Tremblay V, Pumplin N, Harrison MJ. DELLA proteins regulate arbuscule formation in arbuscular mycorrhizal symbiosis. Proc Natl Acad Sci U S A. 2013;110(51):E5025–34.
Yu N, Luo DX, Zhang XW, Liu JZ, Wang WX, Jin Y, Dong WT, Liu JY, Liu H, Yang WB, et al. A DELLA protein complex controls the arbuscular mycorrhizal symbiosis in plants. Cell Res. 2014;24(1):130–3.
Bravo A, Brands M, Wewer V, Doermann P, Harrison MJ. Arbuscular mycorrhiza-specific enzymes FatM and RAM2 fine-tune lipid biosynthesis to promote development of arbuscular mycorrhiza. New Phytol. 2017;214(4):1631–45.
Rich MK, Nouri E, Courty PE, Reinhardt D. Diet of arbuscular mycorrhizal fungi - Bread & Butter? Trends Plant Sci. 2017; in press
Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate phenol chloroform extraction. Anal Biochem. 1987;162(1):156–9.
Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23(2):257–8.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.
We thank Caroline Gutjahr for helpful discussions.
This work was funded by the University of Fribourg, by a grant from the Swiss National Science Foundation to DR (Grant No. 31003A_135778/1), and by a grant from the State Secretariat for Education, Research and Innovation (SERI) of Switzerland to DR in the frame of COST action FA1206 (“STREAM”).
Availability of data and materials
The dataset supporting the conclusions of this article is available in the GEO repository at the following address: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=gzofymqunhctjcd&acc=GSE96896.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary table of RNAseq experiment. (PDF 55 kb)
Summary table of gene expression changes in wild type and ram1. Induction and repression ratios are expressed as the ratio between mycorrhizal roots and their respective non-mycorrhizal controls. (PDF 40 kb)
Table listing the complete gene data set of this study. (XLSX 4836 kb)
Table listing the genes induced at least 5-fold in the wild type and their expession in ram1. All AM-incucible genes were sorted according to their predicted functional category as in Fig. 2. Expression values (RPKM) are given for mycorrhizal wild type (wt-M_RPKM means), wild type control roots (wt-NM_RPKM means), mycorrhizal ram1 (ram1-M_RPKM means), and ram1 control roots (ram1-M_RPKM means). The respective induction ratio (AM vs. control roots) are given for the wildtype (column P) and ram1 (Column BB). (XLSX 542 kb)
Figure showing the global comparison of AM-dependent gene expression in wild type and ram1. Gene induction ratios for AM-induced genes were plotted for wild type (grey), and ram1 (black) for the AM-inducible genes listed in Additional file 4. An arrow indicates the start of the genes belonging to group 7 (Carbohydrate metabolism), which were only moderately induced (<50-fold), relative to the other groups. (PDF 134 kb)
Table showing the global comparison of GO terms induced or repressed in wild type and ram1. Genes were assigned to biochemical and cellular functions based on common GO terms. Only GO terms with significant overrepresentation (p < 0.05) among induced or repressed genes are shown. (XLSX 17 kb)
Table of AM-inducible transcription factors. AM-induced transcription factor genes were classified according to their sequence features. The corresponding induction pattern in Lotus japonicus is given for comparison . (PDF 50 kb)
Figure of phylogenetic tree of GRAS proteins in the AM-specific Pt20/MIG subfamily. AM-induced genes from P. axillaris (Peaxi), M. truncatula (Medtr), and L. japonicus (Lojap) are marked with red circles; the functionally tested MIG1 gene from M. truncatula is marked with a blue circle. The closest homologue in A. thaliana (AT) is highlighted with a red frame. Potri: Populus trichocarpa. The distance bar indicates substitutions per site. (PDF 62 kb)
Figure of phylogenetic tree of GRAS proteins in the AM-specific RAD1 and RAM1 subfamilies. AM-induced genes are marked with red circles; functionally tested homologues from P. axillaris (Peaxi) M. truncatula (Medtr) and L. japonicus (Lojap) are marked with blue circles. The closest homologue in A. thaliana (AT) is highlighted with a red frame. Potri: P. trichocarpa; Sobic: S bicolor; Bradi: Brachipodium distachyon; Os: O. sativa). The distance bar indicates substitutions per site. (PDF 61 kb)
Figure of phylogenetic tree of GRAS proteins in the AM-specific Os19 subfamily. AM-induced genes in L. japonicus (Lojap) and P. axillaris (Peaxi) are marked with red circles. The closest homologue in A. thaliana (AT) is highlighted with a red frame. Potri: P. trichocarpa; Sobic: S. bicolor; Bradi: B. distachyon; Os: O. sativa). The distance bar indicates substitutions per site. (PDF 58 kb)
Figure of phylogenetic tree of GRAS proteins in the SCARECROW subfamily. AM-induced genes from P. axillaris, (Peaxi), M. truncatula (Medtr), and L. japonicus (Lojap) are marked with red circles. The closest homologue in A. thaliana (AT) is highlighted with a red frame. Potri: P. trichocarpa; Sobic: S bicolor; Bradi: B. distachyon; Os: O. sativa). The distance bar indicates substitutions per site. (PDF 59 kb)
Figure of phylogenetic tree of GRAS proteins in the SHORT ROOT subfamily. AM-induced genes from P. axillaris (Peaxi), M. truncatula (Medtr) and L. japonicus (Lojap) are marked with red circles. The functionally characterized NSP1 gene from M. truncatula is marked with a blue circle. The closest homologues in A. thaliana (AT) are highlighted with red frames. Potri: P. trichocarpa; Sobic: S bicolor; Bradi: B. distachyon; Os: O. sativa). The distance bar indicates substitutions per site. (PDF 63 kb)
Figure of phylogenetic tree of GRAS proteins in the SCARECROW-LIKE3 subfamily. AM-induced genes from P. axillaris (Peaxi), M. truncatula (Medtr), and L. japonicus (Lojap) are marked with red circles. The closest homologue in A. thaliana (AT) is highlighted with a red frame. Potri: P. trichocarpa; Sobic: S bicolor; Bradi: B. distachyon; Os: O. sativa). The distance bar indicates substitutions per site. (PDF 62 kb)
Expression of GRAS genes in mycorrhizal roots. Expression analysis by qPCR of various GRAS transcription factor genes in mycorrhizal roots (dark grey columns) and control roots (light grey columns). In all cases, the expression was significantly different. Identities and gene names of GRAS genes can be found in Additional file 3. (PDF 56 kb)
Expression of GRAS genes in various aerial tissues. Global expression analysis by qPCR of various GRAS transcription factor genes in tissues collected from mycorrhizal roots (Root-M), control roots (Root-NM), shoot tips (buds), and various aerial organs (Leaf, Pistil, Stem, Stamen, Petal). Note logarithmic scale of y-axis. Identities and gene names of GRAS genes can be found in Additional file 3. (PDF 96 kb)
Expression of AM-related genes involved at different stages of AM interaction. Expression analysis by qPCR of D27, CCD8, SYMRK, VAPYRIN, PT4, and RAM2 in wild type (dark grey columns) and ram1 mutants (light grey columns) with the AM fungus R. irregularis (Ri) or in nonmycorrhizal controls. Note logarithmic scale of y-axis. Identities and gene names of GRAS genes can be found in Additional file 3. (PDF 40 kb)
About this article
Cite this article
Rich, M.K., Courty, P., Roux, C. et al. Role of the GRAS transcription factor ATA/RAM1 in the transcriptional reprogramming of arbuscular mycorrhiza in Petunia hybrida . BMC Genomics 18, 589 (2017). https://doi.org/10.1186/s12864-017-3988-8
- Petunia hybrida
- Arbuscular mycorrhiza
- GRAS transcription factor