The quantitative genetics of indirect genetic effects: a selective review of modelling issues

Natural selection could determine whether Acropora corals persist under expected climate change

Marine heatwaves are intensifying under climate change, exposing populations of reef-building corals to mass mortality and intense selective pressure. It remains unknown whether adaptation can keep pace with warming and maintain reef functioning. We have developed an eco-evolutionary metapopulation model for Acropora, an ecologically important yet thermally sensitive coral taxon. We found that, although corals have some adaptation capacity, they will suffer severe heatwave-induced declines over the coming decades. For a future in which emissions lead to ~3°C of global warming, natural selection could allow populations to persist, albeit in severely depleted states with elevated extinction risk and potential loss of ecosystem functioning. Yet, for thermally sensitive coral populations to thrive beyond 2050, there must be rapid reductions of greenhouse gas emissions that limit global warming to 2°C.

Review: Opportunities and challenges for the genetic selection of dairy calf disease traits

Interest in dairy cow health continues to grow as we better understand health's relationship with production potential and animal welfare. Over the past decade, efforts have been made to incorporate health traits into national genetic evaluations. However, they have focused on the mature cow, with calf health largely being neglected. Diarrhoea and respiratory disease comprise the main illnesses with regard to calf health. Conventional methods to control calf disease involve early separation of calves from the dam and housing calves individually. However, public concern regarding these methods, and growing evidence that these methods may negatively impact calf development, mean the dairy industry may move away from these practices. Genetic selection may be a promising tool to address these major disease issues. In this review, we examined current literature for enhancing calf health through genetics and discussed alternative approaches to improve calf health via the use of epidemiological modelling approaches, and the potential of indirectly selecting for improved calf health through improving colostrum quality. Heritability estimates on the observed scale for diarrhoea ranged from 0.03 to 0.20, while for respiratory disease, estimates ranged from 0.02 to 0.24. The breadth in these ranges is due, at least in part, to differences in disease prevalence, population structure, data editing and models, as well as data collection practices, which should be all considered when comparing literature values. Incorporation of epidemiological theory into quantitative genetics provides an opportunity to better determine the level of genetic variation in disease traits, as it accounts for disease transmission among contemporaries. Colostrum intake is a major determinant of whether a calf develops either respiratory disease or diarrhoea. Colostrum traits have the advantage of being measured and reported on a continuous scale, which removes the issues classically associated with binary disease traits. Overall, genetic selection for improved calf health is feasible. However, to ensure the maximum response, first steps by any industry members should focus efforts on standardising recording practices and encouragement of uploading information to genetic evaluation centres through herd management software, as high-quality phenotypes are the backbone of any successful breeding programme.

Stability of ecologically scaffolded traits during evolutionary transitions in individuality

Evolutionary transitions in individuality are events in the history of life leading to the emergence of new levels of individuality. Recent studies have described an ecological scaffolding scenario of such transitions focused on the evolutionary consequences of an externally imposed renewing meta-population structure with limited dispersal. One difficulty for such a scenario has been explaining the stability of collective-level traits when scaffolding conditions no longer apply. Here, we show that the stability of scaffolded traits can rely on evolutionary hysteresis: even if the environment is reverted to an ancestral state, collectives do not return to ancestral phenotypes. We describe this phenomenon using a stochastic meta-population model and adaptive dynamics. Further, we show that ecological scaffolding may be limited to Goldilocks zones of the environment. We conjecture that Goldilocks zones-even if they might be rare-could act as initiators of evolutionary transitions and help to explain the near ubiquity of collective-level individuality.

Genetic Variation in Male Aggression Is Influenced by Genotype of Prior Social Partners in Drosophila melanogaster

AbstractSocial behaviors can be influenced by the genotypes of interacting individuals through indirect genetic effects (IGEs) and can also display developmental plasticity. We investigated how developmental IGEs, which describe the effects of a prior social partner's genotype on later behavior, can influence aggression in male Drosophila melanogaster. We predicted that developmental IGEs cannot be estimated by simply extending the effects of contextual IGEs over time and instead have their own unique effects on behavior. On day 1 of the experiment, we measured aggressive behavior in 15 genotypic pairings (n = 600 males). On day 2, each of the males was paired with a new opponent, and aggressive behavior was again measured. We found contextual IGEs on day 1 of the experiment and developmental IGEs on day 2 of the experiment: the influence of the day 1 partner's genotype on the focal individual's day 2 behavior depended on the genotypic identity of both the day 1 partner and the focal male. Importantly, the developmental IGEs in our system produced fundamentally different dynamics than the contextual IGEs, as the presence of IGEs was altered over time. These findings represent some of the first empirical evidence demonstrating developmental IGEs, a first step toward incorporating developmental IGEs into our understanding of behavioral evolution.

Indirect genetic effects are shaped by demographic history and ecology in Arabidopsis thaliana

The phenotype of an individual can be affected by the genes of its conspecifics through indirect genetic effects (IGEs). IGEs have been studied across different organisms including wild and domesticated animals and plants, but little is known about their genetic architecture. Here, in a large-scale intraspecific interaction experiment, we show that the contribution of IGEs to the biomass variation of Arabidopsis thaliana is comparable to values classically reported in animals. Moreover, we identify 11 loci explaining 85.1% of the variability in IGEs. We find that positive IGE alleles (that is, those with positive effects on neighbour biomass) occur both in relict accessions from southern Eurasia and in post-glacial colonizers from northern Scandinavia, and that they are likely to have two divergent origins: for nine loci, they evolved in the post-glacial colonizers independently from the relicts, while the two others were introgressed in the post-glacial colonizer from the relicts. Finally, we find that variation in IGEs probably reflects divergent adaptations to the contrasting environments of the edges and the centre of the native range of the species. These findings reveal a surprisingly tractable genetic basis of IGEs in A. thaliana that is shaped by the ecology and the demographic history of the species.

Microbiability of milk composition and genetic control of microbiota effects in sheep

Recently, high-dimensional omics data are becoming available in larger quantities, and models have been developed that integrate them with genomics to understand in finer detail the relationship between genotype and phenotype, and thus improve the performance of genetic evaluations. Our objectives are to quantify the effect of the inclusion of microbiome data in the genetic evaluation for dairy traits in sheep, through the estimation of the heritability, microbiability, and how the microbiome effect on dairy traits decomposes into genetic and nongenetic parts. In this study we analyzed milk and rumen samples of 795 Lacaune dairy ewes. We included, as phenotype, dairy traits and milk fatty acids and proteins composition; as omics measurements, 16S rRNA rumen bacterial abundances; and as genotyping, 54K SNP chip for all ewes. Two nested genomic models were used: a first model to predict the individual contributions of the genetic and microbial abundances to phenotypes, and a second model to predict the additive genetic effect of the microbial community. In addition, microbiome-wide association studies for all dairy traits were applied using the 2,059 rumen bacterial abundances, and the genetic correlations between microbiome principal components and dairy traits were estimated. Results showed that in general the inclusion of both genetic and microbiome effect did not improve the fit of the model compared with the model with the genetic effect only. In addition, for all dairy traits the total heritability was equal to the direct heritability after fitting microbiota effects, due to a microbiability being almost zero for most dairy traits and heritability of the microbial community was very close to zero. Microbiome-wide association studies did not show operational taxonomic units with major effect for any of the dairy traits evaluated, and the genetic correlations between the first 5 principal components and dairy traits were low to moderate. So far, we can conclude that, using a substantial data set of 795 Lacaune dairy ewes, rumen bacterial abundances do not provide improved genetic evaluation for dairy traits in sheep.

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.

The genetic identity of neighboring plants in intraspecific mixtures modulates disease susceptibility of both wheat and rice

Mixing crop cultivars has long been considered as a way to control epidemics at the field level and is experiencing a revival of interest in agriculture. Yet, the ability of mixing to control pests is highly variable and often unpredictable in the field. Beyond classical diversity effects such as dispersal barrier generated by genotypic diversity, several understudied processes are involved. Among them is the recently discovered neighbor-modulated susceptibility (NMS), which depicts the phenomenon that susceptibility in a given plant is affected by the presence of another healthy neighboring plant. Despite the putative tremendous importance of NMS for crop science, its occurrence and quantitative contribution to modulating susceptibility in cultivated species remains unknown. Here, in both rice and wheat inoculated in greenhouse conditions with foliar fungal pathogens considered as major threats, using more than 200 pairs of intraspecific genotype mixtures, we experimentally demonstrate the occurrence of NMS in 11% of the mixtures grown in experimental conditions that precluded any epidemics. Thus, the susceptibility of these 2 major crops results from indirect effects originating from neighboring plants. Quite remarkably, the levels of susceptibility modulated by plant-plant interactions can reach those conferred by intrinsic basal immunity. These findings open new avenues to develop more sustainable agricultural practices by engineering less susceptible crop mixtures thanks to emergent but now predictable properties of mixtures.

Direct additive genetics and maternal effect contribute to the risk of Tourette disorder

BACKGROUND

Risk for Tourette disorder, and chronic motor or vocal tic disorders (referenced here inclusively as CTD), arise from a combination of genetic and environmental factors. While multiple studies have demonstrated the importance of direct additive genetic variation for CTD risk, little is known about the role of cross-generational transmission of genetic risk, such as maternal effect, which is not transmitted via the inherited parental genomes. Here, we partition sources of variation on CTD risk into direct additive genetic effect (narrow-sense heritability) and maternal effect.

METHODS

The study population consists of 2 522 677 individuals from the Swedish Medical Birth Register, who were born in Sweden between 1 January 1973 and 31 December 2000, and followed for a diagnosis of CTD through 31 December, 2013. We used generalised linear mixed models to partition the liability of CTD into: direct additive genetic effect, genetic maternal effect and environmental maternal effect.

RESULTS

We identified 6227 (0.2%) individuals in the birth cohort with a CTD diagnosis. A study of half-siblings showed that maternal half-siblings had twice higher risk of developing a CTD compared with paternal ones. We estimated 60.7% direct additive genetic effect (95% credible interval, 58.5% to 62.4%), 4.8% genetic maternal effect (95% credible interval, 4.4% to 5.1%) and 0.5% environmental maternal effect (95% credible interval, 0.2% to 7%).

CONCLUSIONS

Our results demonstrate genetic maternal effect contributes to the risk of CTD. Failure to account for maternal effect results in an incomplete understanding of the genetic risk architecture of CTD, as the risk for CTD is impacted by maternal effect which is above and beyond the risk from transmitted genetic effect.

Selection for social genetic effects in purebred pigs improves behaviour and handling of their crossbred progeny

BACKGROUND

In commercial pig production, reduction of harmful social behavioural traits, such as ear manipulation and tail biting, is of major interest. Moreover, farmers prefer animals that are easy to handle. The aim of this experiment was to determine whether selection on social breeding values (SBV) for growth rate in purebred pigs affects behaviour in a weighing crate, lesions from ear manipulation, and tail biting of their crossbred progeny. Data were collected on crossbred F1 pigs allocated to 274 pens, which were progeny of purebred Landrace sows and Yorkshire boars from a DanBred nucleus herd.

RESULTS

Behaviour in the weighing crate scored on a three-level scale showed that groups of pigs with high SBV for growth rate were significantly calmer than groups of pigs with low SBV (P < 0.027). When the mean SBV in the group increased by 1 unit, the proportion of pigs that obtained a calmer score level was increased by 14%. A significant (p = 0.04), favourable effect of SBV was found on both the number of pigs with ear lesions in the group and the mean number of ear lesions per pig. For a 1 unit increase in mean SBV, the mean number of lesions per pig decreased by 0.06 from a mean of 0.98. Individual severity of ear lesions conditional upon the number of ear lesions was also significantly affected (p = 0.05) by the mean SBV in the group. In groups for which the mean SBV increased by 1 unit, the proportion of pigs that were observed with a lower severity score was increased by 20% on a three-level scale. Most pigs received no tail biting injuries and no effect of SBV was observed on the tail injury score.

CONCLUSIONS

After 7 weeks in the finisher unit, crossbred progeny with high SBV were calmer in the weighing crate and had fewer ear lesions. These results indicate that selection of purebred parents for SBV for growth rate will increase welfare in their crossbred progeny by decreasing the number of ear lesions and making them easier to handle.

Using social interaction models for genetic analysis of skin damage in gilts

BACKGROUND

Skin damage is a trait of economic and welfare importance that results from social interactions between animals. These interactions may produce wound signs on the gilt's skin as a result of damage behavior (i.e., fighting), biting syndromes (i.e., tail, vulva, or ear biting), and swine inflammation and necrosis syndrome. Although current selection for traits that are affected by social interactions primarily focuses on improving direct genetic effects, combined selection on direct and social genetic effects could increase genetic gain and avoid a negative response to selection in cases of competitive behavior. The objectives of this study were to (1) estimate variance components for combined skin damage (CSD), with or without accounting for social genetic effects, (2) investigate the impact of including genomic information on the prediction accuracy, bias, and dispersion of CSD estimated breeding values, and (3) perform a single-step genome-wide association study (ssGWAS) of CSD under a classical and a social interaction model.

RESULTS

Our results show that CSD is heritable and affected by social genetic effects. Modeling CSD with social interaction models increased the total heritable variance relative to the phenotypic variance by three-fold compared to the classical model. Including genomic information increased the prediction accuracy of direct, social, and total estimated breeding values for purebred sires by at least 21.2%. Bias and dispersion of estimated breeding values were reduced by including genomic information in classical and social interaction models but remained present. The ssGWAS did not identify any single nucleotide polymorphism that was significantly associated with social or direct genetic effects for CSD.

CONCLUSIONS

Combined skin damage is heritable, and genetic selection against this trait will increase the welfare of animals in the long term. Combined skin damage is affected by social genetic effects, and modeling this trait with a social interaction model increases the potential for genetic improvement. Including genomic information increases the prediction accuracy of estimated breeding values and reduces their bias and dispersion, although some biases persist. The results of the genome-wide association study indicate that CSD has a polygenic architecture and no major quantitative trait locus was detected.

The importance of distinguishing individual differences in 'social impact' versus 'social responsiveness' when quantifying indirect genetic effects on the evolution of social plasticity

Social evolution and the dynamics of social interactions have previously been studied under the frameworks of quantitative genetics and behavioural ecology. In quantitative genetics, indirect genetic effects of social partners on the socially plastic phenotypes of focal individuals typically lack crucial detail already included in treatments of social plasticity in behavioural ecology. Specifically, whilst focal individuals (e.g. receivers) may show variation in their 'responsiveness' to the social environment, individual social partners (e.g. signallers) may have a differential 'impact' on focal phenotypes. Here we propose an integrative framework, that highlights the distinction between responsiveness versus impact in indirect genetic effects for a range of behavioural traits. We describe impact and responsiveness using a reaction norm approach and provide statistical models for the assessment of these effects of focal and social partner identity in different types of social interactions. By providing such a framework, we hope to stimulate future quantitative research investigating the causes and consequences of social interactions on phenotypic evolution.

Direct and indirect phenotypic effects on sociability indicate potential to evolve

The decision to leave or join a group is important as group size influences many aspects of organisms' lives and their fitness. This tendency to socialise with others, sociability, should be influenced by genes carried by focal individuals (direct genetic effects) and by genes in partner individuals (indirect genetic effects), indicating the trait's evolution could be slower or faster than expected. However, estimating these genetic parameters is difficult. Here, in a laboratory population of the cockroach Blaptica dubia, I estimate phenotypic parameters for sociability: repeatability (R) and repeatable influence (RI), that indicate whether direct and indirect genetic effects respectively are likely. I also estimate the interaction coefficient (Ψ), which quantifies how strongly a partner's trait influences the phenotype of the focal individual and is key in models for the evolution of interacting phenotypes. Focal individuals were somewhat repeatable for sociability across a 3-week period (R = 0.080), and partners also had marginally consistent effects on focal sociability (RI = 0.053). The interaction coefficient was non-zero, although in opposite sign for the sexes; males preferred to associate with larger individuals (Ψ_male  = -0.129), while females preferred to associate with smaller individuals (Ψ_female  = 0.071). Individual sociability was consistent between dyadic trials and in social networks of groups. These results provide phenotypic evidence that direct and indirect genetic effects have limited influence on sociability, with perhaps most evolutionary potential stemming from heritable effects of the body mass of partners. Sex-specific interaction coefficients may produce sexual conflict and the evolution of sexual dimorphism in social behaviour.

The social evolution of individual differences: Future directions for a comparative science of personality in social behavior

Personality is essential for understanding the evolution of cooperation and conflict in behavior. However, personality science remains disconnected from the field of social evolution, limiting our ability to explain how personality and plasticity shape phenotypic adaptation in social behavior. Researchers also lack an integrative framework for comparing personality in the contextualized and multifaceted behaviors central to social interactions among humans and other animals. Here we address these challenges by developing a social evolutionary approach to personality, synthesizing theory, methods, and organizing questions in the study of individuality and sociality in behavior. We critically review current measurement practices and introduce social reaction norm models for comparative research on the evolution of personality in social environments. These models demonstrate that social plasticity affects the heritable variance of personality, and that individual differences in social plasticity can further modify the rate and direction of adaptive social evolution. Future empirical studies of frequency- and density-dependent social selection on personality are crucial for further developing this framework and testing adaptive theory of social niche specialization.

Linking quantitative genetics with community-level performance: Are there operational models for plant breeding?

Plant breeding is focused on the genotype and population levels while targeting effects at higher levels of biodiversity, from crop covers to agroecosystems. Making predictions across nested levels of biodiversity is therefore a major challenge for the development of intercropping practices. New prediction tools and concepts are required to design breeding strategies with desirable outcomes at the crop community level. We reviewed theoretical advances in the field of evolutionary ecology to identify potentially operational ways of predicting the effects of artificial selection on community-level performances. We identified three main types of approaches differing in the way they model interspecific indirect genetic effects (IIGEs) at the community level: (1) The community heritability approach estimates the variance for IIGE induced by a focal species at the community level; (2) the joint phenotype approach quantifies genetic constraints between direct genetic effects and IIGE for a set of interacting species; (3) the community-trait genetic gradient approach decomposes the IIGE for a focal species across a multivariate set of its functional traits. We discuss the potential operational capacities of these approaches and stress that each is a special case of a general multitrait and multispecies selection index. Choosing one therefore involves assumptions and goals regarding the breeding target and strategy. Obtaining reliable quantitative, community-level predictions at the genetic level is constrained by the size and complexity of the experimental designs usually required. Breeding strategies should instead be compared using theoretically informed qualitative predictions. The need to estimate genetic covariances between traits measured both within and among species (for IIGE) is another obstacle, as the two are not determined by the exact same biological processes. We suggest future research directions and strategies to overcome these limits. Our synthesis offers an integrative theoretical framework for breeders interested in the genetic improvement of crop communities but also for scientists interested in the genetic bases of plant community functioning.

Kin selection theory and the design of cooperative crops

In agriculture and plant breeding, plant traits may be favoured because they benefit neighbouring plants and ultimately increase total crop yield. This idea of promoting cooperation among crop plants has existed almost as long as W.D. Hamilton's inclusive fitness (kin selection) theory, the leading framework for explaining cooperation in biology. However, kin selection thinking has not been adequately applied to the idea of cooperative crops. Here, I give an overview of modern kin selection theory and consider how it explains three key strategies for designing cooperative crops: (1) selection for a less-competitive plant type (a 'communal ideotype'); (2) group-level selection for yield; and (3) exploiting naturally selected cooperation. The first two strategies, using artificial selection, have been successful in the past but suffer from limitations that could hinder future progress. Instead, I propose an alternative strategy and a new 'colonial ideotype' that exploits past natural selection for cooperation among the modules (e.g., branches or stems) of individual plants. More generally, I suggest that Hamiltonian agriculture-a kin selection view of agriculture and plant breeding-transforms our understanding of how to improve crops of the future.

On the importance of parenting in externalizing disorders: an evaluation of indirect genetic effects in families

BACKGROUND

Theoretical models of the development of childhood externalizing disorders emphasize the role of parents. Empirical studies have not been able to identify specific aspects of parental behaviors explaining a considerable proportion of the observed individual differences in externalizing problems. The problem is complicated by the contribution of genetic factors to externalizing problems, as parents provide both genes and environments to their children. We studied the joint contributions of direct genetic effects of children and the indirect genetic effects of parents through the environment on externalizing problems.

METHODS

The study used genome-wide single nucleotide polymorphism data from 9,675 parent-offspring trios participating in the Norwegian Mother Father and child cohort study. Based on genomic relatedness matrices, we estimated the contribution of direct genetic effects and indirect maternal and paternal genetic effects on ADHD, conduct and disruptive behaviors at 8 years of age.

RESULTS

Models including indirect parental genetic effects were preferred for the ADHD symptoms of inattention and hyperactivity, and conduct problems, but not oppositional defiant behaviors. Direct genetic effects accounted for 11% to 24% of the variance, whereas indirect parental genetic effects accounted for 0% to 16% in ADHD symptoms and conduct problems. The correlation between direct and indirect genetic effects, or gene-environment correlations, decreased the variance with 16% and 13% for conduct and inattention problems, and increased the variance with 6% for hyperactivity problems.

CONCLUSIONS

This study provides empirical support to the notion that parents have a significant role in the development of childhood externalizing behaviors. The parental contribution to decrease in variation of inattention and conduct problems by gene-environment correlations would limit the number of children reaching clinical ranges in symptoms. Not accounting for indirect parental genetic effects can lead to both positive and negative bias when identifying genetic variants for childhood externalizing behaviors.

Niche construction in quantitative traits: heritability and response to selection

A central tenet of niche construction (NC) theory is that organisms can alter their environments in heritable and evolutionarily important ways, often altering selection pressures. We suggest that the physical changes niche constructors make to their environments may also alter trait heritability and the response of phenotypes to selection. This effect might change evolution, over and above the effect of NC acting via selection alone. We develop models of trait evolution that allow us to partition the effects of NC on trait heritability from those on selection to better investigate their distinct effects. We show that the response of a phenotype to selection and so the pace of phenotypic change can be considerably altered in the presence of NC and that this effect is compounded when trans-generational interactions are included. We argue that novel mathematical approaches are needed to describe the simultaneous effects of NC on trait evolution via selection and heritability. Just as indirect genetic effects have been shown to significantly increase trait heritability, the effects of NC on heritability in our model suggest a need for further theoretical development of the concept of heritability.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Capturing indirect genetic effects on phenotypic variability: Competition meets canalization

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is possible with models for inherited variability but may require a two-step analysis.

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.

The Indirect Genetic Effect Interaction Coefficient ψ: Theoretically Essential and Empirically Neglected

The interaction effect coefficient ψ has been a much-discussed, fundamental parameter of indirect genetic effect (IGE) models since its formal mathematical description in 1997. The coefficient simultaneously describes the form of changes in trait expression caused by genes in the social environment and predicts the evolutionary consequences of those IGEs. Here, we report a striking mismatch between theoretical emphasis on ψ and its usage in empirical studies. Surveying all IGE research, we find that the coefficient ψ has not been equivalently conceptualized across studies. Several issues related to its proper empirical measurement have recently been raised, and these may severely distort interpretations about the evolutionary consequences of IGEs. We provide practical advice on avoiding such pitfalls. The majority of empirical IGE studies use an alternative variance-partitioning approach rooted in well-established statistical quantitative genetics, but several hundred estimates of ψ (from 15 studies) have been published. A significant majority are positive. In addition, IGEs with feedback, that is, involving the same trait in both interacting partners, are far more likely to be positive and of greater magnitude. Although potentially challenging to measure without bias, ψ has critically-developed theoretical underpinnings that provide unique advantages for empirical work. We advocate for a shift in perspective for empirical work, from ψ as a description of IGEs, to ψ as a robust predictor of evolutionary change. Approaches that "run evolution forward" can take advantage of ψ to provide falsifiable predictions about specific trait interactions, providing much-needed insight into the evolutionary consequences of IGEs.

Indirect Genetic Effects: A Cross-disciplinary Perspective on Empirical Studies

Indirect genetic effects (IGE) occur when an individual's phenotype is influenced by genetic variation in conspecifics. Opportunities for IGE are ubiquitous, and, when present, IGE have profound implications for behavioral, evolutionary, agricultural, and biomedical genetics. Despite their importance, the empirical study of IGE lags behind the development of theory. In large part, this lag can be attributed to the fact that measuring IGE, and deconvoluting them from the direct genetic effects of an individual's own genotype, is subject to many potential pitfalls. In this Perspective, we describe current challenges that empiricists across all disciplines will encounter in measuring and understanding IGE. Using ideas and examples spanning evolutionary, agricultural, and biomedical genetics, we also describe potential solutions to these challenges, focusing on opportunities provided by recent advances in genomic, monitoring, and phenotyping technologies. We hope that this cross-disciplinary assessment will advance the goal of understanding the pervasive effects of conspecific interactions in biology.

When is Offspring Viability Fitness a Measure of Paternal Fitness and When is it not?

We use population genetics to model the evolution of a gene with an indirect effect owing to paternal care and with a second pleiotropic, direct effect on offspring viability. We use the model to illustrate how the common empirical practice of considering offspring viability as a component of parent fitness can confound a gene's direct and indirect fitness effects. We investigate when this confounding results in a distorted picture of overall evolution and when it does not. We find that the practice has no effect on mean fitness, W, but it does have an effect on the dynamics of gene frequency change, ∆q. We also find that, for some regions of parameter space associated with fitness trade-offs, the distortion is not only quantitative but also qualitative, obscuring the direction of gene frequency change. Because it affects the evolutionary dynamics, it also affects the expected amount of genetic variation at mutation-selection balance, an important consideration in molecular evolution. We discuss empirical techniques for separating direct from indirect effects and how field studies measuring the value of male paternal care might be improved by using them.

A Synthesis of Game Theory and Quantitative Genetic Models of Social Evolution

Two popular approaches for modeling social evolution, evolutionary game theory and quantitative genetics, ask complementary questions but are rarely integrated. Game theory focuses on evolutionary outcomes, with models solving for evolutionarily stable equilibria, whereas quantitative genetics provides insight into evolutionary processes, with models predicting short-term responses to selection. Here we draw parallels between evolutionary game theory and interacting phenotypes theory, which is a quantitative genetic framework for understanding social evolution. First, we show how any evolutionary game may be translated into two quantitative genetic selection gradients, nonsocial and social selection, which may be used to predict evolutionary change from a single round of the game. We show that synergistic fitness effects may alter predicted selection gradients, causing changes in magnitude and sign as the population mean evolves. Second, we show how evolutionary games involving plastic behavioral responses to partners can be modeled using indirect genetic effects, which describe how trait expression changes in response to genes in the social environment. We demonstrate that repeated social interactions in models of reciprocity generate indirect effects and conversely, that estimates of parameters from indirect genetic effect models may be used to predict the evolution of reciprocity. We argue that a pluralistic view incorporating both theoretical approaches will benefit empiricists and theorists studying social evolution. We advocate the measurement of social selection and indirect genetic effects in natural populations to test the predictions from game theory and, in turn, the use of game theory models to aid in the interpretation of quantitative genetic estimates.

Rootstock-Mediated Genetic Variance in Cadmium Uptake by Juvenile Cacao (Theobroma cacao L.) Genotypes, and Its Effect on Growth and Physiology

Grafting typically offers a shortcut to breed tree orchards throughout a multidimensional space of traits. Despite an overwhelming spectrum of rootstock-mediated effects on scion traits observed across several species, the exact nature and mechanisms underlying the rootstock-mediated effects on scion traits in cacao (Theobroma cacao L.) plants often remain overlooked. Therefore, we aimed to explicitly quantify rootstock-mediated genetic contributions in recombinant juvenile cacao plants across target traits, specifically cadmium (Cd) uptake, and its correlation with growth and physiological traits. Content of chloroplast pigments, fluorescence of chlorophyll a, leaf gas exchange, nutrient uptake, and plant biomass were examined across ungrafted saplings and target rootstock × scion combinations in soils with contrasting levels of Cd. This panel considered a total of 320 progenies from open-pollinated half-sib families and reciprocal full-sib progenies (derived from controlled crosses between the reference genotypes IMC67 and PA121). Both family types were used as rootstocks in grafts with two commercial clones (ICS95 and CCN51) commonly grown in Colombia. A pedigree-based best linear unbiased prediction (A-BLUP) mixed model was implemented to quantify rootstock-mediated narrow-sense heritability (h 2) for target traits. A Cd effect measured on rootstocks before grafting was observed in plant biomass, nutrient uptake, and content of chloroplast pigments. After grafting, damage to the Photosystem II (PSII) was also evident in some rootstock × scion combinations. Differences in the specific combining ability for Cd uptake were mostly detected in ungrafted rootstocks, or 2 months after grafting with the clonal CCN51 scion. Moderate rootstock effects (h 2> 0.1) were detected before grafting for five growth traits, four nutrient uptake properties, and chlorophylls and carotenoids content (h 2 = 0.19, 95% CI 0.05-0.61, r = 0.7). Such rootstock effects faded (h 2< 0.1) when rootstock genotypes were examined in soils without Cd, or 4 months after grafting. These results suggest a pervasive genetic conflict between the rootstock and the scion genotypes, involving the triple rootstock × scion × soil interaction when it refers to Cd and nutrient uptake, early growth, and photosynthetic process in juvenile cacao plants. Overall, deepening on these findings will harness early breeding schemes of cacao rootstock genotypes compatible with commercial clonal scions and adapted to soils enriched with toxic levels of Cd.

Contributions of genetic and nongenetic sources to variation in cooperative behavior in a cooperative mammal

The evolution of cooperative behavior is a major area of research among evolutionary biologists and behavioral ecologists, yet there are few estimates of its heritability or its evolutionary potential, and long-term studies of identifiable individuals are required to disentangle genetic and nongenetic components of cooperative behavior. Here, we use long-term data on over 1800 individually recognizable wild meerkats (Suricata suricatta) collected over 30 years and a multigenerational genetic pedigree to partition phenotypic variation in three cooperative behaviors (babysitting, pup feeding, and sentinel behavior) into individual, additive genetic, and other sources, and to assess their repeatability and heritability. In addition to strong effects of sex, age, and dominance status, we found significant repeatability in individual contributions to all three types of cooperative behavior both within and across breeding seasons. Like most other studies of the heritability of social behavior, we found that the heritability of cooperative behavior was low. However, our analysis suggests that a substantial component of the repeatable individual differences in cooperative behavior that we observed was a consequence of additive genetic variation. Our results consequently indicate that cooperative behavior can respond to selection, and suggest scope for further exploration of the genetic basis of social behavior.

The Collective Behavior of Ant Groups Depends on Group Genotypic Composition

Recently, researchers have documented variation between groups in collective behavior. However, how genetic variation within and between groups contributes to population-level variation for collective behavior remains unclear. Understanding how genetic variation of group members relates to group-level phenotypes is evolutionarily important because there is increasing evidence that group-level behavioral variation influences fitness and that the genetic architecture of group-level traits can affect the evolutionary dynamics of traits. Social insects are ideal for studying the complex relationship between individual and group-level variation because they exhibit behavioral variation at multiple scales of organization. To explore how the genetic composition of groups affects collective behavior, we constructed groups of pharaoh ants (Monomorium pharaonis) from 33 genetically distinct colonies of known pedigree. The groups consisted of either all workers from the same single colony or workers from two genetically different colonies, and we assayed the exploration and aggression of the groups. We found that collective exploration, but not aggression, depended on the specific genotypic combination of group members, i.e., we found evidence for genotype-by-genotype epistasis for exploration. Group collective behavior did not depend on the pedigree relatedness between genotypes within groups. Overall, this study highlights that specific combinations of genotypes influence group-level phenotypes, emphasizing the importance of considering nonadditive effects of genotypic interactions between group members.

A neglected conceptual problem regarding phenotypic plasticity's role in adaptive evolution: The importance of genetic covariance and social drive

There is tantalizing evidence that phenotypic plasticity can buffer novel, adaptive genetic variants long enough to permit their evolutionary spread, and this process is often invoked in explanations for rapid adaptive evolution. However, the strength and generality of evidence for it is controversial. We identify a conceptual problem affecting this debate: recombination, segregation, and independent assortment are expected to quickly sever associations between genes controlling novel adaptations and genes contributing to trait plasticity that facilitates the novel adaptations by reducing their indirect fitness costs. To make clearer predictions about this role of plasticity in facilitating genetic adaptation, we describe a testable genetic mechanism that resolves the problem: genetic covariance between new adaptive variants and trait plasticity that facilitates their persistence within populations. We identify genetic architectures that might lead to such a covariance, including genetic coupling via physical linkage and pleiotropy, and illustrate the consequences for adaptation rates using numerical simulations. Such genetic covariances may also arise from the social environment, and we suggest the indirect genetic effects that result could further accentuate the process of adaptation. We call the latter mechanism of adaptation social drive, and identify methods to test it. We suggest that genetic coupling of plasticity and adaptations could promote unusually rapid ‘runaway’ evolution of novel adaptations. The resultant dynamics could facilitate evolutionary rescue, adaptive radiations, the origin of novelties, and other commonly studied processes.

Genetic Consequences of Biologically Altered Environments

Evolvable traits of organisms can alter the environment those organisms experience. While it is well appreciated that those modified environments can influence natural selection to which organisms are exposed, they can also influence the expression of genetic variances and covariances of traits under selection. When genetic variance and covariance change in response to changes in the evolving, modified environment, rates and outcomes of evolution also change. Here we discuss the basic mechanisms whereby organisms modify their environments, review how those modified environments have been shown to alter genetic variance and covariance, and discuss potential evolutionary consequences of such dynamics. With these dynamics, responses to selection can be more rapid and sustained, leading to more extreme phenotypes, or they can be slower and truncated, leading to more conserved phenotypes. Patterns of correlated selection can also change, leading to greater or less evolutionary independence of traits, or even causing convergence or divergence of traits, even when selection on them is consistent across environments. Developing evolutionary models that incorporate changes in genetic variances and covariances when environments themselves evolve requires developing methods to predict how genetic parameters respond to environments-frequently multifactorial environments. It also requires a population-level analysis of how traits of collections of individuals modify environments for themselves and/or others in a population, possibly in spatially explicit ways. Despite the challenges of elucidating the mechanisms and nuances of these processes, even qualitative predictions of how environment-modifying traits alter evolutionary potential are likely to improve projections of evolutionary outcomes.

Genetic parameter estimates for bull prolificacy and its relationship with scrotal circumference in a commercial beef cattle population

The commercial beef cattle industry relies heavily on the use of natural service sires. When artificial insemination is deemed difficult to implement, multisire breeding pastures are used to increase reproductive rates in large breeding herds or to safe-guard against bull injury during the breeding season. Although each bull might be given an equal opportunity to produce offspring, evidence suggest that there is substantial variation in the number of calves sired by each bull in a breeding pasture. With the use of DNA-based paternity testing, correctly assigning calves to their respective sires in multisire pastures is possible and presents an opportunity to investigate the degree to which this trait complex is under genetic control. Field data from a large commercial ranch was used to estimate genetic parameters for calf count (CC; 574 records from 443 sires) and yearling scrotal circumference (SC; n = 1961) using univariate and bivariate animal models. Calf counts averaged 12.2 ± 10.7 and SC averaged 35.4 ± 2.30 cm. Bulls had an average of 1.30 records and there were 23.9 ± 11.1 bulls per contemporary group. The model for CC included fixed effects of age during the breeding season (in years) and contemporary group (concatenation of breeding pasture and year). Random effects included additive genetic and permanent environmental effects, and a residual. The model for SC included fixed effects of age (in days) and contemporary group (concatenation of month and year of measurement). Random effects included an additive genetic effect and a residual. Univariate model heritability estimates for CC and SC were 0.178 ± 0.142 and 0.455 ± 0.072, respectively. Similarly, the bivariate model resulted in heritability estimates for CC and SC of 0.184 ± 0.142 and 0.457 ± 0.072, respectively. Repeatability estimates for CC from univariate and bivariate models were 0.315 ± 0.080 and 0.317 ± 0.080, respectively. The estimate of genetic correlation between CC and SC was 0.268 ± 0.274. Heritability estimates suggest that both CC and SC would respond favorably to selection. Moreover, CC is lowly repeatable and although favorably correlated, SC appears to be weakly associated with CC.

The thermal environment as a moderator of social evolution

Animal sociality plays a crucial organisational role in evolution. As a result, understanding the factors that promote the emergence, maintenance, and diversification of animal societies is of great interest to biologists. Climate is among the foremost ecological factors implicated in evolutionary transitions in social organisation, but we are only beginning to unravel the possible mechanisms and specific climatic variables that underlie these associations. Ambient temperature is a key abiotic factor shaping the spatio-temporal distribution of individuals and has a particularly strong influence on behaviour. Whether such effects play a broader role in social evolution remains to be seen. In this review, we develop a conceptual framework for understanding how thermal effects integrate into pathways that mediate the opportunities, nature, and context of social interactions. We then implement this framework to discuss the capacity for temperature to initiate organisational changes across three broad categories of social evolution: group formation, group maintenance, and group elaboration. For each category, we focus on pivotal traits likely to have underpinned key social transitions and explore the potential for temperature to affect changes in these traits by leveraging empirical examples from the literature on thermal and behavioural ecology. Finally, we discuss research directions that should be prioritised to understand the potentially constructive and/or destructive effects of future warming on the origins, maintenance, and diversification of animal societies.

Dissecting indirect genetic effects from peers in laboratory mice

BACKGROUND

The phenotype of an individual can be affected not only by the individual's own genotypes, known as direct genetic effects (DGE), but also by genotypes of interacting partners, indirect genetic effects (IGE). IGE have been detected using polygenic models in multiple species, including laboratory mice and humans. However, the underlying mechanisms remain largely unknown. Genome-wide association studies of IGE (igeGWAS) can point to IGE genes, but have not yet been applied to non-familial IGE arising from "peers" and affecting biomedical phenotypes. In addition, the extent to which igeGWAS will identify loci not identified by dgeGWAS remains an open question. Finally, findings from igeGWAS have not been confirmed by experimental manipulation.

RESULTS

We leverage a dataset of 170 behavioral, physiological, and morphological phenotypes measured in 1812 genetically heterogeneous laboratory mice to study IGE arising between same-sex, adult, unrelated mice housed in the same cage. We develop and apply methods for igeGWAS in this context and identify 24 significant IGE loci for 17 phenotypes (FDR < 10%). We observe no overlap between IGE loci and DGE loci for the same phenotype, which is consistent with the moderate genetic correlations between DGE and IGE for the same phenotype estimated using polygenic models. Finally, we fine-map seven significant IGE loci to individual genes and find supportive evidence in an experiment with a knockout model that Epha4 gives rise to IGE on stress-coping strategy and wound healing.

CONCLUSIONS

Our results demonstrate the potential for igeGWAS to identify IGE genes and shed light into the mechanisms of peer influence.

A combined GWAS approach reveals key loci for socially-affected traits in Yorkshire pigs

Socially affected traits in pigs are controlled by direct genetic effects and social genetic effects, which can make elucidation of their genetic architecture challenging. We evaluated the genetic basis of direct genetic effects and social genetic effects by combining single-locus and haplotype-based GWAS on imputed whole-genome sequences. Nineteen SNPs and 25 haplotype loci are identified for direct genetic effects on four traits: average daily feed intake, average daily gain, days to 100 kg and time in feeder per day. Nineteen SNPs and 11 haplotype loci are identified for social genetic effects on average daily feed intake, average daily gain, days to 100 kg and feeding speed. Two significant SNPs from single-locus GWAS (SSC6:18,635,874 and SSC6:18,635,895) are shared by a significant haplotype locus with haplotype alleles 'GGG' for both direct genetic effects and social genetic effects in average daily feed intake. A candidate gene, MT3, which is involved in growth, nervous, and immune processes, is identified. We demonstrate the genetic differences between direct genetic effects and social genetic effects and provide an anchor for investigating the genetic architecture underlying direct genetic effects and social genetic effects on socially affected traits in pigs.

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.

Territoriality in Drosophila: indirect effects and covariance with body mass and metabolic rate

Territoriality (i.e., defense of a resource) is the outcome of behavioral interactions that can result in selective advantages in many vertebrates and invertebrates. Since territoriality is expressed in a social context, an individuals’ territoriality may change according to the phenotype of the opponents that they are confronted with (termed “indirect effects”). Defending a territory may also confer energetic costs to individuals, which could be reflected in their standard metabolic rate (SMR), a key component of an ectotherms’ energy budget. Here, we measured territoriality using dyadic contests, body mass, and SMR using flow-through respirometry, twice in each of 192 adult male Drosophila melanogaster. Territoriality, body mass, and (whole-animal) SMR were all significantly repeatable. However, essentially all the among-individual variation in SMR was shared with body mass, as indicated by a very strong among-individual correlation (rind) between body mass and SMR. The among-individual correlation between territoriality and SMR also tended to be positive, suggesting the presence of underlying metabolic costs to territoriality. Although indirect effects on territoriality were present but weak (accounting for 8.4% of phenotypic variance), indirect effects on territoriality were negatively and significantly correlated with body mass. This indicates that larger individuals tended to suppress their opponents territoriality. Variation among individuals in their ability to suppress territoriality in others was not associated with their own territoriality or SMR.

Phenomes: the current frontier in animal breeding

Improvements in genomic technologies have outpaced the most optimistic predictions, allowing industry-scale application of genomic selection. However, only marginal gains in genetic prediction accuracy can now be expected by increasing marker density up to sequence, unless causative mutations are identified. We argue that some of the most scientifically disrupting and industry-relevant challenges relate to ‘phenomics’ instead of ‘genomics’. Thanks to developments in sensor technology and artificial intelligence, there is a wide range of analytical tools that are already available and many more will be developed. We can now address some of the pressing societal demands on the industry, such as animal welfare concerns or efficiency in the use of resources. From the statistical and computational point of view, phenomics raises two important issues that require further work: penalization and dimension reduction. This will be complicated by the inherent heterogeneity and ‘missingness’ of the data. Overall, we can expect that precision livestock technologies will make it possible to collect hundreds of traits on a continuous basis from large numbers of animals. Perhaps the main revolution will come from redesigning animal breeding schemes to explicitly allow for high-dimensional phenomics. In the meantime, phenomics data will definitely enlighten our knowledge on the biological basis of phenotypes.

Why genetic selection to reduce the prevalence of infectious diseases is way more promising than currently believed

Genetic selection for improved disease resistance is an important part of strategies to combat infectious diseases in agriculture. Quantitative genetic analyses of binary disease status, however, indicate low heritability for most diseases, which restricts the rate of genetic reduction in disease prevalence. Moreover, the common liability threshold model suggests that eradication of an infectious disease via genetic selection is impossible because the observed-scale heritability goes to zero when the prevalence approaches zero. From infectious disease epidemiology, however, we know that eradication of infectious diseases is possible, both in theory and practice, because of positive feedback mechanisms leading to the phenomenon known as herd immunity. The common quantitative genetic models, however, ignore these feedback mechanisms. Here, we integrate quantitative genetic analysis of binary disease status with epidemiological models of transmission, aiming to identify the potential response to selection for reducing the prevalence of endemic infectious diseases. The results show that typical heritability values of binary disease status correspond to a very substantial genetic variation in disease susceptibility among individuals. Moreover, our results show that eradication of infectious diseases by genetic selection is possible in principle. These findings strongly disagree with predictions based on common quantitative genetic models, which ignore the positive feedback effects that occur when reducing the transmission of infectious diseases. Those feedback effects are a specific kind of Indirect Genetic Effects; they contribute substantially to the response to selection and the development of herd immunity (i.e., an effective reproduction ratio less than one).

Evidence of horizontal indirect genetic effects in humans

Indirect genetic effects, the effects of the genotype of one individual on the phenotype of other individuals, are environmental factors associated with human disease and complex trait variation that could help to expand our understanding of the environment linked to complex traits. Here, we study indirect genetic effects in 80,889 human couples of European ancestry for 105 complex traits. Using a linear mixed model approach, we estimate partner indirect heritability and find evidence of partner heritability on ~50% of the analysed traits. Follow-up analysis suggests that in at least ~25% of these traits, the partner heritability is consistent with the existence of indirect genetic effects including a wide variety of traits such as dietary traits, mental health and disease. This shows that the environment linked to complex traits is partially explained by the genotype of other individuals and motivates the need to find new ways of studying the environment.

Direct and Indirect Effects of Maternal, Paternal, and Offspring Genotypes: Trio-GCTA

Indirect genetic effects from relatives may result in misleading quantifications of heritability, but can also be of interest in their own right. In this paper we propose Trio-GCTA, a model for separating direct and indirect genetic effects when genome-wide single nucleotide polymorphism data have been collected from parent-offspring trios. The model is applicable to phenotypes obtained from any of the family members. We discuss appropriate parameter interpretations and apply the method to three exemplar phenotypes: offspring birth weight, maternal relationship satisfaction, and paternal body-mass index, using real data from the Norwegian Mother, Father and Child Cohort Study (MoBa).

Breeding Beyond Monoculture: Putting the "Intercrop" Into Crops

Intercropping is both a well-established and yet novel agricultural practice, depending on one's perspective. Such perspectives are principally governed by geographic location and whether monocultural practices predominate. Given the negative environmental effects of monoculture agriculture (loss of biodiversity, reliance on non-renewable inputs, soil degradation, etc.), there has been a renewed interest in cropping systems that can reduce the impact of modern agriculture while maintaining (or even increasing) yields. Intercropping is one of the most promising practices in this regard, yet faces a multitude of challenges if it is to compete with and ultimately replace the prevailing monocultural norm. These challenges include the necessity for more complex agricultural designs in space and time, bespoke machinery, and adapted crop cultivars. Plant breeding for monocultures has focused on maximizing yield in single-species stands, leading to highly productive yet specialized genotypes. However, indications suggest that these genotypes are not the best adapted to intercropping systems. Re-designing breeding programs to accommodate inter-specific interactions and compatibilities, with potentially multiple different intercropping partners, is certainly challenging, but recent technological advances offer novel solutions. We identify a number of such technology-driven directions, either ideotype-driven (i.e., "trait-based" breeding) or quantitative genetics-driven (i.e., "product-based" breeding). For ideotype breeding, plant growth modeling can help predict plant traits that affect both inter- and intraspecific interactions and their influence on crop performance. Quantitative breeding approaches, on the other hand, estimate breeding values of component crops without necessarily understanding the underlying mechanisms. We argue that a combined approach, for example, integrating plant growth modeling with genomic-assisted selection and indirect genetic effects, may offer the best chance to bridge the gap between current monoculture breeding programs and the more integrated and diverse breeding programs of the future.

Optimal definition of contemporary groups for crossbred pigs in a joint purebred and crossbred genetic evaluation

In the pig industry, purebred animals are raised in nucleus herds and selected to produce crossbred progeny to perform in commercial environments. Crossbred and purebred performances are different, correlated traits. All purebreds in a pen have their performance assessed together at the end of a performance test. However, only selected crossbreds are removed (based on visual inspection) and measured at different times creating many small contemporary groups (CGs). This may reduce estimated breeding value (EBV) prediction accuracies. Considering this sequential recording of crossbreds, the objective was to investigate the impact of different CG definitions on genetic parameters and EBV prediction accuracy for crossbred traits. Growth rate (GP) and ultrasound backfat (BFP) records were available for purebreds. Lifetime growth (GX) and backfat (BFX) were recorded on crossbreds. Different CGs were tested: CG_all included farm, sex, birth year, and birth week; CG_week added slaughter week; and CG_day used slaughter day instead of week. Data of 124,709 crossbreds were used. The purebred phenotypes (62,274 animals) included three generations of purebred ancestors of these crossbreds and their CG mates. Variance components for four-trait models with different CG definitions were estimated with average information restricted maximum likelihood. Purebred traits' variance components remained stable across CG definitions and varied slightly for BFX. Additive genetic variances (and heritabilities) for GX fluctuated more: 812 ± 36 (0.28 ± 0.01), 257 ± 15 (0.17 ± 0.01), and 204 ± 13 (0.15 ± 0.01) for CG_all, CG_week, and CG_day, respectively. Age at slaughter (AAS) and hot carcass weight (HCW) adjusted for age were investigated as alternatives for GX. Both have potential for selection but lower heritabilities compared with GX: 0.21 ± 0.01 (0.18 ± 0.01), 0.16 ± 0.02 (0.16 + 0.01), and 0.10 ± 0.01 (0.14 ± 0.01) for AAS (HCW) using CG_all, CG_week, and CG_day, respectively. The predictive ability, linear regression (LR) accuracy, bias, and dispersion of crossbred traits in crossbreds favored CG_day, but correlations with unadjusted phenotypes favored CG_all. In purebreds, CG_all showed the best LR accuracy, while showing small relative differences in bias and dispersion. Different CG scenarios showed no relevant impact on BFX EBV. This study shows that different CG definitions may affect evaluation stability and animal ranking. Results suggest that ignoring slaughter dates in CG is more appropriate for estimating crossbred trait EBV for purebred animals.

Inheritance of Rootstock Effects in Avocado (Persea americana Mill.) cv. Hass

Grafting is typically utilized to merge adapted seedling rootstocks with highly productive clonal scions. This process implies the interaction of multiple genomes to produce a unique tree phenotype. However, the interconnection of both genotypes obscures individual contributions to phenotypic variation (rootstock-mediated heritability), hampering tree breeding. Therefore, our goal was to quantify the inheritance of seedling rootstock effects on scion traits using avocado (Persea americana Mill.) cv. Hass as a model fruit tree. We characterized 240 diverse rootstocks from 8 avocado cv. Hass orchards with similar management in three regions of the province of Antioquia, northwest Andes of Colombia, using 13 microsatellite markers simple sequence repeats (SSRs). Parallel to this, we recorded 20 phenotypic traits (including morphological, biomass/reproductive, and fruit yield and quality traits) in the scions for 3 years (2015-2017). Relatedness among rootstocks was inferred through the genetic markers and inputted in a "genetic prediction" model to calculate narrow-sense heritabilities (h 2) on scion traits. We used three different randomization tests to highlight traits with consistently significant heritability estimates. This strategy allowed us to capture five traits with significant heritability values that ranged from 0.33 to 0.45 and model fits (r) that oscillated between 0.58 and 0.73 across orchards. The results showed significance in the rootstock effects for four complex harvest and quality traits (i.e., total number of fruits, number of fruits with exportation quality, and number of fruits discarded because of low weight or thrips damage), whereas the only morphological trait that had a significant heritability value was overall trunk height (an emergent property of the rootstock-scion interaction). These findings suggest the inheritance of rootstock effects, beyond root phenotype, on a surprisingly wide spectrum of scion traits in "Hass" avocado. They also reinforce the utility of polymorphic SSRs for relatedness reconstruction and genetic prediction of complex traits. This research is, up to date, the most cohesive evidence of narrow-sense inheritance of rootstock effects in a tropical fruit tree crop. Ultimately, our work highlights the importance of considering the rootstock-scion interaction to broaden the genetic basis of fruit tree breeding programs while enhancing our understanding of the consequences of grafting.

Prospects for the Analysis and Reduction of Damaging Behaviour in Group-Housed Livestock, With Application to Pig Breeding

Innovations in the breeding and management of pigs are needed to improve the performance and welfare of animals raised in social groups, and in particular to minimise biting and damage to group mates. Depending on the context, social interactions between pigs can be frequent or infrequent, aggressive, or non-aggressive. Injuries or emotional distress may follow. The behaviours leading to damage to conspecifics include progeny savaging, tail, ear or vulva biting, and excessive aggression. In combination with changes in husbandry practices designed to improve living conditions, refined methods of genetic selection may be a solution reducing these behaviours. Knowledge gaps relating to lack of data and limits in statistical analyses have been identified. The originality of this paper lies in its proposal of several statistical methods for common use in analysing and predicting unwanted behaviours, and for genetic use in the breeding context. We focus on models of interaction reflecting the identity and behaviour of group mates which can be applied directly to damaging traits, social network analysis to define new and more integrative traits, and capture-recapture analysis to replace missing data by estimating the probability of behaviours. We provide the rationale for each method and suggest they should be combined for a more accurate estimation of the variation underlying damaging behaviours.

Prediction of genetic merit for growth rate in pigs using animal models with indirect genetic effects and genomic information

Background

Several studies have found that the growth rate of a pig is influenced by the genetics of the group members (indirect genetic effects). Accounting for these indirect genetic effects in a selection program may increase genetic progress for growth rate. However, indirect genetic effects are small and difficult to predict accurately. Genomic information may increase the ability to predict indirect genetic effects. Thus, the objective of this study was to test whether including indirect genetic effects in the animal model increases the predictive performance when genetic effects are predicted with genomic relationships. In total, 11,255 pigs were phenotyped for average daily gain between 30 and 94 kg, and 10,995 of these pigs were genotyped. Two relationship matrices were used: a numerator relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A) and a combined pedigree and genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H); and two different animal models were used: an animal model with only direct genetic effects and an animal model with both direct and indirect genetic effects. The predictive performance of the models was defined as the Pearson correlation between corrected phenotypes and predicted genetic levels. The predicted genetic level of a pig was either its direct genetic effect or the sum of its direct genetic effect and the indirect genetic effects of its group members (total genetic effect).

Results

The highest predictive performance was achieved when total genetic effects were predicted with genomic information (21.2 vs. 14.7%). In general, the predictive performance was greater for total genetic effects than for direct genetic effects (0.1 to 0.5% greater; not statistically significant). Both types of genetic effects had greater predictive performance when they were predicted with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H rather than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A (5.9 to 6.3%). The difference between predictive performances of total genetic effects and direct genetic effects was smaller when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H was used rather than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A.

Conclusions

This study provides evidence that: (1) corrected phenotypes are better predicted with total genetic effects than with direct genetic effects only; (2) both direct genetic effects and indirect genetic effects are better predicted with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A; (3) using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H rather than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A primarily improves the predictive performance of direct genetic effects.

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

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

Feed competition reduces heritable variation for body weight in Litopenaeus vannamei

Background

Competition is a common social interaction among shrimp and depending on its intensity, it can affect heritable variation and response to selection. Little is known about the variance of indirect genetic effects (IGE) under competitive and non-competitive conditions in shrimp. In this study, we used extended mixed linear models to estimate genetic parameters for the direct genetic effect (DGE) and IGE on body weight in Litopenaeus vannamei raised under ad libitum (AF, non-competitive environment) and restricted (RF, competitive environment) feeding regimes.

Results

Estimates of heritabilities for body weight obtained with a traditional animal model (i.e. without accounting for IGE) were 0.11 ± 0.09 under AF and 0.25 ± 0.11 under RF. With extended animal models that accounted for IGE, the corresponding estimates for body weight were 0.07 ± 0.08 and 0.34 ± 0.11. Thus, heritabilities were higher under the RF regime than under the AF regime, regardless of whether IGE was accounted for or not. The log-likelihood ratio test revealed significant IGE under the RF regime. Although estimates of indirect genetic variance were low (0.0023 ± 0.0013 for AF and 0.0028 ± 0.0012 for RF), they contributed substantially to the total heritable variance: 66.8% for AF and 692.2% for RF. The total heritable variance was smaller under the RF regime (0.7 ± 1.3) than under the AF regime (5.8 ± 2.6) because of the high contribution of the negative covariance between DGE and IGE (− 7.03). Estimates of the correlation between DGE and IGE were 0.32 ± 0.47 under AF and − 0.93 ± 0.15 under RF, those of DGE and IGE for body weight between both regimes were 0.94 ± 0.07 and 0.67 ± 0.20, respectively, and those of IGE for body weight with DGE for survival were − 0.12 ± 0.22 under AF and − 0.58 ± 0.20 under RF.

Conclusions

These results indicate that strong competitive interactions occurred under the RF regime in L. vannamei. Significant reranking and variation in IGE of individuals were observed between the two feeding regimes. Strong competitive interactions reduced the total heritable variation for body weight when food was restricted. These results indicate that the extent of competition among L. vannamei depends on the feeding regime applied and that this competition affects the genetic basis of body weight.

Genomic regions influencing aggressive behavior in honey bees are defined by colony allele frequencies

For social animals, the genotypes of group members affect the social environment, and thus individual behavior, often indirectly. We used genome-wide association studies (GWAS) to determine the influence of individual vs. group genotypes on aggression in honey bees. Aggression in honey bees arises from the coordinated actions of colony members, primarily nonreproductive "soldier" bees, and thus, experiences evolutionary selection at the colony level. Here, we show that individual behavior is influenced by colony environment, which in turn, is shaped by allele frequency within colonies. Using a population with a range of aggression, we sequenced individual whole genomes and looked for genotype-behavior associations within colonies in a common environment. There were no significant correlations between individual aggression and specific alleles. By contrast, we found strong correlations between colony aggression and the frequencies of specific alleles within colonies, despite a small number of colonies. Associations at the colony level were highly significant and were very similar among both soldiers and foragers, but they covaried with one another. One strongly significant association peak, containing an ortholog of the Drosophila sensory gene dpr4 on linkage group (chromosome) 7, showed strong signals of both selection and admixture during the evolution of gentleness in a honey bee population. We thus found links between colony genetics and group behavior and also, molecular evidence for group-level selection, acting at the colony level. We conclude that group genetics dominates individual genetics in determining the fatal decision of honey bees to sting.

Genetic control of non-genetic inheritance in mammals: state-of-the-art and perspectives

Thought to be directly and uniquely dependent from genotypes, the ontogeny of individual phenotypes is much more complicated. Individual genetics, environmental exposures, and their interaction are the three main determinants of individual's phenotype. This picture has been further complicated a decade ago when the Lamarckian theory of acquired inheritance has been rekindled with the discovery of epigenetic inheritance, according to which acquired phenotypes can be transmitted through fertilization and affect phenotypes across generations. The results of Genome-Wide Association Studies have also highlighted a big degree of missing heritability in genetics and have provided hints that not only acquired phenotypes, but also individual's genotypes affect phenotypes intergenerationally through indirect genetic effects. Here, we review available examples of indirect genetic effects in mammals, what is known of the underlying molecular mechanisms and their potential impact for our understanding of missing heritability, phenotypic variation. and individual disease risk.