Alevin efficiently estimates accurate gene abundances from dscRNA-seq data

Novel Porcine Model Reveals Two Distinct LGR5 Cell Types during Lung Development and Homeostasis

Cells expressing leucine-rich repeat-containing G-protein-coupled receptor 5 (LGR5) play a pivotal role in homeostasis, repair, and regeneration in multiple organs, including skin and gastrointestinal tract, but little is known about their role in the lung. Findings from mice, a widely used animal model, suggest that lung LGR5 expression differs from that of humans. In this work, using a new transgenic pig model, we identify two main populations of LGR5+ cells in the lung that are conserved in human but not mouse lungs. Using RNA sequencing, three-dimensional imaging, and organoid models, we determine that in the fetal lung, epithelial LGR5 expression is transient in a subpopulation of SOX9+/ETV5+/SFTPC+ progenitor lung tip cells. In contrast, epithelial LGR5 expression is absent from postnatal lung but is reactivated in bronchioalveolar organoids derived from basal airway cells. We also describe a separate population of mesenchymal LGR5+ cells that surrounds developing and mature airways, lies adjacent to airway basal cells, and is closely associated with nerve fibers. Transcriptionally, mesenchymal LGR5+ cells include a subset of peribronchial fibroblasts that express unique patterns of SHH, FGF, WNT, and TGF-β signaling pathway genes. These results support distinct roles for LGR5+ cells in the lung and describe a physiologically relevant animal model for further studies on the function of these cells in repair and regeneration.

Modeling the dynamics of EMT reveals genes associated with pan-cancer intermediate states and plasticity

Epithelial-mesenchymal transition (EMT) is a cell state transition co-opted by cancer that drives metastasis via stable intermediate states. Here we study EMT dynamics to identify marker genes of highly metastatic intermediate cells via mathematical modeling with single-cell RNA sequencing (scRNA-seq) data. Across multiple tumor types and stimuli, we identified genes consistently upregulated in EMT intermediate states, many previously unrecognized as EMT markers. Bayesian parameter inference of a simple EMT mathematical model revealed tumor-specific transition rates, providing a framework to quantify EMT progression. Consensus analysis of differential expression, RNA velocity, and model-derived dynamics highlighted SFN and NRG1 as key regulators of intermediate EMT. Independent validation confirmed SFN as an intermediate state marker. Our approach integrates modeling and inference to identify genes associated with EMT dynamics, offering biomarkers and therapeutic targets to modulate tumor-promoting cell state transitions driven by EMT.

T-bet+ CXCR3+ B cells drive hyperreactive B-T cell interactions in multiple sclerosis

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system (CNS). Self-peptide-dependent autoproliferation (AP) of B and T cells is a key mechanism in MS. Here, we show that pro-inflammatory B-T cell-enriched cell clusters (BTECs) form during AP and mirror features of a germinal center reaction. T-bet+CXCR3+ B cells are the main cell subset amplifying and sustaining their counterpart Th1 cells via interferon (IFN)-γ and are present in highly inflamed meningeal tissue. The underlying B cell activation signature is reflected by epigenetic modifications and receptor-ligand interactions with self-reactive T cells. AP+ CXCR3+ B cells show marked clonal evolution from memory to somatically hypermutated plasmablasts and upregulation of IFN-γ-related genes. Our data underscore a key role of T-bet+CXCR3+ B cells in the pathogenesis of MS in both the peripheral immune system and the CNS compartment, and thus they appear to be involved in both early relapsing-remitting disease and the chronic stage.

Microbial-mediated carbon metabolism in the subtropical marine mangroves affected by shrimp pond discharge

Mangrove ecosystems exhibit high efficiency in carbon (C) sequestering within the global ecosystem. However, the rapid expansion of the shrimp farming industry poses a significant threat to these delicate ecosystems. The microbial mechanisms driving C metabolism in shrimp-affected sediments remain poorly understood. This study investigates the spatiotemporal dynamics of C metabolism-related microbial communities in shrimp pond and natural mangrove sediments in a subtropical region. Shrimp pond discharge altered soil properties, microbial diversity, and microbial stability, driven by factors such as salinity, sulfide, and total organic C (TOC). Metagenomic analyses reveals shifts in C degradation and oxidation, with a reduction in genes for cellulose and hemicellulose degradation. Microbial markers like Prolixibacteraceae and Nitrosopumilaceae reflect these changes. Co-occurrence network analysis indicates higher connectivity within shrimp pond groups, suggesting nutrient-driven changes in symbiotic relationships. PLS-PM analysis further confirms the interplay between microbial composition, nutrient levels, and C metabolism, with higher 16S rRNA operon copy numbers linked to increased C fixation. These findings demonstrate how shrimp pond discharge alters microbial networks and C metabolism, with implications for ecosystem resilience.

kallisto, bustools and kb-python for quantifying bulk, single-cell and single-nucleus RNA-seq

The term 'RNA-seq' refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, single cells or single nuclei. The kallisto, bustools and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data. Execution of this protocol requires basic familiarity with a command line environment. With this protocol, quantification of a moderately sized RNA-seq dataset can be completed within minutes.

A single-nucleus RNA sequencing atlas of the postnatal retina of the shark Scyliorhinus canicula

The retina, whose basic cellular structure is highly conserved across vertebrates, constitutes an accessible system for studying the central nervous system. In recent years, single-cell RNA sequencing studies have uncovered cellular diversity in the retina of a variety of species, providing new insights on retinal evolution and development. However, similar data in cartilaginous fishes, the sister group to all other extant jawed vertebrates, are still lacking. Here, we present a single-nucleus RNA sequencing atlas of the postnatal retina of the catshark Scyliorhinus canicula, consisting of the expression profiles for 17,438 individual cells from three female, juvenile catshark specimens. Unsupervised clustering revealed 22 distinct cell types comprising all major retinal cell classes, as well as retinal progenitor cells (whose presence reflects the persistence of proliferative activity in postnatal stages in sharks) and oligodendrocytes. Thus, our dataset serves as a foundation for further studies on the development and function of the catshark retina. Moreover, integration of our atlas with data from other species will allow for a better understanding of vertebrate retinal evolution.

Retrotrans-genomics identifies aberrant THE1B endogenous retrovirus fusion transcripts in the pathogenesis of sarcoidosis

Transposon-like human element 1B (THE1B) originates from ancient retroviral sequences integrated into the primate genome approximately 50 million years ago, now accounting for at least 27,233 copies in the human genome, suggesting their extensive influence on human genomic architecture. Here we report identification of 19 THE1B fusion transcripts through short- and long-read RNA-seq analysis, 15 of which are previously unmapped, showing elevated expression in 16 individuals with sarcoid myopathy (SM), as compared to 400 controls with various other muscle diseases. Analysis of publicly available RNA-seq data indicated a correlation between the reduced expression of eight THE1B fusion transcripts and clinical improvement in individuals with cutaneous sarcoidosis receiving tofacitinib treatment. Single-cell or single-nucleus RNA-seq analyses of sarcoidosis not only confirmed these transcripts but also revealed a novel read-through transcript, SIRPB1-SIRPD, and TREM2.1, predominantly in granuloma-associated macrophages. The expression profiles of THE1B fusion transcripts in tuberculosis (TB) significantly differed from SM in single-cell RNA-seq data, suggesting that the differences between TB's caseous granulomas and sarcoidosis's non-caseous granulomas might be linked to disparate expression patterns of THE1B fusion transcripts. Our retrotrans-genomics approach has not only identified the genomic landscape of sarcoidosis but also provided new insights into its etiology.

Phospho-seq: integrated, multi-modal profiling of intracellular protein dynamics in single cells

Cell signaling plays a critical role in neurodevelopment, regulating cellular behavior and fate. While multimodal single-cell sequencing technologies are rapidly advancing, scalable and flexible profiling of cell signaling states alongside other molecular modalities remains challenging. Here we present Phospho-seq, an integrated approach that aims to quantify cytoplasmic and nuclear proteins, including those with post-translational modifications, and to connect their activity with cis-regulatory elements and transcriptional targets. We utilize a simplified benchtop antibody conjugation method to create large custom neuro-focused antibody panels for simultaneous protein and scATAC-seq profiling on whole cells, alongside both experimental and computational strategies to incorporate transcriptomic measurements. We apply our workflow to cell lines, induced pluripotent stem cells, and months-old retinal and brain organoids to demonstrate its broad applicability. We show that Phospho-seq can provide insights into cellular states and trajectories, shed light on gene regulatory relationships, and help explore the causes and effects of diverse cell signaling in neurodevelopment.

Single-cell multiomics to advance cell therapy

Single-cell RNA-sequencing (scRNAseq) was first introduced in 2009 and has evolved with many technological advancements over the last decade. Not only are there several scRNAseq platforms differing in many aspects, but there are also a large number of computational pipelines available for downstream analyses which are being developed at an exponential rate. Such computational data appear in many scientific publications in virtually every field of study; thus, investigators should be able to understand and interpret data in this rapidly evolving field. Here, we discuss key differences in scRNAseq platforms, crucial steps in scRNAseq experiments, standard downstream analyses and introduce newly developed multimodal approaches. We then discuss how single-cell omics has been applied to advance the field of cell therapy.

Genome-wide transcriptomics revealed carbon source-mediated gamma-aminobutyric acid (GABA) production in a probiotic, Lactiplantibacillus pentosus 9D3

GABA-producing probiotics present promising opportunities for developing functional foods. Carbon sources have been identified as a critical influence on GABA production. Therefore, this study investigated the holistic metabolic responses and GABA biosynthesis to various carbon sources of Lactiplantibacillus pentosus 9D3, a proficient GABA producer, using a genome-wide transcriptomic approach. The analysis revealed 414 genes with differential expression responses to altering carbon sources, i.e., glucose, sucrose, and lactose, notably sugar phosphotransferase systems (PTS) (11 genes), indicating carbon source-mediated transcriptional change patterns in L. pentosus 9D3. The integration of transcriptome data with a genome-scale metabolic network (GSMN) revealed that L. pentosus 9D3 displays adaptability by synthesizing GABA as an alternative acid-tolerant mechanism when lactose is used as a carbon source rather than depending on the fatty acid synthesis and the arginine catabolic pathway. The findings of this study offer valuable insights into optimal carbon source utilization and gene expression co-regulation, thereby enhancing the GABA-producing capability of a probiotic and broadening its potential applications in the functional food industry.

Characterizing the omics landscape based on 10,000+ datasets

The characteristics of data produced by omics technologies are pivotal, as they critically influence the feasibility and effectiveness of computational methods applied in downstream analyses, such as data harmonization and differential abundance analyses. Furthermore, variability in these data characteristics across datasets plays a crucial role, leading to diverging outcomes in benchmarking studies, which are essential for guiding the selection of appropriate analysis methods in all omics fields. Additionally, downstream analysis tools are often developed and applied within specific omics communities due to the presumed differences in data characteristics attributed to each omics technology. In this study, we investigate over ten thousand datasets to understand how proteomics, metabolomics, lipidomics, transcriptomics, and microbiome data vary in specific data characteristics. We were able to show patterns of data characteristics specific to the investigated omics types and provide a tool that enables researchers to assess how representative a given omics dataset is for its respective discipline. Moreover, we illustrate how data characteristics can impact analyses at the example of normalization in the presence of sample-dependent proportions of missing values. Given the variability of omics data characteristics, we encourage the systematic inspection of these characteristics in benchmark studies and for downstream analyses to prevent suboptimal method selection and unintended bias.

A Cross-Species Atlas of the Dorsal Vagal Complex Reveals Neural Mediators of Cagrilintide's Effects on Energy Balance

Amylin analogs, including potential anti-obesity therapies like cagrilintide, act on neurons in the brainstem dorsal vagal complex (DVC) that express calcitonin receptors (CALCR). These receptors, often combined with receptor activity-modifying proteins (RAMPs), mediate the suppression of food intake and body weight. To understand the molecular and neural mechanisms of cagrilintide action, we used single-nucleus RNA sequencing to define 89 cell populations across the rat, mouse, and non-human primate caudal brainstem. We then integrated spatial profiling to reveal neuron distribution in the rat DVC. Furthermore, we compared the acute and long-term transcriptional responses to cagrilintide across DVC neurons of rats, which exhibit strong cagrilintide responsiveness, and mice, which respond poorly to cagrilintide over the long term. We found that cagrilintide promoted long-term transcriptional changes, including increased prolactin releasing hormone (Prlh) expression, in the nucleus of the solitary tract (NTS) Calcr/Prlh cells in rats, but not in mice, suggesting the importance of NTS Calcr/Prlh cells for sustained weight loss. Indeed, activating rat area postrema Calcr cells briefly reduced food intake but failed to decrease food intake or body weight over the long term. Overall, these results not only provide a cross-species and spatial atlas of DVC cell populations but also define the molecular and neural mediators of acute and long-term cagrilintide action.

Galaxy as a gateway to bioinformatics: Multi-Interface Galaxy Hands-on Training Suite (MIGHTS) for scRNA-seq

BACKGROUND

Bioinformatics is fundamental to biomedical sciences, but its mastery presents a steep learning curve for bench biologists and clinicians. Learning to code while analyzing data is difficult. The curve may be flattened by separating these two aspects and providing intermediate steps for budding bioinformaticians. Single-cell analysis is in great demand from biologists and biomedical scientists, as evidenced by the proliferation of training events, materials, and collaborative global efforts like the Human Cell Atlas. However, iterative analyses lacking reinstantiation, coupled with unstandardized pipelines, have made effective single-cell training a moving target.

FINDINGS

To address these challenges, we present a Multi-Interface Galaxy Hands-on Training Suite (MIGHTS) for single-cell RNA sequencing (scRNA-seq) analysis, which offers parallel analytical methods using a graphical interface (buttons) or code. With clear, interoperable materials, MIGHTS facilitates smooth transitions between environments. Bridging the biologist-programmer gap, MIGHTS emphasizes interdisciplinary communication for effective learning at all levels. Real-world data analysis in MIGHTS promotes critical thinking and best practices, while FAIR data principles ensure validation of results. MIGHTS is freely available, hosted on the Galaxy Training Network, and leverages Galaxy interfaces for analyses in both settings. Given the ongoing popularity of Python-based (Scanpy) and R-based (Seurat & Monocle) scRNA-seq analyses, MIGHTS enables analyses using both.

CONCLUSIONS

MIGHTS consists of 11 tutorials, including recordings, slide decks, and interactive visualizations, and a demonstrated track record of sustainability via regular updates and community collaborations. Parallel pathways in MIGHTS enable concurrent training of scientists at any programming level, addressing the heterogeneous needs of novice bioinformaticians.

Unique lymphocyte transcriptomic profiles in septic patients with chronic critical illness

Introduction

Despite continued improvement in post-sepsis survival, long term morbidity and mortality remain high. Chronic critical illness (CCI), defined as persistent inflammation and organ injury requiring prolonged intensive care, is a harbinger of poor long-term outcomes in sepsis survivors. Current dogma states that sepsis survivors are immunosuppressed, particularly in CCI. Investigation of this immune suppression in heterogeneous immune populations across distinct clinical trajectories and outcomes, along with limited sampling access, is accessible via single-cell RNA sequencing (scRNA-seq).

Methods

scRNA-seq analysis was performed on healthy subjects (n=12), acutely septic patients at day 4 ± 1 (n=4), and those defined as rapid recovery (n=4) or CCI (n=5) at day 14-21. Differential gene expression and pathway analyses were performed on peripheral blood lymphocytes at both a population and annotated cell subset level. Cellular function was assessed via enzyme-linked immunosorbent spot (ELISpot), cytokine production analysis, and T-cell proliferation assays on an additional cohort of septic patients (19 healthy, 68 acutely septic, 27 rapid recovery and 20 classified as CCI 14-21 days after sepsis onset).

Results

Sepsis survivors that developed CCI exhibited proportional shifts within lymphoid cell populations, with expanded frequency of CD8+ and NK cells. Differential expression and pathway analyses revealed continued activation in T cells and NK cells, with generalized suppression of B-cell function. Both T and NK cell subsets displayed transcriptomic profiles of exhaustion and immunosuppression in CCI, particularly in CD8+ T effector memory (TEM) cells and NK cells. Functional validation of T-cell behavior in an independent cohort demonstrated T cells maintained proliferative responses in vitro yet exhibited a marked loss of cytokine production. IFN-γ production at the acute phase (day 4 ± 1) was significantly reduced in subjects later classified as CCI.

Discussion

Sepsis patients exhibit unique T-, B-, and NK-cell transcriptional patterns that are both time- and clinical trajectory-dependent. These transcriptomic and pathway differences in sepsis patients that develop CCI are associated with exhaustion in CD8+ TEM cells and NK cells. Understanding the specific immune system patterns of different cell subsets after sepsis at a molecular level will be key to the development of personalized immunotherapy and drug-targeting intervention.

Clinical trial registration

https://clinicaltrials.gov/, identifier NCT02276417.

Improved characterization of 3' single-cell RNA-seq libraries with paired-end avidity sequencing

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site or in principle the polyadenylation site. Direct sequencing across this site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of 'avidity base chemistry' DNA sequencing from Element Biosciences to sequence through the primer and enable accurate paired-end read alignment and precise quantification of polyadenylation sites. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The additional sequence enables direct and independent assignment of reads to polyadenylation sites, which bypasses the complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and demonstrate necessary adjustments to adapter trimming and sequence alignment regardless of platform, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site measurement but do not rule out the potential for similar performance from other emerging platforms.

Redistribution of PU.1 partner transcription factor RUNX1 binding secures cell survival during leukemogenesis

Transcription factors (TFs) orchestrating lineage-development often control genes required for cellular survival. However, it is not well understood how cells survive when such TFs are lost, for example in cancer. PU.1 is an essential TF for myeloid fate, and mice with downregulated PU.1 levels develop acute myeloid leukemia (AML). Combining a multi-omics approach with a functional genetic screen, we reveal that PU.1-downregulated cells fundamentally change their survival control from cytokine-driven pathways to overexpression of an autophagy-predominated stem cell gene program, for which we also find evidence in human AML. Control of this program involves redirected chromatin occupancy of the PU.1 partner TF Runx1 to a lineage-inappropriate binding site repertoire. Hence, genomic reallocation of TF binding upon loss of a partner TF can act as a pro-oncogenic failsafe mechanism by sustaining cell survival during leukemogenesis.

A Roadmap for Selecting and Utilizing Optimal Features in scRNA Sequencing Data Analysis for Stem Cell Research: A Comprehensive Review

Stem cells and the cells they produce are unique because they vary from one cell to another. Traditional methods of studying cells often overlook these differences. However, the development of new technologies for studying individual cells has greatly changed biological research in recent years. Among these innovations, single-cell RNA sequencing (scRNA-seq) stands out. This technique allows scientists to examine the activity of genes in each cell, across thousands or even millions of cells. This makes it possible to understand the diversity of cells, identify new types of cells, and see how cells differ across different tissues, individuals, species, times, and conditions. This paper discusses the importance of scRNA-seq and the computational tools and software that are essential for analyzing the vast amounts of data generated by scRNA-seq studies. Our goal is to provide practical advice for bioinformaticians and biologists who are using scRNA-seq to study stem cells. We offer an overview of the scRNA-seq field, including the tools available, how they can be used, and how to present the results of these studies effectively. Our findings include a detailed overview and classification of tools used in scRNA-seq analysis, based on a review of 2,733 scientific publications. This review is complemented by information from the scRNA-tools database, which lists over 1,400 tools for analyzing scRNA-seq data. This database is an invaluable resource for researchers, offering a wide range of options for analyzing their scRNA-seq data.

TET3 regulates terminal cell differentiation at the metabolic level

TET-family members play a critical role in cell fate commitment. Indeed, TET3 is essential to postnatal development due to yet unknown reasons. To define TET3 function in cell differentiation, we have profiled the intestinal epithelium at single-cell level from wild-type and Tet3 knockout mice. We have found that Tet3 is mostly expressed in differentiated enterocytes. In the absence of TET3, enterocytes exhibit an aberrant differentiation trajectory and do not acquire a physiological cell identity due to an impairment in oxidative phosphorylation, specifically due to an ATP synthase assembly deficiency. Moreover, spatial metabolomics analysis has revealed that Tet3 knockout enterocytes exhibit an unphysiological metabolic profile when compared with their wild-type counterparts. In contrast, no metabolic differences have been observed between both genotypes in the stem cell compartment where Tet3 is mainly not expressed. Collectively, our findings suggest a mechanism by which TET3 regulates mitochondrial function and, thus, terminal cell differentiation at the metabolic level.

Uncovering functional lncRNAs by scRNA-seq with ELATUS

Long non-coding RNAs (lncRNAs) play fundamental roles in cellular processes and pathologies, regulating gene expression at multiple levels. Despite being highly cell type-specific, their study at single-cell (sc) level is challenging due to their less accurate annotation and low expression compared to protein-coding genes. Here, we systematically benchmark different preprocessing methods and develop a computational framework, named ELATUS, based on the combination of the pseudoaligner Kallisto with selective functional filtering. ELATUS enhances the detection of functional lncRNAs from scRNA-seq data, detecting their expression with higher concordance than standard methods with the ATAC-seq profiles in single-cell multiome data. Interestingly, the better results of ELATUS are due to its advanced performance with an inaccurate reference annotation such as that of lncRNAs. We independently confirm the expression patterns of cell type-specific lncRNAs exclusively detected with ELATUS and unveil biologically important lncRNAs, such as AL121895.1, a previously undocumented cis-repressor lncRNA, whose role in breast cancer progression is unnoticed by traditional methodologies. Our results emphasize the necessity for an alternative scRNA-seq workflow tailored to lncRNAs that sheds light on the multifaceted roles of lncRNAs.

The Arabidopsis U1 snRNP regulates mRNA 3'-end processing

The removal of introns by the spliceosome is a key gene regulatory mechanism in eukaryotes, with the U1 snRNP subunit playing a crucial role in the early stages of splicing. Studies in metazoans show that the U1 snRNP also conducts splicing-independent functions, but the lack of genetic tools and knowledge about U1 snRNP-associated proteins have limited the study of such splicing-independent functions in plants. Here we describe an RNA-centric approach that identified more than 200 proteins associated with the Arabidopsis U1 snRNP and revealed a tight link to mRNA cleavage and polyadenylation factors. Interestingly, we found that the U1 snRNP protects mRNAs against premature cleavage and polyadenylation within introns-a mechanism known as telescripting in metazoans-while also influencing alternative polyadenylation site selection in 3'-UTRs. Overall, our work provides a comprehensive view of U1 snRNP interactors and reveals novel functions in regulating mRNA 3'-end processing in Arabidopsis, laying the groundwork for understanding non-canonical functions of plant U1 snRNPs.

A global transcriptional atlas of the effect of acute sleep deprivation in the mouse frontal cortex

Sleep deprivation (SD) has negative effects on brain and body function. Sleep problems are prevalent in a variety of disorders, including neurodevelopmental and psychiatric conditions. Thus, understanding the molecular consequences of SD is of fundamental importance in biology. In this study, we present the first simultaneous bulk and single-nuclear RNA sequencing characterization of the effects of SD in the male mouse frontal cortex. We show that SD predominantly affects glutamatergic neurons, specifically in layers 4 and 5, and produces isoform switching of over 1500 genes, particularly those involved in splicing and RNA binding. At both the global and cell-type specific level, SD has a large repressive effect on transcription, downregulating thousands of genes and transcripts. As a resource we provide extensive characterizations of cell-types, genes, transcripts, and pathways affected by SD. We also provide publicly available tutorials aimed at allowing readers adapt analyses performed in this study to their own datasets.

Single-cell RNA sequencing facilitates the elucidation of the complete biosynthesis of the antidepressant hyperforin in St. John's wort

Hyperforin is the compound responsible for the effectiveness of St. John's wort (Hypericum perforatum) as an antidepressant, but its complete biosynthetic pathway remains unknown. Gene discovery based on co-expression analysis of bulk RNA-sequencing data or genome mining failed to discover the missing steps in hyperforin biosynthesis. In this study, we sequenced the 1.54-Gb tetraploid H. perforatum genome assembled into 32 chromosomes with the scaffold N50 value of 42.44 Mb. By single-cell RNA sequencing, we identified a type of cell, "Hyper cells", wherein hyperforin biosynthesis de novo takes place in both the leaves and flowers. Through pathway reconstitution in yeast and tobacco, we identified and characterized four transmembrane prenyltransferases (HpPT1-4) that are localized at the plastid envelope and complete the hyperforin biosynthetic pathway. The hyperforin polycyclic scaffold is created by a reaction cascade involving an irregular isoprenoid coupling and a tandem cyclization. Our findings reveal how and where hyperforin is biosynthesized, enabling synthetic-biology reconstitution of the complete pathway. Thus, this study not only deepens our comprehension of specialized metabolism at the cellular level but also provides strategic guidance for elucidation of the biosynthetic pathways of other specializied metabolites in plants.

ScRNA-seq reveals the spatiotemporal distribution of camptothecin pathway and transposon activity in Camptotheca acuminata shoot apexes and leaves

Camptotheca acuminata Decne., a significant natural source of the anticancer drug camptothecin (CPT), synthesizes CPT through the monoterpene indole alkaloid (MIA) pathway. In this study, we used single-cell RNA sequencing (scRNA-seq) to generate datasets encompassing over 60,000 cells from C. acuminata shoot apexes and leaves. After cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Analysis of MIA pathway gene expression revealed that most of them exhibited heightened expression in proliferating cells (PCs) and vascular cells (VCs). In contrast to MIA biosynthesis in Catharanthus roseus, CPT biosynthesis in C. acuminata did not exhibit multicellular compartmentalization. Some putative genes encoding enzymes and transcription factors (TFs) related to the biosynthesis of CPT and its derivatives were identified through co-expression analysis. These include 19 cytochrome P450 genes, 8 O-methyltransferase (OMT) genes, and 62 TFs. Additionally, these pathway genes exhibited dynamic expression patterns during VC and EC development. Furthermore, by integrating gene and transposable element (TE) expression data, we constructed novel single-cell transcriptome atlases for C. acuminata. This approach significantly facilitated the identification of rare cell types, including peripheral zone cells (PZs). Some TE families displayed cell type specific, tissue specific, or developmental stage-specific expression patterns, suggesting crucial roles for these TEs in cell differentiation and development. Overall, this study not only provides novel insights into CPT biosynthesis and spatial-temporal TE expression characteristics in C. acuminata, but also serves as a valuable resource for further comprehensive investigations into the development and physiology of this species.

TRANSCRIPTOMIC DIFFERENCES IN PERIPHERAL MONOCYTE POPULATIONS IN SEPTIC PATIENTS BASED ON OUTCOME

ABSTRACT

Postsepsis early mortality is being replaced by survivors who experience either a rapid recovery and favorable hospital discharge or the development of chronic critical illness with suboptimal outcomes. The underlying immunological response that determines these clinical trajectories remains poorly defined at the transcriptomic level. As classical and nonclassical monocytes are key leukocytes in both the innate and adaptive immune systems, we sought to delineate the transcriptomic response of these cell types. Using single-cell RNA sequencing and pathway analyses, we identified gene expression patterns between these two groups that are consistent with differences in TNF-α production based on clinical outcome. This may provide therapeutic targets for those at risk for chronic critical illness in order to improve their phenotype/endotype, morbidity, and long-term mortality.

Improved characterization of single-cell RNA-seq libraries with paired-end avidity sequencing

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site, which is expected to be the poly(A) tail or a genomic adenine homopolymer. Direct sequencing across this priming site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of "avidity base chemistry" DNA sequencing from Element Biosciences to sequence through this homopolymer accurately, and the impact of the additional cDNA sequence on read alignment and precise quantification of polyadenylation site usage. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The resulting paired-end alignments enable direct and independent assignment of reads to polyadenylation sites, which bypasses complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and arrive at an adjusted adapter trimming and alignment workflow that significantly improves the alignment of sequence data from Element and Illumina, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site analyses but do not rule out the potential for similar performance from other emerging platforms.

Quantification and statistical modeling of droplet-based single-nucleus RNA-sequencing data

In complex tissues containing cells that are difficult to dissociate, single-nucleus RNA-sequencing (snRNA-seq) has become the preferred experimental technology over single-cell RNA-sequencing (scRNA-seq) to measure gene expression. To accurately model these data in downstream analyses, previous work has shown that droplet-based scRNA-seq data are not zero-inflated, but whether droplet-based snRNA-seq data follow the same probability distributions has not been systematically evaluated. Using pseudonegative control data from nuclei in mouse cortex sequenced with the 10x Genomics Chromium system and mouse kidney sequenced with the DropSeq system, we found that droplet-based snRNA-seq data follow a negative binomial distribution, suggesting that parametric statistical models applied to scRNA-seq are transferable to snRNA-seq. Furthermore, we found that the quantification choices in adapting quantification mapping strategies from scRNA-seq to snRNA-seq can play a significant role in downstream analyses and biological interpretation. In particular, reference transcriptomes that do not include intronic regions result in significantly smaller library sizes and incongruous cell type classifications. We also confirmed the presence of a gene length bias in snRNA-seq data, which we show is present in both exonic and intronic reads, and investigate potential causes for the bias.

Forseti: a mechanistic and predictive model of the splicing status of scRNA-seq reads

MOTIVATION

Short-read single-cell RNA-sequencing (scRNA-seq) has been used to study cellular heterogeneity, cellular fate, and transcriptional dynamics. Modeling splicing dynamics in scRNA-seq data is challenging, with inherent difficulty in even the seemingly straightforward task of elucidating the splicing status of the molecules from which sequenced fragments are drawn. This difficulty arises, in part, from the limited read length and positional biases, which substantially reduce the specificity of the sequenced fragments. As a result, the splicing status of many reads in scRNA-seq is ambiguous because of a lack of definitive evidence. We are therefore in need of methods that can recover the splicing status of ambiguous reads which, in turn, can lead to more accuracy and confidence in downstream analyses.

RESULTS

We develop Forseti, a predictive model to probabilistically assign a splicing status to scRNA-seq reads. Our model has two key components. First, we train a binding affinity model to assign a probability that a given transcriptomic site is used in fragment generation. Second, we fit a robust fragment length distribution model that generalizes well across datasets deriving from different species and tissue types. Forseti combines these two trained models to predict the splicing status of the molecule of origin of reads by scoring putative fragments that associate each alignment of sequenced reads with proximate potential priming sites. Using both simulated and experimental data, we show that our model can precisely predict the splicing status of many reads and identify the true gene origin of multi-gene mapped reads.

AVAILABILITY AND IMPLEMENTATION

Forseti and the code used for producing the results are available at https://github.com/COMBINE-lab/forseti under a BSD 3-clause license.

A temporal developmental map separates human NK cells from noncytotoxic ILCs through clonal and single-cell analysis

ABSTRACT

Natural killer (NK) cells represent the cytotoxic member within the innate lymphoid cell (ILC) family that are important against viral infections and cancer. Although the NK cell emergence from hematopoietic stem and progenitor cells through multiple intermediate stages and the underlying regulatory gene network has been extensively studied in mice, this process is not well characterized in humans. Here, using a temporal in vitro model to reconstruct the developmental trajectory of NK lineage, we identified an ILC-restricted oligopotent stage 3a CD34-CD117+CD161+CD45RA+CD56- progenitor population, that exclusively gave rise to CD56-expressing ILCs in vitro. We also further investigated a previously nonappreciated heterogeneity within the CD56+CD94-NKp44+ subset, phenotypically equivalent to stage 3b population containing both group-1 ILC and RORγt+ ILC3 cells, that could be further separated based on their differential expression of DNAM-1 and CD161 receptors. We confirmed that DNAM-1hi S3b and CD161hiCD117hi ILC3 populations distinctively differed in their expression of effector molecules, cytokine secretion, and cytotoxic activity. Furthermore, analysis of lineage output using DNA-barcode tracing across these stages supported a close developmental relationship between S3b-NK and S4-NK (CD56+CD94+) cells, whereas distant to the ILC3 subset. Cross-referencing gene signatures of culture-derived NK cells and other noncytotoxic ILCs with publicly available data sets validated that these in vitro stages highly resemble transcriptional profiles of respective in vivo ILC counterparts. Finally, by integrating RNA velocity and gene network analysis through single-cell regulatory network inference and clustering we unravel a network of coordinated and highly dynamic regulons driving the cytotoxic NK cell program, as a guide map for future studies on NK cell regulation.

Recapitulation of Perturbed Striatal Gene Expression Dynamics of Donors' Brains With Ventral Forebrain Organoids Derived From the Same Individuals With Schizophrenia

OBJECTIVE

Schizophrenia is a brain disorder that originates during neurodevelopment and has complex genetic and environmental etiologies. Despite decades of clinical evidence of altered striatal function in affected patients, studies examining its cellular and molecular mechanisms in humans are limited. To explore neurodevelopmental alterations in the striatum associated with schizophrenia, the authors established a method for the differentiation of induced pluripotent stem cells (iPSCs) into ventral forebrain organoids (VFOs).

METHODS

VFOs were generated from postmortem dural fibroblast-derived iPSCs of four individuals with schizophrenia and four neurotypical control individuals for whom postmortem caudate genotypes and transcriptomic data were profiled in the BrainSeq neurogenomics consortium. Individuals were selected such that the two groups had nonoverlapping schizophrenia polygenic risk scores (PRSs).

RESULTS

Single-cell RNA sequencing analyses of VFOs revealed differences in developmental trajectory between schizophrenia and control individuals in which inhibitory neuronal cells from the patients exhibited accelerated maturation. Furthermore, upregulated genes in inhibitory neurons in schizophrenia VFOs showed a significant overlap with upregulated genes in postmortem caudate tissue of individuals with schizophrenia compared with control individuals, including the donors of the iPSC cohort.

CONCLUSIONS

The findings suggest that striatal neurons derived from high-PRS individuals with schizophrenia carry abnormalities that originated during early brain development and that the VFO model can recapitulate disease-relevant cell type-specific neurodevelopmental phenotypes in a dish.

Single-cell transcriptome profiling reveals the spatiotemporal distribution of triterpenoid saponin biosynthesis and transposable element activity in Gynostemma pentaphyllum shoot apexes and leaves

Gynostemma pentaphyllum (Thunb.) Makino is an important producer of dammarene-type triterpenoid saponins. These saponins (gypenosides) exhibit diverse pharmacological benefits such as anticancer, antidiabetic, and immunomodulatory effects, and have major potential in the pharmaceutical and health care industries. Here, we employed single-cell RNA sequencing (scRNA-seq) to profile the transcriptomes of more than 50,000 cells derived from G. pentaphyllum shoot apexes and leaves. Following cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Each cell type displayed substantial transcriptomic heterogeneity both within and between tissues. Examining gene expression patterns across various cell types revealed that gypenoside biosynthesis predominantly occurred in mesophyll cells, with heightened activity observed in shoot apexes compared to leaves. Furthermore, we explored the impact of transposable elements (TEs) on G. pentaphyllum transcriptomic landscapes. Our findings the highlighted the unbalanced expression of certain TE families across different cell types in shoot apexes and leaves, marking the first investigation of TE expression at the single-cell level in plants. Additionally, we observed dynamic expression of genes involved in gypenoside biosynthesis and specific TE families during epidermal and vascular cell development. The involvement of TE expression in regulating cell differentiation and gypenoside biosynthesis warrant further exploration. Overall, this study not only provides new insights into the spatiotemporal organization of gypenoside biosynthesis and TE activity in G. pentaphyllum shoot apexes and leaves but also offers valuable cellular and genetic resources for a deeper understanding of developmental and physiological processes at single-cell resolution in this species.

The post-septic peripheral myeloid compartment reveals unexpected diversity in myeloid-derived suppressor cells

Introduction

Sepsis engenders distinct host immunologic changes that include the expansion of myeloid-derived suppressor cells (MDSCs). These cells play a physiologic role in tempering acute inflammatory responses but can persist in patients who develop chronic critical illness.

Methods

Cellular Indexing of Transcriptomes and Epitopes by Sequencing and transcriptomic analysis are used to describe MDSC subpopulations based on differential gene expression, RNA velocities, and biologic process clustering.

Results

We identify a unique lineage and differentiation pathway for MDSCs after sepsis and describe a novel MDSC subpopulation. Additionally, we report that the heterogeneous response of the myeloid compartment of blood to sepsis is dependent on clinical outcome.

Discussion

The origins and lineage of these MDSC subpopulations were previously assumed to be discrete and unidirectional; however, these cells exhibit a dynamic phenotype with considerable plasticity.

Hodge Decomposition of Single-Cell RNA Velocity

RNA velocity has the ability to capture the cell dynamic information in the biological processes; yet, a comprehensive analysis of the cell state transitions and their associated chemical and biological processes remains a gap. In this work, we provide the Hodge decomposition, coupled with discrete exterior calculus (DEC), to unveil cell dynamics by examining the decomposed curl-free, divergence-free, and harmonic components of the RNA velocity field in a low dimensional representation, such as a UMAP or a t-SNE representation. Decomposition results show that the decomposed components distinctly reveal key cell dynamic features such as cell cycle, bifurcation, and cell lineage differentiation, regardless of the choice of the low-dimensional representations. The consistency across different representations demonstrates that the Hodge decomposition is a reliable and robust way to extract these cell dynamic features, offering unique analysis and insightful visualization of single-cell RNA velocity fields.

Split Pool Ligation-based Single-cell Transcriptome sequencing (SPLiT-seq) data processing pipeline comparison

BACKGROUND

Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used.

RESULTS

We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required.

CONCLUSION

Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.

Differences in molecular sampling and data processing explain variation among single-cell and single-nucleus RNA-seq experiments

A mechanistic understanding of the biological and technical factors that impact transcript measurements is essential to designing and analyzing single-cell and single-nucleus RNA sequencing experiments. Nuclei contain the same pre-mRNA population as cells, but they contain a small subset of the mRNAs. Nonetheless, early studies argued that single-nucleus analysis yielded results comparable to cellular samples if pre-mRNA measurements were included. However, typical workflows do not distinguish between pre-mRNA and mRNA when estimating gene expression, and variation in their relative abundances across cell types has received limited attention. These gaps are especially important given that incorporating pre-mRNA has become commonplace for both assays, despite known gene length bias in pre-mRNA capture. Here, we reanalyze public data sets from mouse and human to describe the mechanisms and contrasting effects of mRNA and pre-mRNA sampling on gene expression and marker gene selection in single-cell and single-nucleus RNA-seq. We show that pre-mRNA levels vary considerably among cell types, which mediates the degree of gene length bias and limits the generalizability of a recently published normalization method intended to correct for this bias. As an alternative, we repurpose an existing post hoc gene length-based correction method from conventional RNA-seq gene set enrichment analysis. Finally, we show that inclusion of pre-mRNA in bioinformatic processing can impart a larger effect than assay choice itself, which is pivotal to the effective reuse of existing data. These analyses advance our understanding of the sources of variation in single-cell and single-nucleus RNA-seq experiments and provide useful guidance for future studies.

Cationic cholesterol-dependent LNP delivery to lung stem cells, the liver, and heart

Adding a cationic helper lipid to a lipid nanoparticle (LNP) can increase lung delivery and decrease liver delivery. However, it remains unclear whether charge-dependent tropism is universal or, alternatively, whether it depends on the component that is charged. Here, we report evidence that cationic cholesterol-dependent tropism can differ from cationic helper lipid-dependent tropism. By testing how 196 LNPs delivered mRNA to 22 cell types, we found that charged cholesterols led to a different lung:liver delivery ratio than charged helper lipids. We also found that combining cationic cholesterol with a cationic helper lipid led to mRNA delivery in the heart as well as several lung cell types, including stem cell-like populations. These data highlight the utility of exploring charge-dependent LNP tropism.

Representing and extracting knowledge from single-cell data

Single-cell analysis is currently one of the most high-resolution techniques to study biology. The large complex datasets that have been generated have spurred numerous developments in computational biology, in particular the use of advanced statistics and machine learning. This review attempts to explain the deeper theoretical concepts that underpin current state-of-the-art analysis methods. Single-cell analysis is covered from cell, through instruments, to current and upcoming models. The aim of this review is to spread concepts which are not yet in common use, especially from topology and generative processes, and how new statistical models can be developed to capture more of biology. This opens epistemological questions regarding our ontology and models, and some pointers will be given to how natural language processing (NLP) may help overcome our cognitive limitations for understanding single-cell data.

Phylogenomic curation of Ovate Family Proteins (OFPs) in the U's Triangle of Brassica L. indicates stress-induced growth modulation

The Ovate Family Proteins (OFPs) gene family houses a class of proteins that are involved in regulating plant growth and development. To date, there is no report of the simultaneous functional characterization of this gene family in all members of U's Triangle of Brassica. Here, we retrieved a combined total of 256 OFP protein sequences and analyzed their chromosomal localization, gene structure, conserved protein motif domains, and the pattern of cis-acting regulatory elements. The abundance of light-responsive elements like G-box, MRE, and GT1 motif suggests that OFPs are sensitive to the stimuli of light. The protein-protein interaction network analysis revealed that OFP05 and its orthologous genes were involved in regulating the process of transcriptional repression through their interaction with homeodomain transcription factors like KNAT and BLH. The presence of domains like DNA binding 2 and its superfamily speculated the involvement of OFPs in regulating gene expression. The biotic and abiotic stress, and the tissue-specific expression analysis of the RNA-seq datasets revealed that some of the genes such as BjuOFP30, and BnaOFP27, BolOFP11, and BolOFP10 were highly upregulated in seed coat at the mature stage and roots under various chemical stress conditions respectively which suggests their crucial role in plant growth and development processes. Experimental validation of prominent BnaOFPs such as BnaOFP27 confirmed their involvement in regulating gene expression under salinity, heavy metal, drought, heat, and cold stress. The GO and KEGG pathway enrichment analysis also sheds light on the involvement of OFPs in regulating plant growth and development. These findings have the potential to serve as a forerunner for future studies in terms of functionally diverse analysis of the OFP gene family in Brassica and other plant species.

kallisto, bustools, and kb-python for quantifying bulk, single-cell, and single-nucleus RNA-seq

The term "RNA-seq" refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, from single cells, or from single nuclei. The kallisto, bustools, and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples, or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data.

Expression Atlas update: insights from sequencing data at both bulk and single cell level

Expression Atlas (www.ebi.ac.uk/gxa) and its newest counterpart the Single Cell Expression Atlas (www.ebi.ac.uk/gxa/sc) are EMBL-EBI's knowledgebases for gene and protein expression and localisation in bulk and at single cell level. These resources aim to allow users to investigate their expression in normal tissue (baseline) or in response to perturbations such as disease or changes to genotype (differential) across multiple species. Users are invited to search for genes or metadata terms across species or biological conditions in a standardised consistent interface. Alongside these data, new features in Single Cell Expression Atlas allow users to query metadata through our new cell type wheel search. At the experiment level data can be explored through two types of dimensionality reduction plots, t-distributed Stochastic Neighbor Embedding (tSNE) and Uniform Manifold Approximation and Projection (UMAP), overlaid with either clustering or metadata information to assist users' understanding. Data are also visualised as marker gene heatmaps identifying genes that help confer cluster identity. For some data, additional visualisations are available as interactive cell level anatomograms and cell type gene expression heatmaps.

Vulture: cloud-enabled scalable mining of microbial reads in public scRNA-seq data

The rapidly growing collection of public single-cell sequencing data has become a valuable resource for molecular, cellular, and microbial discovery. Previous studies mostly overlooked detecting pathogens in human single-cell sequencing data. Moreover, existing bioinformatics tools lack the scalability to deal with big public data. We introduce Vulture, a scalable cloud-based pipeline that performs microbial calling for single-cell RNA sequencing (scRNA-seq) data, enabling meta-analysis of host-microbial studies from the public domain. In our benchmarking experiments, Vulture is 66% to 88% faster than local tools (PathogenTrack and Venus) and 41% faster than the state-of-the-art cloud-based tool Cumulus, while achieving comparable microbial read identification. In terms of the cost on cloud computing systems, Vulture also shows a cost reduction of 83% ($12 vs. ${\$}$70). We applied Vulture to 2 coronavirus disease 2019, 3 hepatocellular carcinoma (HCC), and 2 gastric cancer human patient cohorts with public sequencing reads data from scRNA-seq experiments and discovered cell type-specific enrichment of severe acute respiratory syndrome coronavirus 2, hepatitis B virus (HBV), and Helicobacter pylori-positive cells, respectively. In the HCC analysis, all cohorts showed hepatocyte-only enrichment of HBV, with cell subtype-associated HBV enrichment based on inferred copy number variations. In summary, Vulture presents a scalable and economical framework to mine unknown host-microbial interactions from large-scale public scRNA-seq data. Vulture is available via an open-source license at https://github.com/holab-hku/Vulture.

Deep generative modeling of transcriptional dynamics for RNA velocity analysis in single cells

RNA velocity has been rapidly adopted to guide interpretation of transcriptional dynamics in snapshot single-cell data; however, current approaches for estimating RNA velocity lack effective strategies for quantifying uncertainty and determining the overall applicability to the system of interest. Here, we present veloVI (velocity variational inference), a deep generative modeling framework for estimating RNA velocity. veloVI learns a gene-specific dynamical model of RNA metabolism and provides a transcriptome-wide quantification of velocity uncertainty. We show that veloVI compares favorably to previous approaches with respect to goodness of fit, consistency across transcriptionally similar cells and stability across preprocessing pipelines for quantifying RNA abundance. Further, we demonstrate that veloVI's posterior velocity uncertainty can be used to assess whether velocity analysis is appropriate for a given dataset. Finally, we highlight veloVI as a flexible framework for modeling transcriptional dynamics by adapting the underlying dynamical model to use time-dependent transcription rates.

Single-cell characterization of human GBM reveals regional differences in tumor-infiltrating leukocyte activation

Glioblastoma (GBM) harbors a highly immunosuppressive tumor microenvironment (TME) which influences glioma growth. Major efforts have been undertaken to describe the TME on a single-cell level. However, human data on regional differences within the TME remain scarce. Here, we performed high-depth single-cell RNA sequencing (scRNAseq) on paired biopsies from the tumor center, peripheral infiltration zone and blood of five primary GBM patients. Through analysis of >45,000 cells, we revealed a regionally distinct transcription profile of microglia (MG) and monocyte-derived macrophages (MdMs) and an impaired activation signature in the tumor-peripheral cytotoxic-cell compartment. Comparing tumor-infiltrating CD8+ T cells with circulating cells identified CX3CR1high and CX3CR1int CD8+ T cells with effector and memory phenotype, respectively, enriched in blood but absent in the TME. Tumor CD8+ T cells displayed a tissue-resident memory phenotype with dysfunctional features. Our analysis provides a regionally resolved mapping of transcriptional states in GBM-associated leukocytes, serving as an additional asset in the effort towards novel therapeutic strategies to combat this fatal disease.

Model of shrimp pond-mediated spatiotemporal dynamic distribution of antibiotic resistance genes in the mangrove habitat of a subtropical gulf

Aquacultures are the main reason for the environmental selection of antibiotic resistance genes (ARGs), resulting in the enrichment of ARGs. As a filter, a marine mangrove ecosystem can reduce antimicrobial resistance (AMR) or eliminate ARGs; however, its elimination mechanism remains unclear. This study investigated the spatiotemporal dynamic distribution of ARGs in two different types of mangrove habitats (shrimp ponds and virgin forests), within a subtropical gulf located in the Beibu Gulf, China, during dry and wet seasons by using metagenomics and real time quantitative polymerase chain reaction (RT-qPCR) analysis. As the key environmental factors, sulfide, salinity, and mobile genetic elements significantly were found to contribute to ARGs distribution, respectively. Wet and dry seasons influenced the dispersal of ARGs but did not affect the microbial community structure. Three potential biomarkers, TEM-116, smeD, and smeE, played key roles in seasonal differences. The key different genes in the biological relevance of absolute abundance were demonstrated by RT-qPCR. Co-occurrence network analysis indicated that high-abundance ARGs were distributed in a modular manner. For the first time, a risk index weighted by risk rank (RIR) was proposed and used to quantify the human risk of ARGs in the mangrove metagenome. The shrimp ponds during the wet season showed the highest RIR detected. In addition to offering a perspective on reducing AMR in mangrove wetlands, this study constructed the first spatiotemporal dynamic model of ARGs in the Beibu Gulf, China and contributed to revealing the global spread of ARGs. Meanwhile, this study proposes a new pipeline for assessing the risk of ARGs, while also exploring the concept of "One Health."

Automatic quality control of single-cell and single-nucleus RNA-seq using valiDrops

Single-cell and single-nucleus RNA-sequencing (sxRNA-seq) measures gene expression in individual cells or nuclei enabling comprehensive characterization of cell types and states. However, isolation of cells or nuclei for sxRNA-seq releases contaminating RNA, which can distort biological signals, through, for example, cell damage and transcript leakage. Thus, identifying barcodes containing high-quality cells or nuclei is a critical analytical step in the processing of sxRNA-seq data. Here, we present valiDrops, an automated method to identify high-quality barcodes and flag dead cells. In valiDrops, barcodes are initially filtered using data-adaptive thresholding on community-standard quality metrics, and subsequently, valiDrops uses a novel clustering-based approach to identify barcodes with distinct biological signals. We benchmark valiDrops and show that biological signals from cell types and states are more distinct, easier to separate and more consistent after filtering by valiDrops compared to existing tools. Finally, we show that valiDrops can predict and flag dead cells with high accuracy. This novel classifier can further improve data quality or be used to identify dead cells to interrogate the biology of cell death. Thus, valiDrops is an effective and easy-to-use method to improve data quality and biological interpretation. Our method is openly available as an R package at www.github.com/madsen-lab/valiDrops.

Abundant binary promoter switches in lineage-determining transcription factors indicate a digital component of cell fate determination

Previous studies of the murine Ly49 and human KIR gene clusters implicated competing sense and antisense promoters in the control of variegated gene expression. In the current study, an examination of transcription factor genes defines an abundance of convergent and divergent sense/antisense promoter pairs, suggesting that competing promoters may control cell fate determination. Differentiation of CD34+ hematopoietic progenitors in vitro shows that cells with GATA1 antisense transcription have enhanced GATA2 transcription and a mast cell phenotype, whereas cells with GATA2 antisense transcription have increased GATA1 transcripts and an erythroblast phenotype. Detailed analyses of the AHR and RORC genes demonstrate the ability of competing promoters to act as binary switches and the association of antisense transcription with an immature/progenitor cell phenotype. These data indicate that alternative cell fates generated by promoter competition in lineage-determining transcription factors contribute to the programming of cell differentiation.

A Global Transcriptional Atlas of the Effect of Sleep Deprivation in the Mouse Frontal Cortex

Sleep deprivation (SD) has negative effects on brain function. Sleep problems are prevalent in neurodevelopmental, neurodegenerative and psychiatric disorders. Thus, understanding the molecular consequences of SD is of fundamental importance in neuroscience. In this study, we present the first simultaneous bulk and single-nuclear (sn)RNA sequencing characterization of the effects of SD in the mouse frontal cortex. We show that SD predominantly affects glutamatergic neurons, specifically in layers 4 and 5, and produces isoform switching of thousands of transcripts. At both the global and cell-type specific level, SD has a large repressive effect on transcription, down-regulating thousands of genes and transcripts; underscoring the importance of accounting for the effects of sleep loss in transcriptome studies of brain function. As a resource we provide extensive characterizations of cell types, genes, transcripts and pathways affected by SD; as well as tutorials for data analysis.

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

Benchmarking single-cell RNA-seq (scRNA-seq) and single-cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) computational tools demands simulators to generate realistic sequencing reads. However, none of the few read simulators aim to mimic real data. To fill this gap, we introduce scReadSim, a single-cell RNA-seq and ATAC-seq read simulator that allows user-specified ground truths and generates synthetic sequencing reads (in a FASTQ or BAM file) by mimicking real data. At both read-sequence and read-count levels, scReadSim mimics real scRNA-seq and scATAC-seq data. Moreover, scReadSim provides ground truths, including unique molecular identifier (UMI) counts for scRNA-seq and open chromatin regions for scATAC-seq. In particular, scReadSim allows users to design cell-type-specific ground-truth open chromatin regions for scATAC-seq data generation. In benchmark applications of scReadSim, we show that UMI-tools achieves the top accuracy in scRNA-seq UMI deduplication, and HMMRATAC and MACS3 achieve the top performance in scATAC-seq peak calling.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

Microfluidics-free single-cell genomics with templated emulsification

Current single-cell RNA-sequencing approaches have limitations that stem from the microfluidic devices or fluid handling steps required for sample processing. We develop a method that does not require specialized microfluidic devices, expertise or hardware. Our approach is based on particle-templated emulsification, which allows single-cell encapsulation and barcoding of cDNA in uniform droplet emulsions with only a vortexer. Particle-templated instant partition sequencing (PIP-seq) accommodates a wide range of emulsification formats, including microwell plates and large-volume conical tubes, enabling thousands of samples or millions of cells to be processed in minutes. We demonstrate that PIP-seq produces high-purity transcriptomes in mouse-human mixing studies, is compatible with multiomics measurements and can accurately characterize cell types in human breast tissue compared to a commercial microfluidic platform. Single-cell transcriptional profiling of mixed phenotype acute leukemia using PIP-seq reveals the emergence of heterogeneity within chemotherapy-resistant cell subsets that were hidden by standard immunophenotyping. PIP-seq is a simple, flexible and scalable next-generation workflow that extends single-cell sequencing to new applications.

Single-cell RNA-seq of maize meiocytes and pollen grains

RNA-sequencing (RNA-seq) provides invaluable knowledge on developmental pathways and the effects of mutant phenotypes. Plant reproductive cells have traditionally been difficult to isolate for genomics because they are rare and often deeply embedded within somatic tissues. Here, we present a protocol to isolate single maize meiocytes and pollen grains for RNA-seq. We discuss how to identify and isolate each sample type under a microscope, prepare RNA-seq libraries and analyze the data. This technique has several advantages over alternative methods, combining the ability to target specific rare cell types while resolving cell-to-cell heterogeneity with single-cell RNA-seq. The technique is compatible with minute amounts of starting material (e.g., a single anther), making it possible to collect dense time courses. Furthermore, developmentally synchronized anthers are saved for microscopy, allowing staging to be performed in parallel with expression analysis. Up to 200 cells can be collected in 4-5 h by someone proficient in tissue dissection, and library preparation can be completed in 2 d by researchers experienced in molecular biology and genomics. This protocol will facilitate research on plant reproduction, providing insights critical to plant breeding, genetics and agriculture.