Single-cell analysis of human fibrous dysplasia bone reveals a fibrotic transcriptome and GNAS variants in endothelial, perivascular, and stromal cells.

Genetic mosaicism is a leading cause of human disease across the lifespan. Improving the tools to detect somatic mosaicism and applying them to understand the cellular and molecular mechanisms that contribute to disease is of critical importance for improving human health. Fibrous dysplasia (FD) is a prototypical disease of G-GPCR activation caused by somatic, mosaic GNAS variants (c.

Identification and regulation of circulating tumor-TCR-matched cytotoxic CD4+ lymphocytes by KLRG1 in bladder cancer.

While cytotoxic CD4+ tumor-infiltrating lymphocytes have anticancer activity in patients, whether these can be noninvasively monitored and how these are regulated remains obscure. By matching single cells with T cell receptors (TCRs) in tumor and blood of patients with bladder cancer, we identified distinct pools of tumor-matching cytotoxic CD4+ T cells in the periphery directly reflecting the predominant antigenic specificities of intratumoral CD4+ tumor-infiltrating lymphocytes. On one hand, the granzyme B-expressing (GZMB-expressing) cytotoxic CD4+ subset proliferated in blood in response to PD-1 blockade but was separately regulated by the killer cell lectin-like receptor G1 (KLRG1), which inhibited their killing by interacting with E-cadherin.

Reducing batch effects in single cell chromatin accessibility measurements by pooled transposition with MULTI-ATAC.

Large-scale scATAC-seq experiments are challenging because of their costs, lengthy protocols, and confounding batch effects. Several sample multiplexing technologies aim to address these challenges, but do not remove batch effects introduced when performing transposition reactions in parallel. We demonstrate that sample-to-sample variability in nuclei-to-Tn5 ratios is a major cause of batch effects and develop MULTI-ATAC, a multiplexing method that pools samples prior to transposition, as a solution.

Immune Modulation with RANKL Blockade through Denosumab Treatment in Patients with Cancer.

Denosumab is a fully human mAb that binds receptor activator of NFκB ligand (RANKL). It is routinely administered to patients with cancer to reduce the incidence of new bone metastasis. RANK-RANKL interactions regulate bone turnover by controlling osteoclast recruitment, development, and activity.

Compact CRISPR genetic screens enabled by improved guide RNA library cloning.

CRISPR genome editing approaches theoretically enable researchers to define the function of each human gene in specific cell types, but challenges remain to efficiently perform genetic perturbations in relevant models. In this work, we develop a library cloning protocol that increases sgRNA uniformity and greatly reduces bias in existing genome-wide libraries. We demonstrate that our libraries can achieve equivalent or better statistical power compared to previously reported screens using an order of magnitude fewer cells.

D-SPIN constructs gene regulatory network models from multiplexed scRNA-seq data revealing organizing principles of cellular perturbation response.

Gene regulatory networks within cells modulate the expression of the genome in response to signals and changing environmental conditions. Reconstructions of gene regulatory networks can reveal the information processing and control principles used by cells to maintain homeostasis and execute cell-state transitions. Here, we introduce a computational framework, D-SPIN, that generates quantitative models of gene regulatory networks from single-cell mRNA-seq datasets collected across thousands of distinct perturbation conditions.

Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution.

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution.

SCITO-seq: single-cell combinatorial indexed cytometry sequencing.

The development of DNA-barcoded antibodies to tag cell surface molecules has enabled the use of droplet-based single-cell sequencing (dsc-seq) to profile protein abundances from thousands of cells simultaneously. As compared to flow and mass cytometry, the high per cell cost of current dsc-seq-based workflows precludes their use in clinical applications and large-scale pooled screens. Here, we introduce SCITO-seq, a workflow that uses splint oligonucleotides (oligos) to enable combinatorially indexed dsc-seq of DNA-barcoded antibodies from over 10 cells per reaction using commercial microfluidics.

XYZeq: Spatially resolved single-cell RNA sequencing reveals expression heterogeneity in the tumor microenvironment.

Single-cell RNA sequencing (scRNA-seq) of tissues has revealed remarkable heterogeneity of cell types and states but does not provide information on the spatial organization of cells. To better understand how individual cells function within an anatomical space, we developed XYZeq, a workflow that encodes spatial metadata into scRNA-seq libraries. We used XYZeq to profile mouse tumor models to capture spatially barcoded transcriptomes from tens of thousands of cells.

Estimation of Secondary Household Attack Rates for Emergent Spike L452R Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants Detected by Genomic Surveillance at a Community-Based Testing Site in San Francisco.

Background: Sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral genome from patient samples is an important epidemiological tool for monitoring and responding to the pandemic, including the emergence of new mutations in specific communities.

Methods: SARS-CoV-2 genomic sequences were generated from positive samples collected, along with epidemiological metadata, at a walk-up, rapid testing site in the Mission District of San Francisco, California during 22 November to 1 December, 2020, and 10-29 January 2021. Secondary household attack rates and mean sample viral load were estimated and compared across observed variants.

Estimation of secondary household attack rates for emergent SARS-CoV-2 variants detected by genomic surveillance at a community-based testing site in San Francisco.

Background: Sequencing of the SARS-CoV-2 viral genome from patient samples is an important epidemiological tool for monitoring and responding to the pandemic, including the emergence of new mutations in specific communities.

Methods: SARS-CoV-2 genomic sequences were generated from positive samples collected, along with epidemiological metadata, at a walk-up, rapid testing site in the Mission District of San Francisco, California during November 22-December 2, 2020 and January 10-29, 2021. Secondary household attack rates and mean sample viral load were estimated and compared across observed variants.

Macrolide and Nonmacrolide Resistance with Mass Azithromycin Distribution.

Background: Mass distribution of azithromycin to preschool children twice yearly for 2 years has been shown to reduce childhood mortality in sub-Saharan Africa but at the cost of amplifying macrolide resistance. The effects on the gut resistome, a reservoir of antimicrobial resistance genes in the body, of twice-yearly administration of azithromycin for a longer period are unclear.

Methods: We investigated the gut resistome of children after they received twice-yearly distributions of azithromycin for 4 years.

Rapid pathogen detection by metagenomic next-generation sequencing of infected body fluids.

We developed a metagenomic next-generation sequencing (mNGS) test using cell-free DNA from body fluids to identify pathogens. The performance of mNGS testing of 182 body fluids from 160 patients with acute illness was evaluated using two sequencing platforms in comparison to microbiological testing using culture, 16S bacterial PCR and/or 28S-internal transcribed ribosomal gene spacer (28S-ITS) fungal PCR. Test sensitivity and specificity of detection were 79 and 91% for bacteria and 91 and 89% for fungi, respectively, by Illumina sequencing; and 75 and 81% for bacteria and 91 and 100% for fungi, respectively, by nanopore sequencing.

Robust Sequence Determinants of α-Synuclein Toxicity in Yeast Implicate Membrane Binding.

Protein conformations are shaped by cellular environments, but how environmental changes alter the conformational landscapes of specific proteins remains largely uncharacterized, in part due to the challenge of probing protein structures in living cells. Here, we use deep mutational scanning to investigate how a toxic conformation of α-synuclein, a dynamic protein linked to Parkinson's disease, responds to perturbations of cellular proteostasis. In the context of a course for graduate students in the UCSF Integrative Program in Quantitative Biology, we screened a comprehensive library of α-synuclein missense mutants in yeast cells treated with a variety of small molecules that perturb cellular processes linked to α-synuclein biology and pathobiology.

The DNA methylation landscape of advanced prostate cancer.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes that were detectable only with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hypermethylation and somatic mutations in TET2, DNMT3B, IDH1 and BRAF.

ZipSeq: barcoding for real-time mapping of single cell transcriptomes.

Spatial transcriptomics seeks to integrate single cell transcriptomic data within the three-dimensional space of multicellular biology. Current methods to correlate a cell's position with its transcriptome in living tissues have various limitations. We developed an approach, called 'ZipSeq', that uses patterned illumination and photocaged oligonucleotides to serially print barcodes ('zipcodes') onto live cells in intact tissues, in real time and with an on-the-fly selection of patterns.

Deep mutational scanning reveals the structural basis for α-synuclein activity.

Defining the biologically active structures of proteins in their cellular environments remains challenging for proteins with multiple conformations and functions, where only a minor conformer might be associated with a given function. Here, we use deep mutational scanning to probe the structure and dynamics of α-synuclein, a protein known to adopt disordered, helical and amyloid conformations. We examined the effects of 2,600 single-residue substitutions on the ability of intracellularly expressed α-synuclein to slow the growth of yeast.

MULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices.

Sample multiplexing facilitates scRNA-seq by reducing costs and identifying artifacts such as cell doublets. However, universal and scalable sample barcoding strategies have not been described. We therefore developed MULTI-seq: multiplexing using lipid-tagged indices for single-cell and single-nucleus RNA sequencing.

FLASH: a next-generation CRISPR diagnostic for multiplexed detection of antimicrobial resistance sequences.

The growing prevalence of deadly microbes with resistance to previously life-saving drug therapies is a dire threat to human health. Detection of low abundance pathogen sequences remains a challenge for metagenomic Next Generation Sequencing (NGS). We introduce FLASH (Finding Low Abundance Sequences by Hybridization), a next-generation CRISPR/Cas9 diagnostic method that takes advantage of the efficiency, specificity and flexibility of Cas9 to enrich for a programmed set of sequences.

Predicting Splicing from Primary Sequence with Deep Learning.

The splicing of pre-mRNAs into mature transcripts is remarkable for its precision, but the mechanisms by which the cellular machinery achieves such specificity are incompletely understood. Here, we describe a deep neural network that accurately predicts splice junctions from an arbitrary pre-mRNA transcript sequence, enabling precise prediction of noncoding genetic variants that cause cryptic splicing. Synonymous and intronic mutations with predicted splice-altering consequence validate at a high rate on RNA-seq and are strongly deleterious in the human population.

Miniaturization and optimization of 384-well compatible RNA sequencing library preparation.

Preparation of high-quality sequencing libraries is a costly and time-consuming component of metagenomic next generation sequencing (mNGS). While the overall cost of sequencing has dropped significantly over recent years, the reagents needed to prepare sequencing samples are likely to become the dominant expense in the process. Furthermore, libraries prepared by hand are subject to human variability and needless waste due to limitations of manual pipetting volumes.

Pulmonary Metagenomic Sequencing Suggests Missed Infections in Immunocompromised Children.

Background: Despite improved diagnostics, pulmonary pathogens in immunocompromised children frequently evade detection, leading to significant mortality. Therefore, we aimed to develop a highly sensitive metagenomic next-generation sequencing (mNGS) assay capable of evaluating the pulmonary microbiome and identifying diverse pathogens in the lungs of immunocompromised children.

Methods: We collected 41 lower respiratory specimens from 34 immunocompromised children undergoing evaluation for pulmonary disease at 3 children's hospitals from 2014-2016.

Mass Azithromycin Distribution and Community Microbiome: A Cluster-Randomized Trial.

Background: Mass distributions of oral azithromycin have long been used to eliminate trachoma, and they are now being proposed to reduce childhood mortality. The observed benefit appears to be augmented with each additional treatment, suggesting a possible community-level effect. Here, we assess whether 2 biannual mass treatments of preschool children affect the community's gut microbiome at 6 months after the last distribution.