Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches

Adaptive potential in the face of a transmissible cancer in Tasmanian devils

Emerging infectious diseases (EIDs) not only cause catastrophic declines in wildlife populations but also generate selective pressures that may result in rapid evolutionary responses. One such EID is devil facial tumour disease (DFTD) in the Tasmanian devil. DFTD is almost always fatal and has reduced the average lifespan of individuals by around 2 years, likely causing strong selection for traits that reduce susceptibility to the disease, but population decline has also left Tasmanian devils vulnerable to inbreeding depression. We analysed 22 years of data from an ongoing study of a population of Tasmanian devils on Freycinet Peninsula, Tasmania, to (1) identify whether DFTD may be causing selection on body size, by estimating phenotypic and genetic correlations between DFTD and size traits, (2) estimate the additive genetic variance of susceptibility to DFTD, and (3) investigate whether size traits or susceptibility to DFTD were under inbreeding depression. We found a positive phenotypic relationship between head width and susceptibility to DFTD, but this was not underpinned by a genetic correlation. Conversely, we found a negative phenotypic relationship between body weight and susceptibility to DFTD, and there was evidence for a negative genetic correlation between susceptibility to DFTD and body weight. There was additive genetic variance in susceptibility to DFTD, head width and body weight, but there was no evidence for inbreeding depression in any of these traits. These results suggest that Tasmanian devils have the potential to respond adaptively to DFTD, although the realised evolutionary response will critically further depend on the evolution of DFTD itself.

The heritability of fitness in a wild annual plant population with hierarchical size structure

The relative magnitude of additive genetic vs. residual variation for fitness traits is important in models for predicting the rate of evolution and population persistence in response to changes in the environment. In many annual plants, lifetime reproductive fitness is correlated with end-of-season plant biomass, which can vary significantly from plant to plant in the same population. We measured end-of-season plant biomasses and obtained single nucleotide polymorphism (SNP) genotypes of plants in a dense, natural population of the annual plant species Impatiens capensis with hierarchical size structure. These data were used to estimate the amount of heritable variation for position in the size hierarchy and for plant biomass. Additive genetic variance for a position in the size hierarchy and plant biomass were both significantly different from zero. These results are discussed in relationship to the theory for the heritability of fitness in natural populations and ecological factors that potentially influence heritable variation for fitness in this species.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

PERSEUS: an interactive and intuitive web-based tool for pedigree visualization

SUMMARY

Pedigree-based analyses' prime role is to unravel relationships between individuals in breeding programs and germplasms. This is critical information for decoding the genetics underlying main inherited traits of relevance, and unlocking the genotypic variability of a species to carry out genomic selections and predictions. Despite the great interest, current lineage visualizations become quite limiting in terms of public display, exploration, and tracing of traits up to ancestral donors. PERSEUS is a user-friendly, intuitive, and interactive web-based tool for pedigree visualizations represented as directed graph networks distributed using a force-repulsion method. The visualizations do not only showcase individual relationships among accessions, but also facilitate a seamless search and download of phenotypic traits along the pedigrees. PERSEUS is a promising tool for breeders and scientists, advantageous for evolutionary, genealogy, and diversity analyses among related accessions and species.

AVAILABILITY AND IMPLEMENTATION

PERSEUS is freely accessible at https://bioinformatics.cragenomica.es/perseus and GitHub code is available at https://github.com/aranzana-lab/PERSEUS.

Heritability and developmental plasticity of growth in an oviparous lizard

Selective processes act on phenotypic variation although the evolutionary potential of a trait relies on the underlying heritable variation. Developmental plasticity is an important source of phenotypic variation, but it can also promote changes in genetic variation, yet we have a limited understanding of how they are both impacted. Here, we quantified the influence of developmental temperature on growth in delicate skinks (Lampropholis delicata) and partitioned total phenotypic variance using an animal model fitted with a genomic relatedness matrix. We measured mass for 261 individuals (nhot = 125, ncold = 136) over 16 months (nobservations = 3002) and estimated heritability and maternal effects over time. Our results show that lizards reared in cold developmental temperatures had consistently higher mass across development compared to lizards that were reared in hot developmental temperatures. However, developmental temperature did not impact the rate of growth. On average, additive genetic variance, maternal effects and heritability were higher in the hot developmental temperature treatment; however, these differences were not statistically significant. Heritability increased with age, whereas maternal effects decreased upon hatching but increased again at a later age, which could be driven by social competition or intrinsic changes in the expression of variation as an individual's growth. Our work suggests that the evolutionary potential of growth is complex, age-dependent and not overtly affected by extremes in natural nest temperatures.

The estimation of additive genetic variance of body size in a wild passerine is sensitive to the method used to estimate relatedness among the individuals

Assessing additive genetic variance is a crucial step in predicting the evolutionary response of a target trait. However, the estimated genetic variance may be sensitive to the methodology used, e.g., the way relatedness is assessed among the individuals, especially in wild populations where social pedigrees can be inaccurate. To investigate this possibility, we investigated the additive genetic variance in tarsus length, a major proxy of skeletal body size in birds. The model species was the collared flycatcher (Ficedula albicollis), a socially monogamous but genetically polygamous migratory passerine. We used two relatedness matrices to estimate the genetic variance: (1) based solely on social links and (2) a genetic similarity matrix based on a large array of single-nucleotide polymorphisms (SNPs). Depending on the relatedness matrix considered, we found moderate to high additive genetic variance and heritability estimates for tarsus length. In particular, the heritability estimates were higher when obtained with the genetic similarity matrix instead of the social pedigree. Our results confirm the potential for this crucial trait to respond to selection and highlight methodological concerns when calculating additive genetic variance and heritability in phenotypic traits. We conclude that using a social pedigree instead of a genetic similarity matrix to estimate relatedness among individuals in a genetically polygamous wild population may significantly deflate the estimates of additive genetic variation.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Kinship composition in mammals

Understanding the evolution of group-living and cooperation requires information on who animals live and cooperate with. Animals can live with kin, non-kin or both, and kinship structure can influence the benefits and costs of group-living and the evolution of within-group cooperation. One aspect of kinship structure is kinship composition, i.e. a group-level attribute of the presence of kin and/or non-kin dyads in groups. Despite its putative importance, the kinship composition of mammalian groups has yet to be characterized. Here, we use the published literature to build an initial kinship composition dataset in mammals, laying the groundwork for future work in the field. In roughly half of the 18 species in our sample, individuals lived solely with same-sex kin, and, in the other half, individuals lived with related and unrelated individuals of the same sex. These initial results suggest that it is not rare for social mammals to live with unrelated individuals of the same sex, highlighting the importance of considering indirect and direct fitness benefits as co-drivers of the evolution of sociality. We hope that our initial dataset and insights will spur the study of kinship structure and sociality towards new exciting avenues.

Correlates of early reproduction and apparent fitness consequences in male Soay sheep

Life history trade-offs are ubiquitous across species and place constraints on the timing of life history events, including the optimal age at first reproduction. However, studies on lifetime breeding success of male mammals are rare due to sex-biased dispersal and the requirement for genetic paternity inferences. We studied the correlates and apparent fitness consequences of early life reproduction among males in a free-living population of Soay sheep (Ovis aries) on St Kilda, Scotland. We investigated the factors associated with early breeding success and the apparent consequences of early success for survival and future reproduction. We used genetic paternity inferences, population data, and individual morphology measurements collected over 30 years. We found that individuals born in years with low-density population size had the highest early life breeding success and singletons were more likely to be successful than twins. Individuals that bred successfully at 7 months were more likely to survive their first winter. For individuals that survived their first winter, early breeding success was not associated with later breeding success. As individual heterogeneity affects breeding success, we believe that variation in individual quality masks the costs of early reproduction in this population. Our findings provide no evidence for selection for delayed age at reproduction in male Soay sheep.

Quantitative genetic analysis of floral traits shows current limits but potential evolution in the wild

The vast variation in floral traits across angiosperms is often interpreted as the result of adaptation to pollinators. However, studies in wild populations often find no evidence of pollinator-mediated selection on flowers. Evolutionary theory predicts this could be the outcome of periods of stasis under stable conditions, followed by shorter periods of pollinator change that provide selection for innovative phenotypes. We asked if periods of stasis are caused by stabilizing selection, absence of other forms of selection or by low trait ability to respond even if selection is present. We studied a plant predominantly pollinated by one bee species across its range. We measured heritability and evolvability of traits, using genome-wide relatedness in a large wild population, and combined this with estimates of selection on the same individuals. We found evidence for both stabilizing selection and low trait heritability as potential explanations for stasis in flowers. The area of the standard petal is under stabilizing selection, but the variability is not heritable. A separate trait, floral weight, presents high heritability, but is not currently under selection. We show how a simple pollination environment coincides with the absence of current prerequisites for adaptive evolutionary change, while heritable variation remains to respond to future selection pressures.

A polygenic basis for birth weight in a wild population of red deer (Cervus elaphus)

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait.

Genomic heritability and correlation between carcass traits in Japanese Black cattle evaluated under different ceilings of relatedness among individuals

The investigation of carcass traits to produce meat with high efficiency has been in focus on Japanese Black cattle since 1972. To implement a successful breeding program in carcass production, a comprehensive understanding of genetic characteristics and relationships between the traits is of paramount importance. In this study, genomic heritability and genomic correlation between carcass traits, including carcass weight (CW), rib eye area (REA), rib thickness (RT), subcutaneous fat thickness (SFT), yield rate (YI), and beef marbling score (BMS) were estimated using the genomic data of 9,850 Japanese Black cattle (4,142 heifers and 5,708 steers). In addition, we investigated the effect of genetic relatedness degree on the estimation of genetic parameters of carcass traits in sub-populations created based on different GRM-cutoff values. Genome-based restricted maximum likelihood (GREML) analysis was applied to estimate genetic parameters. Using all animal data, the heritability values for carcass traits were estimated as moderate to relatively high magnitude, ranging from 0.338 to 0.509 with standard errors, ranging from 0.014 to 0.015. The genetic correlations were obtained low and negative between SFT and REA [-0.198 (0.034)] and between SFT and BMS [-0.096 (0.033)] traits, and high and negative between SFT and YI [-0.634 (0.022)]. REA trait was genetically highly correlated with YI and BMS [0.811 (0.012) and 0.625 (0.022), respectively]. In sub-populations created based on the genetic-relatedness ceiling, the heritability estimates ranged from 0.212 (0.131) to 0.647 (0.066). At the genetic-relatedness ceiling of 0.15, the correlation values between most traits with low genomic correlation were overestimated while the correlations between the traits with relatively moderate to high correlations, ranging from 0.380 to 0.811, were underestimated. The values were steady at the ceilings of 0.30-0.95 (sample size of 5,443-9,850) for most of the highly correlated traits. The results demonstrated that there is considerable genetic variation and also favorable genomic correlations between carcass traits. Therefore, the genetic improvement for the traits can be simultaneously attained through genomic selection. In addition, we observed that depending on the degree of relationship between individuals and sample size, the genomic heritability and correlation estimates for carcass traits may be different.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

Taking flight: An ecological, evolutionary and genomic perspective on bat telomeres

Over 20% of all living mammals are bats (order Chiroptera). Bats possess extraordinary adaptations including powered flight, laryngeal echolocation and a unique immune system that enables them to tolerate a diversity of viral infections without presenting clinical disease symptoms. They occupy multiple trophic niches and environments globally. Significant physiological and ecological diversity occurs across the order. Bats also exhibit extreme longevity given their body size with many species showing few signs of ageing. The molecular basis of this extended longevity has recently attracted attention. Telomere maintenance potentially underpins bats' extended healthspan, although functional studies are still required to validate the causative mechanisms. In this review, we detail the current knowledge on bat telomeres, telomerase expression, and how these relate to ecology, longevity and life-history strategies. Patterns of telomere shortening and telomerase expression vary across species, and comparative genomic analyses suggest that alternative telomere maintenance mechanisms evolved in the longest-lived bats. We discuss the unique challenges faced when working with populations of wild bats and highlight ways to advance the field including expanding long-term monitoring across species that display contrasting life-histories and occupy different environmental niches. We further review how new high quality, chromosome-level genome assemblies can enable us to uncover the molecular mechanisms governing telomere dynamics and how phylogenomic analyses can reveal the adaptive significance of telomere maintenance and variation in bats.

The impact of SNP density on quantitative genetic analyses of body size traits in a wild population of Soay sheep

Understanding the genetic architecture underpinning quantitative traits in wild populations is pivotal to understanding the processes behind trait evolution. The 'animal model' is a popular method for estimating quantitative genetic parameters such as heritability and genetic correlation and involves fitting an estimate of relatedness between individuals in the study population. Genotypes at genome-wide markers can be used to estimate relatedness; however, relatedness estimates vary with marker density, potentially affecting results. Increasing density of markers is also expected to increase the power to detect quantitative trait loci (QTL). In order to understand how the density of genetic markers affects the results of quantitative genetic analyses, we estimated heritability and performed genome-wide association studies (GWAS) on five body size traits in an unmanaged population of Soay sheep using two different SNP densities: a dataset of 37,037 genotyped SNPs and an imputed dataset of 417,373 SNPs. Heritability estimates did not differ between the two SNP densities, but the high-density imputed SNP dataset revealed four new SNP-trait associations that were not found with the lower density dataset, as well as confirming all previously-found QTL. We also demonstrated that fitting fixed and random effects in the same step as performing GWAS is a more powerful approach than pre-correcting for covariates in a separate model.

Genomic prediction in the wild: A case study in Soay sheep

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

Telomere heritability and parental age at conception effects in a wild avian population

Individual variation in telomere length is predictive of health and mortality risk across a range of species. However, the relative influence of environmental and genetic variation on individual telomere length in wild populations remains poorly understood. Heritability of telomere length has primarily been calculated using parent-offspring regression which can be confounded by shared environments. To control for confounding variables, quantitative genetic "animal models" can be used, but few studies have applied animal models in wild populations. Furthermore, parental age at conception may also influence offspring telomere length, but most studies have been cross-sectional. We investigated within- and between-parental age at conception effects and heritability of telomere length in the Seychelles warbler using measures from birds caught over 20 years and a multigenerational pedigree. We found a weak negative within-paternal age at conception effect (as fathers aged, their offspring had shorter telomeres) and a weak positive between-maternal age at conception effect (females that survived to older ages had offspring with longer telomeres). Animal models provided evidence that heritability and evolvability of telomere length were low in this population, and that variation in telomere length was not driven by early-life effects of hatch period or parental identities. Quantitative polymerase chain reaction plate had a large influence on telomere length variation and not accounting for it in the models would have underestimated heritability. Our study illustrates the need to include and account for technical variation in order to accurately estimate heritability, as well as other environmental effects, on telomere length in natural populations.

The Unusual Value of Long-Term Studies of Individuals: The Example of the Isle of Rum Red Deer Project

Long-term studies of individuals enable incisive investigations of questions across ecology and evolution. Here, we illustrate this claim by reference to our long-term study of red deer on the Isle of Rum, Scotland. This project has established many of the characteristics of social organization, selection, and population ecology typical of large, polygynous, seasonally breeding mammals, with wider implications for our understanding of sexual selection and the evolution of sex differences, as well as for their population dynamics and population management. As molecular genetic techniques have developed, the project has pivoted to investigate evolutionary genetic questions, also breaking new ground in this field. With ongoing advances in genomics and statistical approaches and the development of increasingly sophisticated ways to assay new phenotypic traits, the questions that long-term studies such as the red deer study can answer become both broader and ever more sophisticated. They also offer powerful means of understanding the effects of ongoing climate change on wild populations.

Estimating heritability of song considering within-individual variance in a wild songbird: The collared flycatcher

Heritable genetic variation is a prerequisite for adaptive evolution; however, our knowledge about the heritability of plastic traits, such as behaviors, is scarce, especially in wild populations. In this study, we investigated the heritability of song traits in the collared flycatcher (Ficedula albicollis), a small oscine passerine with complex songs involved in sexual selection. We recorded the songs of 81 males in a natural population and obtained various measures describing the frequency, temporal organization, and complexity of each song. As we had multiple songs from each individual, we were able to statistically account for the first time for the effect of within-individual variance on the heritability of song. Heritability was calculated from the variance estimates of animal models relying on a genetic similarity matrix based on Single Nucleotide Polymorphism screening. Overall, we found small additive genetic variance and heritability values in all song traits, highlighting the role of environmental factors in shaping bird song.

A meta-analysis on the heritability of vertebrate telomere length

Telomere dynamics are linked with both cellular and organismal senescence, and life history, individual quality and health. Telomere dynamics, particularly telomere length, have therefore garnered much research interest in evolutionary biology. To examine the evolution of telomere length, it is important to quantify its heritability, the proportion of total variation explained by additive genetic effects. Many studies have quantified telomere length heritability, but estimates are varied, and no general conclusion has been drawn. Additionally, it is unclear whether biological and methodological factors influence telomere length heritability estimates. We present the first meta-analysis of telomere length heritability, using 104 estimates from 43 studies over 18 vertebrate species. We calculated an overall mean heritability and examined how estimates varied by study, phylogeny, species-specific ecology, environmental setting, age at sampling, laboratory methods, statistical methods, sex and repeated measurements. Overall heritability was moderate (44.9%, 95% CI: 25.2-64.7%), and there was considerable heterogeneity in heritability estimates, in particular among studies and estimates. Laboratory method influenced heritability estimates, with in-gel hybridization TRF yielding higher heritabilities than qPCR and Southern blot TRF. There was also an effect from statistical method, with twin-based and SNP-based estimates lower than correlation-based or pedigree-based estimates. Our results highlight an overall heritable basis of telomere length, and we recommend future research on a wider range of taxa, and the use of variance-partitioning methods with relatedness or SNP data over correlation methods to minimize heritability estimation bias.

Vitamin D status is heritable and under environment-dependent selection in the wild

Vitamin D has a well-established role in skeletal health and is increasingly linked to chronic disease and mortality in humans and companion animals. Despite the clear significance of vitamin D for health and obvious implications for fitness under natural conditions, no longitudinal study has tested whether the circulating concentration of vitamin D is under natural selection in the wild. Here, we show that concentrations of dietary-derived vitamin D2 and endogenously produced vitamin D3  metabolites are heritable and largely polygenic in a wild population of Soay sheep (Ovis aries). Vitamin D2  status was positively associated with female adult survival, and vitamin D3  status predicted female fecundity in particular, good environment years when sheep density and competition for resources was low. Our study provides evidence that vitamin D status has the potential to respond to selection, and also provides new insights into how vitamin D metabolism is associated with fitness in the wild.

Genetic and environmental drivers of colour and pattern in the Australian jacky dragon (Amphibolurus muricatus)

The underlying drivers of variation in the colouration (colour and pattern) of animals can be genetic, non‐genetic, or more likely, a combination of both. Understanding the role of heritable genetic elements, as well as non‐genetic factors such as age, habitat or temperature, in shaping colouration can provide insight into the evolution and function of these traits, as well as the speed of response to changing environments. This project examined the genetic and non‐genetic drivers of continuous variation in colouration in a lizard, the jacky dragon (Amphibolurus muricatus). We leveraged a large captive experiment that manipulated parental and offspring thermal environment to simultaneously estimate the genetic and non‐genetic drivers of variation in colouration. We found that the overall brightness, the elongation of the longitudinal stripes on the dorsum and the contrast between light and dark patches of the pattern were all heritable. Colouration varied according to the age of the hatchling; however, the thermal environment of neither the parents nor offspring contributed significantly to colouration. It appears that developmental plasticity and maternal effects associated with temperature are not important drivers of variation in our measures of colouration.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Using genomic prediction to detect microevolutionary change of a quantitative trait

Detecting microevolutionary responses to natural selection by observing temporal changes in individual breeding values is challenging. The collection of suitable datasets can take many years and disentangling the contributions of the environment and genetics to phenotypic change is not trivial. Furthermore, pedigree-based methods of obtaining individual breeding values have known biases. Here, we apply a genomic prediction approach to estimate breeding values of adult weight in a 35-year dataset of Soay sheep (Ovis aries). Comparisons are made with a traditional pedigree-based approach. During the study period, adult body weight decreased, but the underlying genetic component of body weight increased, at a rate that is unlikely to be attributable to genetic drift. Thus cryptic microevolution of greater adult body weight has probably occurred. Genomic and pedigree-based approaches gave largely consistent results. Thus, using genomic prediction to study microevolution in wild populations can remove the requirement for pedigree data, potentially opening up new study systems for similar research.

Estimating heritability using family-pooled phenotypic and genotypic data: a simulation study applied to aquaculture

Estimating heritability based on individual phenotypic and genotypic measurements can be expensive and labour-intensive in commercial aquaculture breeding. Here, the feasibility of estimating heritability using within-family means of phenotypes and allelic frequencies was investigated. Different numbers of full-sib families and family sizes across ten generations with phenotypic and genotypic information on 10 K SNPs were analysed in ten replicates. Three scenarios, representing differing numbers of pools per family (one, two and five) were considered. The results showed that using one pool per family did not reliably estimate the heritability of family means. Using simulation parameters appropriate for aquaculture, at least 200 families of 60 progeny per family divided equally in two pools per family was required to estimate the heritability of family means effectively. Although application of five pools generated more within- and between- family relationships, it reduced the number of individuals per pool and increased within-family residual variation, hence, decreased the heritability of family means. Moreover, increasing the size of pools resulted in increasing the heritability of family means towards one. In addition, heritability of family mean estimates were higher than family heritabilities obtained from Falconer’s formula due to lower intraclass correlation estimate compared to the coefficient of relationship.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

Inheritance of Yield Components and Morphological Traits in Avocado cv. Hass From "Criollo" "Elite Trees" via Half-Sib Seedling Rootstocks

Grafting induces precocity and maintains clonal integrity in fruit tree crops. However, the complex rootstock × scion interaction often precludes understanding how the tree phenotype is shaped, limiting the potential to select optimum rootstocks. Therefore, it is necessary to assess (1) how seedling progenies inherit trait variation from elite 'plus trees', and (2) whether such family superiority may be transferred after grafting to the clonal scion. To bridge this gap, we quantified additive genetic parameters (i.e., narrow sense heritability-h ², and genetic-estimated breeding values-GEBVs) across landraces, "criollo", "plus trees" of the super-food fruit tree crop avocado (Persea americana Mill.), and their open-pollinated (OP) half-sib seedling families. Specifically, we used a genomic best linear unbiased prediction (G-BLUP) model to merge phenotypic characterization of 17 morpho-agronomic traits with genetic screening of 13 highly polymorphic SSR markers in a diverse panel of 104 avocado "criollo" "plus trees." Estimated additive genetic parameters were validated at a 5-year-old common garden trial (i.e., provenance test), in which 22 OP half-sib seedlings from 82 elite "plus trees" served as rootstocks for the cv. Hass clone. Heritability (h ²) scores in the "criollo" "plus trees" ranged from 0.28 to 0.51. The highest h ² values were observed for ribbed petiole and adaxial veins with 0.47 (CI 95%0.2-0.8) and 0.51 (CI 0.2-0.8), respectively. The h ² scores for the agronomic traits ranged from 0.34 (CI 0.2-0.6) to 0.39 (CI 0.2-0.6) for seed weight, fruit weight, and total volume, respectively. When inspecting yield variation across 5-year-old grafted avocado cv. Hass trees with elite OP half-sib seedling rootstocks, the traits total number of fruits and fruits' weight, respectively, exhibited h ² scores of 0.36 (± 0.23) and 0.11 (± 0.09). Our results indicate that elite "criollo" "plus trees" may serve as promissory donors of seedling rootstocks for avocado cv. Hass orchards due to the inheritance of their outstanding trait values. This reinforces the feasibility to leverage natural variation from "plus trees" via OP half-sib seedling rootstock families. By jointly estimating half-sib family effects and rootstock-mediated heritability, this study promises boosting seedling rootstock breeding programs, while better discerning the consequences of grafting in fruit tree crops.

Genomic analysis reveals a polygenic architecture of antler morphology in wild red deer (Cervus elaphus)

Sexually selected traits show large variation and rapid evolution across the animal kingdom, yet genetic variation often persists within populations despite apparent directional selection. A key step in solving this long-standing paradox is to determine the genetic architecture of sexually selected traits to understand evolutionary drivers and constraints at the genomic level. Antlers are a form of sexual weaponry in male red deer (Cervus elaphus). On the island of Rum, Scotland, males with larger antlers have increased breeding success, yet there has been no evidence of any response to selection at the genetic level. To try and understand the mechanisms underlying this observation, we investigate the genetic architecture of ten antler traits and their principal components using genomic data from >38,000 SNPs. We estimate the heritabilities and genetic correlations of the antler traits using a genomic relatedness approach. We then use genome-wide association and haplotype-based regional heritability to identify regions of the genome underlying antler morphology, and an empirical Bayes approach to estimate the underlying distributions of allele effect sizes. We show that antler morphology is highly repeatable over an individual's lifetime, heritable and has a polygenic architecture and that almost all antler traits are positively genetically correlated with some loci identified as having pleiotropic effects. Our findings suggest that a large mutational target and genetic covariances among antler traits, in part maintained by pleiotropy, are likely to contribute to the maintenance of genetic variation in antler morphology in this population.

MHC class IIa haplotypes derived by high-throughput SNP screening in an isolated sheep population

Investigating the current evolutionary processes acting on a highly polymorphic gene region, such as the major histocompatibility complex (MHC), requires extensive population data for both genotypes and phenotypes. The MHC consists of several tightly linked loci with both allelic and gene content variation, making it challenging to genotype. Eight class IIa haplotypes have previously been identified in the Soay sheep (Ovis aries) of St. Kilda using Sanger sequencing and cloning, but no single locus is representative of all haplotypes. Here, we exploit the closed nature of the island population of Soay sheep and its limited haplotypic variation to identify a panel of SNPs that enable imputation of MHC haplotypes. We compared MHC class IIa haplotypes determined by Sanger sequence-based genotyping of 135 individuals to their SNP profiles generated using the Ovine Infinium HD BeadChip. A panel of 11 SNPs could reliably determine MHC diplotypes, and two additional SNPs within the DQA1 gene enabled detection of a recombinant haplotype affecting only the SNPs downstream of the expressed genes. The panel of 13 SNPs was genotyped in 5951 Soay sheep, of which 5349 passed quality control. Using the Soay sheep pedigree, we were able to trace the origin and inheritance of the recombinant SNP haplotype. This SNP-based method has enabled the rapid generation of locus-specific MHC genotypes for large numbers of Soay sheep. This volume of high-quality genotypes in a well-characterized population of free-living sheep will be valuable for investigating the mechanisms maintaining diversity at the MHC.

How is adaptive potential distributed within species ranges?

Quantitative genetic variation (QGV) represents a major component of adaptive potential and, if reduced toward range-edge populations, could prevent a species' expansion or adaptive response to rapid ecological change. It has been hypothesized that QGV will be lower at the range edge due to small populations-often the result of poor habitat quality-and potentially decreased gene flow. However, whether central populations are higher in QGV is unknown. We used a meta-analytic approach to test for a general QGV-range position relationship, including geographic and climatic distance from range centers. We identified 35 studies meeting our criteria, yielding nearly 1000 estimates of QGV (including broad-sense heritability, narrow-sense heritability, and evolvability) from 34 species. The relationship between QGV and distance from the geographic range or climatic niche center depended on the focal trait and how QGV was estimated. We found some evidence that QGV declines from geographic centers but that it increases toward niche edges; niche and geographic distances were uncorrelated. Nevertheless, few studies have compared QGV in both central and marginal regions or environments within the same species. We call for more research in this area and discuss potential research avenues related to adaptive potential in the context of global change.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

Patterns of MHC-dependent sexual selection in a free-living population of sheep

The MHC is one of the most polymorphic gene clusters in vertebrates and play an essential role in adaptive immunity. Apart from pathogen-mediated selection, sexual selection can also contribute to the maintenance of MHC diversity. MHC-dependent sexual selection could occur via several mechanisms but at present there is no consensus as to which of these mechanisms are involved and their importance. Previous studies have often suffered from limited genetic and behavioural data and small sample size, and were rarely able to examine all the mechanisms together, determine whether signatures of MHC-based non-random mating are independent of genomic effects or differentiate whether MHC-dependent sexual selection takes place at the pre- or post-copulatory stage. In this study, we use Monte Carlo simulation to investigate evidence for non-random MHC-dependent mating patterns by all three mechanisms in a free-living population of Soay sheep. Using 1710 sheep diplotyped at the MHC class IIa region and genome-wide SNPs, together with field observations of consorts, we found sexual selection against a particular haplotype in males at the pre-copulatory stage and sexual selection against female MHC heterozygosity during the rut. We also found MHC-dependent disassortative mating at the post-copulatory stage, along with strong evidence of inbreeding avoidance at both stages. However, results from generalized linear mixed models suggest that the pattern of MHC-dependent disassortative mating could be a by-product of inbreeding avoidance. Our results therefore suggest that while multiple apparent mechanisms of non-random mating with respect to the MHC may occur, some of them have alternative explanations.

Sample identification and pedigree reconstruction in Wolverine (Gulo gulo) using SNP genotyping of non-invasive samples

For conservation genetic studies using non-invasively collected samples, genome-wide data may be hard to acquire. Until now, such studies have instead mostly relied on analyses of traditional genetic markers such as microsatellites (SSRs). Recently, high throughput genotyping of single nucleotide polymorphisms (SNPs) has become available, expanding the use of genomic methods to include non-model species of conservation concern. We have developed a 96-marker SNP array for use in applied conservation monitoring of the Scandinavian wolverine (Gulo gulo) population. By genotyping more than a thousand non-invasively collected samples, we were able to obtain precise estimates of different types of genotyping errors and sample dropout rates. The SNP panel significantly outperforms the SSR markers (and DBY intron markers for sexing) both in terms of precision in genotyping, sex assignment and individual identification, as well as in the proportion of samples successfully genotyped. Furthermore, SNP genotyping offers a simplified laboratory and analysis pipeline with fewer samples needed to be repeatedly genotyped in order to obtain reliable consensus data. In addition, we utilised a unique opportunity to successfully demonstrate the application of SNP genotype data for reconstructing pedigrees in wild populations, by validating the method with samples from wild individuals with known relatedness. By offering a simplified workflow with improved performance, we anticipate this methodology will facilitate the use of non-invasive samples to improve genetic management of many different types of populations that have previously been challenging to survey.

Heritable variation in telomere length predicts mortality in Soay sheep

Telomere length (TL) is considered an important biomarker of whole-organism health and aging. Across humans and other vertebrates, short telomeres are associated with increased subsequent mortality risk, but the processes responsible for this correlation remain uncertain. A key unanswered question is whether TL-mortality associations arise due to positive effects of genes or early-life environment on both an individual's average lifetime TL and their longevity, or due to more immediate effects of environmental stressors on within-individual TL loss and increased mortality risk. Addressing this question requires longitudinal TL and life history data across the entire lifetimes of many individuals, which are difficult to obtain for long-lived species like humans. Using longitudinal data and samples collected over nearly two decades, as part of a long-term study of wild Soay sheep, we dissected an observed positive association between TL and subsequent survival using multivariate quantitative genetic models. We found no evidence that telomere attrition was associated with increased mortality risk, suggesting that TL is not an important marker of biological aging or exposure to environmental stress in our study system. Instead, we find that among-individual differences in average TL are associated with increased lifespan. Our analyses suggest that this correlation between an individual's average TL and lifespan has a genetic basis. This demonstrates that TL has the potential to evolve under natural conditions, and suggests an important role of genetics underlying the widespread observation that short telomeres predict mortality.

Within-trio tests provide little support for post-copulatory selection on major histocompatibility complex haplotypes in a free-living population

Sexual selection has been proposed as a force that could help maintain the diversity of major histocompatibility complex (MHC) genes in vertebrates. Potential selective mechanisms can be divided into pre-copulatory and post-copulatory, and in both cases, the evidence for occurrence is mixed, especially in natural populations. In this study, we used a large number of parent-offspring trios that were diplotyped for MHC class II genes in a wild population of Soay sheep (Ovis aries) to examine whether there was within-trio post-copulatory selection on MHC class II genes at both the haplotype and diplotype levels. We found there was transmission ratio distortion of one of the eight MHC class II haplotypes (E) which was transmitted less than expected by fathers, and transmission ratio distortion of another haplotype (A) which was transmitted more than expected by chance to male offspring. However, in both cases, these deviations were not significant after correction for multiple tests. In addition, we did not find any evidence of post-copulatory selection at the diplotype level. These results imply that, given known parents, there is no strong post-copulatory selection on MHC class II genes in this population.

Inheritance of Rootstock Effects in Avocado (Persea americana Mill.) cv. Hass

Grafting is typically utilized to merge adapted seedling rootstocks with highly productive clonal scions. This process implies the interaction of multiple genomes to produce a unique tree phenotype. However, the interconnection of both genotypes obscures individual contributions to phenotypic variation (rootstock-mediated heritability), hampering tree breeding. Therefore, our goal was to quantify the inheritance of seedling rootstock effects on scion traits using avocado (Persea americana Mill.) cv. Hass as a model fruit tree. We characterized 240 diverse rootstocks from 8 avocado cv. Hass orchards with similar management in three regions of the province of Antioquia, northwest Andes of Colombia, using 13 microsatellite markers simple sequence repeats (SSRs). Parallel to this, we recorded 20 phenotypic traits (including morphological, biomass/reproductive, and fruit yield and quality traits) in the scions for 3 years (2015-2017). Relatedness among rootstocks was inferred through the genetic markers and inputted in a "genetic prediction" model to calculate narrow-sense heritabilities (h 2) on scion traits. We used three different randomization tests to highlight traits with consistently significant heritability estimates. This strategy allowed us to capture five traits with significant heritability values that ranged from 0.33 to 0.45 and model fits (r) that oscillated between 0.58 and 0.73 across orchards. The results showed significance in the rootstock effects for four complex harvest and quality traits (i.e., total number of fruits, number of fruits with exportation quality, and number of fruits discarded because of low weight or thrips damage), whereas the only morphological trait that had a significant heritability value was overall trunk height (an emergent property of the rootstock-scion interaction). These findings suggest the inheritance of rootstock effects, beyond root phenotype, on a surprisingly wide spectrum of scion traits in "Hass" avocado. They also reinforce the utility of polymorphic SSRs for relatedness reconstruction and genetic prediction of complex traits. This research is, up to date, the most cohesive evidence of narrow-sense inheritance of rootstock effects in a tropical fruit tree crop. Ultimately, our work highlights the importance of considering the rootstock-scion interaction to broaden the genetic basis of fruit tree breeding programs while enhancing our understanding of the consequences of grafting.

Heritability Estimates of Antler and Body Traits in White-Tailed Deer (Odocoileus virginianus) From Genomic-Relatedness Matrices

Estimating heritability (h2) is required to predict the response to selection and is useful in species that are managed or farmed using trait information. Estimating h2 in free-ranging populations is challenging due to the need for pedigrees; genomic-relatedness matrices (GRMs) circumvent this need and can be implemented in nearly any system where phenotypic and genome-wide single-nucleotide polymorphism (SNP) data are available. We estimated the heritability of 5 body and 3 antler traits in a free-ranging population of white-tailed deer (Odocoileus virginianus) on Anticosti Island, Quebec, Canada. We generated classic and robust GRMs from >10,000 SNPs: hind foot length, dressed body mass, and peroneus muscle mass had high h2 values of 0.62, 0.44, and 0.55, respectively. Heritability in male-only antler features ranged from 0.07 to 0.33. We explored the influence of filtering by minor allele frequency and data completion on h2: GRMs derived from fewer SNPs had reduced h2 estimates and the relatedness coefficients significantly deviated from those generated with more SNPs. As a corollary, we discussed limitations to the application of GRMs in the wild, notably how skewed GRMs, specifically many unrelated individuals, can increase variance around h2 estimates. This is the first study to estimate h2 on a free-ranging population of white-tailed deer and should be informative for breeding designs and management as these traits could respond to selection.

Polygenic basis for adaptive morphological variation in a threatened Aotearoa | New Zealand bird, the hihi (Notiomystis cincta)

To predict if a threatened species can adapt to changing selective pressures, it is crucial to understand the genetic basis of adaptive traits, especially in species historically affected by severe bottlenecks. We estimated the heritability of three hihi (Notiomystis cincta) morphological traits known to be under selection (nestling tarsus length, body mass and head-bill length) using 523 individuals and 39 699 single nucleotide polymorphisms (SNPs) from a 50 K Affymetrix SNP chip. We then examined the genetic architecture of the traits via chromosome partitioning analyses and genome-wide association scans (GWAS). Heritabilities estimated using pedigree relatedness or genomic relatedness were low. For tarsus length, the proportion of genetic variance explained by each chromosome was positively correlated with its size, and more than one chromosome explained significant variation for body mass and head-bill length. Finally, GWAS analyses suggested many loci of small effect contributing to trait variation for all three traits, although one locus (an SNP within an intron of the transcription factor HEY2) was tentatively associated with tarsus length. Our findings suggest a polygenic nature for the morphological traits, with many small effect size loci contributing to the majority of the variation, similar to results from many other wild populations. However, the small effective population size, polygenic architecture and already low heritabilities suggest that both the total response and rate of response to selection are likely to be limited in hihi.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.

SNP-based heritability and genetic architecture of cranial cruciate ligament rupture in Labrador Retrievers

Cranial cruciate ligament rupture (CCLR) is one of the leading causes of pelvic limb lameness in dogs. About 6% of Labrador Retrievers suffer from this orthopedic problem. The aim of this study was to determine the heritability of CCLR in this breed using SNP array genotyping data. DNA samples were collected from CCLR-affected dogs (n = 190) and unaffected dogs over the age of 8 years (n = 143). All 333 dogs were genotyped directly or imputed up to approximately 710k SNPs on the Affymetrix Axiom CanineHD SNP array. Heritability of CCLR was calculated using multiple methodologies, including linear mixed models, Bayesian models and a model that incorporates LD. The covariates of sex and sterilization status were added to each analysis to assess their impact. Across the algorithms of these models, heritability ranged from 0.550 to 0.886, depending on covariate inclusion. The relatively high heritability for this disease indicates that a substantial genetic component contributes to CCLR in the Labrador Retriever.

Pedigree-Based Estimation of Reproductive Value

How successful an individual or cohort is, in terms of their genetic contribution to the future population, is encapsulated in the concept of reproductive value, and is crucial for understanding selection and evolution. Long-term studies of pedigreed populations offer the opportunity to estimate reproductive values directly. However, the degree to which genetic contributions, as defined by a pedigree, may converge on their long-run values within the time frames of available data sets, such that they may be interpreted as estimates of reproductive value, is unclear. We develop a system for pedigree-based calculation of the expected genetic representation that both individuals and cohorts make to the population in the years following their birth. We apply this system to inference of individual and cohort reproductive values in Soay sheep (Ovis aries) from St Kilda, Outer Hebrides. We observe that these genetic contributions appear to become relatively stable within modest time frames. As such, it may be reasonable to consider pedigree-based calculations of genetic contributions to future generations as estimates of reproductive value. This approach and the knowledge that the estimates can stabilize within decades should offer new opportunities to analyze data from pedigreed wild populations, which will be of value to many fields within evolutionary biology and demography.