Single-Cell and Spatial Transcriptomic Analysis of Human Skin Delineates Intercellular Communication and Pathogenic Cells.

Epidermal homeostasis is governed by a balance between keratinocyte proliferation and differentiation with contributions from cell-cell interactions, but conserved or divergent mechanisms governing this equilibrium across species and how an imbalance contributes to skin disease are largely undefined. To address these questions, human skin single-cell RNA sequencing and spatial transcriptomics data were integrated and compared with mouse skin data. Human skin cell-type annotation was improved using matched spatial transcriptomics data, highlighting the importance of spatial context in cell-type identity, and spatial transcriptomics refined cellular communication inference. In cross-species analyses, we identified a human spinous keratinocyte subpopulation that exhibited proliferative capacity and a heavy metal processing signature, which was absent in mouse and may account for species differences in epidermal thickness. This human subpopulation was expanded in psoriasis and zinc-deficiency dermatitis, attesting to disease relevance and suggesting a paradigm of subpopulation dysfunction as a hallmark of the disease. To assess additional potential subpopulation drivers of skin diseases, we performed cell-of-origin enrichment analysis within genodermatoses, nominating pathogenic cell subpopulations and their communication pathways, which highlighted multiple potential therapeutic targets. This integrated dataset is encompassed in a publicly available web resource to aid mechanistic and translational studies of normal and diseased skin.

A high-resolution flux-matrix model describes the spread of diseases in a spatial network and the effect of mitigation strategies.

Propagation of an epidemic across a spatial network of communities is described by a variant of the SIR model accompanied by an intercommunity infectivity matrix. This matrix is estimated from fluxes between communities, obtained from cell-phone tracking data recorded in the USA between March 2020 and February 2021. We apply this model to the SARS-CoV-2 pandemic by fitting just one global parameter representing the frequency of interaction between individuals. We find that the predicted infections agree reasonably well with the reported cases. We clearly see the effect of "shelter-in-place" policies introduced at the onset of the pandemic. Interestingly, a model with uniform transmission rates produces similar results, suggesting that the epidemic transmission was deeply influenced by air travel. We then study the effect of alternative mitigation policies, in particular restricting long-range travel. We find that this policy is successful in decreasing the epidemic size and slowing down the spread, but less effective than the shelter-in-place policy. This policy can result in a pulled wave of infections. We express its velocity and characterize the shape of the traveling front as a function of the epidemiological parameters. Finally, we discuss a policy of selectively constraining travel based on an edge-betweenness criterion.

Metagenomic prediction of antimicrobial resistance in critically ill patients with lower respiratory tract infections.

BACKGROUND: Antimicrobial resistance (AMR) is rising at an alarming rate and complicating the management of infectious diseases including lower respiratory tract infections (LRTI). Metagenomic next-generation sequencing (mNGS) is a recently established method for culture-independent LRTI diagnosis, but its utility for predicting AMR has remained unclear. We aimed to assess the performance of mNGS for AMR prediction in bacterial LRTI and demonstrate proof of concept for epidemiological AMR surveillance and rapid AMR gene detection using Cas9 enrichment and nanopore sequencing. METHODS: We studied 88 patients with acute respiratory failure between 07/2013 and 9/2018, enrolled through a previous observational study of LRTI. Inclusion criteria were age ≥ 18, need for mechanical ventilation, and respiratory specimen collection within 72 h of intubation. Exclusion criteria were decline of study participation, unclear LRTI status, or no matched RNA and DNA mNGS data from a respiratory specimen. Patients with LRTI were identified by clinical adjudication. mNGS was performed on lower respiratory tract specimens. The primary outcome was mNGS performance for predicting phenotypic antimicrobial susceptibility and was assessed in patients with LRTI from culture-confirmed bacterial pathogens with clinical antimicrobial susceptibility testing (n = 27 patients, n = 32 pathogens). Secondary outcomes included the association between hospital exposure and AMR gene burden in the respiratory microbiome (n = 88 patients), and AMR gene detection using Cas9 targeted enrichment and nanopore sequencing (n = 10 patients). RESULTS: Compared to clinical antimicrobial susceptibility testing, the performance of respiratory mNGS for predicting AMR varied by pathogen, antimicrobial, and nucleic acid type sequenced. For gram-positive bacteria, a combination of RNA + DNA mNGS achieved a sensitivity of 70% (95% confidence interval (CI) 47-87%) and specificity of 95% (CI 85-99%). For gram-negative bacteria, sensitivity was 100% (CI 87-100%) and specificity 64% (CI 48-78%). Patients with hospital-onset LRTI had a greater AMR gene burden in their respiratory microbiome versus those with community-onset LRTI (p = 0.00030), or those without LRTI (p = 0.0024). We found that Cas9 targeted sequencing could enrich for low abundance AMR genes by > 2500-fold and enabled their rapid detection using a nanopore platform. CONCLUSIONS: mNGS has utility for the detection and surveillance of resistant bacterial LRTI pathogens.

Molecular hallmarks of heterochronic parabiosis at single-cell resolution.

The ability to slow or reverse biological ageing would have major implications for mitigating disease risk and maintaining vitality1. Although an increasing number of interventions show promise for rejuvenation2, their effectiveness on disparate cell types across the body and the molecular pathways susceptible to rejuvenation remain largely unexplored. Here we performed single-cell RNA sequencing on 20 organs to reveal cell-type-specific responses to young and aged blood in heterochronic parabiosis. Adipose mesenchymal stromal cells, haematopoietic stem cells and hepatocytes are among those cell types that are especially responsive. On the pathway level, young blood invokes new gene sets in addition to reversing established ageing patterns, with the global rescue of genes encoding electron transport chain subunits pinpointing a prominent role of mitochondrial function in parabiosis-mediated rejuvenation. We observed an almost universal loss of gene expression with age that is largely mimicked by parabiosis: aged blood reduces global gene expression, and young blood restores it in select cell types. Together, these data lay the groundwork for a systemic understanding of the interplay between blood-borne factors and cellular integrity.

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. RESULTS: A total of 12 124 tests were performed yielding 1099 positives. From these, 928 high-quality genomes were generated. Certain viral lineages bearing spike mutations, defined in part by L452R, S13I, and W152C, comprised 54.4% of the total sequences from January, compared to 15.7% in November. Household contacts exposed to the "California" or "West Coast" variants (B.1.427 and B.1.429) were at higher risk of infection compared to household contacts exposed to lineages lacking these variants (0.36 vs 0.29, risk ratio [RR] = 1.28; 95% confidence interval [CI]: 1.00-1.64). The reproductive number was estimated to be modestly higher than other lineages spreading in California during the second half of 2020. Viral loads were similar among persons infected with West Coast versus non-West Coast strains, as was the proportion of individuals with symptoms (60.9% vs 64.3%). CONCLUSIONS: The increase in prevalence, relative household attack rates, and reproductive number are consistent with a modest transmissibility increase of the West Coast variants. Summary: We observed a growing prevalence and modestly elevated attack rate for "West Coast" severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in a community testing setting in San Francisco during January 2021, suggesting its modestly higher transmissibility.

Single-Cell Analysis for Whole-Organism Datasets.

Cell atlases are essential companions to the genome as they elucidate how genes are used in a cell type-specific manner or how the usage of genes changes over the lifetime of an organism. This review explores recent advances in whole-organism single-cell atlases, which enable understanding of cell heterogeneity and tissue and cell fate, both in health and disease. Here we provide an overview of recent efforts to build cell atlases across species and discuss the challenges that the field is currently facing. Moreover, we propose the concept of having a knowledgebase that can scale with the number of experiments and computational approaches and a new feedback loop for development and benchmarking of computational methods that includes contributions from the users. These two aspects are key for community efforts in single-cell biology that will help produce a comprehensive annotated map of cell types and states with unparalleled resolution.

A random-walk-based epidemiological model.

Random walkers on a two-dimensional square lattice are used to explore the spatio-temporal growth of an epidemic. We have found that a simple random-walk system generates non-trivial dynamics compared with traditional well-mixed models. Phase diagrams characterizing the long-term behaviors of the epidemics are calculated numerically. The functional dependence of the basic reproductive number [Formula: see text] on the model's defining parameters reveals the role of spatial fluctuations and leads to a novel expression for [Formula: see text]. Special attention is given to simulations of inter-regional transmission of the contagion. The scaling of the epidemic with respect to space and time scales is studied in detail in the critical region, which is shown to be compatible with the directed-percolation universality class.

Leveraging the Cell Ontology to classify unseen cell types.

Single cell technologies are rapidly generating large amounts of data that enables us to understand biological systems at single-cell resolution. However, joint analysis of datasets generated by independent labs remains challenging due to a lack of consistent terminology to describe cell types. Here, we present OnClass, an algorithm and accompanying software for automatically classifying cells into cell types that are part of the controlled vocabulary that forms the Cell Ontology. A key advantage of OnClass is its capability to classify cells into cell types not present in the training data because it uses the Cell Ontology graph to infer cell type relationships. Furthermore, OnClass can be used to identify marker genes for all the cell ontology categories, regardless of whether the cell types are present or absent in the training data, suggesting that OnClass goes beyond a simple annotation tool for single cell datasets, being the first algorithm capable to identify marker genes specific to all terms of the Cell Ontology and offering the possibility of refining the Cell Ontology using a data-centric approach.

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. RESULTS: A total of 12,124 tests were performed yielding 1,099 positives. From these, 811 high quality genomes were generated. Certain viral lineages bearing spike mutations, defined in part by L452R, S13I, and W152C, comprised 54.9% of the total sequences from January, compared to 15.7% in November. Household contacts exposed to "West Coast" variants were at higher risk of infection compared to household contacts exposed to lineages lacking these variants (0.357 vs 0.294, RR=1.29; 95% CI:1.01-1.64). The reproductive number was estimated to be modestly higher than other lineages spreading in California during the second half of 2020. Viral loads were similar among persons infected with West Coast versus non-West Coast strains, as was the proportion of individuals with symptoms (60.9% vs 64.1%). CONCLUSIONS: The increase in prevalence, relative household attack rates, and reproductive number are consistent with a modest transmissibility increase of the West Coast variants; however, additional laboratory and epidemiological studies are required to better understand differences between these variants. SUMMARY: We observed a growing prevalence and elevated attack rate for "West Coast" SARS-CoV-2 variants in a community testing setting in San Francisco during January 2021, suggesting its modestly higher transmissibility.

A minimal model for household-based testing and tracing in epidemics.

In a previous work (Huberet al.2020Phys. Biol.17065010), we discussed virus transmission dynamics modified by a uniform clustering of contacts in the population: close contacts within households and more distant contacts between households. In this paper, we discuss testing and tracing in such a stratified population. We propose a minimal tracing strategy consisting of random testing of the entire population plus full testing of the households of those persons found positive. We provide estimates of testing frequency for this strategy to work.

Ageing hallmarks exhibit organ-specific temporal signatures.

Ageing is the single greatest cause of disease and death worldwide, and understanding the associated processes could vastly improve quality of life. Although major categories of ageing damage have been identified-such as altered intercellular communication, loss of proteostasis and eroded mitochondrial function1-these deleterious processes interact with extraordinary complexity within and between organs, and a comprehensive, whole-organism analysis of ageing dynamics has been lacking. Here we performed bulk RNA sequencing of 17 organs and plasma proteomics at 10 ages across the lifespan of Mus musculus, and integrated these findings with data from the accompanying Tabula Muris Senis2-or 'Mouse Ageing Cell Atlas'-which follows on from the original Tabula Muris3. We reveal linear and nonlinear shifts in gene expression during ageing, with the associated genes clustered in consistent trajectory groups with coherent biological functions-including extracellular matrix regulation, unfolded protein binding, mitochondrial function, and inflammatory and immune response. Notably, these gene sets show similar expression across tissues, differing only in the amplitude and the age of onset of expression. Widespread activation of immune cells is especially pronounced, and is first detectable in white adipose depots during middle age. Single-cell RNA sequencing confirms the accumulation of T cells and B cells in adipose tissue-including plasma cells that express immunoglobulin J-which also accrue concurrently across diverse organs. Finally, we show how gene expression shifts in distinct tissues are highly correlated with corresponding protein levels in plasma, thus potentially contributing to the ageing of the systemic circulation. Together, these data demonstrate a similar yet asynchronous inter- and intra-organ progression of ageing, providing a foundation from which to track systemic sources of declining health at old age.

Functional Enrichment and Analysis of Antigen-Specific Memory B Cell Antibody Repertoires in PBMCs.

Phenotypic screening of antigen-specific antibodies in human blood is a common diagnostic test for infectious agents and a correlate of protection after vaccination. In addition to long-lived antibody secreting plasma cells residing in the bone marrow, memory B cells are a latent source of antigen-experienced, long-term immunity that can be found at low frequencies in circulating peripheral blood mononuclear cells (PBMCs). Assessing the genotype, clonal frequency, quality, and function of antibodies resulting from an individual's persistent memory B cell repertoire can help inform the success or failure of immune protection. Using in vitro polyclonal stimulation, we functionally expand the memory repertoire from PBMCs and clonally map monoclonal antibodies from this population. We show that combining deep sequencing of stimulated memory B cell repertoires with retrieving single antigen-specific cells is a promising approach in evaluating the latent, functional B cell memory in PBMCs.

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. FLASH-NGS achieves up to 5 orders of magnitude of enrichment and sub-attomolar gene detection with minimal background. We provide an open-source software tool (FLASHit) for guide RNA design. Here we applied it to detection of antimicrobial resistance genes in respiratory fluid and dried blood spots, but FLASH-NGS is applicable to all areas that rely on multiplex PCR.