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.

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.

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.

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.

Discovery of regulatory motifs in 5' untranslated regions using interpretable multi-task learning models.

The sequence in the 5' untranslated regions (UTRs) is known to affect mRNA translation rates. However, the underlying regulatory grammar remains elusive. Here, we propose MTtrans, a multi-task translation rate predictor capable of learning common sequence patterns from datasets across various experimental techniques. The core premise is that common motifs are more likely to be genuinely involved in translation control. MTtrans outperforms existing methods in both accuracy and the ability to capture transferable motifs across species, highlighting its strength in identifying evolutionarily conserved sequence motifs. Our independent fluorescence-activated cell sorting coupled with deep sequencing (FACS-seq) experiment validates the impact of most motifs identified by MTtrans. Additionally, we introduce "GRU-rewiring," a technique to interpret the hidden states of the recurrent units. Gated recurrent unit (GRU)-rewiring allows us to identify regulatory element-enriched positions and examine the local effects of 5' UTR mutations. MTtrans is a powerful tool for deciphering the translation regulatory motifs.

Integrative processing of untargeted metabolomic and lipidomic data using MultiABLER.

Mass spectrometry (MS)-based untargeted metabolomic and lipidomic approaches are being used increasingly in biomedical research. The adoption and integration of these data are critical to the overall multi-omic toolkit. Recently, a sample extraction method called Multi-ABLE has been developed, which enables concurrent generation of proteomic and untargeted metabolomic and lipidomic data from a small amount of tissue. The proteomics field has a well-established set of software for processing of acquired data; however, there is a lack of a unified, off-the-shelf, ready-to-use bioinformatics pipeline that can take advantage of and prepare concurrently generated metabolomic and lipidomic data for joint downstream analyses. Here we present an R pipeline called MultiABLER as a unified and simple upstream processing and analysis pipeline for both metabolomics and lipidomics datasets acquired using liquid chromatography-tandem mass spectrometry. The code is available via an open-source license at https://github.com/holab-hku/MultiABLER.

DCATS: differential composition analysis for flexible single-cell experimental designs.

Differential composition analysis - the identification of cell types that have statistically significant changes in abundance between multiple experimental conditions - is one of the most common tasks in single cell omic data analysis. However, it remains challenging to perform differential composition analysis in the presence of flexible experimental designs and uncertainty in cell type assignment. Here, we introduce a statistical model and an open source R package, DCATS, for differential composition analysis based on a beta-binomial regression framework that addresses these challenges. Our empirical evaluation shows that DCATS consistently maintains high sensitivity and specificity compared to state-of-the-art methods.

Inferring cell diversity in single cell data using consortium-scale epigenetic data as a biological anchor for cell identity.

Methods for cell clustering and gene expression from single-cell RNA sequencing (scRNA-seq) data are essential for biological interpretation of cell processes. Here, we present TRIAGE-Cluster which uses genome-wide epigenetic data from diverse bio-samples to identify genes demarcating cell diversity in scRNA-seq data. By integrating patterns of repressive chromatin deposited across diverse cell types with weighted density estimation, TRIAGE-Cluster determines cell type clusters in a 2D UMAP space. We then present TRIAGE-ParseR, a machine learning method which evaluates gene expression rank lists to define gene groups governing the identity and function of cell types. We demonstrate the utility of this two-step approach using atlases of in vivo and in vitro cell diversification and organogenesis. We also provide a web accessible dashboard for analysis and download of data and software. Collectively, genome-wide epigenetic repression provides a versatile strategy to define cell diversity and study gene regulation of scRNA-seq data.

Altered human gut virome in patients undergoing antibiotics therapy for Helicobacter pylori.

Transient gut microbiota alterations have been reported after antibiotic therapy for Helicobacter pylori. However, alteration in the gut virome after H. pylori eradication remains uncertain. Here, we apply metagenomic sequencing to fecal samples of 44 H. pylori-infected patients at baseline, 6-week (N = 44), and 6-month (N = 33) after treatment. Following H. pylori eradication, we discover contraction of the gut virome diversity, separation of virome community with increased community difference, and shifting towards a higher proportion of core virus. While the gut microbiota is altered at 6-week and restored at 6-month, the virome community shows contraction till 6-month after the treatment with enhanced phage-bacteria interactions at 6-week. Multiple courses of antibiotic treatments further lead to lower virus community diversity when compared with treatment naive patients. Our results demonstrate that H. pylori eradication therapies not only result in transient alteration in gut microbiota but also significantly alter the previously less known gut virome community.

StarmapVis: An interactive and narrative visualisation tool for single-cell and spatial data.

Current single-cell visualisation techniques project high dimensional data into 'map' views to identify high-level structures such as cell clusters and trajectories. New tools are needed to allow the transversal through the high dimensionality of single-cell data to explore the single-cell local neighbourhood. StarmapVis is a convenient web application displaying an interactive downstream analysis of single-cell expression or spatial transcriptomic data. The concise user interface is powered by modern web browsers to explore the variety of viewing angles unavailable to 2D media. Interactive scatter plots display clustering information, while the trajectory and cross-comparison among different coordinates are displayed in connectivity networks. Automated animation of camera view is a unique feature of our tool. StarmapVis also offers a useful animated transition between two-dimensional spatial omic data to three-dimensional single cell coordinates. The usability of StarmapVis is demonstrated by four data sets, showcasing its practical usability. StarmapVis is available at: https://holab-hku.github.io/starmapVis.

Molecular determinants of intrinsic cellular stiffness in health and disease.

In recent years, the role of intrinsic biophysical features, especially cellular stiffness, in diverse cellular and disease processes is being increasingly recognized. New high throughput techniques for the quantification of cellular stiffness facilitate the study of their roles in health and diseases. In this review, we summarized recent discovery about how cellular stiffness is involved in cell stemness, tumorigenesis, and blood diseases. In addition, we review the molecular mechanisms underlying the gene regulation of cellular stiffness in health and disease progression. Finally, we discussed the current understanding on how the cytoskeleton structure and the regulation of these genes contribute to cellular stiffness, highlighting where the field of cellular stiffness is headed.

Incidence of Emergency Department Visits for Sexual Abuse Among Youth in Hong Kong Before and During the COVID-19 Pandemic.

This cohort study assesses the incidence of emergency department (ED) visits in Hong Kong, China, for sexual abuse among youth before and during the COVID-19 pandemic.

Machine learning-coupled combinatorial mutagenesis enables resource-efficient engineering of CRISPR-Cas9 genome editor activities.

The genome-editing Cas9 protein uses multiple amino-acid residues to bind the target DNA. Considering only the residues in proximity to the target DNA as potential sites to optimise Cas9's activity, the number of combinatorial variants to screen through is too massive for a wet-lab experiment. Here we generate and cross-validate ten in silico and experimental datasets of multi-domain combinatorial mutagenesis libraries for Cas9 engineering, and demonstrate that a machine learning-coupled engineering approach reduces the experimental screening burden by as high as 95% while enriching top-performing variants by ∼7.5-fold in comparison to the null model. Using this approach and followed by structure-guided engineering, we identify the N888R/A889Q variant conferring increased editing activity on the protospacer adjacent motif-relaxed KKH variant of Cas9 nuclease from Staphylococcus aureus (KKH-SaCas9) and its derived base editor in human cells. Our work validates a readily applicable workflow to enable resource-efficient high-throughput engineering of genome editor's activity.

MQuad enables clonal substructure discovery using single cell mitochondrial variants.

Mitochondrial mutations are increasingly recognised as informative endogenous genetic markers that can be used to reconstruct cellular clonal structure using single-cell RNA or DNA sequencing data. However, identifying informative mtDNA variants in noisy and sparse single-cell sequencing data is still challenging with few computation methods available. Here we present an open source computational tool MQuad that accurately calls clonally informative mtDNA variants in a population of single cells, and an analysis suite for complete clonality inference, based on single cell RNA, DNA or ATAC sequencing data. Through a variety of simulated and experimental single cell sequencing data, we showed that MQuad can identify mitochondrial variants with both high sensitivity and specificity, outperforming existing methods by a large extent. Furthermore, we demonstrate its wide applicability in different single cell sequencing protocols, particularly in complementing single-nucleotide and copy-number variations to extract finer clonal resolution.

Dynamic changes in antibiotic resistance genes and gut microbiota after Helicobacter pylori eradication therapies.

BACKGROUND: Short-term antibiotics exposure is associated with alterations in microbiota and antibiotic resistance genes (ARGs) in the human gut. While antibiotics are critical in the successful eradication of Helicobacter pylori, the short-term and long-term impacts on the composition and quantity of antibiotics resistance genes after H. pylori eradication are unclear. This study used whole-genome shotgun metagenomic of stool samples to characterize the gut microbiota and ARGs, before and after H. pylori eradication therapy. RESULTS: Forty-four H. pylori-infected patients were recruited, including 21 treatment naïve patients who received clarithromycin-based triple therapy (CLA group) and 23 patients who failed previous therapies, in which 10 received levofloxacin-based quadruple therapy (LEVO group) and 13 received other combinations (OTHER group). Stool samples were collected at baseline (before current treatment), 6 week and 6 month after eradication therapy. At baseline, there was only a slight difference among the three groups on ARGs and gut microbiota. After eradication therapy, there was a transient but significant increase in gut ARGs 6 week post-therapy, among which the LEVO group had the most significant ARGs alteration compared to other two groups. For treatment naïve patients, those with higher ErmF abundance were prone to fail CLA eradication and gain more ARGs after treatment. For gut microbiota, the bacteria richness decreased at 6 week and there was a significant difference in microbiota community among the three groups at 6 week. CONCLUSIONS: Our findings demonstrated the dynamic alterations in gut microbiota and ARGs induced by different eradication therapies, which could influence the choices of antibiotics in eradication therapy.

Automatic flow delay through passive wax valves for paper-based analytical devices.

Microfluidic paper-based analytical devices (µPADs) have been widely explored for point-of-care testing due to their simplicity, low cost, and portability. µPADs with multiple-step reactions usually require precise flow control, especially flow-delay. This paper reports the numerical, mathematical, and experimental studies of flow delay through wax valves surrounded by PDMS walls on paper microfluidics. The predried surfactant in the sample zone diffuses into the liquid sample which can therefore flow through the wax valves. The delay time is automatically regulated by the diffusion of the surfactant after sample loading. The numerical study suggested that both the elevated contact angle and the reduced porosity and pore size in the wax printed region could effectively prevent water but allow liquids with lower contact angles (e.g., surfactant solutions) to flow through. The PDMS walls fabricated using a low-cost liquid dispenser effectively prevented the leakage of surfactant solutions. By controlling the quantity, diffusion distance, and type of the surfactant predried on the chip, the system successfully achieved a delay time ranging from 1.6 to 20 minutes. A mathematical model involving the above parameters was developed based on Fick's second law to predict the delay time. Finally, the flow-delay systems were applied in sequential mixing and distance-based detection of either glucose or alcohol. Linear ranges of 1-100 mg dL-1 and 1-40 mg dL-1 were achieved for glucose and alcohol, respectively. The lower limit detection (LOD) of glucose and alcohol was 1 mg dL-1. The LOD of glucose was only 1/11 of that detected using µPADs without flow control, indicating the advantage of controlling fluid flow. The systematic findings in this study provide critical guidelines for the development and applications of wax valves in automatic flow delay for point-of-care testing.

Generalized and scalable trajectory inference in single-cell omics data with VIA.

Inferring cellular trajectories using a variety of omic data is a critical task in single-cell data science. However, accurate prediction of cell fates, and thereby biologically meaningful discovery, is challenged by the sheer size of single-cell data, the diversity of omic data types, and the complexity of their topologies. We present VIA, a scalable trajectory inference algorithm that overcomes these limitations by using lazy-teleporting random walks to accurately reconstruct complex cellular trajectories beyond tree-like pathways (e.g., cyclic or disconnected structures). We show that VIA robustly and efficiently unravels the fine-grained sub-trajectories in a 1.3-million-cell transcriptomic mouse atlas without losing the global connectivity at such a high cell count. We further apply VIA to discovering elusive lineages and less populous cell fates missed by other methods across a variety of data types, including single-cell proteomic, epigenomic, multi-omics datasets, and a new in-house single-cell morphological dataset.

Genetic screening reveals phospholipid metabolism as a key regulator of the biosynthesis of the redox-active lipid coenzyme Q.

Mitochondrial energy production and function rely on optimal concentrations of the essential redox-active lipid, coenzyme Q (CoQ). CoQ deficiency results in mitochondrial dysfunction associated with increased mitochondrial oxidative stress and a range of pathologies. What drives CoQ deficiency in many of these pathologies is unknown, just as there currently is no effective therapeutic strategy to overcome CoQ deficiency in humans. To date, large-scale studies aimed at systematically interrogating endogenous systems that control CoQ biosynthesis and their potential utility to treat disease have not been carried out. Therefore, we developed a quantitative high-throughput method to determine CoQ concentrations in yeast cells. Applying this method to the Yeast Deletion Collection as a genome-wide screen, 30 genes not known previously to regulate cellular concentrations of CoQ were discovered. In combination with untargeted lipidomics and metabolomics, phosphatidylethanolamine N-methyltransferase (PEMT) deficiency was confirmed as a positive regulator of CoQ synthesis, the first identified to date. Mechanistically, PEMT deficiency alters mitochondrial concentrations of one-carbon metabolites, characterized by an increase in the S-adenosylmethionine to S-adenosylhomocysteine (SAM-to-SAH) ratio that reflects mitochondrial methylation capacity, drives CoQ synthesis, and is associated with a decrease in mitochondrial oxidative stress. The newly described regulatory pathway appears evolutionary conserved, as ablation of PEMT using antisense oligonucleotides increases mitochondrial CoQ in mouse-derived adipocytes that translates to improved glucose utilization by these cells, and protection of mice from high-fat diet-induced insulin resistance. Our studies reveal a previously unrecognized relationship between two spatially distinct lipid pathways with potential implications for the treatment of CoQ deficiencies, mitochondrial oxidative stress/dysfunction, and associated diseases.