Genome Biology最新文献

We introduce SCITO-seq2, an enhanced successor to SCITO-seq that integrates probe-based RNA detection with the established ultra-high-throughput protein profiling. SCITO-seq2 achieves robust quantification of transcripts and surface proteins across more than 100,000 cells, with a shared pool barcoding strategy ensuring precise matching of molecular profiles within multiplexed droplets. SCITO-seq2 is compatible with cell hashing technology, allowing efficient sample multiplexing. We demonstrate its utility in autoimmune diseases, including childhood systemic lupus erythematosus and CTLA4 haploinsufficiency with autoimmune infiltration, enabling the detection of minor immune clusters and disease-specific protein signatures. This platform establishes a scalable, streamlined, and cost-effective next-generation single-cell multi-omics workflow.

Nanopore sequencing enables real-time, long-read analysis by processing raw signals as they are produced. A key step, segmentation of signals into events, is typically handled algorithmically, struggling in noisy regions. We present Campolina, a first deep-learning framework for accurate segmentation of raw nanopore signals. Campolina uses a convolutional model to identify event boundaries and significantly outperforms the traditional Scrappie algorithm on R9.4.1 and R10.4.1 datasets. We introduce a comprehensive evaluation pipeline and show that Campolina aligns better with reference-guided ground-truth segmentation. We show that integrating Campolina segmentation into real-time frameworks, Sigmoni and RawHash2, improves their performance while maintaining time efficiency.

Background: IscB (Insertion sequences Cas9-like OrfB) represents a novel class of RNA-guided nucleases, approximately one-third the size of Cas9 proteins. Despite the limited natural efficiency in eukaryotic cells, recent advances have led to the engineering of several IscBs for mammalian genome editing.

Results: In this study, we screen and identify high-activity IscB variants for rice. A version of pIscB-v3, combining enOgeuIscB and ωRNA-v13, demonstrated superior mutagenesis efficiency compared to other systems. The average editing efficiency of pIscB-v3 is 17.61% from ten endogenous targets, and we obtain edited lines in up to 83.33% of T0 generation with 33.33% of homozygous and bi-allelic mutations. Further analysis reveals that pIscB-v3 exhibits high editing specificity and relaxed target-adjacent motif (TAM) compatibility in rice. Beyond gene knockout systems, we develop cytosine base editors (CBEs) and adenine base editors (ABEs) from pIscB-v3. We find that the ssDNA-targeting SCP1.201 family deaminase Sdd7 outperformed human APOBEC3A in IscB-CBEs for C-to-T conversions in rice. The Sdd7-nIscB achieves precise edits in 22.92% of lines on average, with a maximum frequency of 47.92%. Additionally, TadA8e-nIscB exhibits limited activity. However, fusing an extra copy of TadA-8e to either terminus of TadA8e-nIsc significantly enhances A-to-G conversions.

Conclusions: Collectively, our results demonstrate the robust capabilities of IscB to develop an efficient and versatile miniature plant genome editing toolkit to substantially facilitate crop breeding.

Background: Processing of archaeal 16S and 23S rRNAs is believed to involve excision of individual rRNAs from polycistronic precursors, circularization of excised rRNAs, and re-linearization before the incorporation into ribosomes. However, all the knowledge is derived from several isolated species, leaving open the possibility that different processes may occur in other archaeal groups.

Results: Here, we investigate rRNAs from diverse and mostly uncultivated archaea. Sequencing of total cellular RNA from eight phylum-level lineages indicates that archaeal circular 23S rRNA transcript abundances vastly exceed those of linear counterparts, and linear versions are often undetectable. As the majority of rRNAs derive from mature ribosomes, the data suggest that ribosomes contain circular 23S rRNAs. Thus, we directly sequence RNA extracted from isolated ribosomes of a model archaeon, Methanosarcina acetivorans, and confirm that the 23S rRNAs in the ribosomes are circular. Structural modeling places the 5' and 3' ends of the linear precursors of archaeal 23S rRNAs in close proximity to form a GNRA tetraloop (in which N is A, C, G, or U and R is A or G), consistent with their existence as circular molecules. We also confirm the existence of circular 16S rRNA intermediates in transcriptomes of most archaea, yet a circular form is not evident in some distinct archaeal groups, suggesting that certain archaea do not circularize 16S rRNA during processing.

Conclusions: Our findings uncover unexpected variations in the processing required to generate mature rRNAs and the conformation of functional molecules in archaeal ribosomes.

Background: The viticulture has long suffered from the downy mildew caused by Plasmopara viticola, a strictly obligate biotrophic oomycete. Numerous studies have been performed to reveal how grapevine defends against Plasmopara viticola, but they mainly investigate the plant defense responses on the whole tissue level, not on the cellular level.

Results: Here we employ single-cell RNA sequencing and spatial RNA sequencing to profile approximately 100,000 individual cells (~ 89,000 from scRNA-seq and ~ 11,000 from spRNA-seq), generating the first single-cell transcriptome atlas of grapevine leaves during Plasmopara viticola infection. This high-resolution atlas reveals the dynamic and distinct defense responses of plant cells at early stages of oomycete infection. Notably, we find that Plasmopara viticola reprograms the guard cell transcriptome to facilitate successful invasion, likely by altering the expression of ABA negative regulators and modulating a potassium channel regulatory pathway to influence stomatal movement.

Conclusions: Overall, our work reveals differential and dynamic responses of grapevine to the Plasmopara viticola infection at a single-cell level, providing valuable clues for dissecting the interaction between plants and oomycetes.

Deep mutational scanning (DMS) coupled with fluorescence-activated cell sorting (FACS) provides a high-throughput method to link genetic variants with quantitative molecular phenotypes. Analysis of these experiments is challenging due to measurement variance and the multidimensional FACS readout. However, no statistical method has yet been developed to address these challenges. Here we present Lilace, a Bayesian statistical model to estimate variant effects with uncertainty quantification from FACS-based DMS experiments. We validate Lilace's performance and robustness using simulated data and apply it to OCT1 and Kir2.1 DMS datasets, demonstrating an improved false discovery rate while largely maintaining sensitivity.

Background: Whole genome sequence (WGS) data in multi-ancestry samples supports discovery of low-frequency or population-specific genetic variants associated with chronic obstructive pulmonary disease (COPD) and lung function.

Results: We performed single variant, structural variant, and gene-based analysis of pulmonary function (FEV₁, FVC and FEV₁/FVC) and COPD case-control status in 44,287 multi-ancestry participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We validated findings using the UK Biobank and assessed implicated genes using lung single-cell RNA-seq (scRNA-seq) data sets. Applying a genome-wide significance threshold (P < 5 × 10^-9), we replicated known loci and identified novel associations near LY86, MAGI1, GRK7, and LINC02668. Colocalization with gene expression quantitative trait loci (eQTL) from the Lung Tissue Research Consortium highlighted known candidate genes including ADAM19, THSD4, C4B, and PSMA4, which were not identified through other eQTL sources. Multi-ancestry analysis improved fine-mapping resolution (e.g., HTR4 and RIN3). Gene-based analysis identified and replicated HMCN1. In human lung scRNA-seq data sets, lung epithelial cells and immune cell types showed enriched expression, while fibroblasts showed higher expression for HMCN1. CRISPR targeting HMCN1 in IMR90 demonstrated reduced expression of collagen genes.

Conclusions: Large-scale multi-ancestry WGS analysis improves variant discovery and fine-mapping resolution for lung function and COPD and highlights biologically relevant genes and pathways.

Background: Non-negative matrix factorization is a powerful linear algebra tool used in multiple areas of data analysis, including computational biology. Despite numerous optimization methods devised for non-negative matrix factorization, our understanding of the inherent topological structure within factorizable matrices remains limited.

Results: This study reveals the topological properties of linear mixture data, leading to a remarkable reduction of the non-negative matrix factorization optimization problem to a search for K(K-1) variables, where K represents the number of pure components, regardless of the initial matrix size. This is achieved by revealing complementary simplex structures existing in both feature and sample spaces and leveraging the Sinkhorn transformation to find the relationship between these simplexes. We validate this approach in the context of an unconstrained mixed images scenario and achieve a significant improvement in decomposition accuracy. Furthermore, we successfully applied the proposed approach in the biological context of bulk RNA-seq gene expression deconvolution.

Conclusions: The Dual Simplex unified analytical framework improves robustness to noise and enhances optimization stability, enabling accurate recovery of component proportions and expression profiles. Importantly, the framework naturally accommodates both reference-free and marker-based deconvolution settings, providing a general and efficient solution for analyzing complex biological mixtures such as bulk RNA-seq and single-cell derived data.