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.
Science Journals
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.
Humans comprehend language incrementally, updating the representation of sentence meaning with each incoming word. These updates are guided by the distance between each perceived word and prior expectations--the prediction error. The alignment between large language models (LLMs) and cortical activity inspires the hypothesis that the cortical computation of prediction error is Surface-based, driven by statistical patterns of word form co-occurrence. In contrast, psycholinguistic models propose that prediction error computation is Meaning-based, driven by word semantics. We used polysemic words with ambiguous semantics to distinguish these models: ambiguity would introduce uncertainty into meaning representations and hence the prediction error, if Meaning-based, but would not affect the prediction error, if Surface-based. We examined how ambiguity influenced prediction error signatures in self-paced reading times and magnetoencephalographic (MEG) neural responses during sentence processing. While an LLM-based proxy of prediction error robustly predicted reading times and neural responses to unambiguous words, it failed to predict either under ambiguity. That is, prediction error computation was altered by uncertainty in word meaning, which supports the Meaning-based model and corroborates the essential role of word meaning in predictive language processing. Our findings highlight an important limitation of LLMs as in silico models of the human language faculty.
Transcriptomics has transformed our understanding of the brain, but assigning transcriptomic identities to neurons recorded in vivo remains challenging at scale. Existing platforms can pair transcriptomic identity with two-photon calcium imaging in small populations of approximately 100 neurons, but they require recorded cells to be sparse and therefore cannot be applied to large population recordings. Here, we present coppaFISH 3D, a spatially resolved transcriptomics method, and CASTalign, an in silico alignment framework, which together enable transcriptomic identification of thousands of simultaneously recorded cells. coppaFISH 3D detects hundreds of genes in thick 50m fixed sections while preserving tissue integrity, enabling both 3D registration to in vivo imaging and integration with immunofluorescence labelling. The platform is fully powered by open chemistry and open source software, runs on commodity hardware, and can be performed at very low cost per section. It therefore enables transcriptomic identification of recorded neurons at scale, making it possible to study how transcriptomic identity shapes activity in neural populations.
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.
WD40 domains are major protein-protein interaction (PPI) scaffolds, yet their contributions to fungal pathogenicity remain poorly defined. We systematically analysed 94 canonical WD40 proteins in Cryptococcus neoformans. Conditional knockdown and sporulation identified 36 essential WD40 proteins, while in vitro and in vivo profiling of 103 signature -tagged deletion strains spanning 52 genes uncovered 31 pathogenicity-related WD40 proteins, including epigenetic and post -transcriptional regulators. We identified Wcp1, a dual-domain protein whose WD40-repeat and cyclophilin domains are required for growth at 37{degrees}C under 5% CO2. Its WD40 scaffold and PPIase domain supported CO2/heat tolerance and virulence. Notably, Wcp1 couples these functions to acidic pH adaptation: wcp1{Delta} failed to grow under elevated temperature and CO2 at acidic pH, exhibited enhanced intracellular acidification, reduced macrophage survival and attenuated virulence in Drosophila and mice. Integrated transcriptomic and proteomic analyses place Wcp1 at the centre of intracellular pH homeostasis, coordinating proton transport, metabolic adaptation and stress-buffering networks.
Staphylococcus aureus encounters diverse environmental conditions during colonization and infection, including fluctuations in nutrient availability, oxidative stress, and oxygen limitation. Adaptation to these environments requires regulatory systems that coordinate stress responses with metabolic remodeling. The extracytoplasmic function sigma factor SigS contributes to stress adaptation and virulence in S. aureus and directly activates expression of the sroAB operon, which encodes the small proteins SroA and SroB. While previous work demonstrated that SroA participates in feedback regulation of sigS expression, the broader physiological role of SroA has remained unclear. To define the regulatory functions of SroA, we performed RNA sequencing following inducible overexpression of sroA in S. aureus. Transcriptome analysis revealed extensive remodeling of gene expression, with approximately 200 transcripts significantly altered. Transcriptome analysis revealed coordinated repression of metabolic pathways (including nitrate respiration and nucleotide biosynthesis) alongside activation of stress-response and nutrient acquisition genes. Northern blot and quantitative RT-PCR analysis confirmed repression of narG and narJ transcripts following SroA overexpression. Consistent with these transcriptional changes, nitrate reduction assays demonstrated that SroA overexpression reduces nitrate respiration activity. In addition to repression of nitrate respiration genes, SroA overexpression broadly suppressed genes involved in de novo purine and pyrimidine biosynthesis. In contrast, transcripts associated with stress responses and nutrient acquisition, including the SOS-associated gene sosA and the phosphate transport gene pstS, were upregulated. Together, these findings identify SroA as a regulator that links stress-responsive signaling to metabolic remodeling in S. aureus, particularly through modulation of nitrate respiration pathways.
Quantifying the lipid biosynthesis rate of archaea in hot spring sediments is necessary to interpret the abundance, isotopic patterns, and environmental significance of archaeal lipid biosignatures, with implications for modern biogeochemical cycling and astrobiology. Here, we performed lipid hydrogen stable isotope probing (LH-SIP) experiments on whole sediments collected from two high-temperature, suboxic, circumneutral hot springs in Yellowstone National Park (USA) and El Tatio Geyserfield (Chile). We determined the incorporation of 2H2O into intact polar lipids (IPLs) which provides a taxon- and metabolism-agnostic quantification of biosynthesis under near-natural conditions. We targeted isoprenoid glycerol dialkyl glycerol tetraether lipids (IPL iGDGTs) and recovered structures with 0 to 7 cyclopentyl rings from both springs. We observed minor 2H-uptake into archaeal IPLs in spring sediments in Yellowstone, corresponding to decadal-scale apparent generation times (16 {+/-} 7 years), and no 2H-uptake in El Tatio sediments (consistent with minimum generation times of 35 {+/-} 5 years). We infer that net production of sedimentary IPL-iGDGTs is very slow, consistent with a combination of slow archaeal growth, persistence of older IPLs, lipid recycling, and/or contributions from recently sedimented planktonic biomass. These are the first direct, ex situ estimates of archaeal lipid production rates in terrestrial hydrothermal systems using LH-SIP incubations and provide critical constraints for interpreting archaeal lipids in ancient hot spring deposits. This research establishes a framework for assessing activity by slow-growing extremophilic archaea in hydrothermal environments and provides support for targeting hydrothermal deposits on Mars for biosignature detection efforts.
Collagens are key components of the extracellular matrix (ECM) that play a crucial role in maintaining structure, strength, and function of the lungs. Fibrillar collagens are crosslinked by enzymes such as lysyl oxidases and transglutaminases and organized into networks by proteoglycans and glycoproteins. Collagens are the main load-bearing components and along with elastin may impart a non-linear strain hardening behavior to the lung. In disease, collagen crosslinking and organization can be disrupted, possibly due to abnormal levels of enzymes or ECM components. Few studies have examined collagen crosslinking and organization in healthy and diseased human lungs. In this study, alterations in collagen crosslinking and organization were investigated in human lung control, fibrotic and chronic obstructive pulmonary disease (COPD) tissue sections. Ultra-performance liquid chromatography and second harmonic generation microscopy measured pyridinoline crosslinks and the distribution of mature and immature collagens within the decellularized scaffolds, respectively. Fibrotic scaffolds had higher total collagen but less crosslinking per mole of collagen compared with COPD donors. Image analysis by second harmonic generation microscopy showed mature collagens populated airway or blood vessel walls in all three groups and in the parenchyma of fibrotic scaffolds. Immature collagens, on the other hand, were mainly localized to parenchymal regions in control and COPD scaffolds, with fewer immature collagens in fibrotic parenchyma. Additionally, quantification of the mature to immature collagen ratio in defined regions of control and diseased scaffolds showed increased organized collagen in fibrotic tissue. Our study shows that collagen crosslinking and organization are disrupted in fibrotic and COPD lungs and these changes may be compartment specific and can contribute to aberrant mechanical properties of diseased lungs. Our findings highlight that along with total collagen…
Ocean warming is altering abiotic environments and biotic interactions experienced by marine organisms, where sensitive early developmental windows occur in biologically complex seawater communities. The impact of these interactions on developmental processes and fitness in hosts is not well understood, but likely contingent on the establishment of a host-associated microbiome. Here, we hypothesize that temperature and microbial exposure during embryogenesis influence larval microbiome assembly and host morphology. Strongylocentrotus purpuratus embryos were raised in low microbial richness (LMR) or high microbial richness (HMR) seawater at ambient (14 {ring}C) or elevated (18 {ring}C) temperature, then collected at 2, 4, and 6 days post-fertilization (dpf) following multiple feedings. Higher microbial diversity was observed in larvae that developed in HMR seawater when compared to LMR. Differences in relative abundances of dominant microbial families between seawater and larvae suggest some degree of host selectivity in microbiome assembly. Temperature did not strongly alter microbiome composition, but both temperature and microbial condition led to differences in larval morphology by 6 dpf, potentially due to enrichment of microbes with chemoheterotrophic functions. By linking how temperature and microbial communities interact with host development, we contribute novel insights into how early-life environmental conditions impact holobiont formation and morphology.
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.
Dendritic arbor morphology is shaped in part by interactions with neighboring dendrites, and its geometry strongly influences the spatial distribution and strength of synapses. These observations raise the possibility that local dendritic contacts help determine where synapses accumulate and strengthen. Previous work in cultured hippocampal neurons showed that dendrite-dendrite contact sites are non-random and associated with local synaptic clustering. Here we asked whether a different type of dendritic contact, formed between a dendrite and the soma of a neighboring neuron, behaves similarly. Using dissociated hippocampal cultures, immunofluorescence imaging, time-lapse microscopy, quantitative image analysis, stochastic spatial simulations, and minimal quantitative modeling, we identified three recurrent classes of dendrite-soma interactions (DSIs): dendrites crossing directly over a neighboring soma, growing tangentially along the soma perimeter, or contacting the proximal region where a neighboring dendrite emerges from the soma. These interactions were abundant, occurred exclusively between different neurons, and showed substantial structural persistence over several days. Their overall frequency exceeded stochastic predictions across culture densities, and two configurations - proximal and tangential contacts - were selectively enriched above random expectation, whereas soma-crossing contacts were largely consistent with stochastic overlap. DSI composition also changed over development, with proximal contacts becoming progressively more prevalent. At DSI sites, synaptophysin-positive puncta were significantly denser and more intense than on non-interacting dendritic segments, consistent with local enrichment and strengthening of presynaptic specializations. Minimal modeling further indicated that biased formation together with developmental stabilization explains the observed organization better than stochastic geometry alone. These findings identify DSIs a…
Despite sharing the same genes and the same environment, individuals often develop substantial phenotypic differences. While this pattern has been documented across diverse species and traits, the processes giving rise to this 'stochastic' or non-shared environmental variation remain unclear. Recent mathematical models of development in which phenotypes are gradually constructed may offer some clues. These models show that imperfect environmental cues can generate striking variation in developmental trajectories and adult phenotypes. At the population level, such imperfect cues produce increasing stability of individual differences across ontogeny (e.g. animal personality) and patterned distributions of mature phenotypes (e.g. normal or skewed) that resemble those observed in real organisms. Our paper synthesizes existing models in which stochastic phenotypic variation arises solely as a by-product of mechanisms missing their phenotypic targets because of imperfect cues. We then link these models to related, but independent, mathematical theory exploring the environmental conditions under which stochastic phenotypic variation is favoured by natural selection. Our integration shows that stochastic sampling is often favoured over classic bet-hedging strategies involving non-plastic generalist or specialist strategies. Our findings provide new directions of research on stochastic sampling as a mechanism for adaptive stochastic variation within and across generations.
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.
Mapping the genetic basis of inter-individual heterogeneity in multifactorial diseases opens the door to mechanistic insights and opportunities for targeted intervention. In Alzheimer's disease (AD), clinical and pathological heterogeneity is well recognized, but genetic dissection is limited by a lack of well-powered cohorts with deep phenotypic characterization. Here, we introduce a polygenic score (PGS) analysis strategy to address these limitations by leveraging the inherent pleiotropy in complex trait genetics. We perform a cross-cohort, cross-trait application of pre-trained PGS, integrating 713 UK Biobank-derived PGS with 36 deep AD phenotypes across 1678 ROSMAP participants. We identify 268 statistically significant (FDR<0.1) associations between 12 prioritized PGS and 36 AD phenotypes. Prioritized PGS include blood lipid measurements, inflammatory biomarkers, and cancer traits; observed AD phenotypes include cognition, amyloid, and tangles. Of the 268 associations, 49 persist with APOE-excluded PGS. Predictive models trained on multiple prioritized PGS outperform the AD PGS or APOE alone for predicting amyloid and cognition. Lastly, our approach identifies six individual-level AD polygenic subtypes supported by distinct pathological patterns. Overall, we combine large-scale biobank resources and deeply-phenotyped cohorts using PGS, reveal genetic features underlying AD heterogeneity, and provide a general model for stratifying heterogeneous disease-focused cohorts using genomics.
How post-mitotic neurons maintain precise transcription factor (TF) levels throughout life remains a fundamental open question. Here, we challenge the prevailing model of positive autoregulation by demonstrating that UNC-3 (Collier/EBF1-4), a dosage-sensitive TF continuously required for cholinergic motor neuron identity in C. elegans, negatively regulates its own expression. Using genetics, biochemistry, and inducible protein depletion, we show this self-repression occurs directly at the transcriptional level and persists beyond development. CRISPR/Cas9 disruption of negative autoregulation causes motor neuron identity and locomotion defects, establishing its functional necessity. Mechanistically, the UNC-3 DNA-binding domain is required and sufficient for self-repression, with an AlphaFold2 screen implicating chromatin factors as interaction partners. Critically, UNC-3 self-repression is continuously counterbalanced by positive input from the HOX cofactor CEH-20/PBX, revealing a dynamic "balancing act" between opposing regulatory inputs that stabilize TF dosage over time. Mutations in the unc-3 ortholog EBF3 cause a neurodevelopmental syndrome, and disease-associated variants disrupt UNC-3 self-repression, revealing a key molecular mechanism underlying the disorder. We propose that negative autoregulation continuously counteracted by positive input represents a broadly applicable principle for maintaining dosage-sensitive TF expression to secure post-mitotic cell identity.
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 …