Cell‐mediated immunity and multi‐locus heterozygosity in bluethroat nestlings

Recent evidence suggests that marker‐based heterozygosity‐fitness correlations may be driven by only one or a few markers, indicating local heterozygosity effects caused by linkage disequilibrium with functional genes. In this study, we investigated the relationship between microsatellite heterozygosity and a measure of cell‐mediated immunity (phytohaemagglutinin; PHA) in bluethroat (Luscinia s. svecica) nestlings using a full‐sibling design. We found significant positive associations between PHA response and two different indices of microsatellite heterozygosity, i.e. multi‐locus heterozygosity and mean d². However, model comparisons disclosed that both associations were more likely caused by local effects rather than general effects and that the two local effects appeared to be realized through two different genetic mechanisms. Our results indicate that both the random assortment of parental chromosomes during meiosis as well as inbreeding can drive heterozygosity‐fitness correlations.     

Comparing molecular measures for detecting inbreeding depression

Correlations between heterozygosity and components of fitness have been investigated in natural populations for over 20 years. Positive correlations between a trait of interest and heterozygosity (usually measured at allozyme loci) are generally recognized as evidence of inbreeding depression. More recently, molecular markers such as microsatellites have been employed for the same purpose. A typical study might use around five to ten markers. In this paper we use a panel of 71 microsatellite loci to: (1) Compare the efficacy of heterozygosity and a related microsatellite‐specific variable, mean d², in detecting inbreeding depression; (2) Examine the statistical power of heterozygosity to detect such associations. We performed our analyses in a wild population of red deer (Cervus elaphus) in which inbreeding depression in juvenile traits had previously been detected using a panel of nine markers. We conclude that heterozygosity‐based measures outperform mean d²‐based measures, but that power to detect heterozygosity‐fitness associations is nonetheless low when ten or fewer markers are typed.     

Population differences in behaviour are explained by shared within‐population trait correlations

Correlations in behavioural traits across time, situation and ecological context (i.e. ‘behavioural syndromes’ or ‘personality’) have been documented for a variety of behaviours, and in diverse taxa. Perhaps the most controversial inference from the behavioural syndromes literature is that correlated behaviour may act as an evolutionary constraint and evolutionary change in one’s behaviour may necessarily involve shifts in others. We test the two predictions of this hypothesis using comparative data from eighteen populations of the socially polymorphic spider, Anelosimus studiosus (Araneae, Theriidae). First, we ask whether geographically distant populations share a common syndrome. Second, we test whether population differences in behaviour are correlated similarly to within‐population trait correlations. Our results reveal that populations separated by as much as 36° latitude shared similar syndromes. Furthermore, population differences in behaviour were correlated in the same manner as within‐population trait correlations. That is, population divergence tended to be along the same axes as within‐population covariance. Together, these results suggest a lack of evolutionary independence in the syndrome’s constituent traits.     

Environmental effects on the covariation among pace‐of‐life traits

Pace‐of‐life syndromes (POLSs) are suites of life‐history, physiological and behavioural traits that arise due to trade‐offs between allocation to current and future reproduction. Traits generally show covariation that can arise from genetic and environmental influences on phenotypes and constrain the independent evolution of traits, resulting in fitness consequences and impacts on population dynamics. The notion that correlations among traits may vary among populations along environmental gradients suggests an important role for the environment in shaping and maintaining POLSs. However, no synthesis has been attempted of the myriad ways in which environmental factors should influence POLSs. Here, we formulate a series of hypotheses targeting the critical interfaces of the environment and life‐history ‐ behaviour associations across different organisms. We discuss the hypotheses in light of findings from a systematic review of studies that measured changes in the association between behaviour and life‐history traits as a function of environmental conditions. The review revealed that POLSs are often shaped by environmental variation, where harshness of the environment in early life has the most consistent effects on POLS. However, only partial or no effects of environmental variation were found in a number of studies, which may result from the highly variable study systems, traits and environments studied. We highlight promising directions arising from the available studies and identify knowledge gaps that, if unaddressed, will impede progress in the field.     

Do microsatellites reflect genome‐wide genetic diversity in natural populations? A comment on Väli et al. (2008)

A recent study by Väli et al. (2008) highlights that microsatellites will often provide a poor prediction of the genome‐wide nucleotide diversity of wild populations, but does not fully explain why. To clarify and stress the importance of identity disequilibrium and marker variability for correlations between multilocus heterozygosity and genome‐wide genetic variability, we performed a simple simulation with different types of markers, corresponding to microsatellites and SNPs, in populations with different inbreeding history. The importance of identity disequilibrium was apparent for both markers and there was a clear impact of marker variability.     

Assessment of identity disequilibrium and its relation to empirical heterozygosity fitness correlations: a meta‐analysis

Heterozygosity fitness correlations (HFCs) have frequently been used to detect inbreeding depression, under the assumption that genome‐wide heterozygosity is a good proxy for inbreeding. However, meta‐analyses of the association between fitness measures and individual heterozygosity have shown that often either no correlations are observed or the effect sizes are small. One of the reasons for this may be the absence of variance in inbreeding, a requisite for generating general‐effect HFCs. Recent work has highlighted identity disequilibrium (ID) as a measure that may capture variance in the level of inbreeding within a population; however, no thorough assessment of ID in natural populations has been conducted. In this meta‐analysis, we assess the magnitude of ID (as measured by the g₂ statistic) from 50 previously published HFC studies and its relationship to the observed effect sizes of those studies. We then assess how much power the studies had to detect general‐effect HFCs, and the number of markers that would have been needed to generate a high expected correlation (r² = 0.9) between observed heterozygosity and inbreeding. Across the majority of studies, g₂ values were not significantly different than zero. Despite this, we found that the magnitude of g₂ was associated with the average effect sizes observed in a population, even when point estimates were nonsignificant. These low values of g₂ translated into low expected correlations between heterozygosity and inbreeding and suggest that many more markers than typically used are needed to robustly detect HFCs.     

Seven polymorphic microsatellite DNA loci from the red‐spotted newt (Notophthalmus viridescens)

We describe polymerase chain reaction (PCR) primers and amplification conditions for seven microsatellite DNA loci isolated from the red‐spotted newt (Notophthalmus viridescens). Primers were tested on 16 individuals from two populations on the Savannah River Site in Aiken County, South Carolina. We detected six to 10 alleles per locus and an overall observed heterozygosity range of 0.31–0.81. Despite low heterozygosity at two of the seven loci, the high polymorphic information contents (from 0.54 to 0.85) of these markers render them useful for future studies of the behavioural and population ecology of this common salamander.     

Can a simple algebraic analysis predict markers-genome heterozygosity correlations?

A current algebraic analysis on genome-wide heterozygosity estimates suggests that correlations between molecular markers and genome-wide heterozygosity, rho, depend on the ratio between the number of markers used, r, and the number of genome loci, n; that is: rho approximately square root r/n. Hence, it is unfeasible to obtain reliable estimates of genome-wide heterozygosity in species of large genome using a few markers. We cast some doubts about this analysis as it assumed that the probability that an individual was heterozygous at a locus is equal to the average heterozygosity of this locus in the population. However, we believe that individual heterozygosity at a given locus depends on individual pedigree. Because the pedigree is common for all loci of an individual, their probabilities of heterozygosity are not independent within the genome. We first performed simulations generating random genomes for 100 individuals. Among these individuals, markers and genome-wide heterozygosities correlated as expected from the above equation. However, when we simulated random mating among these individuals and in successive generations including their descendents, as occur in real populations, the correlations between markers and genome-wide heterozygosity were much higher than those predicted from algebraic analyses, and estimates of genome-wide heterozygosity improved slightly with the increment of the number of loci in the genome.     

Evidence of opposing fitness effects of parental heterozygosity and relatedness in a critically endangered marine turtle?

How individual genetic variability relates to fitness is important in understanding evolution and the processes affecting populations of conservation concern. Heterozygosity–fitness correlations (HFCs) have been widely used to study this link in wild populations, where key parameters that affect both variability and fitness, such as inbreeding, can be difficult to measure. We used estimates of parental heterozygosity and genetic similarity (‘relatedness’) derived from 32 microsatellite markers to explore the relationship between genetic variability and fitness in a population of the critically endangered hawksbill turtle, Eretmochelys imbricata. We found no effect of maternal MLH (multilocus heterozygosity) on clutch size or egg success rate, and no single‐locus effects. However, we found effects of paternal MLH and parental relatedness on egg success rate that interacted in a way that may result in both positive and negative effects of genetic variability. Multicollinearity in these tests was within safe limits, and null simulations suggested that the effect was not an artefact of using paternal genotypes reconstructed from large samples of offspring. Our results could imply a tension between inbreeding and outbreeding depression in this system, which is biologically feasible in turtles: female‐biased natal philopatry may elevate inbreeding risk and local adaptation, and both processes may be disrupted by male‐biased dispersal. Although this conclusion should be treated with caution due to a lack of significant identity disequilibrium, our study shows the importance of considering both positive and negative effects when assessing how variation in genetic variability affects fitness in wild systems.     

VNTR loci reveal differentiation between and structure within populations of the eastern barred bandicoot Perameles gunnii

The Eastern Barred Bandicoot Perameles gunnii has declined in abundance within mainland south-eastern Australia, to a relict wild population of less than 100 individuals in Hamilton, Victoria. It is more common, but is also declining in Tasmania. Genomic DN A variability was compared within and between surviving populations of P. gunnii using variable number of tandem repeat (VNTR) markers in one of two ways. First, average percentage differences (APDs) were determined between profiles for two VNTR probe—endonuclease combinations. Secondly, because one of these combinations revealed two multiallelic VNTR loci, genotypes were assigned and analysed for homogeneity of allele frequencies among subpopulations, for deviation of heterozygosity from Hardy-Weinberg equilibrium within populations and for genetic structuring among individuals from different subpopulations. The results of both the APD and defined locus approaches showed consistent trends within and between populations. Genetic variability was higher among mainland P. gunnii than in Tasmanian populations (higher APDs, number of alleles, and heterozygosity at one locus), despite the known decline and subdivision of the Hamilton population. Eleven per cent of the variability detected in Hamilton was attributed to genetic differentiation between east and west subdivisions of the population. Departure from random mating indicating local inbreeding within collecting localities was evident for one locus in both north and south Tasmania, particularly at one locality. AH alleles at both loci were unique to either Hamilton or Tasmanian P. gunnii. The initial captive colony contains high heterozygosity for these loci. It is concluded that VNTR markers can be of benefit for use in studies of population differentiation and for conservation management.     

Population differentiation and behavioural association of the two ‘personality’ genes DRD4 and SERT in dunnocks (Prunella modularis)

Quantifying the variation in behaviour‐related genes within and between populations provides insight into how evolutionary processes shape consistent behavioural traits (i.e. personality). Deliberate introductions of non‐native species offer opportunities to investigate how such genes differ between native and introduced populations and how polymorphisms in the genes are related to variation in behaviour. Here, we compared the genetic variation of the two ‘personality’ genes, DRD4 and SERT, between a native (United Kingdom, UK) and an introduced (New Zealand, NZ) population of dunnocks, Prunella modularis. The NZ population showed a significantly lower number of single nucleotide polymorphisms (SNPs) compared to the UK population. Standardized F’_st estimates of the personality genes and neutral microsatellites indicate that selection (anthropogenic and natural) probably occurred during and post the introduction event. Notably, the largest genetic differentiation was found in the intronic regions of the genes. In the NZ population, we also examined the association between polymorphisms in DRD4 and SERT and two highly repeatable behavioural traits: flight‐initiation distance and mating status (promiscuous females and cobreeding males). We found 38 significant associations (for different allele effect models) between the two behavioural traits and the studied genes. Further, 22 of the tested associations showed antagonistic allele effects for males and females. Our findings illustrate how introduction events and accompanying ecological changes could influence the genetic diversity of behaviour‐related genes.     

Life‐history evolution under fluctuating density‐dependent selection and the adaptive alignment of pace‐of‐life syndromes

We present a novel perspective on life‐history evolution that combines recent theoretical advances in fluctuating density‐dependent selection with the notion of pace‐of‐life syndromes (POLSs) in behavioural ecology. These ideas posit phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits as a continuum from the highly fecund, short‐lived, bold, aggressive and highly dispersive ‘fast’ types at one end of the POLS to the less fecund, long‐lived, cautious, shy, plastic and socially responsive ‘slow’ types at the other. We propose that such variation in life histories and the associated individual differences in behaviour can be explained through their eco‐evolutionary dynamics with population density – a single and ubiquitous selective factor that is present in all biological systems. Contrasting regimes of environmental stochasticity are expected to affect population density in time and space and create differing patterns of fluctuating density‐dependent selection, which generates variation in fast versus slow life histories within and among populations. We therefore predict that a major axis of phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits (i.e. the POLS) should align with these stochastic fluctuations in the multivariate fitness landscape created by variation in density‐dependent selection. Phenotypic plasticity and/or genetic (co‐)variation oriented along this major POLS axis are thus expected to facilitate rapid and adaptively integrated changes in various aspects of life histories within and among populations and/or species. The fluctuating density‐dependent selection POLS framework presented here therefore provides a series of clear testable predictions, the investigation of which should further our fundamental understanding of life‐history evolution and thus our ability to predict natural population dynamics.     

Effect of microsatellite selection on individual and population genetic inferences: an empirical study using cross‐specific and species‐specific amplifications

Although whole‐genome sequencing is becoming more accessible and feasible for nonmodel organisms, microsatellites have remained the markers of choice for various population and conservation genetic studies. However, the criteria for choosing microsatellites are still controversial due to ascertainment bias that may be introduced into the genetic inference. An empirical study of red deer (Cervus elaphus) populations, in which cross‐specific and species‐specific microsatellites developed through pyrosequencing of enriched libraries, was performed for this study. Two different strategies were used to select the species‐specific panels: randomly vs. highly polymorphic markers. The results suggest that reliable and accurate estimations of genetic diversity can be obtained using random microsatellites distributed throughout the genome. In addition, the results reinforce previous evidence that selecting the most polymorphic markers leads to an ascertainment bias in estimates of genetic diversity, when compared with randomly selected microsatellites. Analyses of population differentiation and clustering seem less influenced by the approach of microsatellite selection, whereas assigning individuals to populations might be affected by a random selection of a small number of microsatellites. Individual multilocus heterozygosity measures produced various discordant results, which in turn had impacts on the heterozygosity‐fitness correlation test. Finally, we argue that picking the appropriate microsatellite set should primarily take into account the ecological and evolutionary questions studied. Selecting the most polymorphic markers will generally overestimate genetic diversity parameters, leading to misinterpretations of the real genetic diversity, which is particularly important in managed and threatened populations.     

Environmentally induced dispersal‐related life‐history syndrome in the tropical butterfly,Bicyclus anynana

Dispersal is a key process for understanding the persistence of populations as well as the capacity of organisms to respond to environmental change. Therefore, understanding factors that may facilitate or constrain the evolution of dispersal is of crucial interest. Assessments of phenotypic variation in various behavioural, physiological and morphological traits related to insect dispersal and flight performance are common, yet very little is known about the genetic associations among these traits. We have used experiments on the butterfly Bicyclus anynana to estimate genetic variation and covariation in seven behavioural, physiological and morphological traits related to flight potential and hence dispersal. Our goal was to characterize the heritabilities and genetic correlations among these traits and thus to understand more about the evolution of dispersal‐related life‐history syndromes in butterflies. Using a version of the animal model, we showed that all of the traits varied between the sexes, and most were either positively or negatively (phenotypically and/or genetically) correlated with body size. Heritable variation was present in most traits, with the highest heritabilities estimated for body mass and thorax ratio. The variance in flight activity among multiple measurements for the same individual was high even after controlling for the prevailing environmental conditions, indicating the importance of behavioural switching and/or inherent randomness associated with this type of movement. A number of dispersal‐related traits showed phenotypic correlations among one another, but only a few of these were associated with significant genetic correlations indicating that covariances between these traits in Bicyclus anynana are mainly environmentally induced.     

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

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

Predation mediated population divergence in complex behaviour of nine‐spined stickleback (Pungitius pungitius)

The proximate and ultimate explanations for behavioural syndromes (correlated behaviours – a population trait) are poorly understood, and the evolution of behavioural types (configuration of behaviours – an individual trait) has been rarely studied. We investigated population divergence in behavioural syndromes and types using individually reared, completely predator‐ or conspecific‐naïve adult nine‐spined sticklebacks (Pungitius pungitius) from two marine and two predatory fish free, isolated pond populations. We found little evidence for the existence of behavioural syndromes, but population divergence in behavioural types was profound: individuals from ponds were quicker in feeding, bolder and more aggressive than individuals from marine environments. Our data reject the hypothesis that behavioural syndromes exist as a result of genetic correlations between behavioural traits, and support the contention that different behavioural types can be predominant in populations differing in predation pressure, most probably as a result of repeated independent evolution of separate behavioural traits.     

Genetic variation in the life‐history traits of Epiphyas postvittana: population structure and local adaptation

The light brown apple moth, Epiphyas postvittana (Walker) shows high intraspecific variability in morphological, physiological, demographic and behavioural characters. To gain insight into the extent of adaptation and evolutionary changes in response to environmental heterogeneity in this species, quantitative genetic analyses of life‐history variation were conducted for two natural populations under two thermal conditions (23°C and 28°C). Paternal half‐sib heritability and genetic correlation in six life‐history traits (i.e. development time, adult body weight, adult lifespan, age at first reproduction, the number of eggs laid during the first 5 days after emergence, and total fecundity) were compared. Significant heritabilities were shown consistently in development time; this is also true for adult body weight, except for the Canberra population at 23°C. However, neither population differences nor the effect of temperature were statistically detectable for any of these heritabilities, confirming the genetically determined flexibility. Positive genetic correlations between development time and adult body weight, and negative genetic correlations between the number of eggs laid during the first 5 days and adult lifespan were present for these populations at both temperatures, indicating the presence of genetic constraints. Pairwise comparisons of genetic correlations revealed the heterogeneity of the two populations and across temperatures. These results suggest that the structure of genetic covariance might have changed significantly during the divergence of natural populations and in response to the alteration of environmental conditions in E. postvittana.     

INBREEDING INFLUENCES WITHIN‐BROOD HETEROZYGOSITY‐FITNESS CORRELATIONS (HFCS) IN AN ISOLATED PASSERINE POPULATION

Molecular estimates of inbreeding may be made using genetic markers such as microsatellites, however the interpretation of resulting heterozygosity‐fitness correlations (HFCs) with respect to inbreeding depression is not straightforward. We investigated the relationship between pedigree‐determined inbreeding coefficients (f) and HFCs in a closely monitored, reintroduced population of Stewart Island robins (Petroica australis rakiura) on Ulva Island, New Zealand. Using a full sibling design, we focused on differences in juvenile survival associated specifically with individual sibling variation in standardized multilocus heterozygosity (SH) when expected f was identical. We found that within broods, siblings with higher SH at microsatellite loci experienced a higher probability of juvenile survival. This effect, however, was detected primarily within broods that experienced inbreeding or when inbreeding had occurred in their pedigree histories (i.e., at the parents’ level). Thus we show, for the first time in a wild population, that the strength of an HFC is partially dependent on the presence of inbreeding events in the recent pedigree history. Our results illustrate the importance of realized effects of inbreeding on genetic variation and fitness and the value of full‐sibling designs for the study of HFCs in the context of small, inbred populations.     

Automixis in Artemia: solving a century‐old controversy

Parthenogenesis (reproduction through unfertilized eggs) encompasses a variety of reproduction modes with (automixis) or without (apomixis) meiosis. Different modes of automixis have very different genetic and evolutionary consequences but can be particularly difficult to tease apart. In this study, we propose a new method to discriminate different types of automixis from population‐level genetic data. We apply this method to diploid Artemia parthenogenetica, a crustacean whose reproductive mode remains controversial despite a century of intensive cytogenetic observations. We focus on A. parthenogenetica from two western Mediterranean populations. We show that they are diploid and that markers remain heterozygous in cultures maintained up to ~36 generations in the laboratory. Moreover, parallel patterns of population‐wide heterozygosity levels between the two natural populations strongly support the conclusion that diploid A. parthenogenetica reproduce by automictic parthenogenesis with central fusion and low, but nonzero recombination. This settles a century‐old controversy on Artemia, and, more generally, suggests that many automictic organisms harbour steep within‐chromosome gradients of heterozygosity due to a transition from clonal transmission in centromere‐proximal regions to a form of inbreeding similar to self‐fertilization in centromere‐distal regions. Such systems therefore offer a new avenue for contrasting the genomic consequences of asexuality and inbreeding.     

Heterozygosity at neutral and immune loci is not associated with neonatal mortality due to microbial infection in Antarctic fur seals

Numerous studies have reported correlations between the heterozygosity of genetic markers and fitness. These heterozygosity–fitness correlations (HFCs) play a central role in evolutionary and conservation biology, yet their mechanistic basis remains open to debate. For example, fitness associations have been widely reported at both neutral and functional loci, yet few studies have directly compared the two, making it difficult to gauge the relative contributions of genome‐wide inbreeding and specific functional genes to fitness. Here, we compared the effects of neutral and immune gene heterozygosity on death from bacterial infection in Antarctic fur seal (Arctocephalus gazella) pups. We specifically developed a panel of 13 microsatellites from expressed immune genes and genotyped these together with 48 neutral loci in 234 individuals, comprising 39 pups that were classified at necropsy as having most likely died of bacterial infection together with a five times larger matched sample of healthy surviving pups. Identity disequilibrium quantified from the neutral markers was positive and significant, indicative of variance in inbreeding within the study population. However, multilocus heterozygosity did not differ significantly between healthy and infected pups at either class of marker, and little evidence was found for fitness associations at individual loci. These results support a previous study of Antarctic fur seals that found no effects of heterozygosity at nine neutral microsatellites on neonatal survival and thereby help to refine our understanding of how HFCs vary across the life cycle. Given that nonsignificant HFCs are underreported in the literature, we also hope that our study will contribute toward a more balanced understanding of the wider importance of this phenomenon.