The evolution of antimicrobial peptide resistance in Pseudomonas aeruginosa is severely constrained by random peptide mixtures.

The prevalence of antibiotic-resistant pathogens has become a major threat to public health, requiring swift initiatives for discovering new strategies to control bacterial infections. Hence, antibiotic stewardship and rapid diagnostics, but also the development, and prudent use, of novel effective antimicrobial agents are paramount. Ideally, these agents should be less likely to select for resistance in pathogens than currently available conventional antimicrobials. The usage of antimicrobial peptides (AMPs), key components of the innate immune response, and combination therapies, have been proposed as strategies to diminish the emergence of resistance. Herein, we investigated whether newly developed random antimicrobial peptide mixtures (RPMs) can significantly reduce the risk of resistance evolution in vitro to that of single sequence AMPs, using the ESKAPE pathogen Pseudomonas aeruginosa (P. aeruginosa) as a model gram-negative bacterium. Infections of this pathogen are difficult to treat due the inherent resistance to many drug classes, enhanced by the capacity to form biofilms. P. aeruginosa was experimentally evolved in the presence of AMPs or RPMs, subsequentially assessing the extent of resistance evolution and cross-resistance/collateral sensitivity between treatments. Furthermore, the fitness costs of resistance on bacterial growth were studied and whole-genome sequencing used to investigate which mutations could be candidates for causing resistant phenotypes. Lastly, changes in the pharmacodynamics of the evolved bacterial strains were examined. Our findings suggest that using RPMs bears a much lower risk of resistance evolution compared to AMPs and mostly prevents cross-resistance development to other treatments, while maintaining (or even improving) drug sensitivity. This strengthens the case for using random cocktails of AMPs in favour of single AMPs, against which resistance evolved in vitro, providing an alternative to classic antibiotics worth pursuing.

Modeling Water Interactions with Graphene and Graphite via Force Fields Consistent with Experimental Contact Angles.

Accurate simulation models for water interactions with graphene and graphite are important for nanofluidic applications, but existing force fields produce widely varying contact angles. Our extensive review of the experimental literature reveals extreme variation among reported values of graphene-water contact angles and a clustering of graphite-water contact angles into groups of freshly exfoliated (60° ± 13°) and not-freshly exfoliated graphite surfaces. The carbon-oxygen dispersion energy for a classical force field is optimized with respect to this 60° graphite-water contact angle in the infinite-force-cutoff limit, which in turn yields a contact angle for unsupported graphene of 80°, in agreement with the mean of the experimental results. Interaction force fields for finite cutoffs are also derived. A method for calculating contact angles from pressure tensors of planar equilibrium simulations that is ideally suited to graphite and graphene surfaces is introduced. Our methodology is widely applicable to any liquid-surface combination.

Serial cardiac biomarkers, pulmonary artery pressures and traditional parameters of fluid status in relation to prognosis in patients with chronic heart failure: Design and rationale of the BioMEMS study.

AIMS: Heart failure (HF), a global pandemic affecting millions of individuals, calls for adequate predictive guidance for improved therapy. Congestion, a key factor in HF-related hospitalizations, further underscores the need for timely interventions. Proactive monitoring of intracardiac pressures, guided by pulmonary artery (PA) pressure, offers opportunities for efficient early-stage intervention, since haemodynamic congestion precedes clinical symptoms. METHODS: The BioMEMS study, a substudy of the MONITOR-HF trial, proposes a multifaceted approach integrating blood biobank data with traditional and novel HF parameters. Two additional blood samples from 340 active participants in the MONITOR-HF trial were collected at baseline, 3-, 6-, and 12-month visits and stored for the BioMEMS biobank. The main aims are to identify the relationship between temporal biomarker patterns and PA pressures derived from the CardioMEMS-HF system, and to identify the biomarker profile(s) associated with the risk of HF events and cardiovascular death. CONCLUSION: Since the prognostic value of single baseline measurements of biomarkers like N-terminal pro-B-type natriuretic peptide is limited, with the BioMEMS study we advocate a dynamic, serial approach to better capture HF progression. We will substantiate this by relating repeated biomarker measurements to PA pressures. This design rationale presents a comprehensive review on cardiac biomarkers in HF, and aims to contribute valuable insights into personalized HF therapy and patient risk assessment, advancing our ability to address the evolving nature of HF effectively.

Lipoprotein(a), Oxidized Phospholipids, and Progression to Symptomatic Heart Failure: The CASABLANCA Study.

BACKGROUND: Higher lipoprotein(a) and oxidized phospholipid concentrations are associated with increased risk for coronary artery disease and valvular heart disease. The role of lipoprotein(a) or oxidized phospholipid as a risk factor for incident heart failure (HF) or its complications remains uncertain. METHODS AND RESULTS: A total of 1251 individuals referred for coronary angiography in the Catheter Sampled Blood Archive in Cardiovascular Diseases (CASABLANCA) study were stratified on the basis of universal definition of HF stage; those in stage A/B (N=714) were followed up for an average 3.7 years for incident stage C/D HF or the composite of HF/cardiovascular death. During follow-up, 105 (14.7%) study participants in stage A/B progressed to symptomatic HF and 57 (8.0%) had cardiovascular death. In models adjusted for multiple HF risk factors, including severe coronary artery disease and aortic stenosis, individuals with lipoprotein(a) ≥150 nmol/L were at higher risk for progression to symptomatic HF (hazard ratio [HR], 1.90 [95% CI, 1.15-3.13]; P=0.01) or the composite of HF/cardiovascular death (HR, 1.71 [95% CI, 1.10-2.67]; P=0.02). These results remained significant after further adjustment of the model to include prior myocardial infarction (HF: HR, 1.89, P=0.01; HF/cardiovascular death: HR, 1.68, P=0.02). Elevated oxidized phospholipid concentrations were similarly associated with risk, particularly when added to higher lipoprotein(a). In Kaplan-Meier analyses, individuals with stage A/B HF and elevated lipoprotein(a) had shorter time to progression to stage C/D HF or HF/cardiovascular death (both log-rank P<0.001). CONCLUSIONS: Among individuals with stage A or B HF, higher lipoprotein(a) and oxidized phospholipid concentrations are independent risk factors for progression to symptomatic HF or cardiovascular death. REGISTRATION: URL: https://wwwclinicaltrials.gov; Unique identifier: NCT00842868.

The role of memory-dependent friction and solvent viscosity in isomerization kinetics in viscogenic media.

Molecular isomerization kinetics in liquid solvent depends on a complex interplay between the solvent friction acting on the molecule, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the dihedral free energy profile. Due to the absence of accurate techniques to directly evaluate isomerization friction, it has not been possible to explore these relationships in full. By combining extensive molecular dynamics simulations with friction memory-kernel extraction techniques we consider a variety of small, isomerising molecules under a range of different viscogenic conditions and directly evaluate the viscosity dependence of the friction acting on a rotating dihedral. We reveal that the influence of different viscogenic media on isomerization kinetics can be dramatically different, even when measured at the same viscosity. This is due to the dynamic solute-solvent coupling, mediated by time-dependent friction memory kernels. We also show that deviations from the linear dependence of isomerization rates on solvent viscosity, which are often simply attributed to internal friction effects, are due to the simultaneous violation of two fundamental relationships: the Stokes-Einstein relation and the overdamped Kramers prediction for the barrier-crossing rate, both of which require explicit knowledge of friction.

Psychedelic Experiences Increase Mind Perception but do not Change Atheist-Believer Status: A Prospective Longitudinal Study.

Recent studies suggest psychedelic use may be associated with changes in a variety of beliefs or belief-like states, including increased 1) mind perception, 2) non-naturalistic beliefs, and 3) Atheist-Believer status (e.g. believer, agnostic, or nonbeliever). We conducted a prospective longitudinal study among participants (N = 657) who planned to have a psilocybin experience outside a laboratory setting. We asked participants about their beliefs concerning mind perception of various entities, specific metaphysical positions, and Atheist-Believer status both before (and after their experience. Replicating previous findings, we observed increases in mind perception across a variety of living and non-living targets (e.g. plants, rocks). However, we found little to no change in metaphysical beliefs (e.g. dualism) or Atheist-Believer status. Taken together, these findings contrast with those from cross-sectional studies that psilocybin experiences result in changes to Atheist-Believer status and non-naturalistic beliefs but support the relevance of mind perception and mentalization.

Investigation of Strategies to Block Downstream Effectors of AT1R-Mediated Signalling to Prevent Aneurysm Formation in Marfan Syndrome.

Cardiovascular outcome in Marfan syndrome (MFS) patients most prominently depends on aortic aneurysm progression with subsequent aortic dissection. Angiotensin II receptor blockers (ARBs) prevent aneurysm formation in MFS mouse models. In patients, ARBs only slow down aortic dilation. Downstream signalling from the angiotensin II type 1 receptor (AT1R) is mediated by G proteins and ß-arrestin recruitment. AT1R also interacts with the monocyte chemoattractant protein-1 (MCP-1) receptor, resulting in inflammation. In this study, we explore the targeting of ß-arrestin signalling in MFS mice by administering TRV027. Furthermore, because high doses of the ARB losartan, which has been proven beneficial in MFS, cannot be achieved in humans, we investigate a potential additive effect by combining lower concentrations of losartan (25 mg/kg/day and 5 mg/kg/day) with barbadin, a ß-arrestin blocker, and DMX20, a C-C chemokine receptor type 2 (CCR2) blocker. A high dose of losartan (50 mg/kg/day) slowed down aneurysm progression compared to untreated MFS mice (1.73 ± 0.12 vs. 1.96 ± 0.08 mm, p = 0.0033). TRV027, the combination of barbadin with losartan (25 mg/kg/day), and DMX-200 (90 mg/kg/day) with a low dose of losartan (5 mg/kg/day) did not show a significant beneficial effect. Our results confirm that while losartan effectively halts aneurysm formation in Fbn1C1041G/+ MFS mice, neither TRV027 alone nor any of the other compounds combined with lower doses of losartan demonstrate a notable impact on aneurysm advancement. It appears that complete blockade of AT1R function, achieved by administrating a high dosage of losartan, may be necessary for inhibiting aneurysm progression in MFS.

The "elongate chelicera problem": A virtual approach in an extinct pterygotid sea scorpion from a 3D kinematic point of view.

Chelicerae, distinctive feeding appendages in chelicerates, such as spiders, scorpions, or horseshoe crabs, can be classified based on their orientation relative to the body axis simplified as either orthognathous (parallel) or labidognathous (inclined), exhibiting considerable diversity across various taxa. Among extinct chelicerates, sea scorpions belonging to the Pterygotidae represent the only chelicerates possessing markedly elongated chelicerae relative to body length. Despite various hypotheses regarding the potential ecological functions and feeding movements of these structures, no comprehensive 3D kinematic investigation has been conducted yet to test these ideas. In this study, we generated a comprehensive 3D model of the pterygotid Acutiramus, making the elongated right chelicera movable by equipping it with virtual joint axes for conducting Range of Motion analyses. Due to the absence in the fossil record of a clear indication of the chelicerae orientation and their potential lateral or ventral movements (vertical or horizontal insertion of joint axis 1), we explored the Range of Motion analyses under four distinct kinematic settings with two orientation modes (euthygnathous, klinogathous) analogous to the terminology of the terrestrial relatives. The most plausible kinematic setting involved euthygnathous chelicerae being folded ventrally over a horizontal joint axis. This configuration positioned the chelicera closest to the oral opening. Concerning the maximum excursion angle, our analysis revealed that the chela could open up to 70°, while it could be retracted against the basal element to a maximum of 145°. The maximum excursion in the proximal joint varied between 55° and 120° based on the insertion and orientation. Our findings underscore the utility of applying 3D kinematics to fossilized arthropods for addressing inquiries on functional ecology such as prey capture and handling, enabling insights into their possible behavioral patterns. Pterygotidae likely captured and processed their prey using the chelicerae, subsequently transporting it to the oral opening with the assistance of other prosomal appendages.

Initial experience, safety, and feasibility using remote access or onsite technical support for complex ablation procedures: results of the REMOTE study.

Aims: Electroanatomical mapping (EAM) systems are essential for the treatment of cardiac arrhythmias. The EAM system is usually operated by qualified staff or field technical engineers from the control room. Novel remote support technology allows for remote access of EAM via online services. Remote access increases the flexibility of the electrophysiological lab, reduces travel time, and overcomes hospital access limitations especially during the COVID-19 pandemic. Here, we report on the feasibility and safety of EAM remote access for cardiac ablation procedures. Methods and results: Mapping and ablation were achieved by combining the EnsiteX™ EAM system and the integrated Ensite™ Connect Remote Support software, together with an integrated audiovisual solution system for remote support (Medinbox). Communication between the operator and the remote support was achieved using an incorporated internet-based common communication platform (Zoom™), headphones, and high-resolution cameras. We investigated 50 remote access-assisted consecutive electrophysiological procedures from September 2022 to February 2023 (remote group). The data were compared with matched patients (n = 50) with onsite support from the control room (control group). The median procedure time was 100 min (76, 120; remote) vs. 86 min (60, 110; control), P = 0.090. The procedural success (both groups 100%, P = 0.999) and complication rate (remote: 2%, control: 0%, P = 0.553) were comparable between the groups. Travel burden could be reduced by 11 280 km. Conclusion: Remote access for EAM was feasible and safe in this single-centre study. Procedural data were comparable to procedures with onsite support. In the future, this new solution might have a great impact on facilitating electrophysiological procedures.

Enigmatic host-mite relationships: Unraveling the distribution of quill mites on Birds-of-Paradise.

Mites of the family Syringophilidae (Acariformes: Prostigmata: Cheyletoidea) are permanent and obligatory parasites of birds. This study presents an analysis of mite material collected from 22 avian species belonging to the family Paradisaeidae (Passeriformes), revealing the presence of four mite species belonging to four genera: Syringophiloidus attenboroughi n. sp., Peristerophila regiusi n. comb., Picobia frankei, and Gunabopicobia garylarsoni. In the present work, the genus Neoperisterophila is synonymized with the genus Peristerophila. While the genera Syringophiloidus and Picobia were expectedly found on paradisaeid birds, given their prevalence in passerines, the presence of Peristerophila and Gunabopicobia was intriguing, suggesting potential host-switching events. The specificity of these mites varies, with some showing occurrence on hosts of closely related genera and others infesting phylogenetically distant hosts. Notably, the distribution of specific mite species on the Birds-of-Paradise appears to be influenced by both long coevolutionary histories and incidental contacts between often unrelated or intergeneric hybrid species of paradisaeid birds. Furthermore, our research of 104 specimens from 22 Birds-of-Paradise species shows generally low infestation rates across the studied species, suggesting a nuanced interaction between these mites and their avian hosts. Additionally, our network analysis provides a deeper understanding of these host-parasite interactions, revealing a high level of specialization and complexity in these ecological relationships.

Long-Term Health-Related Quality of Life following Acute Type A Aortic Dissection with a Focus on Male-Female Differences: A Cross Sectional Study.

Objectives: Acute type A aortic dissection (ATAAD) is a life-threatening cardiovascular emergency, of which the long-term impact on health-related quality of life (HRQoL) and male-female-specific insights remain inadequately clarified. Methods: Consecutive adult ATAAD patients who underwent surgery were retrospectively included between 2007 and 2017 in four referral centers in the Netherlands, and baseline data were collected. The 36-Item Short-Form (SF-36) Health Survey was sent to all survivors between 2019 and 2021 and compared to validated SF-36 scores of the Dutch general population stratified by age group and sex. Results: In total, 324/555 surviving patients returned the SF-36 questionnaire (response rate 58%), of which 40.0% were female; the median follow-up was 6.5 years (range: 1.7-13.9, IQR: 4.0-9.4) after surgery for ATAAD. In comparison to the general population, ATAAD patients scored significantly lower on 6/8 SF-36 subdomains and higher on bodily pain. Differences in HRQoL domains compared to the sex-matched data were largely comparable between sexes, apart from bodily pain. In the age-matched subgroups impaired HRQoL was most pronounced in younger patients aged 41-60 (5/8 impaired domains). Female ATAAD patients scored significantly worse on 5/8 SF-36 subdomains and the physical component summary (PCS) scores than male patients. Age at ATAAD, female sex, hypertension, COPD, and prior thoracic aortic aneurysm were associated with worse PCS scores. Conclusions: Long-term HRQoL was impaired in both male and female ATAAD patients when compared to the general population. Further studies on the nature of this impairment and on interventions to improve HRQoL after ATAAD are clearly warranted, with special attention to females and younger patients.

Accurate Memory Kernel Extraction from Discretized Time-Series Data.

Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate, and insightful, particularly when working with high-resolution time series data. However, in experimental systems, high-resolution data are often unavailable, raising questions about the impact of the data resolution on the estimated GLE parameters. This study demonstrates that direct memory extraction from time series data remains accurate when the discretization time is below the memory time. To obtain memory functions reliably, even when the discretization time exceeds the memory time, we introduce a Gaussian Process Optimization (GPO) scheme. This scheme minimizes the deviation of discretized two-point correlation functions between time series data and GLE simulations and is able to estimate accurate memory kernels as long as the discretization time stays below the longest time scale in the data, typically the barrier crossing time.

Unveiling the ventral morphology of a rare early Cambrian great appendage arthropod from the Chengjiang biota of China.

BACKGROUND: The early Cambrian arthropod clade Megacheira, also referred to as great appendage arthropods, comprised a group of diminutive and elongated predators during the early Palaeozoic era, around 518 million years ago. In addition to those identified in the mid-Cambrian Burgess Shale biota, numerous species are documented in the renowned 518-million-year-old Chengjiang biota of South China. Notably, one species, Tanglangia longicaudata, has remained inadequately understood due to limited available material and technological constraints. In this study, we, for the first time, examined eight fossil specimens (six individuals) utilizing state-of-the-art µCT and computer-based 3D rendering techniques to unveil the hitherto hidden ventral and appendicular morphology of this species. RESULTS: We have identified a set of slender endopodites gradually narrowing distally, along with a leaf-shaped exopodite adorned with fringed setae along its margins, and a small putative exite attached to the basipodite. Our techniques have further revealed the presence of four pairs of biramous appendages in the head, aligning with the recently reported six-segmented head in other early euarthropods. Additionally, we have discerned two peduncle elements for the great appendage. These findings underscore that, despite the morphological diversity observed in early euarthropods, there exists similarity in appendicular morphology across various groups. In addition, we critically examine the existing literature on this taxon, disentangling previous mislabelings, mentions, descriptions, and, most importantly, illustrations. CONCLUSIONS: The µCT-based investigation of fossil material of Tanglangia longicaudata, a distinctive early Cambrian euarthropod from the renowned Chengjiang biota, enhances our comprehensive understanding of the evolutionary morphology of the Megacheira. Its overall morphological features, including large cup-shaped eyes, raptorial great appendages, and a remarkably elongated telson, suggest its potential ecological role as a crepuscular predator and adept swimmer in turbid waters.

Nanoscopic Interfacial Hydrogel Viscoelasticity Revealed from Comparison of Macroscopic and Microscopic Rheology.

Deviations between macrorheological and particle-based microrheological measurements are often considered to be a nuisance and neglected. We study aqueous poly(ethylene oxide) (PEO) hydrogels for varying PEO concentrations and chain lengths that contain microscopic tracer particles and show that these deviations reveal the nanoscopic viscoelastic properties of the particle-hydrogel interface. Based on the transient Stokes equation, we first demonstrate that the deviations are not due to finite particle radius, compressibility, or surface-slip effects. Small-angle neutron scattering rules out hydrogel heterogeneities. Instead, we show that a generalized Stokes-Einstein relation, accounting for an interfacial shell around tracers with viscoelastic properties that deviate from bulk, consistently explains our macrorheological and microrheological measurements. The extracted shell diameter is comparable to the PEO end-to-end distance, indicating the importance of dangling chain ends. Our methodology reveals the nanoscopic interfacial rheology of hydrogels and is applicable to different kinds of viscoelastic fluids and particles.

Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments.

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.

The therapeutic alliance between study participants and intervention facilitators is associated with acute effects and clinical outcomes in a psilocybin-assisted therapy trial for major depressive disorder.

We examined if the therapeutic alliance between study participants and intervention facilitators in a psilocybin-assisted therapy (PAT) trial changed over time and whether there were relationships between alliance, acute psilocybin experiences, and depression outcomes. In a randomized, waiting list-controlled clinical trial for major depressive disorder in adults (N = 24), participants were randomized to an immediate (N = 13) or delayed (N = 11) condition with two oral doses of psilocybin (20mg/70kg and 30mg/70kg). Ratings of therapeutic alliance significantly increased from the final preparation session to one-week post-intervention (p = .03, d = .43). A stronger total alliance at the final preparation session predicted depression scores at 4 weeks (r = -.65, p = .002), 6 months (r = -.47, p = .036), and 12 months (r = -.54, p = .014) post-intervention. A stronger total alliance in the final preparation session was correlated with higher peak ratings of mystical experiences (r = .49, p = .027) and psychological insight (r = .52, p = .040), and peak ratings of mystical experience and psychological insight were correlated with depression scores at 4 weeks (r = -.45, p = .030 for mystical; r = -.75, p < .001 for insight). Stronger total alliance one week after the final psilocybin session predicted depression scores at 4 weeks (r = -.85, p < .001), 3 months (r = -.52, p = .010), 6 months (r = -.77, p < .001), and 12 months (r = -.61, p = .001) post-intervention. These findings highlight the importance of the therapeutic relationship in PAT. Future research should explore therapist and participant characteristics which maximize the therapeutic alliance and evaluate its relationship to treatment outcomes. Trial registration: Registration: Clinicaltrials.gov NCT03181529. https://classic.clinicaltrials.gov/ct2/show/NCT03181529.

Data-driven classification of individual cells by their non-Markovian motion.

We present a method to differentiate organisms solely by their motion based on the generalized Langevin equation (GLE) and use it to distinguish two different swimming modes of strongly confined unicellular microalgae Chlamydomonas reinhardtii. The GLE is a general model for active or passive motion of organisms and particles that can be derived from a time-dependent general many-body Hamiltonian and in particular includes non-Markovian effects (i.e., the trajectory memory of its past). We extract all GLE parameters from individual cell trajectories and perform an unbiased cluster analysis to group them into different classes. For the specific cell population employed in the experiments, the GLE-based assignment into the two different swimming modes works perfectly, as checked by control experiments. The classification and sorting of single cells and organisms is important in different areas; our method, which is based on motion trajectories, offers wide-ranging applications in biology and medicine.

How infection-triggered pathobionts influence virulence evolution.

Host-pathogen interactions can be influenced by the host microbiota, as the microbiota can facilitate or prevent pathogen infections. In addition, members of the microbiota can become virulent. Such pathobionts can cause co-infections when a pathogen infection alters the host immune system and triggers dysbiosis. Here we performed a theoretical investigation of how pathobiont co-infections affect the evolution of pathogen virulence. We explored the possibility that the likelihood of pathobiont co-infection depends on the evolving virulence of the pathogen. We found that, in contrast to the expectation from classical theory, increased virulence is not always selected for. For an increasing likelihood of co-infection with increasing pathogen virulence, we found scenario-specific selection for either increased or decreased virulence. Evolutionary changes, however, in pathogen virulence do not always translate into similar changes in combined virulence of the pathogen and the pathobiont. Only in one of the scenarios where pathobiont co-infection is triggered above a specific virulence level we found a reduction in combined virulence. This was not the case when the probability of pathobiont co-infection linearly increased with pathogen virulence. Taken together, our study draws attention to the possibility that host-microbiota interactions can be both the driver and the target of pathogen evolution. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.