scDecouple: decoupling cellular response from infected proportion bias in scCRISPR-seq.

Single-cell clustered regularly interspaced short palindromic repeats-sequencing (scCRISPR-seq) is an emerging high-throughput CRISPR screening technology where the true cellular response to perturbation is coupled with infected proportion bias of guide RNAs (gRNAs) across different cell clusters. The mixing of these effects introduces noise into scCRISPR-seq data analysis and thus obstacles to relevant studies. We developed scDecouple to decouple true cellular response of perturbation from the influence of infected proportion bias. scDecouple first models the distribution of gene expression profiles in perturbed cells and then iteratively finds the maximum likelihood of cell cluster proportions as well as the cellular response for each gRNA. We demonstrated its performance in a series of simulation experiments. By applying scDecouple to real scCRISPR-seq data, we found that scDecouple enhances the identification of biologically perturbation-related genes. scDecouple can benefit scCRISPR-seq data analysis, especially in the case of heterogeneous samples or complex gRNA libraries.

MetaQuad: shared informative variants discovery in metagenomic samples.

Motivation: Strain-level analysis of metagenomic data has garnered significant interest in recent years. Microbial single nucleotide polymorphisms (SNPs) are genomic variants that can reflect strain-level differences within a microbial species. The diversity and emergence of SNPs in microbial genomes may reveal evolutionary history and environmental adaptation in microbial populations. However, efficient discovery of shared polymorphic variants in a large collection metagenomic samples remains a computational challenge. Results: MetaQuad utilizes a density-based clustering technique to effectively distinguish between shared variants and non-polymorphic sites using shotgun metagenomic data. Empirical comparisons with other state-of-the-art methods show that MetaQuad significantly reduces the number of false positive SNPs without greatly affecting the true positive rate. We used MetaQuad to identify antibiotic-associated variants in patients who underwent Helicobacter pylori eradication therapy. MetaQuad detected 7591 variants across 529 antibiotic resistance genes. The nucleotide diversity of some genes is increased 6 weeks after antibiotic treatment, potentially indicating the role of these genes in specific antibiotic treatments. Availability and implementation: MetaQuad is an open-source Python package available via https://github.com/holab-hku/MetaQuad.

StaVia: spatially and temporally aware cartography with higher-order random walks for cell atlases.

Single-cell atlases pose daunting computational challenges pertaining to the integration of spatial and temporal information and the visualization of trajectories across large atlases. We introduce StaVia, a computational framework that synergizes multi-faceted single-cell data with higher-order random walks that leverage the memory of cells' past states, fused with a cartographic Atlas View that offers intuitive graph visualization. This spatially aware cartography captures relationships between cell populations based on their spatial location as well as their gene expression and developmental stage. We demonstrate this using zebrafish gastrulation data, underscoring its potential to dissect complex biological landscapes in both spatial and temporal contexts.

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.

The Effect of the COVID-19 Pandemic on Non-COVID-19 Deaths: Population-Wide Retrospective Cohort Study.

BACKGROUND: Health care avoidance in the COVID-19 pandemic has been widely reported. Yet few studies have investigated the dynamics of hospital avoidance behavior during pandemic waves and inferred its impact on excess non-COVID-19 deaths. OBJECTIVE: This study aimed to measure the impact of hospital avoidance on excess non-COVID-19 deaths in public hospitals in Hong Kong. METHODS: This was a retrospective cohort study involving 11,966,786 patients examined between January 1, 2016, and December 31, 2021, in Hong Kong. All data were linked to service, treatment, and outcomes. To estimate excess mortality, the 2-stage least squares method was used with daily tallies of emergency department (ED) visits and 28-day mortality. Records for older people were categorized by long-term care (LTC) home status, and comorbidities were used to explain the demographic and clinical attributes of excess 28-day mortality. The primary outcome was actual excess death in 2020 and 2021. The 2-stage least squares method was used to estimate the daily excess 28-day mortality by daily reduced visits. RESULTS: Compared with the prepandemic (2016-2019) average, there was a reduction in total ED visits in 2020 of 25.4% (548,116/2,142,609). During the same period, the 28-day mortality of non-COVID-19 ED deaths increased by 7.82% (2689/34,370) compared with 2016-2019. The actual excess deaths in 2020 and 2021 were 3143 and 4013, respectively. The estimated total excess non-COVID-19 28-day deaths among older people in 2020 to 2021 were 1958 (95% CI 1100-2820; no time lag). Deaths on arrival (DOAs) or deaths before arrival (DBAs) increased by 33.6% (1457/4336) in 2020, while non-DOA/DBAs increased only by a moderate 4.97% (1202/24,204). In both types of deaths, the increases were higher during wave periods than in nonwave periods. Moreover, non-LTC patients saw a greater reduction in ED visits than LTC patients across all waves, by more than 10% (non-LTC: 93,896/363,879, 25.8%; LTC: 7,956/67,090, 11.9%). Most of the comorbidity subsets demonstrated an annualized reduction in visits in 2020. Renal diseases and severe liver diseases saw notable increases in deaths. CONCLUSIONS: We demonstrated a statistical method to estimate hospital avoidance behavior during a pandemic and quantified the consequent excess 28-day mortality with a focus on older people, who had high frequencies of ED visits and deaths. This study serves as an informed alert and possible investigational guideline for health care professionals for hospital avoidance behavior and its consequences.

In vitro and in vivo models define a molecular signature reference for human embryonic notochordal cells.

Understanding the emergence of human notochordal cells (NC) is essential for the development of regenerative approaches. We present a comprehensive investigation into the specification and generation of bona fide NC using a straightforward pluripotent stem cell (PSC)-based system benchmarked with human fetal notochord. By integrating in vitro and in vivo transcriptomic data at single-cell resolution, we establish an extended molecular signature and overcome the limitations associated with studying human notochordal lineage at early developmental stages. We show that TGF-ß inhibition enhances the yield and homogeneity of notochordal lineage commitment in vitro. Furthermore, this study characterizes regulators of cell-fate decision and matrisome enriched in the notochordal niche. Importantly, we identify specific cell-surface markers opening avenues for differentiation refinement, NC purification, and functional studies. Altogether, this study provides a human notochord transcriptomic reference that will serve as a resource for notochord identification in human systems, diseased-tissues modeling, and facilitating future biomedical research.