Multivariate genetic architecture reveals testosterone-driven sexual antagonism in contemporary humans.

Sex difference (SD) is ubiquitous in humans despite shared genetic architecture (SGA) between the sexes. A univariate approach, i.e., studying SD in single traits by estimating genetic correlation, does not provide a complete biological overview, because traits are not independent and are genetically correlated. The multivariate genetic architecture between the sexes can be summarized by estimating the additive genetic (co)variance across shared traits, which, apart from the cross-trait and cross-sex covariances, also includes the cross-sex-cross-trait covariances, e.g., between height in males and weight in females. Using such a multivariate approach, we investigated SD in the genetic architecture of 12 anthropometric, fat depositional, and sex-hormonal phenotypes. We uncovered sexual antagonism (SA) in the cross-sex-cross-trait covariances in humans, most prominently between testosterone and the anthropometric traits - a trend similar to phenotypic correlations. 27% of such cross-sex-cross-trait covariances were of opposite sign, contributing to asymmetry in the SGA. Intriguingly, using multivariate evolutionary simulations, we observed that the SGA acts as a genetic constraint to the evolution of SD in humans only when selection is sexually antagonistic and not concordant. Remarkably, we found that the lifetime reproductive success in both the sexes shows a positive genetic correlation with anthropometric traits, but not with testosterone. Moreover, we demonstrated that genetic variance is depleted along multivariate trait combinations in both the sexes but in different directions, suggesting absolute genetic constraint to evolution. Our results indicate that testosterone drives SA in contemporary humans and emphasize the necessity and significance of using a multivariate framework in studying SD.

Response to von Schmalensee et al.

Von Schmalensee et al. present two concerns about our study. While the first stems from a general disagreement about our simulation methodology, the second is a useful observation of a modelling choice we made that affected simulation outcomes, but in ways that do not invalidate our original conclusions.

The genomic consequences of selection across development.

Understanding how natural selection drives diversification in nature has been at the forefront of biological research for over a century. The main idea is simple: natural selection favours individuals best suited to pass on their genes. However, the journey from birth to reproduction is complex as organisms experience multiple developmental stages, each influenced by genetic and environmental factors (Orr, 2009). These complexities compound even further as each stage of development might be governed by a unique underlying set of alleles and genes. In this issue of Molecular Ecology, Goebl et al. (2022) examine the role of natural selection in driving ecotypic divergence across different life history stages of the prairie sunflower Helianthus petiolaris. The authors used reciprocal transplant experiments, demographic models, and genomic sequencing to explore fitness variation across developmental stages. They show how natural selection impacts population divergence across multiple life history stages and evaluate the resulting allele frequency changes. Goebl et al. link these results to the role of chromosomal inversions, thus furthering our understanding of how ecological divergence proceeds in the face of gene flow. Below, we explore these results in detail and complement their interpretation by considering the evolution of genetic correlations amongst traits governing fitness.

Diurnal variation in genetic parameters for locomotor activity in Drosophila melanogaster assessed under natural thermal conditions.

In nature, organisms are exposed to variable and occasionally stressful environmental conditions. Responses to diurnal and seasonal fluctuations, such as temperature and food accessibility, involve adaptive behavioural and physiological changes. While much work has been done on understanding the genetic architecture and evolutionary potential of stress tolerance traits under constant thermal conditions, there has been less focus on the quantitative genetic background in variable environments. In this study, we use the Drosophila Genetic Reference Panel (DGRP) to investigate the locomotor activity, a key behavioural trait, under variable natural thermal conditions during the summer in a temperate environment. Male flies from 100 DGRP lines were exposed to natural thermal and light conditions in Drosophila activity monitors across three experimental days. We found that activity was highly temperature and time dependent and varied between lines both within and between days. Furthermore, we observed variation in genetic and environmental variance components, with low to moderate estimates of the heritability for locomotor activity, consistently peaking in the afternoons. Moreover, we showed that the estimated genetic correlations of locomotor activity between two time points decreased, as the absolute differences in ambient temperature increased. In conclusion, we find that the genetic background for locomotor activity is environment specific, and we conclude that more variable and unpredictable future temperatures will likely have a strong impact on the evolutionary trajectories of behavioural traits in ectotherms.

The fecal microbiota of Holstein cows is heritable and genetically correlated to dairy performances.

The fecal microbiota of ruminants constitutes a diversified community that has been phenotypically associated with a variety of host phenotypes, such as production and health. To gain a better understanding of the complex and interconnected factors that drive the fecal bacterial community, we have aimed to estimate the genetic parameters of the diversity and composition of the fecal microbiota, including heritabilities, genetic correlations among taxa, and genetic correlations between fecal microbiota features and host phenotypes. To achieve this, we analyzed a large population of 1,875 Holstein cows originating from 144 French commercial herds and routinely recorded for production, somatic cell score, and fertility traits. Fecal samples were collected from the animals and subjected to 16S rRNA gene sequencing, with reads classified into Amplicon Sequence Variants (ASVs). The estimated α- and ß-diversity indices (i.e., Observed Richness, Shannon index, Bray-Curtis and Jaccard dissimilarity matrices) and the abundances of ASVs, genera, families and phyla, normalized by centered-log ratio (CLR), were considered as phenotypes. Genetic parameters were calculated using either univariate or bivariate animal models. Heritabilities estimates, ranging from 0.08 to 0.31 for taxa abundances and ß-diversity indices, highlight the influence of the host genetics on the composition of the fecal microbiota. Furthermore, genetic correlations estimated within the microbial community and between microbiota features and host traits reveal the complex networks linking all components of the fecal microbiota together and to their host, thus strengthening the holobiont concept. By estimating the heritabilities of microbiota-associated phenotypes, our study quantifies the impact of the host genetics on the fecal microbiota composition. In addition, genetic correlations between taxonomic groups and between taxa abundances and host performance suggest potential applications for selective breeding to improve host traits or promote a healthier microbiota.

Genetic correlation estimates between calving ease in primiparous cows and economically important traits in Nellore cattle.

This study aimed to estimate (co)variance components and genetic parameters for calving ease (CE) and their genetic correlations with growth, reproductive, carcass, and feed efficiency traits in Nellore cattle. Phenotypes for CE are scored in two categories: normal calving and assisted calving. The traits considered were probability of precocious calving, age at first calving, stayability, adjusted scrotal circumference at 365 days of age, accumulated cow productivity, age at puberty of males, gestation length, birth weight, adjusted weights at 210 and 450 days of age, adult cow weight, frame score, hip height, rib eye area, subcutaneous backfat thickness, rump fat thickness, intramuscular fat percentage, residual feed intake and dry matter intake. The estimation of genetic parameters was performed using a two-trait threshold-linear animal model, except for CE, stayability, and probability of precocious calving, which were evaluated through a two-trait threshold animal model. The direct (0.27) and maternal (0.19) heritability estimates for CE in heifers primiparous Nellore indicated that selecting for this trait is feasible. The selection to improve the female sexual precocity should consider CE during the selection and mating decisions to reduce calving problems. Genetic correlation estimates between CE and BW suggest that selecting low birth weight to reduce calving problems is not an appropriate strategy to improve calving ease in heifers Nellore. Therefore, adopting a multi-trait selection model with CE and BW in the Nellore breed would reduce calving difficulties, particularly in sexually precocious heifers, without impairing the growth, reproductive, feed efficiency conversion, and carcass indicator traits.

Distinguishing happiness and meaning in life from depressive symptoms: A GWAS-by-subtraction study in the UK Biobank.

Hedonic (happiness) and eudaimonic (meaning in life) well-being are negatively related to depressive symptoms. Genetic variants play a role in this association, reflected in substantial genetic correlations. We investigated the overlap and differences between well-being and depressive symptoms, using results of Genome-Wide Association studies (GWAS) in UK Biobank. Subtracting GWAS summary statistics of depressive symptoms from those of happiness and meaning in life, we obtained GWASs of respectively "pure" happiness (neffective = 216,497) and "pure" meaning (neffective = 102,300). For both, we identified one genome-wide significant SNP (rs1078141 and rs79520962, respectively). After subtraction, SNP heritability reduced from 6.3% to 3.3% for pure happiness and from 6.2% to 4.2% for pure meaning. The genetic correlation between the well-being measures reduced from 0.78 to 0.65. Pure happiness and pure meaning became genetically unrelated to traits strongly associated with depressive symptoms, including loneliness, and psychiatric disorders. For other traits, including ADHD, educational attainment, and smoking, the genetic correlations of well-being versus pure well-being changed substantially. GWAS-by-subtraction allowed us to investigate the genetic variance of well-being unrelated to depressive symptoms. Genetic correlations with different traits led to new insights about this unique part of well-being. Our results can be used as a starting point to test causal relationships with other variables, and design future well-being interventions.

Unique drought resistance strategies occur among monkeyflower populations spanning an aridity gradient.

PREMISE: Annual plants often exhibit drought-escape and avoidance strategies to cope with limited water availability. Determining the extent of variation and factors underlying the evolution of divergent strategies is necessary for determining population responses to more frequent and severe droughts. METHODS: We leveraged five Mimulus guttatus populations collected across an aridity gradient within manipulative drought and quantitative genetics experiments to examine constitutive and terminal-drought induced responses in drought resistance traits. RESULTS: Populations varied considerably in drought-escape- and drought-avoidance-associated traits. The most mesic population demonstrated a unique resource conservative strategy. Xeric populations exhibited extreme plasticity when exposed to terminal drought that included flowering earlier at shorter heights, increasing water-use efficiency, and shifting C:N ratios. However, plasticity responses also differed between populations, with two populations slowing growth rates and flowering at earlier nodes and another population increasing growth rate. While nearly all traits were heritable, phenotypic correlations differed substantially between treatments and often, populations. CONCLUSIONS: Our results suggest drought resistance strategies of populations may be finely adapted to local patterns of water availability. Substantial plastic responses suggest that xeric populations can already acclimate to drought through plasticity, but populations not frequently exposed to drought may be more vulnerable.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle.

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Heritability and genetic correlations of total and differential somatic cell count with milk yield and composition traits in Italian Simmental cows.

Costs of production have deeply increased each year in the last decades, breeders are continuously looking for more cost effective and more efficient ways to produce milk. Despite the major signs of progress in productivity, it is fundamental to optimize rather than maximize the performances of the dairy cows. Mastitis is still a highly prevalent disease in the dairy sector which causes several economic losses and environmental effect. Its accurate and early diagnosis is crucial to improve profitability of dairy cows and contribute to a more sustainable dairy industry. Among mastitis reduction strategies, there is the urgent need to implement breeding objectives to select cows displaying mastitis resistance by investigating the genetic mechanisms at the base of the inflammatory response. Therefore, in this study we aimed to further understand the genetic background of the differential somatic cell count (DSCC), which provides thorough insights on the actual inflammatory status of the mammary glands. The objectives of this study were to estimate on a cohort of 20,215 Italian Simmental cows over a 3-yr period: (1) the heritability and repeatability values of somatic cell score (SCS) and DSCC, (2) the genetic and phenotypic correlations between these 2 traits and milk production and milk composition traits, (3) the heritability and repeatability values of SCS and DSCC within class of udder health status. Heritability was low both for SCS (0.06) and DSCC (0.08), whereas the repeatability values for these traits were 0.43 and 0.36, suggesting that the magnitude of cow permanent environmental effect for these traits is remarkable. The genetic and phenotypic correlation of SCS with DSCC was 0.612 and 0.605, respectively. Because both significantly differed from the unit, we must consider those traits as different ones. This latter aspect corroborates the need to consider the DSCC as a further indicator of inflammatory status which might be implemented in the Simmental breed genetic evaluation. It is worthy to mention that heritability estimates for SCS and DSCC were the highest in healthy cows compared with the other udder health classes. This implies that when the udder health status changes, it is most likely due to environmental factors rather than aspects related to the animal's genetics. In contrast, the highest additive genetic variance and heritability found for SCS and DSCC in the healthy group might reveal the potential to further implement breeding strategies to select for healthier animals.

Modularity and selection of nectar traits in the evolution of the selfing syndrome in Ipomoea lacunosa (Convolvulaceae).

Although the evolution of the selfing syndrome often involves reductions in floral size, pollen and nectar, few studies of selfing syndrome divergence have examined nectar. We investigate whether nectar traits have evolved independently of other floral size traits in the selfing syndrome, whether nectar traits diverged due to drift or selection, and the extent to which quantitative trait locus (QTL) analyses predict genetic correlations. We use F5 recombinant inbred lines (RILs) generated from a cross between Ipomoea cordatotriloba and Ipomoea lacunosa. We calculate genetic correlations to identify evolutionary modules, test whether trait divergence was due to selection, identify QTLs and perform correlation analyses to evaluate how well QTL properties reflect genetic correlations. Nectar and floral size traits form separate evolutionary modules. Selection has acted to reduce nectar traits in the selfing I. lacunosa. Genetic correlations predicted from QTL properties are consistent with observed genetic correlations. Changes in floral traits associated with the selfing syndrome reflect independent evolution of at least two evolutionary modules: nectar and floral size traits. We also demonstrate directional selection on nectar traits, which is likely to be independent of selection on floral size traits. Our study also supports the expected mechanistic link between QTL properties and genetic correlations.

How early does the selfing syndrome arise? Associations between selfing ability and flower size within populations of the mixed-mater Collinsia verna.

PREMISE: Widespread associations between selfing rate and floral size within and among taxa suggest that these traits may evolve in concert. Does this association develop immediately because of shared genetic and/or developmental control, or stepwise with selection shaping the evolution of one trait following the other? If the former, then association ought to appear within and across selfing populations. We explore this fundamental question in three populations of the mixed-mater Collinsia verna where autonomous selfing (AS) ability has been shown to be under selection by the pollination environment. METHODS: We grew clonal replicates of C. verna in a controlled environment to characterize broad-sense genetic correlations among traits within populations and to assess whether divergence in mating system and floral traits among these populations is consistent with their previously observed selection pressures. RESULTS: As predicted by their respective pollination environments, we demonstrate significant genetic divergence among populations in AS ability. However, patterns of divergence in floral traits (petal, stamen, and style size, stigmatic receptivity, and stigma-anther distance) were not as expected. Within populations, genetic variation in AS appeared largely independent from floral traits, except for a single weak negative association in one population between flower size and AS rate. CONCLUSIONS: Together, these results suggest that associations between selfing rate and floral traits across Collinsia species are not reflected at microevolutionary scales. If C. verna were to continue evolving toward the selfing syndrome, floral trait evolution would likely follow stepwise from mating system evolution.

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

PREMISE: The evolutionary response of a trait to environmental change depends upon the level of additive genetic variance. It has been long argued that sustained selection will tend to deplete additive genetic variance as favored alleles approach fixation. Non-additive genetic variance, due to interactions among alleles within and between loci, does not immediately contribute to an evolutionary response. However, shifts in the allele frequencies within and between interacting loci may convert non-additive variance into additive variance. Here we consider the possibility that an environmental shift may alter allelic interactions in ways that convert dominance into additive genetic variance. METHODS: We grew a pedigreed population of Brassica rapa in greenhouse and field conditions. The field conditions mimicked agricultural conditions from which the base population was drawn, while the greenhouse featured benign conditions. We used Bayesian models to estimate the additive, dominance, and maternal components of quantitative genetic variance. We also estimated genetic correlations across environments using parental breeding values. RESULTS: Although the additive genetic variance was elevated in the greenhouse condition, no consistent pattens emerged that would indicate a conversion of dominance variance. The unusually low genetic variance and broad confidence intervals for the variance estimates obtained through this analysis preclude definitive interpretations. CONCLUSIONS: Further studies are needed to determine whether between-environment changes in additive genetic variance can be traced to conversion of dominance variance.

Characterisation of the genetic relationship between the domains of sleep and circadian-related behaviours with substance use phenotypes.

Sleep problems and substance use frequently co-occur. While substance use can result in specific sleep deficits, genetic pleiotropy could explain part of the relationship between sleep and substance use and use disorders. Here we use the largest publicly available genome-wide summary statistics of substance use behaviours (N = 79,729-632,802) and sleep/activity phenotypes to date (N = 85,502-449,734) to (1) assess the genetic overlap between substance use behaviours and both sleep and circadian-related activity measures, (2) estimate clusters from genetic correlations and (3) test processes of causality versus genetic pleiotropy. We found 31 genetic correlations between substance use and sleep/activity after Bonferroni correction. These patterns of overlap were represented by two genetic clusters: (1) tobacco use severity (age of first regular tobacco use and smoking cessation) and sleep health (sleep duration, sleep efficiency and chronotype) and (2) substance consumption/problematic use (drinks per day and cigarettes per day, cannabis use disorder, opioid use disorder and problematic alcohol use) and sleep problems (insomnia, self-reported short sleep duration, increased number of sleep episodes, increased sleep duration variability and diurnal inactivity) and measures of circadian-related activity (L5, M10 and sleep midpoint). Latent causal variable analyses determined that horizontal pleiotropy (rather than genetic causality) underlies a majority of the associations between substance use and sleep/circadian related measures, except one plausible genetically causal relationship for opioid use disorder on self-reported long sleep duration. Results show that shared genetics are likely a mechanism that is at least partly responsible for the overlap between sleep and substance use traits.

Re-assessing the importance of linear type traits in predicting genetic merit for survival in an aging Holstein-Friesian dairy cow population.

The cumulative improvement achieved in the genetic merit for reproductive performance in dairy populations will likely improve dairy cow longevity; therefore, it is time to reassess whether linear type traits are still suitable predictors of survival in an aging dairy cow population. The objective of the present study was therefore to estimate the genetic correlations between linear type traits and survival from one parity to the next and, in doing so, evaluate if those genetic correlations change with advancing parity. After edits, 152,894 lactation survival records (first to ninth parity) were available from 52,447 Holstein-Friesian cows, along with linear type trait records from 52,121 Holstein-Friesian cows. A series of bivariate random regression models were used to estimate the genetic covariances between survival in different parities and each linear type trait. Heritability estimates for survival per parity ranged from 0.02 (SE = 0.004; first parity) to 0.05 (SE = 0.01; ninth parity). Pairwise genetic correlations between survival among different parities varied from 0.42 (first and ninth parity) to 1.00 (eighth to ninth parity), with the strength of these genetic correlations being inversely related to the interval between the compared parities. The genetic correlations between survival and the individual linear type traits varied across parities for 9 of the 20 linear type traits examined, but the correlations with only 3 of these linear type traits strengthened as the cows aged; these 3 traits were rear udder height, teat length, and udder depth. Given that linear type traits are frequently scored in first parity and are genetically correlated with survival in older parities, they may be suitable early predictors of survival, especially for later parity cows. Additionally, the direction of the genetic correlations between survival and rear udder height, teat length, and udder depth did not change between parities; hence, selection for survival in older parities using these linear type traits should not hinder genetic improvement for survival in younger parities.

Maternal reactivity of ewes at lambing is genetically linked to their behavioural reactivity in an arena test.

In sheep, the bond between the dam and her lambs is established during the first hours of a lamb's life. Genetic variability for behavioural reactivity of ewes assessed in an arena test performed 24 hr after lambing has already been reported. However, there is no evidence that this reactivity represents the ewe's maternal reactivity at lambing in outdoor conditions. The objective of this study was to investigate whether or not the behavioural reactivity of ewes in the arena test is genetically related to their maternal reactivity measured at lambing. A total of 935 Romane ewes were studied. The maternal reactivity of ewes at the outdoor lambing site was recorded in response to a human approach and to the handling of the lambs. Their behavioural reactivity was also recorded 24 hr post-lambing in the arena test that involved a separation from the litter and a human presence. Flight distance, aggressive reaction, time to restore contact with the litter, maternal behaviour scores and vocalizations recorded at the lambing site were heritable (0.12-0.34). All of these behaviours were genetically correlated with the behavioural reactivity in the arena test. The highest genetic correlations (from 0.60 to 0.90) were found amongst maternal behavioural scores, flight distance and high-pitched bleats. In conclusion, behavioural reactivity in the arena test can be used to assess early maternal reactivity in standardized conditions. Phenotyping of ewes' behavioural reactivity with a simplified arena test can be performed for genetic improvement in maternal behaviour in sheep.

Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout.

BACKGROUND: In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n-3 long-chain polyunsaturated fatty acids (PUFAs)-especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids-in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection. RESULTS: Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait. CONCLUSION: This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Costs of reproduction under experimental climate change across elevations in the perennial forb Boechera stricta.

Investment in current reproduction can reduce future fitness by depleting resources needed for maintenance, particularly under environmental stress. These trade-offs influence life-history evolution. We tested whether climate change alters the future-fitness costs of current reproduction in a large-scale field experiment of Boechera stricta (Brassicaceae). Over 6 years, we simulated climate change along an elevational gradient in the Rocky Mountains through snow removal, which accelerates snowmelt and reduces soil water availability. Costs of reproduction were greatest in arid, lower elevations, where high initial reproductive effort depressed future fitness. At mid-elevations, initial reproduction augmented subsequent fitness in benign conditions, but pronounced costs emerged under snow removal. At high elevation, snow removal dampened costs of reproduction by prolonging the growing season. In most scenarios, failed reproduction in response to resource limitation depressed lifetime fecundity. Indeed, fruit abortion only benefited high-fitness individuals under benign conditions. We propose that climate change could shift life-history trade-offs in an environment-dependent fashion, possibly favouring early reproduction and short lifespans in stressful conditions.

Psychological trauma and the genetic overlap between posttraumatic stress disorder and major depressive disorder.

BACKGROUND: Posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) are commonly reported co-occurring mental health consequences of psychological trauma exposure. The disorders have high genetic overlap. Trauma is a complex phenotype but research suggests that trauma sensitivity has a heritable basis. We investigated whether sensitivity to trauma in those with MDD reflects a similar genetic component in those with PTSD. METHODS: Genetic correlations between PTSD and MDD in individuals reporting trauma and MDD in individuals not reporting trauma were estimated, as well as with recurrent MDD and single-episode MDD, using genome-wide association study (GWAS) summary statistics. Genetic correlations were replicated using PTSD data from the Psychiatric Genomics Consortium and the Million Veteran Program. Polygenic risk scores were generated in UK Biobank participants who met the criteria for lifetime MDD (N = 29 471). We investigated whether genetic loading for PTSD was associated with reporting trauma in these individuals. RESULTS: Genetic loading for PTSD was significantly associated with reporting trauma in individuals with MDD [OR 1.04 (95% CI 1.01-1.07), Empirical-p = 0.02]. PTSD was significantly more genetically correlated with recurrent MDD than with MDD in individuals not reporting trauma (rg differences = ~0.2, p < 0.008). Participants who had experienced recurrent MDD reported significantly higher rates of trauma than participants who had experienced single-episode MDD (χ2 > 166, p < 0.001). CONCLUSIONS: Our findings point towards the existence of genetic variants associated with trauma sensitivity that might be shared between PTSD and MDD, although replication with better powered GWAS is needed. Our findings corroborate previous research highlighting trauma exposure as a key risk factor for recurrent MDD.

Predator- and competitor-induced responses in amphibian populations that evolved different levels of pesticide tolerance.

Exposure to agrochemicals can drive rapid phenotypic and genetic changes in exposed populations. For instance, amphibian populations living far from agriculture (a proxy for agrochemical exposure) exhibit low pesticide tolerance, but they can be induced to possess high tolerance following a sublethal pesticide exposure. In contrast, amphibian populations close to agriculture exhibit high, constitutive tolerance to pesticides. A recent study has demonstrated that induced pesticide tolerance appears to have arisen from plastic responses to predator cues. As a result, we might expect that selection for constitutive pesticide tolerance in populations near agriculture (i.e., genetic assimilation) will lead to the evolution of constitutive responses to natural stressors. Using 15 wood frog (Rana sylvatica) populations from across an agricultural gradient, we conducted an outdoor mesocosm experiment to examine morphological (mass, body length, and tail depth) and behavioral responses (number of tadpoles observed and overall activity) of tadpoles exposed to three stressor environments (no-stressor, competitors, or predator cues). We discovered widespread differences in tadpole traits among populations and stressor environments, but no population-by-environment interaction. Subsequent linear models revealed that population distance to agriculture (DTA) was occasionally correlated with tadpole traits in a given environment and with magnitudes of plasticity, but none of the correlations were significant after Bonferroni adjustment. The magnitudes of predator and competitor plasticity were never correlated with the magnitude of pesticide-induced plasticity that we documented in a companion study. These results suggest that while predator-induced plasticity appears to have laid the foundation for the evolution of pesticide-induced plasticity and its subsequent genetic assimilation, inspection of population-level differences in plastic responses show that the evolution of pesticide-induced plasticity has not had a reciprocal effect on the evolved plastic responses to natural stressors.