The Fossil Calibration Database-A New Resource for Divergence Dating.

Fossils provide the principal basis for temporal calibrations, which are critical to the accuracy of divergence dating analyses. Translating fossil data into minimum and maximum bounds for calibrations is the most important-often least appreciated-step of divergence dating. Properly justified calibrations require the synthesis of phylogenetic, paleontological, and geological evidence and can be difficult for nonspecialists to formulate. The dynamic nature of the fossil record (e.g., new discoveries, taxonomic revisions, updates of global or local stratigraphy) requires that calibration data be updated continually lest they become obsolete. Here, we announce the Fossil Calibration Database (http://fossilcalibrations.org), a new open-access resource providing vetted fossil calibrations to the scientific community. Calibrations accessioned into this database are based on individual fossil specimens and follow best practices for phylogenetic justification and geochronological constraint. The associated Fossil Calibration Series, a calibration-themed publication series at Palaeontologia Electronica, will serve as a key pipeline for peer-reviewed calibrations to enter the database.

A multi-locus inference of the evolutionary diversification of extant flamingos (Phoenicopteridae).

BACKGROUND: Modern flamingos (Phoenicopteridae) occupy a highly specialized ecology unique among birds and represent a potentially powerful model system for informing the mechanisms by which a lineage of birds adapts and radiates. However, despite a rich fossil record and well-studied feeding morphology, molecular investigations of the evolutionary progression among modern flamingos have been limited. Here, using three mitochondrial (mtDNA) markers, we present the first DNA sequence-based study of population genetic variation in the widely distributed Chilean Flamingo and, using two mtDNA and 10 nuclear (nDNA) markers, recover the species tree and divergence time estimates for the six extant species of flamingos. Phylogenetic analyses include likelihood and Bayesian frameworks and account for potential gene tree discordance. Analyses of divergence times are fossil calibrated at the divergence of Mirandornithes (flamingos + grebes) and the divergence of crown grebes. RESULTS: mtDNA sequences confirmed the presence of a single metapopulation represented by two minimally varying mtDNA barcodes in Chilean flamingos. Likelihood and Bayesian methods recovered identical phylogenies with flamingos falling into shallow-keeled (comprising the Greater, American and Chilean Flamingos) and deep-keeled (comprising the Lesser, Andean and James's Flamingos) sub-clades. The initial divergence among flamingos occurred at or shortly after the Mio-Pliocene boundary (6-3 Ma) followed by quick consecutive divergences throughout the Plio-Pleistocene. There is significant incongruence between the ages recovered by the mtDNA and nDNA datasets, likely due to mutational saturation occurring in the mtDNA loci. CONCLUSION: The finding of a single metapopulation in the widespread Chilean Flamingo confirms similar findings in other widespread flamingo species. The robust species phylogeny is congruent with previous classifications of flamingos based on feeding morphology. Modern phoenicopterids likely originated in the New World with each sub-clade dispersing across the Atlantic at least once. Our divergence time estimates place flamingos among the youngest families of birds, counter to the classical notion of flamingos as among the oldest based on biogeography and the fossil record. Finally, we designate 'Phoeniconaias' as a junior synonym of 'Phoenicoparrus' and redefine the latter genus as containing all flamingos more closely related to Phoenicoparrus andinus than Phoenicopterus roseus.

Pathophysiological processes underlying hidden hearing loss revealed in Kcnt1/2 double knockout mice.

Presbycusis is a prevalent condition in older adults characterized by the progressive loss of hearing due to age-related changes in the cochlea, the auditory portion of the inner ear. Many adults also struggle with understanding speech in noise despite having normal auditory thresholds, a condition termed "hidden" hearing loss because it evades standard audiological assessments. Examination of animal models and postmortem human tissue suggests that hidden hearing loss is also associated with age-related changes in the cochlea and may, therefore, precede overt age-related hearing loss. Nevertheless, the pathological mechanisms underlying hidden hearing loss are not understood, which hinders the development of diagnostic biomarkers and effective treatments for age-related hearing loss. To fill these gaps in knowledge, we leveraged a combination of tools, including transcriptomic profiling and morphological and functional assessments, to identify these processes and examine the transition from hidden to overt hearing loss. As a novel approach, we took advantage of a recently characterized model of hidden hearing loss: Kcnt1/2 double knockout mice. Using this model, we find that even before observable morphological pathology, hidden hearing loss is associated with significant alteration in several processes, notably proteostasis, in the cochlear sensorineural structures, and increased susceptibility to overt hearing loss in response to noise exposure and aging. Our findings provide the first insight into the pathophysiology associated with the earliest and, therefore, most treatable stages of hearing loss and provide critical insight directing future investigation of pharmaceutical strategies to slow and possibly prevent overt age-related hearing loss.

Effect of emphysema on AI software and human reader performance in lung nodule detection from low-dose chest CT.

BACKGROUND: Emphysema influences the appearance of lung tissue in computed tomography (CT). We evaluated whether this affects lung nodule detection by artificial intelligence (AI) and human readers (HR). METHODS: Individuals were selected from the "Lifelines" cohort who had undergone low-dose chest CT. Nodules in individuals without emphysema were matched to similar-sized nodules in individuals with at least moderate emphysema. AI results for nodular findings of 30-100 mm3 and 101-300 mm3 were compared to those of HR; two expert radiologists blindly reviewed discrepancies. Sensitivity and false positives (FPs)/scan were compared for emphysema and non-emphysema groups. RESULTS: Thirty-nine participants with and 82 without emphysema were included (n = 121, aged 61 ± 8 years (mean ± standard deviation), 58/121 males (47.9%)). AI and HR detected 196 and 206 nodular findings, respectively, yielding 109 concordant nodules and 184 discrepancies, including 118 true nodules. For AI, sensitivity was 0.68 (95% confidence interval 0.57-0.77) in emphysema versus 0.71 (0.62-0.78) in non-emphysema, with FPs/scan 0.51 and 0.22, respectively (p = 0.028). For HR, sensitivity was 0.76 (0.65-0.84) and 0.80 (0.72-0.86), with FPs/scan of 0.15 and 0.27 (p = 0.230). Overall sensitivity was slightly higher for HR than for AI, but this difference disappeared after the exclusion of benign lymph nodes. FPs/scan were higher for AI in emphysema than in non-emphysema (p = 0.028), while FPs/scan for HR were higher than AI for 30-100 mm3 nodules in non-emphysema (p = 0.009). CONCLUSIONS: AI resulted in more FPs/scan in emphysema compared to non-emphysema, a difference not observed for HR. RELEVANCE STATEMENT: In the creation of a benchmark dataset to validate AI software for lung nodule detection, the inclusion of emphysema cases is important due to the additional number of FPs. KEY POINTS: • The sensitivity of nodule detection by AI was similar in emphysema and non-emphysema. • AI had more FPs/scan in emphysema compared to non-emphysema. • Sensitivity and FPs/scan by the human reader were comparable for emphysema and non-emphysema. • Emphysema and non-emphysema representation in benchmark dataset is important for validating AI.

Functional, Morphological and Molecular Changes Reveal the Mechanisms Associated with Age-Related Vestibular Loss.

Age-related loss of vestibular function and hearing are common disorders that arise from the loss of function of the inner ear and significantly decrease quality of life. The underlying pathophysiological mechanisms are poorly understood and difficult to investigate in humans. Therefore, our study examined young (1.5-month-old) and old (24-month-old) C57BL/6 mice, utilizing physiological, histological, and transcriptomic methods. Vestibular sensory-evoked potentials revealed that older mice had reduced wave I amplitudes and delayed wave I latencies, indicating reduced vestibular function. Immunofluorescence and image analysis revealed that older mice exhibited a significant decline in type I sensory hair cell density, particularly in hair cells connected to dimorphic vestibular afferents. An analysis of gene expression in the isolated vestibule revealed the upregulation of immune-related genes and the downregulation of genes associated with ossification and nervous system development. A comparison with the isolated cochlear sensorineural structures showed similar changes in genes related to immune response, chondrocyte differentiation, and myelin formation. These findings suggest that age-related vestibular hypofunction is linked to diminished peripheral vestibular responses, likely due to the loss of a specific subpopulation of hair cells and calyceal afferents. The upregulation of immune- and inflammation-related genes implies that inflammation contributes to these functional and structural changes. Furthermore, the comparison of gene expression between the vestibule and cochlea indicates both shared and distinct mechanisms contributing to age-related vestibular and hearing impairments. Further research is necessary to understand the mechanistic connection between inflammation and age-related balance and hearing disorders and to translate these findings into clinical treatment strategies.

Transcriptome-Guided Identification of Drugs for Repurposing to Treat Age-Related Hearing Loss.

Age-related hearing loss (ARHL) or presbycusis is a prevalent condition associated with social isolation, cognitive impairment, and dementia. Age-related changes in the cochlea, the auditory portion of the inner ear, are the primary cause of ARHL. Unfortunately, there are currently no pharmaceutical approaches to treat ARHL. To examine the biological processes underlying age-related changes in the cochlea and identify candidate drugs for rapid repurposing to treat ARHL, we utilized bulk RNA sequencing to obtain transcriptomes from the functional substructures of the cochlea-the sensorineural structures, including the organ of Corti and spiral ganglion neurons (OC/SGN) and the stria vascularis and spiral ligament (SV/SL)-in young (6-week-old) and old (2-year-old) C57BL/6 mice. Transcriptomic analyses revealed both overlapping and unique patterns of gene expression and gene enrichment between substructures and with ageing. Based on these age-related transcriptional changes, we queried the protein products of genes differentially expressed with ageing in DrugBank and identified 27 FDA/EMA-approved drugs that are suitable to be repurposed to treat ARHL. These drugs target the protein products of genes that are differentially expressed with ageing uniquely in either the OC/SGN or SV/SL and that interrelate diverse biological processes. Further transcriptomic analyses revealed that most genes differentially expressed with ageing in both substructures encode protein products that are promising drug target candidates but are, nevertheless, not yet linked to approved drugs. Thus, with this study, we apply a novel approach to characterize the druggable genetic landscape for ARHL and propose a list of drugs to test in pre-clinical studies as potential treatment options for ARHL.

Synthesizing and databasing fossil calibrations: divergence dating and beyond.

Divergence dating studies, which combine temporal data from the fossil record with branch length data from molecular phylogenetic trees, represent a rapidly expanding approach to understanding the history of life. National Evolutionary Synthesis Center hosted the first Fossil Calibrations Working Group (3-6 March, 2011, Durham, NC, USA), bringing together palaeontologists, molecular evolutionists and bioinformatics experts to present perspectives from disciplines that generate, model and use fossil calibration data. Presentations and discussions focused on channels for interdisciplinary collaboration, best practices for justifying, reporting and using fossil calibrations and roadblocks to synthesis of palaeontological and molecular data. Bioinformatics solutions were proposed, with the primary objective being a new database for vetted fossil calibrations with linkages to existing resources, targeted for a 2012 launch.

Functional, Morphological, and Evolutionary Characterization of Hearing in Subterranean, Eusocial African Mole-Rats.

Naked mole-rats are highly vocal, eusocial, subterranean rodents with, counterintuitively, poor hearing. The causes underlying their altered hearing are unknown. Moreover, whether altered hearing is degenerate or adaptive to their unique lifestyles is controversial. We used various methods to identify the factors contributing to altered hearing in naked and the related Damaraland mole-rats and to examine whether these alterations result from relaxed or adaptive selection. Remarkably, we found that cochlear amplification was absent from both species despite normal prestin function in outer hair cells isolated from naked mole-rats. Instead, loss of cochlear amplification appears to result from abnormal hair bundle morphologies observed in both species. By exploiting a well-curated deafness phenotype-genotype database, we identified amino acid substitutions consistent with abnormal hair bundle morphology and reduced hearing sensitivity. Amino acid substitutions were found in unique groups of six hair bundle link proteins. Molecular evolutionary analyses revealed shifts in selection pressure at both the gene and the codon level for five of these six hair bundle link proteins. Substitutions in three of these proteins are associated exclusively with altered hearing. Altogether, our findings identify the likely mechanism of altered hearing in African mole-rats, making them the only identified mammals naturally lacking cochlear amplification. Moreover, our findings suggest that altered hearing in African mole-rats is adaptive, perhaps tailoring hearing to eusocial and subterranean lifestyles. Finally, our work reveals multiple, unique evolutionary trajectories in African mole-rat hearing and establishes species members as naturally occurring disease models to investigate human hearing loss.

Tempo and Pattern of Avian Brain Size Evolution.

Relative brain sizes in birds can rival those of primates, but large-scale patterns and drivers of avian brain evolution remain elusive. Here, we explore the evolution of the fundamental brain-body scaling relationship across the origin and evolution of birds. Using a comprehensive dataset sampling> 2,000 modern birds, fossil birds, and theropod dinosaurs, we infer patterns of brain-body co-variation in deep time. Our study confirms that no significant increase in relative brain size accompanied the trend toward miniaturization or evolution of flight during the theropod-bird transition. Critically, however, theropods and basal birds show weaker integration between brain size and body size, allowing for rapid changes in the brain-body relationship that set the stage for dramatic shifts in early crown birds. We infer that major shifts occurred rapidly in the aftermath of the Cretaceous-Paleogene mass extinction within Neoaves, in which multiple clades achieved higher relative brain sizes because of a reduction in body size. Parrots and corvids achieved the largest brains observed in birds via markedly different patterns. Parrots primarily reduced their body size, whereas corvids increased body and brain size simultaneously (with rates of brain size evolution outpacing rates of body size evolution). Collectively, these patterns suggest that an early adaptive radiation in brain size laid the foundation for subsequent selection and stabilization.

Source-sink dynamics structure a common montane mammal.

Assessing the relative role of evolutionary processes on genetic diversity is critical for understanding species response to climatic change. However, many processes, independent of climate, can lead to the same genetic pattern. Because effective population size and gene flow are affected directly by abundance and dispersal, population ecology has the potential to profoundly influence patterns of genetic variation over microevolutionary timescales. Here, we use aDNA data and simulations to explore the influence of population ecology and Holocene climate change on genetic diversity of the Uinta ground squirrel (Spermophilus armatus). We examined phylochronology from three modern and two ancient populations spanning the climate transitions of the last 3000 years. Population genetic analyses based on summary statistics suggest that changes in genetic diversity and structure coincided with the Medieval Warm Period (MWP), c. 1000 years ago. Serial coalescent simulations allowed us to move beyond correlation with climate to statistically compare the likelihoods of alternative population histories given the observed data. The data best fit source-sink models that include large, mid-elevation populations that exchange many migrants and small populations at the elevational extremes. While the MWP is likely to have reduced genetic diversity, our model-testing approach revealed that MWP-driven changes in genetic structure were not better supported for the range of models explored. Our results point to the importance of species ecology in understanding responses to climate, and showcase the use of ancient genetic data and simulation-based inference for unraveling the relative roles of microevolutionary processes.

Sodium-activated potassium channels shape peripheral auditory function and activity of the primary auditory neurons in mice.

Potassium (K+) channels shape the response properties of neurons. Although enormous progress has been made to characterize K+ channels in the primary auditory neurons, the molecular identities of many of these channels and their contributions to hearing in vivo remain unknown. Using a combination of RNA sequencing and single molecule fluorescent in situ hybridization, we localized expression of transcripts encoding the sodium-activated potassium channels KNa1.1 (SLO2.2/Slack) and KNa1.2 (SLO2.1/Slick) to the primary auditory neurons (spiral ganglion neurons, SGNs). To examine the contribution of these channels to function of the SGNs in vivo, we measured auditory brainstem responses in KNa1.1/1.2 double knockout (DKO) mice. Although auditory brainstem response (wave I) thresholds were not altered, the amplitudes of suprathreshold responses were reduced in DKO mice. This reduction in amplitude occurred despite normal numbers and molecular architecture of the SGNs and their synapses with the inner hair cells. Patch clamp electrophysiology of SGNs isolated from DKO mice displayed altered membrane properties, including reduced action potential thresholds and amplitudes. These findings show that KNa1 channel activity is essential for normal cochlear function and suggest that early forms of hearing loss may result from physiological changes in the activity of the primary auditory neurons.

Fire and ice: genetic structure of the Uinta ground squirrel (Spermophilus armatus) across the Yellowstone hotspot.

The range of the Uinta ground squirrel, Spermophilus armatus, is centred over one of the most tectonically active regions today, the Yellowstone hotspot. We document the role of Quaternary tectonic and climatic history on the genetic structure of this species by screening museum and extant individuals throughout its range. Phylogeographic, divergence time, and demographic analyses of partial mitochondrial cytochrome b and control region DNA sequences yield insight into the cadence of evolution across three spatiotemporal scales: (i) a relatively deep intraspecific divergence of S. armatus into three lineages coincident with the last major volcanic eruption in the region and maintained by the Snake River Plain; (ii) demographic expansion in two lineages corresponding to the time of last deglaciation of the region; and (iii) a recent (< 50 years) local extinction of the third lineage coincident with climatic change and conversion of habitat for agricultural purposes in eastern Idaho. Beyond these inferences, our study highlights the unique value of museum material to phylogeography, and shows that small mammal recolonization of previously glaciated montane 'islands' differs from northward postglacial expansion observed in areas previously covered by continental ice sheets. Montane 'islands' may harbour high genetic diversity because of admixture and recurrent expansion/extinction.

Whole-genome resequencing reveals signatures of selection and timing of duck domestication.

Background: The genetic basis of animal domestication remains poorly understood, and systems with substantial phenotypic differences between wild and domestic populations are useful for elucidating the genetic basis of adaptation to new environments as well as the genetic basis of rapid phenotypic change. Here, we sequenced the whole genome of 78 individual ducks, from two wild and seven domesticated populations, with an average sequencing depth of 6.42X per individual. Results: Our population and demographic analyses indicate a complex history of domestication, with early selection for separate meat and egg lineages. Genomic comparison of wild to domesticated populations suggests that genes that affect brain and neuronal development have undergone strong positive selection during domestication. Our FST analysis also indicates that the duck white plumage is the result of selection at the melanogenesis-associated transcription factor locus. Conclusions: Our results advance the understanding of animal domestication and selection for complex phenotypic traits.

Genetic response to climatic change: insights from ancient DNA and phylochronology.

Understanding how climatic change impacts biological diversity is critical to conservation. Yet despite demonstrated effects of climatic perturbation on geographic ranges and population persistence, surprisingly little is known of the genetic response of species. Even less is known over ecologically long time scales pertinent to understanding the interplay between microevolution and environmental change. Here, we present a study of population variation by directly tracking genetic change and population size in two geographically widespread mammal species (Microtus montanus and Thomomys talpoides) during late-Holocene climatic change. We use ancient DNA to compare two independent estimates of population size (ecological and genetic) and corroborate our results with gene diversity and serial coalescent simulations. Our data and analyses indicate that, with population size decreasing at times of climatic change, some species will exhibit declining gene diversity as expected from simple population genetic models, whereas others will not. While our results could be consistent with selection, independent lines of evidence implicate differences in gene flow, which depends on the life history strategy of species.

Potential for bias and low precision in molecular divergence time estimation of the Canopy of Life: an example from aquatic bird families.

Uncertainty in divergence time estimation is frequently studied from many angles but rarely from the perspective of phylogenetic node age. If appropriate molecular models and fossil priors are used, a multi-locus, partitioned analysis is expected to equally minimize error in accuracy and precision across all nodes of a given phylogeny. In contrast, if available models fail to completely account for rate heterogeneity, substitution saturation and incompleteness of the fossil record, uncertainty in divergence time estimation may increase with node age. While many studies have stressed this concern with regard to deep nodes in the Tree of Life, the inference that molecular divergence time estimation of shallow nodes is less sensitive to erroneous model choice has not been tested explicitly in a Bayesian framework. Because of available divergence time estimation methods that permit fossil priors across any phylogenetic node and the present increase in efficient, cheap collection of species-level genomic data, insight is needed into the performance of divergence time estimation of shallow (<10 MY) nodes. Here, we performed multiple sensitivity analyses in a multi-locus data set of aquatic birds with six fossil constraints. Comparison across divergence time analyses that varied taxon and locus sampling, number and position of fossil constraint and shape of prior distribution showed various insights. Deviation from node ages obtained from a reference analysis was generally highest for the shallowest nodes but determined more by temporal placement than number of fossil constraints. Calibration with only the shallowest nodes significantly underestimated the aquatic bird fossil record, indicating the presence of saturation. Although joint calibration with all six priors yielded ages most consistent with the fossil record, ages of shallow nodes were overestimated. This bias was found in both mtDNA and nDNA regions. Thus, divergence time estimation of shallow nodes may suffer from bias and low precision, even when appropriate fossil priors and best available substitution models are chosen. Much care must be taken to address the possible ramifications of substitution saturation across the entire Tree of Life.

Genome skimming identifies polymorphism in tern populations and species.

BACKGROUND: Terns (Charadriiformes: Sterninae) are a lineage of cosmopolitan shorebirds with a disputed evolutionary history that comprises several species of conservation concern. As a non-model system in genetics, previous study has left most of the nuclear genome unexplored, and population-level studies are limited to only 15% of the world's species of terns and noddies. Screening of polymorphic nuclear sequence markers is needed to enhance genetic resolution because of supposed low mitochondrial mutation rate, documentation of nuclear insertion of hypervariable mitochondrial regions, and limited success of microsatellite enrichment in terns. Here, we investigated the phylogenetic and population genetic utility for terns and relatives of a variety of nuclear markers previously developed for other birds and spanning the nuclear genome. Markers displaying a variety of mutation rates from both the nuclear and mitochondrial genome were tested and prioritized according to optimal cross-species amplification and extent of genetic polymorphism between (1) the main tern clades and (2) individual Royal Terns (Thalasseus maxima) breeding on the US East Coast. RESULTS: Results from this genome skimming effort yielded four new nuclear sequence-based markers for tern phylogenetics and 11 intra-specific polymorphic markers. Further, comparison between the two genomes indicated a phylogenetic conflict at the base of terns, involving the inclusion (mitochondrial) or exclusion (nuclear) of the Angel Tern (Gygis alba). Although limited mitochondrial variation was confirmed, both nuclear markers and a short tandem repeat in the mitochondrial control region indicated the presence of considerable genetic variation in Royal Terns at a regional scale. CONCLUSIONS: These data document the value of intronic markers to the study of terns and allies. We expect that these and additional markers attained through next-generation sequencing methods will accurately map the genetic origin and species history of this group of birds.

Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds.

Using standard phylogenetic methods, it can be hard to resolve the order in which speciation events took place when new lineages evolved in the distant past and within a short time frame. As an example, phylogenies of galliform birds (including well-known species such as chicken, turkey, and quail) usually show low bootstrap support values at short internal branches, reflecting the rapid diversification of these birds in the Eocene. However, given the key role of chicken and related poultry species in agricultural, evolutionary, general biological and disease studies, it is important to know their internal relationships. Recently, insertion patterns of transposable elements such as long and short interspersed nuclear element markers have proved powerful in revealing branching orders of difficult phylogenies. Here we decipher the order of speciation events in a group of 27 galliform species based on insertion events of chicken repeat 1 (CR1) transposable elements. Forty-four CR1 marker loci were identified from the draft sequence of the chicken genome, and from turkey BAC clone sequence, and the presence or absence of markers across species was investigated via electrophoretic size separation of amplification products and subsequent confirmation by DNA sequencing. Thirty markers proved possible to type with electrophoresis of which 20 were phylogenetically informative. The distribution of these repeat elements supported a single homoplasy-free cladogram, which confirmed that megapodes, cracids, New World quail, and guinea fowl form outgroups to Phasianidae and that quails, pheasants, and partridges are each polyphyletic groups. Importantly, we show that chicken is an outgroup to turkey and quail, an observation which does not have significant support from previous DNA sequence- and DNA-DNA hybridization-based trees and has important implications for evolutionary studies based on sequence or karyotype data from galliforms. We discuss the potential and limitations of using a genome-based retrotransposon approach in resolving problematic phylogenies among birds.

The shifting baseline of northern fur seal ecology in the northeast Pacific Ocean.

Historical data provide a baseline against which to judge the significance of recent ecological shifts and guide conservation strategies, especially for species decimated by pre-20th century harvesting. Northern fur seals (NFS; Callorhinus ursinus) are a common pinniped species in archaeological sites from southern California to the Aleutian Islands, yet today they breed almost exclusively on offshore islands at high latitudes. Harvest profiles from archaeological sites contain many unweaned pups, confirming the presence of temperate-latitude breeding colonies in California, the Pacific Northwest, and the eastern Aleutian Islands. Isotopic results suggest that prehistoric NFS fed offshore across their entire range, that California populations were distinct from populations to the north, and that populations breeding at temperate latitudes in the past used a different reproductive strategy than modern populations. The extinction of temperate-latitude breeding populations was asynchronous geographically. In southern California, the Pacific Northwest, and the eastern Aleutians, NFS remained abundant in the archaeological record up to the historical period approximately 200 years B.P.; thus their regional collapse is plausibly attributed to historical hunting or some other anthropogenic ecosystem disturbance. In contrast, NFS populations in central and northern California collapsed at approximately 800 years B.P., long before European contact. The relative roles of human hunting versus climatic factors in explaining this ecological shift are unclear, as more paleoclimate information is needed from the coastal zone.