Rhizobium symbiotic efficiency meets CEP signaling peptides.

C-terminally encoded peptides (CEP) signaling peptides are drivers of systemic pathways regulating nitrogen (N) acquisition in different plants, from Arabidopsis to legumes, depending on mineral N availability (e.g. nitrate) and on the whole plant N demand. Recent studies in the Medicago truncatula model legume revealed how root-produced CEP peptides control the root competence for endosymbiosis with N fixing rhizobia soil bacteria through the activity of the Compact Root Architecture 2 (CRA2) CEP receptor in shoots. Among CEP genes, MtCEP7 was shown to be tightly linked to nodulation, and the dynamic temporal regulation of its expression reflects the plant ability to maintain a different symbiotic root competence window depending on the symbiotic efficiency of the rhizobium strain, as well as to reinitiate a new window of root competence for nodulation.

The Medicago SymCEP7 hormone increases nodule number via shoots without compromising lateral root number.

Legumes acquire soil nutrients through nitrogen-fixing root nodules and lateral roots. To balance the costs and benefits of nodulation, legumes negatively control root nodule number by autoregulatory and hormonal pathways. How legumes simultaneously coordinate root nodule and lateral root development to procure nutrients remains poorly understood. In Medicago (Medicago truncatula), a subset of mature C-TERMINALLY ENCODED PEPTIDE (CEP) hormones can systemically promote nodule number, but all CEP hormones tested to date negatively regulate lateral root number. Here we showed that Medicago CEP7 produces a mature peptide, SymCEP7, that promotes nodulation from the shoot without compromising lateral root number. Rhizobial inoculation induced CEP7 in the susceptible root nodulation zone in a Nod factor-dependent manner, and, in contrast to other CEP genes, its transcription level was elevated in the ethylene signaling mutant sickle. Using mass spectrometry, fluorescence microscopy and expression analysis, we demonstrated that SymCEP7 activity requires the COMPACT ROOT ARCHITECTURE 2 receptor and activates the shoot-to-root systemic effector, miR2111. Shoot-applied SymCEP7 rapidly promoted nodule number in the pM to nM range at concentrations up to five orders of magnitude lower than effects mediated by root-applied SymCEP7. Shoot-applied SymCEP7 also promoted nodule number in White Clover (Trifolium repens) and Lotus (Lotus japonicus), which suggests that this biological function may be evolutionarily conserved. We propose that SymCEP7 acts in the Medicago shoot to counter balance the autoregulation pathways induced rapidly by rhizobia to enable nodulation without compromising lateral root growth, thus promoting the acquisition of nutrients other than nitrogen to support their growth.

The Compact Root Architecture 2 systemic pathway is required for the induction of cytokinins and of the miR399 in Medicago truncatula N-satisfied plants.

Plants use a combination of sophisticated local and systemic pathways to optimize growth depending on heterogeneous nutrient availability in the soil. Legume plants can acquire mineral nitrogen (N) either through their roots or via a symbiotic interaction with N-fixing rhizobia bacteria housed in so-called root nodules. To identify shoot-to-root systemic signals acting in Medicago truncatula plants at N-deficit or N-satiety, plants were grown in a split-root experimental design, in which either high or low N was provided to a half of the root system, allowing the analysis of systemic pathways independently of any local N response. Among the plant hormone families analyzed, the cytokinin trans-Zeatin accumulated in plants at N-satiety. Cytokinin application by petiole feeding led to an inhibition of both root growth and nodulation. In addition, an exhaustive analysis of miRNAs revealed that miR2111 accumulates systemically under N-deficit in both shoots and non-treated distant roots, whereas a miRNA related to inorganic Phosphate (Pi)-acquisition, the miR399, does so in plants grown at N-satiety. These two accumulation patterns are dependent on CRA2 (Compact Root Architecture 2), a receptor required for CEP (C-terminally Encoded Peptide) signaling. Constitutive ectopic expression of the miR399 reduced nodule numbers and root biomass depending on Pi availability, suggesting that the miR399-dependent Pi-acquisition regulatory module controlled by N-availability affects the development of the whole legume plant root system.

A CEP Peptide Receptor-like Kinase Regulates Auxin Biosynthesis and Ethylene Signaling to Coordinate Root Growth and Symbiotic Nodulation in Medicago truncatula.

Because of the high energy consumed during symbiotic nitrogen fixation, legumes must balance growth and symbiotic nodulation. Both lateral roots and nodules form on the root system and the developmental coordination of these organs according to reduced nitrogen (N) availability remains elusive. We show that the Compact Root Architecture 2 (MtCRA2) receptor-like kinase is essential to promote the initiation of early symbiotic nodulation and to inhibit root growth in response to low-N. MtCEP1 peptides can activate MtCRA2 under N-starvation conditions, leading to a repression of MtYUC2 auxin biosynthesis gene expression, and therefore of auxin root responses. Accordingly, the compact root architecture phenotype of cra2 can be mimicked by an auxin treatment or by over-expressing MtYUC2, and conversely, a treatment with YUC inhibitors or a MtYUC2 knock-out rescues the cra2 root phenotype. The MtCEP1-activated CRA2 can additionally interact with and phosphorylate the MtEIN2 ethylene signaling component at Ser643 and Ser924, preventing its cleavage and therefore repressing ethylene responses, thus locally promoting the root susceptibility to rhizobia. In agreement, the cra2 low nodulation phenotype is rescued by an ein2 mutation. Overall, by reducing auxin biosynthesis and inhibiting ethylene signaling, the MtCEP1/MtCRA2 pathway balances root and nodule development under low-N conditions.

A CEP Peptide Receptor-Like Kinase Regulates Auxin Biosynthesis and Ethylene Signaling to Coordinate Root Growth and Symbiotic Nodulation in Medicago truncatula.

Because of the large amount of energy consumed during symbiotic nitrogen fixation, legumes must balance growth and symbiotic nodulation. Both lateral roots and nodules form on the root system, and the developmental coordination of these organs under conditions of reduced nitrogen (N) availability remains elusive. We show that the Medicago truncatula COMPACT ROOT ARCHITECTURE2 (MtCRA2) receptor-like kinase is essential to promote the initiation of early symbiotic nodulation and to inhibit root growth in response to low N. C-TERMINALLY ENCODED PEPTIDE (MtCEP1) peptides can activate MtCRA2 under N-starvation conditions, leading to a repression of YUCCA2 (MtYUC2) auxin biosynthesis gene expression, and therefore of auxin root responses. Accordingly, the compact root architecture phenotype of cra2 can be mimicked by an auxin treatment or by overexpressing MtYUC2, and conversely, a treatment with YUC inhibitors or an MtYUC2 knockout rescues the cra2 root phenotype. The MtCEP1-activated CRA2 can additionally interact with and phosphorylate the MtEIN2 ethylene signaling component at Ser643 and Ser924, preventing its cleavage and thereby repressing ethylene responses, thus locally promoting the root susceptibility to rhizobia. In agreement with this interaction, the cra2 low nodulation phenotype is rescued by an ein2 mutation. Overall, by reducing auxin biosynthesis and inhibiting ethylene signaling, the MtCEP1/MtCRA2 pathway balances root and nodule development under low-N conditions.

Nitrate-induced CLE35 signaling peptides inhibit nodulation through the SUNN receptor and miR2111 repression.

Legume plants form nitrogen (N)-fixing symbiotic nodules when mineral N is limiting in soils. As N fixation is energetically costly compared to mineral N acquisition, these N sources, and in particular nitrate, inhibit nodule formation and N fixation. Here, in the model legume Medicago truncatula, we characterized a CLAVATA3-like (CLE) signaling peptide, MtCLE35, the expression of which is upregulated locally by high-N environments and relies on the Nodule Inception-Like Protein (NLP) MtNLP1. MtCLE35 inhibits nodule formation by affecting rhizobial infections, depending on the Super Numeric Nodules (MtSUNN) receptor. In addition, high N or the ectopic expression of MtCLE35 represses the expression and accumulation of the miR2111 shoot-to-root systemic effector, thus inhibiting its positive effect on nodulation. Conversely, ectopic expression of miR2111 or downregulation of MtCLE35 by RNA interference increased miR2111 accumulation independently of the N environment, and thus partially bypasses the nodulation inhibitory action of nitrate. Overall, these results demonstrate that the MtNLP1-dependent, N-induced MtCLE35 signaling peptide acts through the MtSUNN receptor and the miR2111 systemic effector to inhibit nodulation.

A dual legume-rhizobium transcriptome of symbiotic nodule senescence reveals coordinated plant and bacterial responses.

Senescence determines plant organ lifespan depending on aging and environmental cues. During the endosymbiotic interaction with rhizobia, legume plants develop a specific organ, the root nodule, which houses nitrogen (N)-fixing bacteria. Unlike earlier processes of the legume-rhizobium interaction (nodule formation, N fixation), mechanisms controlling nodule senescence remain poorly understood. To identify nodule senescence-associated genes, we performed a dual plant-bacteria RNA sequencing approach on Medicago truncatula-Sinorhizobium meliloti nodules having initiated senescence either naturally (aging) or following an environmental trigger (nitrate treatment or salt stress). The resulting data allowed the identification of hundreds of plant and bacterial genes differentially regulated during nodule senescence, thus providing an unprecedented comprehensive resource of new candidate genes associated with this process. Remarkably, several plant and bacterial genes related to the cell cycle and stress responses were regulated in senescent nodules, including the rhizobial RpoE2-dependent general stress response. Analysis of selected core nodule senescence plant genes allowed showing that MtNAC969 and MtS40, both homologous to leaf senescence-associated genes, negatively regulate the transition between N fixation and senescence. In contrast, overexpression of a gene involved in the biosynthesis of cytokinins, well-known negative regulators of leaf senescence, may promote the transition from N fixation to senescence in nodules.

A Cytokinin Signaling Type-B Response Regulator Transcription Factor Acting in Early Nodulation.

Nitrogen-fixing root nodulation in legumes challenged with nitrogen-limiting conditions requires infection of the root hairs by soil symbiotic bacteria, collectively referred to as rhizobia, and the initiation of cell divisions in the root cortex. Cytokinin hormones are critical for early nodulation to coordinate root nodule organogenesis and the progression of bacterial infections. Cytokinin signaling involves regulation of the expression of cytokinin primary response genes by type-B response regulator (RRB) transcription factors. RNA interference or mutation of MtRRB3, the RRB-encoding gene most strongly expressed in Medicago truncatula roots and nodules, significantly decreased the number of nodules formed, indicating a function of this RRB in nodulation initiation. Fewer infection events were also observed in rrb3 mutant roots associated with a reduced Nod factor induction of the Early Nodulin 11 (MtENOD11) infection marker, and of the cytokinin-regulated Nodulation Signaling Pathway 2 (Mt NSP2) gene. Rhizobial infections correlate with an expansion of the nuclear area, suggesting the activation of endoreduplication cycles linked to the cytokinin-regulated Cell Cycle Switch 52A (Mt CCS52A) gene. Although no significant difference in nucleus size and endoreduplication were detected in rhizobia-infected rrb3 mutant roots, expression of the MtCCS52A endoreduplication marker was reduced. As the MtRRB3 expression pattern overlaps with those of MtNSP2 and MtCCS52A in roots and nodule primordia, chromatin immunoprecipitation-quantitative PCR and protoplast trans-activation assays were used to show that MtRRB3 can interact with and trans-activate MtNSP2 and MtCCS52A promoters. Overall, we highlight that the MtRRB3 cytokinin signaling transcription factor coordinates the expression of key early nodulation genes.

A linkage disequilibrium-based statistical test for Genome-Wide Epistatic Selection Scans in structured populations.

The quest for signatures of selection using single nucleotide polymorphism (SNP) data has proven efficient to uncover genes involved in conserved and/or adaptive molecular functions, but none of the statistical methods were designed to identify interacting alleles as targets of selective processes. Here, we propose a statistical test aimed at detecting epistatic selection, based on a linkage disequilibrium (LD) measure accounting for population structure and heterogeneous relatedness between individuals. SNP-based ([Formula: see text]) and window-based ([Formula: see text]) statistics fit a Student distribution, allowing to test the significance of correlation coefficients. As a proof of concept, we use SNP data from the Medicago truncatula symbiotic legume plant and uncover a previously unknown gene coadaptation between the MtSUNN (Super Numeric Nodule) receptor and the MtCLE02 (CLAVATA3-Like) signaling peptide. We also provide experimental evidence supporting a MtSUNN-dependent negative role of MtCLE02 in symbiotic root nodulation. Using human HGDP-CEPH SNP data, our new statistical test uncovers strong LD between SLC24A5 (skin pigmentation) and EDAR (hairs, teeth, sweat glands development) world-wide, which persists after correction for population structure and relatedness in Central South Asian populations. This result suggests that epistatic selection or coselection could have contributed to the phenotypic make-up in some human populations. Applying this approach to genome-wide SNP data will facilitate the identification of coadapted gene networks in model or non-model organisms.

CEP receptor signalling controls root system architecture in Arabidopsis and Medicago.

Root system architecture (RSA) influences the effectiveness of resources acquisition from soils but the genetic networks that control RSA remain largely unclear. We used rhizoboxes, X-ray computed tomography, grafting, auxin transport measurements and hormone quantification to demonstrate that Arabidopsis and Medicago CEP (C-TERMINALLY ENCODED PEPTIDE)-CEP RECEPTOR signalling controls RSA, the gravitropic set-point angle (GSA) of lateral roots (LRs), auxin levels and auxin transport. We showed that soil-grown Arabidopsis and Medicago CEP receptor mutants have a narrower RSA, which results from a steeper LR GSA. Grafting showed that CEPR1 in the shoot controls GSA. CEP receptor mutants exhibited an increase in rootward auxin transport and elevated shoot auxin levels. Consistently, the application of auxin to wild-type shoots induced a steeper GSA and auxin transport inhibitors counteracted the CEP receptor mutant's steep GSA phenotype. Concordantly, CEP peptides increased GSA and inhibited rootward auxin transport in wild-type but not in CEP receptor mutants. The results indicated that CEP-CEP receptor-dependent signalling outputs in Arabidopsis and Medicago control overall RSA, LR GSA, shoot auxin levels and rootward auxin transport. We propose that manipulating CEP signalling strength or CEP receptor downstream targets may provide means to alter RSA.

Independent Regulation of Symbiotic Nodulation by the SUNN Negative and CRA2 Positive Systemic Pathways.

Plant systemic signaling pathways allow the integration and coordination of shoot and root organ metabolism and development at the whole-plant level depending on nutrient availability. In legumes, two systemic pathways have been reported in the Medicago truncatula model to regulate root nitrogen-fixing symbiotic nodulation. Both pathways involve leucine-rich repeat receptor-like kinases acting in shoots and proposed to perceive signaling peptides produced in roots depending on soil nutrient availability. In this study, we characterized in the M. truncatula Jemalong A17 genotype a mutant allelic series affecting the Compact Root Architecture2 (CRA2) receptor. These analyses revealed that this pathway acts systemically from shoots to positively regulate nodulation and is required for the activity of carboxyl-terminally encoded peptides (CEPs). In addition, we generated a double mutant to test genetic interactions of the CRA2 systemic pathway with the CLAVATA3/EMBRYO SURROUNDING REGION peptide (CLE)/Super Numeric Nodule (SUNN) receptor systemic pathway negatively regulating nodule number from shoots, which revealed an intermediate nodule number phenotype close to the wild type. Finally, we showed that the nitrate inhibition of nodule numbers was observed in cra2 mutants but not in sunn and cra2 sunn mutants. Overall, these results suggest that CEP/CRA2 and CLE/SUNN systemic pathways act independently from shoots to regulate nodule numbers.

Diversification of cytokinin phosphotransfer signaling genes in Medicago truncatula and other legume genomes.

BACKGROUND: Legumes can establish on nitrogen-deprived soils a symbiotic interaction with Rhizobia bacteria, leading to the formation of nitrogen-fixing root nodules. Cytokinin phytohormones are critical for triggering root cortical cell divisions at the onset of nodule initiation. Cytokinin signaling is based on a Two-Component System (TCS) phosphorelay cascade, involving successively Cytokinin-binding Histidine Kinase receptors, phosphorelay proteins shuttling between the cytoplasm and the nucleus, and Type-B Response Regulator (RRB) transcription factors activating the expression of cytokinin primary response genes. Among those, Type-A Response Regulators (RRA) exert a negative feedback on the TCS signaling. To determine whether the legume plant nodulation capacity is linked to specific features of TCS proteins, a genome-wide identification was performed in six legume genomes (Cajanus cajan, pigeonpea; Cicer arietinum, chickpea; Glycine max, soybean; Phaseolus vulgaris, common bean; Lotus japonicus; Medicago truncatula). The diversity of legume TCS proteins was compared to the one found in two non-nodulating species, Arabidopsis thaliana and Vitis vinifera, which are references for functional analyses of TCS components and phylogenetic analyses, respectively. RESULTS: A striking expansion of non-canonical RRBs was identified, notably leading to the emergence of proteins where the conserved phosphor-accepting aspartate residue is replaced by a glutamate or an asparagine. M. truncatula genome-wide expression datasets additionally revealed that only a limited subset of cytokinin-related TCS genes is highly expressed in different organs, namely MtCHK1/MtCRE1, MtHPT1, and MtRRB3, suggesting that this "core" module potentially acts in most plant organs including nodules. CONCLUSIONS: Further functional analyses are required to determine the relevance of these numerous non-canonical TCS RRBs in symbiotic nodulation, as well as of canonical MtHPT1 and MtRRB3 core signaling elements.

Unraveling new molecular players involved in the autoregulation of nodulation in Medicago truncatula.

The number of legume root nodules resulting from a symbiosis with rhizobia is tightly controlled by the plant. Certain members of the CLAVATA3/Embryo Surrounding Region (CLE) peptide family, specifically MtCLE12 and MtCLE13 in Medicago truncatula, act in the systemic autoregulation of nodulation (AON) pathway that negatively regulates the number of nodules. Little is known about the molecular pathways that operate downstream of the AON-related CLE peptides. Here, by means of a transcriptome analysis, we show that roots ectopically expressing MtCLE13 deregulate only a limited number of genes, including three down-regulated genes encoding lysin motif receptor-like kinases (LysM-RLKs), among which are the nodulation factor (NF) receptor NF Perception gene (NFP) and two up-regulated genes, MtTML1 and MtTML2, encoding Too Much Love (TML)-related Kelch-repeat containing F-box proteins. The observed deregulation was specific for the ectopic expression of nodulation-related MtCLE genes and depended on the Super Numeric Nodules (SUNN) AON RLK. Moreover, overexpression and silencing of these two MtTML genes demonstrated that they play a role in the negative regulation of nodule numbers. Hence, the identified MtTML genes are the functional counterpart of the Lotus japonicus TML gene shown to be central in the AON pathway. Additionally, we propose that the down-regulation of a subset of LysM-RLK-encoding genes, among which is NFP, might contribute to the restriction of further nodulation once the first nodules have been formed.

DELLA1-Mediated Gibberellin Signaling Regulates Cytokinin-Dependent Symbiotic Nodulation.

In legume plants, low-nitrogen soils promote symbiotic interactions with rhizobial bacteria, leading to the formation of nitrogen-fixing root nodules. Among critical signals regulating this developmental process are bacterial Nod Factors (NFs) and several plant hormones, including cytokinins (CKs) and gibberellins (GAs). Here, we show in Medicago truncatula that GA signaling mediated by DELLA1 decreases the amount of bioactive CKs in roots and negatively impacts the Cytokinin Response1 (CRE1)-dependent NF activation of a subset of CK-signaling genes as well as of the CK-regulated Nodulation Signaling Pathway2 and Ethylene Response Factor Required for Nodulation1 early nodulation genes. Consistently, a dominant-active DELLA1 protein can partially rescue the reduced nodulation of the cre1 mutant and triggers the formation of nodule-like structures when expressed in the root cortex or in the root epidermis. This suggests a model where the DELLA1-mediated GA signaling interplays with the CRE1-dependent CK pathway to regulate early nodulation in response to both NF and CK signals critical for this symbiotic interaction.

Flavonoids and Auxin Transport Inhibitors Rescue Symbiotic Nodulation in the Medicago truncatula Cytokinin Perception Mutant cre1.

Initiation of symbiotic nodules in legumes requires cytokinin signaling, but its mechanism of action is largely unknown. Here, we tested whether the failure to initiate nodules in the Medicago truncatula cytokinin perception mutant cre1 (cytokinin response1) is due to its altered ability to regulate auxin transport, auxin accumulation, and induction of flavonoids. We found that in the cre1 mutant, symbiotic rhizobia cannot locally alter acro- and basipetal auxin transport during nodule initiation and that these mutants show reduced auxin (indole-3-acetic acid) accumulation and auxin responses compared with the wild type. Quantification of flavonoids, which can act as endogenous auxin transport inhibitors, showed a deficiency in the induction of free naringenin, isoliquiritigenin, quercetin, and hesperetin in cre1 roots compared with wild-type roots 24 h after inoculation with rhizobia. Coinoculation of roots with rhizobia and the flavonoids naringenin, isoliquiritigenin, and kaempferol, or with the synthetic auxin transport inhibitor 2,3,5,-triiodobenzoic acid, rescued nodulation efficiency in cre1 mutants and allowed auxin transport control in response to rhizobia. Our results suggest that CRE1-dependent cytokinin signaling leads to nodule initiation through the regulation of flavonoid accumulation required for local alteration of polar auxin transport and subsequent auxin accumulation in cortical cells during the early stages of nodulation.

KNAT3/4/5-like class 2 KNOX transcription factors are involved in Medicago truncatula symbiotic nodule organ development.

We investigated the role of KNOX genes in legume root nodule organogenesis. Class 1 KNOX homeodomain transcription factors (TFs) are involved in plant shoot development and leaf shape diversity. Class 2 KNOX genes are less characterized, even though an antagonistic function relative to class 1 KNOXs was recently proposed. In silico expression data and further experimental validation identified in the Medicago truncatula model legume three class 2 KNOX genes, belonging to the KNAT3/4/5-like subclass (Mt KNAT3/4/5-like), as expressed during nodulation from early stages. RNA interference (RNAi)-mediated silencing and overexpression studies were used to unravel a function for KNOX TFs in nodule development. Mt KNAT3/4/5-like genes encoded four highly homologous proteins showing overlapping expression patterns during nodule organogenesis, suggesting functional redundancy. Simultaneous reduction of Mt KNAT3/4/5-like genes indeed led to an increased formation of fused nodule organs, and decreased the expression of the MtEFD (Ethylene response Factor required for nodule Differentiation) TF and its direct target MtRR4, a cytokinin response gene. Class 2 KNOX TFs therefore regulate legume nodule development, potentially through the MtEFD/MtRR4 cytokinin-related regulatory module, and may control nodule organ boundaries and shape like class 2 KNOX function in leaf development.

Different Pathways Act Downstream of the CEP Peptide Receptor CRA2 to Regulate Lateral Root and Nodule Development.

C-TERMINALLY ENCODED PEPTIDEs (CEPs) control root system architecture in a non-cell-autonomous manner. In Medicago truncatula, MtCEP1 affects root development by increasing nodule formation and inhibiting lateral root emergence by unknown pathways. Here, we show that the MtCEP1 peptide-dependent increase in nodulation requires the symbiotic signaling pathway and ETHYLENE INSENSITIVE2 (EIN2)/SICKLE (SKL), but acts independently of SUPER NUMERIC NODULES. MtCEP1-dependent inhibition of lateral root development acts through an EIN2-independent mechanism. MtCEP1 increases nodulation by promoting rhizobial infections, the developmental competency of roots for nodulation, the formation of fused nodules, and an increase in frequency of nodule development that initiates at proto-phloem poles. These phenotypes are similar to those of the ein2/skl mutant and support that MtCEP1 modulates EIN2-dependent symbiotic responses. Accordingly, MtCEP1 counteracts the reduction in nodulation induced by increasing ethylene precursor concentrations, and an ethylene synthesis inhibitor treatment antagonizes MtCEP1 root phenotypes. MtCEP1 also inhibits the development of EIN2-dependent pseudonodule formation. Finally, mutants affecting the COMPACT ROOT ARCHITECTURE2 (CRA2) receptor, which is closely related to the Arabidopsis CEP Receptor1, are unresponsive to MtCEP1 effects on lateral root and nodule formation, suggesting that CRA2 is a CEP peptide receptor mediating both organogenesis programs. In addition, an ethylene inhibitor treatment counteracts the cra2 nodulation phenotype. These results indicate that MtCEP1 and its likely receptor, CRA2, mediate nodulation and lateral root development through different pathways.

Opposing Control by Transcription Factors MYB61 and MYB3 Increases Freezing Tolerance by Relieving C-Repeat Binding Factor Suppression.

Cold acclimation is an important process by which plants respond to low temperature and enhance their winter hardiness. C-REPEAT BINDING FACTOR1 (CBF1), CBF2, and CBF3 genes were shown previously to participate in cold acclimation in Medicago truncatula In addition, MtCBF4 is transcriptionally induced by salt, drought, and cold stresses. We show here that MtCBF4, shown previously to enhance drought and salt tolerance, also positively regulates cold acclimation and freezing tolerance. To identify molecular factors acting upstream and downstream of the MtCBF4 transcription factor (TF) in cold responses, we first identified genes that are differentially regulated upon MtCBF4 overexpression using RNAseq Digital Gene Expression Profiling. Among these, we showed that MtCBF4 directly activates the transcription of the COLD ACCLIMATION SPECIFIC15 (MtCAS15) gene. To gain insights into how MtCBF4 is transcriptionally regulated in response to cold, an R2R3-MYB TF, MtMYB3, was identified based on a yeast one-hybrid screen as binding directly to MYB cis-elements in the MtCBF4 promoter, leading to the inhibition of MtCBF4 expression. In addition, another MYB TF, MtMYB61, identified as an interactor of MtMYB3, can relieve the inhibitory effect of MtMYB3 on MtCBF4 transcription. This study, therefore, supports a model describing how MtCBF4 is regulated by antagonistic MtMYB3/MtMYB61 TFs, leading to the up-regulation of downstream targets such as MtCAS15 acting in cold acclimation in M. truncatula.