Evolutionary biology of birds and their associated microbiota in the Madrean Archipelago sky islands

Sky islands are disjunct patches of montane forested habitat in a matrix of desert, grasslands, and scrub. I investigated intraspecific evolutionary biology of two bird species—the White­‐breasted Nuthatch (Sitta carolinensis) and the Brown Creeper (Certhia americana)—in the Madrean Archipelago sky islands (Arizona, USA), a biodiversity and evolution hotspot. In addition, I explored patterns of codiversification of these two birds and their associated gut microbial communities. The two bird studies revealed different patterns of diversification within the sky islands. The White­‐breasted Nuthatch exhibits a pattern of isolation by environment, where genetic differences among populations are related to environmental differences of those localities. In contrast, the Brown Creeper has a strong genetic break between northern and southern populations, with no evidence of gene flow between lineages. When I investigated codiversification of birds and their microbial communities, I found no relationship between host genetic diversity and microbial community alpha diversity, while genetic differentiation between birds was significantly related with beta diversity between microbial communities. This dissertation provides a first step in comparative evolutionary biology of Madrean Archipelago avian taxa, and adds to the knowledge of the factors shaping microbial community diversity in wild animals

Comparison of Target-Capture and Restriction-Site Associated DNA Sequencing for Phylogenomics: A Test in Cardinalid Tanagers (Aves, Genus: Piranga)

This is a pre-copyedited, author-produced version of an article accepted for publication in Systematic Biology following peer review. The version of record, Joseph D. Manthey, Luke C. Campillo, Kevin J. Burns, Robert G. Moyle; Comparison of Target-Capture and Restriction-Site Associated DNA Sequencing for Phylogenomics: A Test in Cardinalid Tanagers (Aves, Genus: Piranga), Systematic Biology, Volume 65, Issue 4, 1 July 2016, Pages 640–650, is available online at: ttps://doi.org/10.1093/sysbio/syw005.Restriction-site associated DNA sequencing (RAD-seq) and target capture of specific genomic regions, such as ultraconserved elements (UCEs), are emerging as two of the most popular methods for phylogenomics using reduced-representation genomic data sets. These two methods were designed to target different evolutionary timescales: RAD-seq was designed for population-genomic level questions and UCEs for deeper phylogenetics. The utility of both data sets to infer phylogenies across a variety of taxonomic levels has not been adequately compared within the same taxonomic system. Additionally, the effects of uninformative gene trees on species tree analyses (for target capture data) have not been explored. Here, we utilize RAD-seq and UCE data to infer a phylogeny of the bird genus Piranga. The group has a range of divergence dates (0.5–6 myr), contains 11 recognized species, and lacks a resolved phylogeny. We compared two species tree methods for the RAD-seq data and six species tree methods for the UCE data. Additionally, in the UCE data, we analyzed a complete matrix as well as data sets with only highly informative loci. A complete matrix of 189 UCE loci with 10 or more parsimony informative (PI) sites, and an approximately 80% complete matrix of 1128 PI single-nucleotide polymorphisms (SNPs) (from RAD-seq) yield the same fully resolved phylogeny of Piranga. We inferred non-monophyletic relationships of Pirangalutea individuals, with all other a priori species identified as monophyletic. Finally, we found that species tree analyses that included predominantly uninformative gene trees provided strong support for different topologies, with consistent phylogenetic results when limiting species tree analyses to highly informative loci or only using less informative loci with concatenation or methods meant for SNPs alone

A genome-wide assessment of stages of elevational parapatry in Bornean passerine birds reveals no introgression: implications for processes and patterns of speciation

Topographically complex regions often contain the close juxtaposition of closely related species along elevational gradients. The evolutionary causes of these elevational replacements, and thus the origin and maintenance of a large portion of species diversity along elevational gradients, are usually unclear because ecological differentiation along a gradient or secondary contact following allopatric diversification can produce the same pattern. We used reduced representation genomic sequencing to assess genetic relationships and gene flow between three parapatric pairs of closely related songbird taxa (Arachnothera spiderhunters, Chloropsis leafbirds, and Enicurus forktails) along an elevational gradient in Borneo. Each taxon pair presents a different elevational range distribution across the island, yet results were uniform: little or no gene flow was detected in any pairwise comparisons. These results are congruent with an allopatric “species-pump” model for generation of species diversity and elevational parapatry of congeners on Borneo, rather than in situ generation of species by “ecological speciation” along an elevational gradient

Tectonic collision and uplift of Wallacea triggered the global songbird radiation

Songbirds (oscine passerines) are the most species-rich and cosmopolitan bird group, comprising almost half of global avian diversity. Songbirds originated in Australia, but the evolutionary trajectory from a single species in an isolated continent to worldwide proliferation is poorly understood. Here, we combine the first comprehensive genome-scale DNA sequence data set for songbirds, fossil-based time calibrations, and geologically informed biogeographic reconstructions to provide a well-supported evolutionary hypothesis for the group. We show that songbird diversification began in the Oligocene, but accelerated in the early Miocene, at approximately half the age of most previous estimates. This burst of diversification occurred coincident with extensive island formation in Wallacea, which provided the first dispersal corridor out of Australia, and resulted in independent waves of songbird expansion through Asia to the rest of the globe. Our results reconcile songbird evolution with Earth history and link a major radiation of terrestrial biodiversity to early diversification within an isolated Australian continent

Data from: Isolation by environment in white-breasted nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach

Understanding landscape processes driving patterns of population genetic differentiation and diversity has been a long-standing focus of ecology and evolutionary biology. Gene flow may be reduced by historical, ecological or geographic factors, resulting in patterns of isolation by distance (IBD) or isolation by environment (IBE). Although IBE has been found in many natural systems, most studies investigating patterns of IBD and IBE in nature have used anonymous neutral genetic markers, precluding inference of selection mechanisms or identification of genes potentially under selection. Using landscape genomics, the simultaneous study of genomic and ecological landscapes, we investigated the processes driving population genetic patterns of White-breasted Nuthatches (Sitta carolinensis) in sky islands (montane forest habitat islands) of the Madrean Archipelago. Using more than 4000 single nucleotide polymorphisms and multiple tests to investigate the relationship between genetic differentiation and geographic or ecological distance, we identified IBE, and a lack of IBD, among sky island populations of S. carolinensis. Using three tests to identify selection, we found 79 loci putatively under selection; of these, seven matched CDS regions in the Zebra Finch. The loci under selection were highly associated with climate extremes (maximum temperature of warmest month and minimum precipitation of driest month). These results provide evidence for IBE – disentangled from IBD – in sky island vertebrates and identify potential adaptive genetic variation

Cryptic diversity in a widespread North American songbird: Phylogeography of the Brown Creeper (Certhia americana)

The identification of species via morphological characteristics has traditionally left cryptic species undescribed in taxa under selection for morphological conservation (or a lack of selection for morphological change). Treecreepers (Genus: Certhia) have a conserved morphological appearance, making it difficult to ascertain relationships in the genus based on morphology alone. Recent genetic and song structure studies of Eurasian Treecreepers identified cryptic species within Old World Certhia that were previously undescribed using morphological characteristics. Here, we use mtDNA to investigate cryptic diversity and patterns of diversification in the Brown Creeper (Certhia americana), the single described Certhia species in the Americas. Phylogenetic analyses identified six well-supported geographically-structured clades; the basal divergence separates a northern and a southern lineage in the Brown Creeper, likely cryptic species previously characterized as many subspecies. Sympatry is prevalent between clades in western North America, where possible contact zones warrant further investigation. Allopatry appears to be the primary driver of deep phylogeographic structure within the Brown Creeper; however, within clade diversity is highly correlated with the life history traits of the populations that comprise the geographically structured phylogroups

Impact of host demography and evolutionary history on endosymbiont molecular evolution: A test in carpenter ants (genus Camponotus) and their Blochmannia endosymbionts

Abstract Obligate endosymbioses are tight associations between symbionts and the hosts they live inside. Hosts and their associated obligate endosymbionts generally exhibit codiversification, which has been documented in taxonomically diverse insect lineages. Host demography (e.g., effective population sizes) may impact the demography of endosymbionts, which may lead to an association between host demography and the patterns and processes of endosymbiont molecular evolution. Here, we used whole‐genome sequencing data for carpenter ants (Genus Camponotus; subgenera Camponotus and Tanaemyrmex) and their Blochmannia endosymbionts as our study system to address whether Camponotus demography shapes Blochmannia molecular evolution. Using whole‐genome phylogenomics, we confirmed previous work identifying codiversification between carpenter ants and their Blochmannia endosymbionts. We found that Blochmannia genes have evolved at a pace ~30× faster than that of their hosts' molecular evolution and that these rates are positively associated with host rates of molecular evolution. Using multiple tests for selection in Blochmannia genes, we found signatures of positive selection and shifts in selection strength across the phylogeny. Host demography was associated with Blochmannia shifts toward increased selection strengths, but not associated with Blochmannia selection relaxation, positive selection, genetic drift rates, or genome size evolution. Mixed support for relationships between host effective population sizes and Blochmannia molecular evolution suggests weak or uncoupled relationships between host demography and Blochmannia population genomic processes. Finally, we found that Blochmannia genome size evolution was associated with genome‐wide estimates of genetic drift and number of genes with relaxed selection pressures