The root nodule symbiosis of plants with nitrogen-fixing bacteria affects global nitrogen cycles and food production but is restricted to a subset of genera within a single clade of flowering plants. To explore the genetic basis for this scattered occurrence, we sequenced the genomes of 10 plant species covering the diversity of nodule morphotypes, bacterial symbionts, and infection strategies. In a genome-wide comparative analysis of a total of 37 plant species, we discovered signatures of multiple independent loss-of-function events in the indispensable symbiotic regulator in 10 of 13 genomes of nonnodulating species within this clade. The discovery that multiple independent losses shaped the present-day distribution of nitrogen-fixing root nodule symbiosis in plants reveals a phylogenetically wider distribution in evolutionary history and a so-far-underestimated selection pressure against this symbiosis.

Frankia strains are nitrogen-fixing soil actinobacteria that can form root symbioses with actinorhizal plants. Phylogenetically, symbiotic frankiae can be divided into three clusters, and this division also corresponds to host specificity groups. The strains of cluster II which form symbioses with actinorhizal Rosales and Cucurbitales, thus displaying a broad host range, show suprisingly low genetic diversity and to date can not be cultured. The genome of the first representative of this cluster, Candidatus Frankia datiscae Dg1 (Dg1), a microsymbiont of Datisca glomerata, was recently sequenced. A phylogenetic analysis of 50 different housekeeping genes of Dg1 and three published Frankia genomes showed that cluster II is basal among the symbiotic Frankia clusters. Detailed analysis showed that nodules of D. glomerata, independent of the origin of the inoculum, contain several closely related cluster II Frankia operational taxonomic units. Actinorhizal plants and legumes both belong to the nitrogen-fixing plant clade, and bacterial signaling in both groups involves the common symbiotic pathway also used by arbuscular mycorrhizal fungi. However, so far, no molecules resembling rhizobial Nod factors could be isolated from Frankia cultures. Alone among Frankia genomes available to date, the genome of Dg1 contains the canonical nod genes nodA, nodB and nodC known from rhizobia, and these genes are arranged in two operons which are expressed in D. glomerata nodules. Furthermore, Frankia Dg1 nodC was able to partially complement a Rhizobium leguminosarum A34 nodC::Tn5 mutant. Phylogenetic analysis showed that Dg1 Nod proteins are positioned at the root of both α- and β-rhizobial NodABC proteins. NodA-like acyl transferases were found across the phylum Actinobacteria, but among Proteobacteria only in nodulators. Taken together, our evidence indicates an Actinobacterial origin of rhizobial Nod factors.

Frankia strains induce the formation of nitrogen-fixing nodules on roots of actinorhizal plants. Phylogenetically, Frankia strains can be grouped in four clusters. The earliest divergent cluster, cluster-2, has a particularly wide host range. The analysis of cluster-2 strains has been hampered by the fact that with two exceptions, they could never be cultured. In this study, 12 Frankia-enriched metagenomes of Frankia cluster-2 strains or strain assemblages were sequenced based on seven inoculum sources. Sequences obtained via DNA isolated from whole nodules were compared with those of DNA isolated from fractionated preparations enhanced in the Frankia symbiotic structures. The results show that cluster-2 inocula represent groups of strains, and that strains not represented in symbiotic structures, that is, unable to perform symbiotic nitrogen fixation, may still be able to colonize nodules. Transposase gene abundance was compared in the different Frankia-enriched metagenomes with the result that North American strains contain more transposase genes than Eurasian strains. An analysis of the evolution and distribution of the host plants indicated that bursts of transposition may have coincided with niche competition with other cluster-2 Frankia strains. The first genome of an inoculum from the Southern Hemisphere, obtained from nodules of Coriaria papuana in Papua New Guinea, represents a novel species, postulated as Candidatus Frankia meridionalis. All Frankia-enriched metagenomes obtained in this study contained homologs of the canonical nod genes nodABC; the North American genomes also contained the sulfotransferase gene nodH, while the genome from the Southern Hemisphere only contained nodC and a truncated copy of nodB.

BackgroundThe ability to establish root nodule symbioses is restricted to four different plant orders. Soil actinobacteria of the genus Frankia can establish a symbiotic relationship with a diverse group of plants within eight different families from three different orders, the Cucurbitales, Fagales and Rosales. Phylogenetically, Frankia strains can be divided into four clusters, three of which (I, II, III) contain symbiotic strains. Members of Cluster II nodulate the broadest range of host plants with species from four families from two different orders, growing on six continents. Two Cluster II genomes were sequenced thus far, both from Asia.ResultsIn this paper we present the first Frankia cluster II genome from North America (California), Dg2, which represents a metagenome of two major and one minor strains. A phylogenetic analysis of the core genomes of 16 Frankia strains shows that Cluster II the ancestral group in the genus, also ancestral to the non-symbiotic Cluster IV. Dg2 contains the canonical nod genes nodABC for the production of lipochitooligosaccharide Nod factors, but also two copies of the sulfotransferase gene nodH. In rhizobial systems, sulfation of Nod factors affects their host specificity and their stability.ConclusionsA comparison with the nod gene region of the previously sequenced Dg1 genome from a Cluster II strain from Pakistan shows that the common ancestor of both strains should have contained nodABC and nodH. Phylogenetically, Dg2 NodH proteins are sister to rhizobial NodH proteins. A glnA-based phylogenetic analysis of all Cluster II strains sampled thus far supports the hypothesis that Cluster II Frankia strains came to North America with Datisca glomerata following the Madrean-Tethyan pattern.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3140-1) contains supplementary material, which is available to authorized users.