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 4 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 2 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) genetrees 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. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

An island 'endemic' born out of hybridization between introduced lineages

Humans have profoundly impacted the distribution of plant and animal species over thousands of years. The most direct example of these effects is human-mediated movement of individuals, either through translocation of individuals within their range or through the introduction of species to new habitats. While human involvement may be suspected in species with obvious range disjunctions, it can be difficult to detect natural versus human-mediated dispersal events for populations at the edge of a species' range, and this uncertainty muddles how we understand the evolutionary history of populations and broad biogeographical patterns. Studies combining genetic data with archaeological, linguistic and historical evidence have confirmed prehistoric examples of human-mediated dispersal; however, it is unclear whether these methods can disentangle recent dispersal events, such as species translocated by European colonizers during the past 500 years. We use genomic DNA from historical museum specimens and historical records to evaluate three hypotheses regarding the timing and origin of Northern Bobwhites (Colinus virginianus) in Cuba, whose status as an endemic or introduced population has long been debated. We discovered that bobwhites from southern Mexico arrived in Cuba between the 12th and 16th centuries, followed by the subsequent introduction of bobwhites from the southeastern USA to Cuba between the 18th and 20th centuries. These dates suggest the introduction of bobwhites to Cuba was human-mediated and concomitant with Spanish colonial shipping routes between Veracruz, Mexico and Havana, Cuba during this period. Our results identify endemic Cuban bobwhites as a genetically distinct population born of hybridization between divergent, introduced lineages.

Implications of headwater contact zones for the riverine barrier hypothesis: a case study of the Blue-capped Manakin (Lepidothrix coronata)

Rivers frequently delimit the geographic ranges of species in the Amazon Basin. These rivers also define the boundaries between genetic clusters within many species, yet river boundaries have been documented to break down in headwater regions where rivers are narrower. To explore the evolutionary implications of headwater contact zones in Amazonia, we examined genetic variation in the Blue-capped Manakin (Lepidothrix coronata), a species previously shown to contain several genetically and phenotypically distinct populations across the western Amazon Basin. We collected restriction site-associated DNA sequence data (RADcap) for 706 individuals and found that spatial patterns of genetic structure indicate several rivers, particularly the Amazon and Ucayali, are dispersal barriers for L. coronata. We also found evidence that genetic connectivity is elevated across several headwater regions, highlighting the importance of headwater gene flow for models of Amazonian diversification. The headwater region of the Ucayali River provided a notable exception to findings of headwater gene flow by harboring non-admixed populations of L. coronata on opposite sides of a < 1-km-wide river channel with a known dynamic history, suggesting that additional prezygotic barriers may be limiting gene flow in this region.

Speciation-by-Extinction

Extinction is a dominant force shaping patterns of biodiversity through time; however its role as a catalyst of speciation through its interaction with intraspecific variation has been overlooked. Here, we synthesize ideas alluded to by Darwin and others into the model of "speciation-by-extinction" in which speciation results from the extinction of intermediate populations within a single geographically variable species. We explore the properties and distinguishing features of speciation-by-extinction with respect to other established speciation models. We demonstrate its plausibility by showing that the experimental extinction of populations within variable species can result in speciation. The prerequisites for speciation-by-extinction, geographically structured intraspecific variation and local extinction, are ubiquitous in nature. We propose that speciation-by-extinction may be a prevalent, but underappreciated, speciation mechanism.

Systematics of the Neopelminae (Aves: Passeriformes: Pipridae) with description of a new genus

N/A.

A high-quality de novo genome assembly for clapper rail (Rallus crepitans)

The clapper rail (Rallus crepitans), of the family Rallidae, is a secretive marsh bird species that is adapted for high salinity habitats. They are very similar in appearance to the closely related king rail (R. elegans), but while king rails are limited primarily to freshwater marshes, clapper rails are highly adapted to tolerate salt marshes. Both species can be found in brackish marshes where they freely hybridize, but the distribution of their respective habitats precludes the formation of a continuous hybrid zone and secondary contact can occur repeatedly. This system, thus, provides unique opportunities to investigate the underlying mechanisms driving their differential salinity tolerance as well as the maintenance of the species boundary between the 2 species. To facilitate these studies, we assembled a de novo reference genome assembly for a female clapper rail. Chicago and HiC libraries were prepared as input for the Dovetail HiRise pipeline to scaffold the genome. The pipeline, however, did not recover the Z chromosome so a custom script was used to assemble the Z chromosome. We generated a near chromosome level assembly with a total length of 994.8 Mb comprising 13,226 scaffolds. The assembly had a scaffold N50 was 82.7 Mb, L50 of four, and had a BUSCO completeness score of 92%. This assembly is among the most contiguous genomes among the species in the family Rallidae. It will serve as an important tool in future studies on avian salinity tolerance, interspecific hybridization, and speciation.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Amazonian birds in more dynamic habitats have less population genetic structure and higher gene flow

Understanding the factors that govern variation in genetic structure across species is key to the study of speciation and population genetics. Genetic structure has been linked to several aspects of life history, such as foraging strategy, habitat association, migration distance, and dispersal ability, all of which might influence dispersal and gene flow. Comparative studies of population genetic data from species with differing life histories provide opportunities to tease apart the role of dispersal in shaping gene flow and population genetic structure. Here, we examine population genetic data from sets of bird species specialized on a series of Amazonian habitat types hypothesized to filter for species with dramatically different dispersal abilities: stable upland forest, dynamic floodplain forest, and highly dynamic riverine islands. Using genome-wide markers, we show that habitat type has a significant effect on population genetic structure, with species in upland forest, floodplain forest, and riverine islands exhibiting progressively lower levels of structure. Although morphological traits used as proxies for individual-level dispersal ability did not explain this pattern, population genetic measures of gene flow are elevated in species from more dynamic riverine habitats. Our results suggest that the habitat in which a species occurs drives the degree of population genetic structuring via its impact on long-term fluctuations in levels of gene flow, with species in highly dynamic habitats having particularly elevated gene flow. These differences in genetic variation across taxa specialized in distinct habitats may lead to disparate responses to environmental change or habitat-specific diversification dynamics over evolutionary time scales.

Phylogenomic analysis of the parrots of the world distinguishes artifactual from biological sources of gene tree discordance

Gene tree discordance is expected in phylogenomic trees and biological processes are often invoked to explain it. However, heterogeneous levels of phylogenetic signal among individuals within datasets may cause artifactual sources of topological discordance. We examined how the information content in tips and subclades impacts topological discordance in the parrots (Order: Psittaciformes), a diverse and highly threatened clade of nearly 400 species. Using ultraconserved elements from 96% of the clade's species-level diversity, we estimated concatenated and species trees for 382 ingroup taxa. We found that discordance among tree topologies was most common at nodes dating between the late Miocene and Pliocene, and often at the taxonomic level of genus. Accordingly, we used two metrics to characterize information content in tips and assess the degree to which conflict between trees was being driven by lower quality samples. Most instances of topological conflict and non-monophyletic genera in the species tree could be objectively identified using these metrics. For subclades still discordant after tip-based filtering, we used a machine learning approach to determine whether phylogenetic signal or noise was the more important predictor of metrics supporting the alternative topologies. We found that when signal favored one of the topologies, noise was the most important variable in poorly performing models that favored the alternative topology. In sum, we show that artifactual sources of gene tree discordance, which are likely a common phenomenon in many datasets, can be distinguished from biological sources by quantifying the information content in each tip and modeling which factors support each topology.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae)

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

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.

The dynamics of introgression across an avian radiation

Hybridization and resulting introgression can play both a destructive and a creative role in the evolution of diversity. Thus, characterizing when and where introgression is most likely to occur can help us understand the causes of diversification dynamics. Here, we examine the prevalence of and variation in introgression using phylogenomic data from a large (1300+ species), geographically widespread avian group, the suboscine birds. We first examine patterns of gene tree discordance across the geographic distribution of the entire clade. We then evaluate the signal of introgression in a subset of 206 species triads using Patterson's D‐statistic and test for associations between introgression signal and evolutionary, geographic, and environmental variables. We find that gene tree discordance varies across lineages and geographic regions. The signal of introgression is highest in cases where species occur in close geographic proximity and in regions with more dynamic climates since the Pleistocene. Our results highlight the potential of phylogenomic datasets for examining broad patterns of hybridization and suggest that the degree of introgression between diverging lineages might be predictable based on the setting in which they occur.

A reference genome for the nectar-robbing Black-throated Flowerpiercer (Diglossa brunneiventris)

Black-throated Flowerpiercers (Diglossa brunneiventris) are one species representing a phenotypically specialized group of tanagers (Thraupidae) that have hooked bills which allow them to feed by stealing nectar from the base of flowers. Members of the genus are widely distributed in montane regions from Mexico to northern Argentina, and previous studies of Diglossa have focused on their systematics, phylogenetics, and interesting natural history. Despite numerous studies of species within the genus, no genome assembly exists to represent these nectivorous tanagers. We described the assembly of a genome sequence representing a museum-vouchered, wild, female D. brunneiventris collected in Peru. By combining Pacific Biosciences Sequel long-read technology with 10× linked-read and reference-based scaffolding, we produced a 1.08 Gbp pseudochromosomal assembly including 600 scaffolds with a scaffold N50 of 67.3 Mbp, a scaffold L50 of 6, and a BUSCO completeness score of 95%. This new assembly improves representation of the diverse species that comprise the tanagers, improves on scaffold lengths and contiguity when compared to existing genomic resources for tanagers, and provides another avenue of research into the genetic basis of adaptations common to a nectivorous lifestyle among vertebrates.

Displaced clines in an avian hybrid zone (Thamnophilidae: Rhegmatorhina) within an Amazonian interfluve

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi.

Multiple species and deep genomic divergences despite little phenotypic differentiation in an ancient Neotropical songbird, Tunchiornis ochraceiceps (Sclater, 1860) (Aves: Vireonidae)

Several bird taxa have been recently described or elevated to full species and almost twice as many bird species than are currently recognized may exist. Defining species is one of the most basic and important issues in biological science because unknown or poorly defined species hamper subsequent studies. Here, we evaluate the species limits and evolutionary history of Tunchiornis ochraceiceps-a widespread forest songbird that occurs in the lowlands of Central America, Chocó and Amazonia-using an integrative approach that includes plumage coloration, morphometrics, vocalization and genomic data. The species has a relatively old crown age (~9 Ma) and comprises several lineages with little, if any, evidence of gene flow among them. We propose a taxonomic arrangement composed of four species, three with a plumage coloration diagnosis and one deeply divergent cryptic species. Most of the remaining lineages have variable but unfixed phenotypic characters despite their relatively old origin. This decoupling of genomic and phenotypic differentiation reveals a remarkable case of phenotypic conservatism, possibly due to strict habitat association. Lineages are geographically delimited by the main Amazonian rivers and the Andes, a pattern observed in studies of other understory upland forest Neotropical birds, although phylogenetic relationships and divergence times among populations are idiosyncratic.

Rethinking Gloger's Rule: Climate, Light Environments, and Color in a Large Family of Tropical Birds (Furnariidae)

AbstractEcogeographic rules provide a framework within which to test evolutionary hypotheses of adaptation. Gloger's rule predicts that endothermic animals should have darker colors in warm/rainy climates. This rule also predicts that animals should be more rufous in warm/dry climates, the so-called complex Gloger's rule. Empirical studies frequently demonstrate that animals are darker in cool/wet climates rather than in warm/wet climates. Furthermore, sensory ecology predicts that, to enhance crypsis, animals should be darker in darker light environments. We aimed to disentangle the effects of climate and light environments on plumage color in the large Neotropical passerine family Furnariidae. We found that birds in cooler and rainier climates had darker plumage even after controlling for habitat type. Birds in darker habitats had darker plumage even after controlling for climate. The effects of temperature and precipitation interact so that the negative effect of precipitation on brightness is strongest in cool temperatures. Finally, birds tended to be more rufous in warm/dry habitats but also, surprisingly, in cool/wet locales. We suggest that Gloger's rule results from complementary selective pressures arising from myriad ecological factors, including crypsis, thermoregulation, parasite deterrence, and resistance to feather abrasion.

The evolution of a tropical biodiversity hotspot

The tropics are the source of most biodiversity yet inadequate sampling obscures answers to fundamental questions about how this diversity evolves. We leveraged samples assembled over decades of fieldwork to study diversification of the largest tropical bird radiation, the suboscine passerines. Our phylogeny, estimated using data from 2389 genomic regions in 1940 individuals of 1283 species, reveals that peak suboscine species diversity in the Neotropics is not associated with high recent speciation rates but rather with the gradual accumulation of species over time. Paradoxically, the highest speciation rates are in lineages from regions with low species diversity, which are generally cold, dry, unstable environments. Our results reveal a model in which species are forming faster in environmental extremes but have accumulated in moderate environments to form tropical biodiversity hotspots.

A simple index to quantify and compare the magnitude of intraspecific geographic plumage colour variation in typical antbirds (Aves: Passeriformes: Thamnophilidae)

Intraspecific geographic phenotypic variation is a crucial theme in evolutionary biology. Comparing its magnitude across species can provide insights into its ecological and genetic correlates. Here, we developed an index, which we dub the V index, to quantify intraspecific plumage colour variation in typical antbirds (Thamnophilidae), a family which has long interested ornithologists due to a high prevalence of intraspecific variation. The V index is based on a bivariate colour space defined by brightness and redness. Its value for each species equals the mean area occupied by each of its subspecies in that colour space, divided by the area of the species. Lower values indicate greater intraspecific geographic variation. Based on this index, Thamnophilus caerulescens (Variable Antshrike) was exceptionally geographically variable compared to other thamnophilids, as previously suggested based on qualitative evidence. In general, we found that the most variable species had disjunct distributions and deep phylogeographic structure, suggesting an effect of historical population dynamics in producing geographic variation. The V index can be adapted for use with other taxa, traits, and taxonomic levels, and we expect it will instigate novel ways of thinking about phenotypic variation in birds and other animals.

Investigating the utility of traditional and genomic multi-locus datasets to resolve relationships in Lipaugus and Tijuca (Cotingidae)

Rapid diversification limits our ability to resolve evolutionary relationships and examine diversification history, as in the case of the Neotropical cotingas. Here we present an analysis with complete taxon sampling for the cotinga genera Lipaugus and Tijuca, which include some of the most range-restricted (e.g., T. condita) and also the most widespread and familiar (e.g., L. vociferans) forest birds in the Neotropics. We used two datasets: (1) Sanger sequencing data sampled from eight loci in 34 individuals across all described taxa and (2) sequence capture data linked to 1,079 ultraconserved elements and conserved exons sampled from one or two individuals per species. Phylogenies estimated from the Sanger sequencing data failed to resolve three nodes, but the sequence capture data produced a well-supported tree. Lipaugus and Tijuca formed a single, highly supported clade, but Tijuca species were not sister and were embedded within Lipaugus. A dated phylogeny confirmed Lipaugus and Tijuca diversified rapidly in the Miocene. Our study provides a detailed evolutionary hypothesis for Lipaugus and Tijuca and demonstrates that increasing genomic sampling can prove instrumental in resolving the evolutionary history of recent radiations.

The dual role of Amazonian rivers in the generation and maintenance of avian diversity

The Amazon River and its major tributaries delimit the distributions of hundreds of terrestrial taxa. It remains unclear whether river-bounded distributions and taxon replacements reflect the historical role of rivers in generating species diversity as vicariant forces, or are the result of their role as secondary barriers, maintaining current levels of species diversity by inhibiting gene flow and population introgression. We use a community-wide comparative phylogeographic and phylogenetic approach to address the roles that the Rio Negro and the Rio Branco play in the avian speciation process in the Guiana Shield. Examining 74 pairs of ecologically similar geographic replacements that turn over across the lower Negro, we found substantial variation in the levels of genetic divergence and the inferred timing of diversification among pairs, ranging from ~0.24 to over 8 million years (Ma ago). The breadth of this variation is inconsistent with a single, shared speciation event. Coalescent simulations also rejected a simultaneous divergence scenario for pairs divided by the Rio Branco but could not reject a single diversification pulse for a subset of 12 pairs of taxa divided by the upper Negro. These results are consistent with recent geomorphological hypotheses regarding the origins of these rivers. Phylogenetically, taxon pairs represent a blend of sister (~40%) and nonsister taxa (~60%), consistent with river-associated allopatric or peripatric speciation and secondary contact, respectively. Our data provide compelling evidence that species turnover across the Rio Negro basin encompasses a mixture of histories, supporting a dual role for Amazonian rivers in the generation and maintenance of biological diversity.

Habitat Association Predicts Genetic Diversity and Population Divergence in Amazonian Birds

The ecological traits of organisms may predict their genetic diversity and population genetic structure and mediate the action of evolutionary processes important for speciation and adaptation. Making these ecological-evolutionary links is difficult because it requires comparable genetic estimates from many species with differing ecologies. In Amazonian birds, habitat association is an important component of ecological diversity. Here, we examine the link between habitat association and genetic parameters using 20 pairs of closely related Amazonian bird species in which one member of the pair occurs primarily in forest edge and floodplains and the other occurs in upland forest interior. We use standardized geographic sampling and data from 2,416 genomic markers to estimate genetic diversity, population genetic structure, and statistics reflecting demographic and evolutionary processes. We find that species of upland forest have greater genetic diversity and divergence across the landscape as well as signatures of older histories and less gene flow than floodplain species. Our results reveal that species ecology in the form of habitat association is an important predictor of genetic diversity and population divergence and suggest that differences in diversity between floodplain and upland avifaunas in the Amazon may be driven by differences in the demographic and evolutionary processes at work in the two habitats.

Positive association between population genetic differentiation and speciation rates in New World birds

An implicit assumption of speciation biology is that population differentiation is an important stage of evolutionary diversification, but its significance as a rate-limiting control on phylogenetic speciation dynamics remains largely untested. If population differentiation within a species is related to its speciation rate over evolutionary time, the causes of differentiation could also be driving dynamics of organismal diversity across time and space. Alternatively, geographic variants might be short-lived entities with rates of formation that are unlinked to speciation rates, in which case the causes of differentiation would have only ephemeral impacts. By pairing population genetics datasets from 173 New World bird species (>17,000 individuals) with phylogenetic estimates of speciation rate, we show that the population differentiation rates within species are positively correlated with their speciation rates over long timescales. Although population differentiation rate explains relatively little of the variation in speciation rate among lineages, the positive relationship between differentiation rate and speciation rate is robust to species-delimitation schemes and to alternative measures of both rates. Population differentiation occurs at least three times faster than speciation, which suggests that most populations are ephemeral. Speciation and population differentiation rates are more tightly linked in tropical species than in temperate species, consistent with a history of more stable diversification dynamics through time in the Tropics. Overall, our results suggest that the processes responsible for population differentiation are tied to those that underlie broad-scale patterns of diversity.

HISTORICAL DIVERSIFICATION OF BIRDS IN NORTHWESTERN SOUTH AMERICA: A MOLECULAR PERSPECTIVE ON THE ROLE OF VICARIANT EVENTS

Studies of South American biodiversity have identified several areas of endemism that may have enhanced historical diversification of South American organisms. Hypotheses concerning the derivation of birds in the Chocó area of endemism in northwestern South America were evaluated using protein electrophoretic data from 14 taxonomically diverse species groups of birds. Nine of these groups demonstrated that the Chocó area of endemism has a closer historical relationship to Central America than to Amazonia, a result that is consistent with phytogeographic evidence. Within species groups, genetic distances between cis-Andean (east of the Andes) and trans-Andean (west of the Andes) taxa are, on average, roughly twice that between Chocó and Central American taxa. The genetic data are consistent with the hypotheses that the divergence of most cis-Andean and trans-Andean taxa was the result of either the Andean uplift fragmenting a once continuous Amazonian-Pacific population (Andean Uplift Hypothesis), the isolation of the two faunas in forest refugia on opposite sides of the Andes during arid climates (Forest Refugia Hypothesis), or dispersal of Amazonian forms directly across the Andes into the trans-Andean region (Across-Andes Dispersal Hypothesis). Disentangling these hypotheses is difficult due to the complexity of the Andean uplift and to the scant geologic and paleoclimatic information that elucidates diversification events in northwestern South America. Regarding the divergence of cis- and trans-Andean taxa, the genetic, geologic, and paleoclimatic data allow weak rejection of the Andean Uplift Hypothesis and weak support for the Forest Refugia and Andean Dispersal Hypotheses. The subsequent diversification of Chocó and Central American taxa was the result of Pleistocene forest refugia, marine transgressions, or parapatric speciation.

A latitudinal phylogeographic diversity gradient in birds

High tropical species diversity is often attributed to evolutionary dynamics over long timescales. It is possible, however, that latitudinal variation in diversification begins when divergence occurs within species. Phylogeographic data capture this initial stage of diversification in which populations become geographically isolated and begin to differentiate genetically. There is limited understanding of the broader implications of intraspecific diversification because comparative analyses have focused on species inhabiting and evolving in restricted regions and environments. Here, we scale comparative phylogeography up to the hemisphere level and examine whether the processes driving latitudinal differences in species diversity are also evident within species. We collected genetic data for 210 New World bird species distributed across a broad latitudinal gradient and estimated a suite of metrics characterizing phylogeographic history. We found that lower latitude species had, on average, greater phylogeographic diversity than higher latitude species and that intraspecific diversity showed evidence of greater persistence in the tropics. Factors associated with species ecologies, life histories, and habitats explained little of the variation in phylogeographic structure across the latitudinal gradient. Our results suggest that the latitudinal gradient in species richness originates, at least partly, from population-level processes within species and are consistent with hypotheses implicating age and environmental stability in the formation of diversity gradients. Comparative phylogeographic analyses scaled up to large geographic regions and hundreds of species can show connections between population-level processes and broad-scale species-richness patterns.

Sequence Capture versus Restriction Site Associated DNA Sequencing for Shallow Systematics

Sequence capture and restriction site associated DNA sequencing (RAD-Seq) are two genomic enrichment strategies for applying next-generation sequencing technologies to systematics studies. At shallow timescales, such as within species, RAD-Seq has been widely adopted among researchers, although there has been little discussion of the potential limitations and benefits of RAD-Seq and sequence capture. We discuss a series of issues that may impact the utility of sequence capture and RAD-Seq data for shallow systematics in non-model species. We review prior studies that used both methods, and investigate differences between the methods by re-analyzing existing RAD-Seq and sequence capture data sets from a Neotropical bird (Xenops minutus). We suggest that the strengths of RAD-Seq data sets for shallow systematics are the wide dispersion of markers across the genome, the relative ease and cost of laboratory work, the deep coverage and read overlap at recovered loci, and the high overall information that results. Sequence capture's benefits include flexibility and repeatability in the genomic regions targeted, success using low-quality samples, more straightforward read orthology assessment, and higher per-locus information content. The utility of a method in systematics, however, rests not only on its performance within a study, but on the comparability of data sets and inferences with those of prior work. In RAD-Seq data sets, comparability is compromised by low overlap of orthologous markers across species and the sensitivity of genetic diversity in a data set to an interaction between the level of natural heterozygosity in the samples examined and the parameters used for orthology assessment. In contrast, sequence capture of conserved genomic regions permits interrogation of the same loci across divergent species, which is preferable for maintaining comparability among data sets and studies for the purpose of drawing general conclusions about the impact of historical processes across biotas. We argue that sequence capture should be given greater attention as a method of obtaining data for studies in shallow systematics and comparative phylogeography.

Comparative Gut Microbiota of 59 Neotropical Bird Species

The gut microbiota of vertebrates are essential to host health. Most non-model vertebrates, however, lack even a basic description of natural gut microbiota biodiversity. Here, we sampled 116 intestines from 59 Neotropical bird species and used the V6 region of the 16S rRNA molecule as a microbial fingerprint (average coverage per bird ~80,000 reads). A core microbiota of Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria was identified, as well as several gut-associated genera. We tested 18 categorical variables associated with each bird for significant correlation to the gut microbiota; host taxonomic categories were most frequently significant and explained the most variation. Ecological variables (e.g., diet, foraging stratum) were also frequently significant but explained less variation. Little evidence was found for a significant influence of geographic space. Finally, we suggest that microbial sampling during field collection of organisms would propel biological understanding of evolutionary history and ecological significance of host-associated microbiota.

Similarity thresholds used in DNA sequence assembly from short reads can reduce the comparability of population histories across species

Comparing inferences among datasets generated using short read sequencing may provide insight into the concerted impacts of divergence, gene flow and selection across organisms, but comparisons are complicated by biases introduced during dataset assembly. Sequence similarity thresholds allow the de novo assembly of short reads into clusters of alleles representing different loci, but the resulting datasets are sensitive to both the similarity threshold used and to the variation naturally present in the organism under study. Thresholds that require high sequence similarity among reads for assembly (stringent thresholds) as well as highly variable species may result in datasets in which divergent alleles are lost or divided into separate loci ('over-splitting'), whereas liberal thresholds increase the risk of paralogous loci being combined into a single locus ('under-splitting'). Comparisons among datasets or species are therefore potentially biased if different similarity thresholds are applied or if the species differ in levels of within-lineage genetic variation. We examine the impact of a range of similarity thresholds on assembly of empirical short read datasets from populations of four different non-model bird lineages (species or species pairs) with different levels of genetic divergence. We find that, in all species, stringent similarity thresholds result in fewer alleles per locus than more liberal thresholds, which appears to be the result of high levels of over-splitting. The frequency of putative under-splitting, conversely, is low at all thresholds. Inferred genetic distances between individuals, gene tree depths, and estimates of the ancestral mutation-scaled effective population size (θ) differ depending upon the similarity threshold applied. Relative differences in inferences across species differ even when the same threshold is applied, but may be dramatically different when datasets assembled under different thresholds are compared. These differences not only complicate comparisons across species, but also preclude the application of standard mutation rates for parameter calibration. We suggest some best practices for assembling short read data to maximize comparability, such as using more liberal thresholds and examining the impact of different thresholds on each dataset.

Whole-genome analyses resolve early branches in the tree of life of modern birds

To better determine the history of modern birds, we performed a genome-scale phylogenetic analysis of 48 species representing all orders of Neoaves using phylogenomic methods created to handle genome-scale data. We recovered a highly resolved tree that confirms previously controversial sister or close relationships. We identified the first divergence in Neoaves, two groups we named Passerea and Columbea, representing independent lineages of diverse and convergently evolved land and water bird species. Among Passerea, we infer the common ancestor of core landbirds to have been an apex predator and confirm independent gains of vocal learning. Among Columbea, we identify pigeons and flamingoes as belonging to sister clades. Even with whole genomes, some of the earliest branches in Neoaves proved challenging to resolve, which was best explained by massive protein-coding sequence convergence and high levels of incomplete lineage sorting that occurred during a rapid radiation after the Cretaceous-Paleogene mass extinction event about 66 million years ago.

Geographic variation and phylogenetic relationships of Myiopagis olallai (Aves: Passeriformes; Tyrannidae), with the description of two new taxa from the Northern Andes

Geographic variation in vocalizations, morphology and plumage patterns in New World flycatchers is little understood, particularly in rare species with disjunct distributions. We discovered a distinct new flycatcher of the genus Myiopagis from cloud forests of the northern Central Andes in Antioquia, Colombia. Comparisons of vocalizations and external morphology, and molecular phylogenetic analyses, demonstrate that the "Antioquia Myiopagis" is a unique lineage of the M. caniceps-olallai group. We show that three specimens collected in 1940-1951 from cloud forests of Serranía de Perijá in Venezuela, and traditionally assigned to M. caniceps, represent another distinct taxon that is closer to the "Antioquia Myiopagis" and M. olallai. Both new taxa, from Antioquia and Perijá, are described as subspecies of M. olallai. We present a phylogenetic hypothesis for the M. caniceps-olallai group, in which M. olallai and the "Antioquia Myiopagis" are phylogenetically nested within the polytypic M. caniceps, which consists of at least four distinct lineages, indicating that species diversity in this group could be underestimated. 

The drivers of tropical speciation

Since the recognition that allopatric speciation can be induced by large-scale reconfigurations of the landscape that isolate formerly continuous populations, such as the separation of continents by plate tectonics, the uplift of mountains or the formation of large rivers, landscape change has been viewed as a primary driver of biological diversification. This process is referred to in biogeography as vicariance. In the most species-rich region of the world, the Neotropics, the sundering of populations associated with the Andean uplift is ascribed this principal role in speciation. An alternative model posits that rather than being directly linked to landscape change, allopatric speciation is initiated to a greater extent by dispersal events, with the principal drivers of speciation being organism-specific abilities to persist and disperse in the landscape. Landscape change is not a necessity for speciation in this model. Here we show that spatial and temporal patterns of genetic differentiation in Neotropical birds are highly discordant across lineages and are not reconcilable with a model linking speciation solely to landscape change. Instead, the strongest predictors of speciation are the amount of time a lineage has persisted in the landscape and the ability of birds to move through the landscape matrix. These results, augmented by the observation that most species-level diversity originated after episodes of major Andean uplift in the Neogene period, suggest that dispersal and differentiation on a matrix previously shaped by large-scale landscape events was a major driver of avian speciation in lowland Neotropical rainforests.

Adaptation to ephemeral habitat may overcome natural barriers and severe habitat fragmentation in a fire-dependent species, the Bachman's Sparrow (Peucaea aestivalis)

Bachman's Sparrow (Peucaea aestivalis) is a fire-dependent species that has undergone range-wide population declines in recent decades. We examined genetic diversity in Bachman's Sparrows to determine whether natural barriers have led to distinct population units and to assess the effect of anthropogenic habitat loss and fragmentation. Genetic diversity was examined across the geographic range by genotyping 226 individuals at 18 microsatellite loci and sequencing 48 individuals at mitochondrial and nuclear genes. Multiple analyses consistently demonstrated little genetic structure and high levels of genetic variation, suggesting that populations are panmictic. Based on these genetic data, separate management units/subspecies designations or translocations to promote gene flow among fragmented populations do not appear to be necessary. Panmixia in Bachman's Sparrow may be a consequence of an historical range expansion and retraction. Alternatively, high vagility in Bachman's Sparrow may be an adaptation to the ephemeral, fire-mediated habitat that this species prefers. In recent times, high vagility also appears to have offset inbreeding and loss of genetic diversity in highly fragmented habitat.

Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater)

Brown-headed Cowbirds (Molothrus ater) are the most widespread avian brood parasite in North America, laying their eggs in the nests of approximately 250 host species that raise the cowbird nestlings as their own. It is currently unknown how these heterospecific hosts influence the cowbird gut microbiota relative to other factors, such as the local environment and genetics. We test a Nature Hypothesis (positing the importance of cowbird genetics) and a Nurture Hypothesis (where the host parents are most influential to cowbird gut microbiota) using the V6 region of 16S rRNA as a microbial fingerprint of the gut from 32 cowbird samples and 16 potential hosts from nine species. We test additional hypotheses regarding the influence of the local environment and age of the birds. We found no evidence for the Nature Hypothesis and little support for the Nurture Hypothesis. Cowbird gut microbiota did not form a clade, but neither did members of the host species. Rather, the physical location, diet and age of the bird, whether cowbird or host, were the most significant categorical variables. Thus, passerine gut microbiota may be most strongly influenced by environmental factors. To put this variation in a broader context, we compared the bird data to a fecal microbiota dataset of 38 mammal species and 22 insect species. Insects were always the most variable; on some axes, we found more variation within cowbirds than across all mammals. Taken together, passerine gut microbiota may be more variable and environmentally determined than other taxonomic groups examined to date.

Species coexistence and the dynamics of phenotypic evolution in adaptive radiation

Interactions between species can promote evolutionary divergence of ecological traits and social signals, a process widely assumed to generate species differences in adaptive radiation. However, an alternative view is that lineages typically interact when relatively old, by which time selection for divergence is weak and potentially exceeded by convergent selection acting on traits mediating interspecific competition. Few studies have tested these contrasting predictions across large radiations, or by controlling for evolutionary time. Thus the role of species interactions in driving broad-scale patterns of trait divergence is unclear. Here we use phylogenetic estimates of divergence times to show that increased trait differences among coexisting lineages of ovenbirds (Furnariidae) are explained by their greater evolutionary age in relation to non-interacting lineages, and that--when these temporal biases are accounted for--the only significant effect of coexistence is convergence in a social signal (song). Our results conflict with the conventional view that coexistence promotes trait divergence among co-occurring organisms at macroevolutionary scales, and instead provide evidence that species interactions can drive phenotypic convergence across entire radiations, a pattern generally concealed by biases in age.