An ecologist’s guide to the animal model

Population Dynamic Consequences of Context-Dependent Trade-Offs across Life Histories

AbstractTrade-offs are central to life history theory and play a role in driving life history diversity. They arise from a finite amount of resources that need to be allocated among different functions by an organism. Yet covariation of demographic rates among individuals frequently do not reflect allocation trade-offs because of variation in resource acquisition. The covariation of traits among individuals can thus vary with the environment and often increases in benign environments. Surprisingly, little is known about how such context-dependent expression of trade-offs among individuals affect population dynamics across species with different life histories. To study their influence on population stability, we develop an individual-based simulation where covariation in demographic rates varies with the environment. We use it to simulate population dynamics for various life histories across the slow-fast pace-of-life continuum. We found that the population dynamics of slower life histories are relatively more sensitive to changes in covariation, regardless of the trade-off considered. Additionally, we found that the impact on population stability depends on which trade-off is considered, with opposite effects of intraindividual and intergenerational trade-offs. Last, the expression of different trade-offs can feed back to influence generation time through selection acting on individual heterogeneity within cohorts, ultimately affecting population dynamics.

Genetically correlated leaf tensile and morphological traits are driven by growing season length in a widespread perennial grass

PREMISE

Leaf tensile resistance, a leaf's ability to withstand pulling forces, is an important determinant of plant ecological strategies. One potential driver of leaf tensile resistance is growing season length. When growing seasons are long, strong leaves, which often require more time and resources to construct than weak leaves, may be more advantageous than when growing seasons are short. Growing season length and other ecological conditions may also impact the morphological traits that underlie leaf tensile resistance.

METHODS

To understand variation in leaf tensile resistance, we measured size-dependent leaf strength and size-independent leaf toughness in diverse genotypes of the widespread perennial grass Panicum virgatum (switchgrass) in a common garden. We then used quantitative genetic approaches to estimate the heritability of leaf tensile resistance and whether there were genetic correlations between leaf tensile resistance and other morphological traits.

RESULTS

Leaf tensile resistance was positively associated with aboveground biomass (a proxy for fitness). Moreover, both measures of leaf tensile resistance exhibited high heritability and were positively genetically correlated with leaf lamina thickness and leaf mass per area (LMA). Leaf tensile resistance also increased with the growing season length in the habitat of origin, and this effect was mediated by both LMA and leaf thickness.

CONCLUSIONS

Differences in growing season length may promote selection for different leaf lifespans and may explain existing variation in leaf tensile resistance in P. virgatum. In addition, the high heritability of leaf tensile resistance suggests that P. virgatum will be able to respond to climate change as growing seasons lengthen.

Fathers and sons, mothers and daughters: Sex-specific genetic architecture for fetal testosterone in a wild mammal

Testosterone plays a critical role in mediating fitness-related traits in many species. Although it is highly responsive to environmental and social conditions, evidence from several species show a heritable component to its individual variation. Despite the known effects that in utero testosterone exposure have on adult fitness, the heritable component of individual testosterone variation in fetuses is mostly unexplored. Furthermore, testosterone has sex-differential effects on fetal development, i.e., a specific level may be beneficial for male fetuses but detrimental for females, producing sexual conflict. Such sexual conflict may be resolved by the evolution of a sex-specific genetic architecture of the trait. Here, we quantified fetal testosterone levels in a wild species, free-ranging nutrias (Myocastor coypus) using hair-testing and estimated testosterone heritability between parent and offspring from the same and opposite sex. We found that in utero accumulated hair testosterone levels were heritable between parents and offspring of the same sex. Moreover, there was a low additive genetic covariance between the sexes, and a low cross-sex genetic correlation, suggesting a potential for sex-specific trait evolution, expressed early on, in utero.

Vertically inherited microbiota and environment modifying behaviours conceal genetic variation in dung beetle life history

Diverse organisms actively manipulate their (sym)biotic and physical environment in ways that feed back on their own development. However, the degree to which these processes affect microevolution remains poorly understood. The gazelle dung beetle both physically modifies its ontogenetic environment and structures its biotic interactions through vertical symbiont transmission. By experimentally eliminating (i) physical environmental modifications and (ii) the vertical inheritance of microbes, we assess how environment modifying behaviour and microbiome transmission shape heritable variation and evolutionary potential. We found that depriving larvae of symbionts and environment modifying behaviours increased additive genetic variance and heritability for development time but not body size. This suggests that larvae's ability to manipulate their environment has the potential to modify heritable variation and to facilitate the accumulation of cryptic genetic variation. This cryptic variation may become released and selectable when organisms encounter environments that are less amenable to organismal manipulation or restructuring. Our findings also suggest that intact microbiomes, which are commonly thought to increase genetic variation of their hosts, may instead reduce and conceal heritable variation. More broadly, our findings highlight that the ability of organisms to actively manipulate their environment may affect the potential of populations to evolve when encountering novel, stressful conditions.

Transposons in the Williams-Beuren Syndrome Critical Region are Associated with Social Behavior in Assistance Dogs

A strong signature of selection in the domestic dog genome is found in a five-megabase region of chromosome six in which four structural variants derived from transposons have previously been associated with human-oriented social behavior, such as attentional bias to social stimuli and social interest in strangers. To explore these genetic associations in more phenotypic detail-as well as their role in training success in a specialized assistance dog program-we genotyped 1001 assistance dogs from Canine Companions for Independence®, including both successful graduates and dogs released from the training program for behaviors incompatible with their working role. We collected phenotypes on each dog using puppy-raiser questionnaires, trainer questionnaires, and both cognitive and behavioral tests. Using Bayesian mixed models, we found strong associations (95% credibility intervals excluding zero) between genotypes and certain behavioral measures, including separation-related problems, aggression when challenged or corrected, and reactivity to other dogs. Furthermore, we found moderate differences in the genotypes of dogs who graduated versus those who did not; insertions in GTF2I showed the strongest association with training success (β = 0.23, CI95% = - 0.04, 0.49), translating to an odds-ratio of 1.25 for one insertion. Our results provide insight into the role of each of these four transposons in canine sociability and may inform breeding and training practices for working dog organizations. Furthermore, the observed importance of the gene GTF2I supports the emerging consensus that variation in GTF2I genotypes and expression have important consequences for social behavior broadly.

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

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

Evolutionary relationships between metabolism and behaviour require genetic correlations

As selection acts on multivariate phenotypes, the evolution of traits within populations not only depends on the genetic basis of each trait, but also on the genetic relationships among traits. As metabolic rate is often related to vital traits such as growth, physiology and behaviour, its variation and evolution is expected to have important repercussions on individual fitness. However, the majority of the correlations between metabolic rate and other traits has been based on phenotypic correlations, while genetic correlations, basis for indirect selection and evolution, have been overlooked. Using a case study, we explore the importance of properly estimating genetic correlations to understand and predict evolution of multivariate phenotypes. We show that selection on metabolic traits could result in indirect selection mainly on growth-related traits, owing to strong genetic correlations, but not on swimming or risk-taking and sociability behaviour even if they covary phenotypically. While phenotypic correlation can inform about genetic correlation direction, caution is needed in predicting the magnitude of genetic correlation. Therefore, even though phenotypic correlations among physiological and behavioural traits could be useful, deriving evolutionary conclusions based purely on them is not robust. In short, proper estimation of genetic correlations is needed when predicting evolutionary consequences. This article is part of the theme issue 'The evolutionary significance of variation in metabolic rates'.

How Can Genomics Help or Hinder Wildlife Conservation?

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

Heritability and developmental plasticity of growth in an oviparous lizard

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

New insight into social relationships in dairy cows and how time of birth, parity, and relatedness affect spatial interactions later in life

Social interactions between cows play a fundamental role in the daily activities of dairy cattle. Real-time location systems provide on a continuous and automated basis information about the position of individual cows inside barns, offering a valuable opportunity to monitor dyadic social contacts. Understanding dyadic social interactions could be applied to enhance the stability of the social structure promoting animal welfare and to model disease transmission in dairy cattle. This study aimed to identify the effect of different cow characteristics on the likelihood of the formation and persistence of social contacts in dairy cattle. The individual position of the lactating cows was automatically collected once per second for 2 wk, using an ultra-wideband system on a Swedish commercial farm consisting of almost 200 dairy cows inside a freestall barn. Social networks were constructed using the position data of 149 cows with available information on all characteristics during the study period. Social contacts were considered as a binary variable indicating whether a cow pair was within 2.5 m of each other for at least 10 min per day. The role of cow characteristics in social networks was studied by applying separable temporal exponential random graph models. Our results revealed that cows of the same parity interacted more consistently, as well as those born within 7 d of each other or closely related by pedigree. The repeatability of the topological parameters indicated a consistent short-term stability of the individual animal roles within the social network structure. Additional research is required to elucidate the underlying mechanisms governing the long-term evolution of social contacts among dairy cattle and to investigate the relationship between these networks and the transmission of diseases in the dairy cattle population.

Estimation of genetic parameters for growth traits and kleiber ratio in dorper sheep breed

The study aimed to estimate genetic and phenotypic parameters for growth traits and Kleiber ratio in Dorper sheep breed. Data used in this study were collected over 12 years (2012-2023) at Debre Birhan Agricultural Research Center sheep research station in Ethiopia. Studied traits were body weights at birth (WT0), weaning (WT3), six month (WT6), and yearling (WT12) age; average daily gains from birth to weaning (ADG0-3), from weaning to six months (ADG3-6), from six months to yearling (ADG6-12); and Kleiber ratios from birth to weaning (KR1) and from weaning to six months (KR2). The (co)variance components were estimated with different animal models using Average Information Restricted Maximum Likelihood (AI-REML) procedure. The best-fitted model for each trait was determined using likelihood ratio tests. Phenotypic performance for WT3, WT6, WT12, ADG0-3 and ADG3-6 showed a decline trend at a rate of 0.216 kg, 0.794 kg, 0.671 kg, 2.601 g and 4.865 g over years respectively. However, WT3, WT6, WT12, ADG0-3 showed a positive genetic improvement trend at a rate of 0.029 kg, 0.043 kg, 0.049 kg and 0.257 g over years respectively. Year of birth had a significant effect (P < 0.001) on all studied traits. Model including direct genetic as well as permanent environmental effect (Model 2) was chosen as the most appropriate model for WT0. Model which included only direct genetic effect (Model 1) was the best-fit model for all other studied traits. Direct heritability estimates based on suitable models were 0.07 ± 0.06, 0.11 ± 0.06, 0.09 ± 0.07, 0.11 ± 0.09 0.11 ± 0.06, 0.00 ± 0.04, 0.15 ± 0.07 and 0.00 ± 0.04 for WT0, WT3, WT6, WT12, ADG0-3, ADG3-6, KR1 and KR2 respectively. The variance ratio for the permanent environmental effect was 0.13 ± 0.04 for WT0. Genetic correlations among the traits ranged from negative (-0.39) for WT0-KR1 to high (0.99) for WT3-ADG0-3 and phenotypic correlations ranged from negative (-0.31) for WT0-KR1 to high (0.98) for WT3-ADG0-3. The low direct heritability estimates for the studied traits indicated that genetic improvement by direct selection might be difficult. Further investigation for the unexpected declined trend of phenotypic performance over years need to be required.

Plant physical defenses contribute to a latitudinal gradient in resistance to insect herbivory within a widespread perennial grass

PREMISE

Herbivore pressure can vary across the range of a species, resulting in different defensive strategies. If herbivory is greater at lower latitudes, plants may be better defended there, potentially driving a latitudinal gradient in defense. However, relationships that manifest across the entire range of a species may be confounded by differences within genetic subpopulations, which may obscure the drivers of these latitudinal gradients.

METHODS

We grew plants of the widespread perennial grass Panicum virgatum in a common garden that included genotypes from three genetic subpopulations spanning an 18.5° latitudinal gradient. We then assessed defensive strategies of these plants by measuring two physical resistance traits-leaf mass per area (LMA) and leaf ash, a proxy for silica-and multiple measures of herbivory by caterpillars of the generalist herbivore fall armyworm (Spodoptera frugiperda).

RESULTS

Across all genetic subpopulations, low-latitude plants experienced less herbivory than high-latitude plants. Within genetic subpopulations, however, this relationship was inconsistent-the most widely distributed and phenotypically variable subpopulation (Atlantic) exhibited more consistent latitudinal trends than either of the other two subpopulations. The two physical resistance traits, LMA and leaf ash, were both highly heritable and positively associated with resistance to different measures of herbivory across all subpopulations, indicating their importance in defense against herbivores. Again, however, these relationships were inconsistent within subpopulations.

CONCLUSIONS

Defensive gradients that occur across the entire species range may not arise within localized subpopulations. Thus, identifying the drivers of latitudinal gradients in herbivory defense may depend on adequately sampling the diversity within a species.

Reproductive biology of Gazella arabica: Predictors of offspring weight and short- and long-term offspring survival

Reproductive traits are central to organismal fitness, and so the factors influencing patterns of reproduction and offspring survival are at the heart of biology. Making use of breeding data collected over 16 years at the King Khalid Wildlife Research Centre in Saudi Arabia, we investigated the reproductive biology of Arabian gazelles Gazella arabica. Offspring survival was mainly a function of birth weight, with heavier offspring having higher survival rates than lighter offspring. However, while sons were heavier than daughters, daughters had higher survival rates. We could not find evidence that giving birth to sons negatively impacts offspring weight in the following year. We uncovered large narrow-sense heritability (h2) in offspring weight at birth, while maternal effects (m2) on birth weight were of lesser importance. However, maternal effects on offspring survival were strong until weaning age, while paternal effects dominated survival to sexual maturity and first reproduction. We propose that variation in maternal postnatal care might overshadow the effects of maternal inheritance of birth weights, while the overall strong heritability of weight at birth and the paternal effects on survival illustrates strong variance in sire fitness based on genetic quality, suggesting a role for sexual selection by female mate choice in wild populations.

Genetic and context-specific effects on individual inhibitory control performance in the guppy (Poecilia reticulata)

Among-individual variation in cognitive traits, widely assumed to have evolved under adaptive processes, is increasingly being demonstrated across animal taxa. As variation among individuals is required for natural selection, characterizing individual differences and their heritability is important to understand how cognitive traits evolve. Here, we use a quantitative genetic study of wild-type guppies repeatedly exposed to a 'detour task' to test for genetic variance in the cognitive trait of inhibitory control. We also test for genotype-by-environment interactions (GxE) by testing related fish under alternative experimental treatments (transparent vs. semi-transparent barrier in the detour-task). We find among-individual variation in detour task performance, consistent with differences in inhibitory control. However, analysis of GxE reveals that heritable factors only contribute to performance variation in one treatment. This suggests that the adaptive evolutionary potential of inhibitory control (and/or other latent variables contributing to task performance) may be highly sensitive to environmental conditions. The presence of GxE also implies that the plastic response of detour task performance to treatment environment is genetically variable. Our results are consistent with a scenario where variation in individual inhibitory control stems from complex interactions between heritable and plastic components.

Repeatability and heritability of inhibitory control performance in wild toutouwai (Petroica longipes)

Despite increasing interest in the evolution of inhibitory control, few studies have examined the validity of widespread testing paradigms, the long-term repeatability and the heritability of this cognitive ability in the wild. We investigated these aspects in the inhibitory control performance of wild toutouwai (North Island robin; Petroica longipes), using detour and reversal learning tasks. We assessed convergent validity by testing whether individual performance correlated across detour and reversal learning tasks. We then further evaluated task validity by examining whether individual performance was confounded by non-cognitive factors. We tested a subset of subjects twice in each task to estimate the repeatability of performance across a 1-year period. Finally, we used a population pedigree to estimate the heritability of task performance. Individual performance was unrelated across detour and reversal learning tasks, indicating that these measured different cognitive abilities. Task performance was not influenced by body condition, boldness or prior experience, and showed moderate between-year repeatability. Yet despite this individual consistency, we found no evidence that task performance was heritable. Our findings suggest that detour and reversal learning tasks measure consistent individual differences in distinct forms of inhibitory control in toutouwai, but this variation may be environmentally determined rather than genetic.

How and Why Does Metabolism Scale with Body Mass?

Most explanations for the relationship between body size and metabolism invoke physical constraints; such explanations are evolutionarily inert, limiting their predictive capacity. Contemporary approaches to metabolic rate and life history lack the pluralism of foundational work. Here, we call for reforging of the lost links between optimization approaches and physiology.

Early-life telomeres are influenced by environments acting at multiple temporal and spatial scales

An individual's telomere length early in life may reflect or contribute to key life-history processes sensitive to environmental variation. Yet, the relative importance of genetic and environmental factors in shaping early-life telomere length is not well understood as it requires samples collected from multiple generations with known developmental histories. We used a confirmed pedigree and conducted an animal model analysis of telomere lengths obtained from nestling house sparrows (Passer domesticus) sampled over a span of 22 years. We found significant additive genetic variation for early-life telomere length, but it comprised a small proportion (9%) of the total biological variation. Three sources of environmental variation were important: among cohorts, among-breeding attempts within years, and among nestmates. The magnitude of variation among breeding attempts and among nestmates also differed by cohort, suggesting that interactive effects of environmental factors across time or spatial scales were important, yet we were unable to identify the specific causes of these interactions. The mean amount of precipitation during the breeding season positively predicted telomere length, but neither weather during a given breeding attempt nor date in the breeding season contributed to an offspring's telomere length. At the level of individual nestlings, offspring sex, size and mass at 10 days of age also did not predict telomere length. Environmental effects appear especially important in shaping early-life telomere length in some species, and more focus on how environmental factors that interact across scales may help to explain some of the variation observed among studies.

Early life drought predicts components of adult body size in wild female baboons

OBJECTIVES

In many taxa, adverse early-life environments are associated with reduced growth and smaller body size in adulthood. However, in wild primates, we know very little about whether, where, and to what degree trajectories are influenced by early adversity, or which types of early adversity matter most. Here, we use parallel-laser photogrammetry to assess inter-individual predictors of three measures of body size (leg length, forearm length, and shoulder-rump length) in a population of wild female baboons studied since birth.

MATERIALS AND METHODS

Using >2000 photogrammetric measurements of 127 females, we present a cross-sectional growth curve of wild female baboons (Papio cynocephalus) from juvenescence through adulthood. We then test whether females exposed to several important sources of early-life adversity-drought, maternal loss, low maternal rank, or a cumulative measure of adversity-were smaller for their age than females who experienced less adversity. Using the "animal model," we also test whether body size is heritable in this study population.

RESULTS

Prolonged early-life drought predicted shorter limbs but not shorter torsos (i.e., shoulder-rump lengths). Our other measures of early-life adversity did not predict variation in body size. Heritability estimates for body size measures were 36%-67%. Maternal effects accounted for 13%-17% of the variance in leg and forearm length, but no variance in torso length.

DISCUSSION

Our results suggest that baboon limbs, but not torsos, grow plastically in response to maternal effects and energetic early-life stress. Our results also reveal considerable heritability for all three body size measures in this study population.

Relative abundance data can misrepresent heritability of the microbiome

BACKGROUND

Host genetics can shape microbiome composition, but to what extent it does, remains unclear. Like any other complex trait, this important question can be addressed by estimating the heritability (h2) of the microbiome-the proportion of variance in the abundance in each taxon that is attributable to host genetic variation. However, unlike most complex traits, microbiome heritability is typically based on relative abundance data, where taxon-specific abundances are expressed as the proportion of the total microbial abundance in a sample.

RESULTS

We derived an analytical approximation for the heritability that one obtains when using such relative, and not absolute, abundances, based on an underlying quantitative genetic model for absolute abundances. Based on this, we uncovered three problems that can arise when using relative abundances to estimate microbiome heritability: (1) the interdependency between taxa can lead to imprecise heritability estimates. This problem is most apparent for dominant taxa. (2) Large sample size leads to high false discovery rates. With enough statistical power, the result is a strong overestimation of the number of heritable taxa in a community. (3) Microbial co-abundances lead to biased heritability estimates.

CONCLUSIONS

We discuss several potential solutions for advancing the field, focusing on technical and statistical developments, and conclude that caution must be taken when interpreting heritability estimates and comparing values across studies. Video Abstract.

Trait-Specific Indirect Effects Underlie Variation in the Response of Spiders to Cannibalistic Social Partners

AbstractIn cannibalistic species, selection to avoid conspecifics may stem from the need to avoid being eaten or to avoid competition. Individuals may thus use conspecific cues to modulate their behavior to such threats. Yet the nature of variation for such cues remains elusive. Here, we use a half-sib/full-sib design to evaluate the contribution of (indirect) genetic or environmental effects to the behavioral response of the cannibalistic wolf spider Lycosa fasciiventris (Dufour, 1835) toward conspecific cues. Spiders showed variation in relative occupancy time, activity, and velocity on patches with or without conspecific cues, but direct genetic variance was found only for occupancy time. These three traits were correlated and could be lumped in a principal component: spiders spending more time in patches with conspecific cues moved less and more slowly in those areas. Genetic and/or environmental components of carapace width and weight loss in the social partner, which may reflect the quality and/or quantity of cues produced, were significantly correlated with this principal component, with larger partners causing focal individuals to move more slowly. Therefore, environmental and genetic trait variation in social partners may maintain trait diversity in focal individuals, even in the absence of direct genetic variation.

Phylogenetic meta-analysis reveals system-specific behavioural type-behavioural predictability correlations

The biological significance of behavioural predictability (environment-independent within-individual behavioural variation) became accepted recently as an important part of an individual's behavioural strategy besides behavioural type (individual mean behaviour). However, we do not know how behavioural type and predictability evolve. Here, we tested different evolutionary scenarios: (i) the two traits evolve independently (lack of correlations) and (ii) the two traits' evolution is constrained (abundant correlations) due to either (ii/a) proximate constraints (direction of correlations is similar) or (ii/b) local adaptations (direction of correlations is variable). We applied a set of phylogenetic meta-analyses based on 93 effect sizes across 44 vertebrate and invertebrate species, focusing on activity and risk-taking. The general correlation between behavioural type and predictability did not differ from zero. Effect sizes for correlations showed considerable heterogeneity, with both negative and positive correlations occurring. The overall absolute (unsigned) effect size was high (Zr = 0.58), and significantly exceeded the null expectation based on randomized data. Our results support the adaptive scenario: correlations between behavioural type and predictability are abundant in nature, but their direction is variable. We suggest that the evolution of these behavioural components might be constrained in a system-specific way.

Early life adversity has sex-dependent effects on survival across the lifespan in rhesus macaques

Exposure to adversity during early life is linked to lasting detrimental effects on evolutionary fitness across many taxa. However, due to the challenges of collecting longitudinal data, especially in species where one sex disperses, direct evidence from long-lived species remains relatively scarce. Here we test the effects of early life adversity on male and female longevity in a free-ranging population of rhesus macaques (Macaca mulatta) at Cayo Santiago, Puerto Rico. We leveraged six decades of data to quantify the relative importance of ten forms of early life adversity for 6,599 macaques (3,230 male, 3,369 female), with a smaller sample size (N=299) for one form of adversity (maternal social isolation) which required high-resolution behavioral data. We found that individuals who experienced more early life adversity died earlier than those who experienced less adversity. Mortality risk was highest during early life, defined as birth to four years old, suggesting acute survival effects of adversity, but heightened mortality risk was also present in macaques who survived to adulthood. Females and males were affected differently by some forms of adversity, and these differences might be driven by varying energetic demands, female philopatry, and male dispersal. By leveraging data on thousands of macaques collected over decades, our results show that the fitness consequences of early life adversity are not uniform across individuals but vary as a function of the type of adversity, timing, and social context, and thus contribute to our limited but growing understanding of the evolution of early life sensitivities in long-lived species.

A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance

Next-generation sequencing (NGS) and metabarcoding approaches are increasingly applied to wild animal populations, but there is a disconnect between the widely applied generalized linear mixed model (GLMM) approaches commonly used to study phenotypic variation and the statistical toolkit from community ecology typically applied to metabarcoding data. Here, we describe the suitability of a novel GLMM-based approach for analyzing the taxon-specific sequence read counts derived from standard metabarcoding data. This approach allows decomposition of the contribution of different drivers to variation in community composition (e.g., age, season, individual) via interaction terms in the model random-effects structure. We provide guidance to implementing this approach and show how these models can identify how responsible specific taxonomic groups are for the effects attributed to different drivers. We applied this approach to two cross-sectional data sets from the Soay sheep population of St. Kilda. GLMMs showed agreement with dissimilarity-based approaches highlighting the substantial contribution of age and minimal contribution of season to microbiota community compositions, and simultaneously estimated the contribution of other technical and biological factors. We further used model predictions to show that age effects were principally due to increases in taxa of the phylum Bacteroidetes and declines in taxa of the phylum Firmicutes. This approach offers a powerful means for understanding the influence of drivers of community structure derived from metabarcoding data. We discuss how our approach could be readily adapted to allow researchers to estimate contributions of additional factors such as host or microbe phylogeny to answer emerging questions surrounding the ecological and evolutionary roles of within-host communities. IMPORTANCE NGS and fecal metabarcoding methods have provided powerful opportunities to study the wild gut microbiome. A wealth of data is, therefore, amassing across wild systems, generating the need for analytical approaches that can appropriately investigate simultaneous factors at the host and environmental scale that determine the composition of these communities. Here, we describe a generalized linear mixed-effects model (GLMM) approach to analyze read count data from metabarcoding of the gut microbiota, allowing us to quantify the contributions of multiple host and environmental factors to within-host community structure. Our approach provides outputs that are familiar to a majority of field ecologists and can be run using any standard mixed-effects modeling packages. We illustrate this approach using two metabarcoding data sets from the Soay sheep population of St. Kilda investigating age and season effects as worked examples.

Rhesus macaque social functioning is paternally, but not maternally, inherited by sons: potential implications for autism

BACKGROUND

Quantitative autistic traits are common, heritable, and continuously distributed across the general human population. Patterns of autistic traits within families suggest that more complex mechanisms than simple Mendelian inheritance-in particular, parent of origin effects-may be involved. The ideal strategy for ascertaining parent of origin effects is by half-sibling analysis, where half-siblings share one, but not both, parents and each individual belongs to a unique combination of paternal and maternal half-siblings. While this family structure is rare in humans, many of our primate relatives, including rhesus macaques, have promiscuous breeding systems that consistently produce paternal and maternal half-siblings for a given index animal. Rhesus macaques, like humans, also exhibit pronounced variation in social functioning.

METHODS

Here we assessed differential paternal versus maternal inheritance of social functioning in male rhesus macaque offspring (N = 407) using ethological observations and ratings on a reverse-translated quantitative autistic trait measurement scale. Restricted Maximum Likelihood mixed models with unbounded variance estimates were used to estimate the variance components needed to calculate the genetic contribution of parents as the proportion of phenotypic variance (σ2P) between sons that could uniquely be attributed to their shared genetics (σ2g), expressed as σ2g/σ2P (or the proportion of phenotypic variance attributable to genetic variance), as well as narrow sense heritability (h2).

RESULTS

Genetic contributions and heritability estimates were strong and highly significant for sons who shared a father but weak and non-significant for sons who shared a mother. Importantly, these findings were detected using the same scores from the same sons in the same analysis, confirmed when paternal and maternal half-siblings were analyzed separately, and observed with two methodologically distinct behavioral measures. Finally, genetic contributions were similar for full-siblings versus half-siblings that shared only a father, further supporting a selective paternal inheritance effect.

LIMITATIONS

These data are correlational by nature. A larger sample that includes female subjects, enables deeper pedigree assessments, and supports molecular genetic analyses is warranted.

CONCLUSIONS

Rhesus macaque social functioning may be paternally, but not maternally, inherited by sons. With continued investigation, this approach may yield important insights into sex differences in autism's genetic liability.

Plastic responses to warmer climates: a semi-natural experiment on lizard populations

Facing warming environments, species can exhibit plastic or microevolutionary changes in their thermal physiology to adapt to novel climates. Here, using semi-natural mesocosms, we experimentally investigated over two successive years whether a 2°C-warmer climate produces selective and inter- and intragenerational plastic changes in the thermal traits (preferred temperature and dorsal coloration) of the lizard Zootoca vivipara. In a warmer climate, the dorsal darkness, dorsal contrast, and preferred temperature of adults plastically decreased and covariances between these traits were disrupted. While selection gradients were overall weak, selection gradients for darkness were slightly different between climates and in the opposite direction to plastic changes. Contrary to adults, male juveniles were darker in warmer climates either through plasticity or selection and this effect was strengthened by intergenerational plasticity when juveniles' mothers also experienced warmer climates. While the plastic changes in adult thermal traits alleviate the immediate overheating costs of warming, its opposite direction to selective gradients and to juveniles' phenotypic responses may slow down evolutionary shifts toward phenotypes that are better adapted to future climates. Our study demonstrates the importance of considering inter- and intragenerational plasticity along with selective processes to better understand adaptation and population dynamics in light of climate change.

Genotypes selected for early and late avian lay date differ in their phenotype, but not fitness, in the wild

Global warming has shifted phenological traits in many species, but whether species are able to track further increasing temperatures depends on the fitness consequences of additional shifts in phenological traits. To test this, we measured phenology and fitness of great tits (Parus major) with genotypes for extremely early and late egg lay dates, obtained from a genomic selection experiment. Females with early genotypes advanced lay dates relative to females with late genotypes, but not relative to nonselected females. Females with early and late genotypes did not differ in the number of fledglings produced, in line with the weak effect of lay date on the number of fledglings produced by nonselected females in the years of the experiment. Our study is the first application of genomic selection in the wild and led to an asymmetric phenotypic response that indicates the presence of constraints toward early, but not late, lay dates.

Genetic rescue from protected areas is modulated by migration, hunting rate, and timing of harvest

In terrestrial and marine ecosystems, migrants from protected areas may buffer the risk of harvest-induced evolutionary changes in exploited populations that face strong selective harvest pressures. Understanding the mechanisms favoring genetic rescue through migration could help ensure evolutionarily sustainable harvest outside protected areas and conserve genetic diversity inside those areas. We developed a stochastic individual-based metapopulation model to evaluate the potential for migration from protected areas to mitigate the evolutionary consequences of selective harvest. We parameterized the model with detailed data from individual monitoring of two populations of bighorn sheep subjected to trophy hunting. We tracked horn length through time in a large protected and a trophy-hunted populations connected through male breeding migrations. We quantified and compared declines in horn length and rescue potential under various combinations of migration rate, hunting rate in hunted areas and temporal overlap in timing of harvest and migrations, which affects the migrants' survival and chances to breed within exploited areas. Our simulations suggest that the effects of size-selective harvest on male horn length in hunted populations can be dampened or avoided if harvest pressure is low, migration rate is substantial, and migrants leaving protected areas have a low risk of being shot. Intense size-selective harvest impacts the phenotypic and genetic diversity in horn length, and population structure through changes in proportions of large-horned males, sex ratio and age structure. When hunting pressure is high and overlaps with male migrations, effects of selective removal also emerge in the protected population, so that instead of a genetic rescue of hunted populations, our model predicts undesirable effects inside protected areas. Our results stress the importance of a landscape approach to management, to promote genetic rescue from protected areas and limit ecological and evolutionary impacts of harvest on both harvested and protected populations.

Tests for associations between sexual dimorphism and patterns of quantitative genetic variation in the water strider, Aquarius remigis

The evolution of sexual dimorphisms requires divergence between sexes in the evolutionary trajectories of the traits involved. Discerning how genetic architecture could facilitate such divergence has proven challenging because of the difficulty in estimating non-additive and sex-linked genetic variances using traditional quantitative genetic designs. Here we use a three-generation, double-first-cousin pedigree design to estimate additive, sex-linked and dominance (co)variances for 12 traits in the water strider, Aquarius remigis. Comparisons among these traits, which have size ratios ranging from 1 to 5 (larger/smaller), allow us to ask if sexual dimorphisms are associated with characteristic patterns of quantitative genetic variation. We frame our analysis around three main questions, derived from existing theory and empirical evidence: Are sexual dimorphisms associated with (1) lower additive inter-sex genetic correlations, (2) higher proportions of sex-linked variance, or (3) differences between sexes in autosomal additive and dominance genetic variances? For questions (1) and (2), we find weak and non-significant trends in the expected directions, which preclude definitive conclusions. However, in answer to question (3), we find strong evidence for a positive relationship between sexual dimorphism and differences between sexes in proportions of autosomal dominance variance. We also find strong interactions among the three genetic components indicating that their relative influence differs among traits and between sexes. These results highlight the need to include all three components of genetic (co)variance in both theoretical evolutionary models and empirical estimations of the genetic architecture of dimorphic traits.

Genetic improvement of speed across distance categories in thoroughbred racehorses in Great Britain

Several studies over recent decades have reported a lack of contemporary improvement in thoroughbred racehorse speed, despite apparent additive genetic variance and putatively strong selection. More recently, it has been shown that some phenotypic improvement is ongoing, but rates are low in general and particularly so over longer distances. Here we used pedigree-based analysis of 692,534 records from 76,960 animals to determine whether these phenotypic trends are underpinned by genetic selection responses, and to evaluate the potential for more rapid improvement. We show that thoroughbred speed in Great Britain is only weakly heritable across sprint (h² = 0.124), middle-distance (h² = 0.122) and long-distance races (h² = 0.074), but that mean predicted breeding values are nonetheless increasing across cohorts born between 1995 and 2012 (and racing from 1997 to 2014). For all three race distance categories, estimated rates of genetic improvement are statistically significant and also greater than can be explained by drift. Taken together our results show genetic improvement for thoroughbred speed is ongoing but slow, likely due to a combination of long generation times and low heritabilities. Additionally, estimates of realised selection intensities raises the possibility that the contemporary selection emerging from the collective actions of horse breeders is weaker than previously assumed, particularly over long distances. We suggest that unmodelled common environment effects may have upwardly biased estimates of heritability, and thus expected selection response, previously.

Shared community effects and the non-genetic maternal environment shape cortisol levels in wild chimpanzees

Mechanisms of inheritance remain poorly defined for many fitness-mediating traits, especially in long-lived animals with protracted development. Using 6,123 urinary samples from 170 wild chimpanzees, we examined the contributions of genetics, non-genetic maternal effects, and shared community effects on variation in cortisol levels, an established predictor of survival in long-lived primates. Despite evidence for consistent individual variation in cortisol levels across years, between-group effects were more influential and made an overwhelming contribution to variation in this trait. Focusing on within-group variation, non-genetic maternal effects accounted for 8% of the individual differences in average cortisol levels, significantly more than that attributable to genetic factors, which was indistinguishable from zero. These maternal effects are consistent with a primary role of a shared environment in shaping physiology. For chimpanzees, and perhaps other species with long life histories, community and maternal effects appear more relevant than genetic inheritance in shaping key physiological traits.

Rapid evolution during climate change: demographic and genetic constraints on adaptation to severe drought

Populations often vary in their evolutionary responses to a shared environmental perturbation. A key hurdle in building more predictive models of rapid evolution is understanding this variation-why do some populations and traits evolve while others do not? We combined long-term demographic and environmental data, estimates of quantitative genetic variance components, a resurrection experiment and individual-based evolutionary simulations to gain mechanistic insights into contrasting evolutionary responses to a severe multi-year drought. We examined five traits in two populations of a native California plant, Clarkia xantiana, at three time points over 7 years. Earlier flowering phenology evolved in only one of the two populations, though both populations experienced similar drought severity and demographic declines and were estimated to have considerable additive genetic variance for flowering phenology. Pairing demographic and experimental data with evolutionary simulations suggested that while seed banks in both populations probably constrained evolutionary responses, a stronger seed bank in the non-evolving population resulted in evolutionary stasis. Gene flow through time via germ banks may be an important, underappreciated control on rapid evolution in response to extreme environmental perturbations.

Variance components of sex determination in the copepod Tigriopus californicus estimated from a pedigree analysis

Extensive theory exists regarding population sex ratio evolution that predicts equal sex ratio (when parental investment is equal). In most animals, sex chromosomes determine the sex of offspring, and this fixed genotype for sex has made theory difficult to test since genotypic variance for the trait (sex) is lacking. It has long been argued that the genotype has become fixed in most animals due to the strong selection for equal sex ratios. The marine copepod Tigriopus californicus has no sex chromosomes, multiple genes affecting female brood sex ratio, and a brood sex ratio that responds to selection. The species thus provides an opportune system in which to test established sex ratio theory. In this paper, we further our exploration of polygenic sex determination in T. californicus using an incomplete diallel crossing design for analysis of the variance components of sex determination in the species. Our data confirm the presence of extra-binomial variance for sex, further confirming that sex is not determined through simple Mendelian trait inheritance. In addition, our crosses and backcrosses of isofemale lines selected for biased brood sex ratios show intermediate phenotypic means, as expected if sex is a threshold trait determined by an underlying "liability" trait controlled by many genes of small effects. Furthermore, crosses between families from the same selection line had similar increases in phenotypic variance as crosses between families from different selection lines, suggesting families from artificial selection lines responded to selection pressure through different underlying genetic bases. Finally, we estimate heritability of an individual to be male or female on the observed binary scale as 0.09 (95% CI: 0.034-0.14). This work furthers our accumulating evidence for polygenic sex determination in T. californicus laying the foundation for this as a model species in future studies of sex ratio evolution theory.

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

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

A generalised approach to the study and understanding of adaptive evolution

Evolutionary theory has made large impacts on our understanding and management of the world, in part because it has been able to incorporate new data and new insights successfully. Nonetheless, there is currently a tension between certain biological phenomena and mainstream evolutionary theory. For example, how does the inheritance of molecular epigenetic changes fit into mainstream evolutionary theory? Is niche construction an evolutionary process? Is local adaptation via habitat choice also adaptive evolution? These examples suggest there is scope (and perhaps even a need) to broaden our views on evolution. We identify three aspects whose incorporation into a single framework would enable a more generalised approach to the understanding and study of adaptive evolution: (i) a broadened view of extended phenotypes; (ii) that traits can respond to each other; and (iii) that inheritance can be non-genetic. We use causal modelling to integrate these three aspects with established views on the variables and mechanisms that drive and allow for adaptive evolution. Our causal model identifies natural selection and non-genetic inheritance of adaptive parental responses as two complementary yet distinct and independent drivers of adaptive evolution. Both drivers are compatible with the Price equation; specifically, non-genetic inheritance of parental responses is captured by an often-neglected component of the Price equation. Our causal model is general and simplified, but can be adjusted flexibly in terms of variables and causal connections, depending on the research question and/or biological system. By revisiting the three examples given above, we show how to use it as a heuristic tool to clarify conceptual issues and to help design empirical research. In contrast to a gene-centric view defining evolution only in terms of genetic change, our generalised approach allows us to see evolution as a change in the whole causal structure, consisting not just of genetic but also of phenotypic and environmental variables.

Atlantic salmon (Salmo salar) age at maturity is strongly affected by temperature, population and age-at-maturity genotype

Age at maturity is a key life history trait involving a trade-off between survival risk and reproductive investment, and is an important factor for population structures. In ectotherms, a warming environment may have a dramatic influence on development and life history, but this influence may differ between populations. While an increasing number of studies have examined population-dependent reactions with temperature, few have investigated this in the context of maturation timing. Atlantic salmon, a species of high conservation relevance, is a good study species for this topic as it displays considerable variation in age at maturity, of which a large proportion has been associated with a genomic region including the strong candidate gene vgll3. Until now, the effect of this gene in the context of different environments and populations has not been studied. Using a large-scale common-garden experiment, we find strong effects of temperature, population-of-origin, and vgll3 genotype on maturation in 2-year-old male Atlantic salmon (Salmo salar). With a temperature difference of 1.8°C, maturation probability was 4.8 times higher in the warm treatment than the cold treatment. This temperature effect was population-specific and was higher in the southern (60.48°N) compared to the northern (65.01°N) population. The early maturation vgll3*E allele was associated with a significantly higher maturation probability, but there was no vgll3 interaction with temperature or population. Both body condition and body mass associated with maturation. The body mass association was only present in the warm treatment. Our findings demonstrate that (i) populations can vary in their response to temperature change in terms of age at maturity, (ii) high intrinsic growth could be associated with higher thermal sensitivity for life history variation and (iii) vgll3 effects on age at maturity might be similar between populations and different thermal environments.

Head Shape Heritability in the Hungarian Meadow Viper Vipera ursinii rakosiensis

Understanding heritability patterns in functionally relevant traits is a cornerstone for evaluating their evolutionary potential and their role in local adaptation. In this study, we investigated patterns of heritability in the head shape of the Hungarian meadow viper (Vipera ursinii rakosiensis). To this end, we used geometric morphometric data from 12 families composed of 8 mothers, 6 fathers and 221 offspring, bred in captivity at the Hungarian Meadow Viper Conservation Centre (Hungary). We separately evaluated maternal and paternal contributions to the offspring phenotype, in addition to additive genetic effects, all determined using a mixed animal model. Our results indicate a strong genetic and maternal contribution to head shape variations. In contrast, the paternal effects-which are rarely evaluated in wild-ranging species-as well as residual environmental variance, were minimal. Overall, our results indicate a high evolutionary potential for head shape in the Hungarian meadow viper, which suggests a strong contribution of this ecologically important trait in shaping the ability of this endangered species to adapt to changing conditions and/or habitats. Furthermore, our results suggest that maternal phenotypes should be carefully considered when designing captive breeding parental pairs for reinforcing the adaptive capacity of threatened populations, whereas the paternal phenotypes seem less relevant.

How longitudinal data can contribute to our understanding of host genetic effects on the gut microbiome

A key component of microbiome research is understanding the role of host genetic influence on gut microbial composition. However, it can be difficult to link host genetics with gut microbial composition because host genetic similarity and environmental similarity are often correlated. Longitudinal microbiome data can supplement our understanding of the relative role of genetic processes in the microbiome. These data can reveal environmentally contingent host genetic effects, both in terms of controlling for environmental differences and in comparing how genetic effects differ by environment. Here, we explore four research areas where longitudinal data could lend new insights into host genetic effects on the microbiome: microbial heritability, microbial plasticity, microbial stability, and host and microbiome population genetics. We conclude with a discussion of methodological considerations for future studies.

Genetic characteristics of colostrum refractive index and its use as a proxy for the concentration of immunoglobulins in Holstein cattle

BACKGROUND

Colostral concentration of immunoglobulins (Ig) is crucial for the passive transfer of antibodies from the cow to the new-born calf. Direct determination of Ig by the gold standard radial immunodiffusion method is demanding in terms of time and costs. For this reason, a refractometer is commonly used at the farm level for an indirect estimation of colostrum quality, which is given as the Ig concentration. In this study, colostrum samples were collected from 548 Italian Holstein cows within 6 h of calving. The refractive index (BRIX, %) of these samples was assessed using a portable optical refractometer, as well as the concentration of total protein, IgG, IgA, and IgM by radial immunodiffusion. A four-trait animal model was used to estimate genetic parameters for BRIX and the different immunoglobulin isotypes. A receiver operating characteristic analysis was carried out to evaluate the BRIX diagnostic accuracy.

RESULTS

Colostral BRIX was moderately heritable (0.26) and its genetic and phenotypic correlations with IgG (0.91, 0.78), IgA (0.57, 0.57), and IgM (0.71, 0.61) were all positive and of similar order, although the genetic correlations were generally higher than the phenotypic correlations. Low-quality colostrum samples, defined as those with an IgG concentration lower than 50 g/L, were accurately identified by the refractive index on the BRIX scale, with an area under the curve of 0.90.

CONCLUSIONS

The use of a refractometer is recommended on dairy farms to produce a proxy for colostral Ig concentration. BRIX is a useful phenotyping tool that can be used in cattle to improve the quality of colostrum for first feeding of calves through both traditional genetic and genomic strategies. Improving colostrum quality will reduce the incidence of failure of passive transfer of immunity in young stock.

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

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

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.

High heritability of telomere length and low heritability of telomere shortening in wild birds

Telomere length and telomere shortening predict survival in many organisms. This raises the question of the contribution of genetic and environmental effects to variation in these traits, which is still poorly known, particularly for telomere shortening. We used experimental (cross-fostering) and statistical (quantitative genetic "animal models") means to disentangle and estimate genetic and environmental contributions to telomere length variation in pedigreed free-living jackdaws (Corvus monedula). Telomere length was measured twice in nestlings, at ages 4 (n = 715) and 29 days (n = 474), using telomere restriction fragment (TRF) analysis, adapted to exclude interstitial telomeric sequences. Telomere length shortened significantly over the nestling period (10.4 ± 0.3 bp day^-1 ) and was highly phenotypically (r_P  = 0.95 ± 0.01) and genetically (r_G  > 0.99 ± 0.01) correlated within individuals. Additive genetic effects explained a major part of telomere length variation among individuals, with its heritability estimated at h²  = 0.74 on average. We note that TRF-based studies reported higher heritabilities than qPCR-based studies, and we discuss possible explanations. Parent-offspring regressions yielded similar heritability estimates for mothers and fathers when accounting for changes in paternal telomere length over life. Year effects explained a small but significant part of telomere length variation. Heritable variation for telomere shortening was low (h²  = 0.09 ± 0.11). The difference in heritability between telomere length (high) and telomere shortening (low) agrees with evolutionary theory, in that telomere shortening has stronger fitness consequences in this population. Despite the high heritability of telomere length, its evolvability, which scales the additive genetic variance by mean telomere length, was on average 0.48%. Hence, evolutionary change of telomere length due to selection is likely to be slow.

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

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

An exploration into the conversion of dominance to additive genetic variance in contrasting environments

PREMISE

The evolutionary response of a trait to environmental change depends upon the level of additive genetic variance. It has been long argued that sustained selection will tend to deplete additive genetic variance as favored alleles approach fixation. Non-additive genetic variance, due to interactions among alleles within and between loci, does not immediately contribute to an evolutionary response. However, shifts in the allele frequencies within and between interacting loci may convert non-additive variance into additive variance. Here we consider the possibility that an environmental shift may alter allelic interactions in ways that convert dominance into additive genetic variance.

METHODS

We grew a pedigreed population of Brassica rapa in greenhouse and field conditions. The field conditions mimicked agricultural conditions from which the base population was drawn, while the greenhouse featured benign conditions. We used Bayesian models to estimate the additive, dominance, and maternal components of quantitative genetic variance. We also estimated genetic correlations across environments using parental breeding values.

RESULTS

Although the additive genetic variance was elevated in the greenhouse condition, no consistent pattens emerged that would indicate a conversion of dominance variance. The unusually low genetic variance and broad confidence intervals for the variance estimates obtained through this analysis preclude definitive interpretations.

CONCLUSIONS

Further studies are needed to determine whether between-environment changes in additive genetic variance can be traced to conversion of dominance variance.

Susceptibility to a sexually transmitted disease in a wild koala population shows heritable genetic variance but no inbreeding depression

The koala, one of the most iconic Australian wildlife species, is facing several concomitant threats that are driving population declines. Some threats are well known and have clear methods of prevention (e.g., habitat loss can be reduced with stronger land-clearing control), whereas others are less easily addressed. One of the major current threats to koalas is chlamydial disease, which can have major impacts on individual survival and reproduction rates and can translate into population declines. Effective management strategies for the disease in the wild are currently lacking, and, to date, we know little about the determinants of individual susceptibility to disease. Here, we investigated the genetic basis of variation in susceptibility to chlamydia using one of the most intensively studied wild koala populations. We combined data from veterinary examinations, chlamydia testing, genetic sampling and movement monitoring. Out of our sample of 342 wild koalas, 60 were found to have chlamydia. Using genotype information on 5007 SNPs to investigate the role of genetic variation in determining disease status, we found no evidence of inbreeding depression, but a heritability of 0.11 (95% CI: 0.06-0.23) for the probability that koalas had chlamydia. Heritability of susceptibility to chlamydia could be relevant for future disease management, as it suggests adaptive potential for the population.

Evaluation of Genetic Diversity and Parasite-Mediated Selection of MHC Class I Genes in Emberiza godlewskii (Passeriformes: Emberizidae)

The major histocompatibility complex (MHC) is a multi-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next-generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles, clustering into nine supertypes. Compared with other passerine birds, the number of MHC class I alleles per individual in Godlewski’s buntings is high (mean 16.1 ± 3.3, median 16). In addition, we demonstrated signatures of historical and contemporary selection on MHC genes. Reflecting historical selection, ten amino acid sites in the antigen-binding domain showed signatures of balancing selection, eight of which exhibit high amino acid polymorphism. In terms of contemporary selection, we found that specific MHC supertypes were nominally associated with the infection of two malaria parasite lineages. These findings indicate the action of historical and possibly also contemporary balancing selection and suggest negative frequency-dependent or fluctuating selection as possible selection mechanisms.

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

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

Heritable and Nonheritable Rumen Bacteria Are Associated with Different Characters of Lactation Performance of Dairy Cows

Recent studies have reported that some rumen microbes are heritable. However, it is necessary to clarify the functions and specific contributions of the heritable rumen microbes to cattle phenotypes (microbiability) in comparison with those that are nonheritable. This study aimed to identify the distribution and predicted functions of heritable and nonheritable bacterial taxa at species level in the rumen of dairy cows and their respective contributions to energy-corrected milk yield, protein content and yield, and fat content and yield in milk. Thirty-two heritable and 674 nonheritable bacterial taxa were identified at species level, and the functional analysis revealed that predicted microbial functions for both groups were mainly enriched for energy, amino acid, and ribonucleotide metabolism. The mean microbiability (to reflect a single taxon's contribution) of heritable bacteria was found to range from 0.16% to 0.33% for the different milk traits, whereas the range for nonheritable bacteria was 0.03% to 0.06%. These findings suggest a strong contribution by host genetics in shaping the rumen microbiota, which contribute significantly to milk production traits. Therefore, there is an opportunity to further improve milk production traits through attention to host genetics and the interaction with the rumen microbiota. IMPORTANCE Rumen bacteria produce volatile fatty acids which exert a far-reaching influence on hepatic metabolism, mammary gland metabolism, and animal production. In the current study, 32 heritable and 674 nonheritable bacterial taxa at species level were identified, and shown to have different microbiability (overall community contribution) and mean microbiability (the average of a single taxon's contribution) for lactation performance. The predicted functions of heritable and nonheritable bacterial taxa also differed, suggesting that targeted nutritional and genetic breeding approaches could be used to manipulate them to improve dairy cow performance.

Heritability and Genetic Parameters for Semen Traits in Australian Sheep

Semen characteristics including volume, gross motility, spermatozoal concentration, and percent post-thaw motility are routinely assessed to determine the quality and quantity of an ejaculate prior to use in artificial breeding programs. Currently, artificial breeding programs in sheep place relatively little emphasis on ram-side factors, such as the fertilising potential of an ejaculate, which may contribute to variability in conception outcomes. Estimating genetic parameters for ejaculate quality and quantity traits could provide insights into whether selective breeding can be used to improve such ram-side traits, improving ovine reproductive performance and farm profitability. Therefore, in this study, a total of 11,470 ejaculate records, including data for ejaculate volume, gross motility, spermatozoal concentration, and percent post-thaw motility, collected over a 20-year period was used to estimate genetic parameters in sheep. Univariate and bivariate mixed model analysis was performed including a variety of fixed effects such as breed, age at collection, centre of collection, collection number, season of collection, and method of collection; and the permanent environmental effects associated with each ram, stud and year of collection, and the breeding value of rams included as random effects. The heritability for ejaculate volume, gross motility, concentration, and percent post-thaw motility was estimated to be 0.161, 0.170, 0.089, and 0.081. Repeatability estimates were moderate, ranging between 0.4126 and 0.5265. Overall, results indicate that semen traits are lowly heritable and moderately repeatable, indicating that these traits are significantly influenced by environmental variables.

Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution

Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.