Rapid radiations outweigh reticulations during the evolution of a 750-million-year-old lineage of cyanobacteria.

Species are a fundamental unit of biodiversity. Yet, the existence of clear species boundaries among bacteria has long been a subject of debate. Here, we studied species boundaries in the context of the phylogenetic history of Nostoc, a widespread genus of photoautotrophic and nitrogen-fixing cyanobacteria that includes many lineages that form symbiotic associations with plants (e.g., cycads and bryophytes) and fungi (e.g., cyanolichens). We found that the evolution of Nostoc was characterized by eight rapid radiations, many of which were associated with major events in the evolution of plants. In addition, incomplete lineage sorting associated with these rapid radiations outweighed reticulations during Nostoc evolution. We then show that the pattern of diversification of Nostoc shapes the distribution of average nucleotide identities (ANIs) into a complex mosaic, wherein some closely related clades are clearly isolated from each other by gaps in genomic similarity, while others form a continuum where genomic species boundaries are expected. Nevertheless, recently diverged Nostoc lineages often form cohesive clades that are maintained by within-clade gene flow. Boundaries to homologous recombination between these cohesive clades persist even when the potential for gene flow is high, i.e., when closely related clades of Nostoc cooccur or are locally found in symbiotic associations with the same lichen-forming fungal species. Our results demonstrate that rapid radiations are major contributors to the complex speciation history of Nostoc. This underscores the need to consider evolutionary information beyond thresholds of genomic similarity to delimit biologically meaningful units of biodiversity for bacteria.