Quantitative genetic methods depending on the nature of the phenotypic trait

Assessing the impact of pedigree attributes on the validity of quantitative genetic parameter estimates

Investigating the evolution of complex traits in nature requires accurate assessment of their genetic basis. Quantitative genetic (QG) modeling is frequently applied to estimate the additive genetic variance (VA) in traits, combining phenotypic and pedigree data from a sample of individuals. Whether reconstructed from social links or molecular markers, empirical pedigrees differ in completeness, genealogical error rates, and other attributes that can impact QG estimation. Here we investigate this impact using human genealogical data for 6 French-Canadian (FC) populations originating from the same genetic founding source but differing in their pedigrees' attributes. First, we simulated phenotypic values along pedigrees and under different trait architectures and "true" parameter values (e.g., VA). Then we fitted mixed effects "animal" models to these simulated data, to assess how QG estimation was impacted by pedigree attributes. Our results show that pedigree size and depth were important determinants of the precision, but not accuracy, of genetic parameter estimates. In contrast, pedigree completeness and entropy, 2 attributes related to the density of genealogical links, were not clearly associated with the performance of parameter estimation. Noticeably, a slight increase in the genealogical error rate was sufficient to cause a detectable underestimation of VA. Including maternal genetic effects in the simulations led to a slight underestimation of VA with pedigrees of smaller size and depth. Despite originating from the same genetic source, the 6 pedigrees yielded wide variations in QG estimates under identical conditions. These findings highlight the importance of sensitivity analyses in pedigree-based genetic studies on natural populations.

Who Keeps the House after Divorcing? Partner and Nest Fidelity in the Long-Lived Alpine Swift

AbstractDespite the advantages of lasting pair bonds and the prevalence of monogamy, at least in avian species, some individuals switch mates (divorce). Divorce is generally considered to be adaptive (i.e., conferring net fitness benefits), although its causes and consequences often remain unclear, most notably regarding the genetic basis of this behavior. Using more than 30 years of data in a long-lived bird with obligate biparental care, the Alpine swift, we first described the overall patterns of mate and nest site fidelity and investigated the predictors of between-year divorce. We show that 16.6% of pairings ended in divorce, with low reproductive success and young age as predictors of divorce, and that males retained the nest site more often than females. By then studying individual repeatability and heritability of divorce, we show moderate repeatability in females and low repeatability in males and little additive genetic variance in either sex. Finally, we assessed the fitness consequences of divorce and report that an active decision to modify the pair bond (divorce) may be more beneficial than reactionary re-pairing following a partner's death. Overall, divorce may provide some reproductive benefits for Alpine swifts, but no microevolutionary potential of this behavior is evident in this population.

Neutral drift upon threshold-like selection promotes variation in antibiotic resistance phenotype

Heritable phenotypic variation plays a central role in evolution by conferring rapid adaptive capacity to populations. Mechanisms that can explain genetic diversity by describing connections between genotype and organismal fitness have been described. However, the difficulty of acquiring comprehensive data on genotype-phenotype-environment relationships has hindered the efforts to explain how the ubiquitously observed phenotypic variation in populations emerges and is maintained. To address this challenge, we establish an experimental system where we can examine the genotype-phenotype relationships in a controlled environment. We perform long-term experimental evolution on VIM-2 β-lactamase, an antibiotic-resistance enzyme, to explore the conditions that promote the emergence and maintenance of phenotypic variation. We found that evolution in a static environment with low antibiotic concentrations can promote and maintain significant phenotypic variation within populations. Notably, evolution of VIM-2 under selection with a low antibiotic concentration led to variants that conferred resistance to over 100-fold higher antibiotic concentrations than used in selection. A model based on the previously described threshold-like relationship between enzyme phenotype and fitness generated using VIM-2's all single amino acid variants, sufficiently explains the emergence of standing phenotypic variation under static environmental conditions. Overall, our approach provides a tractable model for studying phenotypic variation and evolvability at the population level.

Benchmarking performance in Targhee sheep in development of a genetic reference flock

The Targhee breed is important to range sheep production in the Western United States. The objective of this research was to integrate industry sires participating in national genetic evaluation through the National Sheep Improvement Program (NSIP) into the U.S. Sheep Experiment Station (USSES) flock, where the breed originated, to benchmark sire performance. Estimated breeding values (EBV) of industry sires (n = 16) and USSES sires (n = 12) from the 2023 NSIP Targhee genetic evaluation differed (P < 0.05) only for the NSIP Number Born EBV. This difference, and small (nonsignificant) differences favoring industry sires for maternal weaning weight, yearling fiber diameter, and yearling staple length EBV also resulted in greater (P < 0.05) Western Range Index scores for industry sires. The performance of sires' direct progeny born 2016 to 2018 was compared for 17 traits. After data cleaning, there were 664 progeny (417 from industry sires and 247 from USSES sires) with 146 to 664 records per trait. Least squares means differed (P < 0.05) for grease fleece weight, side fiber diameter, and britch fiber diameter in favor of the industry sires. Production traits were measured from the progeny and other retained descendants of both sire genetic groups from 2016 through 2022. After data cleaning, there was a range of 1,138 to 1,493 records per trait. A three-generation pedigree was constructed and the proportional assignment of each lamb to each genetic group (e.g., 0.5 industry, 0.5 USSES) was included in an augmented relationship matrix. A univariate animal model was fitted for each lamb trait and a repeated measures model for each of the three ewe traits was fitted to estimate variance components and predict breeding values. Genetic group solutions did not differ (P > 0.05) for any of these traits but provided a means by which to compare the two groups. Establishment of a Targhee genetic reference flock at USSES has been initiated, which can be used to address issues of importance to the industry.

Reproductive aging weakens offspring survival and constrains the telomerase response to herpesvirus in Pacific oysters

Telomere length (TL) is increasingly recognized as a molecular marker that reflects how reproductive aging affects intergenerational transmissions. Here, we investigated the effects of parental age on offspring survival and the regulation of TL by examining the telomere-elongating activity of telomerase in the Pacific oyster. We assessed the classical hallmarks of aging in parents at three age classes (young, middle-aged, and old) and crossbred them using a split-brood design to examine the consequences of the nine maternal-by-paternal age combinations on their offspring. Reproductive aging leads to increased larval mortality and accelerated telomere shortening in spats, rendering them more susceptible to infection by the Ostreid herpesvirus. Viral exposure stimulates telomerase activity, a response that we identified as adaptive, but weakened by parental aging. While telomerase lengthens a spat's telomere, paradoxically, longer individual TL predicts higher mortality in adults. The telomerase-telomere complex appeared as a conservative biomarker for distinguishing survivors and losers upon exposure to polymicrobial diseases.

The best linear unbiased prediction (BLUP) method as a tool to estimate the lifetime risk of pancreatic ductal adenocarcinoma in high-risk individuals with no known pathogenic germline variants

Pancreatic ductal adenocarcinoma (PDAC) is the fourth leading cause of cancer-related death in the Western world. The number of diagnosed cases and the mortality rate are almost equal as the majority of patients present with advanced disease at diagnosis. Between 4 and 10% of pancreatic cancer cases have an apparent hereditary background, known as hereditary pancreatic cancer (HPC) and familial pancreatic cancer (FPC), when the genetic basis is unknown. Surveillance of high-risk individuals (HRI) from these families by imaging aims to detect PDAC at an early stage to improve prognosis. However, the genetic basis is unknown in the majority of HRIs, with only around 10-13% of families carrying known pathogenic germline mutations. The aim of this study was to assess an individual's genetic cancer risk based on sex and personal and family history of cancer. The Best Linear Unbiased Prediction (BLUP) methodology was used to estimate an individual's predicted risk of developing cancer during their lifetime. The model uses different demographic factors in order to estimate heritability. A reliable estimation of heritability for pancreatic cancer of 0.27 on the liability scale, and 0.07 at the observed data scale as obtained, which is different from zero, indicating a polygenic inheritance pattern of PDAC. BLUP was able to correctly discriminate PDAC cases from healthy individuals and those with other cancer types. Thus, providing an additional tool to assess PDAC risk HRI with an assumed genetic predisposition in the absence of known pathogenic germline mutations.

Limits to the adaptation of herbivorous spider mites to metal accumulation in homogeneous and heterogeneous environments

Metal accumulation is used by some plants as a defence against herbivores. Yet, herbivores may adapt to these defences, becoming less susceptible. Moreover, ecosystems often contain plants that do and do not accumulate metals, but whether such heterogeneity affects herbivore adaptation remains understudied. Here, we performed experimental evolution to test whether the spider mite Tetranychus evansi adapts to plants with high cadmium concentrations, in homogeneous (plants with cadmium) or heterogeneous (plants with or without cadmium) environments. For that we used tomato plants, which accumulate cadmium, thus affecting the performance of these spider mites. We measured mite fecundity, hatching rate, and the number of adult offspring after 12 and 33 generations and habitat choice after 14 and 51 generations, detecting no trait change, which implies the absence of adaptation. We then tested whether this was due to a lack of genetic variation in the traits measured and, indeed, additive genetic variance was low. Interestingly, despite no signs of adaptation, we observed a decrease in fecundity and number of adult offspring produced on cadmium-free plants, in the populations evolving in environments with cadmium. Therefore, evolving in environments with cadmium reduces the growth rate of spider mite populations on non-accumulating plants. Possibly, other traits contributed to population persistence on plants with cadmium. This calls for more studies addressing herbivore adaptation to plant metal accumulation.

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.

Genetic correlations of direct and indirect genetic components of social dominance with fitness and morphology traits in cattle

BACKGROUND

Within the same species, individuals show marked variation in their social dominance. Studies on a handful of populations have indicated heritable genetic variation for this trait, which is determined by both the genetic background of the individual (direct genetic effect) and of its opponent (indirect genetic effect). However, the evolutionary consequences of selection for this trait are largely speculative, as it is not a usual target of selection in livestock populations. Moreover, studying social dominance presents the challenge of working with a phenotype with a mean value that cannot change in the population, as for every winner of an agonistic interaction there will necessarily be a loser. Thus, to investigate what could be the evolutionary response to selection for social dominance, it is necessary to focus on traits that might be correlated with it. This study investigated the genetic correlations of social dominance, both direct and indirect, with several morphology and fitness traits. We used a dataset of agonistic contests involving cattle (Bos taurus): during these contests, pairs of cows compete in ritualized interactions to assess social dominance. The outcomes of 37,996 dominance interactions performed by 8789 cows over 20 years were combined with individual data for fertility, mammary health, milk yield and morphology and analysed using bivariate animal models including indirect genetic effects.

RESULTS

We found that winning agonistic interactions has a positive genetic correlation with more developed frontal muscle mass, lower fertility, and poorer udder health. We also discovered that the trends of changes in the estimated breeding values of social dominance, udder health and more developed muscle mass were consistent with selection for social dominance in the population.

CONCLUSIONS

We present evidence that social dominance is genetically correlated with fitness traits, as well as empirical evidence of the possible evolutionary trade-offs between these traits. We show that it is feasible to estimate genetic correlations involving dyadic social traits.

Visually assessed body condition shows high heritability in a pedigreed great ape population

Body condition, a measure for relative fat mass, is associated with primate health, fitness, and overall welfare. Body condition is often influenced by dietary factors, age, and/or sex, but several body condition measures (body weight, weight-to-height ratios, and so on) also show high heritability across primate species, indicating a role of genetic effects. Although different measures for body condition exist, many require direct handling of animals, which is invasive, time-consuming, and expensive, making them impractical in wild and captive settings. Therefore, noninvasive visual body condition score (BCS) systems were developed for various animal species, including macaques and chimpanzees, to visually assess relative fat mass. Here we evaluate the utility of a visual BCS system in bonobos by assessing (1) inter-rater reliability, (2) links with body mass, a traditional hands-on measure of condition, and (3) the factors driving individual variation in BCS. We adapted the chimpanzee BCS system to rate 76 bonobos in 11 European zoos (92% of the adult population). Inter-rater reliability was high (s* = 0.948), BCSs were positively associated with body mass (β = 0.075) and not predicted by diet, sex, or age, nor were they associated with a higher abundance of obesity-related diseases. Instead, BCSs showed high levels of heritability (h2  = 0.637), indicating that a majority of body condition variation in bonobos is attributable to genetic similarity of the individuals. This is in line with reported h2 -values for traditional body condition measures in primates and provides support for the reliability of visual BCS systems in great apes. The results of this study emphasize an often unanticipated role of genetics in determining primate body fat and health that has implications for the management of captive primates. Application of this tool in wild populations would aid to unravel environmental from genetic drivers of body condition variation in primates.

An investigation of the sex-specific genetic architecture of fitness in Drosophila melanogaster

In dioecious populations, the sexes employ divergent reproductive strategies to maximize fitness and, as a result, genetic variants can affect fitness differently in males and females. Moreover, recent studies have highlighted an important role of the mating environment in shaping the strength and direction of sex-specific selection. Here, we measure adult fitness for each sex of 357 lines from the Drosophila Synthetic Population Resource in two different mating environments. We analyze the data using three different approaches to gain insight into the sex-specific genetic architecture for fitness: classical quantitative genetics, genomic associations, and a mutational burden approach. The quantitative genetics analysis finds that on average segregating genetic variation in this population has concordant fitness effects both across the sexes and across mating environments. We do not find specific genomic regions with strong associations with either sexually antagonistic (SA) or sexually concordant (SC) fitness effects, yet there is modest evidence of an excess of genomic regions with weak associations, with both SA and SC fitness effects. Our examination of mutational burden indicates stronger selection against indels and loss-of-function variants in females than in males.

Same-sex sociosexual behaviour is widespread and heritable in male rhesus macaques

Numerous reports have documented the occurrence of same-sex sociosexual behaviour (SSB) across animal species. However, the distribution of the behaviour within a species needs to be studied to test hypotheses describing its evolution and maintenance, in particular whether the behaviour is heritable and can therefore evolve by natural selection. Here we collected detailed observations across 3 yr of social and mounting behaviour of 236 male semi-wild rhesus macaques, which we combined with a pedigree dating back to 1938, to show that SSB is both repeatable (19.35%) and heritable (6.4%). Demographic factors (age and group structure) explained SSB variation only marginally. Furthermore, we found a positive genetic correlation between same-sex mounter and mountee activities, indicating a common basis to different forms of SSB. Finally, we found no evidence of fitness costs to SSB, but show instead that the behaviour mediated coalitionary partnerships that have been linked to improved reproductive success. Together, our results demonstrate that SSB is frequent in rhesus macaques, can evolve, and is not costly, indicating that SSB may be a common feature of primate reproductive ecology.

Artificial selection for reversal learning reveals limited repeatability and no heritability of cognitive flexibility in great tits (Parus major)

Cognitive flexibility controls how animals respond to changing environmental conditions. Individuals within species vary considerably in cognitive flexibility but the micro-evolutionary potential in animal populations remains enigmatic. One prerequisite for cognitive flexibility to be able to evolve is consistent and heritable among-individual variation. Here we determine the repeatability and heritability of cognitive flexibility among great tits (Parus major) by performing an artificial selection experiment on reversal learning performance using a spatial learning paradigm over three generations. We found low, yet significant, repeatability (R = 0.15) of reversal learning performance. Our artificial selection experiment showed no evidence for narrow-sense heritability of associative or reversal learning, while we confirmed the heritability of exploratory behaviour. We observed a phenotypic, but no genetic, correlation between associative and reversal learning, showing the importance of prior information on reversal learning. We found no correlation between cognitive and personality traits. Our findings emphasize that cognitive flexibility is a multi-faceted trait that is affected by memory and prior experience, making it challenging to retrieve reliable values of temporal consistency and assess the contribution of additive genetic variation. Future studies need to identify what cognitive components underlie variation in reversal learning and study their between-individual and additive genetic components.

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

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

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

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

Genetic correlation between aggressive signals and fighting

Theoretical analyses indicate that aggressive signals should positively correlate with the signallers' willingness and abilities to fight. Few experimental studies, however, have tested this prediction. In two experiments employing distinct, ecologically realistic protocols, we quantified the association between aggressive signals and fighting in fruit fly genotypes and found high positive genetic correlations between threat and fighting (rG = 0.80 and 0.74). Our results add to the growing body of experimental work indicating that aggressive signals have relatively high informational value.

Terpene chemotypes in Gossypium hirsutum (wild cotton) from the Yucatan Peninsula, Mexico

Cultivated plants of Gossypium hirsutum Cav. (cotton) consistently emit low levels of volatile organic compounds, primarily mono- and sesquiterpenoids, which are produced and stored in pigment glands. In this study, we provide a comprehensive evaluation of the terpene profiles of wild G. hirsutum plants sourced from sites located throughout natural distribution of this species, thus providing the first in-depth assessment of the scope of its intraspecific chemotypic diversity. Chemotypic variation can potentially influence resistance to herbivory and diseases, or interact with abiotic stress such as extreme temperatures. Under controlled environmental conditions, plants were grown from seeds of sixteen G. hirsutum populations collected along the coastline of the Yucatan Peninsula, which is its likely centre of origin. We found high levels of intraspecific diversity in the terpene profiles of the plants. Two distinct chemotypes were identified: one chemotype contained higher levels of the monoterpenes γ-terpinene, limonene, α-thujene, α-terpinene, terpinolene, and p-cymene, while the other chemotype was distinguished by higher levels of α- and β-pinene. The distribution of chemotypes followed a geographic gradient from west to east, with an increasing frequency of the former chemotype. Concurrent analysis of maternal plants revealed that chemotypes in wild G. hirsutum are highly heritable.

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.

Maternal effects do not resolve the paradox of stasis in birth weight in a wild red deer populaton

In natural populations, quantitative traits seldom show short-term evolution at the rate predicted by evolutionary models. Resolving this "paradox of stasis" is a key goal in evolutionary biology, as it directly challenges our capacity to predict evolutionary change. One particularly promising hypothesis to explain the lack of evolutionary responses in a key offspring trait, body weight, is that positive selection on juveniles is counterbalanced by selection against maternal investment in offspring growth, given that reproduction is costly for the mothers. Here, we used data from one of the longest individual-based studies of a wild mammal population to test this hypothesis. We first showed that despite positive directional selection on birth weight, and heritable variation for this trait, no genetic change has been observed for birth weight over the past 47 years in the study population. Contrarily to our expectation, we found no evidence of selection against maternal investment in birth weight-if anything, selection favors mothers that produce large calves. Accordingly, we show that genetic change in birth weight over the study period is actually lower than that predicted from models including selection on maternal performance; ultimately our analysis here only deepens rather than resolves the paradox of stasis.

Genetic, maternal, and environmental influences on sociality in a pedigreed primate population

Various aspects of sociality in mammals (e.g., dyadic connectedness) are linked with measures of biological fitness (e.g., longevity). How within- and between-individual variation in relevant social traits arises in uncontrolled wild populations is challenging to determine but is crucial for understanding constraints on the evolution of sociality. We use an advanced statistical method, known as the 'animal model', which incorporates pedigree information, to look at social, genetic, and environmental influences on sociality in a long-lived wild primate. We leverage a longitudinal database spanning 20 years of observation on individually recognized white-faced capuchin monkeys (Cebus capucinus imitator), with a multi-generational pedigree. We analyze two measures of spatial association, using repeat sampling of 376 individuals (mean: 53.5 months per subject, range: 6-185 months per subject). Conditioned on the effects of age, sex, group size, seasonality, and El Niño-Southern Oscillation phases, we show low to moderate long-term repeatability (across years) of the proportion of time spent social (posterior mode [95% Highest Posterior Density interval]: 0.207 [0.169, 0.265]) and of average number of partners (0.144 [0.113, 0.181]) (latent scale). Most of this long-term repeatability could be explained by modest heritability (h²_social: 0.152 [0.094, 0.207]; h²_partners: 0.113 [0.076, 0.149]) with small long-term maternal effects (m²_social: 0.000 [0.000, 0.045]; m²_partners: 0.000 [0.000, 0.041]). Our models capture the majority of variance in our behavioral traits, with much of the variance explained by temporally changing factors, such as group of residence, highlighting potential limits to the evolvability of our trait due to social and environmental constraints.

Intrinsic emergence and modulation of sex-specific dominance reversals in threshold traits

Sex-specific dominance reversals (SSDRs) in fitness-related traits, where heterozygotes' phenotypes resemble those of alternative homozygotes in females versus males, can simultaneously maintain genetic variation in fitness and resolve sexual conflict and thereby shape key evolutionary outcomes. However, the full implications of SSDRs will depend on how they arise and the resulting potential for evolutionary, ecological and environmental modulation. Recent field and laboratory studies have demonstrated SSDRs in threshold(-like) traits with dichotomous or competitive phenotypic outcomes, implying that such traits could promote the emergence of SSDRs. However, such possibilities have not been explicitly examined. I show how phenotypic SSDRs can readily emerge in threshold traits given genetic architectures involving large-effect loci alongside sexual dimorphism in the mean and variance in polygenic liability. I also show how multilocus SSDRs can arise in line-cross experiments, especially given competitive reproductive systems that generate nonlinear fitness outcomes. SSDRs can consequently emerge in threshold(-like) traits as functions of sexual antagonism, sexual dimorphism and reproductive systems, even with purely additive underlying genetic effects. Accordingly, I identify theoretical and empirical advances that are now required to discern the basis and occurrence of SSDRs in nature, probe forms of (co-)evolutionary, ecological and environmental modulation, and evaluate net impacts on sexual conflict.

Reproductive tract size and position score: Estimation of genetic parameters for a novel fertility trait in dairy cows

The dairy industry is moving toward selecting animals with better fertility to decrease the economic losses linked to reproductive issues. The reproductive tract size and position score (SPS) was recently developed in physiological studies as an indicator of pregnancy rate and the number of services to conception. Cows are scored as SPS 1, 2, or 3 based on the size of their reproductive tract and its position in the pelvis, as determined by transrectal palpation. The objective of this study was to estimate genetic parameters for SPS to assess its potential as a novel fertility trait. Phenotypes were collected at the University of British Columbia's research herd from 2017 to 2020, consisting of 3,247 within- and across-lactation SPS records from 490 Holstein cows. A univariate animal model was used to estimate the variance components for SPS. Both threshold and linear models were fit under a Bayesian approach and the results were compared using the Spearman rank correlation (r) between the estimated breeding values. The 2 models ranked the animals very similarly (r = 0.99), and the linear model was selected for further analysis. Genetic correlations with other currently evaluated traits were estimated using a bivariate animal model. The posterior means (± posterior standard deviation) for heritability and repeatability within- and across-lactation were 0.113 (± 0.013), 0.242 (± 0.012), and 0.134 (± 0.014), respectively. The SPS showed null correlations with production traits and favorable correlations with traditional fertility traits, varying from -0.730 (nonreturn rate) to 0.931 (number of services). Although preliminary, these results are encouraging because SPS seems to be more heritable than and strongly genetically correlated with number of services, nonreturn rate, and first service to conception, indicating potential for effective indirect selection response on these traits from SPS genetic selection. Therefore, further studies with larger data sets to validate these findings are warranted.

The quantitative genetics of fitness in a wild seabird

Additive genetic variance in fitness is a prerequisite for adaptive evolution, as a trait must be genetically correlated with fitness to evolve. Despite its relevance, additive genetic variance in fitness has not often been estimated in nature. Here, we investigate additive genetic variance in lifetime and annual fitness components in common terns (Sterna hirundo). Using 28 years of data comprising approximately 6000 pedigreed individuals, we find that additive genetic variances in the zero-inflated and Poisson components of lifetime fitness were effectively zero but estimated with high uncertainty. Similarly, additive genetic variances in adult annual reproductive success and survival did not differ from zero but were again associated with high uncertainty. Simulations suggested that we would be able to detect additive genetic variances as low as 0.05 for the zero-inflated component of fitness but not for the Poisson component, for which adequate statistical power would require approximately two more decades (four tern generations) of data collection. As such, our study suggests heritable variance in common tern fitness to be rather low if not zero, shows how studying the quantitative genetics of fitness in natural populations remains challenging, and highlights the importance of maintaining long-term individual-based studies of natural populations.

Interpreting Mendelian-randomization estimates of the effects of categorical exposures such as disease status and educational attainment

BACKGROUND

Mendelian randomization has been previously used to estimate the effects of binary and ordinal categorical exposures-e.g. Type 2 diabetes or educational attainment defined by qualification-on outcomes. Binary and categorical phenotypes can be modelled in terms of liability-an underlying latent continuous variable with liability thresholds separating individuals into categories. Genetic variants influence an individual's categorical exposure via their effects on liability, thus Mendelian-randomization analyses with categorical exposures will capture effects of liability that act independently of exposure category.

METHODS AND RESULTS

We discuss how groups in which the categorical exposure is invariant can be used to detect liability effects acting independently of exposure category. For example, associations between an adult educational-attainment polygenic score (PGS) and body mass index measured before the minimum school leaving age (e.g. age 10 years), cannot indicate the effects of years in full-time education on this outcome. Using UK Biobank data, we show that a higher educational-attainment PGS is strongly associated with lower smoking initiation and higher odds of glasses use at age 15 years. These associations were replicated in sibling models. An orthogonal approach using the raising of the school leaving age (ROSLA) policy change found that individuals who chose to remain in education to age 16 years before the reform likely had higher liability to educational attainment than those who were compelled to remain in education to age 16 years after the reform, and had higher income, lower pack-years of smoking, higher odds of glasses use and lower deprivation in adulthood. These results suggest that liability to educational attainment is associated with health and social outcomes independently of years in full-time education.

CONCLUSIONS

Mendelian-randomization studies with non-continuous exposures should be interpreted in terms of liability, which may affect the outcome via changes in exposure category and/or independently.

Age-dependent genetic architecture across ontogeny of body size in sticklebacks

Heritable variation in traits under natural selection is a prerequisite for evolutionary response. While it is recognized that trait heritability may vary spatially and temporally depending on which environmental conditions traits are expressed under, less is known about the possibility that genetic variance contributing to the expected selection response in a given trait may vary at different stages of ontogeny. Specifically, whether different loci underlie the expression of a trait throughout development and thus providing an additional source of variation for selection to act on in the wild, is unclear. Here we show that body size, an important life-history trait, is heritable throughout ontogeny in the nine-spined stickleback (Pungitius pungitius). Nevertheless, both analyses of quantitative trait loci and genetic correlations across ages show that different chromosomes/loci contribute to this heritability in different ontogenic time-points. This suggests that body size can respond to selection at different stages of ontogeny but that this response is determined by different loci at different points of development. Hence, our study provides important results regarding our understanding of the genetics of ontogeny and opens an interesting avenue of research for studying age-specific genetic architecture as a source of non-parallel evolution.

Identifying Developmental Patterns in Structured Plant Phenotyping Data

Technological breakthroughs concerning both sensors and robotized plant phenotyping platforms have totally renewed the plant phenotyping paradigm in the last two decades. This has impacted both the nature and the throughput of data with the availability of data at high-throughput from the tissular to the whole plant scale. Sensor outputs often take the form of 2D or 3D images or time series of such images from which traits are extracted while organ shapes, shoot or root system architectures can be deduced. Despite this change of paradigm, many phenotyping studies often ignore the structure of the plant and therefore loose the information conveyed by the temporal and spatial patterns emerging from this structure. The developmental patterns of plants often take the form of succession of well-differentiated phases, stages or zones depending on the temporal, spatial or topological indexing of data. This entails the use of hierarchical statistical models for their identification.The objective here is to show potential approaches for analyzing structured plant phenotyping data using state-of-the-art methods combining probabilistic modeling, statistical inference and pattern recognition. This approach is illustrated using five different examples at various scales that combine temporal and topological index parameters, and development and growth variables obtained using prospective or retrospective measurements.

Additive genetic variance for traits least related to fitness increases with environmental stress in the desert locust, Schistocerca gregaria

Under environmental stress, previously hidden additive genetic variation can be unmasked and exposed to selection. The amount of hidden variation is expected to be higher for life history traits, which strongly correlate to individual fitness, than for morphological traits, in which fitness effects are more ambiguous. However, no consensual pattern has been recovered yet, and this idea is still debated in the literature. Here, we hypothesize that the classical categorization of traits (i.e., life history and morphology) may fail to capture their proximity to fitness. In the desert locust, Schistocerca gregaria, a model organism for the study of insect polyphenism, we quantified changes in additive genetic variation elicited by lifetime thermal stress for ten traits, in which evolutionary significance is known. Irrespective of their category, traits under strong stabilizing selection showed genetic invariance with environmental stress, while traits more loosely associated with fitness showed a marked increase in additive genetic variation in the stressful environment. Furthermore, traits involved in adaptive phenotypic plasticity (growth compensation) showed either no change in additive genetic variance or a change of moderate magnitude across thermal environments. We interpret this mitigated response of plastic traits in the context of integrated evolution to adjust the entire phenotype in heterogeneous environments (i.e., adaptiveness of initial plasticity, compromise of phenotypic compensation with stress, and shared developmental pathway). Altogether, our results indicate, in agreement with theoretical expectations, that environmental stress can increase available additive genetic variance in some desert locust traits, but those closely linked to fitness are largely unaffected. Our study also highlights the importance of assessing the proximity to fitness of a trait on a case-by-case basis and in an ecologically relevant context, as well as considering the processes of canalization and plasticity, involved in the control of phenotypic variation.

Plasticity and artificial selection for developmental mode in a poecilogonous sea slug

The contribution of phenotypically plastic traits to evolution depends on the degree of environmental influence on the target of selection (the phenotype) as well as the underlying genetic structure of the trait and plastic response. Likewise, maternal effects can help or hinder evolution through affects to the response to selection. The sacoglossan sea slug Alderia willowi exhibits intraspecific variation for developmental mode (= poecilogony) that is environmentally modulated with populations producing more yolk-feeding (lecithotrophic) larvae during the summer, and more planktonic-feeding (planktotrophic) larvae in the winter. I found significant family-level variation in the reaction norms between 17 maternal families of A. willowi when reared in a split-brood design in low (16 ppt) versus high (32 ppt) salinity, conditions which mimic seasonal variation in salinity of natural populations. I documented a significant response to selection for lecithotrophic larvae in high and low salinity. The slope of the reaction norm was maintained following one generation of selection for lecithotrophy. When the maternal environment was controlled in the laboratory, I found significant maternal effects, which reduced the response to selection. These results suggest there is standing genetic variation for egg-mass type in A. willowi, but the ability of selection to act on that variation may depend on the environment in which the phenotype is expressed in preceding generations.

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.

How does contemporary selection shape oak phenotypes?

Most existing forests are subjected to natural and human-mediated selection pressures, which have increased due to climate change and the increasing needs of human societies for wood, fibre and fuel resources. It remains largely unknown how these pressures trigger evolutionary changes. We address this issue here for temperate European oaks (Quercus petraea and Q. robur), which grow in mixed stands, under even-aged management regimes. We screened numerous functional traits for univariate selection gradients and for expected and observed genetic changes over two successive generations. In both species, growth, leaf morphology and physiology, and defence-related traits displayed significant selection gradients and predicted shifts, whereas phenology, water metabolism, structure and resilience-related traits did not. However, the direction of the selection response and the potential for adaptive evolution differed between the two species. Quercus petraea had a much larger phenotypic and genetic variance of fitness than Q. robur. This difference raises concerns about the adaptive response of Q. robur to contemporary selection pressures. Our investigations suggest that Q. robur will probably decline steadily, particularly in mixed stands with Q. petraea, consistent with the contrasting demographic dynamics of the two species.

Genetic and nongenetic profiling of milk β-hydroxybutyrate and acetone and their associations with ketosis in Holstein cows

Ketosis is a metabolic disorder of increasing importance in high-yielding dairy cows, but accurate population-wide binary health trait recording is difficult to implement. Against this background, proper Gaussian indicator traits, which can be routinely measured in milk, are needed. Consequently, we focused on the ketone bodies acetone and β-hydroxybutyrate (BHB), measured via Fourier-transform infrared spectroscopy (FTIR) in milk. In the present study, 62,568 Holstein cows from large-scale German co-operator herds were phenotyped for clinical ketosis (KET) according to a veterinarian diagnosis key. A sub-sample of 16,861 cows additionally had first test-day observations for FTIR acetone and BHB. Associations between FTIR acetone and BHB with KET and with test-day traits were studied phenotypically and quantitative genetically. Furthermore, we estimated SNP marker effects for acetone and BHB (application of genome-wide association studies) based on 40,828 SNP markers from 4,384 genotyped cows, and studied potential candidate genes influencing body fat mobilization. Generalized linear mixed models were applied to infer the influence of binary KET on Gaussian-distributed acetone and BHB (definition of an identity link function), and vice versa, such as the influence of acetone and BHB on KET (definition of a logit link function). Additionally, linear models were applied to study associations between BHB, acetone and test-day traits (milk yield, fat percentage, protein percentage, fat-to-protein ratio and somatic cell score) from the first test-day after calving. An increasing KET incidence was statistically significant associated with increasing FTIR acetone and BHB milk concentrations. Acetone and BHB concentrations were positively associated with fat percentage, fat-to-protein ratio and somatic cell score. Bivariate linear animal models were applied to estimate genetic (co)variance components for KET, acetone, BHB and test-day traits within parities 1 to 3, and considering all parities simultaneously in repeatability models. Pedigree-based heritabilities were quite small (i.e., in the range from 0.01 in parity 3 to 0.07 in parity 1 for acetone, and from 0.03-0.04 for BHB). Heritabilites from repeatability models were 0.05 for acetone, and 0.03 for BHB. Genetic correlations between acetone and BHB were moderate to large within parities and considering all parities simultaneously (0.69-0.98). Genetic correlations between acetone and BHB with KET from different parities ranged from 0.71 to 0.99. Genetic correlations between acetone across parities, and between BHB across parities, ranged from 0.55 to 0.66. Genetic correlations between KET, acetone, and BHB with fat-to-protein ratio and with fat percentage were large and positive, but negative with milk yield. In genome-wide association studies, we identified SNP on BTA 4, 10, 11, and 29 significantly influencing acetone, and on BTA 1 and 16 significantly influencing BHB. The identified potential candidate genes NRXN3, ACOXL, BCL2L11, HIBADH, KCNJ1, and PRG4 are involved in lipid and glucose metabolism pathways.

In situ estimation of genetic variation of functional and ecological traits in Quercus petraea and Q.robur

BACKGROUND

Predicting the evolutionary potential of natural tree populations requires the estimation of heritability and genetic correlations among traits on which selection acts, as differences in evolutionary success between species may rely on differences for these genetic parameters. In situ estimates are expected to be more accurate than measures done under controlled conditions which do not reflect the natural environmental variance.

AIMS

The aim of the current study was to estimate three genetic parameters (i.e. heritability, evolvability and genetic correlations) in a natural mixed oak stand composed of Quercus petraea and Quercus robur about 100 years old, for 58 traits of ecological and functional relevance (growth, reproduction, phenology, physiology, resilience, structure, morphology and defence).

METHODS

First we estimated genetic parameters directly in situ using realized genomic relatedness of adult trees and parentage relationships over two generations to estimate the traits additive variance. Secondly, we benefited from existing ex situ experiments (progeny tests and conservation collection) installed with the same populations, thus allowing comparisons of in situ heritability estimates with more traditional methods.

RESULTS

Heritability and evolvability estimates obtained with different methods varied substantially and showed large confidence intervals, however we found that in situ were less precise than ex situ estimates, and assessments over two generations (with deeper relatedness) improved estimates of heritability while large sampling sizes are needed for accurate estimations. At the biological level, heritability values varied moderately across different ecological and functional categories of traits, and genetic correlations among traits were conserved over the two species.

CONCLUSION

We identified limits for using realized genomic relatedness in natural stands to estimate the genetic variance, given the overall low variance of genetic relatedness and the rather low sampling sizes of currently used long term genetic plots in forestry. These limits can be overcome if larger sample sizes are considered, or if the approach is extended over the next generation.

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.

Estimation of Genetic Variance in Fitness, and Inference of Adaptation, When Fitness Follows a Log-Normal Distribution

Additive genetic variance in relative fitness (σA2(w)) is arguably the most important evolutionary parameter in a population because, by Fisher’s fundamental theorem of natural selection (FTNS; Fisher RA. 1930. The genetical theory of natural selection. 1st ed. Oxford: Clarendon Press), it represents the rate of adaptive evolution. However, to date, there are few estimates of σA2(w) in natural populations. Moreover, most of the available estimates rely on Gaussian assumptions inappropriate for fitness data, with unclear consequences. “Generalized linear animal models” (GLAMs) tend to be more appropriate for fitness data, but they estimate parameters on a transformed (“latent”) scale that is not directly interpretable for inferences on the data scale. Here we exploit the latest theoretical developments to clarify how best to estimate quantitative genetic parameters for fitness. Specifically, we use computer simulations to confirm a recently developed analog of the FTNS in the case when expected fitness follows a log-normal distribution. In this situation, the additive genetic variance in absolute fitness on the latent log-scale (σA2(l)) equals (σA2(w)) on the data scale, which is the rate of adaptation within a generation. However, due to inheritance distortion, the change in mean relative fitness between generations exceeds σA2(l) and equals (exp⁡(σA2(l))−1). We illustrate why the heritability of fitness is generally low and is not a good measure of the rate of adaptation. Finally, we explore how well the relevant parameters can be estimated by animal models, comparing Gaussian models with Poisson GLAMs. Our results illustrate 1) the correspondence between quantitative genetics and population dynamics encapsulated in the FTNS and its log-normal-analog and 2) the appropriate interpretation of GLAM parameter estimates.

Alterations in Life-History Associated With Non-target-site Herbicide Resistance in Alopecurus myosuroides

The evolution of resistance to herbicides is a classic example of rapid contemporary adaptation in the face of a novel environmental stress. Evolutionary theory predicts that selection for resistance will be accompanied by fitness trade-offs in environments where the stress is absent. Alopecurus myosuroides, an autumn-germinating grass weed of cereal crops in North-West Europe, has evolved resistance to seven herbicide modes-of-action, making this an ideal species to examine the presence and magnitudes of such fitness costs. Here, we use two contrasting A. myosuroides phenotypes derived from a common genetic background, one with enhanced metabolism resistance to a commercial formulation of the sulfonylurea (ALS) actives mesosulfuron and iodosulfuron, and the other with susceptibility to these actives (S). Comparisons of plant establishment, growth, and reproductive potential were made under conditions of intraspecific competition, interspecific competition with wheat, and over a gradient of nitrogen deprivation. Herbicide dose response assays confirmed that the two lines had contrasting resistance phenotypes, with a 20-fold difference in resistance between them. Pleiotropic effects of resistance were observed during plant development, with R plants having a greater intraspecific competitive effect and longer tiller lengths than S plants during vegetative growth, but with S plants allocating proportionally more biomass to reproductive tissues during flowering. Direct evidence of a reproductive cost of resistance was evident in the nitrogen deprivation experiment with R plants producing 27% fewer seed heads per plant, and a corresponding 23% reduction in total seed head length. However, these direct effects of resistance on fecundity were not consistent across experiments. Our results demonstrate that a resistance phenotype based on enhanced herbicide metabolism has pleiotropic impacts on plant growth, development and resource partitioning but does not support the hypothesis that resistance is associated with a consistent reproductive fitness cost in this species. Given the continued difficulties associated with unequivocally detecting costs of herbicide resistance, we advocate future studies that adopt classical evolutionary quantitative genetics approaches to determine genetic correlations between resistance and fitness-related plant life history traits.