iMeta

Microbiome analysis has become a pivotal tool in understanding the role of microbiota in human health and disease. However, the lack of standardized workflows, together with the limitations of proprietary software solutions, hampers reproducibility and flexibility. Here, we present Mbiome, an open-source, user-friendly and automated workflow designed to streamline amplicon-based microbiome analysis. Built upon QIIME2, Mbiome supports both bacterial (16S rRNA) and fungal (ITS) profiling, and is compatible with raw fastq files generated by Ion Torrent (IT) and Illumina (IL) sequencing platforms. The workflow guides users through an interactive setup process via a simple configuration file, enabling researchers with minimal bioinformatics experience to perform comprehensive analyses without writing code. Once configured, Mbiome automates major steps including quality control, taxonomic assignment, and {beta}-diversity analyses, functional predictions (via q2-metnet), and customizable visualizations and statistical analyses. Mbiome has been validated using real-world datasets from multiple sclerosis research projects, performing a comparison between different microbiome analysis approaches, including 16S hypervariable region reconstruction, amplicon-based strategies, and cross-platform sequencing (IT and IL), as well as against results obtained with Ion Reporter (IR) commercial software. This evaluation demonstrated its versatility and effectiveness across different sequencing platforms. Moreover, Mbiome provided improved flexibility, transparency, and taxonomic resolution compared to IR. By combining accessibility, reproducibility, and cross-platform compatibility, Mbiome lowers the barrier to microbiome data analysis and facilitates high-quality, standardized workflows in both research and applied settings. Mbiome is publicly available at https://github.com/MGorostidi/mbiome.

Motivation: Horizontal gene transfer (HGT) shapes bacterial evolution and microbial ecosystems, yet detecting HGT within microbiomes remains a challenge due to fragmented metagenomic assemblies, reference bias, reliance on gene boundaries, and limited ability to model structural mosaicism and patterns across genomes. Methods: We present Kente, a novel pangenome graph-based framework designed for HGT detection that aligns metagenomic assembly contigs to a curated database of >600 genus-level bacterial pangenome graphs constructed using minigraph. Kente infers local taxonomic composition along contigs using alignment evidence and classifies candidate transfers using structured clade-transition topologies (e.g., A-B-A sandwich, open tips, and mosaic patterns). A complementary intra-genus module detects inter-species transfers within a single genus graph using segment-level clade annotations. Results: Across simulated intra- and inter-genus transfer scenarios, Kente achieves higher precision and comparable recall relative to existing gene-centric microbiome HGT detection approaches while reducing false positives from fragmented assemblies. Application to real human gut metagenomes (HMP2, n = 26) demonstrates Kente's ability to detect candidate cross-lineage transfer regions in complex microbial communities. Runtime profiling shows near-linear scaling with input size, enabling efficient analysis of large metagenomic assemblies. Availability and Implementation: https://github.com/treangenlab/Kente

Viruses play indispensable roles in ecosystems and human health, yet deciphering their molecular functions remains challenging. Many viral protein annotations are incomplete or poorly characterized. Existing tools typically predict functional categories without linking to verifiable evidence, hindering the credibility of functional interpretation. Here, we present VirProtRAG, a viral protein function annotation framework that integrates information retrieval with evidence-grounded knowledge generation. It introduces three task-adapted components: a hybrid retrieval module combining keyword-based and semantic dense retrieval to maximize literature coverage, synonym-expanded and rank-aware retrieval with reciprocal rank fusion for improved search effectiveness, and literature quality and evidence-oriented re-ranking to enhance reliability and interpretability. Results show that hybrid retrieval strategy performed best, with quality and evidence features further enhancing re-ranking. Compared with direct LLM prompting without retrieved literature, it consistently improves generation performance, underscoring the critical role of external knowledge. Finally, we built a searchable database comprising all 17,484 reviewed Swiss-Prot viral proteins, supporting both sequence- and text-based queries. VirProtRAG introduced 32.53% non-overlapping function annotations beyond existing expert curation, and independently supported 56.34% of sequence-inferred function points with retrieved literature. Case studies further demonstrate its capability to augment and refine the characterization of previously unannotated or poorly understood viral proteins.

Genome-scale metabolic models (GSMMs) are important aids towards system-level understanding of the metabolic physiology of the gut microbes and for rational microbiome engineering. While large-scale repositories of GSMMs for gut-associated bacteria are available, strain-level variability and the continuous discovery of novel taxa through metagenomics and culturomics underscore the need for scalable, ab initio reconstruction tools. Here, we present CarveMe-GutMicrobes, a client-side framework for rapid reconstruction of metabolic models directly from (meta)genomic input. Building upon the original CarveMe framework, CarveMe-GutMicrobes incorporates an expanded, gut-microbe-centric biochemical database that includes reactions, metabolites, and gene-protein-reaction (GPR) associations curated specifically for Bacteria and Archaea inhabiting the human gut. The tool supports taxonomic restriction of the reference database to improve context-specific accuracy. To test the CarveMe-GutMicrobes and to address the paucity of experimental data for non-model gut taxa, we generated new experimental datasets on metabolite secretion profiles and gene essentiality. CarveMe-GutMicrobes models demonstrated high predictive performance performance against these as well as previously available datasets. By integrating curated resources, extending reaction coverage, and offering new empirical datasets, CarveMe-GutMicrobes provides a scalable platform for high-resolution metabolic reconstruction towards broader adoption of GSMMs in gut microbiome research.

Genomic epidemiology workflows typically require expert curation of multiple specialized tools, extensive manual parameter tuning, and access to heterogeneous compute infrastructure. While standard generative AI models often hallucinate in complex biological domains, we introduce biomeStat: an autonomous AI agent that functions as a strict deterministic orchestrator. By automatically writing code to execute established bioinformatics tools in sandboxed environments, biomeStat dynamically provisions compute resources (CPU and GPU) and guarantees reproducibility, making it immediately useful for scientists without requiring command-line expertise. To demonstrate the platform, we performed a fully autonomous genomic epidemiology and structural analysis of 1,000 Dengue virus (DENV) genomes sampled from 16 Asian countries between 2000 and 2025. The agent seamlessly orchestrated phylogenetic reconstruction (IQ-TREE, TreeTime), Bayesian phylodynamics (BEAST2 via NVIDIA H200 GPU), selection pressure analysis (HyPhy), and structural mapping (PyMOL). The analysis was completed in under 24 hours of wall-clock time, revealing endemic stability (R_e [~]1.0) and identifying 1,869 candidate immune escape sites structurally colocalized with B-cell and T-cell epitopes. Furthermore, the agent validated 176 highly conserved drug target residues across the viral replication complex, confirming that resistance-associated positions for emerging antivirals JNJ-1802 and NITD-688 remain absolutely conserved across all four serotypes. By bridging the gap between natural language intent and deterministic computational execution, biomeStat reduces weeks of expert effort into a single-session analysis with full methodological transparency.

Here we introduce OpenEvo, a fully open-source, low-cost turbidostat platform for automated continuous culture and directed evolution experiments. Existing tools are expensive, complex, or lack open-source hardware; OpenEvo addresses this gap. OpenEvo is a complete, fully automated evolution platform with detailed, illustrated construction instructions for beginners, open-source software and firmware, and a single device priced around $300. An optional PC-based version offers enhanced functionality, including remote access, programmable evolution cycles, programmable LED stimulation, and a data visualization tool. OpenEvo can cycle through three types of media for positive, negative, and neutral selection conditions, supporting a wide range of experimental designs. We validate the use of OpenEvo by evolving H. volcanii to grow from 15% to 12% salt over ~150 cycles, ~1,000 hours. Evolved cells grew 36% faster than wild-type at 12% salt. Whole-genome sequencing of adapted cells found SNPs and large deletions. We also demonstrate positive and negative selection using the OpenEvo LEDs to drive optogenetics via a Phytochrome B-based optogenetic tool, with light as the selection stimulus during over 4000 hours of growth. OpenEvo lowers the technical and cost barriers for continuous evolution experiments, serves as a teaching tool, and is designed to grow an open community of users who share modifications.

The International Committee on Taxonomy of Viruses (ICTV) is responsible for developing and maintaining a universal virus taxonomy. As the reference framework for organising the viral world, it is essential for virology and related fields. Despite its widespread use in research and public health, programmatic access to ICTV taxonomy has remained limited, posing challenges for integration, versioning, and interoperability across databases and bioinformatics resources requiring up-to-date virus taxonomy. To address this, we developed a public and sustainable solution leveraging ontology-based APIs. Successive ICTV Master Species List (MSL) releases were transformed into a structured ontology and deployed as a unified representation through the Ontology Lookup Service (OLS). The framework also provides ICTV-NCBI mappings and helper libraries for integration into downstream systems. This enables, for the first time, public programmatic retrieval of current and historical virological taxon names, taxonomic relationships, metadata, and persistent identifiers through stable endpoints. More broadly, this work illustrates a general strategy for transforming structured biological datasets into semantically enriched graph resources exposed through scalable public APIs. These developments enhance interoperability, reduce manual curation, and support FAIR-aligned taxonomic data management in virology and pandemic preparedness. Key pointsO_LIICTV provides the official taxonomy for classifying viruses and naming virus taxa, but lacks standardised programmatic access. C_LIO_LITransforming ICTV data into an ontology enables semantic, machine-actionable access across releases via ontology-based APIs. C_LIO_LIICTV-NCBI mappings support interoperability across bioinformatics resources. C_LIO_LIThe framework enables programmatic resolution of current and historical viral taxa. C_LIO_LIThis approach provides a reusable model for exposing biological datasets through public APIs. C_LI

The human gut, containing 100 trillion microbes, is also considered the "second brain," having control over the different functions of the physiological system. With advancements in bioinformatics and the development of sequencing technologies, researchers are able to explore the diversity and functional implications of gut microbiota (GM), which have become strongly associated with a variety of diseases. Microbial imbalance, or dysbiosis, acts as a biomarker for early detection and prognosis of a disease. Artificial Intelligence and Machine Learning (AI/ML) methods, although extensively used in predicting GM associated diseases, are seldom translated to having practical real-world outcomes, necessitating the design of robust AI/ML models applicable in real-world scenario. We have therefore come up with designing stacking-based ensemble architectures (EM1 and EM2), developed by integrating multiple ML-based learning algorithms for improving disease prediction accuracy. The GM datasets, after split into training and test sets, were eventually fed into the proposed two-layer ensemble models, which combines the output from standardized base learners via a meta-classifier, strengthening classification robustness as well as ensuring consistency in optimized performance across diverse datasets. Both the proposed hybrid ensemble models have emerged to be superior performers over all baseline and deep learning models, with an average accuracy of 0.87 and 0.84 respectively. By combining multiple learners, the proposed ensemble models outperform traditional single-algorithm-based approaches to attain higher accuracy and robustness on complex GM datasets. Key messagesO_LIDevelopment of stacking-based hybrid ensemble models (EM), which can be employed to integrate different AI/ML algorithms with better prediction accuracy of gut microbiome (GM)-associated diseases. C_LIO_LIUse of independent GM datasets with preprocessing methods such as SMOTE and PCA to address class imbalance and high dimensionality. C_LIO_LIAll the proposed EM architectures are mostly superior to the existing state-of-the-art AI/ML methods (highest prediction accuracy: 0.87 and 0.84 with EM1 and EM2 models respectively) for GM diseases predictions. C_LIO_LIThe cross-cohort validation demonstrates high prediction accuracy and robustness, (AUC values close to 0.98 and 0.99, for EM1 and EM2). C_LIO_LIThese therefore demonstrate the effectiveness of EM frameworks for GM associated disease prediction, paving the way for corresponding applications in precision medicine. C_LI

Accurate coding sequence prediction in short prokaryotic metagenomic reads remains challenging due to sequence fragmentation, unknown sequence origins, and sequencing errors. Here we introduce DeepCDS, a deep learning-based ab initio coding sequence predictor trained on short prokaryotic sequences with and without simulated Illumina-like sequencing errors. DeepCDS integrates ESM-2 protein language model embeddings with nucleotide-level information to predict complete and fragmented coding sequence regions. Benchmarking on 215 phylogenetically diverse prokaryotic organisms demonstrates that DeepCDS consistently outperforms current state-of-the-art methods in coding sequence detection, start and stop codon localization, and robustness to different sequencing error profiles, while remaining operational at shorter sequence lengths than existing tools support. These findings demonstrate that protein language models capture distinct signals relevant for nucleotide-level coding sequence detection, especially at very short lengths. Ultimately, DeepCDS may help uncover the functional potential of the vast microbial diversity that remains genomically uncharacterized.

Organisms have continuously evolved in response to environmental conditions. Pathogenic bacteria evolve under host and environmental pressures, reshaping their genomes through insertions, inversions, deletions, and duplications during infection. In clinical settings, phenotypic traits of pathogenic bacteria such as virulence or antimicrobial resistance directly affect disease severity, transmission, and treatment. Conventional genotyping provides insights into genomic relatedness but does not always align with these clinically relevant traits, limiting its utility for phenotype-driven interventions. Here, we develop ABComp (Assembly polishing and Bacterial whole-genome Comparison for multi-group clinical isolates), a modular and Snakemake-based workflow for phenotype-driven comparative genomics. ABComp automates assembly polishing, group-wise pangenome analysis, and enables flexible pathogenic marker discovery through user-defined comparisons. We validated ABComp using a Klebsiella pneumoniae ground truth dataset stratified by yersiniabactin presence and successfully recovered the entire locus as a group-specific core marker. By applying ABComp to another dataset of clinical isolates with experimentally measured virulence, we discovered the ferric citrate (Fec) uptake system as a potential marker specific to a hypervirulent group. These results demonstrate ABComps utility in uncovering phenotype-linked genomic markers with clinical significance, supporting targeted treatments and rapid diagnosis.

The rapid rate of virus discovery renders manual curation by taxonomy experts increasingly impractical, creating a need for reliable software that can reproducibly assign viral contigs to taxa at all fifteen ranks of the virus taxonomy. We led an open community challenge for the computational taxonomic classification of viruses and assembled a dataset of virus sequences combining expert-curated and metagenomic sequences. Seventeen teams contributed a total of thirty-four automated, fully reproducible classification pipelines. Most tools correctly assigned viruses belonging to established species, genera, or families, but viruses that are unclassified at those lower ranks remain challenging. This study provides datasets, open-source software, novel approaches, and recommendations to benchmark computational taxonomic classification of viruses, and support organizing the many viruses discovered in big omics data.

Flavins are organic redox cofactors central to metabolism and uniquely capable of acting as extracellular electron shuttles. For life to have emerged, it must have disengaged itself from its stationary geochemical environment, a step requiring mobile redox-active components. The role of flavins at life's origin has been debated for decades, centered on their capacity for both one- and two-electron chemistry, distinguishing them from nicotinamides and iron-sulfur clusters. Here we chart the abiotic reduction of flavin mononucleotide (FMN), flavin adenine dinucleotide (FAD), and riboflavin under hydrothermal conditions (40 {degrees}C, 1 bar N2 or 5 bar H2, pH 6, 8, and 10) by nickel (Ni) and iron (Fe). Flavins show greater environmental versatility than hydride carriers such as NAD and can harvest electrons from metals that would otherwise reduce water's protons to H2. Reduction is favoured under acidic conditions, while increasing molecular charge at higher pH impedes electron transfer. Ni acts as a hydrogenation catalyst, reducing deprotonated flavins via hydride transfer, suggesting mineral composition could have influenced geochemical selection of early electron carriers. Reduced FMNH2 and FADH2 were tested as electron shuttles toward Fe3+-containing minerals, revealing that FMNH2 enables faster mineral dissolution than FADH2. We further demonstrate complete redox cycling of FMN through Ni-assisted H2 reduction and subsequent oxidation by magnetite (Fe3O4) under inert atmosphere, releasing Fe2+. This study highlights the versatility, stability and redox chemical capabilities of flavins in prebiotic context.

GenBank and GISAID are the primary repositories for viral genomic data, but integrating records across them remains a challenge. The same sequence could be made available in both databases without any cross-reference linking the two entries. Consequently, there is no systematic way to identify this redundancy, which compromises the compilation of representative, non-redundant large-scale datasets. In parallel, the growth of viral genomic data has increased the risk of systematic technical artifacts introduced during sequencing or assembly. These artifacts can inflate substitution rate estimates and degrade temporal signal, biasing evolutionary rate estimates. To address both challenges, here we present a formal, reproducible workflow integrating two newly developed complementary tools: G2G matcher for cross-repository harmonization and Lab-Specific Bias FILTer (LSBFILT) for masking of laboratory-specific artifacts. Using the Eastern/Central/South African (ECSA) chikungunya virus lineage as a proof-of-concept, we demonstrate that our integrated workflow restores temporal signal and provides a robust, curated dataset for downstream phylodynamic analyses. Critically, restricting masking of homoplastic sites to specific sequences reduces the substitution rate estimate from an inflated 8.517 x 10-4 to 5.078 x 10-4 substitutions/site/year and increases the coefficient of determination (R2) of the root-to-tip regression analysis from 0.353 to 0.677. By enabling systematic cross-repository harmonization and source-specific artifact masking, we provide the molecular epidemiological community with scalable tools to reconcile fragmented genomic data and reduce technical biases, fostering more accurate and reproducible phylogenetic analysis. G2G matcher is available at https://github.com/andrezaleite/G2G-Matcher, and LSBFILT at https://github.com/khourious/LSBFILT.

Motivation: Bacteria carry a large repertoire of antiviral defence systems, our knowledge of which is expanding rapidly. Several bioinformatics tools now exist to identify them. Though powerful, these tools can differ in the models they use and the nomenclature they return, thus a single tool could both miss an annotation and disagree with its peers. Results: Here we describe Ptolemaea, a pipeline for harmonising phage-defence annotations across multiple tools by reconciling PADLOC, DefenseFinder, and a bidirectional BLAST. Over a common predicted set of proteins, Ptolemaea provides a consensus annotation list per genome. The pipeline is not intended to outperform or replace its component tools; its purpose is to maximise the number of defence systems recovered from a genome and to make disagreements between tools explicit and resolvable. We demonstrate the pipeline on 700 complete genomes spanning the ESKAPE pathogens and Escherichia coli, recovering 32,509 defence annotations, of which 50.6% were supported by more than one annotation source. Availability: The Ptolemaea pipeline is freely available at https://github.com/ecampbell50/Ptolemaea. Supplementary information: Genome accessions used in this analysis can be found in S1, and script for genome retrieval in S2. Collated consensus annotation counts can be found in S3, while all raw tool outputs and curated decisions for each species can be found in S4-10

Intestinal organoids are powerful in vitro models, but their use in large-scale analyses remains constrained by the low throughput, labor-intensive handling, and high reagent consumption of conventional Matrigel dome culture. Here, we present Organoid-in-Bead (OrB), a vortex-based compartmentalization workflow that partitions organoid fragments into thousands of discrete Matrigel microbeads, enabling scalable, high-density culture from a single batch preparation. OrB maintains dome-comparable organoid growth and epithelial polarity, supports passaging-based culture expansion, yields more than 5,000 organoids in the final 10 cm dish format, and reduces Matrigel and medium consumption by approximately 70% on a per-organoid basis. OrB therefore provides a practical and scalable upstream workflow for generating screening-scale intestinal organoids. HighlightsO_LIOrB generates Matrigel microcompartments by vortexing without microfluidics C_LIO_LIOrB enables scalable, high-density intestinal organoid culture in one batch C_LIO_LIOrB maintains dome-comparable growth and epithelial polarity and supports passaging C_LIO_LIOrB yields >5,000 organoids per batch with [~]70% less Matrigel/medium per organoid C_LI

Quantifying short-term evolutionary rates of microbial genomes is essential for understanding the processes that shape within-host evolution and for establishing thresholds needed to track transmission. In studies of short-term evolutionary rates, samples are often collected from closely related clusters (e.g. longitudinally from the same host or from transmission pairs), with substantial time intervals separating genomes between clusters. Distinguishing strain replacement from persistence presents is also difficult in these studies. In addition, many public health and metagenomic bacterial strain tracking pipelines output pairwise SNP distances rather than the multiple sequence alignments required by common substitution rate estimation pipelines. This makes it hard to estimate within-host evolutionary rates in many commensal bacterial species that are difficult to culture and isolate. To address these challenges, we introduce MxSure, a tool for estimating substitution rates and transmission thresholds while accounting for strain replacement from pairwise SNP distance data, as commonly generated by transmission tracking and metagenomic analysis pipelines. We demonstrate the accuracy of MxSure through extensive simulations and by analysing species with previously estimated substitution rates from longitudinal metagenomic datasets. Using MxSure, we estimated within-host substitution rates and transmission SNP thresholds for multiple commensal bacterial species including Bifidobacterium longum and Bifidobacterium bifidum from a longitudinal study of the infant gut microbiome.

Integrity of scientific data is critical in biomedical research, where images often serve as primary evidence for experimental observations and conclusions. Advances in image-editing technologies and generative artificial intelligence (AI) have increased the accessibility and realism of visual manipulation, making detection through manual review increasingly challenging. To empower our laboratory researchers to continuously monitor and uphold scientific rigor and data integrity, and serve the global scientific community, we developed VrySure, an easy-to-deploy, AI-driven multi-task platform for automated image-integrity screening in biomedical research. VrySure integrates four detection modules: cross-image transformation detection, within-image copy-move detection, splicing detection in blot and gel images, and AI-generated image detection. The system identifies potentially manipulated images and, when possible, localizes suspicious regions using bounding-box outputs to support downstream verification. To support development and evaluation, we constructed task-specific datasets by combining public biomedical image resources, curated manipulated examples, and synthetic images generated by multiple generative AI systems. We evaluated VrySure using region-level F1 score, recall, precision, false negative rate (FNR), and false discovery rate (FDR) across multiple manipulation categories and compared its performance with two commonly used commercial image-integrity screening platforms under a predefined benchmark protocol. Under the tested conditions, VrySure achieved a higher F1 score and recall, lower FNR, and maintained a low FDR for within-image copy-move detection, splicing detection, and AI-generated image detection, while showing comparable performance in transformation detection. Beyond automated screening, VrySure is designed to support source-data comparison and evidence-based assessment in scientific integrity investigations. By integrating multiple detection capabilities into a unified and scalable workflow, VrySure provides a practical framework to improve the efficiency and consistency of image-integrity screening in biomedical research.

Cysteine and serine proteases act as an immune hub in the plant apoplast to provide robust extracellular immunity during microbial colonisation. Microbial pathogens counteract these immune proteases by inhibiting their activity using small secreted proteins (SSPs). Traditionally, SSPs with protease-inhibitory activity are predicted using sequence-dependent database searches. However, in recent years, fungal SSPs have been shown to exhibit protease-inhibitory functions despite lacking the inhibitor domain that is annotated through sequence similarity searches. Hence, a large number of these novel SSPs with putative protease inhibitor functions are missed during detection and filtered out during sequence similarity searches. This necessitates the development of newer approaches to predict SSPs lacking an annotated inhibitor domain. Machine learning approaches, such as protein language models, have emerged as powerful tools for predicting protein functions. To date, no machine learning models have been developed to predict the protease-inhibitory activities of SSPs lacking an annotated inhibitor domain. Here, we introduce a protease inhibitor prediction pipeline, PINPOINT (Protease INhibitor PredictiOn at plant-pathogen INTerface). The PINPOINT pipeline combines fine-tuned protein language model classifiers, a structure-aware autoencoder, and effector prediction into a multi-level framework for identifying SSPs with predicted protease inhibitor functions. PINPOINT predicts protease inhibitors using SSPs sequences and monomeric structures with pre-computed structures obtained from the AlphaFold Protein Structure Database or predicted using the ESMFold public API. We successfully validated the PINPOINT platform using SSPs from the plant fungal pathogen Macrophomina phaseolina. Notably, the PINPOINT platform robustly predicted several of these SSPs as protease inhibitors including Sequence-unrelated but structurally similar (SUSS) effectors. We further validated the inhibitory potential of these predicted M. phaseolina SSPs using AlphaFold Multimer (AFM) screening against candidate apoplastic soybean cysteine and serine proteases. Additionally, this platform can be used as a pre-filtering step in AFM screening approaches to reduce the number of candidates for discovering novel SSPs with protease inhibitor function for cross-kingdom plant-microbe interaction studies. The PINPOINT platform will accelerate the prediction of novel SSPs including SUSS effectors with protease inhibitor functions in proteomes of any organisms. We made the PINPOINT pipeline accessible to the research community as a web-based notebook environment for interactive computing in Google Colab, available at https://github.com/iitj-mpg-lab/PINPOINT

Brain organoids offer an invaluable model system for studying human brain development and disease. However, the establishment of high-fidelity brain organoids with multiple cell lineages including vasculature and immune cells remains a huge challenge. Here, we present a new strategy to generate human cerebral organoids with vasculature and microglia-like cells using genetic code expansion technology (GCE-T) via site-specific protein engineering. The strategy integrates orthogonal genetic translation machinery in hPSCs via PiggyBac transposon system, enabling temporally control of ETV2 expression and endothelial differentiation in hPSC-derived cerebral organoids. The vascularized human cerebral organoids (vhCOs) exhibit coordinated development of multiple cell lineages and blood-brain barrier (BBB) features. Moreover, vhCOs form perfusable vascular network after transplanted in the immune-deficient mice. Single-nucleus RNA sequencing reveals enhanced neurovascular interactions, multi-brain-regional identities, diverse neuronal subtypes and specialized endothelial subclusters in vhCOs, closely resembling human fetal brain. Strikingly, we identify enriched microglia-like cells comprising three distinct subtypes in vhCOs, which contribute to microglia-vascular interactions and synergistically modulate vascular development. Upon Zika virus (ZIKV) infection, vhCOs show neurovascular dysfunction and impaired microglia development, offering new insights into viral-induced neurodevelopmental disorders. This study offers a unique platform for producing more valuable brain organoids with vasculature and immune components, opening a new avenue to advance organoid research and applications.

Drug-induced liver injury (DILI) remains a major cause of clinical attrition and postmarketing withdrawal, but structure only DILI predictors are difficult to compare because public benchmarks are vulnerable to compound overlap, scaffold similarity and shared label provenance. We present OakuloidTM, an open DILI prediction framework that pairs a leakage audited structure based model with an optional iBAC 3D primary human hepatocyte IC50/Cmax confirmation signal. The structural model integrates gradient boosted descriptor backbones, fingerprint random forests and LivTox proxy-DILI features through a logistic meta learner. Its evaluation is designed as part of the contribution: internal DILIrank, strict external TDC, scaffold disjoint TDC and independent Geci provenance checks are reported with released per compound predictions. Oakuloid reaches AUROC 0.811 on the strict external TDC benchmark and remains competitive under scaffold and fully clean TDC filtering. A channel attribution ablation shows that the external benchmark lead is driven by descriptor based gradient boosted trees rather than by DILIPredictor derived proxy features, reducing a potential circularity concern. The wet lab IC50/Cmax signal is largely orthogonal to structure and supports a confirmation mode that shifts the internal operating point toward higher specificity without claiming a universal AUROC gain. Oakuloid is released with code, model artifacts, calibration analysis, a 122 compound wet lab benchmark and a model card under the Apache License 2.0, supporting reproducible DILI screening and benchmark auditing.