The evolution of reproductive isolation is central to speciation, yet the earliest stages of this process remain poorly understood. In particular, it is unclear how rapidly barriers to mating arise during adaptation, whether they accumulate predictably, and how they depend on ecological context. Here, we investigate the evolution of mating efficiency during prolonged asexual adaptation in diploid Saccharomyces cerevisiae. Twelve replicate populations were evolved for 1200 generations in two distinct carbon environments, glucose and galactose, under strictly asexual conditions. At regular intervals, we induced sporulation and quantified mating efficiency using three complementary assays: within-population crosses, crosses between populations evolved in different environments, and crosses between evolved populations and the ancestral strain. We find that mating efficiency evolves during asexual adaptation, with outcomes that depend strongly on the environment. While glucose-evolved populations remain largely stable, galactose-evolved populations exhibit a reversible decline. Overall, changes in mating efficiency are dynamic, heterogeneous, and often transient, with evidence for both intrinsic reductions in mating competence and context-dependent incompatibilities between populations. Together, these results show that asexual adaptation can generate rapid but non-monotonic changes in mating compatibility. Early reductions in mating efficiency are heterogeneous, environment-dependent, and often reversible, and do not accumulate into stable reproductive isolation over the timescale examined. Our findings suggest that the initial stages of divergence are characterized by dynamic and contingent perturbations of reproductive traits, rather than a steady progression toward speciation.
Science Journals
Strong population contractions can leave a persistent genomic legacy that can influence populations long after their demographic recovery. While bottlenecks facilitate the removal of strongly deleterious mutations, the effectiveness of purging may be limited in historically small populations. The Kirtlands warbler (Setophaga kirtlandii) is a rare North American songbird with an ancestrally small population. After narrowly evading extinction, they are one of few species that have been delisted from federal protections in the USA. Despite their recovery, a previous study showed evidence for recent inbreeding and a high burden of deleterious mutations that may have not been purged despite strong bottlenecks. Historical DNA offers a unique opportunity to understand the consequences of recent demographic declines on genetic diversity. Here, we use DNA from over 100-year-old museum specimens to estimate changes in genetic load in the Kirtlands warblers pre- and post-bottleneck. We validate our results with forward-in-time genetic simulations and explore how sample size and missing data can affect estimates. Both empirical data and simulations suggest a reduced ability to purge deleterious mutations in this historically small population. Our simulations also highlight that limited sampling design and data quality can constrain the ability to detect changes.
Biodiversity is commonly summarized by macroecological mean patterns, most prominently the species-area relationship (SAR) linking habitat area to expected species richness. Yet conservation, policy, and economic decisions increasingly require risk metrics: probabilities of rare but consequential biodiversity shortfalls, including local collapse. Such tail risks are central in finance and insurance but remain difficult to quantify in ecology because the data needed to estimate full richness distributions are rarely available at decision scales. Here we provide a mechanistic route from species-area relationships to biodiversity risk metrics. We show that when regional species abundances are well approximated by Fisher's log-series, a minimal immigration-extinction mechanism yields a closed-form stationary distribution for local richness whose structure tightly couples the mean SAR to richness variability and lower-tail probabilities. This coupling implies exact fluctuation-response identities and an explicit integral transform that reconstructs collapse probabilities and other tail risk measures directly from the mean SAR. These results define ecological analogues of financial risk metrics---such as collapse probability and lower-tail quantiles---without requiring direct estimation of the full richness distribution. Using high-resolution ForestGEO tree censuses spanning tropical, subtropical, and temperate forests, we find empirical support for these predictions across spatial scales. Together, our results show how widely measurable species-area relationships can be elevated from descriptive averages to operational tools for biodiversity risk assessment and reliability-based conservation planning.
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.
Elucidating how habitat degradation facilitates extinction is critical for effective conservation efforts. Here, we propose integrating physiologically-structured population models into stochastic population viability analyses to assess how differing consequences of habitat degradation interact to drive extinction dynamics in a focal population. Using the isolated spectacled caiman Caiman crocodilus population/ecomorph from the Apaporis River as a case study, we find that threatening the resource base, which individuals increasingly rely upon, to outgrow vulnerable size ranges and mature accelerates extinction. We also found that when habitat degradation impacts both the primary adult and juvenile resource bases, this can have marked synergistic effects on threatening population viability. By contrast, destroying nesting sites has only a small effect on accelerating the impact of deteriorating prey availability. Through integrating community-level feedback between habitat degradation/change and population dynamics/structure, our approach provides a comparative framework for assessing the relative importance of distinct mechanisms through which habitat degradation ultimately drives extinction risk.
Wildlife vaccination could become a powerful strategy to mitigate disease-induced biodiversity losses, yet many vaccines for wildlife diseases provide only limited protection. Notably, tools to control the fungal pathogen Batrachochytrium dendrobatidis (Bd) are urgently needed for amphibian conservation. Laboratory experiments have demonstrated that prophylactic exposure to Bd metabolites increases host resistance, significantly reducing infection intensity in amphibians subsequently challenged with live Bd. Because Bd metabolites are non-infectious and applied topically, this treatment has potential to be administered to waterbodies to vaccinate and protect amphibians. We developed an agent-based model that indicated imperfect vaccination could reduce or amplify Bd infections at the population level, depending on degree of enhanced resistance or tolerance. Utilizing a Before-After-Control-Impact design with ten years of data, we conducted an ecosystem-level trial where we applied low levels of Bd metabolites or a sham control treatment to ponds in California and subsequently quantified Bd prevalence and infection intensity in metamorphosing Pacific chorus frogs (Pseudacris regilla). Unexpectedly, infection intensity was significantly greater in treated ponds relative to control ponds following metabolite addition. Additional model simulations indicated that this could occur via two mechanisms: (1) if treatment greatly increased tolerance alone or in combination with smaller increases in resistance, or (2) if a deleterious environmental interaction caused the treatment to increase susceptibility, rather than promote resistance. Future research is needed to determine whether tolerance or environmental factors drove heightened Bd infection intensities in this field trial to identify contexts in which this treatment can be used as a conservation tool.
Cliffs are environmentally extreme yet biodiversity-rich ecosystems that harbour specialist plants, many endemic and threatened. Plant persistence in these nutrient-poor substrates may depend on tightly linked soil- and root-associated microbial communities, which remain poorly understood. These interactions may become increasingly important with the global expansion of recreational climbing. While physical climbing impacts on vegetation are documented, potential chemical effects, from the use of climbing chalk (magnesium carbonate), on soil properties and plant-associated microbiota remain unknown. We sampled soils and roots beneath cliff-specialist and generalist plants, and unvegetated soils, across climbed and unclimbed routes in northern, central, and southern Spain. Soil physicochemical properties were quantified, fungal communities were characterized using ITS-metabarcoding, and structural equation modelling was used to disentangle direct and indirect effects. Climbing increased soil pH and altered soil chemical properties, driving shifts in fungal diversity and functional composition in soil and roots. The relative read abundance of root-associated symbiotrophic fungi declined, whereas arbuscular mycorrhizal fungi and pathogens increased in climbed cliffs. Overall effects were consistent, with cliff-specialist plants mediating nutrient and fungal shifts. Our findings show that climbing can reshape cliff soil chemistry and fungal communities, with potential cascading consequences for plant functional performance, nutrient dynamics, and ecosystem resilience.
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…
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.
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.
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.
Genetic variation is the raw material for evolution. One source of variation is chromosomal rearrangements, which can bring genes together and form genetic linkage. Rearrangements can also suppress recombination and gene flow, as in the case of sex chromosome evolution. We conducted the first population genomic study of the red harvester ant Pogonomyrmex barbatus to investigate genomic rearrangements that differentiate the lineages J1 and J2 in the "dependent-lineage system" (also known as "social hybridogenesis"). In this unusual reproductive system, males and females from different lineages mate to create hybrids, yet these hybrids develop into sterile offspring (workers), and so the two lineages remain reproductively isolated. We sequenced high-quality reference genomes for the two lineages to search for a potential explanation of the suppression of gene flow between them. Comparison of the two genome assemblies revealed multiple large-scale genomic rearrangements, all of which occurred in the J1 lineage. The rearrangements formed some of the largest J1 chromosomes, including the largest scaffold in the assembly that was formed by at least two translocation events and additional intra-chromosomal rearrangements. The translocations brought together 118 odorant receptor (OR) genes on this rearranged chromosome, 44 of which are 9-exon ORs, which are implicated in chemical communication in ants. We also identified an enrichment of transposable elements in a large synteny gap between the translocated segments. The discovery of multiple translocations that formed large rearranged chromosomes provides a potential explanation for the reproductive isolation between the pair of dependent lineages in this system, and opens the way for the study of the molecular genetic basis of an intriguing evolutionary phenomenon in these and in other ant lineages.
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.
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.
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.
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.
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.
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.
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.
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.
Double-stranded RNA (dsRNA) is a potent immunogenic impurity and its detection is a critical quality attribute in characterizing mRNA therapeutics. Standard analytical methods (e.g., sandwich ELISA) are only able to resolve the bulk presence of dsRNA and cannot characterize the different sub-species that may be present within a mRNA sample.. In this study, we use mass photometry (MP) as a single-molecule analytical platform for the simultaneous detection and characterization of dsRNA impurities in mRNA samples. We demonstrate how ionic strength can interfere with the stability of the mAb/dsRNA complex and measure the binding affinity (1 nM) under a set of parameters for reproducible characterization of the complex. We then leverage the J2 antibody to identify antibody/dsRNA complexes that then resolve dsRNA-positive species within an mRNA sample based on discrete molecular weight profiles. Furthermore, we introduce a novel MP assay that harnesses the repulsive surface chemistry of uncoated glass to exclude the bulk mRNA analyte to enable the use of higher loading concentrations to sensitively profile trace dsRNA impurities as antibody-bound species. This work establishes MP as a valuable next generation mRNA analytical tool for analyzing dsRNA byproducts within mRNA samples.