How genomics can help biodiversity conservation.

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians.

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Understanding the evolution of immune genes in jawed vertebrates.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

King or royal family? Testing for species boundaries in the King Cobra, Ophiophagus hannah (Cantor, 1836), using morphology and multilocus DNA analyses.

In widespread species, the diverse ecological conditions in which the populations occur, and the presence of many potential geographical barriers through their range are expected to have created ample opportunities for the evolution of distinct, often cryptic lineages. In this work, we tested for species boundaries in one such widespread species, the king cobra, Ophiophagus hannah (Cantor, 1836), a largely tropical elapid snake distributed across the Oriental realm. Based on extensive geographical sampling across most of the range of the species, we initially tested for candidate species (CS) using Maximum-Likelihood analysis of mitochondrial genes. We then tested the resulting CS using both morphological data and sequences of three single-copy nuclear genes. We used snapclust to determine the optimal number of clusters in the nuclear dataset, and Bayesian Phylogenetics and Phylogeography (BPP) to test for likely species status. We used non-metric multidimensional scaling (nMDS) analysis for discerning morphological separation. We recovered four independently evolving, geographically separated lineages that we consider Confirmed Candidate Species: (1) Western Ghats lineage; (2) Indo-Chinese lineage (3) Indo-Malayan lineage; (4) Luzon Island lineage, in the Philippine Archipelago. We discuss patterns of lineage divergence, particularly in the context of low morphological divergence, and the conservation implications of recognizing several endemic king cobra lineages.

Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish.

The use of genetic markers in the context of conservation is largely being outcompeted by whole-genome data. Comparative studies between the two are sparse, and the knowledge about potential effects of this methodology shift is limited. Here, we used whole-genome sequencing data to assess the genetic status of peripheral populations of the wels catfish (Silurus glanis), and discuss the results in light of a recent microsatellite study of the same populations. The Swedish populations of the wels catfish have suffered from severe declines during the last centuries and persists in only a few isolated water systems. Fragmented populations generally are at greater risk of extinction, for example due to loss of genetic diversity, and may thus require conservation actions. We sequenced individuals from the three remaining native populations (Båven, Emån, and Möckeln) and one reintroduced population of admixed origin (Helge å), and found that genetic diversity was highest in Emån but low overall, with strong differentiation among the populations. No signature of recent inbreeding was found, but a considerable number of short runs of homozygosity were present in all populations, likely linked to historically small population sizes and bottleneck events. Genetic substructure within any of the native populations was at best weak. Individuals from the admixed population Helge å shared most genetic ancestry with the Båven population (72%). Our results are largely in agreement with the microsatellite study, and stresses the need to protect these isolated populations at the northern edge of the distribution of the species.

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation.

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene-environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

Small-scale population divergence is driven by local larval environment in a temperate amphibian.

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi.

BACKGROUND: The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai-Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history. RESULTS: Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche. CONCLUSIONS: The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans.

BACKGROUND: While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns. RESULTS: We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values. CONCLUSION: Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Occurrence of Batrachochytrium dendrobatidis in Sweden: higher infection prevalence in southern species.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused worldwide declines in amphibian populations. While Bd is widespread in southern and central Europe, its occurrence and distribution in northernmost Europe is mostly unknown. We surveyed for Bd in breeding anurans in Sweden by sampling 1917 amphibians from 101 localities and 3 regions in Sweden (southern, northern and central). We found that Bd was widespread in southern and central Sweden, occurring in all 9 investigated species and in 45.5% of the 101 localities with an overall prevalence of 13.8%. No infected individuals were found in the 4 northern sites sampled. The records from central Sweden represent the northernmost records of Bd in Europe. While the proportion of sites positive for Bd was similar between the southern and central regions, prevalence was much higher in the southern region. This was because southern species with a distribution mainly restricted to southernmost Sweden had a higher prevalence than widespread generalist species. The nationally red-listed green toad Bufotes variabilis and the fire-bellied toad Bombina bombina had the highest prevalence (61.4 and 48.9%, respectively). Across species, Bd prevalence was strongly positively, correlated with water temperature at the start of egg laying. However, no individuals showing visual signs of chytridiomycosis were found in the field. These results indicate that Bd is widespread and common in southern and central Sweden with southern species, breeding in higher temperatures and with longer breeding periods, having higher prevalence. However, the impact of Bd on amphibian populations in northernmost Europe remains unknown.

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation.

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long-term persistence of populations, understanding how these processes affect variation over large-scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700-km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction-site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post-glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate-mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.

Genetic basis of amphibian larval development along a latitudinal gradient: Gene diversity, selection and links with phenotypic variation in transcription factor C/EBP-1.

Ectotherm development rates often show adaptive divergence along climatic gradients, but the genetic basis for this variation is rarely studied. Here, we investigated the genetic basis for phenotypic variation in larval development in the moor frog Rana arvalis from five regions along a latitudinal gradient from Germany to northern Sweden. We focused on the C/EBP-1 gene, a transcription factor associated with larval development time. Allele frequencies at C/EBP-1 varied strongly among geographical regions. Overall, the distribution of alleles along the gradient was in concordance with the dual post-glacial colonization routes into Scandinavia, with a large number of alleles exclusively present along the southern colonization route. Only three of 38 alleles were shared between the routes. Analysis of contemporary selection on C/EBP-1 showed divergent selection among the regions, probably reflecting adaptation to the local environmental conditions, although this was especially strong between southern and northern regions coinciding also with lineages from different colonization routes. Overall, the C/EBP-1 gene has historically been under purifying selection, but two specific amino acid positions showed significant signals of positive selection. These positions showed divergence between southern and northern regions, and we suggest that they are functionally involved in the climatic adaptation of larval development. Using phenotypic data from a common garden experiment, we found evidence for specific C/EBP-1 alleles being correlated with larval development time, suggesting a functional role in adaptation of larval development to large-scale climatic variation.

Post-glacial colonization routes coincide with a life-history breakpoint along a latitudinal gradient.

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post-glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical QST -FST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The QST estimates between the colonization routes were four times higher than FST , indicating a prominent role for natural selection. QST within the colonization routes did not generally differ from FST , but we found temperature-dependent adaptive divergence close to the contact zone. These results indicate that lineage-specific variation can account for much of the adaptive divergence along a latitudinal gradient.

Urban environment and reservoir host species are associated with Batrachochytrium dendrobatidis infection prevalence in the common toad.

Human-induced changes of the environment, including landscape alteration and habitat loss, may affect wildlife disease dynamics and have important ramifications for wildlife conservation. Amphibians are among the vertebrate taxa most threatened by anthropogenic habitat change. The emerging fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused extinctions and population declines in hundreds of anuran species globally. We studied how the urban landscape is associated with the prevalence of Bd infections by sampling 655 anurans of 3 species (mainly the common toad Bufo bufo) in 42 ponds surrounded by different amounts of urban habitat (defined as towns, cities or villages). We also examined the association between Bd infections and a potential reservoir host species (the moor frog Rana arvalis). We found that 38% of the sites were positive for Bd with an infection prevalence of 4.4%. The extent of urban landscape was negatively correlated with Bd infection prevalence. However, the positive association of Bd with the presence of the possible reservoir species was substantially stronger than the urban effects. The body condition index of B. bufo was negatively associated with Bd infection. This Bd effect was stronger than the negative effect of urban landscape on body condition. Our results suggest that urban environments in Sweden have a negative impact on Bd infections, while the presence of the reservoir species has a positive impact on Bd prevalence. Our study also highlights the potential importance of Bd infection on host fitness, especially in rural landscapes.

Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian.

BACKGROUND: Past events like fluctuations in population size and post-glacial colonization processes may influence the relative importance of genetic drift, migration and selection when determining the present day patterns of genetic variation. We disentangle how drift, selection and migration shape neutral and adaptive genetic variation in 12 moor frog populations along a 1700 km latitudinal gradient. We studied genetic differentiation and variation at a MHC exon II locus and a set of 18 microsatellites. RESULTS: Using outlier analyses, we identified the MHC II exon 2 (corresponding to the ß-2 domain) locus and one microsatellite locus (RCO8640) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. STRUCTURE and DAPC analyses on the neutral microsatellites assigned populations to a northern and a southern cluster, reflecting two different post-glacial colonization routes found in previous studies. Genetic variation overall was lower in the northern cluster. The signature of selection on MHC exon II was weaker in the northern cluster, possibly as a consequence of smaller and more fragmented populations. CONCLUSION: Our results show that historical demographic processes combined with selection and drift have led to a complex pattern of differentiation along the gradient where some loci are more divergent among populations than predicted from drift expectations due to diversifying selection, while other loci are more uniform among populations due to stabilizing selection. Importantly, both overall and MHC genetic variation are lower at northern latitudes. Due to lower evolutionary potential, the low genetic variation in northern populations may increase the risk of extinction when confronted with emerging pathogens and climate change.

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media).

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds.

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum.

Looking into the past - the reaction of three grouse species to climate change over the last million years using whole genome sequences.

Tracking past population fluctuations can give insight into current levels of genetic variation present within species. Analysing population dynamics over larger timescales can be aligned to known climatic changes to determine the response of species to varying environments. Here, we applied the Pairwise Sequentially Markovian Coalescent (psmc) model to infer past population dynamics of three widespread grouse species; black grouse, willow grouse and rock ptarmigan. This allowed the tracking of the effective population size (Ne ) of all three species beyond 1 Mya, revealing that (i) early Pleistocene cooling (~2.5 Mya) caused an increase in the willow grouse and rock ptarmigan populations, (ii) the mid-Brunhes event (~430 kya) and following climatic oscillations decreased the Ne of willow grouse and rock ptarmigan, but increased the Ne of black grouse and (iii) all three species reacted differently to the last glacial maximum (LGM) - black grouse increased prior to it, rock ptarmigan experienced a severe bottleneck and willow grouse was maintained at large population size. We postulate that the varying psmc signal throughout the LGM depicts only the local history of the species. Nevertheless, the large population fluctuations in willow grouse and rock ptarmigan indicate that both species are opportunistic breeders while black grouse tracks the climatic changes more slowly and is maintained at lower Ne . Our results highlight the usefulness of the psmc approach in investigating species' reaction to climate change in the deep past, but also that caution should be taken in drawing general conclusions about the recent past.

Whole genome sequencing of the black grouse (Tetrao tetrix): reference guided assembly suggests faster-Z and MHC evolution.

BACKGROUND: The different regions of a genome do not evolve at the same rate. For example, comparative genomic studies have suggested that the sex chromosomes and the regions harbouring the immune defence genes in the Major Histocompatability Complex (MHC) may evolve faster than other genomic regions. The advent of the next generation sequencing technologies has made it possible to study which genomic regions are evolutionary liable to change and which are static, as well as enabling an increasing number of genome studies of non-model species. However, de novo sequencing of the whole genome of an organism remains non-trivial. In this study, we present the draft genome of the black grouse, which was developed using a reference-guided assembly strategy. RESULTS: We generated 133 Gbp of sequence data from one black grouse individual by the SOLiD platform and used a combination of de novo assembly and chicken reference genome mapping to assemble the reads into 4572 scaffolds with a total length of 1022 Mb. The draft genome well covers the main chicken chromosomes 1 ~ 28 and Z which have a total length of 1001 Mb. The draft genome is fragmented, but has a good coverage of the homologous chicken genes. Especially, 33.0% of the coding regions of the homologous genes have more than 90% proportion of their sequences covered. In addition, we identified ~1 M SNPs from the genome and identified 106 genomic regions which had a high nucleotide divergence between black grouse and chicken or between black grouse and turkey. CONCLUSIONS: Our results support the hypothesis that the chromosome X (Z) evolves faster than the autosomes and our data are consistent with the MHC regions being more liable to change than the genome average. Our study demonstrates how a moderate sequencing effort can be combined with existing genome references to generate a draft genome for a non-model species.

Limited indirect fitness benefits of male group membership in a lekking species.

In group living species, individuals may gain the indirect fitness benefits characterizing kin selection when groups contain close relatives. However, tests of kin selection have primarily focused on cooperatively breeding and eusocial species, whereas its importance in other forms of group living remains to be fully understood. Lekking is a form of grouping where males display on small aggregated territories, which females then visit to mate. As females prefer larger aggregations, territorial males might gain indirect fitness benefits if their presence increases the fitness of close relatives. Previous studies have tested specific predictions of kin selection models using measures such as group-level relatedness. However, a full understanding of the contribution of kin selection in the evolution of group living requires estimating individuals' indirect fitness benefits across multiple sites and years. Using behavioural and genetic data from the black grouse (Tetrao tetrix), we show that the indirect fitness benefits of group membership were very small because newcomers joined leks containing few close relatives who had limited mating success. Males' indirect fitness benefits were higher in yearlings during increasing population density but marginally changed the variation in male mating success. Kin selection acting through increasing group size is therefore unlikely to contribute substantially to the evolution and maintenance of lekking in this black grouse population.