Biology

Abstract Background Prostate cancer (PCa) diagnosis remains challenged by the limited specificity of prostate-specific antigen (PSA) testing, which cannot reliably distinguish malignancy from benign prostatic hyperplasia (BPH). MicroRNAs (miRNAs) are emerging candidates for liquid biopsy-based diagnostics, but most studies assess expression in isolation within a single compartment (biological source - Tissue, blood, serum, urine etc.), overlooking both compartment-specific behavior and the coordinated relationships among miRNAs. Methods We profiled four candidate miRNAs --- miR-19b-3p, miR-21-5p, miR-101-3p and miR-375-3p, across four biological compartments (prostate tumor tissue, urine, serum, and blood) in 179 patients undergoing prostate biopsy for clinical suspicion of PCa (104 PCa, 75 BPH) using qRT-PCR. Urinary exosomal RNA was isolated with a commercial exosome isolation kit so from here onwards this compartment will be referred to as urine. Differential expression was quantified using Cohen's d; inter-miRNA coordination was assessed via Spearman correlation and differential correlation ({delta} r) analysis; and a compartment-level network rewiring score was derived as the sum of {delta} r| across miRNA pairs. Cross-compartment structural alignment was evaluated by comparing correlation patterns at the population level. Diagnostic models combining PSA, age, and urinary exosomal-miRNA features were evaluated using Logistic Regression, Elastic Net Logistic Regression and Naive Bayes classifiers under leave-one-out cross-validation (LOOCV). Results Effect sizes were largest and most consistent in urine, with miR-101-3p showing the strongest separation between PCa and BPH (d = -1.01), followed by miR-21-5p (d {approx}-0.72$) and miR-19b-3p (d {approx}-0.64). Two markers (miR-19b-3p, miR-375-3p) showed directional reversals across compartments, indicating that disease-associated signals are compartment-specific rather than uniformly conserved. In tumor tissue, PCa was associated with substantial reorganization of inter-miRNA coordination (network rewiring score = 2.46), including the emergence of a strong miR-21-5p--miR-375-3p co-regulatory axis ({delta} r = +0.87$) and decoupling of the miR-21-5p--miR-19b-3p relationship ({delta}r = -0.64$). Urine showed a structurally distinct coordination pattern (rewiring score = 1.77), dominated by a miR-101-3p--miR-19b-3p axis (r = +0.56) absent from tissue; cross-compartment comparison showed concordance in only 1 of 5 miRNA pairs, indicating that urine's architecture is largely independent of tissue's. For diagnostic translation, the conventional PSA cutoff (4 ng/mL) achieved 100% sensitivity but only 23.5% specificity. In urine, miR-101-3p performs better than other miRNAs, with AUC of 0.77 (95% CI: 0.62--0.90). Adding PSA and age to the urinary miR-101-3p further improved discrimination to an AUC of 0.91 (95% CI: 0.82--0.99), with 70% specificity at 92% sensitivity; this pattern was consistent across Elastic Net and Logistic Regression classifiers. Expanding the model to include all urinary miRNAs, age, and pair-derived coordination features did not improve on this result (AUC = 0.88), indicating that population-level coordination changes did not translate into additional individual-level diagnostic value in this cohort. Conclusions miRNA signals in extracellular compartments do not represent direct surrogates of tumor-level molecular architecture; each compartment harbors a distinct, transformed coordination structure reflecting its biological context. While these coordination-level changes are mechanistically informative, the most direct translational gain in this study came from a parsimonious model combining PSA, age with a single urinary marker, miR-101-3p, which improved AUC from 0.77 to 0.91, with specificity 70.5% at 90% sensitivity criteria. This combination represents a promising, interpretable candidate for reducing unnecessary prostate biopsies, pending validation in larger, independent cohorts. Keywords: MicroRNA, Compartment-Specific Biomarkers, Urinary Exosomes, Differential Correlation, Liquid Biopsy, Machine learning, PSA, Early diagnosis

Seasonal breeding is a remarkable adaptive trait, but it constrains efficient production in the sheep industry. Recent studies have shown that seasonal breeding is associated with endogenous circannual rhythms, which are regulated in part by the circadian clock system. FBXL3, a pivotal component of the SCF (SKP1 - CUL1 - F-box) E3 ubiquitin ligase complex, is a known determinant of the mammalian circadian period. In this study, we identified a missense mutation, T183M, in FBXL3 through selective sweep analysis. The allele frequency of this mutation differed significantly between sheep breeds exhibiting year-round estrus and those showing seasonal breeding patterns. The association between the T183M mutation and seasonal breeding was further validated using an ovariectomized, estradiol-implanted sheep model. We then generated mice carrying the homologous T183M mutation and found that they exhibited significantly lengthened circadian periods, accompanied by reduced CRY1 expression and increased CLOCK expression. Co-immunoprecipitation assays confirmed that the mutation reduced the interaction between FBXL3 and CRY1. These findings demonstrate an evolutionarily conserved role of FBXL3 in the circadian clock system. We propose that the T183M mutation disrupts day-length recognition, thereby influencing seasonal estrus in sheep. Author summarySeasonal breeding limits sheep productivity and is regulated by circadian rhythms, yet the key genetic determinants remain poorly understood. Here, we identified an FBXL3 T183M missense mutation whose allele frequency differed markedly between year-round-estrous and seasonally breeding sheep, and validated its association with seasonal reproduction in a sheep population. Functional analyses in mutant mice showed that this variant lengthened the circadian period, disrupted the expression of core clock genes, and weakened the interaction between FBXL3 and CRY1. These findings suggest that the FBXL3 T183M variant impairs day-length perception, thereby modulating seasonal estrus in sheep. Our study reveals a conserved circadian mechanism underlying seasonal breeding and highlights FBXL3 T183M as a promising genetic target for improving reproductive performance in sheep.

This study employed a "one-versus-many" (a single pregnant woman compared with multiple non-pregnant women) urinary proteomic comparative framework to examine whether fetal development-related signals can be captured through changes in the urinary proteome under conditions of limited sample size. The experimental group consisted of urinary proteomic data from three women with normal pregnancies (R6, R15, R16) at three gestational time points ([~]6-8 weeks, 22-24 weeks, and 32-34 weeks; Wang et al., 2022), while the control group consisted of urinary proteomic data from six healthy non-pregnant women (Bao & Gao, ChinaXiv: 202302.00108v2). A total of nine independent "1 vs. 6" differential protein analyses and DAVID GO Biological Process enrichment analyses (P < 0.05) were performed. The results showed that all three pregnant women exhibited a large number of differentially expressed proteins at each time point, and all enriched GO BP terms highly relevant to concurrent fetal organ development (nervous system, lung, eye, ear, kidney, etc.). This suggests that the pregnancy urinary proteome can reflect fetal development signals, corroborating the findings reported by Wang et al. (2025) in a rat model. This study demonstrates that the one-versus-many comparative approach maintains high sensitivity under small-sample conditions and can provide a methodological reference for personalized pregnancy medicine.

Abstract Background: Acute Lymphoblastic Leukemia (ALL) is a highly heterogeneous pediatric malignancy. Despite high survival rates, relapse and the involvement of central nervous system (CNS) remains a significant clinical challenge. Traditional clinical parameters often lack the precision required for early detection and risk stratification. This study utilizes high-throughput proteomics and machine learning to identify molecular signatures in cerebrospinal fluid (CSF) that characterize disease effect and treatment response. Methods: 82 CSF samples from 41 pediatric ALL patients at diagnosis (VD) and remission (VR) were analyzed. Proteomic profiling of 276 proteins was performed using Olink Proximity Extension Assay. Differentially abundant proteins were identified (q-value< 0.05, |Log_2FC| > 0.5) using the Wilcoxon rank-sum test. Three machine-learning algorithms - Random Forest, LASSO, and SVM-RFE - were integrated to select the differentially abundant proteins in VR and VD and between CNS involvement levels. To validate the data Pan-Cancer Atlas analysis was done using two different platforms. Results: In the remission phase, we observed significant alterations in the expression of key proteins compared to diagnosis, with ADGRG1 and KYNU showing a marked increase, while CCL17, CD5, CD27, CXCL9, CXCL11, FASLG, GZMA, and TNFRSF9 were significantly downregulated. Furthermore, our analysis identified distinct protein signatures associated with CNS involvement: CCL4, CTSC, CXCL10, CXCL9, and MMP7 were differentially abundant at the VD stage, whereas CAIX, CASP-8, HAGH, CXCL9, MMP7, MCP-2, and VWC2 at the VR stage. Conclusion: Integrating Olink proteomics with machine learning identified molecular signatures in ALL that have the potential to be further developed to a biomarker panel for monitoring treatment response and guiding personalized therapeutic strategies shifting the focus toward the Precision One Health approaches.

Background: Bladder cancer has entered an era in which immune checkpoint blockade (ICB) and antibody-drug conjugate (ADC)-based combinations are reshaping clinical management. However, transcriptomic scores that connect prognosis, tumor microenvironment state, and treatment response are incompletely defined. Methods: Open-access TCGA-BLCA RNA-seq, clinical, mutation, copy-number, and RPPA data were downloaded from the Genomic Data Commons (GDC). Tumor-normal differential expressions, survival screening, LASSO-Cox modeling, train-test validation, GEO validation, pathway enrichment, immune signature scoring, mutation/CNV/RPPA support, drug sensitivity prediction, single-cell/spatial localization, and ICB validation were performed using reproducible Python and R scripts. A reduced model was derived using only genes shared by TCGA, GSE13507, and GSE31684. The fixed formula was then applied without refitting to IMvigor210 and GSE176307. Results: A five-gene model composed of EMP1, AHNAK, TNFRSF14, CLEC2D, and GSDMB retained TCGA internal prognostic value (train C-index 0.693, test C-index 0.605, all-sample C-index 0.667; TCGA test log-rank p = 0.015), although GEO survival validation in GSE13507 and GSE31684 was modest. High-risk tumors were enriched for epithelial-mesenchymal transition (EMT), TNF-alpha/NF-kB signaling, inflammatory response, hypoxia, complement, CAF, macrophage, checkpoint, and cytotoxic programs. Single-cell and spatial analyses localized the score to basal tumor, endothelial, fibroblast, and perivascular compartments. In IMvigor210, risk scores were higher in ICB non-responders than responders (Wilcoxon p = 0.044; AUC for non-response = 0.580), high-risk tumors had a lower responder rate (17.6% vs. 28.0%), and high risk predicted poorer OS (log-rank p = 0.016; multivariate continuous risk HR = 3.15, p = 0.044). GSE176307 showed directionally consistent but non-significant response results (AUC = 0.576). Conclusions: The five-gene score is best interpreted not as a standalone universal prognostic classifier, but as a compact stromal-EMT and immune-suppression phenotype associated with inferior ICB response. These findings support a framework linking prognosis, microenvironment biology, immunotherapy resistance, and therapeutic hypotheses in bladder cancer.

Mate preference learning, where individuals use past social experience to choose mates, is prevalent in many species. Yet its consequences on reproductive investment, and whether individuals differentially invest in reproduction based on their learned perception of the mates phenotype, are unknown. We addressed these questions using the butterfly Bicyclus anynana, where males acquire individual preferences for artificially painted 0-UV dorsal hindwing spotted (DHSN) females, who then lay more eggs. We measured the spermatophore proteins transferred by naive and experienced males to females manipulated to have preferred (0-DHSN) and unpreferred (2-DHSN) wing patterns. Using data independent acquisition (DIA), we identified 2144 proteins in the B. anynana spermatophore transferred to the female during mating. Experienced males that mated with preferred and unpreferred females had more differentially abundant (DA) proteins in their spermatophores that have functions in circadian rhythms, oogenesis, and neural signalling than naive males. Proteins associated with oogenesis were DA in spermatophores transferred to preferred females, which may contribute to why preferred females lay more eggs. Overall, our study provides evidence for the role of experience-induced behavioural plasticity in tailoring male ejaculates in butterflies and identifies proteins that influence female physiology and behaviours which directly affect the pairs overall reproductive fitness.

Modeling the circadian impact of seasonal clock changing requires precise synchronization between solar and social time. This report critiques a recent study that associated disease prevalence in the United States with seasonal clock exposure. We identify a fundamental computational error in which a sign reversal of the longitudinal offset effectively inverted the US East-West axis, cross-correlating local health data with the circadian burden of hypothetical locations on the opposite side of a time zone. We outline the methodology for a correct modelization of the circadian process in the context of US geography.

Motivation: With the scaling frequency and intensity of extreme heat events across the globe, it is critical for public institutions to develop early detection systems and continuous monitoring of these events and their impacts. In Brazil, Rio de Janeiro was the first city to publish its heat protocol, with the Rio Heat Dashboard as a central component of this system. Implementation: The dashboard was implemented using R/Shiny and integrates climatic and health data from multiple sources. General features: The application comprises real-time heat exposure monitoring and automatic alert level classification, which is monitored daily by multiple municipal actors and supports activation of actions specified in the heat protocol. It also features a health impact module, which lists each heat event and its impact on mortality, and primary care and emergency visits. Availability: The source for full reproducibility is available through https://github.com/joaohmorais/RioHeatDashboard.

Background: Biochemical recurrence (BCR) occurs in 20-40% of men after radical prostatectomy. Existing postoperative recurrence risk tools based on PSA and pathology are clinically useful but show only moderate and variable discrimination, highlighting the need for biomarkers that improve risk stratification and consequent treatment decisions. We hypothesized that the prostate microenvironment, including both the tumor and non-cancerous adjacent tissue, may contain prognostic features associated with adverse postoperative PSA outcomes. Methods: We assembled a cohort of matched tumor-adjacent benign and tumor prostate tissue from 243 men across three institutions to establish a discovery cohort (n=123; 43 postoperative PSA events, 35%) and validation cohort (n=120; 46 events, 38%). For primary binary analyses, a postoperative PSA event included BCR, defined as two consecutive postoperative PSA values >=0.2 ng/mL, or PSA persistence. We performed RNA sequencing of matched tumor-adjacent benign and tumor tissues, quantified immune signatures, and developed an integrated model combining the adjacent-tissue B-cell signature, preoperative PSA, and radical prostatectomy Gleason score (BRIGADE). CAPRA-S-adjusted Cox analyses excluding recurrence-time-0 cases evaluated time to BCR, and CD19 multiplex immunofluorescence provided tissue-level confirmation (n=10). Results: In prostatectomy specimens, tumors from patients without a postoperative PSA event were enriched for B-cell transcriptional programs, whereas tumors from event-positive patients showed elevated proliferation signatures. B-cell-related transcriptional programs were correlated between tumor and adjacent tissue. Tumor-adjacent benign B-cell scores were higher in no-event cases and discriminated postoperative PSA-event status in PCBN discovery (AUC 0.63) and BM validation (AUC 0.81) cohorts, outperforming numerous other immune-related signatures. In CAPRA-S-adjusted Cox sensitivity analyses excluding recurrence-time-0 cases, higher adjacent-tissue B-cell activity was associated with reduced recurrence risk in PCBN (HR 0.42, 95% CI 0.19-0.94; BH-adjusted p=0.035) and BM (HR 0.54, 95% CI 0.30-0.95; BH-adjusted p=0.034). Tissue-based validation showed that CD19+ B-cell density in adjacent benign tissue was higher in no-event than event-positive patients (median 0.1145 vs 0.0471; p=0.008). BRIGADE achieved an AUC of 0.68 in cross-validation and 0.83 in independent validation, compared to AUCs of 0.54-0.63 and 0.44-0.78 for the tested clinical predictors, respectively. At the fixed classification threshold, the validation-cohort odds ratio for BRIGADE was 2.75. The adjacent B-cell score remained associated with lower odds of a postoperative PSA event after adjustment for PSA and Gleason score. Conclusions: B-cell infiltration in tumor-adjacent benign prostate tissue may complement existing clinicopathologic models for stratifying adverse postoperative PSA outcomes and subsequent BCR after radical prostatectomy. The transcriptomic signal was recapitulated by CD19-based tissue staining, supporting further development of a pathology-based assay.

BackgroundIgA nephropathy (IgAN) is increasingly recognized as a systemic immune-mediated disease characterized by aberrant IgA1 glycosylation, circulating immune complex formation, complement activation, and emerging metabolic perturbations. However, clinical diagnosis still relies on invasive renal biopsy, and non-invasive biomarkers capable of capturing both systemic immune activation and kidney-specific alterations remain lacking. Extracellular vesicles (EVs), as biologically active carriers of proteins and metabolites, provide a unique opportunity to interrogate compartment-specific molecular signatures underlying IgAN pathophysiology. MethodsWe performed an integrated, untargeted multi-omics analysis of plasma- and urine-derived EVs from 60 individuals (24 IgAN, 21 chronic kidney disease [CKD], and 15 controls). Differentially expressed proteins (DEPs) and metabolite features (DEFs) discriminating IgAN from CKD and controls were identified using Venn diagram analysis, followed by pathway enrichment and receiver operating characteristic (ROC) evaluation. ResultsVenn analysis identified 22 and 3 candidate DEPs in plasma EVs (pEVs) and urinary EVs (uEVs), respectively, revealing broader systemic proteomic alterations relative to renal EV cargo. Notably, complement and coagulation regulators, including C4b-binding protein alpha chain (C4BPA) and vitamin K-dependent protein S (PROS1), demonstrated strong discriminatory performance between IgAN and CKD (AUC = 0.826 and 0.795), suggesting EV-associated complement-coagulation crosstalk in IgAN. Metabolomic profiling revealed 1,006 and 540 candidate DEFs in pEVs and uEVs, respectively. Enrichment analyses highlighted steroid biosynthesis and fatty acid metabolism pathways in both compartments, indicating immune-metabolic reprogramming. Three metabolite features (C27H44O, C30H50O, and C28H46O) distinguished IgAN from CKD with high accuracy (AUC = 0.942-0.877). ConclusionsThis study provides the first compartment-resolved, plasma- and urine-derived EV multi-omics landscape of IgAN. Our findings suggest that EV cargo reflects coordinated complement dysregulation and metabolic alterations, extending current understanding of IgAN beyond glomerular immune complex deposition. These EV-associated proteins and metabolites offer a mechanistically informed framework for non-invasive biomarker development and for exploring immune-metabolic pathways involved in IgAN progression.

Brook trout (Salvelinus fontinalis) exhibit strain-level variation in growth performance, environmental tolerance, and survival, yet the biochemical mechanisms underlying these differences remain poorly understood. We developed and applied a ratiometric biochemical framework integrating the pentose-phosphate pathway (PPP) and glutathione metabolism to characterize strain-specific hepatic metabolic organization in brook trout. Five strains reared under standardized conditions differed significantly in hepatic soluble protein density, glutathione pool size, total NADP(H) concentration, and activities of glucose-6-phosphate dehydrogenase (G6PDH), glutathione reductase (GR), and transketolase (TKT). These differences were not uniformly coordinated across pathways, demonstrating that metabolic phenotype cannot be inferred from individual biomarkers alone. Derived ratios describing oxidative-to-non-oxidative PPP capacity (G6PDH/TKT) and glutathione buffering relative to recycling capacity ((GSH+GSSG)/GR) resolved distinct patterns of metabolic allocation among strains. Despite shared ancestry, the Temiscamie (TEM) strain and its domestic x TEM hybrid (TXD) exhibited markedly divergent metabolic phenotypes, demonstrating that closely related strains can differ substantially in hepatic metabolic organization. Together, these findings identify relative allocation among interconnected metabolic pathways as an axis of physiologic diversity and establish a ratiometric approach for comparing metabolic organization across populations and species. Graphical abstractHepatic metabolic phenotypes of brook trout strains were characterized by integrating pentose phosphate pathway enzyme capacities, glutathione metabolism, NADP(H) availability, and soluble protein into a ratiometric framework. Ratios distinguish investment in oxidative versus non-oxidative PPP capacity (G6PDH/TKT), antioxidant buffering versus glutathione recycling capacity (total glutathione/GR), and hepatic protein density (soluble protein/liver mass), revealing distinct metabolic organization among strains. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=88 SRC="FIGDIR/small/743818v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@1694676org.highwire.dtl.DTLVardef@90f2d4org.highwire.dtl.DTLVardef@365327org.highwire.dtl.DTLVardef@8d56ca_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LIA ratiometric framework was developed to characterize hepatic metabolic organization in brook trout C_LIO_LIGlutathione buffering and recycling capacity distinguish alternative redox phenotypes C_LIO_LIInvestment in oxidative and non-oxidative PPP capacity varies independently among strains C_LIO_LIG6PDH/TKT and total glutathione (GSH+GSSG)/GR reveal distinct metabolic phenotypes C_LIO_LIRatiometric indices provide a framework for interpreting redox metabolism and carbon allocation C_LI

The relationship between genotypic and phenotypic variation is determined by the complex interaction of genetic and environmental factors. While statistical methods capable of detecting such interactions exist, an axiomatic mathematical framework that seamlessly describes the combined effects of genetic modifications and environmental exposures on a common scale is lacking. In this report, buffering concepts are used to construct a measurement system that enables the geometric representation of both gene-by-gene and gene-by-environment interactions on the extended complex plane (i.e., as projections on the Riemann sphere). In this manner, any such interaction, or combination thereof, can be precisely defined and quantified as the deviation from the neutral value calculated through the applicable complex transformation. When thus conceptualized, the framework's parameterization defines the "state space" of a given measurable phenotype along both the real and imaginary dimensions, thus establishing an unambiguous and broadly applicable method for determining the phenotypic value expected upon combinatorial changes in genetic and/or environmental variables. Remarkably, by applying these methods, it is possible to quantify the effects of any gene-by-environment interaction using the equation, AGxE=Im([z]obs*zexp)/2, where zobs and zexp are complex numbers representing the observed and expected phenotypes of a given genotype expressed in terms of the buffering parameters, and b.

BackgroundSepsis-associated acute kidney injury (S-AKI) is a major contributor to morbidity and mortality in critically ill patients. However, the molecular mechanisms underlying its temporal progression remain poorly understood because conventional biomarkers primarily reflect renal dysfunction rather than disease pathogenesis. Urinary extracellular vesicles (uEVs), which carry kidney-derived molecular cargo, provide a promising platform for monitoring renal-specific biological alterations during disease progression. MethodsWe conducted a longitudinal multi-omics study of uEVs collected from 81 patients with sepsis, including 48 patients with S-AKI and 33 sepsis-only controls. Patients were randomly assigned to a discovery cohort (n = 52) and an independent validation cohort (n = 29). Urine samples were collected at Day 1, Day 4, and Day 8 after AKI diagnosis. High-resolution proteomic and metabolomic profiling was performed to characterize temporal molecular alterations. Enriched pathways identified in the discovery cohort were evaluated in the validation cohort using pathway-level concordance analysis. ResultsComparative analysis between S-AKI and sepsis-only patients identified distinct stage-specific molecular alterations throughout disease progression. At the early stage (Day 1), validated pathways included complement and coagulation cascades, ferroptosis, HIF-1 signaling, sphingolipid metabolism, and arachidonic acid metabolism, highlighting coordinated inflammatory, hypoxic, and lipid metabolic responses. During the mid-stage (Day 4), persistent activation of complement and coagulation cascades, ferroptosis, and HIF-1 signaling was accompanied by metabolic reprogramming involving alanine, aspartate and glutamate metabolism and tyrosine metabolism. Although limited sample availability reduced statistical power at Day 8, phenylalanine metabolism remained validated in the metabolomic analysis, suggesting persistent metabolic dysregulation during late-stage disease progression. ConclusionsThis study provides the first longitudinal, independently validated multi-omics characterization of human uEVs in S-AKI. By integrating proteomic and metabolomic profiling, we reveal the temporal evolution of renal-specific molecular pathways from early inflammatory and hypoxic responses to subsequent metabolic reprogramming. These findings establish uEV-based multi-omics as a promising strategy for molecular phenotyping of S-AKI beyond conventional clinical biomarkers and provide a valuable resource for future biomarker discovery and therapeutic target identification.

Animal venoms are a rich source of bioactive molecules, yet their diversity remains incompletely characterized in many species. Here we present the first long-read transcriptomic analysis of venom glands from Philippine tarantulas (Theraphosidae), a highly endemic but understudied group. Using Oxford Nanopore sequencing, we reconstructed near full-length venom gland transcriptomes across multiple species and identified extensive repertoires of toxin-encoding peptides. Venom glands were enriched in cysteine-rich inhibitor cystine knot (ICK) peptides, which dominated the toxin landscape and are known modulators of ion channels. Cross-species comparative analyses revealed a distinct transcriptional signature separating venom from non-venom tissues, driven by coordinated expression of toxin-associated and regulatory gene families. Phylogenomic reconstruction based on orthologous peptides recovered expected taxonomic relationships while revealing potential lineage-specific diversification and potential cryptic taxa. Despite a conserved core set of toxin families, substantial variation in toxin composition was observed among species, consistent with rapid evolution driven by gene duplication and functional divergence. Analysis of highly expressed ICK peptides showed a conserved cysteine framework alongside marked sequence variability in inter-cysteine regions, supporting a model in which structural stability is maintained while functional diversification proceeds. Together, these findings establish the first long-read transcriptomic resource for Philippine theraphosid spiders, reveal a conserved molecular signature underlying venom gland specialization, and provide new insights into the diversification of ICK toxin repertoires that may facilitate future evolutionary and functional studies, including the discovery and characterization of bioactive venom peptides.

Summary: Agentic AI platforms enable the engineering of autonomous workflows but are not designed for experimentation and hypothesis testing. ASAREE (Analytical Sandbox for Agentic AI Research, Engineering, and Experimentation), is an open-source platform to address this gap. ASAREE creates agents, connects to MCP servers and tools, and designs factorial experiments through a visual interface or Python SDK. It records a full provenance trace for every run and routes all model calls through a provider-agnostic bridge that supports local deployments, ensuring data privacy. As a use-case, we use ASAREE to evaluate key design choices in a mutli-agent machine learning pipeline. Across a 2 x 2 x 2 factorial design, more advanced models, greater reasoning effort, and critic agent use significantly increased compute time, token use, cost, and feature count without improving predictive performance. The lowest-cost baseline, Claude Sonnet 5 with medium effort and no critic, achieved the highest mean PR AUC while Claude Opus 5 with extra high effort and a critic agent cost 15.5x more (USD) and ran 13.1x longer while performing worse on average. These findings highlight ASAREE as a robust framework for evaluating agentic system performance and resource efficiency.

Transcriptomic profiling has become a cornerstone of modern biomedical research. To make transcriptomic analyses accessible to a broader scientific community, specifically including researchers with limited bioinformatics expertise, we introduced ArrayAnalysis in 2013 as a user-friendly web-based application for microarray data analysis. We now present a major update (https://arrayanalysis.org), introducing a strongly interactive platform that facilitates the dedicated exploration and analysis of both microarray and RNA-seq data, and allows for the generation of publication-ready outputs. Users can perform key analysis steps, including data pre-processing and quality control, differential expression analysis, and gene set analysis, via a sequential, interactive workflow. At each step, the application provides interactive visualizations accompanied by information pages to support interpretation. Users can dynamically adjust figure layouts and colour palettes and export figures as vector graphics and high-resolution raster images. For non-expert users, ArrayAnalysis offers step-by-step guidance to support correct usage and facilitate learning, while for experienced bioinformaticians, it provides a streamlined and flexible workflow ideal for large-scale analyses requiring efficient and consistent processing. ArrayAnalysis is available both as a web application and for local deployment as a desktop application, Docker image, or R package, making it suitable for diverse computational environments, user groups, and analytical purposes. Together, ArrayAnalysis empowers a broad community of biomedical researchers to unlock the full potential of transcriptomic data. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=96 SRC="FIGDIR/small/738193v1_ufig1.gif" ALT="Figure 1"> View larger version (39K): org.highwire.dtl.DTLVardef@1072123org.highwire.dtl.DTLVardef@11095e0org.highwire.dtl.DTLVardef@1dfaee7org.highwire.dtl.DTLVardef@53d31e_HPS_FORMAT_FIGEXP M_FIG C_FIG

Prolonged low-temperature exposure may extend the preservation window of mammalian cells but can also disrupt cellular homeostasis and ultimately compromise cell viability. This study investigated the time-dependent phenotypic and transcriptomic responses of primary canine dermal fibroblasts to sustained hypothermic stress. Passage-three fibroblasts were continuously maintained at 15 for up to 15 days, with samples collected on Days 0, 3, 6, 9, 12, and 15. Cellular morphology, metabolic activity and viability, and apoptosis were evaluated using bright-field microscopy, Cell Counting Kit-8 assays, and Annexin V-FITC/propidium iodide flow cytometry, respectively. RNA sequencing was performed to characterize dynamic transcriptional changes throughout the exposure period. Early low-temperature exposure was associated with relatively preserved cellular morphology and viability, suggesting a transient adaptive response. With increasing exposure duration, fibroblasts exhibited progressive morphological deterioration, reduced metabolic activity, loss of adhesion, and increased apoptosis. Time-series transcriptomic analysis further revealed temporally coordinated and stage-dependent gene-expression programs associated with metabolic regulation, cellular stress responses, structural homeostasis, and cell survival. Integration of phenotypic and transcriptomic data demonstrated that the response of primary canine dermal fibroblasts to 15 was dynamic rather than linear, progressing from early adaptation to cumulative dysfunction during prolonged exposure. These findings provide a framework for defining the low-temperature tolerance of primary canine dermal fibroblasts and may inform the optimization of protocols for their short- to medium-term preservation and transportation.

Survival analysis in head and neck squamous cell carcinoma (HNSCC) is traditionally performed using Cox proportional hazards models, alongside some exploration into black-box machine learning methods. The Super Learner (SL) algorithm addresses this model selection dilemma by combining diverse candidate algorithms into a weighted ensemble to perform comparably to the best candidate method. This study evaluates the performance of SL in HNSCC. Proteomic features as well as clinical covariates from 96 CPTAC HNSCC samples were modeled with three candidate algorithms (Cox LASSO, Cox Ridge, and Random Survival Forest) as well as the ensemble SL method. Models were optimized via Unos time-dependent Concordance Index (C-index) and tested at 1- and 3-year time horizons using 2000 bootstrap resamples. The Cox Ridge regression model achieved the highest predictive accuracy among the four total methods. However, the SL demonstrated stable performance over both time horizons (1-year C-index: 0.985; 3-year C-index: 0.960). Variable importance analysis of the Cox Ridge model successfully identified malignant proteins (ATR, MAML1, MIEN1) alongside novel potential prognostic indicators (ZNF800, KERA). This analysis emphasizes the statistical necessity for larger cohorts for ensemble learning, while providing a benchmark of proteomic indicators in HNSCC.

Background The placenta serves a critical role in nutrient uptake and waste elimination for the developing fetus. The placenta is also responsible for the uptake and/or exchange of xenobiotics, including medications, between the maternal and fetal bloodstreams. An estimated 40-80% of women take medications or drugs during pregnancy for a variety of conditions. Very little is understood about fetal drug and nutrient exposure during pregnancy and how it may change over the course of fetal development. Objective This study aimed to characterize the abundance of transport proteins in placental tissue, which are important in modulating fetal nutrient and drug exposure, over the duration of pregnancy. Mass spectrometry-based global proteomic analysis revealed trends in the expression of thousands of proteins throughout gestation. Focusing on the membrane-associated proteome enabled an increased emphasis on the solute carrier and ATP-binding cassette families of transporter proteins that are critical for nutrient and xenobiotic transport across the maternal-fetal barrier. Study Design Using data-independent acquisition proteomics, relative abundance of proteins in placental tissue samples was profiled across all three trimesters of pregnancy (Trimester 1 = 16, Trimester 2 = 9, and Term = 9). Membrane fractions were generated to enrich membrane-associated proteins for proteomic analysis. Placental samples were grouped into randomized batches for membrane fraction generation and mass spectrometry analysis. Proteomic search results from each batch were imported into the R programming environment from Skyline, concatenated, and normalized as one data set for downstream analysis. Results A total of 6,331 proteins were detected across all samples with 4,210 proteins identified in every sample. Pathway analysis revealed that as gestational age increases, membrane-associated proteins involved in more complex metabolic pathways increase in relative abundance while those involved in extracellular remodeling events and simple organic ion transport tended to decrease. A total of 139 solute carrier and ATP-binding cassette transport proteins were identified in all samples, and 80 were identified in every sample. In general, membrane-associated proteins, including solute carrier and ATP-binding cassette transport proteins, were significantly enriched in placental tissue collected during early gestation compared to term placental tissue. Conclusion This study presents a comprehensive profiling of membrane-associated proteomic changes during gestation and identifies significant gestational age associated abundance changes at the protein level in several transport protein families. The application of data-independent acquisition global proteomic techniques enabled in-depth analysis of thousands of proteomic changes across pregnancy in a single experiment. These data provide critical information to support future studies into the understanding of fetal exposure to xenobiotics and nutrients circulating in the maternal bloodstream.

Dung beetles can use a variety of orientation cues to maintain a consistent bearing during ball-rolling. Where several cues are available, they appear to learn the spatial relationship between them, providing redundancy if some cues are removed. Mounting evidence indicates that such a learning process is implemented in the insect head direction circuit; specifically, in the plastic substrate between sensory input neurons and compass neurons in the central complex. This plasticity appears to be driven by rotational movements, providing a clear link with observed beetle 'dance' behaviour. Here, we extend our functional model of this circuit and use it on a robot platform, to test it in the same behavioural assay as was used for the beetles. The robot was able to replicate the beetle's ability to substitute a directional wind cue for a point source light cue in guiding straight-line movement. However, it also revealed significant biasing coupled to dance direction. This biasing appears to be caused by inherent conflict between recurrent and instantaneous inputs to the compass circuit. We predict that the real insect should experience similar issues unless it has evolved a neural mechanism to compensate.