Genetic depletion at adaptive but not neutral loci in an endangered bird species

Differences between mates at the TLR1Lb locus are associated with lower reproductive success in a long-lived seabird

Genetic diversity can influence fitness components such as survival and reproductive success. Yet the association between genetic diversity and fitness based on neutral loci is sometime very weak and inconsistent, with relationships varying among taxa due to confounding effects of population demography and life history. Fitness-diversity relationships are likely to be stronger and more consistent for genes known to influence phenotypic traits, such as immunity-related genes, and may also depend on the genetic differences between breeding partners. We recorded breeding success of individuals and breeding pairs over 20 years to evaluate the relationships between reproductive success and both neutral genetic variation (using 7,830 single nucleotide polymorphisms) and functional variation (four toll-like receptor [TLRs] loci) with reproductive success in thick-billed murres (Uria lomvia). Individual genetic diversity (both neutral and functional) was unrelated to reproductive success, but surprisingly, successful multi-year reproductive success decreased with the genetic difference between breeding partners at TLR1Lb. This result may be due to an advantage of specific alleles at TLR1Lb. This study is one of few addressing both individual genetic variation and genetic similarity between mates at both neutral and functional variation in a long-lived bird.

Evaluating genotyping-in-thousands by sequencing as a genetic monitoring tool for a climate sentinel mammal using non-invasive and archival samples

Genetic tools for wildlife monitoring can provide valuable information on spatiotemporal population trends and connectivity, particularly in systems experiencing rapid environmental change. Multiplexed targeted amplicon sequencing techniques, such as genotyping-in-thousands by sequencing (GT-seq), can provide cost-effective approaches for collecting genetic data from low-quality and quantity DNA samples, making them potentially useful for long-term wildlife monitoring using non-invasive and archival samples. Here, we developed a GT-seq panel as a potential monitoring tool for the American pika (Ochotona princeps) and evaluated its performance when applied to traditional, non-invasive, and archival samples, respectively. Specifically, we optimized a GT-seq panel (307 single nucleotide polymorphisms (SNPs)) that included neutral, sex-associated, and putatively adaptive SNPs using contemporary tissue samples (n = 77) from the Northern Rocky Mountains lineage of American pikas. The panel demonstrated high genotyping success (94.7%), low genotyping error (0.001%), and excellent performance identifying individuals, sex, relatedness, and population structure. We subsequently applied the GT-seq panel to archival tissue (n = 17) and contemporary fecal pellet samples (n = 129) collected within the Canadian Rocky Mountains to evaluate its effectiveness. Although the panel demonstrated high efficacy with archival tissue samples (90.5% genotyping success, 0.0% genotyping error), this was not the case for the fecal pellet samples (79.7% genotyping success, 28.4% genotyping error) likely due to the exceptionally low quality/quantity of recovered DNA using the approaches implemented. Overall, our study reinforced GT-seq as an effective tool using contemporary and archival tissue samples, providing future opportunities for temporal applications using historical specimens. Our results further highlight the need for additional optimization of sample and genetic data collection techniques prior to broader-scale implementation of a non-invasive genetic monitoring tool for American pikas.

Urbanization is associated with non-coding polymorphisms in candidate behavioural genes in the Eurasian coot

Extensive transformation of natural land cover into urbanized areas enhances accumulation of phenotypic differences between animals from urban and nonurban populations, but there is little information on whether these changes, especially in terms of animal behaviour and circadian rhythm, have a genetic basis. The aim of this study was to investigate genetic background of behavioural differences between four pairs of urban and nonurban populations of a common waterbird, the Eurasian coot Fulica atra. For this purpose, we quantified polymorphisms in personality-related candidate genes, previously reported to be associated with avian circadian rhythms and behavioural traits that may be crucial for urban life. We found general associations between landscape urbanization level and polymorphisms in 3'UTR region of CREB1 gene encoding transcriptional factor, which participates in development of cognitive functions and regulation of circadian rhythm. We also found significant differentiation between urban and nonurban populations in the intronic region of CKIɛ gene responsible for regulation of circadian clock. Although we lacked evidence for linkage of this intronic variation with coding polymorphisms, genetic differentiation between urban populations was significantly stronger at CKIɛ intron compared with neutral microsatellite markers, suggesting possible local adaptations of CKIɛ expression regulation to specific urban sites. Our results indicate that behavioural differentiation between urban and nonurban coot populations may be the effect of habitat-specific selective pressure resulting in genetic adaptations to urban environment and supporting the microevolutionary scenario. These adaptations, however, prevailed in non-coding regulatory rather than coding gene regions and showed either general or local patterns, revealing high complexity of associations between behaviour and landscape urbanization in birds.

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

Immune genotypes, immune responses, and survival in a wild bird population

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

Leveraging an existing whole-genome resequencing population data set to characterize toll-like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll-like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species-specific primers and PCR. Here, we leverage an existing short-read reference genome, whole-genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat'/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat'/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Using temporal genomics to understand contemporary climate change responses in wildlife

Monitoring the evolutionary responses of species to ongoing global climate change is critical for informing conservation. Population genomic studies that use samples from multiple time points ("temporal genomics") are uniquely able to make direct observations of change over time. Consequently, only temporal studies can show genetic erosion or spatiotemporal changes in population structure. Temporal genomic studies directly examining climate change effects are currently rare but will likely increase in the coming years due to their high conservation value. Here, we highlight four key genetic indicators that can be monitored using temporal genomics to understand how species are responding to climate change. All indicators crucially rely on having a suitable baseline that accurately represents the past condition of the population, and we discuss aspects of study design that must be considered to achieve this.

Restoring faith in conservation action: Maintaining wild genetic diversity through the Tasmanian devil insurance program

Conservation breeding programs aim to maintain 90% wild genetic diversity, but rarely assess functional diversity. Here, we compare both genome-wide and functional diversity (in over 500 genes) of Tasmanian devils (Sarcophilus harrisii) within the insurance metapopulation and across the species’ range (64,519 km²). Populations have declined by 80% since 1996 due to a contagious cancer, devil facial tumor disease (DFTD). However, predicted local extinctions have not occurred. Recent suggestions of selection for “resistance” alleles in the wild precipitated concerns that insurance population devils may be unsuitable for translocations. Using 830 wild samples collected at 31 locations between 2012 and 2021, and 553 insurance metapopulation devils, we show that the insurance metapopulation is representative of current wild genetic diversity. Allele frequencies at DFTD-associated loci were not substantially different between captive and wild devils. Methods presented here are valuable for others investigating evolutionary potential in threatened species, particularly ones under significant selective pressures.

•

Developed target capture to assess functional diversity at over 500 genes

•

Fine-scale structure exists in the genetically depauperate Tasmanian devil

•

Insurance metapopulation is representative of wild genetic diversity

•

Allele frequencies at disease-associated loci were similar in captivity to the wild

Selection balancing at innate immune genes: adaptive polymorphism maintenance in Toll-like receptors

Balancing selection is a classic mechanism for maintaining variability in immune genes involved in host–pathogen interactions. However, it remains unclear how widespread the mechanism is across immune genes other than the major histocompatibility complex (MHC). Although occasional reports suggest that balancing selection (heterozygote advantage, negative frequency-dependent selection, and fluctuating selection) may act on other immune genes, the current understanding of the phenomenon in non-MHC immune genes is far from solid. In this review, we focus on Toll-like receptors (TLRs), innate immune genes directly involved in pathogen recognition and immune response activation, as there is a growing body of research testing the assumptions of balancing selection in these genes. After reviewing infection- and fitness-based evidence, along with evidence based on population allelic frequencies and heterozygosity levels, we conclude that balancing selection maintains variation in TLRs, though it tends to occur under specific conditions in certain evolutionary lineages rather than being universal and ubiquitous. Our review also identifies key gaps in current knowledge and proposes promising areas for future research. Improving our understanding of host–pathogen interactions and balancing selection in innate immune genes are increasingly important, particularly regarding threats from emerging zoonotic diseases.

Contemporary evolution of the innate immune receptor gene TLR3 in an isolated vertebrate population

Understanding where genetic variation exists, and how it influences fitness within populations is important from an evolutionary and conservation perspective. Signatures of past selection suggest that pathogen-mediated balancing selection is a key driver of immunogenetic variation, but studies tracking contemporary evolution are needed to help resolve the evolutionary forces and mechanism at play. Previous work in a bottlenecked population of Seychelles warblers (Acrocephalus sechellensis) show that functional variation has been maintained at the viral-sensing Toll-like receptor 3 (TLR3) gene, including one nonsynonymous SNP, resulting in two alleles. Here, we characterise evolution at this TLR3 locus over a 25-year period within the original remnant population of the Seychelles warbler, and in four other derived, populations. Results show a significant and consistent temporal decline in the frequency of the TLR3^C allele in the original population, and that similar declines in the TLR3^C allele frequency occurred in all the derived populations. Individuals (of both sexes) with the TLR3^CC genotype had lower survival, and males - but not females - that carry the TLR3^C allele had significantly lower lifetime reproductive success than those with only the TLR3^A allele. These results indicate that positive selection on the TLR3^A allele, caused by an as yet unknown agent, is driving TLR3 evolution in the Seychelles warbler. No evidence of heterozygote advantage was detected. However, whether the positive selection observed is part of a longer-term pattern of balancing selection (through fluctuating selection or rare-allele advantage) cannot be resolved without tracking the TLR3^C allele over an extended time period.

Genetic drift and bottleneck do not influence diversity in Toll-like receptor genes at a small spatial scale in a Himalayan passerine

Genetic diversity is important for long-term viability of a population. Low genetic diversity reduces persistence and survival of populations and increases susceptibility to diseases. Comparisons of the neutral markers with functional loci such as immune genes [Toll-like receptors; TLR] can provide useful insights into evolutionary potential of a species and how the diversity of pathogens and selection pressures on their hosts are directly linked to their environment. In this study, we compare genetic diversity in neutral (eleven microsatellite loci) and adaptive (seven TLR loci) loci to determine genetic variation in a nonmigratory western Himalayan passerine, the black-throated tit (Aegithalos concinnus), distributed across an elevation gradient with varying degree of pathogen-mediated selection pressure. We further compare the diversity in TLR loci with a high-elevation sister species, the white-throated tit (Aegithalos niveogularis). Our results indicate a lack of population genetic structure in the black-throated tit and signatures of a past bottleneck. In contrast, we found high diversity in TLR loci and locus-specific (TLR7) signatures of pathogen-mediated selection, which was comparable to diversity in the white-throated tit. Levels of diversity at TLR5 locus corresponded very closely with neutral microsatellite variation. We found evidence of positive selection in TLR1LA, TLR5, and TLR7 loci highlighting the importance in pathogen recognition. Our finding demonstrates that reduction in neutral variation does not necessarily lead to reduction in functional genetic diversity and probably helps in revival of population in a widespread species.

Evidence of Pathogen-Induced Immunogenetic Selection across the Large Geographic Range of a Wild Seabird

Over evolutionary time, pathogen challenge shapes the immune phenotype of the host to better respond to an incipient threat. The extent and direction of this selection pressure depend on the local pathogen composition, which is in turn determined by biotic and abiotic features of the environment. However, little is known about adaptation to local pathogen threats in wild animals. The Gentoo penguin (Pygoscelis papua) is a species complex that lends itself to the study of immune adaptation because of its circumpolar distribution over a large latitudinal range, with little or no admixture between different clades. In this study, we examine the diversity in a key family of innate immune genes-the Toll-like receptors (TLRs)-across the range of the Gentoo penguin. The three TLRs that we investigated present varying levels of diversity, with TLR4 and TLR5 greatly exceeding the diversity of TLR7. We present evidence of positive selection in TLR4 and TLR5, which points to pathogen-driven adaptation to the local pathogen milieu. Finally, we demonstrate that two positively selected cosegregating sites in TLR5 are sufficient to alter the responsiveness of the receptor to its bacterial ligand, flagellin. Taken together, these results suggest that Gentoo penguins have experienced distinct pathogen-driven selection pressures in different environments, which may be important given the role of the Gentoo penguin as a sentinel species in some of the world's most rapidly changing environments.

Genetic diversity of toll-like receptor genes in the vulnerable Chinese egret (Egretta eulophotes)

Toll-like receptor (TLR) genes have recently been employed to assess genetic diversity, as they can be used to infer both demographic history and adaptation to environments with different pathogen pressure. Here, we sampled 120 individuals of the Chinese egret (Egretta eulophotes), a globally vulnerable species, from four breeding populations across China. We assessed the levels of genetic diversity, selection pressure, and population differentiation at seven TLR loci (TLR1LB, TLR2A, TLR3, TLR4, TLR5, TLR7, and TLR15). Using a variety of metrics (SNPs, heterozygosity, nucleotides, haplotypes), our analyses showed that genetic diversity was lower at 4 of the 7 TLR loci in the vulnerable Chinese egret compared to the more common little egret (Egretta garzetta). The selection test indicated TLRs, except for TLR5, were under purifying selection in TLR evolution, suggesting that low TLR genetic diversity in the Chinese egret may be caused by purifying selection. Moreover, analysis of molecular variance indicated low but significant population differentiation among four populations at all of the TLR loci in this egret. However, some comparisons based on fixation index analyses did not show significant population differentiation, and Bayesian clustering showed admixture. Our finding suggested that these four populations of the Chinese egret in China may be considered a single unit for conservation planning. These results, the new report of TLR genetic diversity in a long-distance migratory vulnerable Ardeid species, will provide fundamental TLR information for further studies on the conservation genetics of the Chinese egret and other Ardeids.

Adaptive genetic diversity and evidence of population genetic structure in the endangered Sierra Madre Sparrow (Xenospiza baileyi)

The magnitude and distribution of genetic diversity through space and time can provide useful information relating to evolutionary potential and conservation status in threatened species. In assessing genetic diversity in species that are of conservation concern, several studies have focused on the use of Toll-like receptors (TLRs). TLRs are innate immune genes related to pathogen resistance, and polymorphisms may reflect not only levels of functional diversity, but may also be used to assess genetic diversity within and among populations. Here, we combined four potentially adaptive markers (TLRs) with one mitochondrial (COI) marker to evaluate genetic variation in the endangered Sierra Madre Sparrow (Xenospiza baileyi). This species offers an ideal model to investigate population and evolutionary genetic processes that may be occurring in a habitat restricted endangered species with disjunct populations (Mexico City and Durango), the census sizes of which differ by an order of magnitude. TLRs diversity in the Sierra Madre Sparrow was relatively high, which was not expected given its two small, geographically isolated populations. Genetic diversity was different (but not significantly so) between the two populations, with less diversity seen in the smaller Durango population. Population genetic structure between populations was due to isolation and different selective forces acting on different TLRs; population structure was also evident in COI. Reduction of genetic diversity in COI was observed over 20 years in the Durango population, a result likely caused by habitat loss, a factor which may be the main cause of diversity decline generally. Our results provide information related to the ways in which adaptive variation can be altered by demographic changes due to human-mediated habitat alterations. Furthermore, our findings may help to guide conservation schemes for both populations and their restricted habitat.

Anthropogenic pressures negatively impact genomic diversity of the vulnerable seagrass Zostera capensis

Zostera capensis is a keystone species providing essential ecosystem services to southern African coastal systems. Like most seagrasses globally, Z. capensis is declining and under threat from anthropogenic pressures, and indicators of seagrass health and resilience may be of interest in preventing further declines. As intraspecific diversity is an important component of resilience, we used a pooled RADseq approach to generate genome-wide measures of variation across the entire South African distribution of Z. capensis. Using nucleotide diversity, heterozygosity and allelic richness we tested for associations with fine-scale anthropogenic pressure data compiled by the South African National Biodiversity Assessment using generalised linear models. Increased fishing effort, habitat loss, sand mining and a change in estuary flow dynamics were found to play an important role in decreasing nucleotide diversity and expected heterozygosity, most likely due to the loss of less resilient genotypes as a result of direct physical damage or indirect consequences. As the building block for adaptation, nucleotide diversity is particularly important for resilience. Because of this, as well as the fact that nucleotide diversity displayed the most distinct difference between the west and east coast, and responded most strongly to anthropogenic pressures, we suggest that this may be a useful measure for monitoring genetic or genomic variation. As genomic diversity influences resilience and resistance to disturbances, the remaining diversity in South African seagrass beds urgently needs to be conserved through restoration efforts and careful management of pressures.

Purifying Selection in the Toll-Like Receptors of Song Sparrows Melospiza melodia

Variation in immune gene sequences is known to influence resistance to infectious diseases and parasites, and hence survival and mate choice, across animal taxa. Toll-like receptors (TLRs) comprise one essential gene family in the vertebrate innate immune system and recognize evolutionarily conserved structures from all major microorganism classes. However, the causes and consequences of TLR variation in passerine birds remain largely unexplored. We examined 7 TLR genes in song sparrows (Melospiza melodia), a species that is studied across North America. We then examined sequences from 4 unduplicated TLRs (TLR1LB, TLR3, TLR4, and TLR15) from birds in 2 parts of the species' range (N = 27, N = 6), tested for evidence of selection, and conducted pilot analyses of the role of TLR heterozygosity in survival. We identified 45 SNPs: 19 caused changes in amino acid sequences and 2 of these were likely deleterious. We found no evidence of codon-level episodic positive selection but detected purifying selection at codons in all TLRs. Contrary to expectations we found no strong correlation between heterozygosity at TLRs and inbreeding coefficient f (estimate ± standard error [SE] = -0.68 ± 0.37, Radj2 = 0.08, F1,25 = 3.38, P = 0.08). In addition, pilot analyses revealed no relationship between TLR heterozygosity and survival (β ± SE: 0.09 ± 2.00, P = 0.96), possibly due to small sample size. Further analyses of genetic diversity in TLRs are likely to advance understanding of the effects of innate immune gene diversity on the fitness and persistence of wild populations.

Estimated six per cent loss of genetic variation in wild populations since the industrial revolution

Genetic variation is fundamental to population fitness and adaptation to environmental change. Human activities are driving declines in many wild populations and could have similar effects on genetic variation. Despite the importance of estimating such declines, no global estimate of the magnitude of ongoing genetic variation loss has been conducted across species. By combining studies that quantified recent changes in genetic variation across a mean of 27 generations for 91 species, we conservatively estimate a 5.4%-6.5% decline in within-population genetic diversity of wild organisms since the industrial revolution. This loss has been most severe for island species, which show a 27.6% average decline. We identified taxonomic and geographical gaps in temporal studies that must be urgently addressed. Our results are consistent with single time-point meta-analyses, which indicated that genetic variation is likely declining. However, our results represent the first confirmation of a global decline and provide an estimate of the magnitude of the genetic variation lost from wild populations.

Neutral variation does not predict immunogenetic variation in the European grayling (Thymallus thymallus)-implications for management

Preservation of genetic diversity is critical to successful conservation, and there is increasing demand for the inclusion of ecologically meaningful genetic information in management decisions. Supportive breeding programmes are increasingly implemented to combat declines in many species, yet their effect on adaptive genetic variation is understudied. This is despite the fact that supportive breeding may interfere with natural evolutionary processes. Here, we assessed the performance of neutral and adaptive markers (major histocompatibility complex; MHC) to inform management of European grayling (Thymallus thymallus), which routinely involves supplementation of natural populations with hatchery-reared fish (stocking). This study is the first to characterize MH II DAA and DAB loci in grayling and to investigate immune genetic variation in relation to management practice in this species. High-throughput Illumina sequencing of "introduced," "stocked native" and "non-stocked native" populations revealed significantly higher levels of allelic richness and heterozygosity for MH markers than microsatellites exclusively in non-stocked native populations. Likewise, significantly lower differentiation at the MH II than for microsatellites was apparent when considering non-stocked native populations, but not stocked populations. We developed a simulation model to test the effects of relaxation of selection during the early life stage within captivity. Dependent on the census population size and stocking intensity, there may be long-term effects of stocking on MH II, but not neutral genetic diversity. This is consistent with our empirical results. This study highlights the necessity for considering adaptive genetic variation in conservation decisions and raises concerns about the efficiency of stocking as a management practice.

A High-Quality, Long-Read De Novo Genome Assembly to Aid Conservation of Hawaii's Last Remaining Crow Species

Genome-level data can provide researchers with unprecedented precision to examine the causes and genetic consequences of population declines, which can inform conservation management. Here, we present a high-quality, long-read, de novo genome assembly for one of the world’s most endangered bird species, the ʻAlalā (Corvus hawaiiensis; Hawaiian crow). As the only remaining native crow species in Hawaiʻi, the ʻAlalā survived solely in a captive-breeding program from 2002 until 2016, at which point a long-term reintroduction program was initiated. The high-quality genome assembly was generated to lay the foundation for both comparative genomics studies and the development of population-level genomic tools that will aid conservation and recovery efforts. We illustrate how the quality of this assembly places it amongst the very best avian genomes assembled to date, comparable to intensively studied model systems. We describe the genome architecture in terms of repetitive elements and runs of homozygosity, and we show that compared with more outbred species, the ʻAlalā genome is substantially more homozygous. We also provide annotations for a subset of immunity genes that are likely to be important in conservation management, and we discuss how this genome is currently being used as a roadmap for downstream conservation applications.

Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

Adaptive genetic variation at three loci in South African vervet monkeys (Chlorocebus pygerythrus) and the role of selection within primates

Vervet monkeys (Chlorocebus pygerythrus) are one of the most widely distributed non-human primate species found in South Africa. They occur across all the South African provinces, inhabiting a large variety of habitats. These habitats vary sufficiently that it can be assumed that various factors such as pathogen diversity could influence populations in different ways. In turn, these factors could lead to varied levels of selection at specific fitness linked loci. The Toll-like receptor (TLR) gene family, which play an integral role in vertebrate innate immunity, is a group of fitness linked loci which has been the focus of much research. In this study, we assessed the level of genetic variation at partial sequences of two TLR loci (TLR4 and 7) and a reproductively linked gene, acrosin (ACR), across the different habitat types within the vervet monkey distribution range. Gene variation and selection estimates were also made among 11-21 primate species. Low levels of genetic variation for all three gene regions were observed within vervet monkeys, with only two polymorphic sites identified for TLR4, three sites for TLR7 and one site for ACR. TLR7 variation was positively correlated with high mean annual rainfall, which was linked to increased pathogen abundance. The observed genetic variation at TLR4 might have been influenced by numerous factors including pathogens and climatic conditions. The ACR exonic regions showed no variation in vervet monkeys, which could point to the occurrence of a selective sweep. The TLR4 and TLR7 results for the among primate analyses was mostly in line with previous studies, indicating a higher rate of evolution for TLR4. Within primates, ACR coding regions also showed signs of positive selection, which was congruent with previous reports on mammals. Important additional information to the already existing vervet monkey knowledge base was gained from this study, which can guide future research projects on this highly researched taxon as well as help conservation agencies with future management planning involving possible translocations of this species.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Decline in territory size and fecundity as a response to carrying capacity in an endangered songbird

Density-dependent processes are fundamental mechanisms for the regulation of populations. Ecological theories differ in their predictions on whether increasing population density leads to individual adjustments of survival and reproductive output or to dominance and monopolization of resources. Here, we use a natural experiment to examine which factors limit population growth in the only remaining population of the endangered pale-headed brush finch (Atlapetes pallidiceps). For three distinct phases (a phase of population suppression, 2001-2002; expansion due to conservation management, 2003-2008; and equilibrium phase, 2009-2014), we estimated demographic parameters with an integrated population model using population size, the proportion of successfully breeding pairs and their productivity, territory size, and mark-recapture data of adult birds. A low proportion of successful breeders due to brood parasitism (0.42, 95% credible interval 0.26-0.59) limited population growth before 2003; subsequent culling of the brood parasite resulted in a two-fold increase of the proportion of successful breeders during the 'expansion phase'. When the population approached the carrying capacity of its habitat, territory size declined by more than 50% and fecundity declined from 1.9 (1.54-2.27) to 1.3 (1.12-1.53) chicks per breeding pair, but the proportion of successful breeders remained constant (expansion phase: 0.85; 0.76-0.93; equilibrium phase: 0.86; 0.79-0.92). This study demonstrates that limiting resources can lead to individual adjustments instead of despotic behavior, and the individual reduction of reproductive output at high population densities is consistent with the slow life-history of many tropical species.

Lack of diversity at innate immunity Toll-like receptor genes in the Critically Endangered White-winged Flufftail (Sarothrura ayresi)

The White-winged Flufftail (Sarothrura ayresi) population is listed as globally Critically Endangered. White-winged Flufftails are only known to occur, with any regularity, in the high-altitude wetlands of South Africa and Ethiopia. Threats to the species include the limited number of suitable breeding sites in Ethiopia and severe habitat degradation and loss both in Ethiopia and South Africa. Toll-like receptors (TLRs) are increasingly being studied in a variety of taxa as a broader approach to determine functional genetic diversity. In this study, we confirm low genetic diversity in the innate immune regions of the White-winged Flufftail similar to that observed in other bird species that have undergone population bottlenecks. Low TLR diversity in White-winged Flufftail indicates that this species is more likely to be threatened by changes to the environment that would potentially expose the species to new diseases. Thus, conservation efforts should be directed towards maintaining pristine habitat for White-winged Flufftail in its current distribution range. To date, no studies on immunogenetic variation in White-winged Flufftail have been conducted and to our knowledge, this is the first study of TLR genetic diversity in a critically endangered species.

Drift, not selection, shapes toll-like receptor variation among oceanic island populations

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll‐like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen‐mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.

Genetic diversity does not explain variation in extra-pair paternity in multiple populations of a songbird

Many songbirds are socially monogamous but genetically polyandrous, mating with individuals outside their pair bonds. Extra-pair paternity (EPP) varies within and across species, but reasons for this variation remain unclear. One possible source of variation is population genetic diversity, which has been shown in interspecific meta-analyses to correlate with EPP but which has limited support from intraspecific tests. Using eight populations of the genetically polyandrous red-winged blackbird (Agelaius phoeniceus), including an island population, we investigated whether population-level differences in genetic diversity led to differences in EPP. We first measured genetic diversity over 10 microsatellite loci and found, as predicted, low genetic diversity in the island population. Additional structure analyses with multilocus genotypes and mtDNA showed the island population to be distinct from the continental populations. However, the island population's EPP rate fell in the middle of the continental populations' distribution, whereas the continental populations themselves showed significant variation in EPP. This result suggests that genetic diversity by itself is not a predictor of EPP rate. We discuss reasons for the departure from previous results, including hypotheses for EPP that do not solely implicate female-driven behaviour.