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.

Multigenerational Fitness Effects of Natural Immigration Indicate Strong Heterosis and Epistatic Breakdown in a Wild Bird Population.

AbstractThe fitness of immigrants and their descendants produced within recipient populations fundamentally underpins the genetic and population dynamic consequences of immigration. Immigrants can in principle induce contrasting genetic effects on fitness across generations, reflecting multifaceted additive, dominance, and epistatic effects. Yet full multigenerational and sex-specific fitness effects of regular immigration have not been quantified within naturally structured systems, precluding inference on underlying genetic architectures and population outcomes. We used four decades of song sparrow (Melospiza melodia) life history and pedigree data to quantify fitness of natural immigrants, natives, and their F1, F2, and backcross descendants and test for evidence of nonadditive genetic effects. Values of key fitness components (including adult lifetime reproductive success and zygote survival) of F1 offspring of immigrant-native matings substantially exceeded their parent mean, indicating strong heterosis. Meanwhile, F2 offspring of F1-F1 matings had notably low values, indicating surprisingly strong epistatic breakdown. Furthermore, magnitudes of effects varied among fitness components and differed between female and male descendants. These results demonstrate that strong nonadditive genetic effects on fitness can arise within weakly structured and fragmented populations experiencing frequent natural immigration. Such effects will substantially affect the net degree of effective gene flow and resulting local genetic introgression and adaptation.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation.

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

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.

Pedigree-based inbreeding coefficient explains more variation in fitness than heterozygosity at 160 microsatellites in a wild bird population.

Although the pedigree-based inbreeding coefficient F predicts the expected proportion of an individual's genome that is identical-by-descent (IBD), heterozygosity at genetic markers captures Mendelian sampling variation and thereby provides an estimate of realized IBD. Realized IBD should hence explain more variation in fitness than their pedigree-based expectations, but how many markers are required to achieve this in practice remains poorly understood. We use extensive pedigree and life-history data from an island population of song sparrows (Melospiza melodia) to show that the number of genetic markers and pedigree depth affected the explanatory power of heterozygosity and F, respectively, but that heterozygosity measured at 160 microsatellites did not explain more variation in fitness than F This is in contrast with other studies that found heterozygosity based on far fewer markers to explain more variation in fitness than F Thus, the relative performance of marker- and pedigree-based estimates of IBD depends on the quality of the pedigree, the number, variability and location of the markers employed, and the species-specific recombination landscape, and expectations based on detailed and deep pedigrees remain valuable until we can routinely afford genotyping hundreds of phenotyped wild individuals of genetic non-model species for thousands of genetic markers.

Balancing selection and genetic drift create unusual patterns of MHCIIß variation in Galápagos mockingbirds.

The extracellular subunit of the major histocompatibility complex MHCIIß plays an important role in the recognition of pathogens and the initiation of the adaptive immune response of vertebrates. It is widely accepted that pathogen-mediated selection in combination with neutral micro-evolutionary forces (e.g. genetic drift) shape the diversity of MHCIIß, but it has proved difficult to determine the relative effects of these forces. We evaluated the effect of genetic drift and balancing selection on MHCIIß diversity in 12 small populations of Galápagos mockingbirds belonging to four different species, and one larger population of the Northern mockingbird from the continental USA. After genotyping MHCIIß loci by high-throughput sequencing, we applied a correlational approach to explore the relationships between MHCIIß diversity and population size by proxy of island size. As expected when drift predominates, we found a positive effect of population size on the number of MHCIIß alleles present in a population. However, the number of MHCIIß alleles per individual and number of supertypes were not correlated with population size. This discrepancy points to an interesting feature of MHCIIß diversity dynamics: some levels of diversity might be shaped by genetic drift while others are independent and possibly maintained by balancing selection.

Mother-offspring and nest-mate resemblance but no heritability in early-life telomere length in white-throated dippers.

Telomeres are protective DNA-protein complexes located at the ends of eukaryotic chromosomes, whose length has been shown to predict life-history parameters in various species. Although this suggests that telomere length is subject to natural selection, its evolutionary dynamics crucially depends on its heritability. Using pedigree data for a population of white-throated dippers (Cinclus cinclus), we test whether and how variation in early-life relative telomere length (RTL, measured as the number of telomeric repeats relative to a control gene using qPCR) is transmitted across generations. We disentangle the relative effects of genes and environment and test for sex-specific patterns of inheritance. There was strong and significant resemblance among offspring sharing the same nest and offspring of the same cohort. Furthermore, although offspring resemble their mother, and there is some indication for an effect of inbreeding, additive genetic variance and heritability are close to zero. We find no evidence for a role of either maternal imprinting or Z-linked inheritance in generating these patterns, suggesting they are due to non-genetic maternal and common environment effects instead. We conclude that in this wild bird population, environmental factors are the main drivers of variation in early-life RTL, which will severely bias estimates of heritability when not modelled explicitly.

Reverse attenuation in interaction terms due to covariate measurement error.

Covariate measurement error may cause biases in parameters of regression coefficients in generalized linear models. The influence of measurement error on interaction parameters has, however, only rarely been investigated in depth, and if so, attenuation effects were reported. In this paper, we show that also reverse attenuation of interaction effects may emerge, namely when heteroscedastic measurement error or sampling variances of a mismeasured covariate are present, which are not unrealistic scenarios in practice. Theoretical findings are illustrated with simulations. A Bayesian approach employing integrated nested Laplace approximations is suggested to model the heteroscedastic measurement error and covariate variances, and an application shows that the method is able to reveal approximately correct parameter estimates.

Hybrid ancestry of an island subspecies of Galápagos mockingbird explains discordant gene trees.

Introgression of genes through hybridization has been proposed to be an important driver of speciation, but in animals this has been shown only in relatively few cases until recently. Additionally, introgressive hybridization among non-sister species leads to a change in the gene tree topology of the concerned loci and thus complicates phylogenetic reconstruction. However, such cases of ancient introgression have been very difficult to demonstrate in birds. Here, we present such an example in an island bird subspecies, the Genovesa mockingbird (Mimus parvulus bauri). We assessed phylogenetic relationships and population structure among mockingbirds of the Galápagos archipelago using mitochondrial and nuclear DNA sequences, autosomal microsatellites, and morphological measurements. Mitochondrial haplotypes of Genovesa mockingbirds clustered closely with the haplotypes from two different species, San Cristóbal (M. melanotis) and Española (M. macdonaldi) mockingbirds. The same pattern was found for some haplotypes of two nuclear gene introns, while the majority of nuclear haplotypes of Genovesa mockingbirds were shared with other populations of the same species (M. parvulus). At 26 autosomal microsatellites, Genovesa mockingbirds grouped with other M. parvulus populations. This pattern shows that Genovesa mockingbirds contain mitochondria and some autosomal alleles that have most likely introgressed from M. melanotis into a largely M. parvulus background, making Genovesa mockingbirds a lineage of mixed ancestry, possibly undergoing speciation. Consistent with this hypothesis, mockingbirds on Genovesa are more clearly differentiated morphologically from other M. parvulus populations than M. melanotis is from M. parvulus.

Are there indirect fitness benefits of female extra-pair reproduction? Lifetime reproductive success of within-pair and extra-pair offspring.

The forces driving extra-pair reproduction by socially monogamous females, and the resulting genetic polyandry, remain unclear. A testable prediction of the hypothesis that extra-pair reproduction partly reflects indirect selection on females is that extra-pair young (EPY) will be fitter than their within-pair young (WPY) maternal half-siblings. This prediction has not been comprehensively tested in a wild population, requiring data on the lifetime reproductive success (LRS) of maternal half-sib EPY and WPY. We used 17 years of genetic parentage data from song sparrows, Melospiza melodia, to compare the LRS of hatched EPY and WPY maternal half-siblings measured as their lifetime number of hatched offspring, recruited offspring, and hatched grandoffspring. EPY hatchlings were not significantly fitter than WPY hatchlings for any of three measures of LRS. Furthermore, opposite to prediction, EPY hatchlings tended to have lower LRS than their maternal half-sibling WPY hatchlings on average. EPY also tended to be less likely to survive to hatch than their maternal half-sibling WPY. Taken together, these results fail to support one key hypothesis explaining the evolution of genetic polyandry by socially monogamous females and suggest there may be weak indirect selection against female extra-pair reproduction in song sparrows.

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.

A hitchhikers guide to the Galápagos: co-phylogeography of Galápagos mockingbirds and their parasites.

BACKGROUND: Parasites are evolutionary hitchhikers whose phylogenies often track the evolutionary history of their hosts. Incongruence in the evolutionary history of closely associated lineages can be explained through a variety of possible events including host switching and host independent speciation. However, in recently diverged lineages stochastic population processes, such as retention of ancestral polymorphism or secondary contact, can also explain discordant genealogies, even in fully co-speciating taxa. The relatively simple biogeographic arrangement of the Galápagos archipelago, compared with mainland biomes, provides a framework to identify stochastic and evolutionary informative components of genealogic data in these recently diverged organisms. RESULTS: Mitochondrial DNA sequences were obtained for four species of Galápagos mockingbirds and three sympatric species of ectoparasites--two louse and one mite species. These data were complemented with nuclear EF1α sequences in selected samples of parasites and with information from microsatellite loci in the mockingbirds. Mitochondrial sequence data revealed differences in population genetic diversity between all taxa and varying degrees of topological congruence between host and parasite lineages. A very low level of genetic variability and lack of congruence was found in one of the louse parasites, which was excluded from subsequent joint analysis of mitochondrial data. The reconciled multi-species tree obtained from the analysis is congruent with both the nuclear data and the geological history of the islands. CONCLUSIONS: The gene genealogies of Galápagos mockingbirds and two of their ectoparasites show strong phylogeographic correlations, with instances of incongruence mostly explained by ancestral genetic polymorphism. A third parasite genealogy shows low levels of genetic diversity and little evidence of co-phylogeny with their hosts. These differences can mostly be explained by variation in life-history characteristics, primarily host specificity and dispersal capabilities. We show that pooling genetic data from organisms living in close ecological association reveals a more accurate phylogeographic history for these taxa. Our results have implications for the conservation and taxonomy of Galápagos mockingbirds and their parasites.

Additive genetic variance, heritability, and inbreeding depression in male extra-pair reproductive success.

The hypothesis that female extra-pair reproduction in socially monogamous animals reflects indirect genetic benefits requires that there be additive and/or nonadditive genetic variance in fitness. However, the specific hypotheses that male extra-pair reproductive success (EPRS) shows additive genetic variance (V(A)), heritability (h2), or inbreeding depression, and hence that females could acquire indirect genetic benefits through increased EPRS of sons, have not been explicitly tested. We used comprehensive genetic pedigree data from song sparrows (Melospiza melodia) to estimate V(A), h2, and inbreeding depression in the number of extra-pair offspring a male sired per year and the probability that a male would sire any extra-pair offspring per year. Inbreeding depression was substantial: more inbred males sired fewer extra-pair offspring and were less likely to sire any extra-pair offspring. In contrast, estimates of V(A) and h2 were close to 0, although 95% credible intervals were relatively wide. These data suggest that females could accrue indirect genetic benefits, in terms of increased EPRS of outbred sons, by mating with unrelated social or extra-pair mates. In contrast, any indirect benefit of extra-pair reproduction in terms of producing sons with high additive genetic value for EPRS is most likely to be small.

Heritability of female extra-pair paternity rate in song sparrows (Melospiza melodia).

The forces driving the evolution of extra-pair reproduction in socially monogamous animals remain widely debated and unresolved. One key hypothesis is that female extra-pair reproduction evolves through indirect genetic benefits, reflecting increased additive genetic value of extra-pair offspring. Such evolution requires that a female's propensity to produce offspring that are sired by an extra-pair male is heritable. However, additive genetic variance and heritability in female extra-pair paternity (EPP) rate have not been quantified, precluding accurate estimation of the force of indirect selection. Sixteen years of comprehensive paternity and pedigree data from socially monogamous but genetically polygynandrous song sparrows (Melospiza melodia) showed significant additive genetic variance and heritability in the proportion of a female's offspring that was sired by an extra-pair male, constituting major components of the genetic architecture required for extra-pair reproduction to evolve through indirect additive genetic benefits. However, estimated heritabilities were moderately small (0.12 and 0.18 on the observed and underlying latent scales, respectively). The force of selection on extra-pair reproduction through indirect additive genetic benefits may consequently be relatively weak. However, the additive genetic variance and non-zero heritability observed in female EPP rate allow for multiple further genetic mechanisms to drive and constrain mating system evolution.

Sex-specific differential survival of extra-pair and within-pair offspring in song sparrows, Melospiza melodia.

It is widely hypothesized that the evolution of female extra-pair reproduction in socially monogamous species reflects indirect genetic benefits to females. However, a critical prediction of this hypothesis, that extra-pair young (EPY) are fitter than within-pair young (WPY), has rarely been rigorously tested. We used 18 years of data from free-living song sparrows, Melospiza melodia, to test whether survival through major life-history stages differed between EPY and WPY maternal half-siblings. On average, survival of hatched chicks to independence from parental care and recruitment, and their total lifespan, did not differ significantly between EPY and WPY. However, EPY consistently tended to be less likely to survive, and recruited EPY survived for significantly fewer years than recruited WPY. Furthermore, the survival difference between EPY and WPY was sex-specific; female EPY were less likely to survive to independence and recruitment and lived fewer years than female WPY, whereas male EPY were similarly or slightly more likely to survive and to live more years than male WPY. These data indicate that extra-pair paternity may impose an indirect cost on females via their female offspring and that sex-specific genetic, environmental or maternal effects may shape extra-pair reproduction.

Disentangling the effect of genes, the environment and chance on sex ratio variation in a wild bird population.

Sex ratio theory proposes that the equal sex ratio typically observed in birds and mammals is the result of natural selection. However, in species with chromosomal sex determination, the same 1 : 1 sex ratio is expected under random Mendelian segregation. Here, we present an analysis of 14 years of sex ratio data for a population of song sparrows (Melospiza melodia) on Mandarte Island, at the nestling stage and at independence from parental care. We test for the presence of variance in sex ratio over and above the binomial variance expected under Mendelian segregation, and thereby quantify the potential for selection to shape sex ratio. Furthermore, if sex ratio variation is to be shaped by selection, we expect some of this extra-binomial variation to have a genetic basis. Despite ample statistical power, we find no evidence for the existence of either genetic or environmentally induced variation in sex ratio, in the nest or at independence. Instead, the sex ratio variation observed matches that expected under random Mendelian segregation. Using one of the best datasets of its kind, we conclude that female song sparrows do not, and perhaps cannot, adjust the sex of their offspring. We discuss the implications of this finding and make suggestions for future research.

Immigration counter-acts local micro-evolution of a major fitness component: Migration-selection balance in free-living song sparrows.

Ongoing adaptive evolution, and resulting "evolutionary rescue" of declining populations, requires additive genetic variation in fitness. Such variation can be increased by gene flow resulting from immigration, potentially facilitating evolution. But, gene flow could in fact constrain rather than facilitate local adaptive evolution if immigrants have low additive genetic values for local fitness. Local migration-selection balance and micro-evolutionary stasis could then result. However, key quantitative genetic effects of natural immigration, comprising the degrees to which gene flow increases the total local additive genetic variance yet counteracts local adaptive evolutionary change, have not been explicitly quantified in wild populations. Key implications of gene flow for population and evolutionary dynamics consequently remain unclear. Our quantitative genetic analyses of long-term data from free-living song sparrows (Melospiza melodia) show that mean breeding value for local juvenile survival to adulthood, a major component of fitness, increased across cohorts more than expected solely due to drift. Such micro-evolutionary change should be expected given nonzero additive genetic variance and consistent directional selection. However, this evolutionary increase was counteracted by negative additive genetic effects of recent immigrants, which increased total additive genetic variance but prevented a net directional evolutionary increase in total additive genetic value. These analyses imply an approximate quantitative genetic migration-selection balance in a major fitness component, and hence demonstrate a key mechanism by which substantial additive genetic variation can be maintained yet decoupled from local adaptive evolutionary change.

Comprehensive paternity assignment: genotype, spatial location and social status in song sparrows, Melospiza Melodia.

Comprehensive, accurate paternity assignment is critical to answering numerous questions in evolutionary ecology. Yet, most studies of species with extra-pair paternity (EPP) fail to assign sires to all offspring. Common limitations include incomplete and biased sampling of offspring and males, particularly with respect to male location and social status, potentially biasing estimated patterns of paternity. Studies that achieve comprehensive sampling and paternity assignment are therefore required. Accordingly, we genotyped virtually all males and >99% of 6-day-old offspring over 16 years in a song sparrow (Melospiza melodia) population and used three complementary statistical methodologies to attempt complete paternity assignment for all 2207 offspring. Assignments were highly consistent across maximum likelihood methods that used solely genotype data, and heuristic and integrated Bayesian analyses that included data describing individual locations. Sires were assigned to >99% of all genotyped offspring with ≥95% confidence, revealing an EPP rate of c. 28%. Extra-pair sires primarily occupied territories neighbouring their extra-pair offspring; spatial location was therefore highly informative for paternity assignment. EPP was biased towards paired territorial males, although unpaired territorial and floater males sired c. 13% of extra-pair offspring. Failing to sample and include unpaired males as candidate sires would therefore substantially reduce assignment rates. These analyses demonstrate the integration of genetic and ecological information to achieve comprehensive paternity assignment and direct biological insight, illustrate the potential biases that common forms of incomplete sampling could have on estimated patterns of EPP, and provide an essential basis for understanding the evolutionary causes and consequences of EPP.

Inbreeding coefficient and heterozygosity-fitness correlations in unhatched and hatched song sparrow nestmates.

Heterozygosity-fitness correlations use molecular measures of heterozygosity as proxy estimates of individual inbreeding coefficients (f) to examine relationships between inbreeding and fitness traits. Heterozygosity-fitness correlations partly depend on the assumption that individual heterozygosity and f are strongly and negatively correlated. Although theory predicts that this relationship will be strongest when mean f and variance in f are high, few studies of heterozygosity-fitness correlations include estimates of f based on pedigrees, which allow for more thorough examinations of the relationship between f, heterozygosity and fitness in nature. We examined relationships between pedigree-based estimates of f, multilocus heterozygosity (MLH) and the probability of survival to hatch in song sparrow nestmates. f and MLH were weakly, but significantly negatively correlated. Inbreeding coefficient predicted the probability of survival to hatch. In contrast, MLH did not predict the probability of survival to hatch nor did it account for residual variation in survival to hatch after statistically controlling for the effects of f. These results are consistent with the expectation that heterozygosity-f correlations will be weak when mean and variance in f are low. Our results also provide empirical support for recent simulation studies, which show that variation in MLH among siblings with equal f can be large and may obscure MLH-fitness relationships.

Saving Darwin's muse: evolutionary genetics for the recovery of the Floreana mockingbird.

The distribution of mockingbird species among the Galápagos Islands prompted Charles Darwin to question, for the first time in writing, the 'stability of species'. Some 50 years after Darwin's visit, however, the endemic Floreana mockingbird (Mimus trifasciatus) had become extinct on Floreana Island and, today, only two small populations survive on two satellite islets. As Darwin noted, rarity often precedes extinction. To avert extinction, plans are being developed to reintroduce M. trifasciatus to Floreana. Here, we integrate evolutionary thinking and conservation practice using coalescent analyses and genetic data from contemporary and museum samples, including two collected by Darwin and Robert Fitzroy on Floreana in 1835. Our microsatellite results show substantial differentiation between the two extant populations, but our coalescence-based modelling does not indicate long, independent evolutionary histories. One of the populations is highly inbred, but both harbour unique alleles present on Floreana in 1835, suggesting that birds from both islets should be used to establish a single, mixed population on Floreana. Thus, Darwin's mockingbird specimens not only revealed to him a level of variation that suggested speciation following geographical isolation but also, more than 170 years later, return important information to their place of origin for the conservation of their conspecifics.