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.
Science Journals
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.
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.
Gene duplication is a major driver of evolution, yet how it generates fundamentally new molecular functions remains poorly understood. Here, we show how such novelty arose in KLMT-1, a selfish toxin that causes genetic incompatibilities in Caenorhabditis tropicalis. KLMT-1 evolved via duplication of an essential tRNA synthetase but, strikingly, lost its ancestral role in tRNA biology and translation. Instead, KLMT-1 localizes to centrosomes, where it targets Aurora kinase A (AIR-1). This innovation is mediated by a three-amino acid insertion that extends a beta-hairpin loop, enabling electrostatic interaction with a regulatory interface on the kinase. Our results demonstrate how changes in selective pressure, combined with minimal modifications in neutrally evolving regions, allow duplicated proteins to access new functional space and evolve entirely new molecular activities.
LD Score Regression (LDSC) is a prominent method, which estimates whole-genome SNP heritability from summary statistics via the slope of a linear regression of GWAS test statistics corresponding to a trait of interest against LD scores. It was claimed by the LDSC authors that the free intercept in the regression accounts for confounding bias such as population stratification. In this study, we argue that the intercept in LDSC must be fixed to 1 for accurate SNP heritability estimation. We show both theoretically and with simulations that the estimated intercept does not accurately capture population stratification effects, and that it adversely affects the accuracy of the heritability estimate introducing bias and increasing variance. Fixing the intercept to 1 eliminates bias and reduces variance when no population stratification is present. On the other hand, under population stratification, LDSC is biased with both the free and the fixed intercept. Additionally, we show that estimated standard errors in LDSC are underestimated, potentially leading to false-positives in downstream GWAS analyses.
Maximum utilization of existing genetic variability in a breeding program depends on the efficient classification of the inbred lines into heterotic groups, particularly under stress conditions. This study applied practical breeding approaches to determine the mode of genetic inheritance for Striga resistance and proposes a weighted heterotic grouping method based on the general combining ability of multiple traits (WHGCAMT) and compares its effectiveness with other existing methods in classifying the inbred lines into heterotic groups in Striga-infested and optimum environments. Using Diallel design IV, 300 crosses were generated from 21 inbred lines and 4 standard testers. The crosses, along with six checks, were evaluated in an 18 x 17 alpha lattice design with two replications at two locations, in both artificial Striga-infested and Striga-free environments. The inbred lines were genotyped using DArTtag SNP markers. Phenotypic and genotypic data were analyzed using R. Analysis of variance revealed significant mean squares for hybrid, general combining ability (GCA), specific combining ability (SCA) and their interactions with environment. Significant positive and negative GCA and SCA effects were detected for grain yield and other measured traits. However, a larger proportion of additive gene action than non-additive gene action was observed for grain yield and most measured traits. The analysis of molecular variance also showed substantial genetic differences within and between clusters. Except for HSCA, the mean grain yield between the inter-group and intra-group hybrids was significant for each method. Pairwise comparison of the inter- and intra-group hybrids of all the methods showed significant differences between the WHGCAMT and all other methods in most cases. WHGCAMT consistently produced higher-yielding inter-group hybrids and lower-yielding intra-group hybrids, achieving breeding efficiency improvements of 55.8%, 4.3%, 15.7%, and 11.4% over the HSCA…
Gross chromosomal rearrangements are a hallmark of many diseases and cancers. The study of their biogenesis and the mechanisms underlying their formation is greatly facilitated by the availability of genetic reporter assays in model organisms. We present here a novel GCR assay developed in fission yeast, a highly relevant model for understanding genome instability related to human biology. The reporter employs canavanine counter-selection to detect GCRs within a chromosomal context. Using this assay, we identified natural hotspots for GCRs, including inverted long terminal repeats (IR-LTRs). Structural analysis of GCR events showed that IR-LTR-induced GCRs mainly result in either terminal deletions with adjacent inverted duplications or repair via long-range break-induced replication (BIR). Deleting IR-LTRs reduces the GCR rate and reveals another hotspot driven by BIR between homeologous aldo/keto reductase genes on opposite arms of chromosome I. This is the first evidence that BIR can occur in S. pombe on long tracks reaching up to 600 kb. Besides highlighting genome rearrangement hotspots, the assay also identifies regulators of genome instability in fission yeast. Loss of Nup132, a component of the nuclear pore complex, increases IR-LTRs-induced GCRs, while the budding yeast homolog Nup133 has no effect on the stability of a structurally similar IR. In contrast, disrupting djc9, which encodes a conserved histone H3-H4 binding protein, decreases GCR rates. Overall, this sensitive GCR assay enables the identification of factors that control spontaneous and fragile motif-induced chromosomal instability, including those conserved in humans but lost through evolution in other organisms.
LINE-1 retrotransposons are the only autonomous mobile elements still active in human genomes and remain a potent source of mutation, genome remodeling, and disease risk. However, young, full-length, potentially active copies (the elements most likely to shape present-day genomes) have been largely inaccessible to population-scale analysis because they are long, repetitive, and poorly resolved by short-read sequencing. Here, we use 47 phased long-read assemblies from the Human Pangenome Reference Consortium, representing 94 haplotypes, to build an allele-resolved view of recent human LINE-1 evolution. We identify 13,617 LINE-1 alleles with intact ORF1 and ORF2 across 683 unique insertion sites, revealing that every genome carries a distinct repertoire of potentially active source elements. These intact LINE-1 profiles recapitulate broad human population structure while exposing a large, rare, and population-enriched reservoir of mobile-element diversity missed by single-reference approaches. We also resolve a structurally variable chromosome 11 LINE-1 array, demonstrating that local duplication and rearrangement can amplify LINE-1 sequence independently of canonical retrotransposition. By comparing full-length LINE-1 sequences, we define activity signatures that separate ancient remnants from recently expanding lineages and uncover young LINE-1 groups whose activity is not fully explained by canonical subfamily labels. Sequence-network analyses further reveal a dynamic history of lineage turnover, in which successful source elements rise, seed new insertions, and are replaced by descendants marked by specific nucleotide changes. Together, these data transform human LINE-1s from a repetitive background into a resolved evolutionary system, linking insertion polymorphism, coding potential, population history, and recent retrotransposon adaptation. Our findings establish the human pangenome as a framework for discovering active source elements and for testing how mob…
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.
Abstract As part of preparedness activities supporting pathogens classified under the UK High Consequence Infectious Diseases (HCID) framework, we previously evaluated both a whole-genome tiling amplicon sequencing scheme and a pan-viral hybridisation capture approach (TWIST-CVRP) for sequencing Andes virus (ANDV). In light of the recent outbreak, we make available viral sequencing datasets generated using a historical ANDV isolate (Chile, 1997). In addition, we provide an evaluation of tiling amplicon scheme performance and present recommended primer updates informed by in silico comparison with the recently released outbreak genome. These datasets are intended to support benchmarking, validation, and optimisation of bioinformatic pipelines across the community.
Phyllanthus niruri (Phyllanthaceae) is a medicinally important herb known for producing phyllanthin, a bioactive dibenzylbutane lignan with reported hepatoprotective and antioxidant properties. However, the biosynthetic basis of phyllanthin production remains unresolved, largely due to the absence of a reference genome for the species. We report a chromosome-scale assembly of P. niruri generated by integrating PacBio HiFi long reads and Illumina short reads, followed by reference-guided scaffolding against Phyllanthus cochinchinensis. The assembly has an L50 of 7 and 97.6% BUSCO completeness. Annotation predicted 19,254 protein-coding genes, of which 91.1% were functionally annotated, with phenylpropanoid biosynthesis emerging as the most enriched specialized-metabolism pathway in the genome. Using pathway-guided genome mining, structural similarity analysis, and comparative metabolic reconstruction, we propose a putative biosynthetic pathway for phyllanthin originating from the phenylpropanoid-lignan branch through secoisolariciresinol-like intermediates followed by terminal O-methylation reactions. A total of 305 unique candidate genes associated with the proposed pathway were identified, including expanded families of dirigent proteins, peroxidases, secoisolariciresinol dehydrogenases, and O-methyltransferases. Comparative transcriptomic analyses across related Phyllanthus species further supported the proposed pathway through coordinated expression of lignan-associated genes and tissue-specific enrichment of O-methyltransferases. This work provides the first reference genome for P. niruri and a prioritized candidate gene set for functional characterization of phyllanthin biosynthesis.
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.
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.
Measles virus remains a significant global health threat, and despite the availability of an effective vaccine, measles cases continue to increase worldwide in recent years. Genomic surveillance has become an essential tool for monitoring virus circulation and investigating outbreaks. Here, we describe a wet laboratory method for whole genome sequencing of measles virus using a tiled amplicon approach and Illumina sequencing technology. A previously published Oxford Nanopore based tiled primer scheme was adapted to include both circulating measles genotypes and for use on the Illumina platform. Two Illumina library preparation kits, Illumina DNA Prep (IDP) and Nextera XT (XT), were evaluated for performance. The IDP kit demonstrated more complete genomes and consistent genome coverage compared with XT. Using quantified reference genomes, the limit of detection was determined to be 10,000 genome copies for genotype B3 and D8. Sequence accuracy was evaluated using previously characterized clinical samples and showed high concordance. This method provides a reliable and sensitive approach for measles virus whole-genome sequencing using Illumina platforms and is suitable for genomic surveillance applications.
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.
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.
Conformational plasticity of RNAs plays important roles in recognizing RNA-binding proteins, and is often modulated by their binding partners. Here, we investigate RNA conformational preferences in a non-redundant dataset of 263 protein-RNA complexes to characterize the structural landscape associated with protein recognition. RNA dinucleotide segments are analyzed using seven backbone torsion angles ({delta}1, {varepsilon}1, {zeta}1, 2, {beta}2, {gamma}2, and {delta}2), two glycosidic torsion angles ({chi}1 and {chi}2) and the pseudo-torsion angle . Focusing on dinucleotide steps present in both interface and non-interface regions, we performed density-based clustering using selected backbone torsion angles to identify recurrent conformational states. We identify 28 distinct RNA dinucleotide conformers containing at least ten members each. Among these, eight conformers represent previously unreported nucleotide conformers (NtCs), including the transitional and the non-canonical states AB06, AB07, BB21, BB22, OP32, OP33, IC08 and IC09. Several of these conformers are preferentially enriched at protein-binding interfaces, suggesting their involvement in local conformational adaptation during protein-RNA recognition. The newly identified conformers span transitional A-B geometries, distorted B-like states, open conformations and compact intercalated structures, highlighting the remarkable structural plasticity of RNA in ribonucleoprotein complexes. Overall, this study expands the current understanding of RNA conformational space and provides a refined RNA dinucleotide conformer library for protein-RNA complexes. These findings will facilitate the identification of novel RNA structural motifs and improved RNA structural modeling, docking protein-RNA complexes and deep learning-based prediction frameworks for describing RNA tertiary structures.
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 …
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.
How do humans store sequences that far exceed working memory capacity? Using visuo-spatial and binary auditory sequences, we previously showed that a Language of Thought (LoT) architecture, in which simple primitives are recursively combined into hierarchical programs, enables efficient storage of structured sequences. Here we ask whether this principle extends to purely ordinal structure: sequences defined by how items repeat and in what order, as in AABBCCAABBCC, independently of their spatial content. Across three experiments, participants reproduced 12-item sequences of spatial locations with various ordinal structures. The minimal description length derived from the LoT model predicted recall accuracy with remarkable precision (r = .96), substantially outperforming Shannon entropy, Lempel-Ziv complexity, chunking models and subjective complexity ratings. Critically, fine-grained analyses of participants' inter-click intervals during reproduction revealed systematic slowdowns at the hierarchical boundaries predicted by the LoT programs, providing a behavioral signature of the underlying mental syntax. These results identify a compact vocabulary of mental primitives, repetition, mirroring, and interleaving, whose composition accounts for the symbolic compression of ordinal structures. For ordinal regularities, human sequence memory operates as a form of program induction, leveraging a domain-general capacity for hierarchical compression to encode complex structured information.
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.
Short-read amplicon sequencing is widely used for fungal surveys but can limit taxonomic resolution. Long-read sequencing enables recovery of the full internal transcribed spacer (ITS) region and may improve ecological and taxonomic inference. Here, we conducted a paired comparison of Illumina ITS2 and PacBio HiFi full-length ITS sequencing using identical DNA extracts from built-environmental air and surface samples (n = 68) collected across homes, a dormitory, and laboratories. Both datasets were taxonomically assigned using the same algorithm and reference database. We performed paired statistics, in-silico ITS2 trimming of long-read sequences, and cross-platform mapping at multiple identity thresholds. Full-length ITS provided higher taxonomic resolution, assigning a greater fraction of ASVs at the family (98% vs. 88%) and species (42% vs. 32%) ranks than ITS2 (paired Wilcoxon q=0.002). Alpha-diversity comparisons showed similar Shannon diversity across pipelines, whereas richness metrics were consistently higher for full-length ITS. Beta-diversity analyses indicated broadly comparable community-level patterns, although full-length ITS revealed stronger sample-type- and location-associated structure (PERMANOVA R{superscript 2} 0.06, p=0.0001). In-silico ITS2 trimming reduced these differences, indicating that amplicon length is a major contributor to enhanced taxonomic resolution and ecological inference. Cross-platform mapping further showed extensive one-to-many relationships between ITS2 and full-length ITS ASVs, consistent with increased sequence resolution in long-read data.Together, these results show that ITS2 sequencing provides robust community-level profiling, while full-length ITS enables improved richness estimates and finer ecological and taxonomic resolution. This paired, bias-aware framework provides a practical template for selecting fungal amplicon sequencing strategies in built-environment mycobiome studies.
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…
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.
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.