Phylogenomics of the extinct Heath Hen provides support for sex-biased introgression among extant prairie grouse.

Rapid divergence and subsequent reoccurring patterns of gene flow can complicate our ability to discern phylogenetic relationships among closely related species. To what degree such patterns may differ across the genome can provide an opportunity to extrapolate better how life history constraints may influence species boundaries. By exploring differences between autosomal and Z (or X) chromosomal-derived phylogenetic patterns, we can better identify factors that may limit introgression despite patterns of incomplete lineage sorting among closely related taxa. Here, using a whole-genome resequencing approach coupled with an exhaustive sampling of subspecies within the recently divergent prairie grouse complex (genus: Tympanuchus), including the extinct Heath Hen (T. cupido cupido), we show that their phylogenomic history differs depending on autosomal or Z-chromosome partitioned SNPs. Because the Heath Hen was allopatric relative to the other prairie grouse taxa, its phylogenetic signature should not be influenced by gene flow. In contrast, all the other extant prairie grouse taxa, except Attwater's Prairie-chicken (T. c. attwateri), possess overlapping contemporary geographic distributions and have been known to hybridize. After excluding samples that were likely translocated prairie grouse from the Midwest to the eastern coastal states or their resulting hybrids with mainland Heath Hens, species tree analyses based on autosomal SNPs consistently identified a paraphyletic relationship with regard to the Heath Hen with Lesser Prairie-chicken (T. pallidicinctus) sister to Greater Prairie-chicken (T. c. pinnatus) regardless of genic or intergenic partitions. In contrast, species trees based on the Z-chromosome were consistent with Heath Hen sister to a clade that included its conspecifics, Greater and Attwater's Prairie-chickens (T. c. attwateri). These results were further explained by historic gene flow, as shown with an excess of autosomal SNPs shared between Lesser and Greater Prairie-chickens but not with the Z-chromosome. Phylogenetic placement of Sharp-tailed Grouse (T. phasianellus), however, did not differ among analyses and was sister to a clade that included all other prairie grouse despite low levels of autosomal gene flow with Greater Prairie-chicken. These results, along with strong sexual selection (i.e., male hybrid behavioral isolation) and a lek breeding system (i.e., high variance in male mating success), are consistent with a pattern of female-biased introgression between prairie grouse taxa with overlapping geographic distributions. Additional study is warranted to explore how genomic components associated with the Z-chromosome influence the phenotype and thereby impact species limits among prairie grouse taxa despite ongoing contemporary gene flow.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies.

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

A supergene underlies linked variation in color and morphology in a Holarctic songbird.

The genetic architecture of a phenotype can have considerable effects on the evolution of a trait or species. Characterizing genetic architecture provides insight into the complexity of a given phenotype and, potentially, the role of the phenotype in evolutionary processes like speciation. We use genome sequences to investigate the genetic basis of phenotypic variation in redpoll finches (Acanthis spp.). We demonstrate that variation in redpoll phenotype is broadly controlled by a ~55-Mb chromosomal inversion. Within this inversion, we find multiple candidate genes related to melanogenesis, carotenoid coloration, and bill shape, suggesting the inversion acts as a supergene controlling multiple linked traits. A latitudinal gradient in ecotype distribution suggests supergene driven variation in color and bill morphology are likely under environmental selection, maintaining supergene haplotypes as a balanced polymorphism. Our results provide a mechanism for the maintenance of ecotype variation in redpolls despite a genome largely homogenized by gene flow.

Migratory movements of Atlantic puffins Fratercula arctica naumanni from high Arctic Greenland.

Although the Atlantic puffin Fratercula arctica is well studied throughout its temperate and low Arctic breeding range, few have studied the species in its far northern distribution. This study is the first to present data on the migratory movements of the "large-billed" subspecies, F. a. naumanni, that breeds in the high Arctic and which has significantly larger body size than those farther south. During 2013-2015, migration tracks were collected from nine adult puffins (6 males and 3 females) tagged with geolocators in northwest Greenland. Overall, female puffins traveled farther than males on their annual migration, with one female puffin traveling over 13,600 km, which was nearly a third farther than any tagged male in our study. Differential migration was observed in migratory phenology and route, with males using a form of chain migration with acute synchrony between individuals while females appeared to largely use leap-frog migration and showed little synchrony between individuals. Extreme sexual segregation in wintering areas was evidenced by two females that migrated to the southern limit of the species' range while the six males remained at the northern limit, and wintered along the sea ice edge during portions of the non-breeding season. Male puffins thus wintered in regions with sea surface temperatures up to 10° C cooler than female puffins, and in areas with generally colder sea surface temperatures when compared to previously known wintering areas of temperate and low Arctic puffin breeding populations. The degree to which body size enables male F. a. naumanni to remain in colder waters likely reflects differing life history constraints between sexes and populations (i.e., subspecies). Further study is warranted to investigate how recent changes in climate have further exacerbated the observed differences between sexes in high Arctic puffins and possibly other marine avian species.

Bird feathers are potential biomonitors for airborne elemental carbon.

Birds can serve as effective biomonitors of air pollution, yet few studies have quantified external particulate matter accumulation on bird feathers. Biomonitoring of airborne elemental carbon (EC) is of critical significance because EC is a component of particulate matter with adverse effects on air quality and human health. To assess their effectiveness for use in EC monitoring, we compared EC accumulation on bird feathers at two sites that differed in vehicular traffic volume in an urban environment within the Dallas-Fort Worth Metropolitan Area, USA. Moulted flight feathers from domestic chickens were experimentally exposed to ambient EC pollution for 5 days in two urban microenvironments 1.5 km distant from each other that differed in traffic volume--adjacent to an interstate highway and a university campus bus stop. Feathers near the highway accumulated approximately eight times more EC (307 ± 34 µg m-2 day-1), on average, than feathers near the bus stop (40 ± 9 µg m-2 day-1). These findings indicate that EC accumulation on feathers varies over short distances within urban areas and that bird feathers potentially can be used for biomonitoring airborne EC.

Correction to: Status and trends of circumpolar peregrine falcon and gyrfalcon populations.

While collating contributions and comments from 36 researchers, the coordinating authors accidentally omitted Dr. Suzanne Carrière from the list of contributing co-authors. Dr. Carrière's data are described in Tables 1 and 3, Figure 2 and several places in the narrative.The new author list is thus updated in this article.

Status and trends of circumpolar peregrine falcon and gyrfalcon populations.

The peregrine falcon (Falco peregrinus) and the gyrfalcon (Falco rusticolus) are top avian predators of Arctic ecosystems. Although existing monitoring efforts are well established for both species, collaboration of activities among Arctic scientists actively involved in research of large falcons in the Nearctic and Palearctic has been poorly coordinated. Here we provide the first overview of Arctic falcon monitoring sites, present trends for long-term occupancy and productivity, and summarize information describing abundance, distribution, phenology, and health of the two species. We summarize data for 24 falcon monitoring sites across the Arctic, and identify gaps in coverage for eastern Russia, the Arctic Archipelago of Canada, and East Greenland. Our results indicate that peregrine falcon and gyrfalcon populations are generally stable, and assuming that these patterns hold beyond the temporal and spatial extents of the monitoring sites, it is reasonable to suggest that breeding populations at broader scales are similarly stable. We have highlighted several challenges that preclude direct comparisons of Focal Ecosystem Components (FEC) attributes among monitoring sites, and we acknowledge that methodological problems cannot be corrected retrospectively, but could be accounted for in future monitoring. Despite these drawbacks, ample opportunity exists to establish a coordinated monitoring program for Arctic-nesting raptor species that supports CBMP goals.

Evaluation of a Chicken 600K SNP genotyping array in non-model species of grouse.

The use of single nucleotide polymorphism (SNP) arrays to generate large SNP datasets for comparison purposes have recently become an attractive alternative to other genotyping methods. Although most SNP arrays were originally developed for domestic organisms, they can be effectively applied to wild relatives to obtain large panels of SNPs. In this study, we tested the cross-species application of the Affymetrix 600K Chicken SNP array in five species of North American prairie grouse (Centrocercus and Tympanuchus genera). Two individuals were genotyped per species for a total of ten samples. A high proportion (91%) of the total 580 961 SNPs were genotyped in at least one individual (73-76% SNPs genotyped per species). Principal component analysis with autosomal SNPs separated the two genera, but failed to clearly distinguish species within genera. Gene ontology analysis identified a set of genes related to morphogenesis and development (including genes involved in feather development), which may be primarily responsible for large phenotypic differences between Centrocercus and Tympanuchus grouse. Our study provided evidence for successful cross-species application of the chicken SNP array in grouse which diverged ca. 37 mya from the chicken lineage. As far as we are aware, this is the first reported application of a SNP array in non-passerine birds, and it demonstrates the feasibility of using commercial SNP arrays in research on non-model bird species.

Major histocompatibility complex variation and blood parasites in resident and migratory populations of the common yellowthroat.

Genes of the major histocompatibility complex (MHC) are a critical part of the adaptive immune response, and the most polymorphic genes in the vertebrate genome, especially in passerine birds. This diversity is thought to be influenced by exposure to pathogens which can vary in relation to numerous factors. Migratory behaviour may be a particularly important trait to consider because migratory birds are exposed to a greater number of different pathogens and parasites at both breeding (i.e. temperate) and overwintering (i.e. tropical and subtropical) areas, as well as at stopover sites during migration. Thus, migrants are predicted to have greater MHC diversity than residents. We compared MHC variation, at both class I and II, and levels of haemosporidian infection between one resident and two migratory populations of the common yellowthroat (Geothlypis trichas). We found that residents were less likely to be infected with haemosporidian parasites and had lower MHC diversity at class I; however, variation at MHC class II was greater in residents than migrants, contrary to our prediction. These patterns were not likely to be caused by differences in population demography as genomewide heterozygosity (based on 9225 single nucleotide polymorphisms) was high in all three populations and not correlated with MHC variation. Our different results for MHC class I and II suggest that studies of immune gene variation in relation to life history need to consider that there could be different selection pressures arising from intracellular (class I) and extracellular (class II) pathogens in different populations.

Extensive shared polymorphism at non-MHC immune genes in recently diverged North American prairie grouse.

Gene polymorphisms shared between recently diverged species are thought to be widespread and most commonly reflect introgression from hybridization or retention of ancestral polymorphism through incomplete lineage sorting. Shared genetic diversity resulting from incomplete lineage sorting is usually maintained for a relatively short period of time, but under strong balancing selection it may persist for millions of years beyond species divergence (balanced trans-species polymorphism), as in the case of the major histocompatibility complex (MHC) genes. However, balancing selection is much less likely to act on non-MHC immune genes. The aim of this study was to investigate the patterns of shared polymorphism and selection at non-MHC immune genes in five grouse species from Centrocercus and Tympanuchus genera. For this purpose, we genotyped five non-MHC immune genes that do not interact directly with pathogens, but are involved in signaling and regulate immune cell growth. In contrast to previous studies with MHC, we found no evidence for balancing selection or balanced trans-species polymorphism among the non-MHC immune genes. No haplotypes were shared between genera and in most cases more similar allelic variants sorted by genus. Between species within genera, however, we found extensive shared polymorphism, which was most likely attributable to introgression or incomplete lineage sorting following recent divergence and large ancestral effective population size (i.e., weak genetic drift). Our study suggests that North American prairie grouse may have attained relatively low degree of reciprocal monophyly at nuclear loci and reinforces the rarity of balancing selection in non-MHC immune genes.

Spatial Microbial Composition Along the Gastrointestinal Tract of Captive Attwater's Prairie Chicken.

Gastrointestinal microbiota is increasingly recognized as an important component of individual health, and therefore, our ability to quantify its diversity accurately is central for exploring different ways to improve health. Non-invasive sampling methods, such as cloaca swabs, are often used to measure gastrointestinal microbiota diversity within an individual. However, few studies have addressed to what degree differences exist in microbial community composition along the gastrointestinal tract, and measures obtained from the cloaca may not actually represent the diversity present elsewhere in the gastrointestinal tract. In this study, we systematically characterized the gastrointestinal microbial community of the critically endangered Attwater's Prairie chicken (Tympanuchus cupido attwateri) by opportunistically sampling four different locations (ileum, cecum, large intestine, and cloaca) along the gastrointestinal tract of eight individuals. Spatial variation of microbial community was observed at different sampling locations within the gastrointestinal tract. The cecum harbored the most diverse and significantly different microbiota from the other locations, while the microbial α- and ß-diversities were similar in the ileum, large intestine, and cloaca. The results of this study provide evidence that microbiota diversity can differ depending on sampling location and metric used to quantify diversity. As shown here, non-invasive cloacal sampling strategies may reflect microbiota diversity elsewhere in the gastrointestinal tract, yet caution is warranted when making generalizations in terms of the microbiota diversity correlations when samples are obtained from a single location within the gastrointestinal tract.

Multi-locus phylogenetic inference among New World Vultures (Aves: Cathartidae).

New World Vultures are large-bodied carrion feeding birds in the family Cathartidae, currently consisting of seven species from five genera with geographic distributions in North and South America. No study to date has included all cathartid species in a single phylogenetic analysis. In this study, we investigated the phylogenetic relationships among all cathartid species using five nuclear (nuc; 4060bp) and two mitochondrial (mt; 2165bp) DNA loci with fossil calibrated gene tree (27 outgroup taxa) and coalescent-based species tree (2 outgroup taxa) analyses. We also included an additional four nuclear loci (2578bp) for the species tree analysis to explore changes in nodal support values. Although the stem lineage is inferred to have originated ∼69 million years ago (Ma; 74.5-64.9 credible interval), a more recent basal split within Cathartidae was recovered at ∼14Ma (17.1-11.1 credible interval). Two primary clades were identified: (1) Black Vulture (Coragyps atratus) together with the three Cathartes species (Lesser C. burrovianus and Greater C. melambrotus Yellow-headed Vultures, and Turkey Vulture C. aura), and (2) King Vulture (Sarcoramphus papa), California (Gymnogyps californianus) and Andean (Vultur gryphus) Condors. Support for taxon relationships within the two basal clades were inconsistent between analyses with the exception of Black Vulture sister to a monophyletic Cathartes clade. Increased support for a yellow-headed vulture clade was recovered in the species tree analysis using the four additional nuclear loci. Overall, these results are in agreement with cathartid life history (e.g. olfaction ability and behavior) and contrasting habitat affinities among sister taxa with overlapping geographic distributions. More research is needed using additional molecular loci to further resolve the phylogenetic relationships within the two basal cathartid clades, as speciation appeared to have occurred in a relatively short period of time.

Specific alleles at immune genes, rather than genome-wide heterozygosity, are related to immunity and survival in the critically endangered Attwater's prairie-chicken.

The negative effects of inbreeding on fitness are serious concerns for populations of endangered species. Reduced fitness has been associated with lower genome-wide heterozygosity and immune gene diversity in the wild; however, it is rare that both types of genetic measures are included in the same study. Thus, it is often unclear whether the variation in fitness is due to the general effects of inbreeding, immunity-related genes or both. Here, we tested whether genome-wide heterozygosity (20 990 SNPs) and diversity at nine immune genes were better predictors of two measures of fitness (immune response and survival) in the endangered Attwater's prairie-chicken (Tympanuchus cupido attwateri). We found that postrelease survival of captive-bred birds was related to alleles of the innate (Toll-like receptors, TLRs) and adaptive (major histocompatibility complex, MHC) immune systems, but not to genome-wide heterozygosity. Likewise, we found that the immune response at the time of release was related to TLR and MHC alleles, and not to genome-wide heterozygosity. Overall, this study demonstrates that immune genes may serve as important genetic markers when monitoring fitness in inbred populations and that in some populations specific functional genes may be better predictors of fitness than genome-wide heterozygosity.

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens.

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

Mercury Concentrations in Birds from Two Atmospherically Contaminated Sites in North Texas, USA.

Mercury (Hg) is a ubiquitous and highly toxic contaminant that can have negative effects on wildlife. Only a few studies have measured Hg concentrations in birds from the south central United States, and the potential threat of Hg contamination to birds in this region is largely unknown. In the present study, we assess Hg concentrations in blood and feathers from five bird species [eastern bluebird (Sialis sialis), Carolina wren (Thryothorus ludovicianus), wood duck (Aix sponsa), great egret (Ardea alba), and great blue heron (Ardea herodias)] that occupy different trophic levels at Caddo Lake and Lewisville Lake, located in northeast and north central Texas, respectively. Both sites are contaminated with Hg from the atmosphere. Adult passerines had higher Hg concentrations in their blood than conspecific nestlings. Mercury concentrations in feathers differed between species by more than an order of magnitude with large piscivorous species having higher concentrations than smaller insectivorous species. Mercury concentrations in eastern bluebirds were higher at Caddo Lake than Lewisville Lake. The present study represents one of the first studies of Hg concentrations in multiple bird species in north Texas and suggests that Hg concentrations in birds from atmospherically polluted sites in this region may be high enough to compromise fitness in those species.

Differential introgression and effective size of marker type influence phylogenetic inference of a recently divergent avian group (Phasianidae: Tympanuchus).

Life history strategies can influence the effective population size (Ne) of loci differently based on their mode of inheritance. Recognizing how this may affect the rate of lineage sorting among marker types is important for studies focused on resolving phylogenetic relationships among recently divergent taxa. In this study, we use gene tree, coalescent-based species tree, and isolation-with-migration analyses to explore the differences between marker types (autosomal, Z-linked, and mitochondrial) in resolving phylogenetic relationships among North American prairie grouse (Tympanuchus). We found that Z-linked loci were more likely to identify monophyletic relationships among prairie grouse species compared to autosomal and mtDNA loci in both species and gene tree analyses, with species tree analyses outperforming gene trees. These results were further supported with isolation-with-migration analyses, where Z-linked loci largely followed a strict isolation model while autosomal loci were more likely to fit a model with gene flow between species following population divergence. While accounting for differences in inheritance pattern (or Ne) for marker type, results suggest that additional factors, such as strong sexual selection and sex-biased introgression (i.e., male-biased postzygotic hybrid behavioral isolation or "unsexy son"), may further explain the decreased diversity levels and increased rate of lineage sorting observed with the Z-linked loci relative to autosomal and mtDNA loci. In fact, to our knowledge no hybrid male prairie grouse have been observed breeding in the wild, yet hybrid females along with backcross females are known to produce viable offspring. Overall, this study highlights that more work is needed to determine how complex models of gene flow (i.e., sex biased introgression) and differences in the effective size among marker types based on differing life history strategies influence divergence date estimation and species delimitation.

Rapid diversification of falcons (Aves: Falconidae) due to expansion of open habitats in the Late Miocene.

Understanding how and why lineages diversify is central to understanding the origins of biological diversity. The avian family Falconidae (caracaras, forest-falcons, falcons) has an uneven distribution of species among multiple well-supported clades, and provides a useful system for testing hypotheses about diversification rate and correlation with environmental changes. We analyzed eight independent loci for 1-7 individuals from each of the 64 currently recognized Falconidae species, together with two fossil falconid temporal calibrations, to assess phylogeny, absolute divergence times and potential shifts in diversification rate. Our analyses supported similar diversification ages in the Early to Middle Miocene for the three traditional subfamilies, Herpetotherinae, Polyborinae and Falconinae. We estimated that divergences within the subfamily Falconinae began about 16mya and divergences within the most species-rich genus, Falco, including about 60% of all Falconidae species, began about 7.5mya. We found evidence for a significant increase in diversification rate at the basal phylogenetic node for the genus Falco, and the timing for this rate shift correlates generally with expansion of C4 grasslands beginning around the Miocene/Pliocene transition. Concomitantly, Falco lineages that are distributed primarily in grassland or savannah habitats, as opposed to woodlands, and exhibit migratory, as opposed to sedentary, behavior experienced a higher diversification rate.

A comparison of pedigree- and DNA-based measures for identifying inbreeding depression in the critically endangered Attwater's Prairie-chicken.

The primary goal of captive breeding programmes for endangered species is to prevent extinction, a component of which includes the preservation of genetic diversity and avoidance of inbreeding. This is typically accomplished by minimizing mean kinship in the population, thereby maintaining equal representation of the genetic founders used to initiate the captive population. If errors in the pedigree do exist, such an approach becomes less effective for minimizing inbreeding depression. In this study, both pedigree- and DNA-based methods were used to assess whether inbreeding depression existed in the captive population of the critically endangered Attwater's Prairie-chicken (Tympanuchus cupido attwateri), a subspecies of prairie grouse that has experienced a significant decline in abundance and concurrent reduction in neutral genetic diversity. When examining the captive population for signs of inbreeding, variation in pedigree-based inbreeding coefficients (f(pedigree)) was less than that obtained from DNA-based methods (f(DNA)). Mortality of chicks and adults in captivity were also positively correlated with parental relatedness (r(DNA)) and f(DNA), respectively, while no correlation was observed with pedigree-based measures when controlling for additional variables such as age, breeding facility, gender and captive/release status. Further, individual homozygosity by loci (HL) and parental rDNA values were positively correlated with adult mortality in captivity and the occurrence of a lethal congenital defect in chicks, respectively, suggesting that inbreeding may be a contributing factor increasing the frequency of this condition among Attwater's Prairie-chickens. This study highlights the importance of using DNA-based methods to better inform management decisions when pedigrees are incomplete or errors may exist due to uncertainty in pairings.

Genetics of plumage color in the Gyrfalcon (Falco rusticolus): analysis of the melanocortin-1 receptor gene.

Genetic variation at the melanocortin-1 receptor (MC1R) gene is correlated with melanin color variation in a few reported vertebrates. In Gyrfalcon (Falco rusticolus), plumage color variation exists throughout their arctic and subarctic circumpolar distribution, from white to gray and almost black. Multiple color variants do exist within the majority of populations; however, a few areas (e.g., northern Greenland and Iceland) possess a single color variant. Here, we show that the white/melanic color pattern observed in Gyrfalcons is explained by allelic variation at MC1R. Six nucleotide substitutions in MC1R resulted in 9 alleles that differed in geographic frequency with at least 2 MC1R alleles observed in almost all sampled populations in Greenland, Iceland, Canada, and Alaska. In north Greenland, where white Gyrfalcons predominate, a single MC1R allele was observed at high frequency (>98%), whereas in Iceland, where only gray Gyrfalcons are known to breed, 7 alleles were observed. Of the 6 nucleotide substitutions, 3 resulted in amino acid substitutions, one of which (Val(128)Ile) was perfectly associated with the white/melanic polymorphism. Furthermore, the degree of melanism was correlated with number of MC1R variant alleles, with silver Gyrfalcons all heterozygous and the majority of dark gray individuals homozygous (Ile(128)). These results provide strong support that MC1R is associated with plumage color in this species.

Drift and selection influence geographic variation at immune loci of prairie-chickens.

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.