Synergistic effects of plant genotype and soil microbiome on growth in Lotus japonicus.

The biological interactions between plants and their root microbiomes are essential for plant growth, and even though plant genotype (G), soil microbiome (M), and growth conditions (environment; E) are the core factors shaping root microbiome, their relationships remain unclear. In this study, we investigated the effects of G, M, and E and their interactions on the Lotus root microbiome and plant growth using an in vitro cross-inoculation approach, which reconstructed the interactions between nine Lotus accessions and four soil microbiomes under two different environmental conditions. Results suggested that a large proportion of the root microbiome composition is determined by M and E, while G-related (G, G × M, and G × E) effects were significant but small. In contrast, the interaction between G and M had a more pronounced effect on plant shoot growth than M alone. Our findings also indicated that most microbiome variations controlled by M have little effect on plant phenotypes, whereas G × M interactions have more significant effects. Plant genotype-dependent interactions with soil microbes warrant more attention to optimize crop yield and resilience.

Differential amplification of the subtelomeric satellite DNA JcSAT1 in the genus Jatropha L. (Euphorbiaceae).

Satellite DNAs (satDNAs) are highly repetitive sequences that occur in virtually all eukaryotic genomes and can undergo rapid copy number and nucleotide sequence variation among relatives. After chromosomal mapping of the satDNA JcSAT1, it was found a large accumulation at subtelomeres of Jatropha curcas (subgenus Curcas), but an absence of these monomers in J. integerrima (subgenus Jatropha). This fact suggests a dynamic scenario for this satellite repeat in Jatropha genomes. Here, we used a multitasking approach (sequence analysis, DNA blotting and chromosomal mapping) to investigate the molecular organization and chromosomal abundance and distribution of JcSAT1 in a broader group of species from the subgenus Jatropha (J. gossypiifolia, J. mollissima, J. podagrica, and J. multifida) in addition to J. curcas, with the aiming of understanding the evolution of this satDNA. Based on the analysis of BAC clone sequences of J. curcas, a large array (~ 30 kb) of 80 homogeneous monomers of JcSAT1 was identified in BAC 23J11. The monomer size was conserved (~ 358 bp) and contained a telomeric motif at the 5' end. PCR amplification coupled with a Southern blot revealed the presence of JcSAT1-like sequences in all species examined. However, a large set of genome copies was identified only in J. curcas, where a ladder-like pattern with multimers of different sizes was observed. In situ hybridization of BAC 23J11 confirmed the subtelomeric pattern for J. curcas, but showed no signals on chromosomes of species from the subgenus Jatropha. Our data indicate that JcSAT1 is a highly homogeneous satDNA that originated from a region near the telomeres and spread throughout the chromosomal subtermini, possibly due to frequent ectopic recombination between these regions. The abundance of JcSAT1 in the genome of J. curcas suggests that an amplification event occurred either at the base of the subgenus Curcas or at least in this species, although the repeat is shared by all species of the genus studied so far.

Bradyrhizobium ottawaense efficiently reduces nitrous oxide through high nosZ gene expression.

N2O is an important greenhouse gas influencing global warming, and agricultural land is the predominant (anthropogenic) source of N2O emissions. Here, we report the high N2O-reducing activity of Bradyrhizobium ottawaense, suggesting the potential for efficiently mitigating N2O emission from agricultural lands. Among the 15 B. ottawaense isolates examined, the N2O-reducing activities of most (13) strains were approximately five-fold higher than that of Bradyrhizobium diazoefficiens USDA110T under anaerobic conditions. This robust N2O-reducing activity of B. ottawaense was confirmed by N2O reductase (NosZ) protein levels and by mitigation of N2O emitted by nodule decomposition in laboratory system. While the NosZ of B. ottawaense and B. diazoefficiens showed high homology, nosZ gene expression in B. ottawaense was over 150-fold higher than that in B. diazoefficiens USDA110T, suggesting the high N2O-reducing activity of B. ottawaense is achieved by high nos expression. Furthermore, we examined the nos operon transcription start sites and found that, unlike B. diazoefficiens, B. ottawaense has two transcription start sites under N2O-respiring conditions, which may contribute to the high nosZ expression. Our study indicates the potential of B. ottawaense for effective N2O reduction and unique regulation of nos gene expression towards the high performance of N2O mitigation in the soil.

Self-grafting-induced epigenetic changes leading to drought stress tolerance in tomato plants.

Grafting is widely used as a method to increase stress tolerance in good fruiting lines of Solanaceae plants. However, little is known about how grafting, affects epigenetic modifications and leads to stress tolerance, especially within the same line. Here, we studied the effects of self-grafting in tomato plants on histone and DNA modifications and changes in gene expression related to drought stress. We found that at the three-leaf stage, 1 week after self-grafting, histone H3 K4 trimethylation and K27 trimethylation changes were observed in more than 500 genes each, and DNA methylation changes in more than 5,000 gene regions at the shoot apex compared to the non-grafted control. In addition, two weeks after the epigenomic changes, global expression changes continued to be observed at the shoot apex in several genes related to the metabolic process of nitrogen compounds, responses to stimulus, chromosome organization, cell cycle-related genes, and regulation of hormone levels. Finally, these grafted seedlings acquired remarkable drought tolerance, suggesting that epigenomic modifications during the wound-healing process mitigate stress tolerance in tomato plants.

FER and LecRK show haplotype-dependent cold-responsiveness and mediate freezing tolerance in Lotus japonicus.

Many plant species have succeeded in colonizing a wide range of diverse climates through local adaptation, but the underlying molecular genetics remain obscure. We previously found that winter survival was a direct target of selection during colonization of Japan by the perennial legume Lotus japonicus and identified associated candidate genes. Here, we show that two of these, FERONIA-receptor like kinase (LjFER) and a S-receptor-like kinase gene (LjLecRK), are required for non-acclimated freezing tolerance and show haplotype-dependent cold-responsive expression. Our work suggests that recruiting a conserved growth regulator gene, FER, and a receptor-like kinase gene, LecRK, into the set of cold-responsive genes has contributed to freezing tolerance and local climate adaptation in L. japonicus, offering functional genetic insight into perennial herb evolution.

Construction of prediction models for growth traits of soybean cultivars based on phenotyping in diverse genotype and environment combinations.

As soybean cultivars are adapted to a relatively narrow range of latitude, the effects of climate changes are estimated to be severe. To address this issue, it is important to improve our understanding of the effects of climate change by applying the simulation model including both genetic and environmental factors with their interactions (G×E). To achieve this goal, we conducted the field experiments for soybean core collections using multiple sowing times in multi-latitudinal fields. Sowing time shifts altered the flowering time (FT) and growth phenotypes, and resulted in increasing the combinations of genotypes and environments. Genome-wide association studies for the obtained phenotypes revealed the effects of field and sowing time to the significance of detected alleles, indicating the presence of G×E. By using accumulated phenotypic and environmental data in 2018 and 2019, we constructed multiple regression models for FT and growth pattern. Applicability of the constructed models was evaluated by the field experiments in 2020 including a novel field, and high correlation between the predicted and measured values was observed, suggesting the robustness of the models. The models presented here would allow us to predict the phenotype of the core collections in a given environment.

Development of a genome-wide marker design workflow for onions and its application in target amplicon sequencing-based genotyping.

Onions are one of the most widely cultivated vegetables worldwide; however, the development and utilization of molecular markers have been limited because of the large genome of this plant. We present a genome-wide marker design workflow for onions and its application in a high-throughput genotyping method based on target amplicon sequencing. The efficiency of the method was evaluated by genotyping of F2 populations. In the marker design workflow, unigene and genomic sequence data sets were constructed, and polymorphisms between parental lines were detected through transcriptome sequence analysis. The positions of polymorphisms detected in the unigenes were mapped onto the genome sequence, and primer sets were designed. In total, 480 markers covering the whole genome were selected. By genotyping an F2 population, 329 polymorphic sites were obtained from the estimated positions or the flanking sequences. However, missing or sparse marker regions were observed in the resulting genetic linkage map. We modified the markers to cover these regions by genotyping the other F2 populations. The grouping and order of markers on the linkages were similar across the genetic maps. Our marker design workflow and target amplicon sequencing are useful for genome-wide genotyping of onions owing to their reliability, cost effectiveness, and flexibility.

Comparative Transcriptome Analysis of Grafted Tomato with Drought Tolerance.

Grafting is a method used in agriculture to improve crop production and tolerance to biotic and abiotic stress. This technique is widely used in tomato, Solanum lycopersicum L.; however, the effects of grafting on changes in gene expression associated with stress tolerance in shoot apical meristem cells are still under-discovered. To clarify the effect of grafting, we performed a transcriptomic analysis between non-grafted and grafted tomatoes using the tomato variety Momotaro-scion and rootstock varieties, TD1, GS, and GF. Drought tolerance was significantly improved not only by a combination of compatible resistant rootstock TD1 but also by self-grafted compared to non-grafted lines. Next, we found the differences in gene expression between grafted and non-grafted plants before and during drought stress treatment. These altered genes are involved in the regulation of plant hormones, stress response, and cell proliferation. Furthermore, when comparing compatible (Momo/TD1 and Momo/Momo) and incompatible (Momo/GF) grafted lines, the incompatible line reduced gene expression associated with phytohormones but increased in wounding and starvation stress-response genes. These results conclude that grafting generates drought stress tolerance through several gene expression changes in the apical meristem.

Widespread and transgenerational retrotransposon activation in inter- and intraspecies recombinant inbred populations of Lotus japonicus.

Transposable elements (TEs) constitute a large proportion of genomes of multicellular eukaryotes, including flowering plants. TEs are normally maintained in a silenced state and their transpositions rarely occur. Hybridization between distant species has been regarded as a 'shock' that stimulates genome reorganization, including TE mobilization. However, whether crosses between genetically close parents that result in viable and fertile offspring can induce TE transpositions has remained unclear. Here, we investigated the activation of long terminal repeat (LTR) retrotransposons in three Lotus japonicus recombinant inbred line (RIL) populations. We found that at least six LTR retrotransposon families were activated and transposed in 78% of the RILs investigated. LORE1a, one of the transposed LTR retrotransposons, showed transgenerational epigenetic activation, indicating the long-term effects of epigenetic instability induced by hybridization. Our study highlights TE activation as an unexpectedly common event in plant reproduction.

In Vivo Evidence of Single 13C and 15N Isotope-Labeled Methanotrophic Nitrogen-Fixing Bacterial Cells in Rice Roots.

Methane-oxidizing bacteria (methanotrophs) play an ecological role in methane and nitrogen fluxes because they are capable of nitrogen fixation and methane oxidation, as indicated by genomic and cultivation-dependent studies. However, the chemical relationships between methanotrophy and diazotrophy and aerobic and anaerobic reactions, respectively, in methanotrophs remain unclear. No study has demonstrated the cooccurrence of both bioactivities in a single methanotroph bacterium in its natural environment. Here, we demonstrate that both bioactivities in type II methanotrophs occur at the single-cell level in the root tissues of paddy rice (Oryza sativa L. cv. Nipponbare). We first verified that difluoromethane, an inhibitor of methane monooxygenase, affected methane oxidation in rice roots. The results indicated that methane assimilation in the roots mostly occurred due to oxygen-dependent processes. Moreover, the results indicated that methane oxidation-dependent and methane oxidation-independent nitrogen fixation concurrently occurred in bulk root tissues. Subsequently, we performed fluorescence in situ hybridization and NanoSIMS analyses, which revealed that single cells of type II methanotrophs (involving six amplicon sequence variants) in paddy rice roots simultaneously and logarithmically fixed stable isotope gases 15N2 and 13CH4 during incubation periods of 0, 23, and 42 h, providing in vivo functional evidence of nitrogen fixation in methanotrophic cells. Furthermore, 15N enrichment in type II methanotrophs at 42 h varied among cells with an increase in 13C accumulation, suggesting that either the release of fixed nitrogen into root systems or methanotroph metabolic specialization is dependent on different microenvironmental niches in the root. IMPORTANCE Atmospheric methane concentrations have been continually increasing, causing methane to become a considerable environmental concern. Methanotrophy may be the key to regulating methane fluxes. Although research suggests that type II methanotrophs are involved in methane oxidation aerobically and nitrogen fixation anaerobically, direct evidence of simultaneous aerobic and anaerobic bioreactions of methanotrophs in situ is still lacking. In this study, a single-cell isotope analysis was performed to demonstrate these in vivo parallel functions of type II methanotrophs in the root tissues of paddy rice (Oryza sativa L. cv. Nipponbare). The results of this study indicated that methanotrophs might provide fixed nitrogen to root systems or depend on cells present in the spatially localized niche of the root tissue. Furthermore, our results suggested that single type II methanotrophic cells performed simultaneous methane oxidation and nitrogen fixation in vivo. Under natural conditions, however, nitrogen accumulation varied at the single-cell level.

Phenolic Acids Induce Nod Factor Production in Lotus japonicus-Mesorhizobium Symbiosis.

In legume-rhizobia symbiosis, partner recognition and the initiation of symbiosis processes require the mutual exchange of chemical signals. Chemicals, generally (iso)flavonoids, in the root exudates of the host plant induce the expression of nod genes in rhizobia, and, thus, are called nod gene inducers. The expression of nod genes leads to the production of lipochitooligosaccharides (LCOs) called Nod factors. Natural nod gene inducer(s) in Lotus japonicus-Mesorhizobium symbiosis remain unknown. Therefore, we developed an LCO detection method based on ultra-high-performance liquid chromatography-tandem-quadrupole mass spectrometry (UPLC-TQMS) to identify these inducers and used it herein to screen 40 phenolic compounds and aldonic acids for their ability to induce LCOs in Mesorhizobium japonicum MAFF303099. We identified five phenolic acids with LCO-inducing activities, including p-coumaric, caffeic, and ferulic acids. The induced LCOs caused root hair deformation, and nodule numbers in L. japonicus inoculated with M. japonicum were increased by these phenolic acids. The three phenolic acids listed above induced the expression of the nodA, nodB, and ttsI genes in a strain harboring a multicopy plasmid encoding NodD1, but not that encoding NodD2. The presence of p-coumaric and ferulic acids in the root exudates of L. japonicus was confirmed by UPLC-TQMS, and the induction of ttsI::lacZ in the strain harboring the nodD1 plasmid was detected in the rhizosphere of L. japonicus. Based on these results, we propose that phenolic acids are a novel type of nod gene inducer in L. japonicus-Mesorhizobium symbiosis.

Evolution of rhizobial symbiosis islands through insertion sequence-mediated deletion and duplication.

Symbiosis between organisms influences their evolution via adaptive changes in genome architectures. Immunity of soybean carrying the Rj2 allele is triggered by NopP (type III secretion system [T3SS]-dependent effector), encoded by symbiosis island A (SymA) in B. diazoefficiens USDA122. This immunity was overcome by many mutants with large SymA deletions that encompassed T3SS (rhc) and N2 fixation (nif) genes and were bounded by insertion sequence (IS) copies in direct orientation, indicating homologous recombination between ISs. Similar deletion events were observed in B. diazoefficiens USDA110 and B. japonicum J5. When we cultured a USDA122 strain with a marker gene sacB inserted into the rhc gene cluster, most sucrose-resistant mutants had deletions in nif/rhc gene clusters, similar to the mutants above. Some deletion mutants were unique to the sacB system and showed lower competitive nodulation capability, indicating that IS-mediated deletions occurred during free-living growth and the host plants selected the mutants. Among 63 natural bradyrhizobial isolates, 2 possessed long duplications (261-357 kb) harboring nif/rhc gene clusters between IS copies in direct orientation via homologous recombination. Therefore, the structures of symbiosis islands are in a state of flux via IS-mediated duplications and deletions during rhizobial saprophytic growth, and host plants select mutualistic variants from the resultant pools of rhizobial populations. Our results demonstrate that homologous recombination between direct IS copies provides a natural mechanism generating deletions and duplications on symbiosis islands.

Insights from the first genome assembly of Onion (Allium cepa).

Onion is an important vegetable crop with an estimated genome size of 16 Gb. We describe the de novo assembly and ab initio annotation of the genome of a doubled haploid onion line DHCU066619, which resulted in a final assembly of 14.9 Gb with an N50 of 464 Kb. Of this, 2.4 Gb was ordered into eight pseudomolecules using four genetic linkage maps. The remainder of the genome is available in 89.6 K scaffolds. Only 72.4% of the genome could be identified as repetitive sequences and consist, to a large extent, of (retro) transposons. In addition, an estimated 20% of the putative (retro) transposons had accumulated a large number of mutations, hampering their identification, but facilitating their assembly. These elements are probably already quite old. The ab initio gene prediction indicated 540,925 putative gene models, which is far more than expected, possibly due to the presence of pseudogenes. Of these models, 47,066 showed RNASeq support. No gene rich regions were found, genes are uniformly distributed over the genome. Analysis of synteny with Allium sativum (garlic) showed collinearity but also major rearrangements between both species. This assembly is the first high-quality genome sequence available for the study of onion and will be a valuable resource for further research.

A Dual-Color Tyr-FISH Method for Visualizing Genes/Markers on Plant Chromosomes to Create Integrated Genetic and Cytogenetic Maps.

In situ imaging of molecular markers on a physical chromosome is an indispensable tool for refining genetic maps and validation genome assembly at the chromosomal level. Despite the tremendous progress in genome sequencing, the plant genome assembly at the chromosome level remains a challenge. Recently developed optical and Hi-C mapping are aimed at assistance in genome assembly. For high confidence in the genome assembly at chromosome level, more independent approaches are required. The present study is aimed at refining an ultrasensitive Tyr-FISH technique and developing a reliable and simple method of in situ mapping of a short unique DNA sequences on plant chromosomes. We have carefully analyzed the critical steps of the Tyr-FISH to find out the reasons behind the flaws of this technique. The accurate visualization of markers/genes appeared to be significantly dependent on the means of chromosome slide preparation, probe design and labeling, and high stringency washing. Appropriate adjustment of these steps allowed us to detect a short DNA sequence of 1.6 Kb with a frequency of 51.6%. Based on our results, we developed a more reliable and simple protocol for dual-color Tyr-FISH visualization of unique short DNA sequences on plant chromosomes. This new protocol can allow for more accurate determination of the physical distance between markers and can be applied for faster integration of genetic and cytogenetic maps.

Mesorhizobium sp. J8 can establish symbiosis with Glycyrrhiza uralensis, increasing glycyrrhizin production.

Licorice (Glycyrrhiza uralensis) is a medicinal plant that contains glycyrrhizin (GL), which has various pharmacological activities. Because licorice is a legume, it can establish a symbiotic relationship with nitrogen-fixing rhizobial bacteria. However, the effect of this symbiosis on GL production is unknown. Rhizobia were isolated from root nodules of Glycyrrhiza glabra, and a rhizobium that can form root nodules in G. uralensis was selected. Whole-genome analysis revealed a single circular chromosome of 6.7 Mbp. This rhizobium was classified as Mesorhizobium by phylogenetic analysis and was designated Mesorhizobium sp. J8. When G. uralensis plants grown from cuttings were inoculated with J8, root nodules formed. Shoot biomass and SPAD values of inoculated plants were significantly higher than those of uninoculated controls, and the GL content of the roots was 3.2 times that of controls. Because uninoculated plants from cuttings showed slight nodule formation, we grew plants from seeds in plant boxes filled with sterilized vermiculite, inoculated half of the seedlings with J8, and grew them with or without 100 µM KNO3. The SPAD values of inoculated plants were significantly higher than those of uninoculated plants. Furthermore, the expression level of the CYP88D6 gene, which is a marker of GL synthesis, was 2.5 times higher than in inoculated plants. These results indicate that rhizobial symbiosis promotes both biomass and GL production in G. uralensis.

Construction of a high-density linkage map and graphical representation of the arrangement of transcriptome-based unigene markers on the chromosomes of onion, Allium cepa L.

BACKGROUND: Genomic information for Allium cepa L. is limited as it is heterozygous and its genome is very large. To elucidate potential SNP markers obtained by NGS, we used a complete set of A. fistulosum L.-A. cepa monosomic addition lines (MALs) and doubled haploids (DHs). These were the parental lines of an A. cepa mapping population for transcriptome-based SNP genotyping. RESULTS: We mapped the transcriptome sequence reads from a series of A. fistulosum-A. cepa MALs onto the unigene sequence of the doubled haploid shallot A. cepa Aggregatum group (DHA) and compared the MAL genotype call for parental bunching onion and shallot transcriptome mapping data. We identified SNP sites with at least four reads on 25,462 unigenes. They were anchored on eight A. cepa chromosomes. A single SNP site was identified on 3,278 unigenes and multiple SNPs were identified on 22,184 unigenes. The chromosome marker information was made public via the web database Allium TDB ( http://alliumtdb.kazusa.or.jp/ ). To apply transcriptome based genotyping approach for genetic mapping, we gathered RNA sequence data from 96 lines of a DHA × doubled haploid bulb onion A. cepa common onion group (DHC) mapping population. After selecting co-dominant SNP sites, 16,872 SNPs were identified in 5,339 unigenes. Of these, at least two SNPs with identical genotypes were found in 1,435 unigenes. We developed a linkage map using genotype information from these unigenes. All unigene markers mapped onto the eight chromosomes and graphical genotyping was conducted based on the unigene order information. Another 2,963 unigenes were allocated onto the eight chromosomes. To confirm the accuracy of this transcriptome-based genetic linkage map, conventional PCR-based markers were used for linkage analysis. All SNP - and PCR-based markers were mapped onto the expected linkage groups and no inconsistency was found among these chromosomal locations. CONCLUSIONS: Effective transcriptome analysis with unique Allium resources successfully associated numerous chromosome markers with unigene information and a high-density A. cepa linkage map. The information on these unigene markers is valuable in genome sequencing and useful trait detection in Allium.

Identification of a Flavin Monooxygenase-Like Flavonoid 8-Hydroxylase with Gossypetin Synthase Activity from Lotus japonicus.

Lotus japonicus is a model legume that accumulates 8-hydroxyflavonol derivatives, such as gossypetin (8-hydroxyquercetin) 3-O-glycoside, which confer the yellow color to its petals. An enzyme, flavonoid 8-hydroxylase (F8H; LjF8H), is assumed to be involved in the biosynthesis, but the specific gene is yet to be identified. The LjF8H cDNA was isolated as a flavin adenine dinucleotide (FAD)-binding monooxygenase-like protein using flower buds and flower-specific EST data of L. japonicus. LjF8H is a single copy gene on chromosome III consisting of six exons. The conserved FAD- and NAD(P)H-dependent oxidase motifs were found in LjF8H. Phylogenetic analysis suggested that LjF8H is a member of the flavin monooxygenase group but distinctly different from other known flavonoid oxygenases. Analysis of recombinant yeast microsome expressing LjF8H revealed that the enzyme catalyzed the 8-hydroxylation of quercetin. Other flavonoids, such as naringenin, eriodictyol, apigenin, luteolin, taxifolin and kaempferol, also acted as substrates of LjF8H. This broad substrate acceptance was unlike known F8Hs in other plants. Interestingly, flavanone and flavanonol, which have saturated C-C bond at positions 2 and 3 of the flavonoid C-ring, produced 6-hyroxylflavonoids as a by-product of the enzymatic reaction. Furthermore, LjF8H only accepted the 2S-isomer of naringenin, suggesting that the conformational state of the substrates might affect product specificity. The overexpression of LjF8H in Arabidopsis thaliana and Petunia hybrida synthesized gossypetin and 8-hydroxykaempferol, respectively, indicating that LjF8H was functional in plant cells. In conclusion, this study represents the first instance of cloning and identification of F8Hs responsible for gossypetin biosynthesis.

Rhizobia use a pathogenic-like effector to hijack leguminous nodulation signalling.

Legume plants form a root-nodule symbiosis with rhizobia. This symbiosis establishment generally relies on rhizobium-produced Nod factors (NFs) and their perception by leguminous receptors (NFRs) that trigger nodulation. However, certain rhizobia hijack leguminous nodulation signalling via their type III secretion system, which functions in pathogenic bacteria to deliver effector proteins into host cells. Here, we report that rhizobia use pathogenic-like effectors to hijack legume nodulation signalling. The rhizobial effector Bel2-5 resembles the XopD effector of the plant pathogen Xanthomonas campestris and could induce nitrogen-fixing nodules on soybean nfr mutant. The soybean root transcriptome revealed that Bel2-5 induces expression of cytokinin-related genes, which are important for nodule organogenesis and represses ethylene- and defense-related genes that are deleterious to nodulation. Remarkably, Bel2-5 introduction into a strain unable to nodulate soybean mutant affected in NF perception conferred nodulation ability. Our findings show that rhizobia employ and have customized pathogenic effectors to promote leguminous nodulation signalling.

Lotus japonicus Nuclear Factor YA1, a nodule emergence stage-specific regulator of auxin signalling.

Organogenesis of legume root nodules begins with the nodulation factor-dependent stimulation of compatible root cells to initiate divisions, signifying an early nodule primordium formation event. This is followed by cellular differentiation, including cell expansion and vascular bundle formation, and we previously showed that Lotus japonicus NF-YA1 is essential for this process, presumably by regulating three members of the SHORT INTERNODES/STYLISH (STY) transcription factor gene family. In this study, we used combined genetics, genomics and cell biology approaches to characterize the role of STY genes during root nodule formation and to test a hypothesis that they mediate nodule development by stimulating auxin signalling. We show here that L. japonicus STYs are required for nodule emergence. This is attributed to the NF-YA1-dependent regulatory cascade, comprising STY genes and their downstream targets, YUCCA1 and YUCCA11, involved in a local auxin biosynthesis at the post-initial cell division stage. An analogous NF-YA1/STY regulatory module seems to operate in Medicago truncatula in association with the indeterminate nodule patterning. Our data define L. japonicus and M. truncatula NF-YA1 genes as important nodule emergence stage-specific regulators of auxin signalling while indicating that the inductive stage and subsequent formation of early nodule primordia are mediated through an independent mechanism(s).

Massive rhizobial genomic variation associated with partner quality in Lotus-Mesorhizobium symbiosis.

Variation in partner quality is commonly observed in diverse cooperative relationships, despite the theoretical prediction that selection favoring high-quality partners should eliminate such variation. Here, we investigated how genetic variation in partner quality could be maintained in the nitrogen-fixing mutualism between Lotus japonicus and Mesorhizobium bacteria. We reconstructed de novo assembled full-genome sequences from nine rhizobial symbionts, finding massive variation in the core genome and the similar symbiotic islands, indicating recent horizontal gene transfer (HGT) of the symbiosis islands into diverse Mesorhizobium lineages. A cross-inoculation experiment using 9 sequenced rhizobial symbionts and 15 L. japonicus accessions revealed extensive quality variation represented by plant growth phenotypes, including genotype-by-genotype interactions. Variation in quality was not associated with the presence/absence variation in known symbiosis-related genes in the symbiosis island; rather, it showed significant correlation with the core genome variation. Given the recurrent HGT of the symbiosis islands into diverse Mesorhizobium strains, local Mesorhizobium communities could serve as a major source of variation for core genomes, which might prevent variation in partner quality from fixing, even in the presence of selection favoring high-quality partners. These findings highlight the novel role of HGT of symbiosis islands in maintaining partner quality variation in the legume-rhizobia symbiosis.