Evolution

Rapid evolution allows populations to persist in environments where they would otherwise go extinct. This phenomenon, known as evolutionary rescue, is typically studied in the framework of biological evolution, yet adaptive traits can also arise and spread through cultural evolution. The present study developed a stochastic eco-evolutionary model to compare rescue probabilities through biological and cultural evolution. Transmission bias governed the rescue probability under cultural evolution by setting how readily a rare adaptive trait was copied. Conformity bias suppressed population persistence because a rare trait was the least likely to be copied. Content bias toward the adaptive trait enabled evolutionary rescue when social learning was rapid, but it typically yielded a lower rescue probability than biological evolution. Only anticonformity bias, together with a high social learning rate, exceeded the rescue probability of biological evolution by enabling the adaptive trait to be established more rapidly. These results demonstrate that transmission bias alters the demographic consequences of cultural evolution and highlight the importance of transmission processes in evolutionary rescue theory. Understanding how adaptive behaviours are socially transmitted may also improve predictions of animal population persistence and inform conservation efforts in rapidly changing environments.

Heritable bacterial symbionts are pervasive in terrestrial arthropods, often imposing reproductive manipulations to promote their own spread within host populations. Co-infections are common, potentially allowing symbiont co-infectors to hitchhike through a host population. However, adverse thermal conditions can disrupt these communities, particularly when co-infectors vary in their thermal sensitivity. We used a multi-generation experiment to test whether warm (29 {degrees}C) conditions disrupted spread of heritable symbionts through uninfected populations of the spider, Mermessus fradeorum. We tested two common infection combinations: a single infection with a cytoplasmic incompatibility (CI) inducing Rickettsiella or a feminizing co-infection that included a feminizing Wolbachia, the same Rickettsiella, and up to three apparent hitchhikers (two additional Wolbachia strains and Tisiphia). We initiated replicate populations with 1/3 of one infection type and 2/3 uninfected spiders, evaluating population infection rate over 5 spider generations under different temperature regimes. Under cool (21{degrees}C) conditions, Wolbachia feminization drove co-infection to 88% and Rickettsiella CI drove single infection to 83% of host populations. Vertical transmission for all symbionts was high (97-99%) and hitchhiking symbionts also spread effectively. Under warm conditions, feminization and CI efficacy were reduced, and symbionts suffered variably reduced vertical transmission. Warm conditions ultimately destroyed the co-infecting symbiont consortium and impeded symbiont spread. On its own, though, Rickettsiella was still able to increase, despite reduced strength of CI. We hypothesize that contrasting tensions between feminizing spread of the symbiont consortium versus environmentally driven loss of function and transmission may explain observed patterns of mixed infections in field populations of this spider.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

The shape of density-dependence governs species persistence, and ecosystem stability. Yet, whether per-capita growth declines sublinearily, or superlinearily with density remains hotly debated. Growth rates across the tree of life have been shown to decline sublinearly with density, whereas theory founded on resource competition predicts the opposite. Here, we resolve this discrepancy and show that sublinearity can readily emerge from geometric constraints on consumer interactions. By linking inter individual spacing, movement and interference rates, we derive two limiting-interference regimes, one of which the well-mixed limit recovers the form of classic Beddington DeAngelis interference response. We then developed an individual-based model from first principles which reproduces the derived sublinearity response, and further use empirical data from published consumer-resource experiments that also bears the signature of sublinear density-dependence. Further, embedding the interference mechanisms underlying the emergence of sublinear density-dependence in coexistence theory opens a new regime for species coexistence where classical theory fails to predict. Our framework indicates that non-consumptive interactions are not merely a correction to resource competition but might be a distinct axis along which diverse communities may potentially coexist.

Standard methods for ancestral sequence reconstruction (ASR) rely on substitution models for the residues in a biological sequence and assume independent evolution across these sites, ignoring the epistatic interactions that shape molecular evolution. In contrast, deep learning models like variational autoencoders (VAEs) can learn low-dimensional representations ("embeddings") of sequences in a protein family that may implicitly handle these dependencies, raising the possibility of performing more accurate ASR by interpolating between extant sequence embeddings within the VAE's latent space. In this study, we test this hypothesis by developing and evaluating a VAE-based ASR pipeline. Benchmarking this approach against established likelihood-based and parsimony methods using various simulations of protein evolution, including scenarios with and without epistasis, we find that the VAE-based approach is consistently and significantly outperformed by standard methods, even in epistatic regimes where it was hypothesized to have an advantage. We further show that this failure is not due to a lack of phylogenetic structure in the latent space, which does contain evolutionary signal. Rather, the primary limitation is the information loss inherent to the autoencoding process: the VAE's decoder cannot generate sequences with sufficient fidelity for the precise demands of ASR.

Infections with bacteria which have evolved multi-drug resistance (MDR) cause millions of deaths worldwide. Large-scale efforts are gathering genotypic and phenotypic data on MDR bacteria, but methods for learning the structure, diversity, and predictors of evolutionary pathways to MDR have yet to take full advantage of these data. Here, we use evolutionary accumulation modelling (EvAM), an emerging class of machine learning methods with roots in cancer progression, to infer these evolutionary pathways across ESKAPEE pathogens (seven bacterial species that dominate health burdens), using a database of over 635k genotyped phenotypic observations from around the world. We identify global patterns in MDR evolutionary pathways, remarkably shared across multiple ESKAPEE species. Species-specific deviations from these stereotypical pathways are connected with geographical and demographic covariates, facilitating predictions of future MDR evolution. We verify these predictions with several hundred new phenotypes from ESKAPEE samples spanning decades of clinical infections in sub-Saharan Africa, demonstrating the capacity to forecast future MDR evolution from these inferred shared pathways.

Across different modalities, signals play a core role in attracting mates and influencing mating success. In several non-human primate species, females produce calls during mating that are thought to promote male competition over receptive females. The extent to which social system characteristics modulate the function of copulation calls remains less clear. We studied copulation calls in wild Guinea baboons (Papio papio), who live in a multilevel society structured around units in which females associate and mate almost exclusively with a single male. We hypothesised that females use copulation calls as an indirect form of mate competition, with competition increasing in larger units. In addition, we hypothesised that females are more likely to mate again after calling. We analysed 6116 copulations between 2014 and 2025, involving 99 reproductively active females and 78 subadult and adult males. Females produced copulation calls in 72.7% of copulations, with large inter-individual variation. Neither unit size nor its interaction with the female's swelling size or the presence of simultaneously receptive females affected the probability of calling. A survival analysis with a subset of the data (2353 copulations) revealed no effect of calling on the latency to the next mating. Our results render the hypothesis that female Guinea baboons use calls in indirect mate competition unlikely. Yet, the probability of calling varied with sexual swelling size, suggesting that calls signal female fertility. Possibly, Guinea baboon copulation calls represent an evolutionary remnant, no longer under selective pressure, and can be considered index signals of female fertility.

In vertebrates, top-down visual attention is a cognitive process where internal goals modulate the tuning of peripheral sensory systems. This leads to increased perceived contrast to both goal-relevant objects and areas of the visual field that are attended. Such a system would also be beneficial to bees, enabling them to detect and recognise the most profitable flowers in their environment. We tested whether bumblebees possess a top-down attentional system resembling that seen in vertebrates. We trained two groups of bees to collect rewards under high contrast targets. To potentially induce a difference in attention while searching for the targets, one group received a higher concentration of sucrose rewards compared to the other. During tests, the targets were presented with a series of lower contrasts to measure the contrast sensitivity curves of the bees induced by the different learnt reward levels. We predicted a stronger effect of any attention-like process on contrast sensitivity in the high reward group. We also repeated this experiment with the neonicotinoid pesticide imidacloprid dissolved in the sucrose rewards to test whether this affects bee attention. Across all test contrasts, higher rewards significantly increased bee accuracy when locating targets, lowered contrast thresholds and reduced the latency to make first choices. Imidacloprid reduced bee accuracy but did not influence first choice latency. These results suggest that learnt floral rewards can influence bee behavioural contrast sensitivity in a manner resembling vertebrate top-down attention and that imidacloprid may modulate this through effects on their nervous system.

Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.

Developmental regression is a severe but poorly understood complication of several neurodevelopmental disorders. In Kleefstra syndrome (KLEFS1), caused by EHMT1 haploinsufficiency, regression often emerges during adolescence or early adulthood and is frequently preceded by marked sleep disturbance. Experimental work implicating EHMT1/G9a in metabolic regulation and stress responses raises the possibility that impaired metabolic resilience contributes to this vulnerability. Here, we aimed to investigate whether altered glucose metabolism is a feature of KLEFS1 and whether it relates to clinical variability, including regression. Through [18F]FDG-PET/CT, individuals with KLEFS1 who had experienced regression (n=4) exhibited a hypometabolic brain profile, whereas one individual who had not experienced regression showed globally elevated metabolic activity. In parallel, G9a mutant flies exhibited increased baseline metabolic rate and neuronal ATP levels together with sleep fragmentation resembling the clinical phenotype. Providing flies with oxidative stress to model KLEFS1 regression further exacerbated sleep disruption and was associated with a reduction in metabolic output. Importantly, adult high sugar feeding in flies prevented oxidative stress-induced worsening of sleep and maintained metabolic stability under challenge. Together, these findings suggest that regression in KLEFS1 and associated sleep disturbances are linked to underlying metabolic vulnerability and impaired maintenance of energy homeostasis under stress.

Background/Objectives: Childhood appetitive traits are heritable behavioural phenotypes hypothesized to link genetic susceptibility to obesity risk. Yet their genetic architecture and role in mediating polygenic adiposity risk remain poorly understood. Methods: We conducted the largest survey of childhood eating behaviour to date, allowing us to perform genome-wide association studies of six appetitive domains derived from 18 items of the parent-reported Children's Eating Behaviour Questionnaire in up to 31,018 eight-year-old children from the Norwegian Mother, Father and Child Cohort Study (MoBa). A trio-based design enabled decomposition of direct and indirect genetic effects on appetite and BMI. Results: We identified ten independent genome-wide significant loci for childhood eating behaviour, primarily across Food Responsiveness, Satiety Responsiveness, and Food Fussiness, eight of which lie at established childhood or adult BMI loci. Food Responsiveness and Satiety Responsiveness showed both phenotypic and genetic correlations with BMI trajectories from early childhood through adolescence. Statistical mediation analyses indicated that 22.1% and 10.4% of the aggregated genetic association with BMI at age 8 could be decomposed through these traits, respectively. Locus-specific patterns further suggested mechanistic pathways, with the FTO locus acting predominantly via Food Responsiveness, and the ADCY3 locus via Satiety Responsiveness. Trio analyses demonstrated that both BMI and eating behaviour associations were predominantly explained by children's inherited alleles, with minimal contribution from indirect effect from parental adiposity, although parental genetic liability influenced reporting of Satiety Responsiveness. Conclusions: Childhood appetitive traits capture a substantial proportion of genetic susceptibility to adiposity through distinct eating behaviour pathways (under standard mediation assumptions). These effects are primarily driven by the child's own genotype rather than indirect parental influences, positioning appetite as a plausible, biologically grounded target for early obesity prevention.

Ecologists remotely track movement steps of animals (e.g., via global positioning systems) and use step selection functions to study the effect of environmental factors upon their movement decisions. Constructing such functions requires pairing each observed step with some number of unobserved but feasible comparison steps. Larger numbers of comparison steps generally yield better estimates but also incur potentially challenging computational demands. Thus, it is important to determine an appropriate number of comparison steps. No established guidance exists for this decision. Here, we use simulated tracks to assess how many comparison steps are needed, fitting each set of steps to a conditional logistic regression model. We monitor errors in estimated effects for several classes of tracks, identifying the number of comparison steps for which mean relative absolute error in estimated effects is consistently low. By this criterion, 32 comparison steps per observed step are needed for our primary class of simulated tracks. Tracks in more homogeneous landscapes, tracks with shorter mean step lengths, or shorter tracks generally require more comparison steps (ranging from 64 to 128 per observed step) to achieve the same level of accuracy. Longer tracks generally require fewer comparison steps (16 per observed step). These results clearly demonstrate that the number of comparison steps influences how well step selection functions estimate covariate effects and provides initial direction in a research area that currently lacks quantitative guidance. Movement ecologists should take care when selecting the number of comparison steps paired with each observed step because those decisions matter.

Three new species of sun orchid (Thelymitra) endemic to Aotearoa New Zealand are here described. These are T. palustris, T. scabrifolia and T. semaphora. The morphological distinctiveness of these three species has been acknowledged for decades; however, their taxonomic status has remained unresolved. Evidence from existing karyological data, recently generated DNA sequence data (LFY and ycf1) and morphological studies from historical and fresh collections are used here to support their formal description. Both, T. palustris and T. semaphora are restricted to wet habitats north of Auckland (North Island). Thelymitra scabrifolia inhabits mostly scrub, and it has a similar northern North Island distribution, but is has been found also in Manawat[a]whi / Three Kings Islands and historically in Otago (South Island). All three species are polyploids and are of conservation concern.

Plant cold acclimation emerges from coordinated adjustments in photosynthesis, primary metabolism, and intracellular carbon allocation. Yet, the regulatory role of subcellular metabolite compartmentation in natural variation of cold acclimation remains insufficiently understood. Here, we investigated four Arabidopsis thaliana accessions grown either individually or in bulk to determine how growth configuration and genotype shape the metabolism of sugars and organic acids during cold exposure. Using non-aqueous fractionation, we quantified plastidial, cytosolic, and vacuolar sugar pools alongside whole-cell carbohydrates, organic acids, enzyme activities, photosynthetic parameters, and stress markers. A neural-network classifier revealed that subcellular sugar distribution together with sugar amounts and organic acids provided the strongest discriminatory power among accessions, surpassing photosynthetic traits and enzyme activities. Our findings demonstrate that natural variation in cold acclimation is strongly determined by genotype-specific subcellular metabolite architectures, and that the cultivation strategy modulates these intracellular signatures. We conclude that subcellular compartmentation of metabolites represents a cellular control layer for natural variation of cold acclimation and resilience in Arabidopsis thaliana.

Unmeasured confounding threatens causal inference and replicability in observational multi-omic studies across variable environments. Genetic instrumental variables (Mendelian randomization) and negative-control calibration each address complementary sources of unmeasured confounding, yet no existing framework unifies them for omics-scale mediation analysis. We introduce ICONIC, an R package that embeds genetic instruments and negative controls within a proximal causal inference framework for total-effect and mediation analysis. ICONIC implements eight estimators spanning five confounding-control strategies, supports continuous, binary, and time-to-event outcomes, and provides extensive diagnostics including sensitivity analyses that map estimator performance across plausible assumptions. Ground-truth benchmarks are calibrated to real-omics covariance structures via a hybrid generative model (GAN + feature-level Gaussian copula) rather than parametric simulation, and a companion planning tool predicts performance gains from collecting additional omic data. We demonstrate ICONIC in two case studies: identifying placental transcriptomic mediators of gestational diabetes on birth weight (n = 164), and tumor-expression mediators of smoking intensity on lung cancer survival (n = 494). Notably, ICONIC's diagnostics recommended different estimation strategies across the two scenarios, reflecting differences in the likely influence of unmeasured confounding. ICONIC is freely available at https://github.com/sbresnahan/iconic/.

Circadian clocks enable mosquitoes to anticipate recurring environmental variations and coordinate behaviors critical for survival and disease transmission, such as locomotion, reproduction, host-seeking, and blood-feeding, with times of day when performance is maximal. In Aedes aegypti, locomotor activity follows a robust diurnal rhythm shaped by endogenous circadian clocks and environmental cues, among which light has been shown to be the primary source of temporal information. While early studies established the role of light in regulating locomotor activity, behavior, oviposition and pupation, it remains unclear which features of a light cycle drive changes in circadian rhythms. This question is increasingly relevant as Ae. aegypti is frequently exposed to artificial and dynamic lighting conditions in urban environments. Here, we investigated how transient changes in light schedules influence circadian rhythms in locomotor activity by systematically manipulating the timing, duration, and direction of light exposure. Using a high-throughput assay, we tested over 1900 individuals, including wild-type and timeless knockout mutants, and showed that a single day of al tered lighting is sufficient to induce robust phase shifts, with no evidence of masking effects. A 6-hour light pulse was sufficient to re-entrain mosquitoes regardless of the timing of the pulse, and phase shifts were primarily driven by the offset time of the light pulse, indicating that light-offset acts as a major zeitgeber. Together, these findings challenge conventional assumptions about the timescale of circadian synchronization and highlight the remarkable plasticity of mosquito behavior in response to anthropogenic light. Eventually, these effects could explain the rapid adaptation of the species to urban environments and have potential consequences for disease transmission dynamics.

Standard polygenic risk scores (PRSs) are constructed based on additive genome-wide association study (GWAS) summary statistics. Nonlinear machine learning methods have been increasingly applied to construct PRSs directly from individual-level data, with the aim of improving predictive performance over standard PRSs through their ability to model non-additive genetic effects. However, their superiority across studies has been inconsistent, and the conditions under which they provide meaningful improvements remain unclear. We combined theoretical analysis, simulations and a real-world application to investigate when two widely used nonlinear machine learning methods, random forest and XGBoost, outperform standard PRSs. Theoretical analysis showed that standard PRSs can implicitly capture part of the genetic variance attributable to nonadditive genetic effects through their contributions to marginal SNP effects, thereby losing less information than commonly assumed. Although nonlinear models have a higher theoretical potential, their greater flexibility incurs a bias-variance trade-off that can limit predictive gains at finite sample sizes. Simulations showed that XGBoost outperformed the standard PRS only when the genetic architecture involves a sufficiently large proportion of interaction genetic variance concentrated across relatively few interaction effects and large training samples were available. Random forest consistently underperformed the standard PRS. In an application to ischemic heart disease prediction using UK Biobank data, XGBoost showed no meaningful improvement in predictive performance over the standard PRS, whereas random forest again performed worse. Together, these findings suggest that nonlinear machine learning do not uniformly outperform standard PRSs; rather, their relative performance depends jointly on genetic architecture and training sample size. Our study helps to reconcile the inconsistent results reported across previous studies and provides a framework for identifying settings in which more complex PRS models are likely to be beneficial.

Brettanomyces bruxellensis is a yeast species associated with diverse fermentation environments and characterized by extensive genetic diversity, including diploid, autotriploid, and allotriploid lineages resulting from independent hybridization events. These lineages are associated with distinct ecological niches and provide a framework for studying metabolic trait evolution in complex genomes. Nitrate assimilation is a relatively uncommon trait among yeasts and has been reported in B. bruxellensis, but its distribution and evolutionary history within the species remain poorly understood. Here, we combined phenotypic characterization of 151 strains with genomic analyses of 946 whole-genome sequences to investigate nitrate assimilation. Growth assays revealed that nitrate assimilation is widespread but unevenly distributed across genetic lineages, with some populations largely retaining the trait whereas others have frequently lost it. Genomic analyses identified extensive variation affecting the nitrate assimilation gene cluster composed of YNR1, YNI1, and YNT1. Nitrate assimilation was strongly associated with both gene copy number and predicted gene functionality, with nitrate-assimilating strains generally carrying more functional copies of the cluster. Leveraging the complex genomic architecture of the species, we independently analyzed primary and acquired genomes in allotriploid lineages and uncovered contrasting evolutionary trajectories following hybridization. While nitrate assimilation genes were generally maintained in primary genomes, acquired genomes showed a higher prevalence of gene loss and predicted loss-of-function variants, revealing asymmetric dynamics between subgenomes. Altogether, our results suggest that nitrate assimilation represents an ancestral trait that has been differentially maintained across B. bruxellensis lineages through a combination of copy number variation, gene degeneration, and genome-specific evolutionary dynamics. These findings provide new insights into how genome architecture and polyploid evolution shape the maintenance and loss of metabolic traits in an industrially relevant yeast species.

Prions, once mainly studied for their pathogenic roles, are now gaining recognition as adaptive elements in microbial physiology. Over one-third of wild yeast isolates harbor prion proteins, yet their impact on host-microbe interactions remains poorly characterized. Given the ecological dominance of yeasts in the Drosophila mycobiome, we leveraged the Drosophila melanogaster-Saccharomyces cerevisiae system to investigate how the mycobiome-derived prion, [MRPL10+], modulates host physiology. We show that flies exposed to [MRPL10+] yeast exhibit significantly enhanced cold tolerance and increased locomotor activity. This effect persists with heat-killed yeast and diluted culture, suggesting a stable, potent bioactive factor. Using the genetically diverse Drosophila Global Diversity Lines (GDL), we identified natural variation in responsiveness to [MRPL10+] yeast. Genome-wide association and functional RNAi screening revealed a gut-brain signaling axis involving genes critical for digestion, intercellular communication, transcription regulation, and neural transmission. Notably, serotonin and octopamine pathways were essential for [MRPL10+]-induced changes in cold tolerance and locomotion, implicating neuromodulatory circuits in prion-mediated microbial signaling. Our findings establish a mechanistic link between a fungal prion and host metabolic and neural adaptation. This work provides the first genetic dissection of a prion-mediated host-microbe interaction, laying the groundwork for investigating beneficial prions in complex microbial communities and highlighting a new dimension of the mycobiomes influence on animal physiology.

Objective: To determine whether the most cited and the most publicly discussed benign prostatic hyperplasia (BPH) literature describe the same body of work, and whether clinicians and patients read different evidence. Methods: Four Boolean Web of Science searches and four matched Altmetric Explorer searches were run in February 2024, restricted to literature indexed with urologic terminology, and screened in duplicate. Two arms were assembled: the 50 most-cited articles (citation census May 2024) and the 50 with the highest Altmetric Attention Scores. Funding source and intervention focus were hand-coded from the full text of all 100 articles; open-access status came from Unpaywall. Arms were compared with Mann-Whitney U and Fisher exact tests Results: Eight of 50 articles (16%) appeared in both arms. Citation selected articles were older (median 2008 versus 2018, p<0.001), more often randomized trials (58% versus 22%, p<0.001), and more often published in urology-specific journals (96% versus 66%, p<0.001). Industry funded 50% versus 18% of articles (p=0.001) and non-industry sources 8% versus 36% (p=0.001). Only 10% of citation-selected articles were open access versus 56% (p<0.001). Attention data were recoverable for only 13 citation-selected articles (median score 9 versus 17). The two arms were cited in the 2026 AUA BPH Guideline at indistinguishable rates (22% versus 20%, p=1.00). Conclusions: These are largely distinct bodies of work, separated most decisively by whether they can be read without a subscription, yet both inform guideline development equally. Patients arrive with evidence systematically different from, and no less guideline-relevant than, that underpinning their urologist's training.