q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics.

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

SCING: Inference of robust, interpretable gene regulatory networks from single cell and spatial transcriptomics.

Gene regulatory network (GRN) inference is an integral part of understanding physiology and disease. Single cell/nuclei RNA-seq (scRNA-seq/snRNA-seq) data has been used to elucidate cell-type GRNs; however, the accuracy and speed of current scRNAseq-based GRN approaches are suboptimal. Here, we present Single Cell INtegrative Gene regulatory network inference (SCING), a gradient boosting and mutual information-based approach for identifying robust GRNs from scRNA-seq, snRNA-seq, and spatial transcriptomics data. Performance evaluation using Perturb-seq datasets, held-out data, and the mouse cell atlas combined with the DisGeNET database demonstrates the improved accuracy and biological interpretability of SCING compared to existing methods. We applied SCING to the entire mouse single cell atlas, human Alzheimer's disease (AD), and mouse AD spatial transcriptomics. SCING GRNs reveal unique disease subnetwork modeling capabilities, have intrinsic capacity to correct for batch effects, retrieve disease relevant genes and pathways, and are informative on spatial specificity of disease pathogenesis.

Tiered sympathetic control of cardiac function revealed by viral tracing and single cell transcriptome profiling.

The cell bodies of postganglionic sympathetic neurons innervating the heart primarily reside in the stellate ganglion (SG), alongside neurons innervating other organs and tissues. Whether cardiac-innervating stellate ganglionic neurons (SGNs) exhibit diversity and distinction from those innervating other tissues is not known. To identify and resolve the transcriptomic profiles of SGNs innervating the heart, we leveraged retrograde tracing techniques using adeno-associated virus (AAV) expressing fluorescent proteins (GFP or Td-tomato) with single cell RNA sequencing. We investigated electrophysiologic, morphologic, and physiologic roles for subsets of cardiac-specific neurons and found that three of five adrenergic SGN subtypes innervate the heart. These three subtypes stratify into two subpopulations; high (NA1a) and low (NA1b and NA1c) neuropeptide-Y (NPY) -expressing cells, exhibit distinct morphological, neurochemical, and electrophysiologic characteristics. In physiologic studies in transgenic mouse models modulating NPY signaling, we identified differential control of cardiac responses by these two subpopulations to high and low stress states. These findings provide novel insights into the unique properties of neurons responsible for cardiac sympathetic regulation, with implications for novel strategies to target specific neuronal subtypes for sympathetic blockade in cardiac disease.

Tiered Sympathetic Control of Cardiac Function Revealed by Viral Tracing and Single Cell Transcriptome Profiling.

The cell bodies of postganglionic sympathetic neurons innervating the heart primarily reside in the stellate ganglion (SG), alongside neurons innervating other organs and tissues. Whether cardiac-innervating stellate ganglionic neurons (SGNs) exhibit diversity and distinction from those innervating other tissues is not known. To identify and resolve the transcriptomic profiles of SGNs innervating the heart we leveraged retrograde tracing techniques using adeno-associated virus (AAV) expressing fluorescent proteins (GFP or Td-tomato) with single cell RNA sequencing. We investigated electrophysiologic, morphologic, and physiologic roles for subsets of cardiac-specific neurons and found that three of five adrenergic SGN subtypes innervate the heart. These three subtypes stratify into two subpopulations; high (NA1a) and low (NA1b and NA1c) Npy-expressing cells, exhibit distinct morphological, neurochemical, and electrophysiologic characteristics. In physiologic studies in transgenic mouse models modulating NPY signaling, we identified differential control of cardiac responses by these two subpopulations to high and low stress states. These findings provide novel insights into the unique properties of neurons responsible for cardiac sympathetic regulation, with implications for novel strategies to target specific neuronal subtypes for sympathetic blockade in cardiac disease.

A comprehensive benchmarking of WGS-based deletion structural variant callers.

Advances in whole-genome sequencing (WGS) promise to enable the accurate and comprehensive structural variant (SV) discovery. Dissecting SVs from WGS data presents a substantial number of challenges and a plethora of SV detection methods have been developed. Currently, evidence that investigators can use to select appropriate SV detection tools is lacking. In this article, we have evaluated the performance of SV detection tools on mouse and human WGS data using a comprehensive polymerase chain reaction-confirmed gold standard set of SVs and the genome-in-a-bottle variant set, respectively. In contrast to the previous benchmarking studies, our gold standard dataset included a complete set of SVs allowing us to report both precision and sensitivity rates of the SV detection methods. Our study investigates the ability of the methods to detect deletions, thus providing an optimistic estimate of SV detection performance as the SV detection methods that fail to detect deletions are likely to miss more complex SVs. We found that SV detection tools varied widely in their performance, with several methods providing a good balance between sensitivity and precision. Additionally, we have determined the SV callers best suited for low- and ultralow-pass sequencing data as well as for different deletion length categories.

Joint cell segmentation and cell type annotation for spatial transcriptomics.

RNA hybridization-based spatial transcriptomics provides unparalleled detection sensitivity. However, inaccuracies in segmentation of image volumes into cells cause misassignment of mRNAs which is a major source of errors. Here, we develop JSTA, a computational framework for joint cell segmentation and cell type annotation that utilizes prior knowledge of cell type-specific gene expression. Simulation results show that leveraging existing cell type taxonomy increases RNA assignment accuracy by more than 45%. Using JSTA, we were able to classify cells in the mouse hippocampus into 133 (sub)types revealing the spatial organization of CA1, CA3, and Sst neuron subtypes. Analysis of within cell subtype spatial differential gene expression of 80 candidate genes identified 63 with statistically significant spatial differential gene expression across 61 (sub)types. Overall, our work demonstrates that known cell type expression patterns can be leveraged to improve the accuracy of RNA hybridization-based spatial transcriptomics while providing highly granular cell (sub)type information. The large number of newly discovered spatial gene expression patterns substantiates the need for accurate spatial transcriptomic measurements that can provide information beyond cell (sub)type labels.

Single-cell transcriptomic profiling of satellite glial cells in stellate ganglia reveals developmental and functional axial dynamics.

Stellate ganglion neurons, important mediators of cardiopulmonary neurotransmission, are surrounded by satellite glial cells (SGCs), which are essential for the function, maintenance, and development of neurons. However, it remains unknown whether SGCs in adult sympathetic ganglia exhibit any functional diversity, and what role this plays in modulating neurotransmission. We performed single-cell RNA sequencing of mouse stellate ganglia (n = 8 animals), focusing on SGCs (n = 11,595 cells). SGCs were identified by high expression of glial-specific transcripts, S100b and Fabp7. Microglia and Schwann cells were identified by expression of C1qa/C1qb/C1qc and Ncmap/Drp2, respectively, and excluded from further analysis. Dimensionality reduction and clustering of SGCs revealed six distinct transcriptomic subtypes, one of which was characterized the expression of pro-inflammatory markers and excluded from further analyses. The transcriptomic profiles and corresponding biochemical pathways of the remaining subtypes were analyzed and compared with published astrocytic transcriptomes. This revealed gradual shifts of developmental and functional pathways across the subtypes, originating from an immature and pluripotent subpopulation into two mature populations of SGCs, characterized by upregulated functional pathways such as cholesterol metabolism. As SGCs aged, these functional pathways were downregulated while genes and pathways associated with cellular stress responses were upregulated. These findings were confirmed and furthered by an unbiased pseudo-time analysis, which revealed two distinct trajectories involving the five subtypes that were studied. These findings demonstrate that SGCs in mouse stellate ganglia exhibit transcriptomic heterogeneity along maturation or differentiation axes. These subpopulations and their unique biochemical properties suggest dynamic physiological adaptations that modulate neuronal function.

Benchmarking of computational error-correction methods for next-generation sequencing data.

BACKGROUND: Recent advancements in next-generation sequencing have rapidly improved our ability to study genomic material at an unprecedented scale. Despite substantial improvements in sequencing technologies, errors present in the data still risk confounding downstream analysis and limiting the applicability of sequencing technologies in clinical tools. Computational error correction promises to eliminate sequencing errors, but the relative accuracy of error correction algorithms remains unknown. RESULTS: In this paper, we evaluate the ability of error correction algorithms to fix errors across different types of datasets that contain various levels of heterogeneity. We highlight the advantages and limitations of computational error correction techniques across different domains of biology, including immunogenomics and virology. To demonstrate the efficacy of our technique, we apply the UMI-based high-fidelity sequencing protocol to eliminate sequencing errors from both simulated data and the raw reads. We then perform a realistic evaluation of error-correction methods. CONCLUSIONS: In terms of accuracy, we find that method performance varies substantially across different types of datasets with no single method performing best on all types of examined data. Finally, we also identify the techniques that offer a good balance between precision and sensitivity.

Therapeutic IDOL Reduction Ameliorates Amyloidosis and Improves Cognitive Function in APP/PS1 Mice.

Brain lipoprotein receptors have been shown to regulate the metabolism of ApoE and ß-amyloid (Aß) and are potential therapeutic targets for Alzheimer's disease (AD). Previously, we identified E3 ubiquitin ligase IDOL as a negative regulator of brain lipoprotein receptors. Genetic ablation of Idol increases low-density lipoprotein receptor protein levels, which facilitates Aß uptake and clearance by microglia. In this study, we utilized an antisense oligonucleotide (ASO) to reduce IDOL expression therapeutically in the brains of APP/PS1 male mice. ASO treatment led to decreased Aß pathology and improved spatial learning and memory. Single-cell transcriptomic analysis of hippocampus revealed that IDOL inhibition upregulated lysosomal/phagocytic genes in microglia. Furthermore, clustering of microglia revealed that IDOL-ASO treatment shifted the composition of the microglia population by increasing the prevalence of disease-associated microglia. Our results suggest that reducing IDOL expression in the adult brain promotes the phagocytic clearance of Aß and ameliorates Aß-dependent pathology. Pharmacological inhibition of IDOL activity in the brain may represent a therapeutic strategy for the treatment of AD.

Challenges and recommendations to improve the installability and archival stability of omics computational tools.

Developing new software tools for analysis of large-scale biological data is a key component of advancing modern biomedical research. Scientific reproduction of published findings requires running computational tools on data generated by such studies, yet little attention is presently allocated to the installability and archival stability of computational software tools. Scientific journals require data and code sharing, but none currently require authors to guarantee the continuing functionality of newly published tools. We have estimated the archival stability of computational biology software tools by performing an empirical analysis of the internet presence for 36,702 omics software resources published from 2005 to 2017. We found that almost 28% of all resources are currently not accessible through uniform resource locators (URLs) published in the paper they first appeared in. Among the 98 software tools selected for our installability test, 51% were deemed "easy to install," and 28% of the tools failed to be installed at all because of problems in the implementation. Moreover, for papers introducing new software, we found that the number of citations significantly increased when authors provided an easy installation process. We propose for incorporation into journal policy several practical solutions for increasing the widespread installability and archival stability of published bioinformatics software.