Spatial proteogenomics reveals distinct and evolutionarily conserved hepatic macrophage niches.

The liver is the largest solid organ in the body, yet it remains incompletely characterized. Here we present a spatial proteogenomic atlas of the healthy and obese human and murine liver combining single-cell CITE-seq, single-nuclei sequencing, spatial transcriptomics, and spatial proteomics. By integrating these multi-omic datasets, we provide validated strategies to reliably discriminate and localize all hepatic cells, including a population of lipid-associated macrophages (LAMs) at the bile ducts. We then align this atlas across seven species, revealing the conserved program of bona fide Kupffer cells and LAMs. We also uncover the respective spatially resolved cellular niches of these macrophages and the microenvironmental circuits driving their unique transcriptomic identities. We demonstrate that LAMs are induced by local lipid exposure, leading to their induction in steatotic regions of the murine and human liver, while Kupffer cell development crucially depends on their cross-talk with hepatic stellate cells via the evolutionarily conserved ALK1-BMP9/10 axis.

The Genetics of Transcription Factor DNA Binding Variation.

Most complex trait-associated variants are located in non-coding regulatory regions of the genome, where they have been shown to disrupt transcription factor (TF)-DNA binding motifs. Variable TF-DNA interactions are therefore increasingly considered as key drivers of phenotypic variation. However, recent genome-wide studies revealed that the majority of variable TF-DNA binding events are not driven by sequence alterations in the motif of the studied TF. This observation implies that the molecular mechanisms underlying TF-DNA binding variation and, by extrapolation, inter-individual phenotypic variation are more complex than originally anticipated. Here, we summarize the findings that led to this important paradigm shift and review proposed mechanisms for local, proximal, or distal genetic variation-driven variable TF-DNA binding. In addition, we discuss the biomedical implications of these findings for our ability to dissect the molecular role(s) of non-coding genetic variants in complex traits, including disease susceptibility.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

Live-seq enables temporal transcriptomic recording of single cells.

Single-cell transcriptomics (scRNA-seq) has greatly advanced our ability to characterize cellular heterogeneity1. However, scRNA-seq requires lysing cells, which impedes further molecular or functional analyses on the same cells. Here, we established Live-seq, a single-cell transcriptome profiling approach that preserves cell viability during RNA extraction using fluidic force microscopy2,3, thus allowing to couple a cell's ground-state transcriptome to its downstream molecular or phenotypic behaviour. To benchmark Live-seq, we used cell growth, functional responses and whole-cell transcriptome read-outs to demonstrate that Live-seq can accurately stratify diverse cell types and states without inducing major cellular perturbations. As a proof of concept, we show that Live-seq can be used to directly map a cell's trajectory by sequentially profiling the transcriptomes of individual macrophages before and after lipopolysaccharide (LPS) stimulation, and of adipose stromal cells pre- and post-differentiation. In addition, we demonstrate that Live-seq can function as a transcriptomic recorder by preregistering the transcriptomes of individual macrophages that were subsequently monitored by time-lapse imaging after LPS exposure. This enabled the unsupervised, genome-wide ranking of genes on the basis of their ability to affect macrophage LPS response heterogeneity, revealing basal Nfkbia expression level and cell cycle state as important phenotypic determinants, which we experimentally validated. Thus, Live-seq can address a broad range of biological questions by transforming scRNA-seq from an end-point to a temporal analysis approach.

Chromatin modules and their implication in genomic organization and gene regulation.

Regulation of gene expression is a complex but highly guided process. While genomic technologies and computational approaches have allowed high-throughput mapping of cis-regulatory elements (CREs) and their interactions in 3D, their precise role in regulating gene expression remains obscure. Recent complementary observations revealed that interactions between CREs frequently result in the formation of small-scale functional modules within topologically associating domains. Such chromatin modules likely emerge from a complex interplay between regulatory machineries assembled at CREs, including site-specific binding of transcription factors. Here, we review the methods that allow identifying chromatin modules, summarize possible mechanisms that steer CRE interactions within these modules, and discuss outstanding challenges to uncover how chromatin modules fit in our current understanding of the functional 3D genome.

Mammalian adipogenesis regulator (Areg) cells use retinoic acid signalling to be non- and anti-adipogenic in age-dependent manner.

Adipose stem and precursor cells (ASPCs) give rise to adipocytes and determine the composition and plasticity of adipose tissue. Recently, several studies have demonstrated that ASPCs partition into at least three distinct cell subpopulations, including the enigmatic CD142+ cells. An outstanding challenge is to functionally characterise this population, as discrepant properties, from adipogenic to non- and anti-adipogenic, have been reported for these cells. To resolve these phenotypic ambiguities, we characterised mammalian subcutaneous CD142+ ASPCs across various experimental conditions, demonstrating that CD142+ ASPCs exhibit high molecular and phenotypic robustness. Specifically, we find these cells to be firmly non- and anti-adipogenic both in vitro and in vivo, with their inhibitory signals also impacting adipogenic human cells. However, these CD142+ ASPC-specific properties exhibit surprising temporal phenotypic alterations, and emerge only in an age-dependent manner. Finally, using multi-omic and functional assays, we show that the inhibitory nature of these adipogenesis-regulatory CD142+ ASPCs (Aregs) is driven by specifically expressed secretory factors that cooperate with the retinoic acid signalling pathway to transform the adipogenic state of CD142- ASPCs into a non-adipogenic, Areg-like state.

Pseudomonas aeruginosa type IV pili actively induce mucus contraction to form biofilms in tissue-engineered human airways.

The opportunistic pathogen Pseudomonas aeruginosa causes antibiotic-recalcitrant pneumonia by forming biofilms in the respiratory tract. Despite extensive in vitro experimentation, how P. aeruginosa forms biofilms at the airway mucosa is unresolved. To investigate the process of biofilm formation in realistic conditions, we developed AirGels: 3D, optically accessible tissue-engineered human lung models that emulate the airway mucosal environment. AirGels recapitulate important factors that mediate host-pathogen interactions including mucus secretion, flow and air-liquid interface (ALI), while accommodating high-resolution live microscopy. With AirGels, we investigated the contributions of mucus to P. aeruginosa biofilm biogenesis in in vivo-like conditions. We found that P. aeruginosa forms mucus-associated biofilms within hours by contracting luminal mucus early during colonization. Mucus contractions facilitate aggregation, thereby nucleating biofilms. We show that P. aeruginosa actively contracts mucus using retractile filaments called type IV pili. Our results therefore suggest that, while protecting epithelia, mucus constitutes a breeding ground for biofilms.

Deterministic scRNA-seq captures variation in intestinal crypt and organoid composition.

Single-cell RNA sequencing (scRNA-seq) approaches have transformed our ability to resolve cellular properties across systems, but are currently tailored toward large cell inputs (>1,000 cells). This renders them inefficient and costly when processing small, individual tissue samples, a problem that tends to be resolved by loading bulk samples, yielding confounded mosaic cell population read-outs. Here, we developed a deterministic, mRNA-capture bead and cell co-encapsulation dropleting system, DisCo, aimed at processing low-input samples (<500 cells). We demonstrate that DisCo enables precise particle and cell positioning and droplet sorting control through combined machine-vision and multilayer microfluidics, enabling continuous processing of low-input single-cell suspensions at high capture efficiency (>70%) and at speeds up to 350 cells per hour. To underscore DisCo's unique capabilities, we analyzed 31 individual intestinal organoids at varying developmental stages. This revealed extensive organoid heterogeneity, identifying distinct subtypes including a regenerative fetal-like Ly6a+ stem cell population that persists as symmetrical cysts, or spheroids, even under differentiation conditions, and an uncharacterized 'gobloid' subtype consisting predominantly of precursor and mature (Muc2+) goblet cells. To complement this dataset and to demonstrate DisCo's capacity to process low-input, in vivo-derived tissues, we also analyzed individual mouse intestinal crypts. This revealed the existence of crypts with a compositional similarity to spheroids, which consisted predominantly of regenerative stem cells, suggesting the existence of regenerating crypts in the homeostatic intestine. These findings demonstrate the unique power of DisCo in providing high-resolution snapshots of cellular heterogeneity in small, individual tissues.

Single-cell biology: what does the future hold?

In this Editorial, our Chief Editor and members of our Advisory Editorial Board discuss recent breakthroughs, current challenges, and emerging opportunities in single-cell biology and share their vision of "where the field is headed."

Multi-well plate lid for single-step pooling of 96 samples for high-throughput barcode-based sequencing.

High-throughput transcriptomics is of increasing fundamental biological and clinical interest. The generation of molecular data from large collections of samples, such as biobanks and drug libraries, is boosting the development of new biomarkers and treatments. Focusing on gene expression, the transcriptomic market exploits the benefits of next-generation sequencing (NGS), leveraging RNA sequencing (RNA-seq) as standard for measuring genome-wide gene expression in biological samples. The cumbersome sample preparation, including RNA extraction, conversion to cDNA and amplification, prevents high-throughput translation of RNA-seq technologies. Bulk RNA barcoding and sequencing (BRB-seq) addresses this limitation by enabling sample preparation in multi-well plate format. Sample multiplexing combined with early pooling into a single tube reduces reagents consumption and manual steps. Enabling simultaneous pooling of all samples from the multi-well plate into one tube, our technology relies on smart labware: a pooling lid comprising fluidic features and small pins to transport the liquid, adapted to standard 96-well plates. Operated with standard fluidic tubes and pump, the system enables over 90% recovery of liquid in a single step in less than a minute. Large scale manufacturing of the lid is demonstrated with the transition from a milled polycarbonate/steel prototype into an injection molded polystyrene lid. The pooling lid demonstrated its value in supporting high-throughput barcode-based sequencing by pooling 96 different DNA barcodes directly from a standard 96-well plate, followed by processing within the single sample pool. This new pooling technology shows great potential to address medium throughput needs in the BRB-seq workflow, thereby addressing the challenge of large-scale and cost-efficient sample preparation for RNA-seq.

Disparate temperature-dependent virus-host dynamics for SARS-CoV-2 and SARS-CoV in the human respiratory epithelium.

Since its emergence in December 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread globally and become a major public health burden. Despite its close phylogenetic relationship to SARS-CoV, SARS-CoV-2 exhibits increased human-to-human transmission dynamics, likely due to efficient early replication in the upper respiratory epithelium of infected individuals. Since different temperatures encountered in the human upper and lower respiratory tract (33°C and 37°C, respectively) have been shown to affect the replication kinetics of several respiratory viruses, as well as host innate immune response dynamics, we investigated the impact of temperature on SARS-CoV-2 and SARS-CoV infection using the primary human airway epithelial cell culture model. SARS-CoV-2, in contrast to SARS-CoV, replicated to higher titers when infections were performed at 33°C rather than 37°C. Although both viruses were highly sensitive to type I and type III interferon pretreatment, a detailed time-resolved transcriptome analysis revealed temperature-dependent interferon and pro-inflammatory responses induced by SARS-CoV-2 that were inversely proportional to its replication efficiency at 33°C or 37°C. These data provide crucial insight on pivotal virus-host interaction dynamics and are in line with characteristic clinical features of SARS-CoV-2 and SARS-CoV, as well as their respective transmission efficiencies.

A stromal cell population that inhibits adipogenesis in mammalian fat depots.

Adipocyte development and differentiation have an important role in the aetiology of obesity and its co-morbidities1,2. Although multiple studies have investigated the adipogenic stem and precursor cells that give rise to mature adipocytes3-14, our understanding of their in vivo origin and properties is incomplete2,15,16. This is partially due to the highly heterogeneous and unstructured nature of adipose tissue depots17, which has proven difficult to molecularly dissect using classical approaches such as fluorescence-activated cell sorting and Cre-lox lines based on candidate marker genes16,18. Here, using the resolving power of single-cell transcriptomics19 in a mouse model, we reveal distinct subpopulations of adipose stem and precursor cells in the stromal vascular fraction of subcutaneous adipose tissue. We identify one of these subpopulations as CD142+ adipogenesis-regulatory cells, which can suppress adipocyte formation in vivo and in vitro in a paracrine manner. We show that adipogenesis-regulatory cells are refractory to adipogenesis and that they are functionally conserved in humans. Our findings point to a potentially critical role for adipogenesis-regulatory cells in modulating adipose tissue plasticity, which is linked to metabolic control, differential insulin sensitivity and type 2 diabetes.

MarrowQuant 2.0: A Digital Pathology Workflow Assisting Bone Marrow Evaluation in Experimental and Clinical Hematology.

Bone marrow (BM) cellularity assessment is a crucial step in the evaluation of BM trephine biopsies for hematologic and nonhematologic disorders. Clinical assessment is based on a semiquantitative visual estimation of the hematopoietic and adipocytic components by hematopathologists, which does not provide quantitative information on other stromal compartments. In this study, we developed and validated MarrowQuant 2.0, an efficient, user-friendly digital hematopathology workflow integrated within QuPath software, which serves as BM quantifier for 5 mutually exclusive compartments (bone, hematopoietic, adipocytic, and interstitial/microvasculature areas and other) and derives the cellularity of human BM trephine biopsies. Instance segmentation of individual adipocytes is realized through the adaptation of the machine-learning-based algorithm StarDist. We calculated BM compartments and adipocyte size distributions of hematoxylin and eosin images obtained from 250 bone specimens, from control subjects and patients with acute myeloid leukemia or myelodysplastic syndrome, at diagnosis and follow-up, and measured the agreement of cellularity estimates by MarrowQuant 2.0 against visual scores from 4 hematopathologists. The algorithm was capable of robust BM compartment segmentation with an average mask accuracy of 86%, maximal for bone (99%), hematopoietic (92%), and adipocyte (98%) areas. MarrowQuant 2.0 cellularity score and hematopathologist estimations were highly correlated (R2 = 0.92-0.98, intraclass correlation coefficient [ICC] = 0.98; interobserver ICC = 0.96). BM compartment segmentation quantitatively confirmed the reciprocity of the hematopoietic and adipocytic compartments. MarrowQuant 2.0 performance was additionally tested for cellularity assessment of specimens prospectively collected from clinical routine diagnosis. After special consideration for the choice of the cellularity equation in specimens with expanded stroma, performance was similar in this setting (R2 = 0.86, n = 42). Thus, we conclude that these validation experiments establish MarrowQuant 2.0 as a reliable tool for BM cellularity assessment. We expect this workflow will serve as a clinical research tool to explore novel biomarkers related to BM stromal components and may contribute to further validation of future digitalized diagnostic hematopathology workstreams.

LncRNA Ctcflos orchestrates transcription and alternative splicing in thermogenic adipogenesis.

The recruitment of thermogenic brite adipocytes within white adipose tissue attenuates obesity and metabolic comorbidities, arousing interest in understanding the underlying regulatory mechanisms. The molecular network of brite adipogenesis, however, remains largely unresolved. In this light, long noncoding RNAs (lncRNAs) emerged as a versatile class of modulators that control many steps within the differentiation machinery. Leveraging the naturally varying propensities of different inbred mouse strains for white adipose tissue browning, we identify the nuclear lncRNA Ctcflos as a pivotal orchestrator of thermogenic gene expression during brite adipocyte differentiation. Mechanistically, Ctcflos acts as a pleiotropic regulator, being essential for the transcriptional recruitment of the early core thermogenic regulatory program and the modulation of alternative splicing to drive brite adipogenesis. This is showcased by Ctcflos-regulated gene transcription and splicing of the key browning factor Prdm16 toward the isoform that is specific for the thermogenic gene program. Conclusively, our findings emphasize the mechanistic versatility of lncRNAs acting at several independent levels of gene expression for effective regulation of key differentiation factors to direct cell fate and function.

A parallelized, automated platform enabling individual or sequential ChIP of histone marks and transcription factors.

Despite its popularity, chromatin immunoprecipitation followed by sequencing (ChIP-seq) remains a tedious (>2 d), manually intensive, low-sensitivity and low-throughput approach. Here, we combine principles of microengineering, surface chemistry, and molecular biology to address the major limitations of standard ChIP-seq. The resulting technology, FloChIP, automates and miniaturizes ChIP in a beadless fashion while facilitating the downstream library preparation process through on-chip chromatin tagmentation. FloChIP is fast (<2 h), has a wide dynamic range (from 106 to 500 cells), is scalable and parallelized, and supports antibody- or sample-multiplexed ChIP on both histone marks and transcription factors. In addition, FloChIP's interconnected design allows for straightforward chromatin reimmunoprecipitation, which allows this technology to also act as a microfluidic sequential ChIP-seq system. Finally, we ran FloChIP for the transcription factor MEF2A in 32 distinct human lymphoblastoid cell lines, providing insights into the main factors driving collaborative DNA binding of MEF2A and into its role in B cell-specific gene regulation. Together, our results validate FloChIP as a flexible and reproducible automated solution for individual or sequential ChIP-seq.

ASAP 2020 update: an open, scalable and interactive web-based portal for (single-cell) omics analyses.

Single-cell omics enables researchers to dissect biological systems at a resolution that was unthinkable just 10 years ago. However, this analytical revolution also triggered new demands in 'big data' management, forcing researchers to stay up to speed with increasingly complex analytical processes and rapidly evolving methods. To render these processes and approaches more accessible, we developed the web-based, collaborative portal ASAP (Automated Single-cell Analysis Portal). Our primary goal is thereby to democratize single-cell omics data analyses (scRNA-seq and more recently scATAC-seq). By taking advantage of a Docker system to enhance reproducibility, and novel bioinformatics approaches that were recently developed for improving scalability, ASAP meets challenging requirements set by recent cell atlasing efforts such as the Human (HCA) and Fly (FCA) Cell Atlas Projects. Specifically, ASAP can now handle datasets containing millions of cells, integrating intuitive tools that allow researchers to collaborate on the same project synchronously. ASAP tools are versioned, and researchers can create unique access IDs for storing complete analyses that can be reproduced or completed by others. Finally, ASAP does not require any installation and provides a full and modular single-cell RNA-seq analysis pipeline. ASAP is freely available at https://asap.epfl.ch.

Multilevel regulation of the glass locus during Drosophila eye development.

Development of eye tissue is initiated by a conserved set of transcription factors termed retinal determination network (RDN). In the fruit fly Drosophila melanogaster, the zinc-finger transcription factor Glass acts directly downstream of the RDN to control identity of photoreceptor as well as non-photoreceptor cells. Tight control of spatial and temporal gene expression is a critical feature during development, cell-fate determination as well as maintenance of differentiated tissues. The molecular mechanisms that control expression of glass, however, remain largely unknown. We here identify complex regulatory mechanisms controlling expression of the glass locus. All information to recapitulate glass expression are contained in a compact 5.2 kb cis-acting genomic element by combining different cell-type specific and general enhancers with repressor elements. Moreover, the immature RNA of the locus contains an alternative small open reading frame (smORF) upstream of the actual glass translation start, resulting in a small peptide instead of the three possible Glass protein isoforms. CRISPR/Cas9-based mutagenesis shows that the smORF is not required for the formation of functioning photoreceptors, but is able to attenuate effects of glass misexpression. Furthermore, editing the genome to generate glass loci eliminating either one or two isoforms shows that only one of the three proteins is critical for formation of functioning photoreceptors, while removing the two other isoforms did not cause defects in developmental or photoreceptor function. Our results show that eye development and function is largely unaffected by targeted manipulations of critical features of the glass transcript, suggesting a strong selection pressure to allow the formation of a functioning eye.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Inter-embryo gene expression variability recapitulates the hourglass pattern of evo-devo.

BACKGROUND: The evolution of embryological development has long been characterized by deep conservation. In animal development, the phylotypic stage in mid-embryogenesis is more conserved than either early or late stages among species within the same phylum. Hypotheses to explain this hourglass pattern have focused on purifying the selection of gene regulation. Here, we propose an alternative-genes are regulated in different ways at different stages and have different intrinsic capacities to respond to perturbations on gene expression. RESULTS: To eliminate the influence of natural selection, we quantified the expression variability of isogenetic single embryo transcriptomes throughout fly Drosophila melanogaster embryogenesis. We found that the expression variability is lower at the phylotypic stage, supporting that the underlying regulatory architecture in this stage is more robust to stochastic variation on gene expression. We present evidence that the phylotypic stage is also robust to genetic variations on gene expression. Moreover, chromatin regulation appears to play a key role in the variation and evolution of gene expression. CONCLUSIONS: We suggest that a phylum-level pattern of embryonic conservation can be explained by the intrinsic difference of gene regulatory mechanisms in different stages.

SMiLE-seq identifies binding motifs of single and dimeric transcription factors.

Resolving the DNA-binding specificities of transcription factors (TFs) is of critical value for understanding gene regulation. Here, we present a novel, semiautomated protein-DNA interaction characterization technology, selective microfluidics-based ligand enrichment followed by sequencing (SMiLE-seq). SMiLE-seq is neither limited by DNA bait length nor biased toward strong affinity binders; it probes the DNA-binding properties of TFs over a wide affinity range in a fast and cost-effective fashion. We validated SMiLE-seq by analyzing 58 full-length human, mouse, and Drosophila TFs from distinct structural classes. All tested TFs yielded DNA-binding models with predictive power comparable to or greater than that of other in vitro assays. De novo motif discovery on all JUN-FOS heterodimers and several nuclear receptor-TF complexes provided novel insights into partner-specific heterodimer DNA-binding preferences. We also successfully analyzed the DNA-binding properties of uncharacterized human C2H2 zinc-finger proteins and validated several using ChIP-exo.