Museum Genomics.

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

Boom-bust population dynamics drive rapid genetic change.

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.

Genomes, fossils, and the concurrent rise of modern birds and flowering plants in the Late Cretaceous.

The phylogeny and divergence timing of the Neoavian radiation remain controversial despite recent progress. We analyzed the genomes of 124 species across all Neoavian orders, using data from 25,460 loci spanning four DNA classes, including 5,756 coding sequences, 12,449 conserved nonexonic elements, 4,871 introns, and 2,384 intergenic segments. We conducted a comprehensive sensitivity analysis to account for the heterogeneity across different DNA classes, leading to an optimal tree of Neoaves with high resolution. This phylogeny features a novel Neoavian dichotomy comprising two monophyletic clades: a previously recognized Telluraves (land birds) and a newly circumscribed Aquaterraves (waterbirds and relatives). Molecular dating analyses with 20 fossil calibrations indicate that the diversification of modern birds began in the Late Cretaceous and underwent a constant and steady radiation across the KPg boundary, concurrent with the rise of angiosperms as well as other major Cenozoic animal groups including placental and multituberculate mammals. The KPg catastrophe had a limited impact on avian evolution compared to the Paleocene-Eocene Thermal Maximum, which triggered a rapid diversification of seabirds. Our findings suggest that the evolution of modern birds followed a slow process of gradualism rather than a rapid process of punctuated equilibrium, with limited interruption by the KPg catastrophe. This study places bird evolution into a new context within vertebrates, with ramifications for the evolution of the Earth's biota.

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.

Historic samples reveal loss of wild genotype through domestic chicken introgression during the Anthropocene.

Human activities have precipitated a rise in the levels of introgressive gene flow among animals. The investigation of conspecific populations at different time points may shed light on the magnitude of human-mediated introgression. We used the red junglefowl Gallus gallus, the wild ancestral form of the chicken, as our study system. As wild junglefowl and domestic chickens readily admix, conservationists fear that domestic introgression into junglefowl may compromise their wild genotype. By contrasting the whole genomes of 51 chickens with 63 junglefowl from across their natural range, we found evidence of a loss of the wild genotype across the Anthropocene. When comparing against the genomes of junglefowl from approximately a century ago using rigorous ancient-DNA protocols, we discovered that levels of domestic introgression are not equal among and within modern wild populations, with the percentage of domestic ancestry around 20-50%. We identified a number of domestication markers in which chickens are deeply differentiated from historic junglefowl regardless of breed and/or geographic provenance, with eight genes under selection. The latter are involved in pathways dealing with development, reproduction and vision. The wild genotype is an allelic reservoir that holds most of the genetic diversity of G. gallus, a species which is immensely important to human society. Our study provides fundamental genomic infrastructure to assist in efforts to prevent a further loss of the wild genotype through introgression of domestic alleles.

Convergence in hearing-related genes between echolocating birds and mammals.

Echolocation, the detection of objects by means of sound waves, has evolved independently in diverse animals. Echolocators include not only mammals such as toothed whales and yangochiropteran and rhinolophoid bats but also Rousettus fruit bats, as well as two bird lineages, oilbirds and swiftlets. In whales and yangochiropteran and rhinolophoid bats, positive selection and molecular convergence has been documented in key hearing-related genes, such as prestin (SLC26A5), but few studies have examined these loci in other echolocators. Here, we examine patterns of selection and convergence in echolocation-related genes in echolocating birds and Rousettus bats. Fewer of these loci were under selection in Rousettus or birds compared with classically recognized echolocators, and elevated convergence (compared to outgroups) was not evident across this gene set. In certain genes, however, we detected convergent substitutions with potential functional relevance, including convergence between Rousettus and classic echolocators in prestin at a site known to affect hair cell electromotility. We also detected convergence between Yangochiroptera, Rhinolophidea, and oilbirds in TMC1, an important mechanosensory transduction channel in vertebrate hair cells, and observed an amino acid change at the same site within the pore domain. Our results suggest that although most proteins implicated in echolocation in specialized mammals may not have been recruited in birds or Rousettus fruit bats, certain hearing-related loci may have undergone convergent functional changes. Investigating adaptations in diverse echolocators will deepen our understanding of this unusual sensory modality.

Tadpoles rely on mechanosensory stimuli for communication when visual capabilities are poor.

The ways in which animals sense the world changes throughout development. For example, young of many species have limited visual capabilities, but still make social decisions, likely based on information gathered through other sensory modalities. Poison frog tadpoles display complex social behaviors that have been suggested to rely on vision despite a century of research indicating tadpoles have poorly-developed visual systems relative to adults. Alternatively, other sensory modalities, such as the lateral line system, are functional at hatching in frogs and may guide social decisions while other sensory systems mature. Here, we examined development of the mechanosensory lateral line and visual systems in tadpoles of the mimic poison frog (Ranitomeya imitator) that use vibrational begging displays to stimulate egg feeding from their mothers. We found that tadpoles hatch with a fully developed lateral line system. While begging behavior increases with development, ablating the lateral line system inhibited begging in pre-metamorphic tadpoles, but not in metamorphic tadpoles. We also found that the increase in begging and decrease in reliance on the lateral line co-occurs with increased retinal neural activity and gene expression associated with eye development. Using the neural tracer neurobiotin, we found that axonal innervations from the eye to the brain proliferate during metamorphosis, with few retinotectal connections in recently-hatched tadpoles. We then tested visual function in a phototaxis assay and found tadpoles prefer darker environments. The strength of this preference increased with developmental stage, but eyes were not required for this behavior, possibly indicating a role for the pineal gland. Together, these data suggest that tadpoles rely on different sensory modalities for social interactions across development and that the development of sensory systems in socially complex poison frog tadpoles is similar to that of other frog species.

Genetic Basis and Evolution of Structural Color Polymorphism in an Australian Songbird.

Island organisms often evolve phenotypes divergent from their mainland counterparts, providing a useful system for studying adaptation under differential selection. In the white-winged fairywren (Malurus leucopterus), subspecies on two islands have a black nuptial plumage whereas the subspecies on the Australian mainland has a blue nuptial plumage. The black subspecies have a feather nanostructure that could in principle produce a blue structural color, suggesting a blue ancestor. An earlier study proposed independent evolution of melanism on the islands based on the history of subspecies divergence. However, the genetic basis of melanism and the origin of color differentiation in this group are still unknown. Here, we used whole-genome resequencing to investigate the genetic basis of melanism by comparing the blue and black M. leucopterus subspecies to identify highly divergent genomic regions. We identified a well-known pigmentation gene ASIP and four candidate genes that may contribute to feather nanostructure development. Contrary to the prediction of convergent evolution of island melanism, we detected signatures of a selective sweep in genomic regions containing ASIP and SCUBE2 not in the black subspecies but in the blue subspecies, which possesses many derived SNPs in these regions, suggesting that the mainland subspecies has re-evolved a blue plumage from a black ancestor. This proposed re-evolution was likely driven by a preexisting female preference. Our findings provide new insight into the evolution of plumage coloration in island versus continental populations, and, importantly, we identify candidate genes that likely play roles in the development and evolution of feather structural coloration.

A phylogenetic method linking nucleotide substitution rates to rates of continuous trait evolution.

Genomes contain conserved non-coding sequences that perform important biological functions, such as gene regulation. We present a phylogenetic method, PhyloAcc-C, that associates nucleotide substitution rates with changes in a continuous trait of interest. The method takes as input a multiple sequence alignment of conserved elements, continuous trait data observed in extant species, and a background phylogeny and substitution process. Gibbs sampling is used to assign rate categories (background, conserved, accelerated) to lineages and explore whether the assigned rate categories are associated with increases or decreases in the rate of trait evolution. We test our method using simulations and then illustrate its application using mammalian body size and lifespan data previously analyzed with respect to protein coding genes. Like other studies, we find processes such as tumor suppression, telomere maintenance, and p53 regulation to be related to changes in longevity and body size. In addition, we also find that skeletal genes, and developmental processes, such as sprouting angiogenesis, are relevant.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae).

BACKGROUND: Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary. RESULTS: Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression. CONCLUSIONS: Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

PhyloAcc-GT: A Bayesian Method for Inferring Patterns of Substitution Rate Shifts on Targeted Lineages Accounting for Gene Tree Discordance.

An important goal of evolutionary genomics is to identify genomic regions whose substitution rates differ among lineages. For example, genomic regions experiencing accelerated molecular evolution in some lineages may provide insight into links between genotype and phenotype. Several comparative genomics methods have been developed to identify genomic accelerations between species, including a Bayesian method called PhyloAcc, which models shifts in substitution rate in multiple target lineages on a phylogeny. However, few methods consider the possibility of discordance between the trees of individual loci and the species tree due to incomplete lineage sorting, which might cause false positives. Here, we present PhyloAcc-GT, which extends PhyloAcc by modeling gene tree heterogeneity. Given a species tree, we adopt the multispecies coalescent model as the prior distribution of gene trees, use Markov chain Monte Carlo (MCMC) for inference, and design novel MCMC moves to sample gene trees efficiently. Through extensive simulations, we show that PhyloAcc-GT outperforms PhyloAcc and other methods in identifying target lineage-specific accelerations and detecting complex patterns of rate shifts, and is robust to specification of population size parameters. PhyloAcc-GT is usually more conservative than PhyloAcc in calling convergent rate shifts because it identifies more accelerations on ancestral than on terminal branches. We apply PhyloAcc-GT to two examples of convergent evolution: flightlessness in ratites and marine mammal adaptations, and show that PhyloAcc-GT is a robust tool to identify shifts in substitution rate associated with specific target lineages while accounting for incomplete lineage sorting.

Whole-genome Analyses Reveal Past Population Fluctuations and Low Genetic Diversities of the North Pacific Albatrosses.

Throughout the Plio-Pleistocene, climate change has impacted tropical marine ecosystems substantially, with even more severe impacts predicted in the Anthropocene. Although many studies have clarified demographic histories of seabirds in polar regions, the history of keystone seabirds of the tropics is unclear, despite the prominence of albatrosses (Diomedeidae, Procellariiformes) as the largest and most threatened group of oceanic seabirds. To understand the impact of climate change on tropical albatrosses, we investigated the evolutionary and demographic histories of all four North Pacific albatrosses and their prey using whole-genome analyses. We report a striking concordance in demographic histories among the four species, with a notable dip in effective population size at the beginning of the Pleistocene and a population expansion in the Last Glacial Period when sea levels were low, which resulted in increased potential coastal breeding sites. Abundance of the black-footed albatross dropped again during the Last Glacial Maximum, potentially linked to climate-driven loss of breeding sites and concordant genome-derived decreases in its major prey. We find very low genome-wide (π < 0.001) and adaptative genetic diversities across the albatrosses, with genes of the major histocompatibility complex close to monomorphic. We also identify recent selective sweeps at genes associated with hyperosmotic adaptation, longevity, and cognition and memory. Our study has shed light on the evolutionary and demographic histories of the largest tropical oceanic seabirds and provides evidence for their large population fluctuations and alarmingly low genetic diversities.

Testosterone Coordinates Gene Expression Across Different Tissues to Produce Carotenoid-Based Red Ornamentation.

Carotenoid pigments underlie most of the red, orange, and yellow visual signals used in mate choice in vertebrates. However, many of the underlying processes surrounding the production of carotenoid-based traits remain unclear due to the complex nature of carotenoid uptake, metabolism, and deposition across tissues. Here, we leverage the ability to experimentally induce the production of a carotenoid-based red plumage patch in the red-backed fairywren (Malurus melanocephalus), a songbird in which red plumage is an important male sexual signal. We experimentally elevated testosterone in unornamented males lacking red plumage to induce the production of ornamentation and compared gene expression in both the liver and feather follicles between unornamented control males, testosterone-implanted males, and naturally ornamented males. We show that testosterone upregulates the expression of CYP2J19, a gene known to be involved in ketocarotenoid metabolism, and a putative carotenoid processing gene (ELOVL6) in the liver, and also regulates the expression of putative carotenoid transporter genes in red feather follicles on the back, including ABCG1. In black feathers, carotenoid-related genes are downregulated and melanin genes upregulated, but we find that carotenoids are still present in the feathers. This may be due to the activity of the carotenoid-cleaving enzyme BCO2 in black feathers. Our study provides a first working model of a pathway for carotenoid-based trait production in free-living birds, implicates testosterone as a key regulator of carotenoid-associated gene expression, and suggests hormones may coordinate the many processes that underlie the production of these traits across multiple tissues.

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum.

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.

Evolution of the DAN gene family in vertebrates.

The DAN gene family (DAN, Differential screening-selected gene Aberrant in Neuroblastoma) is a group of genes that is expressed during development and plays fundamental roles in limb bud formation and digitation, kidney formation and morphogenesis and left-right axis specification. During adulthood the expression of these genes are associated with diseases, including cancer. Although most of the attention to this group of genes has been dedicated to understanding its role in physiology and development, its evolutionary history remains poorly understood. Thus, the goal of this study is to investigate the evolutionary history of the DAN gene family in vertebrates, with the objective of complementing the already abundant physiological information with an evolutionary context. Our results recovered the monophyly of all DAN gene family members and divide them into five main groups. In addition to the well-known DAN genes, our phylogenetic results revealed the presence of two new DAN gene lineages; one is only retained in cephalochordates, whereas the other one (GREM3) was only identified in cartilaginous fish, holostean fish, and coelacanth. According to the phyletic distribution of the genes, the ancestor of gnathostomes possessed a repertoire of eight DAN genes, and during the radiation of the group GREM1, GREM2, SOST, SOSTDC1, and NBL1 were retained in all major groups, whereas, GREM3, CER1, and DAND5 were differentially lost.

Whole-genome phylogeography of the blue-faced honeyeater (Entomyzon cyanotis) and discovery and characterization of a neo-Z chromosome.

Whole-genome surveys of genetic diversity and geographic variation often yield unexpected discoveries of novel structural variation, which long-read DNA sequencing can help clarify. Here, we report on whole-genome phylogeography of a bird exhibiting classic vicariant geographies across Australia and New Guinea, the blue-faced honeyeater (Entomyzon cyanotis), and the discovery and characterization of a novel neo-Z chromosome by long-read sequencing. Using short-read genome-wide SNPs, we inferred population divergence events within E. cyanotis across the Carpentarian and other biogeographic barriers during the Pleistocene (~0.3-1.7 Ma). Evidence for introgression between nonsister populations supports a hypothesis of reticulate evolution around a triad of dynamic barriers around Pleistocene Lake Carpentaria between Australia and New Guinea. During this phylogeographic survey, we discovered a large (134 Mbp) neo-Z chromosome and we explored its diversity, divergence and introgression landscape. We show that, as in some sylvioid passerine birds, a fusion occurred between chromosome 5 and the Z chromosome to form a neo-Z chromosome; and in E. cyanotis, the ancestral pseudoautosomal region (PAR) appears nonrecombinant between Z and W, along with most of the fused chromosome 5. The added recombination-suppressed portion of the neo-Z (~37.2 Mbp) displays reduced diversity and faster population genetic differentiation compared with the ancestral-Z. Yet, the new PAR (~17.4 Mbp) shows elevated diversity and reduced differentiation compared to autosomes, potentially resulting from introgression. In our case, long-read sequencing helped clarify the genomic landscape of population divergence on autosomes and sex chromosomes in a species where prior knowledge of genome structure was still incomplete.

Episodic population fragmentation and gene flow reveal a trade-off between heterozygosity and allelic richness.

In episodic environments like deserts, populations of some animal species exhibit irregular fluctuations such that populations are alternately large and connected or small and isolated. Such dynamics are typically driven by periodic resource pulses due, for example, to large but infrequent rainfall events. The repeated population bottlenecks resulting from fragmentation should lower genetic diversity over time, yet species undergoing these fluctuations appear to maintain high levels of genetic diversity. To resolve this apparent paradox, we simulated a metapopulation of constant size undergoing repeat episodes of fragmentation and change in gene flow to mimic outcomes experienced by mammals in an Australian desert. We show that episodic fragmentation and gene flow have contrasting effects on two measures of genetic diversity: heterozygosity and allelic richness. Specifically, fragmentation into many, small subpopulations, coupled with periods of infrequent gene flow, preserves allelic richness at the expense of heterozygosity. In contrast, fragmentation into a few, large subpopulations maintains heterozygosity at the expense of allelic richness. The strength of the trade-off between heterozygosity and allelic richness depends on the amount of gene flow and the frequency of gene flow events. Our results imply that the type of genetic diversity maintained among species living in strongly fluctuating environments will depend on the way populations fragment, with our results highlighting different mechanisms for maintaining allelic richness and heterozygosity in small, fragmented populations.

A Unifying Framework for Estimating Generation Time in Age-Structured Populations: Implications for Phylogenetics and Conservation Biology.

AbstractGeneration time is a measure of the pace of life and is used to describe processes in population dynamics and evolution. We show that three commonly used mathematical definitions of generation time in age-structured populations can produce different estimates of up to several years for the same set of life history data. We present and prove a mathematical theorem that reveals a general order relation among the definitions. Furthermore, the exact population growth rate at the time of sampling influences estimates of generation time, which calls for attention. For phylogenetic estimates of divergence times between species, included demographic data should be collected when the population growth rate for each species is most common and typical. In conservation biology, demographic data should be collected during phases of population decline in declining species, contrary to common recommendations to use predisturbance data. The results can be used to improve the International Union for Conservation of Nature's recommendation in parameterizing models for evaluating threat categories of threatened species and to avoid underestimating extinction risk.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population.

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.