LATEST IN SCIENCE

Abstract: High-throughput sequencing has generated protein datasets whose scale increasingly exceeds the practical limits of conventional functional annotation workflows. We present Sma3s v3, a scalable reimplementation of the Sma3s three-step annotation strategy, which combines transfer from highly similar homologs, orthology-based inference, and functional enrichment among homologous proteins. Sma3s v3 replaces BLAST-based searches with MMseqs2 and introduces parallel processing, reusable SQLite caches, taxonomic filtering, and traceable outputs that retain the evidence underlying each assignment. We evaluated the method on a Vibrio cholerae pangenome comprising 50,415 gene clusters from 11,295 quality-filtered genomes and on a metagenomic catalogue containing 843,935 proteins. After excluding non-informative assignments, Sma3s v3 annotated 30,662 pangenome clusters (60.8%), comparable to InterProScan (60.2%) and exceeding eggNOG-mapper (41.4%), while providing 5,747 annotations not recovered by either comparator. Gene Ontology comparisons showed broad semantic agreement between methods, with Sma3s v3 frequently contributing more non-redundant information in Molecular Function and Biological Process. Within the pangenome, annotation coverage reached 97.1% for core clusters and approximately 59% for accessory and unique clusters. Exact protein matches to non-Vibrio genera identified 1,838 candidate horizontally transferred clusters enriched in genetic mobility, antimicrobial resistance, and metal tolerance functions. In the metagenomic catalogue, Sma3s v3 annotated 728,014 proteins (86.3%), compared with 616,895 (73.1%) using InterProScan 2026, and recovered approximately 20,000 unique functional terms. These results establish Sma3s v3 as a scalable and interpretable tool for functional annotation and re-annotation of proteomes, pangenomes, and metagenomic protein catalogues.

Abstract: Separating somatic from germline variants and reconstructing somatic haplotypes are the two central problems of tumor-only cancer genome analysis. Long reads carry the linkage needed to solve both, but chromosome-scale loss of heterozygosity (LOH) and an unknown degree of normal-cell admixture blur the distinction between somatic and germline haplotypes. Here we present LongPhase-TO, the first method to reconstruct somatic haplotypes from a tumor sample alone. Rather than mapping somatic variants onto germline haplotypes, LongPhase-TO co-phases germline and somatic alleles in a unified graph, in which LOH and tumor DNA fraction are resolved internally from heterozygosity depletion and haplotype imbalance rather than a copy-number and ploidy model. Across eight datasets from six cancer cell lines, LongPhase-TO increased haplotype block N50 by a median of 2.9-fold relative to germline phasers. It also consistently improved somatic single-nucleotide variant (SNV) and indel calls from ClairS-TO and DeepSomatic-TO, raising mean F1 from 0.55 to 0.62 and 0.65 for SNVs and from 0.19 to 0.23 for indels, with the largest gains at low tumor DNA fraction. Across breast, melanoma and lung cancer cell lines, LongPhase-TO improves the accuracy of existing somatic callers and reconstructs megabase-scale somatic haplotypes.

Abstract: Identification of transcription factor target gene interactions and construction of the gene regulatory networks (GRNs) are essential for understanding the molecular mechanisms underlying transcriptional gene regulation. Large scale single cell transcriptomics across different tissues offers unprecedented resolution of cellular diversity and regulatory dynamics by capturing gene expression heterogeneity. However, existing methods often lack effective multimodal integration and fail to fully exploit the hierarchical structure in Gene Ontology (GO) and gene sequence level representations, which limits their ability for predictive performance and biological interpretability. We present scMGFGRN, a multi-model deep learning framework that integrates single-cell transcriptomic profiles with GO hierarchical relationships, gene sequences by leveraging denoising auto encoders, graph attention feature extraction and pertained DNA language model to capture multi-source dependencies within multi-model biological knowledge, while its gated multi head attention module effectively identifies informative regulatory signatures and integrate complementary features from different sources to predict accurate gene regulatory networks. Benchmarking on the seven datasets of human and mouse demonstrates that scMGFGRN outperforms state of the art methods in identifying GRNs. Further analyses reveal that scMGFGRN effectively identifies novel TF gene interactions (TGIs) and reconstructs cell type specific GRNs. Interpretability analysis reveals the contribution patterns of heterogeneous biological sources, demonstrating the ability of scMGFGRN to integrate transcriptomic profiles with multi model structure information.

Abstract: Immune checkpoint inhibitors (ICI) have revolutionized the treatment of brain metastasis (BM), demonstrating intracranial response rates and durability not previously seen with other therapies. However, approximately 50% of patients do not respond to ICI, and the brain-specific interactions that shape anti-tumor immunity remain poorly understood. We profiled five syngeneic BM models using spatial transcriptomics, identifying myeloid-rich niches enriched for interferon-responsive and disease-associated microglia (DAM)-like programs. Ligand-receptor colocalization analyses identified progranulin (PGRN) as a candidate mediator of these niches. Host- or tumor-cell Grn loss reduced BM burden, and Grn-deficient macrophages exhibited altered metabolic programs and tumor-cell engulfment in vitro. Across independent human BM datasets, GRN expression was associated with conserved lysosomal and DAM-like myeloid programs, and GRN-high myeloid regions colocalized with immunosuppressive signatures and dysfunctional CD8+ T cell states. These cross-species findings identify PGRN as a candidate mediator of the BM immune niche and support further investigation of its therapeutic relevance.