Abstract Root nodule symbiosis (RNS) allows plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. RNS is a complex trait requiring coordination from both the plant host and the bacterial symbiont, and pinpointing the evolutionary origins of root nodules is critical for understanding the genetic basis of RNS. This endeavor is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for all major lineages of the NFC and an enhanced root nodule trait database to reconstruct the evolution of RNS. Through identification of the evolutionary pathway to RNS gain, we show that shifts among heterogeneous evolutionary rates can explain how a complex trait such as RNS can arise many times across a large phylogeny. Our analysis identifies a two-step process in which an ancestral precursor state gave rise to a more labile state from which RNS was quickly gained at specific points in the NFC. Our rigorous reconstruction of ancestral states illustrates how a two-step pathway could have led to multiple independent gains and losses of RNS, contrary to recent hypotheses invoking just a single gain and numerous losses. RNS may be an example of multi-level convergent evolution, thus requiring a broader phylogenetic and genetic scope for genome-phenome mapping to elucidate mechanisms enabling fully functional RNS.

Two shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants