Multi-generation genetic contributions of immigrants reveal cryptic elevated and sex-biased effective gene flow within a natural meta-population

Impacts of immigration on micro-evolution and population dynamics fundamentally depend on net rates and forms of resulting gene flow into recipient populations. Yet, the degrees to which observed rates and sex ratios of physical immigration translate into multi-generational genetic legacies have not been explicitly quantified in natural meta-populations, precluding inference on how movements translate into effective gene flow and eco-evolutionary outcomes. Our analyses of three decades of complete song sparrow (Melospiza melodia) pedigree data show that multi-generational genetic contributions from regular natural immigrants substantially exceeded those from contemporary natives, consistent with heterosis-enhanced introgression. However, while contributions from female immigrants exceeded those from female natives by up to three-fold, male immigrants' lineages typically went locally extinct soon after arriving. Both the overall magnitude, and the degree of female bias, of effective gene flow therefore greatly exceeded those which would be inferred from observed physical arrivals, altering multiple eco-evolutionary implications of immigration.

Candidate genes under selection in song sparrows co-vary with climate and body mass in support of Bergmann's Rule

Ecogeographic rules denote spatial patterns in phenotype and environment that may reflect local adaptation as well as a species' capacity to adapt to change. To identify genes underlying Bergmann's Rule, which posits that spatial correlations of body mass and temperature reflect natural selection and local adaptation in endotherms, we compare 79 genomes from nine song sparrow (Melospiza melodia) subspecies that vary ~300% in body mass (17 - 50 g). Comparing large- and smaller-bodied subspecies revealed 9 candidate genes in three genomic regions associated with body mass. Further comparisons to the five smallest subspecies endemic to California revealed eight SNPs within four of the candidate genes (GARNL3, RALGPS1, ANGPTL2, and COL15A1) associated with body mass and varying as predicted by Bergmann's Rule. Our results support the hypothesis that co-variation in environment, body mass and genotype reflect the influence of natural selection on local adaptation and a capacity for contemporary evolution in this diverse species.

Adaptation to climate change through seasonal migration revealed by climatic versus demographic niche models

Predicting the geographic range of species and their response to climatic variation and change are entwined goals in conservation and evolutionary ecology. Species distribution models (SDMs) are foundational in this effort and used to visualize the geographic range of species as the spatial representation of its realized niche. SDMs are also used to forecast range shifts under climate change, but often in the absence of empirical evidence that climate limits population growth. We explored the influence of climate on demography, seasonal migration, and the extent of the geographic range in song sparrows (Melospiza melodia), a species thought to display marked local adaptation to regional climate. To do so, we developed SDMs to predict the demographic and climate niches of migratory and resident song sparrows across our study area in western North America from California to Alaska, using 48 years of demographic data from a focal population in British Columbia and 1.2 million continental-scale citizen science observations. Spatial agreement of our demographic and climate niche models in the region of our focal population was strong (76%), supporting the hypothesis that demographic performance and the occurrence of seasonal migration varied predictably with climatic conditions. In contrast, agreement at the northern (58%) and southern (40%) extents of our study area was lower, as expected if the factors limiting population growth vary regionally. Our results support the hypothesis that local climate drives spatial variation in the occurrence of seasonal migration in song sparrows by limiting the fitness of year-round residents, and suggest that climate warming has favored range expansions and facilitated an upward shift in elevational range song sparrows that forgo seasonal migration. Our work highlights the potential role of seasonal migration in climate adaptation and limits on the reliability of climate niche models not validated with demographic data.

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.

Understanding the Social Dynamics of Breeding Phenology: Indirect Genetic Effects and Assortative Mating in a Long-Distance Migrant

AbstractPhenological traits, such as the timing of reproduction, are often influenced by social interactions between paired individuals. Such partner effects may occur when pair members affect each other's prebreeding environment. Partner effects can be environmentally and/or genetically determined, and quantifying direct and indirect genetic effects is important for understanding the evolutionary dynamics of phenological traits. Here, using 26 years of data from a pedigreed population of a migratory seabird, the common tern (Sterna hirundo), we investigate male and female effects on female laying date. We find that female laying date harbors both genetic and environmental variation and is additionally influenced by the environmental and, to a lesser extent, genetic component of its mate. We demonstrate this partner effect to be largely explained by male arrival date. Interestingly, analyses of mating patterns with respect to arrival date show mating to be strongly assortative, and using simulations we show that assortative mating leads to overestimation of partner effects. Our study provides evidence for partner effects on breeding phenology in a long-distance migrant while uncovering the potential causal pathways underlying the observed effects and raising awareness for confounding effects resulting from assortative mating or other common environmental effects.

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

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

Heritabilities, social environment effects and genetic correlations of social behaviours in a cooperatively breeding vertebrate

Social animals interact frequently with conspecifics, and their behaviour is influenced by social context, environmental cues and the behaviours of interaction partners, allowing for adaptive, flexible adjustments to social encounters. This flexibility can be limited by part of the behavioural variation being genetically determined. Furthermore, behaviours can be genetically correlated, potentially constraining independent evolution. Understanding social behaviour thus requires carefully disentangling genetic, environmental, maternal and social sources of variations as well as the correlation structure between behaviours. Here, we assessed heritability, maternal, common environment and social effects of eight social behaviours in Neolamprologus pulcher, a cooperatively breeding cichlid. We bred wild-caught fish in a paternal half-sibling design and scored ability to defend a resource against conspecifics, to integrate into a group and the propensity to help defending the group territory ("helping behaviour"). We assessed genetic, social and phenotypic correlations within clusters of behaviours predicted to be functionally related, namely "competition," "aggression," "aggression-sociability," "integration" and "integration-help." Helping behaviour and two affiliative behaviours were heritable, whereas there was little evidence for a genetic basis in all other traits. Phenotypic social effects explained part of the variation in a sociable and a submissive behaviour, but there were no maternal or common environment effects. Genetic and phenotypic correlation within clusters was mostly positive. A group's social environment influenced covariances of social behaviours. Genetic correlations were similar in magnitude but usually exceeding the phenotypic ones, indicating that conclusions about the evolution of social behaviours in this species could be provisionally drawn from phenotypic data in cases where data for genetic analyses are unobtainable.

Animal models with group-specific additive genetic variances: extending genetic group models

BACKGROUND

The animal model is a key tool in quantitative genetics and has been used extensively to estimate fundamental parameters, such as additive genetic variance or heritability. An implicit assumption of animal models is that all founder individuals derive from a single population. This assumption is commonly violated, for instance in crossbred livestock or when a meta-population is split into genetically differentiated subpopulations. Ignoring that base populations are genetically heterogeneous and thus split into different 'genetic groups' may lead to biased parameter estimates, especially for additive genetic variance. To avoid such biases, genetic group animal models, which account for the presence of more than one genetic group, have been proposed. Unfortunately, the method to date is only computationally feasible when the breeding values of the groups are allowed to differ in their means, but not in their variances.

RESULTS

We present an extension of the animal model that permits estimation of group-specific additive genetic variances. This is achieved by employing group-specific relatedness matrices for the breeding value components to different genetic groups. We derive these matrices by decomposing the full relatedness matrix via the generalized Cholesky decomposition, and by scaling the respective matrix components for each group. We propose a computationally convenient approximation for the matrix component that encodes for the Mendelian sampling variance, and show that this approximation is not critical. In addition, we explain why segregation variances are often negligible when analyzing the complex polygenic traits that are frequently the focus of evolutionary ecologists and animal breeders. Simulations and an example from an insular meta-population of house sparrows in Norway with three distinct genetic groups illustrate that the method is successful in estimating group-specific additive genetic variances, and that segregation variances are indeed negligible in the empirical example.

CONCLUSIONS

Quantifying differences in additive genetic variance within and among populations is of major biological interest in ecology, evolution, and animal and plant breeding. The proposed method allows to estimate such differences for subpopulations that form a connected set of populations, and may thus also be useful to study temporal or spatial variation of additive genetic variances.

The Contemporary Evolution of Fitness

The rate of evolution of population mean fitness informs how selection acting in contemporary populations can counteract environmental change and genetic degradation (mutation, gene flow, drift, recombination). This rate influences population increases (e.g., range expansion), population stability (e.g., cryptic eco-evolutionary dynamics), and population recovery (i.e., evolutionary rescue). We review approaches for estimating such rates, especially in wild populations. We then review empirical estimates derived from two approaches: mutation accumulation (MA) and additive genetic variance in fitness (IAw). MA studies inform how selection counters genetic degradation arising from deleterious mutations, typically generating estimates of <1% per generation. IAw studies provide an integrated prediction of proportional change per generation, nearly always generating estimates of <20% and, more typically, <10%. Overall, considerable, but not unlimited, evolutionary potential exists in populations facing detrimental environmental or genetic change. However, further studies with diverse methods and species are required for more robust and general insights.

Individual repeatability and heritability of divorce in a wild population

Understanding micro-evolutionary responses of mating systems to contemporary selection requires estimating sex-specific additive genetic variances and cross-sex genetic covariances in key reproductive strategy traits. One key trait comprises the occurrence of divorce versus mate fidelity across sequential reproductive attempts. If divorce represents an evolving behavioural strategy that responds to selection it must have non-zero individual repeatability and heritability, but quantitative estimates from wild populations are scarce. We used 39 years of individual breeding records and pedigree data from free-living song sparrows (Melospiza melodia) to quantify sex-specific permanent individual and additive genetic variances, and hence estimate repeatability and heritability, in liability for divorce. We estimated moderate repeatability among females, but little repeatability among males. Estimates of additive genetic variance were small in both sexes, and the cross-sex genetic covariance was close to zero. Consequently, the total heritability was small but likely non-zero, indicating low potential for micro-evolutionary response to selection. Rapid micro-evolutionary change of divorce rate, therefore, appears unlikely, even if there were substantial fitness benefits of divorce and resulting selection.

The Consequences of Polyandry for Sibship Structures, Distributions of Relationships and Relatedness, and Potential for Inbreeding in a Wild Population

The evolutionary benefits of simultaneous polyandry (female multiple mating within a single reproductive event) remain elusive. One potential benefit could arise if polyandry alters sibship structures and consequent relationships and relatedness among females' descendants, thereby intrinsically reducing future inbreeding risk (the indirect inbreeding avoidance hypothesis). However such effects have not been quantified in naturally complex mating systems that also encompass iteroparity, overlapping generations, sequential polyandry, and polygyny. We used long-term social and genetic pedigree data from song sparrows (Melospiza melodia) to quantify cross-generational consequences of simultaneous polyandry for offspring sibship structures and distributions of relationships and relatedness among possible mates. Simultaneous polyandry decreased full sibships and increased half-sibships, on average, but such effects varied among females and were smaller than would occur in the absence of sequential polyandry or polygyny. Further, while simultaneous polyandry decreased the overall frequencies of possible matings among adult full sibs, it increased the frequencies of possible matings among adult half-sibs and more distant relatives. These results imply that the intrinsic consequences of simultaneous polyandry for inbreeding risk could cause weak indirect selection on polyandry, but the magnitude and direction of such effects will depend on complex interactions with other mating system components and the form of inbreeding depression.

No evidence of inbreeding depression in sperm performance traits in wild song sparrows

Inbreeding is widely hypothesized to shape mating systems and population persistence, but such effects will depend on which traits show inbreeding depression. Population and evolutionary consequences could be substantial if inbreeding decreases sperm performance and hence decreases male fertilization success and female fertility. However, the magnitude of inbreeding depression in sperm performance traits has rarely been estimated in wild populations experiencing natural variation in inbreeding. Further, the hypothesis that inbreeding could increase within-ejaculate variation in sperm traits and thereby further affect male fertilization success has not been explicitly tested. We used a wild pedigreed song sparrow (Melospiza melodia) population, where frequent extrapair copulations likely create strong postcopulatory competition for fertilization success, to quantify effects of male coefficient of inbreeding (f) on key sperm performance traits. We found no evidence of inbreeding depression in sperm motility, longevity, or velocity, and the within-ejaculate variance in sperm velocity did not increase with male f. Contrary to inferences from highly inbred captive and experimental populations, our results imply that moderate inbreeding will not necessarily constrain sperm performance in wild populations. Consequently, the widely observed individual-level and population-level inbreeding depression in male and female fitness may not stem from reduced sperm performance in inbred males.

Accounting for Sampling Error in Genetic Eigenvalues Using Random Matrix Theory

The distribution of genetic variance in multivariate phenotypes is characterized by the empirical spectral distribution of the eigenvalues of the genetic covariance matrix. Empirical estimates of genetic eigenvalues from random effects linear models are known to be overdispersed by sampling error, where large eigenvalues are biased upward, and small eigenvalues are biased downward. The overdispersion of the leading eigenvalues of sample covariance matrices have been demonstrated to conform to the Tracy-Widom (TW) distribution. Here we show that genetic eigenvalues estimated using restricted maximum likelihood (REML) in a multivariate random effects model with an unconstrained genetic covariance structure will also conform to the TW distribution after empirical scaling and centering. However, where estimation procedures using either REML or MCMC impose boundary constraints, the resulting genetic eigenvalues tend not be TW distributed. We show how using confidence intervals from sampling distributions of genetic eigenvalues without reference to the TW distribution is insufficient protection against mistaking sampling error as genetic variance, particularly when eigenvalues are small. By scaling such sampling distributions to the appropriate TW distribution, the critical value of the TW statistic can be used to determine if the magnitude of a genetic eigenvalue exceeds the sampling error for each eigenvalue in the spectral distribution of a given genetic covariance matrix.

Accounting for genetic differences among unknown parents in microevolutionary studies: how to include genetic groups in quantitative genetic animal models

Quantifying and predicting microevolutionary responses to environmental change requires unbiased estimation of quantitative genetic parameters in wild populations. 'Animal models', which utilize pedigree data to separate genetic and environmental effects on phenotypes, provide powerful means to estimate key parameters and have revolutionized quantitative genetic analyses of wild populations. However, pedigrees collected in wild populations commonly contain many individuals with unknown parents. When unknown parents are non-randomly associated with genetic values for focal traits, animal model parameter estimates can be severely biased. Yet, such bias has not previously been highlighted and statistical methods designed to minimize such biases have not been implemented in evolutionary ecology. We first illustrate how the occurrence of non-random unknown parents in population pedigrees can substantially bias animal model predictions of breeding values and estimates of additive genetic variance, and create spurious temporal trends in predicted breeding values in the absence of local selection. We then introduce 'genetic group' methods, which were developed in agricultural science, and explain how these methods can minimize bias in quantitative genetic parameter estimates stemming from genetic heterogeneity among individuals with unknown parents. We summarize the conceptual foundations of genetic group animal models and provide extensive, step-by-step tutorials that demonstrate how to fit such models in a variety of software programs. Furthermore, we provide new functions in r that extend current software capabilities and provide a standardized approach across software programs to implement genetic group methods. Beyond simply alleviating bias, genetic group animal models can directly estimate new parameters pertaining to key biological processes. We discuss one such example, where genetic group methods potentially allow the microevolutionary consequences of local selection to be distinguished from effects of immigration and resulting gene flow. We highlight some remaining limitations of genetic group models and discuss opportunities for further development and application in evolutionary ecology. We suggest that genetic group methods should no longer be overlooked by evolutionary ecologists, but should become standard components of the toolkit for animal model analyses of wild population data sets.

Direct and indirect genetic and fine-scale location effects on breeding date in song sparrows

Quantifying direct and indirect genetic effects of interacting females and males on variation in jointly expressed life-history traits is central to predicting microevolutionary dynamics. However, accurately estimating sex-specific additive genetic variances in such traits remains difficult in wild populations, especially if related individuals inhabit similar fine-scale environments. Breeding date is a key life-history trait that responds to environmental phenology and mediates individual and population responses to environmental change. However, no studies have estimated female (direct) and male (indirect) additive genetic and inbreeding effects on breeding date, and estimated the cross-sex genetic correlation, while simultaneously accounting for fine-scale environmental effects of breeding locations, impeding prediction of microevolutionary dynamics. We fitted animal models to 38 years of song sparrow (Melospiza melodia) phenology and pedigree data to estimate sex-specific additive genetic variances in breeding date, and the cross-sex genetic correlation, thereby estimating the total additive genetic variance while simultaneously estimating sex-specific inbreeding depression. We further fitted three forms of spatial animal model to explicitly estimate variance in breeding date attributable to breeding location, overlap among breeding locations and spatial autocorrelation. We thereby quantified fine-scale location variances in breeding date and quantified the degree to which estimating such variances affected the estimated additive genetic variances. The non-spatial animal model estimated nonzero female and male additive genetic variances in breeding date (sex-specific heritabilities: 0·07 and 0·02, respectively) and a strong, positive cross-sex genetic correlation (0·99), creating substantial total additive genetic variance (0·18). Breeding date varied with female, but not male inbreeding coefficient, revealing direct, but not indirect, inbreeding depression. All three spatial animal models estimated small location variance in breeding date, but because relatedness and breeding location were virtually uncorrelated, modelling location variance did not alter the estimated additive genetic variances. Our results show that sex-specific additive genetic effects on breeding date can be strongly positively correlated, which would affect any predicted rates of microevolutionary change in response to sexually antagonistic or congruent selection. Further, we show that inbreeding effects on breeding date can also be sex specific and that genetic effects can exceed phenotypic variation stemming from fine-scale location-based variation within a wild population.

Variation in parent-offspring kinship in socially monogamous systems with extra-pair reproduction and inbreeding

Female extra‐pair reproduction in socially monogamous systems is predicted to cause cuckolded socially‐paired males to conditionally reduce paternal care, causing selection against extra‐pair reproduction and underlying polyandry. However, existing models and empirical studies have not explicitly considered that cuckolded males might be related to their socially‐paired female and/or to her extra‐pair mate, and therefore be related to extra‐pair offspring that they did not sire but could rear. Selection against paternal care, and hence against extra‐pair reproduction, might then be weakened. We derive metrics that quantify allele‐sharing between within‐pair and extra‐pair offspring and their mother and her socially‐paired male in terms of coefficients of kinship and inbreeding. We use song sparrow (Melospiza melodia) paternity and pedigree data to quantify these metrics, and thereby quantify the joint effects of extra‐pair reproduction and inbreeding on a brood's total allelic value to its socially‐paired parents. Cuckolded male song sparrows were almost always detectably related to extra‐pair offspring they reared. Consequently, although brood allelic value decreased substantially following female extra‐pair reproduction, this decrease was reduced by within‐pair and extra‐pair reproduction among relatives. Such complex variation in kinship within nuclear families should be incorporated into models considering coevolutionary dynamics of extra‐pair reproduction, parental care, and inbreeding.