A guide to using a multiple-matrix animal model to disentangle genetic and nongenetic causes of phenotypic variance

Effects of the Prenatal Maternal Health Status on Calf Disease Prevalences and Respective Genetic Parameter Estimates in German Holstein Cattle

The aim of the present study was to infer phenotypic responses and genetic parameters of the F1 calf diseases diarrhoea (DIAR) and pneumonia (PNEU) in dependency of the prenatal maternal health status (PMHS) of the dam and of the herd-calving year. The PMHS considered diagnoses for the cow disease mastitis (MAST) and claw disorders (CD) during gestation of F0 dams. Furthermore, 305-d milk production traits of F1 offspring from either healthy or diseased dam groups were compared. The study comprised 20,045 female calves (F1 = generation 1) and their corresponding dams (F0 = parental generation 0), kept in 41 large-scale herds. All F1 calves were from their dams' 2nd parity, implying that all dam (maternal) diseases were recorded during the first lactation and dry period of the dams. The F1 calves were phenotyped for DIAR up to 30 days post-partum, and for PNEU up to 180 days of age. At least one entry for the respective disease implied a score = 1 = sick, otherwise, a score = 0 = healthy, was assigned. Production records of the 10,129 F1 cows comprised 305-d records in first lactation for milk yield (MY), protein yield (PY) and fat yield (FY). Linear and generalised linear mixed models were applied to infer phenotypic responses of F1 traits in dependency of the PMHS for CD and MAST. A diagnosis for MAST or CD in F0 cows during gestation was significantly (p ≤ 0.05) associated with an increased prevalence for DIAR and PNEU, with pairwise differences of least-squares-means between calves from healthy and diseased cow groups up to 3.61%. The effects of PMHS on 305-d production traits in offspring were non-significant (p > 0.05). In bivariate genetic analyses, DIAR and PNEU were defined as different traits according to the PMHS, i.e., DIAR-MASThealthy and DIAR-MASTdiseased, DIAR-CDhealthy and DIAR-CDdiseased, PNEU-MASThealthy and PNEU-MASTdiseased, and PNEU-CDhealthy and PNEU-CDdiseased. The direct heritabilities for DIAR and PNEU were quite similar in the healthy and respective diseased dam group. Slightly larger direct heritabilities in the diseased dam groups were due to increased genetic variances. Maternal heritabilities were quite stable and smaller than the direct heritabilities. In random regression models, genetic parameters for DIAR and PNEU were estimated along the continuous herd-calving-year prevalence scale, considering a prevalence for MAST and CD (based on the 20,045 dam records plus 16,193 herd contemporary records) in the range from 0% to 30%. Direct heritabilities for PNEU were quite stable along the herd-calving-year gradient for MAST and CD. For DIAR, we observed stronger estimate fluctuations, especially increasing direct heritabilities in dependency of the herd-calving-year prevalence for MAST from 0.13 (at a MAST prevalence of 0%) to 0.30 (at a MAST prevalence of 30%). Consequently, obvious genotype x herd-calving-year PMHS interactions were observed for DIAR on the prenatal MAST scale, with a minimal correlation of 0.48 between direct genetic effects at 0% MAST prevalence and at 30% MAST prevalence. The correlations between direct genetic and maternal genetic effects were antagonistic at all herd-calving-year prevalence levels, displaying strongest fluctuations for "DIAR-MAST." The genotype x herd-calving-year PMHS interactions for DIAR suggest consideration of specific sires according to the herd health status for CD and for MAST.

Shared environmental similarity between relatives influences heritability of reproductive timing in wild great tits

Intraspecific variation is necessary for evolutionary change and population resilience, but the extent to which it contributes to either depends on the causes of this variation. Understanding the causes of individual variation in traits involved with reproductive timing is important in the face of environmental change, especially in systems where reproduction must coincide with seasonal resource availability. However, separating the genetic and environmental causes of variation is not straightforward, and there has been limited consideration of how small-scale environmental effects might lead to similarity between individuals that occupy similar environments, potentially biasing estimates of genetic heritability. In ecological systems, environments are often complex in spatial structure, and it may therefore be important to account for similarities in the environments experienced by individuals within a population beyond considering spatial distances alone. Here, we construct multi-matrix quantitative genetic animal models using over 11,000 breeding records (spanning 35 generations) of individually-marked great tits (Parus major) and information about breeding proximity and habitat characteristics to quantify the drivers of variability in two key seasonal reproductive timing traits. We show that the environment experienced by related individuals explains around a fifth of the variation seen in reproductive timing, and accounting for this leads to decreased estimates of heritability. Our results thus demonstrate that environmental sharing between relatives can strongly affect estimates of heritability and therefore alter our expectations of the evolutionary response to selection.

Genetic, natal and spatial drivers of social phenotypes in wild great tits

In social animals, group dynamics profoundly influence collective behaviours, vital in processes like information sharing and predator vigilance. Disentangling the causes of individual-level variation in social behaviours is crucial for understanding the evolution of sociality. This requires the estimation of the genetic and environmental basis of these behaviours, which is challenging in uncontrolled wild populations. In this study, we partitioned genetic, developmental and spatial environmental variation in repeatable social network traits derived from foraging events. We used a multi-generational pedigree and social data for 1823 individuals with over 800,000 observations from a long-term monitored great tit population. Animal models indicated minimal narrow-sense heritability (2%-3%) in group size choice, further reduced when the spatial location was considered, which itself explains a substantial 30% of the observed variation. Individual gregariousness also had a small genetic component, with a low heritability estimate for degree (<5%). Centrality showed heritability up to 10% in one of 3 years sampled, whereas betweenness showed none, indicating modest genetic variation in individual sociability, but not group-switching tendencies. These findings suggest a small, albeit detectable, genetic basis for individual sociality but pronounced spatial effects. Furthermore, our study highlights the importance of common environment effects (natal origin and brood identity), which essentially negated genetic effects when explicitly accounted for. In addition, we demonstrate that phenotypic resemblance can be a result of similarities beyond shared genes; spatial proximity at birth and natal environmental similarity explained up to 8% of the variation in individual sociability. Our results thus emphasise the role of non-genetic factors, particularly developmental and spatial variation, in shaping individual social behavioural tendencies.

A Recursive Model Approach to Include Epigenetic Effects in Genetic Evaluations Using Simulated DNA Methylation Effects

The advancement of epigenetics has highlighted DNA methylation as an intermediate-omic influencing gene regulation and phenotypic expression. With emerging technologies enabling the large-scale and affordable capture of methylation data, there is growing interest in integrating this information into genetic evaluation models for animal breeding. This study used methylome information from six dairy cows to simulate the methylation profile of 13,183 genotyped animals. The liability to methylation was treated as an additive trait, while a trait moderated by methylation effects was also simulated. A multiomic model (GOBLUP) was adapted to incorporate methylation data in genomic and genetic evaluations, using the traditional BLUP method as a benchmark. The GOBLUP accurately recovered heritability estimates for the liability to methylation in all low, medium and high heritability scenarios and was consistent at estimating the heritability for the epigenetics-moderated trait of interest at a low-medium heritability of 0.14. The genetic variance recovered by the BLUP model was influenced by the h2 of the liability to methylation, and a part of the methylation variance for the phenotypic trait was captured as additive. The h2 of the phenotypic trait partially relies on the h2 value for the methylation windows in the traditional model. A newly proposed estimated epigenetic value (EEV) combines the traditional additive genetic information from genotyping arrays with epigenetic information. The correlation between the traditional estimated breeding value (EBV) and EEV was high (0.92-0.99 depending on the scenario), but the correlation of the EEV with the true breeding value was higher than the correlation between the traditional EBV and the TBV (0.85 vs. 0.75, 0.71 vs. 0.66 and 0.61 vs. 0.62 depending on the scenario). This study demonstrates that the GOBLUP multiomic recursive model can effectively separates additive and epigenetic variances, enabling improved breeding decisions by accounting for genetic liability to DNA methylation. This enables more informed breeding decisions, optimising selection for desired traits. Emerging sequencing techniques offer new opportunities for cost-effective simultaneous acquisition of genetic and epigenetic data, further enhancing breeding accuracy.

Choose your mother wisely: the familial resemblance of bone adaptation

This study demonstrates how complex bone microarchitectural features can be summarized to describe bone adaptations seen with aging in women, which are consistent with the stages of osteoporosis. Additionally, we showed familial resemblance in these bone microarchitectural traits between mothers and daughters that can be used to predict bone adaptations.

INTRODUCTION

Patient-specific characterization of bone quality can reduce complex microarchitectural features to common combinations of bone characteristics, known as bone phenotypes. This study investigated whether there is a general trend in bone phenotype change over time seen with aging in females and whether there is a familial resemblance to phenotype membership between mothers and daughters.

METHODS

Bone phenotype membership was calculated on biological mother and daughter pairs (Participants = 101), scanned using high resolution peripheral quantitative computed tomography, to the three pre-defined phenotypes (healthy, low volume, and low density). The trajectory of bone phenotype with age was explored using all participant's data. Linear regression models were used to assess the familial resemblance of phenotyping in the mother-daughter pairs.

RESULTS

When stratified for age, the trajectory of the phenotype membership transitioned from healthy (20-40 years), to low volume (40-60 years), to low density (60-80 years), which similarly aligns with the stages of osteoporosis observed in females. Familial resemblance (½h2) was observed in the healthy phenotype (β = 0.432, p < 0.01). Predictive modelling showed a significant association in phenotype membership between mothers and daughters in the healthy (R2 = 0.347, p = 0.04) and low volume (R2 = 0.416, p < 0.01) phenotypes, adjusted for age, height, and weight.

CONCLUSION

Our results suggest that phenotype membership in females changes with age in a pattern that is consistent with the stages of osteoporosis. Additionally, we showed familial resemblance in bone phenotype, which can be used to predict bone adaptations between mothers and daughters that are associated with bone loss with aging.

Ontogeny shapes individual dietary specialization in female European brown bears (Ursus arctos)

Individual dietary specialization, where individuals occupy a subset of a population's wider dietary niche, is a key factor determining a species resilience against environmental change. However, the ontogeny of individual specialization, as well as associated underlying social learning, genetic, and environmental drivers, remain poorly understood. Using a multigenerational dataset of female European brown bears (Ursus arctos) followed since birth, we discerned the relative contributions of environmental similarity, genetic heritability, maternal effects, and offspring social learning from the mother to individual specialization. Individual specialization accounted for 43% of phenotypic variation and spanned half a trophic position, with individual diets ranging from omnivorous to carnivorous. The main determinants of dietary specialization were social learning during rearing (13%), environmental similarity (5%), maternal effects (11%), and permanent between-individual effects (9%), whereas the contribution of genetic heritability (3%) was negligible. The trophic position of offspring closely resembled the trophic position of their mothers during the first 3-4 years of independence, but waned with increasing time since separation. Our study shows that social learning and maternal effects were more important for individual dietary specialization than environmental composition. We propose a tighter integration of social effects into studies of range expansion and habitat selection under global change.

Family resemblance in color-patch size is not affected by stress experience in a cichlid fish

Animal body coloration is often linked to social dominance and mating success. This is because it can carry information on an animal's body condition and competitive ability by reflecting the genetic quality of individuals or by responding to their current or past living conditions. The present study investigates genetic and environmental effects on a conspicuous color pattern of the cichlid fish Tropheus sp. black "Ikola," in which the size of a carotenoid-based yellow area on the body co-varies with social dominance. To examine environmental plasticity of the color pattern, we tested for effects of early-life stress, induced by reduced feeding of juveniles prior to color pattern formation, as well as effects of a stress treatment administered to fully colored adult fish. None of the stress treatments affected the color pattern as quantified by the width of the yellow bar. However, offspring bar width was correlated to parental values in mid-parent-mid-offspring regression analyses, and animal models estimated significant additive genetic effects on bar width, indicating heritability of the trait. Depending on the random effects structure of the animal models (i.e., whether including or excluding maternal and brood effects), narrow-sense heritability estimates for bar width ranged between 0.2 and 0.8, with the strongest statistical support for the highest estimate. In each of the alternative models, a large proportion of the total variance in bar width was explained by the included random effects, suggesting that bar width is strongly determined by genetic factors or shared maternal and brood environments, with limited scope for environmental influences later in life.

Perspective: the evolution of hormones and person perception-a quantitative genetic framework

Evolutionary biology provides a unifying theory for testing hypotheses about the relationship between hormones and person perception. Person perception usually receives attention from the perspective of sexual selection. However, because person perception is one trait in a suite regulated by hormones, univariate approaches are insufficient. In this Perspectives article, quantitative genetics is presented as an important but underutilized framework for testing evolutionary hypotheses within this literature. We note tacit assumptions within the current literature on psychiatric genetics, which imperil the interpretation of findings thus far. As regulators of a diverse manifold of traits, hormones mediate tradeoffs among an array of functions. Hormonal pleiotropy also provides the basis of correlational selection, a process whereby selection on one trait in a hormone-mediated suite generates selection on the others. This architecture provides the basis for conflicts between sexual and natural selection within hormone-mediated suites. Due to its role in person perception, psychiatric disorders, and reproductive physiology, the sex hormone estrogen is highlighted as an exemplar here. The implications of this framework for the evolution of person perception are discussed. Empirical quantification of selection on traits within hormone-mediated suites remains an important gap in this literature with great potential to illuminate the fundamental nature of psychiatric disorders.

Vocal convergence and social proximity shape the calls of the most basal Passeriformes, New Zealand Wrens

Despite extensive research on avian vocal learning, we still lack a general understanding of how and when this ability evolved in birds. As the closest living relatives of the earliest Passeriformes, the New Zealand wrens (Acanthisitti) hold a key phylogenetic position for furthering our understanding of the evolution of vocal learning because they share a common ancestor with two vocal learners: oscines and parrots. However, the vocal learning abilities of New Zealand wrens remain unexplored. Here, we test for the presence of prerequisite behaviors for vocal learning in one of the two extant species of New Zealand wrens, the rifleman (Acanthisitta chloris). We detect the presence of unique individual vocal signatures and show how these signatures are shaped by social proximity, as demonstrated by group vocal signatures and strong acoustic similarities among distantly related individuals in close social proximity. Further, we reveal that rifleman calls share similar phenotypic variance ratios to those previously reported in the learned vocalizations of the zebra finch, Taeniopygia guttata. Together these findings provide strong evidence that riflemen vocally converge, and though the mechanism still remains to be determined, they may also suggest that this vocal convergence is the result of rudimentary vocal learning abilities.

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.

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.

Early-life behavior, survival, and maternal personality in a wild marsupial

Individual behavior varies for many reasons, but how early in life are such differences apparent, and are they under selection? We investigated variation in early-life behavior in a wild eastern gray kangaroo (Macropus giganteus) population, and quantified associations of behavior with early survival. Behavior of young was measured while still in the pouch and as subadults, and survival to weaning was monitored. We found consistent variation between offspring of different mothers in levels of activity at the pouch stage, in flight initiation distance (FID) as subadults, and in subadult survival, indicating similarity between siblings. There was no evidence of covariance between the measures of behavior at the pouch young versus subadult stages, nor of covariance of the early-life behavioral traits with subadult survival. However, there was a strong covariance between FIDs of mothers and those of their offspring tested at different times. Further, of the total repeatability of subadult FID (51.5%), more than half could be attributed to differences between offspring of different mothers. Our results indicate that 1) behavioral variation is apparent at a very early stage of development (still in the pouch in the case of this marsupial); 2) between-mother differences can explain much of the repeatability (or "personality") of juvenile behavior; and 3) mothers and offspring exhibit similar behavioral responses to stimuli. However, 4) we found no evidence of selection via covariance between early-life or maternal behavioral traits and juvenile survival in this wild marsupial.

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.

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.

Position of Hungarian Merino among other Merinos, within-breed genetic similarity network and markers associated with daily weight gain

OBJECTIVE

In this study, we aimed to position the Hungarian Merino among other Merino-derived sheep breeds, explore the characteristics of our sampled animals' genetic similarity network within the breed, and highlight single nucleotide polymorphisms (SNPs) associated with daily weight-gain.

METHODS

Hungarian Merino (n = 138) was genotyped on Ovine SNP50 Bead Chip (Illumina, San Diego, CA, USA) and positioned among 30 Merino and Merino-derived breeds (n = 555). Population characteristics were obtained via PLINK, SVS, Admixture, and Treemix software, within-breed network was analysed with python networkx 2.3 library. Daily weight gain of Hungarian Merino was standardised to 60 days and was collected from the database of the Association of Hungarian Sheep and Goat Breeders. For the identification of loci associated with daily weight gain, a multi-locus mixed-model was used.

RESULTS

Supporting the breed's written history, the closest breeds to Hungarian Merino were Estremadura and Rambouillet (pairwise FST values are 0.035 and 0.036, respectively). Among Hungarian Merino, a highly centralised connectedness has been revealed by network analysis of pairwise values of identity-by-state, where the animal in the central node had a betweenness centrality value equal to 0.936. Probing of daily weight gain against the SNP data of Hungarian Merinos revealed five associated loci. Two of them, OAR8_17854216.1 and s42441.1 on chromosome 8 and 9 (-log10P>22, false discovery rate<5.5e-20) and one locus on chromosome 20, s28948.1 (-log10P = 13.46, false discovery rate = 4.1e-11), were close to the markers reported in other breeds concerning daily weight gain, six-month weight, and post-weaning gain.

CONCLUSION

The position of Hungarian Merino among other Merino breeds has been determined. We have described the similarity network of the individuals to be applied in breeding practices and highlighted several markers useful for elevating the daily weight gain of Hungarian Merino.

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.

Social animal models for quantifying plasticity, assortment, and selection on interacting phenotypes

Both assortment and plasticity can facilitate social evolution, as each may generate heritable associations between the phenotypes and fitness of individuals and their social partners. However, it currently remains difficult to empirically disentangle these distinct mechanisms in the wild, particularly for complex and environmentally responsive phenotypes subject to measurement error. To address this challenge, we extend the widely used animal model to facilitate unbiased estimation of plasticity, assortment and selection on social traits, for both phenotypic and quantitative genetic (QG) analysis. Our social animal models (SAMs) estimate key evolutionary parameters for the latent reaction norms underlying repeatable patterns of phenotypic interaction across social environments. As a consequence of this approach, SAMs avoid inferential biases caused by various forms of measurement error in the raw phenotypic associations between social partners. We conducted a simulation study to demonstrate the application of SAMs and investigate their performance for both phenotypic and QG analyses. With sufficient repeated measurements, we found desirably high power, low bias and low uncertainty across model parameters using modest sample and effect sizes, leading to robust predictions of selection and adaptation. Our results suggest that SAMs will readily enhance social evolutionary research on a variety of phenotypes in the wild. We provide detailed coding tutorials and worked examples for implementing SAMs in the Stan statistical programming language.

Permutation tests for hypothesis testing with animal social network data: Problems and potential solutions

Permutation tests are widely used to test null hypotheses with animal social network data, but suffer from high rates of type I and II error when the permutations do not properly simulate the intended null hypothesis.Two common types of permutations each have limitations. Pre-network (or datastream) permutations can be used to control 'nuisance effects' like spatial, temporal or sampling biases, but only when the null hypothesis assumes random social structure. Node (or node-label) permutation tests can test null hypotheses that include nonrandom social structure, but only when nuisance effects do not shape the observed network.We demonstrate one possible solution addressing these limitations: using pre-network permutations to adjust the values for each node or edge before conducting a node permutation test. We conduct a range of simulations to estimate error rates caused by confounding effects of social or non-social structure in the raw data.Regressions on simulated datasets suggest that this 'double permutation' approach is less likely to produce elevated error rates relative to using only node permutations, pre-network permutations or node permutations with simple covariates, which all exhibit elevated type I errors under at least one set of simulated conditions. For example, in scenarios where type I error rates from pre-network permutation tests exceed 30%, the error rates from double permutation remain at 5%.The double permutation procedure provides one potential solution to issues arising from elevated type I and type II error rates when testing null hypotheses with social network data. We also discuss alternative approaches that can provide robust inference, including fitting mixed effects models, restricted node permutations, testing multiple null hypotheses and splitting large datasets to generate replicated networks. Finally, we highlight ways that uncertainty can be explicitly considered and carried through the analysis.

Gut microbiome heritability is nearly universal but environmentally contingent

Relatives have more similar gut microbiomes than nonrelatives, but the degree to which this similarity results from shared genotypes versus shared environments has been controversial. Here, we leveraged 16,234 gut microbiome profiles, collected over 14 years from 585 wild baboons, to reveal that host genetic effects on the gut microbiome are nearly universal. Controlling for diet, age, and socioecological variation, 97% of microbiome phenotypes were significantly heritable, including several reported as heritable in humans. Heritability was typically low (mean = 0.068) but was systematically greater in the dry season, with low diet diversity, and in older hosts. We show that longitudinal profiles and large sample sizes are crucial to quantifying microbiome heritability, and indicate scope for selection on microbiome characteristics as a host phenotype.

Social competition as a driver of phenotype-environment correlations: implications for ecology and evolution

While it is universally recognised that environmental factors can cause phenotypic trait variation via phenotypic plasticity, the extent to which causal processes operate in the reverse direction has received less consideration. In fact individuals are often active agents in determining the environments, and hence the selective regimes, they experience. There are several important mechanisms by which this can occur, including habitat selection and niche construction, that are expected to result in phenotype–environment correlations (i.e. non‐random assortment of phenotypes across heterogeneous environments). Here we highlight an additional mechanism – intraspecific competition for preferred environments – that may be widespread, and has implications for phenotypic evolution that are currently underappreciated. Under this mechanism, variation among individuals in traits determining their competitive ability leads to phenotype–environment correlation; more competitive phenotypes are able to acquire better patches. Based on a concise review of the empirical evidence we argue that competition‐induced phenotype–environment correlations are likely to be common in natural populations before highlighting the major implications of this for studies of natural selection and microevolution. We focus particularly on two central issues. First, competition‐induced phenotype–environment correlation leads to the expectation that positive feedback loops will amplify phenotypic and fitness variation among competing individuals. As a result of being able to acquire a better environment, winners gain more resources and even better phenotypes – at the expense of losers. The distinction between individual quality and environmental quality that is commonly made by researchers in evolutionary ecology thus becomes untenable. Second, if differences among individuals in competitive ability are underpinned by heritable traits, competition results in both genotype–environment correlations and an expectation of indirect genetic effects (IGEs) on resource‐dependent life‐history traits. Theory tells us that these IGEs will act as (partial) constraints, reducing the amount of genetic variance available to facilitate evolutionary adaptation. Failure to recognise this will lead to systematic overestimation of the adaptive potential of populations. To understand the importance of these issues for ecological and evolutionary processes in natural populations we therefore need to identify and quantify competition‐induced phenotype–environment correlations in our study systems. We conclude that both fundamental and applied research will benefit from an improved understanding of when and how social competition causes non‐random distribution of phenotypes, and genotypes, across heterogeneous environments.

Indirect genetic and environmental effects on behaviors, morphology, and life-history traits in a wild Eastern chipmunk population

Additive genetic variance in a trait reflects its potential to respond to selection, which is key for adaptive evolution in the wild. Social interactions contribute to this genetic variation through indirect genetic effects-the effect of an individual's genotype on the expression of a trait in a conspecific. However, our understanding of the evolutionary importance of indirect genetic effects in the wild and of their strength relative to direct genetic effects is limited. In this study, we assessed how indirect genetic effects contribute to genetic variation of behavioral, morphological, and life-history traits in a wild Eastern chipmunk population. We also compared the contribution of direct and indirect genetic effects to traits evolvabilities and related these effects to selection strength across traits. We implemented a novel approach integrating the spatial structure of social interactions in quantitative genetic analyses, and supported the reliability of our results with power analyses. We found indirect genetic effects for trappability and relative fecundity, little direct genetic effects in all traits and a large role for direct and indirect permanent environmental effects. Our study highlights the potential evolutionary role of social permanent environmental effects in shaping phenotypes of conspecifics through adaptive phenotypic plasticity.

Genome-Wide Variation in DNA Methylation Predicts Variation in Leaf Traits in an Ecosystem-Foundational Oak Species

Epigenetic modifications such as DNA methylation are a potential mechanism for trees to respond to changing environments. However, it remains controversial the extent to which DNA methylation impacts ecologically important traits that influence fitness. In this study, we used reduced-representation bisulfite sequencing to associate genomic and epigenomic variation with seven phenotypic traits related to growth, leaf function, and disease susceptibility in 160 valley oak (Quercus lobata) saplings planted across two common gardens in California. We found that DNA methylation was associated with a significant fraction of phenotypic variance in plant height, leaf lobedness, powdery mildew infection, and trichome density. Two of the seven traits were significantly associated with DNA methylation in the CG context, three traits were significantly associated with CHG methylation, and two traits were significantly associated with CHH methylation. Notably, controlling for genomic variation in SNPs generally reduced the amount of trait variation explained by DNA methylation. Our results suggest that DNA methylation may serve as a useful biomarker to predict phenotypic variation in trees, though it remains unclear the degree to which DNA methylation is a causal mechanism driving phenotypic variation in forest tree species.

Multiple spatial behaviours govern social network positions in a wild ungulate

The structure of wild animal social systems depends on a complex combination of intrinsic and extrinsic drivers. Population structuring and spatial behaviour are key determinants of individuals' observed social behaviour, but quantifying these spatial components alongside multiple other drivers remains difficult due to data scarcity and analytical complexity. We used a 43-year dataset detailing a wild red deer population to investigate how individuals' spatial behaviours drive social network positioning, while simultaneously assessing other potential contributing factors. Using Integrated Nested Laplace Approximation (INLA) multi-matrix animal models, we demonstrate that social network positions are shaped by two-dimensional landscape locations, pairwise space sharing, individual range size, and spatial and temporal variation in population density, alongside smaller but detectable impacts of a selection of individual-level phenotypic traits. These results indicate strong, multifaceted spatiotemporal structuring in this society, emphasising the importance of considering multiple spatial components when investigating the causes and consequences of sociality.

Phenotypic Response to Light Versus Shade Associated with DNA Methylation Changes in Snapdragon Plants (Antirrhinum majus)

The phenotypic plasticity of plants in response to change in their light environment, and in particularly, to shade is a schoolbook example of ecologically relevant phenotypic plasticity with evolutionary adaptive implications. Epigenetic variation is known to potentially underlie plant phenotypic plasticity. Yet, little is known about its role in ecologically and evolutionary relevant mechanisms shaping the diversity of plant populations in nature. Here we used a reference-free reduced representation bisulfite sequencing method for non-model organisms (epiGBS) to investigate changes in DNA methylation patterns across the genome in snapdragon plants (Antirrhinum majus L.). We exposed plants to sunlight versus artificially induced shade in four highly inbred lines to exclude genetic confounding effects. Our results showed that phenotypic plasticity in response to light versus shade shaped vegetative traits. They also showed that DNA methylation patterns were modified under light versus shade, with a trend towards global effects over the genome but with large effects found on a restricted portion. We also detected the existence of a correlation between phenotypic and epigenetic variation that neither supported nor rejected its potential role in plasticity. While our findings imply epigenetic changes in response to light versus shade environments in snapdragon plants, whether these changes are directly involved in the phenotypic plastic response of plants remains to be investigated. Our approach contributed to this new finding but illustrates the limits in terms of sample size and statistical power of population epigenetic approaches in non-model organisms. Pushing this boundary will be necessary before the relationship between environmentally induced epigenetic changes and phenotypic plasticity is clarified for ecologically relevant mechanisms with evolutionary implications.

Different perspectives on non-genetic inheritance illustrate the versatile utility of the Price equation in evolutionary biology

The diversity of genetic and non-genetic processes that make offspring resemble their parents are increasingly well understood. In addition to genetic inheritance, parent-offspring similarity is affected by epigenetic, behavioural and cultural mechanisms that collectively can be referred to as non-genetic inheritance. Given the generality of the Price equation as a description of evolutionary change, is it not surprising that the Price equation has been adopted to model the evolutionary implications of non-genetic inheritance. In this paper, we briefly introduce the heredity perspectives on which those models rely, discuss the extent to which these perspectives make different assumptions and place different emphases on the roles of heredity and development in evolution, and the types of empirical research programmes they motivate. The existence of multiple perspectives and explanatory aims highlight, on the one hand, the versatility of the Price equation and, on the other hand, the importance of understanding how heredity and development can be conceptualized in evolutionary studies. This article is part of the theme issue 'Fifty years of the Price equation'.

Pedigree-free quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer

Assessing the evolutionary potential of animal populations in the wild is crucial to understanding how they may respond to selection mediated by rapid environmental change (e.g. habitat loss and fragmentation). A growing number of studies have investigated the adaptive role of behaviour, but assessments of its genetic basis in a natural setting remain scarce. We combined intensive biologging technology with genome-wide data and a pedigree-free quantitative genetic approach to quantify repeatability, heritability and evolvability for a suite of behaviours related to the risk avoidance-resource acquisition trade-off in a wild roe deer (Capreolus capreolus) population inhabiting a heterogeneous, human-dominated landscape. These traits, linked to the stress response, movement and space-use behaviour, were all moderately to highly repeatable. Furthermore, the repeatable among-individual component of variation in these traits was partly due to additive genetic variance, with heritability estimates ranging from 0.21 ± 0.08 to 0.70 ± 0.11 and evolvability ranging from 1.1% to 4.3%. Changes in the trait mean can therefore occur under hypothetical directional selection over just a few generations. To the best of our knowledge, this is the first empirical demonstration of additive genetic variation in space-use behaviour in a free-ranging population based on genomic relatedness data. We conclude that wild animal populations may have the potential to adjust their spatial behaviour to human-driven environmental modifications through microevolutionary change.

The role of common ancestry and gene flow in the evolution of human-directed play behaviour in dogs

Among-population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human-directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among-breed variance of human-directed play behaviour while also accounting for within-breed variance. We found that 75% of the among-breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human-directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within-species/among-population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among-breed variation of a behavioural phenotype.

Intergenerational Transmission of Characters Through Genetics, Epigenetics, Microbiota, and Learning in Livestock

Evolutionary biologists studying wild species have demonstrated that genetic and non-genetic sources of information are inherited across generations and are therefore responsible for phenotypic resemblance between relatives. Although it has been postulated that non-genetic sources of inheritance are important in natural selection, they are not taken into account for livestock selection that is based on genetic inheritance only. According to the natural selection theory, the contribution of non-genetic inheritance may be significant for the transmission of characters. If this theory is confirmed in livestock, not considering non-genetic means of transmission in selection schemes might prevent achieving maximum progress in the livestock populations being selected. The present discussion paper reviews the different mechanisms of genetic and non-genetic inheritance reported in the literature as occurring in livestock species. Non-genetic sources of inheritance comprise information transmitted via physical means, such as epigenetic and microbiota inheritance, and those transmitted via learning mechanisms: behavioral, cultural and ecological inheritance. In the first part of this paper we review the evidence that suggests that both genetic and non-genetic information contribute to inheritance in livestock (i.e. transmitted from one generation to the next and causing phenotypic differences between individuals) and discuss how the environment may influence non-genetic inherited factors. Then, in a second step, we consider methods for favoring the transmission of non-genetic inherited factors by estimating and selecting animals on their extended transmissible value and/or introducing favorable non-genetic factors via the animals’ environment.

Anthropometric heritability and child growth in a Caribbean village: A quantitative genetic analysis of longitudinal height, weight, and body mass index in Bwa Mawego, Dominica

OBJECTIVES

Body size and composition vary widely among individuals and populations, and long-term research in diverse contexts informs our understanding of genetic, cultural, and environmental impacts on this variation. We analyze longitudinal measures of height, weight, and body mass index (BMI) from a Caribbean village, estimating the extent to which these anthropometrics are shaped by genetic variance in a small-scale population of mixed ancestry.

MATERIALS AND METHODS

Longitudinal data from a traditionally horticultural village in Dominica document height and weight in a non-Western population that is transitioning to increasingly Westernized lifestyles, and an 11-generation pedigree enables us to estimate the proportions of phenotypic variation in height, weight, and BMI attributed to genetic variation. We assess within-individual variation across growth curves as well as heritabilities of these traits for 260 individuals using Bayesian variance component estimation.

RESULTS

Age, sex, and secular trends account for the majority of anthropometric variation in these longitudinal data. Independent of age, sex, and secular trends, our analyses show high repeatabilities for the remaining variation in height, weight, and BMI growth curves (>0.75), and moderate heritabilities (h² _height = 0.68, h² _weight = 0.64, h² _BMI = 0.49) reveal clear genetic signals that account for large proportions of the variation in body size observed between families. Secular trends show increases of 6.5% in height and 16.0% in weight from 1997 to 2017.

DISCUSSION

This horticultural Caribbean population has transitioned to include more Westernized foods and technologies over the decades captured in this analysis. BMI varies widely between individuals and is significantly shaped by genetic variation, warranting future exploration with other physiological correlates and associated genetic variants.

Adaptive Alignment of Plasticity With Genetic Variation and Selection

Theoretical research has outlined how selection may shape both genetic variation and the expression of phenotypic plasticity in multivariate trait space. Specifically, research regarding the evolution of patterns of additive genetic variances and covariances (summarized in matrix form as G) and complementary research into how selection may shape adaptive plasticity lead to the general prediction that G, plasticity, and selection surfaces are all expected to align with each other. However, less well discussed is how this prediction might be assessed and how the modeled theoretical processes are expected to manifest in actual populations. Here, we discuss the theoretical foundations of the overarching prediction of alignment, what alignment mathematically means, how researchers might test for alignment and important caveats to this testing. The most important caveat concerns the fact that, for plasticity, the prediction of alignment only applies to cases where plasticity is adaptive, whereas organisms express considerable plasticity that may be neutral or even maladaptive. We detail the ramifications of these alternative expressions of plasticity vis-à-vis predictions of alignment. Finally, we briefly highlight some important interpretations of deviations from the prediction of alignment and what alignment might mean for populations experiencing environmental and selective changes.