Specimen collection is essential for modern science

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

Predictors of genomic diversity within North American squamates

Comparisons of intraspecific genetic diversity across species can reveal the roles of geography, ecology, and life history in shaping biodiversity. The wide availability of mitochondrial DNA (mtDNA) sequences in open-access databases makes this marker practical for conducting analyses across several species in a common framework, but patterns may not be representative of overall species diversity. Here, we gather new and existing mtDNA sequences and genome-wide nuclear data (genotyping-by-sequencing; GBS) for 30 North American squamate species sampled in the Southeastern and Southwestern United States. We estimated mtDNA nucleotide diversity for two mtDNA genes, COI (22 species alignments; average 16 sequences) and cytb (22 species; average 58 sequences), as well as nuclear heterozygosity and nucleotide diversity from GBS data for 118 individuals (30 species; four individuals and 6,820-44,309 loci per species). We showed that nuclear genomic diversity estimates were highly consistent across individuals for some species, while other species showed large differences depending on the locality sampled. Range size was positively correlated with both cytb diversity (Phylogenetically Independent Contrasts: R 2 = 0.31, p = 0.007) and GBS diversity (R 2 = 0.21; p = 0.006), while other predictors differed across the top models for each dataset. Mitochondrial and nuclear diversity estimates were not correlated within species, although sampling differences in the data available made these datasets difficult to compare. Further study of mtDNA and nuclear diversity sampled across species' ranges is needed to evaluate the roles of geography and life history in structuring diversity across a variety of taxonomic groups.

Phylogenetics And The Evolution Of Terrestriality In Mudskippers (Gobiidae: Oxudercinae)

The initial vertebrate conquest of land by stegocephalians (Sarcopterygia) allowed access to new resources and exploitation of untapped niches precipitating a major phylogenetic diversification. However, a paucity of fossils has left considerable uncertainties about phylogenetic relationships and the eco-morphological stages in this key transition in Earth history. Among extant actinopterygians, three genera of mudskippers (Gobiidae: Oxudercinae), Boleophthalmus, Periophthalmus and Periophthalmodon are the most terrestrialized, with vertebral, appendicular, locomotory, respiratory, and epithelial specializations enabling overland excursions up to 14 hours. Unlike early stegocephalians, the ecologies and morphologies of the 45 species of oxudercines are well known, making them viable analogs for the initial vertebrate conquest of land. Nevertheless, they have received little phylogenetic attention. We compiled the largest molecular dataset to date, with 29 oxudercine species, and 5 nuclear and mitochondrial loci. Phylogenetic and comparative analyses revealed strong support for two independent terrestrial transitions, and a complex suit of ecomorphological forms in estuarine environments. Furthermore, neither Oxudercinae nor their presumed sister-group the eel gobies (Amblyopinae, a group of elongated gobies) were monophyletic with respect to each other, requiring a merging of these two subfamilies and revealing an expansion of phenotypic variation within the "mudskipper" clade. We did not find support for the expected linear model of ecomorphological and locomotory transition from fully aquatic, to mudswimming, to pectoral-aided mudswimming, to lobe-finned terrestrial locomotion proposed by earlier morphological studies. This high degree of convergent or parallel transitions to terrestriality, and apparent divergent directions of estuarine adaptation, promises even greater potential for this clade to illuminate the conquest of land. Future work should focus on these less-studied species with "transitional" and other mud-habitat specializations to fully resolve the dynamics of this diversification.

Signatures of north-eastern expansion and multiple refugia: genomic phylogeography of the Pine Barrens tree frog, Hyla andersonii (Anura: Hylidae)

Range fragmentation poses challenges for species persistence over time and can be caused by both historical and contemporary processes. We combined genomic data, phylogeographical model testing and palaeoclimatic niche modelling to infer the evolutionary history of the Pine Barrens tree frog (Hyla andersonii), a seepage bog specialist, in eastern North America to gain a better understanding of the historical context of its fragmented distribution. We sampled H. andersonii populations across the three disjunct regions of the species range: Alabama/Florida (AF), the Carolinas (CL) and New Jersey (NJ). Phylogenetic relationships within H. andersonii were consistent between the nuclear species tree and mitochondrial analyses, indicating divergence between AF and CL/NJ (Atlantic clade) ~0.9 Mya and divergence of the NJ clade ~0.15 Mya. Several predictions of north-eastern expansion along the Atlantic coast were supported by phylogeographical analyses. Model testing using genome-wide single nucleotide polymorphism data and species distribution models both provided evidence for multiple disjunct refugia. This comprehensive phylogeographical study of H. andersonii demonstrates a long history of range fragmentation within an endemic coastal plain species and highlights the influence of historical climate change on the current distribution of species and their genetic diversity.

Detecting turnover among complex communities using null models: a case study with sky-island haemosporidian parasites

Turnover in species composition between sites, or beta diversity, is a critical component of species diversity that is typically influenced by geography, environment, and biotic interactions. Quantifying turnover is particularly challenging, however, in multi-host, multi-parasite assemblages where undersampling is unavoidable, resulting in inflated estimates of turnover and uncertainty about its spatial scale. We developed and implemented a framework using null models to test for community turnover in avian haemosporidian communities of three sky islands in the southwestern United States. We screened 776 birds for haemosporidian parasites from three genera (Parahaemoproteus, Plasmodium, and Leucocytozoon) by amplifying and sequencing a mitochondrial DNA barcode. We detected infections in 280 birds (36.1%), sequenced 357 infections, and found a total of 99 parasite haplotypes. When compared to communities simulated from a regional pool, we observed more unique, single-mountain haplotypes and fewer haplotypes shared among three mountain ranges than expected, indicating that haemosporidian communities differ to some degree among adjacent mountain ranges. These results were robust even after pruning datasets to include only identical sets of host species, and they were consistent for two of the three haemosporidian genera. The two more distant mountain ranges were more similar to each other than the one located centrally, suggesting that the differences we detected were due to stochastic colonization-extirpation dynamics. These results demonstrate that avian haemosporidian communities of temperate-zone forests differ on relatively fine spatial scales between adjacent sky islands. Null models are essential tools for testing the spatial scale of turnover in complex, undersampled, and poorly known systems.

Predicting amphibian intraspecific diversity with machine learning: Challenges and prospects for integrating traits, geography, and genetic data

The growing availability of genetic data sets, in combination with machine learning frameworks, offers great potential to answer long-standing questions in ecology and evolution. One such question has intrigued population geneticists, biogeographers, and conservation biologists: What factors determine intraspecific genetic diversity? This question is challenging to answer because many factors may influence genetic variation, including life history traits, historical influences, and geography, and the relative importance of these factors varies across taxonomic and geographic scales. Furthermore, interpreting the influence of numerous, potentially correlated variables is difficult with traditional statistical approaches. To address these challenges, we analysed repurposed data using machine learning and investigated predictors of genetic diversity, focusing on Nearctic amphibians as a case study. We aggregated species traits, range characteristics, and >42,000 genetic sequences for 299 species using open-access scripts and various databases. After identifying important predictors of nucleotide diversity with random forest regression, we conducted follow-up analyses to examine the roles of phylogenetic history, geography, and demographic processes on intraspecific diversity. Although life history traits were not important predictors for this data set, we found significant phylogenetic signal in genetic diversity within amphibians. We also found that salamander species at northern latitudes contained low genetic diversity. Data repurposing and machine learning provide valuable tools for detecting patterns with relevance for conservation, but concerted efforts are needed to compile meaningful data sets with greater utility for understanding global biodiversity.

Inter- and intra-archipelago dynamics of population structure and gene flow in a Polynesian bird

Islands are separated by natural barriers that prevent gene flow between terrestrial populations and promote allopatric diversification. Birds in the South Pacific are an excellent model to explore the interplay between isolation and gene flow due to the region's numerous archipelagos and well-characterized avian communities. The wattled honeyeater complex (Foulehaio spp.) comprises three allopatric species that are widespread and common across Fiji, Tonga, Samoa, and Wallis and Futuna. Here, we explored patterns of diversification within and among these lineages using genomic and morphometric data. We found support for three clades of Foulehaio corresponding to three recognized species. Within F. carunculatus, population genetic analyses identified nine major lineages, most of which were composed of sub-lineages that aligned nearly perfectly to individual island populations. Despite genetic structure and great geographic distance between populations, we found low levels of gene flow between populations in adjacent archipelagos. Additionally, body size of F. carunculatus varied randomly with respect to evolutionary history (as Ernst Mayr predicted), but correlated negatively with island size, consistent with the island rule. Our findings support a hypothesis that widespread taxa can show population structure between immediately adjacent islands, and likely represent many independent lineages loosely connected by gene flow.

Building Natural History Collections for the Twenty-First Century and Beyond

Natural history collections (NHCs) are important resources for a diverse array of scientific fields. Recent digitization initiatives have broadened the user base of NHCs, and new technological innovations are using materials generated from collections to address novel scientific questions. Simultaneously, NHCs are increasingly imperiled by reductions in funding and resources. Ensuring that NHCs continue to serve as a valuable resource for future generations will require the scientific community to increase their contribution to and acknowledgement of collections. We provide recommendations and guidelines for scientists to support NHCs, focusing particularly on new users that may be unfamiliar with collections. We hope that this perspective will motivate debate on the future of NHCs and the role of the scientific community in maintaining and improving biological collections.

Genomic sequence capture of haemosporidian parasites: Methods and prospects for enhanced study of host-parasite evolution

Avian malaria and related haemosporidians (Plasmodium, [Para]Haemoproteus and Leucocytoozoon) represent an exciting multihost, multiparasite system in ecology and evolution. Global research in this field accelerated after the publication in 2000 of PCR protocols to sequence a haemosporidian mitochondrial (mtDNA) barcode and the development in 2009 of an open-access database to document the geographic and host ranges of parasite mtDNA haplotypes. Isolating haemosporidian nuclear DNA from bird hosts, however, has been technically challenging, slowing the transition to genomic-scale sequencing techniques. We extend a recently developed sequence capture method to obtain hundreds of haemosporidian nuclear loci from wild bird samples, which typically have low levels of infection, or parasitemia. We tested 51 infected birds from Peru and New Mexico and evaluated locus recovery in light of variation in parasitemia, divergence from reference sequences and pooling strategies. Our method was successful for samples with parasitemia as low as ~0.02% (2 of 10,000 blood cells infected) and mtDNA divergence as high as 15.9% (one Leucocytozoonsample), and using the most cost-effective pooling strategy tested. Phylogenetic relationships estimated with >300 nuclear loci were well resolved, providing substantial improvement over the mtDNA barcode. We provide protocols for sample preparation and sequence capture including custom probe sequences and describe our bioinformatics pipeline using atram 2.0, phyluce and custom Perl/Python scripts. This approach can be applied to thousands of avian samples that have already been found to have haemosporidian infections of at least moderate intensity, greatly improving our understanding of parasite speciation, biogeography and evolutionary dynamics.

Deeply conserved susceptibility in a multi-host, multi-parasite system

Variation in susceptibility is ubiquitous in multi-host, multi-parasite assemblages, and can have profound implications for ecology and evolution in these systems. The extent to which susceptibility to parasites is phylogenetically conserved among hosts can be revealed by analysing diverse regional communities. We screened for haemosporidian parasites in 3983 birds representing 40 families and 523 species, spanning ~ 4500 m elevation in the tropical Andes. To quantify the influence of host phylogeny on infection status, we applied Bayesian phylogenetic multilevel models that included a suite of environmental, spatial, temporal, life history and ecological predictors. We found evidence of deeply conserved susceptibility across the avian tree; host phylogeny explained substantial variation in infection status, and results were robust to phylogenetic uncertainty. Our study suggests that susceptibility is governed, in part, by conserved, latent aspects of anti-parasite defence. This demonstrates the importance of deep phylogeny for understanding present-day ecological interactions.

Isolation and characterization of species-specific microsatellite markers for blue and black wildebeest (Connochaetes taurinus and C. gnou)

The bluewildebeest (Connochaetes taurinus) is distributed throughout southern and east Africa while the black wildebeest (Connochaetes gnou) is endemic to South Africa and was driven to near extinction in the early 1900s due to hunting pressure and disease outbreaks. Extensive translocation of both species throughout South Africa is threatening the genetic integrity of blue and blackwilde beest. To effectively manage these species, genetic tools that can be used to detect hybrid individuals, identify genetically unique subpopulations and determine the levels of genetic diversity are required. In this study, 11 microsatellite markers were developed for wildebeest through next-generation sequencing. The microsatellite loci displayed 2.00-4.14 alleles, unbiased heterozygosity values ranged from 0.32 to 0.60 and observed heterozygosity values ranged from 0.26 to 0.52. The comparatively high level of polymorphism observed in the microsatellite markers indicates that these markers can contribute significantly to our knowledge of population genetic structure, relatedness, genetic diversity and hybridization in these species.

Phylogenetic patterns of trait and trait plasticity evolution: Insights from amphibian embryos

Environmental variation favors the evolution of phenotypic plasticity. For many species, we understand the costs and benefits of different phenotypes, but we lack a broad understanding of how plastic traits evolve across large clades. Using identical experiments conducted across North America, we examined prey responses to predator cues. We quantified five life-history traits and the magnitude of their plasticity for 23 amphibian species/populations (spanning three families and five genera) when exposed to no cues, crushed-egg cues, and predatory crayfish cues. Embryonic responses varied considerably among species and phylogenetic signal was common among the traits, whereas phylogenetic signal was rare for trait plasticities. Among trait-evolution models, the Ornstein-Uhlenbeck (OU) model provided the best fit or was essentially tied with Brownian motion. Using the best fitting model, evolutionary rates for plasticities were higher than traits for three life-history traits and lower for two. These data suggest that the evolution of life-history traits in amphibian embryos is more constrained by a species' position in the phylogeny than is the evolution of life history plasticities. The fact that an OU model of trait evolution was often a good fit to patterns of trait variation may indicate adaptive optima for traits and their plasticities.

Diversity, abundance, and host relationships of avian malaria and related haemosporidians in New Mexico pine forests

Avian malaria and related haemosporidian parasites (genera Haemoproteus, Plasmodium, and Leucocytozoon) affect bird demography, species range limits, and community structure, yet they remain unsurveyed in most bird communities and populations. We conducted a community-level survey of these vector-transmitted parasites in New Mexico, USA, to describe their diversity, abundance, and host associations. We focused on the breeding-bird community in the transition zone between piñon-juniper woodland and ponderosa pine forests (elevational range: 2,150–2,460 m). We screened 186 birds representing 49 species using both standard PCR and microscopy techniques to detect infections of all three avian haemosporidian genera. We detected infections in 68 out of 186 birds (36.6%), the highest proportion of which were infected with Haemoproteus (20.9%), followed by Leucocytozoon (13.4%), then Plasmodium (8.0%). We sequenced mtDNA for 77 infections representing 43 haplotypes (25 Haemoproteus, 12 Leucocytozoon, 6 Plasmodium). When compared to all previously known haplotypes in the MalAvi and GenBank databases, 63% (27) of the haplotypes we recovered were novel. We found evidence for host specificity at the avian clade and species level, but this specificity was variable among parasite genera, in that Haemoproteus and Leucocytozoon were each restricted to three avian groups (out of six), while Plasmodium occurred in all groups except non-passerines. We found striking variation in infection rate among host species, with nearly universal infection among vireos and no infection among nuthatches. Using rarefaction and extrapolation, we estimated the total avian haemosporidian diversity to be 70 haplotypes (95% CI [43–98]); thus, we may have already sampled ∼60% of the diversity of avian haemosporidians in New Mexico pine forests. It is possible that future studies will find higher diversity in microhabitats or host species that are under-sampled or unsampled in the present study. Fortunately, this study is fully extendable via voucher specimens, frozen tissues, blood smears, parasite images, and documentation provided in open-access databases (MalAvi, GenBank, and ARCTOS).

Quantifying the spatiotemporal dynamics in a chorus frog (Pseudacris) hybrid zone over 30 years

Although theory suggests that hybrid zones can move or change structure over time, studies supported by direct empirical evidence for these changes are relatively limited. We present a spatiotemporal genetic study of a hybrid zone between Pseudacris nigrita and P. fouquettei across the Pearl River between Louisiana and Mississippi. This hybrid zone was initially characterized in 1980 as a narrow and steep “tension zone,” in which hybrid populations were inferior to parentals and were maintained through a balance between selection and dispersal. We reanalyzed historical tissue samples and compared them to samples of recently collected individuals using microsatellites. Clinal analyses indicate that the cline has not shifted in roughly 30 years but has widened significantly. Anthropogenic and natural changes may have affected selective pressure or dispersal, and our results suggest that the zone may no longer best be described as a tension zone. To the best of our knowledge, this study provides the first evidence of significant widening of a hybrid cline but stasis of its center. Continued empirical study of dynamic hybrid zones will provide insight into the forces shaping their structure and the evolutionary potential they possess for the elimination or generation of species.

Reevaluation of a classic phylogeographic barrier: new techniques reveal the influence of microgeographic climate variation on population divergence

We evaluated the mtDNA divergence and relationships within Geomys pinetis to assess the status of formerly recognized Geomys taxa. Additionally, we integrated new hypothesis-based tests in ecological niche models (ENM) to provide greater insight into causes for divergence and potential barriers to gene flow in Southeastern United States (Alabama, Florida, and Georgia). Our DNA sequence dataset confirmed and strongly supported two distinct lineages within G. pinetis occurring east and west of the ARD. Divergence date estimates showed that eastern and western lineages diverged about 1.37 Ma (1.9 Ma-830 ka). Predicted distributions from ENMs were consistent with molecular data and defined each population east and west of the ARD with little overlap. Niche identity and background similarity tests were statistically significant suggesting that ENMs from eastern and western lineages are not identical or more similar than expected based on random localities drawn from the environmental background. ENMs also support the hypothesis that the ARD represents a ribbon of unsuitable climate between more suitable areas where these populations are distributed. The estimated age of divergence between eastern and western lineages of G. pinetis suggests that the divergence was driven by climatic conditions during Pleistocene glacial-interglacial cycles. The ARD at the contact zone of eastern and western lineages of G. pinetis forms a significant barrier promoting microgeographic isolation that helps maintain ecological and genetic divergence.

Increased liver weights in severely burned children: comparison of ultrasound and autopsy measurements

Hepatomegaly is a common finding at autopsy in severely burned children surviving less than 6 months. This study validates a reliable ultrasound method which can be used to identify changes in liver size in severely burned children during acute hospitalization. Thirty-eight children, age 0.5-17 years with burns covering over 40% of their total surface area were studied at autopsy. Liver weight was measured at autopsy and compared to predicted liver weight for age and height. Eighteen had liver size measured by ultrasound within 10 days of death while five had ultrasound liver measures after death just prior to autopsy. All burned children who survived 7 days or more (n = 33) had liver weights at autopsy that were greater than predicted for age and height while all 23 livers measured by ultrasound were greater than predicted. Autopsy weights correlated well with weights estimated by ultrasound, R = 0.824. At autopsy, those who survived 7 days or more had enlarged livers ranging from 142 to 406% of their predicted normal age and height. Common histologic findings include large and small-droplet fat deposits, and cholestasis. The degree of these histologic abnormalities correlated with the increase in liver weight, R = 0.652. Ultrasound is a valid, noninvasive method for measuring liver weight changes in severely burned children during acute hospitalization. Ninety-five percent of the severely burned children from this institute had significant hepatomegaly identified at autopsy.

Chest radiographic appearances in severely burned adults. A comparison of early radiographic and extravascular lung thermal volume changes

Chest radiographs (CXRs) have previously been used as a diagnostic tool to detect changes in lung water. In this study CXR changes in severely burned adults, in the absence of an inhalation injury, preceded detectable increases in extravascular lung thermal volume (ELTV) by 3 to 5 days. The hypothesis that early CXR density changes in burned patients have an infectious cause, not related to changes in ELTV, was tested. Blood cultures, CXRs, and ELTV were evaluated during the first 15 days after injury in severely burned adults who had no identified inhalation injury. Chest radiographs were scored daily on a 1 to 5 scale, with 1 = normal, 2 = peribronchial cuffing, 3 = mild interstitial infiltrates, 4 = severe interstitial infiltrates, and 5 = alveolar infiltrates. In all patients, except those who were septic, increases in their CXR density scores correlated well with increases in ELTV. The ELTV/CXR score ratios for septic burn patients on days 1 to 6 postburn was 1.7 +/- 0.2 compared with 4.2 +/- 0.4, (means +/- SEM) for nonseptic (P < .001), whereas the ELTV/CXR score ratios for septic and nonseptic patients, 7 to 15 days postburn, were 3.8 +/- 0.4 and 3.4 +/- 0.5, respectively. We suggest that before any measurable change in ELTV early increases in CXR density scores in burned patients without a concomitant inhalation injury are caused by intraalveolar pneumonitis or hyaline membrane atelectasis and not increased ELTV.