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.

Cannabaceae are a relatively small family of angiosperms, but they include several species of huge economic and cultural significance: marijuana or hemp (Cannabis sativa) and hops (Humulus lupulus). Previous phylogenetic studies have clarified the most deep relationships in Cannabaceae, but relationships remain ambiguous among several major lineages. Here, we sampled 82 species representing all genera of Cannabaceae and utilized a new dataset of 90 nuclear genes and 82 chloroplast loci from Hyb‐Seq to investigate the phylogenomics of Cannabaceae. Nuclear phylogenetic analyses revealed a robust and consistent backbone for Cannabaceae. We observed nuclear gene‐tree conflict at several deep nodes in inferred species trees, also cyto‐nuclear discordance concerning the relationship between Gironniera and Lozanella and the relationships among Trema s.l. (including Parasponia), Cannabis + Humulus, and Chaetachme + Pteroceltis. Coalescent simulations and network analyses suggest that observed deep cyto‐nuclear discordances were most likely to stem from incomplete lineage sorting (ILS); nuclear gene‐tree conflict might be caused by both ILS and gene flow between species. All genera of Cannabaceae were recovered as monophyletic, except for Celtis, which consisted of two distinct clades: Celtis I (including most Celtis species) and Celtis II (including Celtis gomphophylla and Celtis schippii). We suggest that Celtis II should be recognized as the independent genus Sparrea based on both molecular and morphological evidence. Our work provides the most comprehensive and reliable phylogeny to date for Cannabaceae, enabling further exploration of evolutionary patterns across this family and highlighting the necessity of comparing nuclear with chloroplast data to examine the evolutionary history of plant groups.