A rapid and cost-effective pipeline for digitization of museum specimens with 3D photogrammetry.

Natural history collections are yielding more information as digitization brings specimen data to researchers, connects specimens across museums, and as new technologies allow for more large-scale data collection. Therefore, a key goal in specimen digitization is developing methods that both increase access and allow for the highest yield of phenomic data. 3D digitization is increasingly popular because it has the potential to meet both aspects of that key goal. However, current methods overlook or do not prioritize some of the most sought-after phenotypic traits, those involving the external appearance of specimens, especially color. Here, we introduce an efficient and cost-effective pipeline for 3D photogrammetry to capture the external appearance of natural history specimens and other museum objects. 3D photogrammetry aligns and compares sets of dozens, hundreds, or even thousands of photos to create 3D models. The hardware set-up requires little physical space and around $3,000 in initial investment, while the software pipeline requires $1,400/year for proprietary software subscriptions (with open-source alternatives). The creation of each 3D model takes 1-2 hours/specimen and much of the software pipeline is automated with minimal supervision required, including the onerous step of mesh processing. We showcase the method by creating 3D models for most of the type specimens in the Moore Laboratory of Zoology bird collection and show that digital bill measurements are comparable to hand-taken measurements. Color data, while not included as part of this pipeline, is easily extractable from the models and one of the most promising areas of data collection. Future advances can adapt the method for ultraviolet reflectance capture and increased efficiency and model quality. Combined with genomic data, phenomic data from 3D models including photogrammetry will open new doors to understanding organismal evolution.

Comparative Analyses of Vertebrate Gut Microbiomes Reveal Convergence between Birds and Bats.

Diet and host phylogeny drive the taxonomic and functional contents of the gut microbiome in mammals, yet it is unknown whether these patterns hold across all vertebrate lineages. Here, we assessed gut microbiomes from ∼900 vertebrate species, including 315 mammals and 491 birds, assessing contributions of diet, phylogeny, and physiology to structuring gut microbiomes. In most nonflying mammals, strong correlations exist between microbial community similarity, host diet, and host phylogenetic distance up to the host order level. In birds, by contrast, gut microbiomes are only very weakly correlated to diet or host phylogeny. Furthermore, while most microbes resident in mammalian guts are present in only a restricted taxonomic range of hosts, most microbes recovered from birds show little evidence of host specificity. Notably, among the mammals, bats host especially bird-like gut microbiomes, with little evidence for correlation to host diet or phylogeny. This suggests that host-gut microbiome phylosymbiosis depends on factors convergently absent in birds and bats, potentially associated with physiological adaptations to flight. Our findings expose major variations in the behavior of these important symbioses in endothermic vertebrates and may signal fundamental evolutionary shifts in the cost/benefit framework of the gut microbiome.IMPORTANCE In this comprehensive survey of microbiomes of >900 species, including 315 mammals and 491 birds, we find a striking convergence of the microbiomes of birds and animals that fly. In nonflying mammals, diet and short-term evolutionary relatedness drive the microbiome, and many microbial species are specific to a particular kind of mammal, but flying mammals and birds break this pattern with many microbes shared across different species, with little correlation either with diet or with relatedness of the hosts. This finding suggests that adaptation to flight breaks long-held relationships between hosts and their microbes.

More than skin and bones: Comparing extraction methods and alternative sources of DNA from avian museum specimens.

Next-generation sequencing has greatly expanded the utility and value of museum collections by revealing specimens as genomic resources. As the field of museum genomics grows, so does the need for extraction methods that maximize DNA yields. For avian museum specimens, the established method of extracting DNA from toe pads works well for most specimens. However, for some specimens, especially those of birds that are very small or very large, toe pads can be a poor source of DNA. In this study, we apply two DNA extraction methods (phenol-chloroform and silica column) to three different sources of DNA (toe pad, skin punch and bone) from 10 historical avian museum specimens. We show that a modified phenol-chloroform protocol yielded significantly more DNA than a silica column protocol (e.g., Qiagen DNeasy Blood & Tissue Kit) across all tissue types. However, extractions using the silica column protocol contained longer fragments on average than those using the phenol-chloroform protocol, probably as a result of loss of small fragments through the silica column. While toe pads yielded more DNA than skin punches and bone fragments, skin punches proved to be a reliable alternative source of DNA and might be especially appealing when toe pad extractions are impractical. Overall, we found that historical bird museum specimens contain substantial amounts of DNA for genomic studies under most extraction scenarios, but that a phenol-chloroform protocol consistently provides the high quantities of DNA required for most current genomic protocols.

Museum genomics reveals the speciation history of Dendrortyx wood-partridges in the Mesoamerican highlands.

Natural history collections are increasingly valued as genomic resources. Their specimens reflect the combined efforts of collectors and curators over hundreds of years. For many rare or endangered species, specimens are the only readily available source of DNA. We leveraged specimens from a historical collection to study the evolutionary history of wood-partridges in the genus Dendrortyx. The three Dendrortyx species are found in the highlands of central Mexico and Central America south to Costa Rica. One of these species is endangered, and in general, Dendrortyx are secretive and poorly represented in tissue collections. We extracted DNA from historical museum specimens and sequenced ultraconserved elements (UCEs) and mitochondrial DNA (mtDNA) to assess their phylogeny and divergence times. Phylogenies built from hundreds to thousands of nuclear markers were well resolved and largely congruent with an mtDNA phylogeny. The divergence times revealed an unusually old avian divergence across the Isthmus of Tehuantepec in the Pliocene around 3.6 million years ago. Combined with other recent studies, our results challenge the general pattern that highland bird divergences in Mesoamerica are relatively young and influenced by the Pleistocene glacial cycles compared to the older divergences of reptiles and plants, which are thought to overlap more with periods of mountain formation. We also found evidence for monophyletic genetic lineages in mountain ranges within the widespread D. macroura, which should be investigated further with integrative taxonomic methods. Our study demonstrates the power of museum genomics to provide insight into the evolutionary histories of groups where modern samples are lacking.

Combining ultraconserved elements and mtDNA data to uncover lineage diversity in a Mexican highland frog (Sarcohyla; Hylidae).

Molecular studies have uncovered significant diversity in the Mexican Highlands, leading to the description of many new endemic species. DNA approaches to this kind of species discovery have included both mitochondrial DNA (mtDNA) sequencing and multilocus genomic methods. While these marker types have often been pitted against one another, there are benefits to deploying them together, as linked mtDNA data can provide the bridge between uncovering lineages through rigorous multilocus genomic analysis and identifying lineages through comparison to existing mtDNA databases. Here, we apply one class of multilocus genomic marker, ultraconserved elements (UCEs), and linked mtDNA data to a species complex of frogs (Sarcohyla bistincta, Hylidae) found in the Mexican Highlands. We generated data from 1,891 UCEs, which contained 1,742 informative SNPs for S. bistincta and closely related species and captured mitochondrial genomes for most samples. Genetic analyses based on both whole loci and SNPs agree there are six to seven distinct lineages within what is currently described as S. bistincta. Phylogenies from UCEs and mtDNA mostly agreed in their topologies, and the few differences suggested a more complex evolutionary history of the mtDNA marker. Our study demonstrates that the Mexican Highlands still hold substantial undescribed diversity, making their conservation a particularly urgent goal. The Trans-Mexican Volcanic Range stands out as a significant geographic feature in Sarcohyla and may have acted as a dispersal corridor for S. bistincta to spread to the north. Combining multilocus genomic data with linked mtDNA data is a useful approach for identifying potential new species and associating them with already described taxa, which will be especially important in groups with undescribed subadult phenotypes and cryptic species.

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.

Speciation in Western Scrub-Jays, Haldane's rule, and genetic clines in secondary contact.

BACKGROUND: Haldane's Rule, the tendency for the heterogametic sex to show reduced fertility in hybrid crosses, can obscure the signal of gene flow in mtDNA between species where females are heterogametic. Therefore, it is important when studying speciation and species limits in female-heterogametic species like birds to assess the signature of gene flow in the nuclear genome as well. We studied introgression of microsatellites and mtDNA across a secondary contact zone between coastal and interior lineages of Western Scrub-Jays (Aphelocoma californica) to test for a signature of Haldane's Rule: a narrower cline of introgression in mtDNA compared to nuclear markers. RESULTS: Our initial phylogeographic analysis revealed that there is only one major area of contact between coastal and interior lineages and identified five genetic clusters with strong spatial structuring: Pacific Slope, Interior US, Edwards Plateau (Texas), Northern Mexico, and Southern Mexico. Consistent with predictions from Haldane's Rule, mtDNA showed a narrower cline than nuclear markers across a transect through the hybrid zone. This result is not being driven by female-biased dispersal because neutral diffusion analysis, which included estimates of sex-specific dispersal rates, also showed less diffusion of mtDNA. Lineage-specific plumage traits were associated with nuclear genetic profiles for individuals in the hybrid zone, indicating that these differences are under genetic control. CONCLUSIONS: This study adds to a growing list of studies that support predictions of Haldane's Rule using cline analysis of multiple loci of differing inheritance modes, although alternate hypotheses like selection on different mtDNA types cannot be ruled out. That Haldane's Rule appears to be operating in this system suggests a measure of reproductive isolation between the Pacific Slope and interior lineages. Based on a variety of evidence from the phenotype, ecology, and genetics, we recommend elevating three lineages to species level: A. californica (Pacific Slope); A. woodhouseii (Interior US plus Edwards Plateau plus Northern Mexico); A. sumichrasti (Southern Mexico). The distinctive Edwards Plateau population in Texas, which was monophyletic in mtDNA except for one individual, should be studied in greater detail given habitat threat.

HZAR: hybrid zone analysis using an R software package.

We present a new software package (HZAR) that provides functions for fitting molecular genetic and morphological data from hybrid zones to classic equilibrium cline models using the Metropolis-Hastings Markov chain Monte Carlo (MCMC) algorithm. The software applies likelihood functions appropriate for different types of data, including diploid and haploid genetic markers and quantitative morphological traits. The modular design allows flexibility in fitting cline models of varying complexity. To facilitate hypothesis testing, an autofit function is included that allows automated model selection from a set of nested cline models. Cline parameter values, such as cline centre and cline width, are estimated and may be compared statistically across clines. The package is written in the R language and is available through the Comprehensive R Archive Network (CRAN; http://cran.r-project.org/). Here, we describe HZAR and demonstrate its use with a sample data set from a well-studied hybrid zone in western Panama between white-collared (Manacus candei) and golden-collared manakins (M. vitellinus). Comparisons of our results with previously published results for this hybrid zone validate the hzar software. We extend analysis of this hybrid zone by fitting additional models to molecular data where appropriate.

Next-generation sequencing reveals phylogeographic structure and a species tree for recent bird divergences.

Next generation sequencing (NGS) technologies are revolutionizing many biological disciplines but have been slow to take root in phylogeography. This is partly due to the difficulty of using NGS to sequence orthologous DNA fragments for many individuals at low cost. We explore cases of recent divergence in four phylogenetically diverse avian systems using a method for quick and cost-effective generation of primary DNA sequence data using pyrosequencing. NGS data were processed using an analytical pipeline that reduces many reads into two called alleles per locus per individual. Using single nucleotide polymorphisms (SNPs) mined from the loci, we detected population differentiation in each of the four bird systems, including: a case of ecological speciation in rails (Rallus); a rapid postglacial radiation in the genus Junco; recent in situ speciation among hummingbirds (Trochilus) in Jamaica; and subspecies of white-crowned sparrows (Zonotrichia leucophrys) along the Pacific coast. The number of recovered loci aligning closely to chromosomal locations on the zebra finch (Taeniopygia guttata) genome was highly correlated to the size of the chromosome, suggesting that loci are randomly distributed throughout the genome. Using eight loci found in Zonotrichia and Junco lineages, we were also able to generate a species tree of these sparrow sister genera, demonstrating the potential of this method for generating data amenable to coalescent-based analysis. We discuss improvements that should enhance the method's utility for primary data generation.

Diversification at high latitudes: speciation of buntings in the genus Plectrophenax inferred from mitochondrial and nuclear markers.

High-latitude diversification is a process characterized by speciation and extinction due to climatically driven vicariance and dispersal events. McKay's buntings (Plectrophenax hyperboreus) are high-latitude island endemic songbirds, and their global range is restricted to Beringia. Snow buntings (P. nivalis), their closest relatives, are distributed throughout the Holarctic, breeding in available habitat surrounding the island range of McKay's buntings. We sequenced 1123 base pairs of mitochondrial DNA for 40 individuals of each species and analysed a total of 913 AFLPs for 57 individuals. Both marker types suggested weak but significant genetic differentiation. Analysis of sequence data indicated divergence occurring when the current breeding range of McKay's buntings was a hill on the Beringian steppe (approximately 18,400 to approximately 73,700 years before present), suggesting that snow buntings were restricted to lower latitudes by ice sheets. Ancestral effective population size estimates indicate a founder event in McKay's buntings followed by an expansion and then a reduction in effective size. Rising sea levels and asymmetric hybridization from McKay's buntings into the postglacially-colonizing population of snow buntings could account for this reduction. Reproductive isolation is likely maintained through differential arrival dates on breeding grounds and the high breeding density of McKay's buntings. This recent, high-latitude divergence best fits a model of founder event speciation driven by vicariance and oscillations in habitat due to climate change.

Paraphyly of Cinclodes fuscus (Aves: Passeriformes: Furnariidae): implications for taxonomy and biogeography.

The Andes are a hotspot of global avian diversity, but studies on the historical diversification of Andean birds remain relatively scarce. Evolutionary studies on avian lineages with Andean-Patagonian distributions have focused on reconstructing species-level phylogenies, whereas no detailed phylogeographic studies on widespread species have been conducted. Here, we describe phylogeographic patterns in the Bar-winged Cinclodes (Cinclodes fuscus), a widespread and common species of ovenbird (Furnariidae) that breeds from Tierra del Fuego to the northern Andes. Traditionally, C. fuscus has been considered a single species composed of nine subspecies, but its long and narrow range suggests the possibility of considerable genetic variation among populations. Sequences of two mitochondrial genes revealed three discrete and geographically coherent groups of C. fuscus, occupying the southern, central, and northern Andes. Surprisingly, phylogenetic analyses indicated that these groups were more closely related to other species of Cinclodes than to each other. Relationships of the southern and northern C. fuscus clades to other species of Cinclodes were straightforward; in combination with available information on plumage, behavioral, and vocal variation, this suggests that each should be recognized as a distinct biological species. The central Andean group was paraphyletic with respect to C. oustaleti, and relationships among these taxa and C. olrogi were poorly resolved. We suggest that the central Andean C. fuscus should also be considered a different species, pending new information to clarify species limits in this group. These new phylogenetic data, along with recently developed methods, allowed us to review the biogeography of the genus, confirming southern South America and the central Andes as important areas for the diversification of these birds.