Geogenomic predictors of genetree heterogeneity explain phylogeographic and introgression history: a case study in an Amazonian bird (Thamnophilus aethiops).

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of four populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered two primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) gene-trees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification.

Evidence of introgression, ecological divergence and adaptation in Asterias sea stars.

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

River network rearrangements promote speciation in lowland Amazonian birds.

Large Amazonian rivers impede dispersal for many species, but lowland river networks frequently rearrange, thereby altering the location and effectiveness of river barriers through time. These rearrangements may promote biotic diversification by facilitating episodic allopatry and secondary contact among populations. We sequenced genome-wide markers to evaluate the histories of divergence and introgression in six Amazonian avian species complexes. We first tested the assumption that rivers are barriers for these taxa and found that even relatively small rivers facilitate divergence. We then tested whether species diverged with gene flow and recovered reticulate histories for all species, including one potential case of hybrid speciation. Our results support the hypothesis that river rearrangements promote speciation and reveal that many rainforest taxa are micro-endemic, unrecognized, and thus threatened with imminent extinction. We propose that Amazonian hyper-diversity originates partly from fine-scale barrier displacement processes-including river dynamics-which allow small populations to differentiate and disperse into secondary contact.

Why is Amazonia a 'source' of biodiversity? Climate-mediated dispersal and synchronous speciation across the Andes in an avian group (Tityrinae).

Amazonia is a 'source' of biodiversity for other Neotropical ecosystems, but which conditions trigger in situ speciation and emigration is contentious. Three hypotheses for how communities have assembled include (1) a stochastic model wherein chance dispersal events lead to gradual emigration and species accumulation, (2) diversity-dependence wherein successful dispersal events decline through time due to ecological limits, and (3) barrier displacement wherein environmental change facilitates dispersal to other biomes via transient habitat corridors. We sequenced thousands of molecular markers for the Neotropical Tityrinae (Aves) and applied a novel filtering protocol to identify loci with high utility for dated phylogenomics. We used these loci to estimate divergence times and model Tityrinae's evolutionary history. We detected a prominent role for speciation driven by barriers including synchronous speciation across the Andes and found that dispersal increased toward the present. Because diversification was continuous but dispersal was non-random over time, we show that barrier displacement better explains Tityrinae's history than stochasticity or diversity-dependence. We propose that Amazonia is a source of biodiversity because (1) it is a relic of a biome that was once more extensive, (2) environmentally mediated corridors facilitated emigration and (3) constant diversification is attributed to a spatially heterogeneous landscape that is perpetually dynamic through time.

Verbalizing phylogenomic conflict: Representation of node congruence across competing reconstructions of the neoavian explosion.

Phylogenomic research is accelerating the publication of landmark studies that aim to resolve deep divergences of major organismal groups. Meanwhile, systems for identifying and integrating the products of phylogenomic inference-such as newly supported clade concepts-have not kept pace. However, the ability to verbalize node concept congruence and conflict across multiple, in effect simultaneously endorsed phylogenomic hypotheses, is a prerequisite for building synthetic data environments for biological systematics and other domains impacted by these conflicting inferences. Here we develop a novel solution to the conflict verbalization challenge, based on a logic representation and reasoning approach that utilizes the language of Region Connection Calculus (RCC-5) to produce consistent alignments of node concepts endorsed by incongruent phylogenomic studies. The approach employs clade concept labels to individuate concepts used by each source, even if these carry identical names. Indirect RCC-5 modeling of intensional (property-based) node concept definitions, facilitated by the local relaxation of coverage constraints, allows parent concepts to attain congruence in spite of their differentially sampled children. To demonstrate the feasibility of this approach, we align two recent phylogenomic reconstructions of higher-level avian groups that entail strong conflict in the "neoavian explosion" region. According to our representations, this conflict is constituted by 26 instances of input "whole concept" overlap. These instances are further resolvable in the output labeling schemes and visualizations as "split concepts", which provide the labels and relations needed to build truly synthetic phylogenomic data environments. Because the RCC-5 alignments fundamentally reflect the trained, logic-enabled judgments of systematic experts, future designs for such environments need to promote a culture where experts routinely assess the intensionalities of node concepts published by our peers-even and especially when we are not in agreement with each other.

Phylogenomics and species delimitation of a complex radiation of Neotropical suboscine birds (Pachyramphus).

Phylogeographic studies within the Neotropics continue to uncover hidden diversity, the extent of which remains poorly known. In birds, molecular studies are producing evidence that species-level diversity is substantially underestimated. Many avian taxa comprise large complexes of subspecies that often represent species-level taxa by various criteria. One such group of Neotropical suboscine birds, the becards (Pachyramphus), ranges from Argentina through northern Mexico. Their taxonomic limits have been complex and controversial as the genus has bounced around a number of suboscine families. Additionally, the phylogenetic relationships within Pachyramphus are unresolved due to insufficient sampling of taxa and populations across species' ranges. We used target capture of ultraconserved elements for 62 individuals representing 42 taxa, and sequenced two mitochondrial genes and two nuclear introns covering 265 individuals of 51 taxa, including all recognized species, resulting in the most densely and completely sampled phylogenetic hypothesis for Pachyramphus to date. We delimited species using a traditional taxonomic approach and then tested them under a Bayesian multi-species coalescent framework. In doing so, we provide evidence for multiple young, previously undetected evolutionary lineages within Pachyramphus. Deep, well-supported branches and a high number of intraspecific lineages across the tree suggest that at least 50% of species diversity may be unrecognized.