Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria.

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Moving beyond the "Diversity Paradox": The Limitations of Competition-Based Frameworks in Understanding Species Diversity.

AbstractOver the past century, ecologists have attempted to understand patterns of species diversity by showing stable coexistence arising from a baseline expectation of competitive exclusion. This expectation stems from an explicit assumption of resource scarcity and implicit assumptions of Malthusian struggle and winner-takes-all dynamics. Fidelity to the competitive exclusion principle (CEP) presents species diversity as a paradox: if species compete for limited resources, how can they coexist? In this article, we investigate the contradiction between the theoretical expectation of competitive exclusion and the empirical prevalence of multispecies communities. We trace the persistence of the CEP in ecological research despite numerous challenges and explore publishing trends suggesting that this framework has resulted in a disproportionate focus on competition and exclusion in contemporary research. From a critical science studies perspective, we analyze the sociopolitical factors that have contributed to these patterns. We argue that we must excavate the ideological foundation on which competition-based coexistence research has been built to move beyond the current perceived "diversity paradox." To that end, we propose shifting the baseline expectation of coexistence research, introducing the notion of a coexistence principle, which positions the persistence of multispecies communities as the rule rather than the exception in nature.

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks.

Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α-diversity. To calculate the PSS index, we then compare observed partner phylogenetic α-diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host-parasite, avian seed-dispersal, lichenized fungi-cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much-needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss.

Turnover of Lecanoroid Mycobionts and Their Trebouxia Photobionts Along an Elevation Gradient in Bolivia Highlights the Role of Environment in Structuring the Lichen Symbiosis.

Shifts in climate along elevation gradients structure mycobiont-photobiont associations in lichens. We obtained mycobiont (lecanoroid Lecanoraceae) and photobiont (Trebouxia alga) DNA sequences from 89 lichen thalli collected in Bolivia from a ca. 4,700 m elevation gradient encompassing diverse natural communities and environmental conditions. The molecular dataset included six mycobiont loci (ITS, nrLSU, mtSSU, RPB1, RPB2, and MCM7) and two photobiont loci (ITS, rbcL); we designed new primers to amplify Lecanoraceae RPB1 and RPB2 with a nested PCR approach. Mycobionts belonged to Lecanora s.lat., Bryonora, Myriolecis, Protoparmeliopsis, the "Lecanora" polytropa group, and the "L." saligna group. All of these clades except for Lecanora s.lat. occurred only at high elevation. No single species of Lecanoraceae was present along the entire elevation gradient, and individual clades were restricted to a subset of the gradient. Most Lecanoraceae samples represent species which have not previously been sequenced. Trebouxia clade C, which has not previously been recorded in association with species of Lecanoraceae, predominates at low- to mid-elevation sites. Photobionts from Trebouxia clade I occur at the upper extent of mid-elevation forest and at some open, high-elevation sites, while Trebouxia clades A and S dominate open habitats at high elevation. We did not find Trebouxia clade D. Several putative new species were found in Trebouxia clades A, C, and I. These included one putative species in clade A associated with Myriolecis species growing on limestone at high elevation and a novel lineage sister to the rest of clade C associated with Lecanora on bark in low-elevation grassland. Three different kinds of photobiont switching were observed, with certain mycobiont species associating with Trebouxia from different major clades, species within a major clade, or haplotypes within a species. Lecanoraceae mycobionts and Trebouxia photobionts exhibit species turnover along the elevation gradient, but with each partner having a different elevation threshold at which the community shifts completely. A phylogenetically defined sampling of a single diverse family of lichen-forming fungi may be sufficient to document regional patterns of Trebouxia diversity and distribution.