Characterizing morphological (co)variation using structural equation models: Body size, allometric relationships and evolvability in a house sparrow metapopulation

Long-distance dispersal in the short-distance dispersing house sparrow (Passer domesticus)

The house sparrow (Passer domesticus) is a small passerine known to be highly sedentary. Throughout a 30-year capture-mark-recapture study, we have obtained occasional reports of recoveries far outside our main metapopulation study system, documenting unusually long dispersal distances. Our records constitute the highest occurrence of long-distance dispersal events recorded for this species in Scandinavia. Such long-distance dispersals radically change the predicted distribution of dispersal distances and connectedness for our study metapopulation. Moreover, it reveals a much greater potential for colonization than formerly recorded for the house sparrow, which is an invasive species across four continents. These rare and occasional long-distance dispersal events are challenging to document but may have important implications for the genetic composition of small and isolated populations and for our understanding of dispersal ecology and evolution.

The genetic basis of dispersal in a vertebrate metapopulation

Dispersal affects evolutionary processes by changing population size and genetic composition, influencing the viability and persistence of populations. Investigating which mechanisms underlie variation in dispersal phenotypes and whether populations harbour adaptive potential for dispersal is crucial to understanding the eco-evolutionary dynamics of this important trait. Here, we investigate the genetic architecture of dispersal among successfully recruited individuals in an insular metapopulation of house sparrows. We use an extensive long-term individual-based ecological data set and high-density single-nucleotide polymorphism (SNP) genotypes for over 2500 individuals. We conducted a genome-wide association study (GWAS), and found a relationship between dispersal probability and a SNP located near genes known to regulate circadian rhythm, glycogenesis and exercise performance, among other functions. However, this SNP only explained 3.8% of variance, suggesting that dispersal is a polygenic trait. We then used an animal model to estimate heritable genetic variation (σA 2 ), which composes 10% of the total variation in dispersal probability. Finally, we investigated differences in σA 2 across populations occupying ecologically relevant habitat types (farm vs. non-farm) using a genetic groups animal model. We found different adaptive potentials across habitats, with higher mean breeding value, σA 2 , and heritability for the habitat presenting lower dispersal rates, suggesting also different roles of environmental variation. Our results suggest a complex genetic architecture of dispersal and demonstrate that adaptive potential may be environment dependent in key eco-evolutionary traits. The eco-evolutionary implications of such environment dependence and consequent spatial variation are likely to become ever more important with the increased fragmentation and loss of suitable habitats for many natural populations.

Spatial variation in the evolutionary potential and constraints of basal metabolic rate and body mass in a wild bird

An organism's energy budget is strongly related to resource consumption, performance, and fitness. Hence, understanding the evolution of key energetic traits, such as basal metabolic rate (BMR), in natural populations is central for understanding life-history evolution and ecological processes. Here we used quantitative genetic analyses to study evolutionary potential of BMR in two insular populations of the house sparrow (Passer domesticus). We obtained measurements of BMR and body mass (M_b ) from 911 house sparrows on the islands of Leka and Vega along the coast of Norway. These two populations were the source populations for translocations to create an additional third, admixed 'common garden' population in 2012. With the use of a novel genetic group animal model concomitant with a genetically determined pedigree, we differentiate genetic and environmental sources of variation, thereby providing insight into the effects of spatial population structure on evolutionary potential. We found that the evolutionary potential of BMR was similar in the two source populations, whereas the Vega population had a somewhat higher evolutionary potential of M_b than the Leka population. BMR was genetically correlated with M_b in both populations, and the conditional evolutionary potential of BMR (independent of body mass) was 41% (Leka) and 53% (Vega) lower than unconditional estimates. Overall, our results show that there is potential for BMR to evolve independently of M_b , but that selection on BMR and/or M_b may have different evolutionary consequences in different populations of the same species.

Genetic architecture and heritability of early-life telomere length in a wild passerine

Early-life telomere length (TL) is associated with fitness in a range of organisms. Little is known about the genetic basis of variation in TL in wild animal populations, but to understand the evolutionary and ecological significance of TL it is important to quantify the relative importance of genetic and environmental variation in TL. In this study, we measured TL in 2746 house sparrow nestlings sampled across 20 years and used an animal model to show that there is a small heritable component of early-life TL (h²  = 0.04). Variation in TL among individuals was mainly driven by environmental (annual) variance, but also brood and parental effects. Parent-offspring regressions showed a large maternal inheritance component in TL ( h maternal 2  = 0.44), but no paternal inheritance. We did not find evidence for a negative genetic correlation underlying the observed negative phenotypic correlation between TL and structural body size. Thus, TL may evolve independently of body size and the negative phenotypic correlation is likely to be caused by nongenetic environmental effects. We further used genome-wide association analysis to identify genomic regions associated with TL variation. We identified several putative genes underlying TL variation; these have been inferred to be involved in oxidative stress, cellular growth, skeletal development, cell differentiation and tumorigenesis in other species. Together, our results show that TL has a low heritability and is a polygenic trait strongly affected by environmental conditions in a free-living bird.

Artificial size selection experiment reveals telomere length dynamics and fitness consequences in a wild passerine

Telomere dynamics could underlie life-history trade-offs among growth, size and longevity, but our ability to quantify such processes in natural, unmanipulated populations is limited. We investigated how 4 years of artificial selection for either larger or smaller tarsus length, a proxy for body size, affected early-life telomere length (TL) and several components of fitness in two insular populations of wild house sparrows over a study period of 11 years. The artificial selection was expected to shift the populations away from their optimal body size and increase the phenotypic variance in body size. Artificial selection for larger individuals caused TL to decrease, but there was little evidence that TL increased when selecting for smaller individuals. There was a negative correlation between nestling TL and tarsus length under both selection regimes. Males had longer telomeres than females and there was a negative effect of harsh weather on TL. We then investigated whether changes in TL might underpin fitness effects due to the deviation from the optimal body size. Mortality analyses indicated disruptive selection on TL because both short and long early-life telomeres tended to be associated with the lowest mortality rates. In addition, there was a tendency for a negative association between TL and annual reproductive success, but only in the population where body size was increased experimentally. Our results suggest that natural selection for optimal body size in the wild may be associated with changes in TL during growth, which is known to be linked to longevity in some bird species.

Inbreeding is associated with shorter early-life telomere length in a wild passerine

Inbreeding can have negative effects on survival and reproduction, which may be of conservation concern in small and isolated populations. However, the physiological mechanisms underlying inbreeding depression are not well-known. The length of telomeres, the DNA sequences protecting chromosome ends, has been associated with health or fitness in several species. We investigated effects of inbreeding on early-life telomere length in two small island populations of wild house sparrows (Passer domesticus) known to be affected by inbreeding depression. Using genomic measures of inbreeding we found that inbred nestling house sparrows (n = 371) have significantly shorter telomeres. Using pedigree-based estimates of inbreeding we found a tendency for inbred nestling house sparrows to have shorter telomeres (n = 1195). This negative effect of inbreeding on telomere length may have been complemented by a heterosis effect resulting in longer telomeres in individuals that were less inbred than the population average. Furthermore, we found some evidence of stronger effects of inbreeding on telomere length in males than females. Thus, telomere length may reveal subtle costs of inbreeding in the wild and demonstrate a route by which inbreeding negatively impacts the physiological state of an organism already at early life-history stages.

Evolvability: A Quantitative-Genetics Perspective

The concept of evolvability emerged in the early 1990s and soon became fashionable as a label for different streams of research in evolutionary biology. In evolutionary quantitative genetics, evolvability is defined as the ability of a population to respond to directional selection. This differs from other fields by treating evolvability as a property of populations rather than organisms or lineages and in being focused on quantification and short-term prediction rather than on macroevolution. While the term evolvability is new to quantitative genetics, many of the associated ideas and research questions have been with the field from its inception as biometry. Recent research on evolvability is more than a relabeling of old questions, however. New operational measures of evolvability have opened possibilities for understanding adaptation to rapid environmental change, assessing genetic constraints, and linking micro- and macroevolution.

Spatial structure and dispersal dynamics in a house sparrow metapopulation

The effects of spatial structure on metapopulation dynamics depend upon the interaction between local population dynamics and dispersal, and how this relationship is affected by the geographical isolation and spatial heterogeneity in habitat characteristics. Our aim is to examine how emigration and immigration of house sparrows Passer domesticus in a Norwegian archipelagic metapopulation are affected by key factors predicted by classic metapopulation models to affect dispersal-spatial and temporal variation in population size, inter-island distance, local demography and habitat characteristics. This metapopulation can be divided into two major habitat types: (a) islands closer to the mainland where sparrows breed in colonies on farms, and (b) islands without farms, situated farther away from the mainland where sparrows are exposed to harsher environmental conditions. Dispersal was spatially structured within the metapopulation; there was proportionally and numerically less emigration and immigration involving farm islands, as compared to non-farm islands. Furthermore, emigration and immigration occurred mostly between nearby islands. Moreover, emigration in response to spatial differences in mean population size differed between the habitat types, but populations with large mean received more immigrants in both habitat types. The number of emigrants and immigrants was negatively related to long-term recruit production, which was not the case in non-farm islands. The proportion and number of emigrants was positively related to temporal increases in recruit production on farm islands, however not on non-farm islands. Our results demonstrate that spatial heterogeneity in environmental conditions influences how spatial variation in long-term mean population size, and temporal and spatial variation in recruit production, affects dispersal dynamics. The spatial structure of this metapopulation is therefore best described by a spatially explicit model in which the exchange of individuals within each habitat type is strongly affected by the degree of geographical isolation, population size and recruit production. However, these relationships differed between the two habitat types; non-farm islands showing similarities to a mainland-island model type of structure, whereas farm islands showed features more associated with source-sink or balanced dispersal models. Such differential dispersal dynamics between habitat types are expected to have important consequences for the ecological and evolutionary dynamics within this metapopulation.

Conceptualizing and quantifying body condition using structural equation modelling: A user guide

Body condition is an important concept in behaviour, evolution and conservation, commonly used as a proxy of an individual's performance, for example in the assessment of environmental impacts. Although body condition potentially encompasses a wide range of health state dimensions (nutritional, immune or hormonal status), in practice most studies operationalize body condition using a single (univariate) measure, such as fat storage. One reason for excluding additional axes of variation may be that multivariate descriptors of body condition impose statistical and analytical challenges. Structural equation modelling (SEM) is used in many fields to study questions relating multidimensional concepts, and we here explain how SEM is a useful analytical tool to describe the multivariate nature of body condition. In this 'Research Methods Guide' paper, we show how SEM can be used to resolve different challenges in analysing the multivariate nature of body condition, such as (a) variable reduction and conceptualization, (b) specifying the relationship of condition to performance metrics, (c) comparing competing causal hypothesis and (d) including many pathways in a single model to avoid stepwise modelling approaches. We illustrated the use of SEM on a real-world case study and provided R-code of worked examples as a learning tool. We compared the predictive power of SEM with conventional statistical approaches that integrate multiple variables into one condition variable: multiple regression and principal component analyses. We show that model performance on our dataset is higher when using SEM and led to more accurate and precise estimates compared to conventional approaches. We encourage researchers to consider SEM as a flexible framework to describe the multivariate nature of body condition and thus understand how it affects biological processes, thereby improving the value of body condition proxies for predicting organismal performance. Finally, we highlight that it can be useful for other multidimensional ecological concepts as well, such as immunocompetence, oxidative stress and environmental conditions.

Dispersal in a house sparrow metapopulation: An integrative case study of genetic assignment calibrated with ecological data and pedigree information

Dispersal has a crucial role determining ecoevolutionary dynamics through both gene flow and population size regulation. However, to study dispersal and its consequences, one must distinguish immigrants from residents. Dispersers can be identified using telemetry, capture-mark-recapture (CMR) methods, or genetic assignment methods. All of these methods have disadvantages, such as high costs and substantial field efforts needed for telemetry and CMR surveys, and adequate genetic distance required in genetic assignment. In this study, we used genome-wide 200K Single Nucleotide Polymorphism data and two different genetic assignment approaches (GSI_SIM, Bayesian framework; BONE, network-based estimation) to identify the dispersers in a house sparrow (Passer domesticus) metapopulation sampled over 16 years. Our results showed higher assignment accuracy with BONE. Hence, we proceeded to diagnose potential sources of errors in the assignment results from the BONE method due to variation in levels of interpopulation genetic differentiation, intrapopulation genetic variation and sample size. We show that assignment accuracy is high even at low levels of genetic differentiation and that it increases with the proportion of a population that has been sampled. Finally, we highlight that dispersal studies integrating both ecological and genetic data provide robust assessments of the dispersal patterns in natural populations.

Contrary to vertebrates, less aggressive and more consistent individuals are common in disturbed habitats in the colonial spider Metabus gravidus (Araneae: Araneidae)

Habitat disturbance may affect average behavioural types and consistency/plasticity of behaviour. Studies with solitary vertebrates suggest that human-modified habitats may favour bolder, more aggressive and more plastic individuals. We evaluated whether wild colonial spiders, Metabus gravidus, vary in the magnitude, consistency and plasticity of boldness and aggressiveness between an undisturbed forest and an adjacent urban area in Monteverde, Costa Rica. Repeatability of aggressiveness was high at the disturbed site but moderate at the undisturbed site; repeatability of boldness was low at both sites. Individual and population plasticity was similar between sites for both behaviours. Aggressiveness decreases with increasing colony size at the disturbed site; this trend was not observed at the undisturbed site. Boldness did not change with colony size. In contrast to solitary animals, our results indicate that less aggressive and more consistent colonies may have an advantage living in human-disturbed habitats.

Consistent scaling of inbreeding depression in space and time in a house sparrow metapopulation

Inbreeding may increase the extinction risk of small populations. Yet, studies using modern genomic tools to investigate inbreeding depression in nature have been limited to single populations, and little is known about the dynamics of inbreeding depression in subdivided populations over time. Natural populations often experience different environmental conditions and differ in demographic history and genetic composition, characteristics that can affect the severity of inbreeding depression. We utilized extensive long-term data on more than 3,100 individuals from eight islands in an insular house sparrow metapopulation to examine the generality of inbreeding effects. Using genomic estimates of realized inbreeding, we discovered that inbred individuals had lower survival probabilities and produced fewer recruiting offspring than noninbred individuals. Inbreeding depression, measured as the decline in fitness-related traits per unit inbreeding, did not vary appreciably among populations or with time. As a consequence, populations with more resident inbreeding (due to their demographic history) paid a higher total fitness cost, evidenced by a larger variance in fitness explained by inbreeding within these populations. Our results are in contrast to the idea that effects of inbreeding generally depend on ecological factors and genetic differences among populations, and expand the understanding of inbreeding depression in natural subdivided populations.

Ontogenetic shape changes in the pelvis of the Greater Rhea (Aves, Palaeognathae) and their relationships with cursorial locomotion: a geometric morphometric approach

Knowledge of the ontogenetic pattern of morphological features is essential to improve biological interpretations. The study of morphological features of the pelvic girdle and hind limb apparatus throughout growth is an excellent approach to understand how the skeletal morphology and muscles are interrelated during growth in a bird with a specialized mode of locomotion. The Greater Rhea (Rhea americana) is a large cursorial palaeognathous bird with long legs and powerful musculature. The postnatal shape changes of the pelvis of this bird were studied with geometric morphometric techniques, using landmarks and semilandmarks. In addition, regression analyses were used to explore the association between pelvic shape changes with muscle and body mass. The pelvises of 16 specimens of Rhea americana from 1 month old to adulthood were studied in dorsal and lateral views. Noticeable differences in pelvic shape were noted between ages, particularly in lateral view. In young birds, the pre- and post-acetabular ilium was subequal in length, whereas in adults the pre-acetabular ilium became shorter. In dorsal view, the main shape changes observed were the progressive thinning of both ilium portions and the elongation of the vertex craniolateralis ilii from chicks to adulthood. In this view, the only clear differentiation was between young and adult birds. Shape differences were influenced by body mass and pelvic muscles; the post-acetabular muscle mass explained the highest percentage of the variation. The specialized locomotion of Greater Rhea is reflected in their pelvic musculoskeletal system, in which the change to a longer post-acetabular ilium correlates with the growth of the powerful post-acetabular muscles. The actions of these muscles provide the necessary strength to support the body mass, minimize the body swinging movements and propel the body forward during locomotion. Bone morphology is affected by the forces produced by body mass and the muscle activity, demonstrating the presence of common growth mechanisms, which are primordial and gave rise to a functional and properly proportioned adult.

Multi-generational genetic consequences of reinforcement in a bird metapopulation

Translocation of conspecific individuals to reduce extinction risk of small, isolated populations and prevent genetic depletion is a powerful tool in conservation biology. An important question is how the translocated individuals influence the long-term genetic composition of the recipient population. Here, we experimentally reinforced a house sparrow (Passer domesticus) population, and examined the impact of this translocation on allele frequencies, levels of heterozygosity and genetic differentiation over six cohorts. We found no permanent increase in the mean number of alleles across loci or levels of observed heterozygosity, but a few alleles private to the translocated individuals remained in the population and we found a short-term increase in heterozygosity. Consequently, genetic differentiation of the recipient population compared to the genetic composition prior to reinforcement was small. The limited genetic impact was due to combined effects of a small probability of establishment and low mating success for the translocated individuals, together with increased genetic drift in the recipient population. Our findings emphasize the importance of selection and genetic drift as forces that may decrease the genetic contribution of reinforcement, especially in small populations. Conservation managers should aim to improve habitat quality in the recipient population to reduce genetic drift following translocation and thereby avoid the need for continued reinforcement. Furthermore, by facilitating establishment success and selecting individuals expected to have high mating success, possibly indicated by sexually selected traits, genetic contribution of released individuals is increased which in turn will decrease reproductive skew and genetic drift.

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.