Ecological specialization emerges when adaptation to a focal context increases fitness in that context relative to others. Experimental evolution has been widely used to study microbial specialization in abiotic environments but not predatory specialization. Here we demonstrate evolutionary specialization by a bacterial predator and characterize associated diversification of its predation profile. Populations of the bacterium Myxococcus xanthus evolving on single, non-evolving prey species diverged in predatory performance, showing increased performance on their home prey relative to their common ancestor, and relative to foreign prey not encountered during adaptation. Adaptation to the single-prey environments resulted in striking radiation of performance profiles across a diverse panel of foreign prey that was shaped interactively by selection, chance, and indirect effects, with home-prey identity modulating the degree of stochastic indirect diversification. Despite a great diversity of indirect evolutionary effects, correlated evolution was net-positive, yielding positive predatory specialization as the general outcome. Genomic evolution mirrored phenotypic evolution in that degrees of genomic parallelism differed as a function of home-prey identity. These findings show that adaptation to even simple biotic conditions can generate great ecological and behavioral diversity, linking direct selection, deterministic indirect effects of adaptation, stochasticity, and the origins of predator specialization and diversification.
Science Journals
In modular organisms, where growth and fragmentation blur the boundaries between individuals, the interplay between asexual and sexual reproduction creates complex fitness trade - offs. Life - history theory predicts that resources allocated to one fitness component necessarily reduce investment in others, yet detecting these trade - offs in wild populations of clonal organisms remains challenging. Phenotypic plasticity can enhance survival, yet its influence on reproductive capacity and life - history trade - offs remains poorly understood. Using a fully crossed reciprocal transplant design, we tracked 263 colonies of the branching coral Acropora cervicornis across nine reef sites over 42 months, investigating relationships between fragmentation, morphological plasticity, and the capacity for sexual reproduction. Breakage patterns reflected both environmental and genetic factors. Primary branch breaks created a "double negative" effect - simultaneously more than doubling mortality risk and delaying attainment of a validated reproductive size class by ~40%. Conversely, higher morphological plasticity in surface area - to - volume ratio accelerated sexual maturation up to 6 - fold, counteracting the negative effects of fragmentation. In parallel, a simple demographic model parameterized with published fecundity data estimated that primary breakage reduces expected cumulative reproductive output by ~58%, a result robust across a wide range of parameter assumptions. These results demonstrate a fundamental reproductive trade-off in which asexual reproduction through fragmentation undermines sexual reproductive potential by reducing colony size. Moreover, our findings reveal that fragmentation susceptibility is broadly heritable and subject to selection, and identify a compensatory mechanism through which plasticity enhances fitness beyond immediate survival.
Allostery enables proteins to couple environmental signals to functional outputs, yet how allosteric mechanisms diversify during evolution remains poorly understood. Here, we address this question in the ubiquitous and functionally diverse arsenic repressor (ArsR) superfamily by integrating information-theoretic bioinformatics, structural characterization of DNA recognition and NMR measurements of fast internal dynamics. We identify conserved residues that define the structural scaffold of ArsR proteins and subfamily-specific positions that encode inducer and DNA specificity. In the persulfide sensor SqrR, the crystal structure of the DNA-bound complex reveals how operator specificity is encoded by a limited set of residues, consistent with sequence-derived predictions functionally validated by in vitro transcription assays across divergent ArsR regulators. We further show that allosteric inhibition of DNA binding in SqrR occurs without large-scale conformational rearrangements and is instead associated with changes in internal dynamics, as previously observed for the zinc sensor CzrA. Together, these results support a model in which conformational entropy preserves allosteric connectivity while relaxing sequence constraints, thereby enabling functional diversification within a protein superfamily.
Holcosus orcesi, the Orces Blue Whiptail, is a Critically Endangered lizard endemic to the upper Jubones River basin in southern Ecuador. Restricted to a narrow elevational range within semi-arid Andean shrublands, it represents one of the few montane members of a predominantly lowland lineage. Here we present the first high-quality reference genome for H. orcesi, generated using Oxford Nanopore Technologies long-read sequencing. The assembly spans 1.68 Gb across only 91 contigs, with an N50 of 76.2 Mb and a BUSCO completeness of 96.8%, making it among the most contiguous and complete squamate genomes to date. Structural annotation predicted 25,682 genes, of which 85% showed homology to known proteins and 45% were assigned Gene Ontology terms. Repetitive elements accounted for 46.3% of the genome, with LINEs representing the predominant class. This genome provides a foundational resource for future evolutionary, comparative and conservation-genomic research of H. orcesi and other mountain reptiles, enabling studies of population genomics, local adaptation, and genomic erosion in isolated populations. By expanding the genomic representation of tropical montane reptiles, this work helps address longstanding phylogenetic and geographic gaps in global biodiversity genomics and provides a foundation for evidence-based conservation of H. orcesi and related taxa.
Biological sequences are known to be not random. Thus, the comparison of in silico restriction fragment distributions of random and biological sequences may be an indicator of this non-randomness. Our analyses show that for most of the tested combinations of restriction enzyme and genome sequence the fragments per Megabase of the biological sequence deviate at least more then 10% from the corresponding random sequence. This deviation goes into both directions, i.e. clearly increased values are as common as clearly decreased values. Although there is no species- or restriction-enzyme-specific effect, a clear impact of the GC content both of the restriction site and of the genome sequence can be seen. In contrast to the random sequences, the genome sequences show distinct peaks in their fragment length distributions, hinting to repetitive elements such as transposons.
Cancer progression is increasingly understood as an evolutionary process shaped not only by competition but also by cooperative interactions including those mediated through diffusible ``public goods'' (PGs). Classical evolutionary game theory predicts that PG-producing (altruistic) subclones cannot invade well-mixed populations of non-producers, creating a paradox given their observed emergence in tumors. Here, we resolve this contradiction by combining stochastic spatial simulations with an analytically tractable Moran model to study the invasion dynamics of PG-producing cells in structured populations. Starting from a single producer cell, we explicitly model stochastic PG secretion, diffusion, binding/unbinding, and cell proliferation across biologically relevant parameter ranges. We demonstrate that spatial structure fundamentally alters invasion dynamics, enabling PG producers to invade and establish even when production incurs a fitness cost. Both numerical and analytical approaches converge on a key unifying parameter, a characteristic length scale {delta}, that captures the combined effects of diffusivity, binding kinetics, and degradation. This length scale determines the spatial extent of PG availability and thus the selective advantage of producers. We identify distinct regimes: when PGs are localized (small {delta}), producers preferentially benefit and invasion is likely; when PGs are widely dispersed (large {delta}), benefits are shared and invasion approaches neutrality or is suppressed by costs. Our results highlight that invasion of cooperative traits is governed by spatially mediated resource localization rather than intrinsic fitness alone. This framework provides a mechanistic basis for understanding the emergence of cooperative subclones in tumors and suggests that modulating biophysical transport properties of signaling molecules could influence tumor evolution, metastasis, and therapeutic resistance.
As one of the earliest-diverging multicellular eukaryotic lineages, the bladed Bangiales (Rhodophyta) possess a deep evolutionary history with a central role in the multi-billion-dollar global seaweed aquaculture industry. Although North Atlantic representatives are emerging candidates for regional mariculture, the scarcity of high-quality genomic resources for these taxa hinders both fundamental research and commercial optimization. To address this, we present the first chromosome-level genome assemblies for two native European species: Porphyra dioica (150.44 Mbp) and Porphyra linearis (95.22 Mbp). By integrating Oxford Nanopore Technologies (ONT) long-read sequencing with Hi-C proximity ligation, we generated highly contiguous nuclear genomes resolved into five chromosomes. Structural gene models were predicted through the BRAKER3 pipeline, identifying 12,548 and 10,382 protein-coding genes for P. dioica and P. linearis, respectively. Subsequent homology-based functional annotation characterized 57.4% and 59.8% of these predicted proteins. Supplemented by circularized organellar genomes, these reference genomes provide a critical framework for future research, enabling comparative studies of Atlantic-Pacific divergence and facilitating the development of selective breeding programs for sustainable European aquaculture.
The deployment of clothianidin-based insecticide formulations in malaria vector control has highlighted the capacity of Anopheles funestus to displace more susceptible mosquito species in treated areas and to rapidly evolve resistance under selection pressure. Metabolic detoxification, together with structural and genetic changes in nicotinic acetylcholine receptors (nAChRs), the primary molecular targets of neonicotinoids, can reduce insecticide efficacy. Here, we characterized amino acid substitutions across all 11 nAChR subunits in An. funestus to assess standing variation that may facilitate adaptive responses to chemical exposure. Using whole-genome sequencing data from 656 mosquitoes sampled in 13 African countries, we found marked contrasts in the distribution of nonsynonymous variants among nAChR subunits. Most subunits are strongly constrained and carry no missense variants, whereas two loci (3 and 7) display three geographically widespread amino acid substitutions across the continent. In contrast, 9 and {beta}2 accumulate dozens of nonsynonymous mutations occurring at intermediate to high frequencies, including within domains involved in orthosteric ligand binding and channel gating. Genetic differentiation at nAChR loci among populations from different countries is low to moderate, although several nonsynonymous mutations display high FST values consistent with geographic structuring. These results highlight relaxed constraint on two subunits that may provide opportunities for evolutionary diversification within a conserved family of multimeric receptor assemblies. Such diversification has not been observed in vector species displaced by An. funestus in indoor residual spraying areas, and the potential implications for reduced sensitivity to neonicotinoids are discussed.
Birds and mammals are shrinking and shapeshifting as global temperatures rise. Ecogeographic rules predict that such changes should ease heat stress by increasing surface-area-to-volume ratios, and thus, the capacity for heat exchange. This has led to the hypothesis that body size reductions are driven by thermoregulatory selection or adaptive plasticity, although recent syntheses point to more complex, multifactorial causes. Crucially, recent theoretical models predict that thermoregulatory benefits of smaller body size only emerge at extreme deviations from average phenotypes. Here, we exploit agricultural selection in Japanese quail to directly test this hypothesis, using three breeds spanning extreme differences in body mass, surface area, and relative appendage lengths. Evaporative cooling capacity and the scope for evaporative water loss broadly followed allometric predictions when contrasting small and larger breeds. As expected, this allowed the smallest breed to tolerate higher air temperatures. However, differences in heat tolerance limits between breeds were consistently much smaller than predicted. Additionally, the breadth of thermoneutral zones overlapped in full, and upper critical temperatures were remarkably similar, between breeds. Together, these results show that heat tolerance is only weakly linked to surface-area-to-volume relationships and cannot be explained by size alone. Thus, although smaller bodies may modestly enhance heat dissipation when size variation in a population is substantial, our findings suggest that recent body size reductions and morphological shifts are unlikely to be driven primarily by thermoregulatory benefits.
Although resources are typically distributed continuously in space, species distributions often organize into discrete clusters. In his seminal paper, Turing demonstrated that such clusters can spontaneously arise in population densities, even when populations evolve in environments with continuously varying conditions. This phenomenon is known as Turing instability. In this work, we focus on two models grounded in population dynamics: a one-dimensional model based on the nonlocal Fisher-KPP equation, and a two-dimensional model involving an environmental gradient. We show that phenotypic clusters (sometimes referred to as "species") emerge in these models. We prove that they do not emerge because of Turing instability, but because of stochasticity, and that they disappear when stochasticity is reduced. First, for both models, we start our simulations with initial populations uniformly distributed in the state space. We show that phenotypic clusters quickly emerge and that the distances between them depend on the population size, that is, on the degree of stochasticity. Next, we start from already clearly defined phenotypic clusters. We identify three regimes in the connection between population size, the initial distances between clusters, and the distances between clusters at equilibrium. Last, on the two-dimensional model, we relax the hypothesis of complete clonality by varying the effective recombination rate, explore its effect on phenotypic clustering, and show that phenotypic clustering decays drastically with slight recombination.
Accurate and reproducible assessment of foliar disease severity is essential for evaluating the performance of heterogeneous plant communities and understanding host-pathogen interactions. However, traditional visual scoring methods remain subjective, with limited precision, and difficult to scale in large phenotyping experiments. Here, we present a semi-automated image analysis workflow designed to quantify multiple foliar disease symptoms simultaneously on wheat flag leaves sampled from varietal mixtures. The workflow combines three methodological components: (i) a standardized protocol for leaf sampling and imaging, (ii) supervised machine learning segmentation using Random Forest implemented in Ilastik to classify multiple symptoms (powdery mildew and yellow rust), and (iii) a graphical user interface facilitating pipeline deployment by non-specialist operators. To evaluate the influence of image representation on classification performance, four color spaces (RGB, HSV, HLS, LAB) were systematically compared. The approach was validated using images of durum wheat flag leaves collected from a field experiment assessing eight-way varietal mixtures under natural fungal pressure. Cross-validation against manually annotated images demonstrated high segmentation accuracy across all symptom. Comparison among color spaces revealed only minor differences in performance. Overall, this workflow offers a cost-effective, annotation-efficient and reproducible alternative to deep learning approaches, leveraging open-source and actively maintained tools while requiring limited training data and enabling objective, reproducible and scalable disease phenotyping.
RNA-seq experiments routinely identify thousands of differentially expressed genes, but translating these into biological insights and therapeutic hypotheses often requires integrating multiple tools. Existing web platforms such as iDEP, NetworkAnalyst, and GEPIA2 address individual steps, differential expression, network visualization, or TCGA queries, but lack a unified environment spanning raw data processing to clinical and pharmacological interpretation. TransXplorer (https://www.transxplorer.org) is a freely available web platform that addresses this limitation by integrating the complete RNA-seq analytical workflow. It supports processing from raw FASTQ files using HISAT2 or Salmon, as well as direct GEO dataset import with automated metadata handling. Differential expression analysis is implemented via DESeq2, edgeR, and limma-voom, followed by functional enrichment across more than 1,800 species using Bioconductor resources. Batch effects are automatically detected and corrected using a composite of PVCA, kBET, and Silhouette metrics without requiring predefined batch annotations. Downstream analyses include co-expression network construction (WGCNA), protein-protein interaction mapping (STRING), cell-type deconvolution, and transcription factor inference using integrated DoRothEA and TFLink resources. The platform further links gene signatures to drug candidates through DGIdb and OpenTargets and enables survival and tumour-normal comparisons across TCGA cohorts. Application to cardiac endothelial differentiation (GSE151427) and kidney renal papillary cell carcinoma (TCGA-KIRP) datasets demonstrates accurate batch correction, biologically consistent pathway enrichment, recovery of expected cell-type proportions, and identification of clinically relevant genes and drug candidates. TransXplorer is freely available without a login.
Anxiety has been extensively studied in relation to memory, yet its dynamic association with spatial episodic memory in naturalistic clinical settings remains largely unexplored. We developed an anxiety-spatial-memory EMA protocol (asm-EMA) and deployed it in 30 epilepsy patients undergoing inpatient EEG monitoring, delivering combined momentary anxiety ratings and a validated spatial memory task pseudo-randomly every 90-150 minutes across multiple days. Subject-level asm-EMA means and session-to-session variability both correlated significantly with standard neuropsychological assessments, supporting the clinical validity of our design. Elevated within-person STAI-6 was selectively associated with faster retrieval responses, yet spatial memory accuracy was independent of all three anxiety measures, suggesting a shift in response strategy rather than memory impairment. Within-day anxiety showed short-term carryover between consecutive sessions, with little persistence beyond the next session. The asm-EMA protocol provides a feasible, autonomous framework for capturing moment-to-moment anxiety-memory dynamics in naturalistic settings.
Chromosome-scale assemblies are increasingly available for non-model organisms, but functional annotation remains limited when deep evolutionary divergence erodes primary amino-acid sequence identity even though protein structural similarity can remain conserved. We present a hybrid annotation framework that decouples gene-model discovery from cross-species similarity assignment by combining Evo2-based ab initio prediction of exon-intron structures with ESM-2 protein-embedding-based structural similarity mapping. Applied to the sea lamprey, the framework derives high- or medium-confidence cross-species similarity assignments for 73,485 Evo2-derived translated protein models, including 35,395 high-confidence calls, and expands the deduplicated structural catalog to 31,286 loci, including 20,871 additions absent from the Ensembl baseline. A joint alignment-structure classification identifies 21,391 structurally supported catalog loci that a fixed human DIAMOND protein search does not confidently assign on its own, including 21,184 loci with no detectable human protein-sequence match and 207 loci with only low-confidence matches in the classical 20-30% amino-acid-identity twilight zone. These rescue-space totals describe catalog loci rather than validated one-to-one human-absent genes. In a single-cell RNA sequencing application, a stricter UTR-aware Ensembl+Evo2 reference improves gene recovery and expands the interpretable feature space of the lamprey immune compartment relative to the Ensembl baseline. This enables more resolved annotation of four transcriptionally defined immune cell states, including VLRA+-associated T-like and VLRB+-associated B-like programs together with oxidative iron-handling and iron-associated VLR-linked states. Together, these results show that structural protein signal often persists beyond the limits of pairwise sequence alignment and that an embedding-based annotation layer can extend that signal to improve downstream comparative and…
Terpenes constitute the largest and most structurally diverse class of plant secondary metabolites, with critical roles in plant-environment interactions and broad industrial applications. Although nuclear genome engineering of terpene pathways has been extensively explored, chloroplast genome engineering remains largely undeveloped, with all reported studies restricted to the model plant Nicotiana. Here we report successful chloroplast genome engineering for diterpene production in the crop plant potato (Solanum tuberosum). First, we identified the trnT/trnL plastomic locus as optimal for minimizing integration-associated growth penalties. Insertion of a bifunctional diterpene synthase gene into this plastomic site yielded transplastomic plants with successful diterpene production, but with reduced growth. The co-expression of a geranylgeranyl diphosphate synthase gene to enhance precursor supply restored normal growth while elevating diterpene accumulation. Transplastomic plants were otherwise agronomically comparable to wild-type. This work expands chloroplast engineering as a viable strategy for terpene pathway engineering in crop improvement and high-value terpene production.
The hippocampus is organized along dorsal-ventral and left-right axes, but whether and how these axes interact within defined neuronal populations across behavioral states remains unresolved. Here, we combined within-animal slice electrophysiology with dual-site fiber photometry to compare dorsal and ventral CA1 activity across contralateral hemispheric configurations in mice expressing CaMKII-jGCaMP8s and SynI-jRCaMP1b at distinct longitudinal sites. Ventral CA1 pyramidal neurons exhibited greater intrinsic excitability and stronger AMPAR-mediated synaptic responses than dorsal CA1 neurons. In vivo, CaMKII-defined pyramidal recordings during home cage rest revealed a left-biased event-rate asymmetry within dorsal but not ventral CA1, with no comparable asymmetry in pan-neuronal SynI recordings. Apparent dorsal-ventral differences in spontaneous event rate were therefore configuration-dependent and resolved into a hemispheric, cell-type-specific effect restricted to the CaMKII-defined population. Lead-lag analysis showed that dorsal-ventral temporal coordination was likewise reorganized across configurations and was restricted to pyramidal-cell-biased recordings. During open-field center entries, dorsal CA1 was preferentially recruited before entry across both configurations, whereas non-coordinated entries revealed a relative post-entry suppression of contralateral ventral CA1. Together, these findings suggest that dorsal-ventral CA1 organization cannot be inferred from hemisphere-pooled designs and identify a pyramidal-cell-specific left dorsal CA1 asymmetry as a structural feature that shapes both spontaneous activity and behaviorally driven recruitment along the longitudinal hippocampal axis.
Rapid connectivity alterations of thalamic nuclei during initial learning of goal-directed behaviour
The thalamus is essential for learning, dynamically engaging with other subcortical and cerebral cortex regions throughout the learning process. Here, the thalamus serves as a critical connector hub and synchroniser within the thalamocortical system of the brain. However, whilst higher order thalamic nuclei are known to be particularly important for this process, the exact contributions of individual higher order and first order thalamic nuclei, alongside their individual involvement with cortical networks and subcortical regions, remains unexplored within the initial phase of learning. In light of this, we analysed fMRI data obtained within a paradigm which is designed to examine initial learning processes within feedback-driven stimulus-response learning, in order to explore thalamic contributions. We investigated dynamic learning-related functional connectivity alterations between various thalamic nuclei with other subcortical regions and cortical networks. Our results show that the initial phase of learning was associated with: (1) decreasing functional connectivity between thalamic nuclei and frontoparietal and cingulo-opercular networks, (2) increasing functional connectivity between thalamic nuclei with default mode and salience networks, (3) decreasing functional connectivity between thalamic nuclei and the putamen, and (4) decreasing functional connectivity amongst higher order thalamic nuclei. Furthermore (5) these dynamic alterations were associated primarily by mediodorsal thalamus. Altogether, these results indicate that higher order thalamic nuclei play a crucial role within initial learning and in the generation of novel goal-directed behaviour. This was demonstrated through enhanced functional connectivity with selected cortical networks which drive goal-directed behaviour, alongside decreased functional connectivity with striatal regions which drive motor selectivity.
Coxiella burnetii is the only member of the order Legionellales known to primarily infect vertebrates. The Q fever pathogen is also unusual in that it replicates within an acidified phagolysosome-like vacuole. The evolutionary origins of the virulence determinants underlying this lifestyle remain unclear. More broadly, little is known about how virulence-related traits arise in specialized intracellular lineages, where access to foreign-origin DNA may be more episodic. To address this question, we used Legionellales-wide comparative phylogenomics to reconstruct the gain and loss of traits affecting host interaction, immune evasion, intracellular survival, and metabolism. We found that many virulence-associated traits in C. burnetii predate the modern pathogen and were assembled stepwise in ancestors that likely occupied niches distinct from the acidified vacuolar niche of modern C. burnetii. The common ancestor shared with soft-tick Coxiella endosymbionts likely encoded most C. burnetii type IVB secretion system effectors, indicating that much of the host-manipulation repertoire in C. burnetii was already present before the emergence of the modern pathogen. Distinctive lipopolysaccharide features associated with immune evasion also appear to have accumulated progressively within the Coxiella lineage, including genes implicated in synthesis of virenose, a unique O-antigen sugar critical for C. burnetii virulence. Traits likely to support replication in the acidic Coxiella-containing vacuole likewise accumulated gradually, with generalized stress-tolerance functions predating acquisition of an Mrp cation/proton antiporter that may have further supported pH homeostasis. Additional changes in sugar transport and catabolism, glycolytic control, and respiratory metabolism may have enhanced metabolic flexibility and access to diverse substrates in this nutrient-rich niche. Together, these findings support a model in which vertebrate pathogenicity in C. burnetii emerged …
Microbes that remain uncultivated occupy nearly every ecosystem on the planet; this is particularly true in soils, where despite their prevalence, the roles of rarely cultivated microbes in driving biogeochemical cycles and ecosystem function remain poorly explored. We combine metagenome-informed substrate selection with enrichment sub-communities to generate reduced-complexity communities that preserve co-occurrence and expand experimental access to underrepresented soil lineages without requiring prior isolation of each member. Carbohydrate-active enzyme (CAZyme) profiles from soil-derived genomes were used to select carbon compounds predicted to enrich difficult to culture taxa, including members of the phylum Acidobacteriota. Based on 16S rRNA amplicon sequencing, we reproducibly enriched Terriglobus (Acidobacteriota) on multiple metagenome-guided substrates. Select communities with consistent presence and varying abundance of Terriglobus were passaged in a longitudinal design to generate 89 metagenomes; genus-level profiling revealed that community composition varied between biological replicates but remained consistent within replicates over time, providing diverse Acidobacteriota-containing configurations for downstream analysis. Association network inference identified a core set of co-occurring taxa that positively tracked with Terriglobus across the longitudinal series. In parallel, the substrate-guided approach led to isolation of a novel Terriglobus species, the first cultured representative of its GTDB species cluster. Together, these results establish a generalizable strategy for generating communities enriched with rarely cultivated taxa, yielding tractable systems for studying microbial interactions and community assembly in soil.
Course-based undergraduate research experiences (CUREs) can expand undergraduates' access to research and motivate students to stay in science. Yet, little research has examined how CURE instruction shapes student motivation. We leveraged a motivation-related characterization of non-content talk of 48 CURE and non-CURE instructors to predict the motivation-related outcomes of 462 students. We fit a series of multi-level models (MLM) in which we regressed students' post-course scientific self-efficacy, task values, scientific identity, and science-related intentions onto instructors' self-efficacy and task values-related talk, controlling for students' pre-course levels. We also fit an MLM to explore whether instructors' relationship-building talk (immediacy talk) was associated with students' rapport with their instructor. Instructors' self-efficacy talk did not affect students' self-efficacy, and instructors' immediacy talk had a marginally positive but non-significant association with students' rapport ratings. Instructors' task values talk positively influenced students' scientific identity and some but not all of their task values. Instructors' task values talk also positively influenced students' intentions to pursue a science career, but not graduate education or research careers. Collectively, these results suggest that instructors' task values talk may underpin some of the motivational effects of CURE instruction, but that task values talk need not be limited to CUREs.
Metatranscriptomic (MTX) sequencing enables profiling of gene expression across microbial communities, providing a framework for linking genetic potential with functional activity. However, standard pipelines report normalized abundances rather than raw counts, limiting the use of count-based RNA-seq methods, while Gaussian-based alternatives rely on transformations and assumptions that are often poorly suited to MTX data. We propose a new modeling framework for differential expression analysis of MTX data, built on a scale mixture of exponential distributions, that incorporates DNA abundance to adjust for genomic potential, accommodates subject-specific random effects, treats zeros as left-censored, and employs a mixture prior to handle extreme sparsity. Applied to the IBDMDB multi-omics cohort, differential expression results vary substantially across models, including among Gaussian approaches with different pseudocount choices. Our approach identifies a distinct subset of candidate genes not detected by existing Gaussian methods; these may provide useful leads toward a novel understanding of transcriptomic patterns associated with dysbiosis in inflammatory bowel disease. Estimated dysbiosis effect directions are consistent between our model and Gaussian-based approaches, while effect sizes from our model tend to be larger in absolute value.
Advancing the utility of plant synthetic biology requires the continued development of protein engineering tools. Self-assembling protein compartments, such as virus-like particles (VLPs), provide versatile scaffolds for synthetic biology. However, few plant-expressed VLPs have demonstrated broad amenability to protein engineering, restricting their applications to specific contexts. Here, the Salmonella typhimurium bacteriophage P22 VLP is explored as a novel protein scaffold for plant synthetic biology, demonstrating its application in a eukaryote for the first time. Through transient expression in the biofactory plant Nicotiana benthamiana, the capacity for P22 VLPs to correctly assemble and selectively encapsulate recombinant protein cargo is demonstrated. The durability of this protein scaffold is explored, through co-encapsulation of multiple cargo protein species and by encapsulation through direct fusion to the P22 coat protein. Finally, the ability to simultaneously program cargo encapsulation and external protein display on P22 VLPs in vivo is demonstrated through SpyTag/SpyCatcher-mediated protein conjugation. This work demonstrates the broad utility of P22 VLPs as nanoscale protein scaffolds for plant synthetic biology. Keywords: protein scaffolds, cargo encapsulation, protein display, SpyTag/SpyCatcher, transient expression, Nicotiana benthamiana.
Protein labelling by covalent attachment of a specific substrate to a self-labelling protein tag has become a regular in the life sciences. Herein, we report the design of a two-component labelling system, comprised of a non-fluorescent difluorinated xanthene, called F2X, and a HaloTag mutant engineered for targeted reactivity towards F2X. Upon primary covalent locking of the ligand at the canonical aspartate residue, two proximal lysine residues located at the protein surface can undergo nucleophilic aromatic substitution with the F2X core, building a fluorescent rhodamine via triple-covalent fusion. We used a generalizable in silico pipeline for heuristic conformational sampling of covalent protein-ligand complexes to find suitable mutation sites, culminating in the curation of 7 double-lysine HaloTag mutants for targeted in vitro testing. Reaction with the best-performing mutant, HTPL161K_Q165K, is characterized by full protein mass spectrometry, fluorescence polarization fluorescence lifetime, and fluorescence anisotropy and rationalized by computational modelling. We showcase the system in single molecule microscopy, where obviation of post-labelling purification is a prime advantage when targeting recombinant proteins that may not be expressed in larger quantities, and employ F2X in living cells with reduced photobleaching. Lastly, a cell-impermeable version was obtained by means of sulfonation, exclusively targeting extracellularly exposed HTPKK fused to the neuromodulatory G protein-coupled receptor metabotropic glutamate receptor 2.
WD40 domains share a widespread {beta}-propeller fold, and often act as versatile scaffold proteins. Despite their central role in organizing dynamic cellular complexes, the molecular and structural mechanisms of many WD40 proteins remain poorly understood. Among them, DCAF7, an ubiquitously expressed and essential gene in human, also encodes a highly conserved WD40 protein in eukaryotic organisms. It is known to interact with multiple and functionally diverse partners to coordinates cellular activity of several protein kinases as well as transcriptional regulators, thereby modulating key cellular processes such as cell growth, differentiation, and transcriptional regulation. However, the precise mode of action of DCAF7 is unknown and its important divergence in sequence from better characterize WD40 prevent information transfer by similarity. Structural interactomic can reveal how protein-protein interactions (PPIs) occur within an organism and are essential for understanding biological functions and developing new therapeutic strategies. Using SLiMAn2, AlphaFold2/3 and PSSMsearch, we identified a conserved -helical short linear motif (SLiM) in several well known DCAF7 partners that binds to the top surface of its {beta}-propeller. This motif was subsequently used to generate a regular expression, to identify potential new direct binders across the DCAF7 meta-interactome and the human proteome. Domain-domain interactions were also predicted for some other partners. Finally, modeling of oligomeric complexes with such new hits reveals the structural basis of DCAF7 scaffolding, with links to neurodevelopmental disorders such as autism.
Conventional diagonal stride skiing traditionally includes a glide phase, characterised by a period of relatively passive gliding on one ski. While the glide phase may take advantage of low ski-snow friction, it does not exhibit the same whole-cycle mechanical energy fluctuations seen in running or walking on foot. A new sub-technique, known as running style, substantially reduces the glide phase and may alter the role of elastic tissues, making the movement pattern more similar to uphill running on foot in its temporal organisation. We examined knee extensor and plantar flexor muscle-tendon behaviour in eight competitive skiers performing conventional diagonal and running techniques on a treadmill inclined at 10 deg. Using synchronised ultrasonography, 3D kinematics, ski forces and EMG, we quantified gastrocnemius medialis and vastus lateralis fascicle and muscle-tendon unit (MTU) dynamics in both the running (RUN) and conventional (CON) styles. Shorter glide and total cycle durations during RUN shifted MTU peak length and velocity earlier during the kick phase. Fascicles in both muscles operated at similar velocities across techniques, showing MTU-fascicle decoupling. Vastus lateralis fascicles shortened at higher absolute peak velocities than gastrocnemius in both conditions, while normalised velocities were similar. RUN increased preactivation and advanced EMG timing, while integrated EMG during the kick was lower compared to CON. These findings suggest that, despite large shifts in external mechanics between glide-based and more running-like skiing, elastic tissues may help stabilise fascicle behaviour and preserve a similar contractile strategy across muscles and techniques.