Bento: a toolkit for subcellular analysis of spatial transcriptomics data.

The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell-cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene-gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

Identifying the core genome of the nucleus-forming bacteriophage family and characterization of Erwinia phage RAY.

We recently discovered that some bacteriophages establish a nucleus-like replication compartment (phage nucleus), but the core genes that define nucleus-based phage replication and their phylogenetic distribution were unknown. By studying phages that encode the major phage nucleus protein chimallin, including previously sequenced yet uncharacterized phages, we discovered that chimallin-encoding phages share a set of 72 highly conserved genes encoded within seven distinct gene blocks. Of these, 21 core genes are unique to this group, and all but one of these unique genes encode proteins of unknown function. We propose that phages with this core genome comprise a novel viral family we term Chimalliviridae. Fluorescence microscopy and cryo-electron tomography studies of Erwinia phage vB_EamM_RAY confirm that many of the key steps of nucleus-based replication encoded in the core genome are conserved among diverse chimalliviruses, and reveal that non-core components can confer intriguing variations on this replication mechanism. For instance, unlike previously studied nucleus-forming phages, RAY doesn't degrade the host genome, and its PhuZ homolog appears to form a five-stranded filament with a lumen. This work expands our understanding of phage nucleus and PhuZ spindle diversity and function, providing a roadmap for identifying key mechanisms underlying nucleus-based phage replication.

Identifying the core genome of the nucleus-forming bacteriophage family and characterization of Erwinia phage RAY.

We recently discovered that some bacteriophages establish a nucleus-like replication compartment (phage nucleus), but the core genes that define nucleus-based phage replication and their phylogenetic distribution were still to be determined. Here, we show that phages encoding the major phage nucleus protein chimallin share 72 conserved genes encoded within seven gene blocks. Of these, 21 core genes are unique to nucleus-forming phage, and all but one of these genes encode proteins of unknown function. We propose that these phages comprise a novel viral family we term Chimalliviridae. Fluorescence microscopy and cryoelectron tomography studies of Erwinia phage vB_EamM_RAY confirm that many of the key steps of nucleus-based replication are conserved among diverse chimalliviruses and reveal variations on this replication mechanism. This work expands our understanding of phage nucleus and PhuZ spindle diversity and function, providing a roadmap for identifying key mechanisms underlying nucleus-based phage replication.

Breakage fusion bridge cycles drive high oncogene copy number, but not intratumoral genetic heterogeneity or rapid cancer genome change.

Oncogene amplification is a major driver of cancer pathogenesis. Breakage fusion bridge (BFB) cycles, like extrachromosomal DNA (ecDNA), can lead to high copy numbers of oncogenes, but their impact on intratumoral heterogeneity, treatment response, and patient survival are not well understood due to difficulty in detecting them by DNA sequencing. We describe a novel algorithm that detects and reconstructs BFB amplifications using optical genome maps (OGMs), called OM2BFB. OM2BFB showed high precision (>93%) and recall (92%) in detecting BFB amplifications in cancer cell lines, PDX models and primary tumors. OM-based comparisons demonstrated that short-read BFB detection using our AmpliconSuite (AS) toolkit also achieved high precision, albeit with reduced sensitivity. We detected 371 BFB events using whole genome sequences from 2,557 primary tumors and cancer lines. BFB amplifications were preferentially found in cervical, head and neck, lung, and esophageal cancers, but rarely in brain cancers. BFB amplified genes show lower variance of gene expression, with fewer options for regulatory rewiring relative to ecDNA amplified genes. BFB positive (BFB (+)) tumors showed reduced heterogeneity of amplicon structures, and delayed onset of resistance, relative to ecDNA(+) tumors. EcDNA and BFB amplifications represent contrasting mechanisms to increase the copy numbers of oncogene with markedly different characteristics that suggest different routes for intervention.

Circular extrachromosomal DNA promotes tumor heterogeneity in high-risk medulloblastoma.

Circular extrachromosomal DNA (ecDNA) in patient tumors is an important driver of oncogenic gene expression, evolution of drug resistance and poor patient outcomes. Applying computational methods for the detection and reconstruction of ecDNA across a retrospective cohort of 481 medulloblastoma tumors from 465 patients, we identify circular ecDNA in 82 patients (18%). Patients with ecDNA-positive medulloblastoma were more than twice as likely to relapse and three times as likely to die within 5 years of diagnosis. A subset of tumors harbored multiple ecDNA lineages, each containing distinct amplified oncogenes. Multimodal sequencing, imaging and CRISPR inhibition experiments in medulloblastoma models reveal intratumoral heterogeneity of ecDNA copy number per cell and frequent putative 'enhancer rewiring' events on ecDNA. This study reveals the frequency and diversity of ecDNA in medulloblastoma, stratified into molecular subgroups, and suggests copy number heterogeneity and enhancer rewiring as oncogenic features of ecDNA.