Nodule number regulation in legumes is controlled by a feedback loop that integrates nutrient and rhizobia symbiont status signals to regulate nodule development. Signals from the roots are perceived by shoot receptors, including a CLV1-like receptor-like kinase known as SUNN in Medicago truncatula. In the absence of functional SUNN, the autoregulation feedback loop is disrupted, resulting in hypernodulation. To elucidate early autoregulation mechanisms disrupted in SUNN mutants, we searched for genes with altered expression in the loss-of-function sunn-4 mutant and included the rdn1-2 autoregulation mutant for comparison. We identified constitutively altered expression of small groups of genes in sunn-4 roots and in sunn-4 shoots. All genes with verified roles in nodulation that were induced in wild-type roots during the establishment of nodules were also induced in sunn-4, including autoregulation genes TML2 and TML1. Only an isoflavone-7-O-methyltransferase gene was induced in response to rhizobia in wild-type roots but not induced in sunn-4. In shoot tissues of wild-type, eight rhizobia-responsive genes were identified, including a MYB family transcription factor gene that remained at a baseline level in sunn-4; three genes were induced by rhizobia in shoots of sunn-4 but not wild-type. We cataloged the temporal induction profiles of many small secreted peptide (MtSSP) genes in nodulating root tissues, encompassing members of twenty-four peptide families, including the CLE and IRON MAN families. The discovery that expression of TML2 in roots, a key factor in inhibiting nodulation in response to autoregulation signals, is also triggered in sunn-4 in the section of roots analyzed, suggests that the mechanism of TML regulation of nodulation in M. truncatula may be more complex than published models.

In response to colonization by rhizobia bacteria, legumes are able to form nitrogen-fixing nodules in their roots, allowing the plants to grow efficiently in nitrogen-depleted environments. Legumes utilize a complex, long-distance signaling pathway to regulate nodulation that involves signals in both roots and shoots. We measured the transcriptional response to treatment with rhizobia in both the shoots and roots of Medicago truncatula over a 72-h time course. To detect temporal shifts in gene expression, we developed GeneShift, a novel computational statistics and machine learning workflow that addresses the time series replicate the averaging issue for detecting gene expression pattern shifts under different conditions. We identified both known and novel genes that are regulated dynamically in both tissues during early nodulation including leginsulin, defensins, root transporters, nodulin-related, and circadian clock genes. We validated over 70% of the expression patterns that GeneShift discovered using an independent M. truncatula RNA-Seq study. GeneShift facilitated the discovery of condition-specific temporally differentially expressed genes in the symbiotic nodulation biological system. In principle, GeneShift should work for time-series gene expression profiling studies from other systems.

Legumes can establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules after root exposure to lipo-chito-oligosaccharides secreted by the bacteria. Nodule development initiates with anticlinal mitotic divisions in the pericycle and endodermal and inner cortical cells, establishing cell lineages that ultimately form each nodule compartment. We characterized these lineages by isolating and sequencing the transcriptome ofMedicago truncatulasingle nuclei derived from uninoculated roots and roots undergoing early nodule development at 24, 48, and 96 hours after inoculation. To enrich samples for cells responding to the rhizobia, we complemented the analysis of theMedicagowild-type genotype A17 with a mutant for the autoregulation of nodulation,sunn-4. Analysis of cell lineage trajectories derived from the cortex indicates that their transcriptome is initially enriched for cytokinin perception and signaling while repressing auxin accumulation. As these cells differentiate to form nodules, expression of genes related to auxin biosynthesis, transport, and signaling was enhanced, while genes involved in cytokinin degradation were activated as lineages bifurcated to form the nodule meristem and infection zones. While the contribution of auxin and cytokinin in nodule development has been recognized, this single-cell resource quantifies the expression of each of their regulators, receptors and targets as cells divide and differentiate to form each nodule compartment.

We report a public resource for examining the spatiotemporal RNA expression of 54,893 M. truncatula genes during the first 72 hours of response to rhizobial inoculation. Using a methodology that allows synchronous inoculation and growth of over 100 plants in a single media container, we harvested the same segment of each root responding to rhizobia in the initial inoculation over a time course, collected individual tissues from these segments with laser capture microdissection, and created and sequenced RNA libraries generated from these tissues. We demonstrate the utility of the resource by examining the expression patterns of a set of genes induced very early in nodule signaling, as well as two gene families (CLE peptides and nodule specific PLAT-domain proteins) and show that despite similar whole root expression patterns, there are tissue differences in expression between the genes. Using a rhizobial response data set generated from transcriptomics on intact root segments, we also examined differential temporal expression patterns and determined that, after nodule tissue, the epidermis and cortical cells contained the most temporally patterned genes. We circumscribed gene lists for each time and tissue examined and developed an expression pattern visualization tool. Finally, we explored transcriptomic differences between the inner cortical cells that become nodules and those that do not, confirming that the expression of ACC synthases distinguishes inner cortical cells that become nodules and provide and describe potential downstream genes involved in early nodule cell division.