Single-cell spatial reconstruction reveals global division of labour in the mammalian liver. Nature. 2017;542:352-356. [PMID: 28166538 PMCID: PMC5321580 DOI: 10.1038/nature21065]

Self-Assembled "Off/On" Nanopomegranate for In Vivo Photoacoustic and Fluorescence Imaging: Strategic Arrangement of Kupffer Cells in Mouse Hepatic Lobules

Kupffer cells (KCs), potent scavenger cells located in hepatic sinusoids, constantly phagocytize and degrade foreign materials to maintain metabolism and clearance. Understanding the strategic KC arrangement which links to their spatial location and function in hepatic lobules, the basic functional unit in the liver, is highly valuable for characterizing liver function. However, selectively labeling KCs and characterizing their function in vivo remains challenging. Herein, a fast self-assembled pomegranate structure-like nanoparticle with "nanopomegranate seeds" of dye aggregates has been developed, which has dual-modality "off/on" capability. This nanopomegranate shows good photostability, a high extinction coefficient, a high KC labeling efficiency (98.8%), and better visualization of KC morphology than commercial FluoSpheres. In vivo photoacoustic (PA) and fluorescence imaging consistently visualize that KCs are strategically distributed along the central vein (CV)-portal triad (PT) axis in each liver lobule: more and larger KCs exist in areas closer to the PTs. The high-resolution PA quantitative data further revealed that the density of KCs was linearly dependent on the r _n/ r_max ratio (their relative location along the CV-PT axis) ( R² = 0.7513), and the KC density at the outermost layer is almost 246-fold that at the innermost layer (each layer is 8 μm). Notably, the phagocytic ability of KCs located in layers with r _n/ r_max ratios of 0.167-0.3 varies in a zigzag pattern, as evidenced by their different PA intensities. Additionally, the fluorescence imaging quantitation suggests similar fluorescence activation of nanopomegranate in KCs. Nanopomegranates combined with dual-modality imaging reveal the strategic arrangement of KCs in vivo, greatly extending our understanding of liver physiology.

Stem-cell derived hepatocyte-like cells for the assessment of drug-induced liver injury

Drug-induced liver injury is a major cause of drug discovery failure in clinical trials and a leading cause of liver disease. Current preclinical drug testing does not predict hepatotoxicity which highlights the importance of developing highly predictive cell-based models. The use of stem cell technology and differentiation into hepatocyte-like cells (HLCs) could provide a stable source of hepatocytes for multiple applications, including drug screening. HLCs derived from both embryonic and induced pluripotent stem cells have been used to accurately predict hepatotoxicity as well as to test individual-specific toxicity. Although there are still many limitations, mainly related to the lack of fully maturity of the HLCs derived from pluripotent stem cells, they could provide a relative unlimited and consistent supply of cells with stable phenotype, that could be obtained from different donors, enabling the generation of a library of HLCs representative of the variability of human population.

Lineage Tracing Reveals the Bipotency of SOX9+ Hepatocytes during Liver Regeneration

Elucidation of the role of different cell lineages in the liver could offer avenues to drive liver regeneration. Previous studies showed that SOX9⁺ hepatocytes can differentiate into ductal cells after liver injuries. It is unclear whether SOX9⁺ hepatocytes are uni- or bipotent progenitors at a single-cell level during liver injury. Here, we developed a genetic tracing system to delineate the lineage potential of SOX9⁺ hepatocytes during liver homeostasis and regeneration. Fate-mapping data showed that these SOX9⁺ hepatocytes respond specifically to different liver injuries, with some contributing to a substantial number of ductal cells. Clonal analysis demonstrated that a single SOX9⁺ hepatocyte gives rise to both hepatocytes and ductal cells after liver injury. This study provides direct evidence that SOX9⁺ hepatocytes can serve as bipotent progenitors after liver injury, producing both hepatocytes and ductal cells for liver repair and regeneration.

•

SOX9⁺ hepatocytes respond distinctly to different liver injuries

•

Generation of a Confetti reporter responsive to dual recombinases

•

SOX9⁺ hepatocytes can serve as bipotent progenitors after liver injury

[Dissecting early development of the kidney by single cell transcriptomics]

Each of the billions of the cells in our body exhibits their identity with unique gene expression profile. Recent advances in single cell transcriptomics enable to conduct cell taxonomy identifying new cell types and to re-arrange cells in order of pseudo-time course describing differentiation status of each cell. Even though the cost is still high, the single cell transcriptomics now becomes one of the conventional assays. We have applied the single cell gene expression analysis to dissect human development. In this article, we show our recent progress on a study describing early development of the kidney using human iPS cells by the single cell transcriptomics.

β-catenin-activated hepatocellular carcinomas are addicted to fatty acids

OBJECTIVES

CTNNB1-mutated hepatocellular carcinomas (HCCs) constitute a major part of human HCC and are largely inaccessible to target therapy. Yet, little is known about the metabolic reprogramming induced by β-catenin oncogenic activation in the liver. We aimed to decipher such reprogramming and assess whether it may represent a new avenue for targeted therapy of CTNNB1-mutated HCC.

DESIGN

We used mice with hepatocyte-specific oncogenic activation of β-catenin to evaluate metabolic reprogramming using metabolic fluxes on tumourous explants and primary hepatocytes. We assess the role of Pparα in knock-out mice and analysed the consequences of fatty acid oxidation (FAO) using etomoxir. We explored the expression of the FAO pathway in an annotated human HCC dataset.

RESULTS

β-catenin-activated HCC were not glycolytic but intensively oxidised fatty acids. We found that Pparα is a β-catenin target involved in FAO metabolic reprograming. Deletion of Pparα was sufficient to block the initiation and progression of β-catenin-dependent HCC development. FAO was also enriched in human CTNNB1-mutated HCC, under the control of the transcription factor PPARα.

CONCLUSIONS

FAO induced by β-catenin oncogenic activation in the liver is the driving force of the β-catenin-induced HCC. Inhibiting FAO by genetic and pharmacological approaches blocks HCC development, showing that inhibition of FAO is a suitable therapeutic approach for CTNNB1-mutated HCC.

Modeling Human Digestive Diseases With CRISPR-Cas9-Modified Organoids

Insights into the stem cell niche have allowed researchers to cultivate adult tissue stem cells as organoids that display structural and phenotypic features of healthy and diseased epithelial tissues. Organoids derived from patients' tissues are used as models of disease and to test drugs. CRISPR-Cas9 technology can be used to genetically engineer organoids for studies of monogenic diseases and cancer. We review the derivation of organoids from human gastrointestinal tissues and how CRISPR-Cas9 technology can be used to study these organoids. We discuss burgeoning technologies that are broadening our understanding of diseases of the digestive system.

Identification of Enteroendocrine Regulators by Real-Time Single-Cell Differentiation Mapping

Homeostatic regulation of the intestinal enteroendocrine lineage hierarchy is a poorly understood process. We resolved transcriptional changes during enteroendocrine differentiation in real time at single-cell level using a novel knockin allele of Neurog3, the master regulator gene briefly expressed at the onset of enteroendocrine specification. A bi-fluorescent reporter, Neurog3Chrono, measures time from the onset of enteroendocrine differentiation and enables precise positioning of single-cell transcriptomes along an absolute time axis. This approach yielded a definitive description of the enteroendocrine hierarchy and its sub-lineages, uncovered differential kinetics between sub-lineages, and revealed time-dependent hormonal plasticity in enterochromaffin and L cells. The time-resolved map of transcriptional changes predicted multiple novel molecular regulators. Nine of these were validated by conditional knockout in mice or CRISPR modification in intestinal organoids. Six novel candidate regulators (Sox4, Rfx6, Tox3, Myt1, Runx1t1, and Zcchc12) yielded specific enteroendocrine phenotypes. Our time-resolved single-cell transcriptional map presents a rich resource to unravel enteroendocrine differentiation.

Allele-specific RNA imaging shows that allelic imbalances can arise in tissues through transcriptional bursting

Extensive cell-to-cell variation exists even among putatively identical cells, and there is great interest in understanding how the properties of transcription relate to this heterogeneity. Differential expression from the two gene copies in diploid cells could potentially contribute, yet our ability to measure from which gene copy individual RNAs originated remains limited, particularly in the context of tissues. Here, we demonstrate quantitative, single molecule allele-specific RNA FISH adapted for use on tissue sections, allowing us to determine the chromosome of origin of individual RNA molecules in formaldehyde-fixed tissues. We used this method to visualize the allele-specific expression of Xist and multiple autosomal genes in mouse kidney. By combining these data with mathematical modeling, we evaluated models for allele-specific heterogeneity, in particular demonstrating that apparent expression from only one of the alleles in single cells can arise as a consequence of low-level mRNA abundance and transcriptional bursting.

In mammals, most cells of the body contain two genetic datasets: one from the mother and one from the father, and—in theory—these two sets of information could contribute equally to produce the molecules in a given cell. In practice, however, this is not always the case, which can have dramatic implications for many traits, including visible features (such as fur color) and even disease outcomes. However, it remains difficult to measure the parental origin of individual molecules in a given cell and thus to assess what processes lead to an imbalance of the maternal and paternal contribution. We adapted a microscopy technique—called allele-specific single molecule RNA FISH—that uses a combination of fluorescent tags to specifically label one type of molecule, RNA, depending on its origin, for use in mouse kidney sections. Focusing on RNAs that were previously reported to show imbalance, we performed measurements and combined these with mathematical modeling to quantify imbalance in tissues and explain how these can arise. We found that we could recapitulate the observed imbalances using models that only take into account the random processes that produce RNA, without the need to invoke special regulatory mechanisms to create unequal contributions.

Potential application of cell reprogramming techniques for cancer research

The ability to control the transition from an undifferentiated stem cell to a specific cell fate is one of the key techniques that are required for the application of interventional technologies to regenerative medicine and the treatment of tumors and metastases and of neurodegenerative diseases. Reprogramming technologies, which include somatic cell nuclear transfer, induced pluripotent stem cells, and the direct reprogramming of specific cell lineages, have the potential to alter cell plasticity in translational medicine for cancer treatment. The characterization of cancer stem cells (CSCs), the identification of oncogene and tumor suppressor genes for CSCs, and the epigenetic study of CSCs and their microenvironments are important topics. This review summarizes the application of cell reprogramming technologies to cancer modeling and treatment and discusses possible obstacles, such as genetic and epigenetic alterations in cancer cells, as well as the strategies that can be used to overcome these obstacles to cancer research.

Single-Cell Genomics: A Stepping Stone for Future Immunology Discoveries

The immunology field has invested great efforts and ingenuity to characterize the various immune cell types and elucidate their functions. However, accumulating evidence indicates that current technologies and classification schemes are limited in their ability to account for the functional heterogeneity of immune processes. Single-cell genomics hold the potential to revolutionize the way we characterize complex immune cell assemblies and study their spatial organization, dynamics, clonal distribution, pathways, function, and crosstalks. In this Perspective, we consider recent and forthcoming technological and analytical advances in single-cell genomics and the potential impact of those advances on the future of immunology research and immunotherapy.

Molecular regulation of mammalian hepatic architecture

The essential liver exocrine and endocrine functions require a precise spatial arrangement of the hepatic lobule consisting of the central vein, portal vein, hepatic artery, intrahepatic bile duct system, and hepatocyte zonation. This allows blood to be carried through the liver parenchyma sampled by all hepatocytes and bile produced by the hepatocytes to be carried out of the liver through the intrahepatic bile duct system composed of cholangiocytes. The molecular orchestration of multiple signaling pathways and epigenetic factors is required to set up lineage restriction of the bipotential hepatoblast progenitor into the hepatocyte and cholangiocyte cell lineages, and to further refine cell fate heterogeneity within each cell lineage reflected in the functional heterogeneity of hepatocytes and cholangiocytes. In addition to the complex molecular regulation, there is a complicated morphogenetic choreography observed in building the refined hepatic epithelial architecture. Given the multifaceted molecular and cellular regulation, it is not surprising that impairment of any of these processes can result in acute and chronic hepatobiliary diseases. To enlighten the development of potential molecular and cellular targets for therapeutic options, an understanding of how the intricate hepatic molecular and cellular interactions are regulated is imperative. Here, we review the signaling pathways and epigenetic factors regulating hepatic cell lineages, fates, and epithelial architecture.

Liver Regeneration: Different Sub-Populations of Parenchymal Cells at Play Choreographed by an Injury-Specific Microenvironment

Liver regeneration is crucial for the maintenance of liver functional mass during homeostasis and diseases. In a disease context-dependent manner, liver regeneration is contributed to by hepatocytes or progenitor cells. As long as they are replicatively competent, hepatocytes are the main cell type responsible for supporting liver size homeostasisand regeneration. The concept that all hepatocytes within the lobule have the same proliferative capacity but are differentially recruited according to the localization of the wound, or whether a yet to be defined sub-population of hepatocytes supports regeneration is still debated. In a chronically or severely injured liver, hepatocytes may enter a state of replicative senescence. In such conditions, small biliary cells activate and expand, a process called ductular reaction (DR). Work in the last few decades has demonstrated that DR cells can differentiate into hepatocytes and thereby contribute to parenchymal reconstitution. In this study we will review the molecular mechanisms supporting these two processes to determine potential targets that would be amenable for therapeutic manipulation to enhance liver regeneration.

Understanding tumor ecosystems by single-cell sequencing: promises and limitations

Cellular heterogeneity within and across tumors has been a major obstacle in understanding and treating cancer, and the complex heterogeneity is masked if bulk tumor tissues are used for analysis. The advent of rapidly developing single-cell sequencing technologies, which include methods related to single-cell genome, epigenome, transcriptome, and multi-omics sequencing, have been applied to cancer research and led to exciting new findings in the fields of cancer evolution, metastasis, resistance to therapy, and tumor microenvironment. In this review, we discuss recent advances and limitations of these new technologies and their potential applications in cancer studies.

ALK3 undergoes ligand-independent homodimerization and BMP-induced heterodimerization with ALK2

The bone morphogenetic protein (BMP) type I receptors ALK2 and ALK3 are essential for expression of hepcidin, a key iron regulatory hormone. In mice, hepatocyte-specific Alk2 deficiency leads to moderate iron overload with periportal liver iron accumulation, while hepatocyte-specific Alk3 deficiency leads to severe iron overload with centrilobular liver iron accumulation and a more marked reduction of basal hepcidin levels. The objective of this study was to investigate whether the two receptors have additive roles in hepcidin regulation. Iron overload in mice with hepatocyte-specific Alk2 and Alk3 (Alk2/3) deficiency was characterized and compared to hepatocyte-specific Alk3 deficient mice. Co-immunoprecipitation studies were performed to detect the formation of ALK2 and ALK3 homodimer and heterodimer complexes in vitro in the presence and absence of ligands. The iron overload phenotype of hepatocyte-specific Alk2/3-deficient mice was more severe than that of hepatocyte-specific Alk3-deficient mice. In vitro co-immunoprecipitation studies in Huh7 cells showed that ALK3 can homodimerize in absence of BMP2 or BMP6. In contrast, ALK2 did not homodimerize in either the presence or absence of BMP ligands. However, ALK2 did form heterodimers with ALK3 in the presence of BMP2 or BMP6. ALK3-ALK3 and ALK2-ALK3 receptor complexes induced hepcidin expression in Huh7 cells. Our data indicate that: (I) ALK2 and ALK3 have additive functions in vivo, as Alk2/3 deficiency leads to a greater degree of iron overload than Alk3 deficiency; (II) ALK3, but not ALK2, undergoes ligand-independent homodimerization; (III) the formation of ALK2-ALK3 heterodimers is ligand-dependent and (IV) both receptor complexes functionally induce hepcidin expression in vitro.

High-throughput single-cell transcriptomics on organoids

Three-dimensional (3D) tissues grown in culture from human stem cells offer the incredible opportunity to analyze and manipulate human development, and to generate patient-specific models of disease. Methods to sequence DNA and RNA in single cells are being used to analyze these so-called 'organoid' systems in high-resolution. Single-cell transcriptomics has been used to quantitate the similarity of organoid cells to primary tissue counterparts in the brain, intestine, liver, and kidney, as well as identify cell-specific responses to environmental variables and disease conditions. The merging of these two technologies, single-cell genomics and organoids, will have profound impact on personalized medicine in the near future.

Transcriptional Heterogeneity of Beta Cells in the Intact Pancreas

Pancreatic beta cells have been shown to be heterogeneous at multiple levels. However, spatially interrogating transcriptional heterogeneity in the intact tissue has been challenging. Here, we developed an optimized protocol for single-molecule transcript imaging in the intact pancreas and used it to identify a sub-population of "extreme" beta cells with elevated mRNA levels of insulin and other secretory genes. Extreme beta cells contain higher ribosomal and proinsulin content but lower levels of insulin protein in fasted states, suggesting they may be tuned for basal insulin secretion. They exhibit a distinctive intra-cellular polarization pattern, with elevated mRNA concentrations in an apical ER-enriched compartment, distinct from the localization of nascent and mature proteins. The proportion of extreme cells increases in db/db diabetic mice, potentially facilitating the required increase in basal insulin. Our results thus highlight a sub-population of beta cells that may carry distinct functional roles along physiological and pathological timescales.

Machine learning based classification of cells into chronological stages using single-cell transcriptomics

Age-associated deterioration of cellular physiology leads to pathological conditions. The ability to detect premature aging could provide a window for preventive therapies against age-related diseases. However, the techniques for determining cellular age are limited, as they rely on a limited set of histological markers and lack predictive power. Here, we implement GERAS (GEnetic Reference for Age of Single-cell), a machine learning based framework capable of assigning individual cells to chronological stages based on their transcriptomes. GERAS displays greater than 90% accuracy in classifying the chronological stage of zebrafish and human pancreatic cells. The framework demonstrates robustness against biological and technical noise, as evaluated by its performance on independent samplings of single-cells. Additionally, GERAS determines the impact of differences in calorie intake and BMI on the aging of zebrafish and human pancreatic cells, respectively. We further harness the classification ability of GERAS to identify molecular factors that are potentially associated with the aging of beta-cells. We show that one of these factors, junba, is necessary to maintain the proliferative state of juvenile beta-cells. Our results showcase the applicability of a machine learning framework to classify the chronological stage of heterogeneous cell populations, while enabling detection of candidate genes associated with aging.

From single cells to tissue self-organization

Self-organization is a process by which interacting cells organize and arrange themselves in higher order structures and patterns. To achieve this, cells must have molecular mechanisms to sense their complex local environment and interpret it to respond accordingly. A combination of cell-intrinsic and cell-extrinsic cues are decoded by the single cells dictating their behaviour, their differentiation and symmetry-breaking potential driving development, tissue remodeling and regenerative processes. A unifying property of these self-organized pattern-forming systems is the importance of fluctuations, cell-to-cell variability, or noise. Cell-to-cell variability is an inherent and emergent property of populations of cells that maximize the population performance instead of the individual cell, providing tissues the flexibility to develop and maintain homeostasis in diverse environments. In this review, we will explore the role of self-organization and cell-to-cell variability as fundamental properties of multicellularity-and the requisite of single-cell resolution for its understanding. Moreover, we will analyze how single cells generate emergent multicellular dynamics observed at the tissue level 'travelling' across different scales: spatial, temporal and functional.

Spatial Analysis of Single Fiber Cells of the Developing Ocular Lens Reveals Regulated Heterogeneity of Gene Expression

The developing eye lens presents an exceptional paradigm for spatial transcriptomics. It is composed of highly organized long, slender transparent fiber cells, which differentiate from the edges of the anterior epithelium of the lens (equator), attended by high expression of crystallins, which generates transparency. Every fiber cell, therefore, is an optical unit whose refractive properties derive from its gene activity. Here, we probe this tangible relationship between the gene activity and the phenotype by studying the expression of all known 17 crystallins and 77 other non-crystallin genes in single fiber cells isolated from three states/regions of differentiation, allowing us to follow molecular progression at the single-cell level. The data demonstrate highly variable gene activity in cortical fibers, interposed between the nascent and the terminally differentiated fiber cell transcription. These data suggest that the so-called stochastic, highly heterogeneous gene activity is a regulated intermediate in the realization of a functional phenotype.

Paired-cell sequencing enables spatial gene expression mapping of liver endothelial cells

Spatially resolved single-cell RNA sequencing (scRNAseq) is a powerful approach for inferring connections between a cell's identity and its position in a tissue. We recently combined scRNAseq with spatially mapped landmark genes to infer the expression zonation of hepatocytes. However, determining zonation of small cells with low mRNA content, or without highly expressed landmark genes, remains challenging. Here we used paired-cell sequencing, in which mRNA from pairs of attached mouse cells were sequenced and gene expression from one cell type was used to infer the pairs' tissue coordinates. We applied this method to pairs of hepatocytes and liver endothelial cells (LECs). Using the spatial information from hepatocytes, we reconstructed LEC zonation and extracted a landmark gene panel that we used to spatially map LEC scRNAseq data. Our approach revealed the expression of both Wnt ligands and the Dkk3 Wnt antagonist in distinct pericentral LEC sub-populations. This approach can be used to reconstruct spatial expression maps of non-parenchymal cells in other tissues.

Large-scale reconstruction of cell lineages using single-cell readout of transcriptomes and CRISPR-Cas9 barcodes by scGESTALT

Lineage relationships among the large number of heterogeneous cell types generated during development are difficult to reconstruct in a high-throughput manner. We recently established a method, scGESTALT, that combines cumulative editing of a lineage barcode array by CRISPR-Cas9 with large-scale transcriptional profiling using droplet-based single-cell RNA sequencing (scRNA-seq). The technique generates edits in the barcode array over multiple timepoints using Cas9 and pools of single-guide RNAs (sgRNAs) introduced during early and late zebrafish embryonic development, which distinguishes it from similar Cas9 lineage-tracing methods. The recorded lineages are captured, along with thousands of cellular transcriptomes, to build lineage trees with hundreds of branches representing relationships among profiled cell types. Here, we provide details for (i) generating transgenic zebrafish; (ii) performing multi-timepoint barcode editing; (iii) building scRNA-seq libraries from brain tissue; and (iv) concurrently amplifying lineage barcodes from captured single cells. Generating transgenic lines takes 6 months, and performing barcode editing and generating single-cell libraries involve 7 d of hands-on time. scGESTALT provides a scalable platform to map lineage relationships between cell types in any system that permits genome editing during development, regeneration, or disease.

Identification of spatially associated subpopulations by combining scRNAseq and sequential fluorescence in situ hybridization data

How intrinsic gene-regulatory networks interact with a cell's spatial environment to define its identity remains poorly understood. We developed an approach to distinguish between intrinsic and extrinsic effects on global gene expression by integrating analysis of sequencing-based and imaging-based single-cell transcriptomic profiles, using cross-platform cell type mapping combined with a hidden Markov random field model. We applied this approach to dissect the cell-type- and spatial-domain-associated heterogeneity in the mouse visual cortex region. Our analysis identified distinct spatially associated, cell-type-independent signatures in the glutamatergic and astrocyte cell compartments. Using these signatures to analyze single-cell RNA sequencing data, we identified previously unknown spatially associated subpopulations, which were validated by comparison with anatomical structures and Allen Brain Atlas images.

Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations

The liver is the largest solid organ in the body and is critical for metabolic and immune functions. However, little is known about the cells that make up the human liver and its immune microenvironment. Here we report a map of the cellular landscape of the human liver using single-cell RNA sequencing. We provide the transcriptional profiles of 8444 parenchymal and non-parenchymal cells obtained from the fractionation of fresh hepatic tissue from five human livers. Using gene expression patterns, flow cytometry, and immunohistochemical examinations, we identify 20 discrete cell populations of hepatocytes, endothelial cells, cholangiocytes, hepatic stellate cells, B cells, conventional and non-conventional T cells, NK-like cells, and distinct intrahepatic monocyte/macrophage populations. Together, our study presents a comprehensive view of the human liver at single-cell resolution that outlines the characteristics of resident cells in the liver, and in particular provides a map of the human hepatic immune microenvironment.

The development of single cell RNA sequencing technologies has been instrumental in advancing our understanding of tissue biology. Here, MacParland et al. performed single cell RNA sequencing of human liver samples, and identify distinct populations of intrahepatic macrophages that may play specific roles in liver disease.

Causality Analysis and Cell Network Modeling of Spatial Calcium Signaling Patterns in Liver Lobules

Dynamics as well as localization of Ca²⁺ transients plays a vital role in liver function under homeostatic conditions, repair, and disease. In response to circulating hormonal stimuli, hepatocytes exhibit intracellular Ca²⁺ responses that propagate through liver lobules in a wave-like fashion. Although intracellular processes that control cell autonomous Ca²⁺ spiking behavior have been studied extensively, the intra- and inter-cellular signaling factors that regulate lobular scale spatial patterns and wave-like propagation of Ca²⁺ remain to be determined. To address this need, we acquired images of cytosolic Ca²⁺ transients in 1300 hepatocytes situated across several mouse liver lobules over a period of 1600 s. We analyzed this time series data using correlation network analysis, causal network analysis, and computational modeling, to characterize the spatial distribution of heterogeneity in intracellular Ca²⁺ signaling components as well as intercellular interactions that control lobular scale Ca²⁺ waves. Our causal network analysis revealed that hepatocytes are causally linked to multiple other co-localized hepatocytes, but these influences are not necessarily aligned uni-directionally along the sinusoids. Our computational model-based analysis showed that spatial gradients of intracellular Ca²⁺ signaling components as well as intercellular molecular exchange are required for lobular scale propagation of Ca²⁺ waves. Additionally, our analysis suggested that causal influences of hepatocytes on Ca²⁺ responses of multiple neighbors lead to robustness of Ca²⁺ wave propagation through liver lobules.

Epigenomic map of human liver reveals principles of zonated morphogenic and metabolic control

A deeper epigenomic understanding of spatial organization of cells in human tissues is an important challenge. Here we report the first combined positional analysis of transcriptomes and methylomes across three micro-dissected zones (pericentral, intermediate and periportal) of human liver. We identify pronounced anti-correlated transcriptional and methylation gradients including a core of 271 genes controlling zonated metabolic and morphogen networks and observe a prominent porto-central gradient of DNA methylation at binding sites of 46 transcription factors. The gradient includes an epigenetic and transcriptional Wnt signature supporting the concept of a pericentral hepatocyte regeneration pathway under steady-state conditions. While donors with non-alcoholic fatty liver disease show consistent gene expression differences corresponding to the severity of the disease across all zones, the relative zonated gene expression and DNA methylation patterns remain unchanged. Overall our data provide a wealth of new positional insights into zonal networks controlled by epigenetic and transcriptional gradients in human liver.

Spatial mapping of genomic programs in tissue cells is an important step in the understanding of organ function and disease. Here, the authors provide a spatially resolved epigenomic and transcriptomic map of human liver and show porto-central gradients in metabolic and morphogen networks and transcription factor binding sites as a basis to better understand liver regeneration and function.

The 'TranSeq' 3'-end sequencing method for high-throughput transcriptomics and gene space refinement in plant genomes

High-throughput RNA sequencing has proven invaluable not only to explore gene expression but also for both gene prediction and genome annotation. However, RNA sequencing, carried out on tens or even hundreds of samples, requires easy and cost-effective sample preparation methods using minute RNA amounts. Here, we present TranSeq, a high-throughput 3'-end sequencing procedure that requires 10- to 20-fold fewer sequence reads than the current transcriptomics procedures. TranSeq significantly reduces costs and allows a greater increase in size of sample sets analyzed in a single experiment. Moreover, in comparison with other 3'-end sequencing methods reported to date, we demonstrate here the reliability and immediate applicability of TranSeq and show that it not only provides accurate transcriptome profiles but also produces precise expression measurements of specific gene family members possessing high sequence similarity. This is difficult to achieve in standard RNA-seq methods, in which sequence reads cover the entire transcript. Furthermore, mapping TranSeq reads to the reference tomato genome facilitated the annotation of new transcripts improving >45% of the existing gene models. Hence, we anticipate that using TranSeq will boost large-scale transcriptome assays and increase the spatial and temporal resolution of gene expression data, in both model and non-model plant species. Moreover, as already performed for tomato (ITAG3.0; www.solgenomics.net), we strongly advocate its integration into current and future genome annotations.

Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris

Here we present a compendium of single-cell transcriptomic data from the model organism Mus musculus that comprises more than 100,000 cells from 20 organs and tissues. These data represent a new resource for cell biology, reveal gene expression in poorly characterized cell populations and enable the direct and controlled comparison of gene expression in cell types that are shared between tissues, such as T lymphocytes and endothelial cells from different anatomical locations. Two distinct technical approaches were used for most organs: one approach, microfluidic droplet-based 3'-end counting, enabled the survey of thousands of cells at relatively low coverage, whereas the other, full-length transcript analysis based on fluorescence-activated cell sorting, enabled the characterization of cell types with high sensitivity and coverage. The cumulative data provide the foundation for an atlas of transcriptomic cell biology.

Spatial Reconstruction of Single Enterocytes Uncovers Broad Zonation along the Intestinal Villus Axis

The intestinal epithelium is a highly structured tissue composed of repeating crypt-villus units. Enterocytes perform the diverse tasks of absorbing a wide range of nutrients while protecting the body from the harsh bacterium-rich environment. It is unknown whether these tasks are spatially zonated along the villus axis. Here, we extracted a large panel of landmark genes characterized by transcriptomics of laser capture microdissected villus segments and utilized it for single-cell spatial reconstruction, uncovering broad zonation of enterocyte function along the villus. We found that enterocytes at villus bottoms express an anti-bacterial gene program in a microbiome-dependent manner. They next shift to sequential expression of carbohydrates, peptides, and fat absorption machineries in distinct villus compartments. Finally, they induce a Cd73 immune-modulatory program at the villus tips. Our approach can be used to uncover zonation patterns in other organs when prior knowledge of landmark genes is lacking.

Patient-Derived Induced Pluripotent Stem Cells for Alpha-1 Antitrypsin Deficiency Disease Modeling and Therapeutic Discovery

PIZZ alpha-1 antitrypsin deficiency (AATD) is an autosomal recessive disease affecting approximately 100,000 individuals in the United States and one of the most common hereditary causes of liver disease.¹ The most common form of the disease results from a single base pair mutation (Glu342Lys), known as the "Z" mutation, that encodes a mutant protein (Z alpha-1 antritypsin [AAT]) that is prone to misfolding and is retained in the endoplasmic reticulum (ER) rather than appropriately secreted. Some of the retained mutant protein attains an unusual aggregated or polymerized conformation. Retained polymeric ZAAT aggregates are hepatotoxic and lead to downstream liver disease in a subset of PiZZ neonates and adults through a gain-of-function mechanism. PiZZ individuals are likewise highly predisposed to developing chronic obstructive pulmonary disease (COPD)/emphysema as a result of low circulating levels of AAT protein and associated protease-antiprotease imbalance. Much of our understanding of the molecular pathogenesis of AATD is based on studies employing either transgenic mice that express the mutant human Z allele or immortalized cell lines transduced to overexpress ZAAT. While they have been quite informative, these models fail to capture the patient-to-patient variability in disease phenotype that clinicians observe in their AATD patients, raising the question of whether alternative models might provide new insight. Induced pluripotent stem cells (iPSCs), first described in 2006, have the capacity to differentiate into a broad array of cell types from all 3 germ layers, including hepatocytes. Disease-specific iPSCs have been derived from patients with a variety of monogenic disorders and have been found to faithfully recapitulate features of such diseases as spinal muscular atrophy, familial dysautonomia, Rett syndrome, polycythemia vera, type 1A glycogen storage disease, familial hypercholesterolemia, long QT syndrome, and others. This discussion reviews the potential applications of iPSCs for understanding AATD-associated liver disease as well as for development of potential therapeutic strategies.

The contributions of mesoderm-derived cells in liver development

The liver is an indispensable organ for metabolism and drug detoxification. The liver consists of endoderm-derived hepatobiliary lineages and various mesoderm-derived cells, and interacts with the surrounding tissues and organs through the ventral mesentery. Liver development, from hepatic specification to liver maturation, requires close interactions with mesoderm-derived cells, such as mesothelial cells, hepatic stellate cells, mesenchymal cells, liver sinusoidal endothelial cells and hematopoietic cells. These cells affect liver development through precise signaling events and even direct physical contact. Through the use of new techniques, emerging studies have recently led to a deeper understanding of liver development and its related mechanisms, especially the roles of mesodermal cells in liver development. Based on these developments, the current protocols for in vitro hepatocyte-like cell induction and liver-like tissue construction have been optimized and are of great importance for the treatment of liver diseases. Here, we review the roles of mesoderm-derived cells in the processes of liver development, hepatocyte-like cell induction and liver-like tissue construction.

Cholangiopathies - Towards a molecular understanding

Liver diseases constitute an important medical problem, and a number of these diseases, termed cholangiopathies, affect the biliary system of the liver. In this review, we describe the current understanding of the causes of cholangiopathies, which can be genetic, viral or environmental, and the few treatment options that are currently available beyond liver transplantation. We then discuss recent rapid progress in a number of areas relevant for decoding the disease mechanisms for cholangiopathies. This includes novel data from analysis of transgenic mouse models and organoid systems, and we outline how this information can be used for disease modeling and potential development of novel therapy concepts. We also describe recent advances in genomic and transcriptomic analyses and the importance of such studies for improving diagnosis and determining whether certain cholangiopathies should be viewed as distinct or overlapping disease entities.

Single molecule approaches for studying gene regulation in metabolic tissues

Gene expression in metabolic tissues can be regulated at multiple levels, ranging from the control of promoter accessibilities, transcription rates, mRNA degradation rates and mRNA localization. Modulating these processes can differentially affect important performance criteria of cells. These include precision, cellular economy, rapid response and maintenance of DNA integrity. In this review we will describe how distinct strategies of gene regulation impact the trade-offs between the cells' performance criteria. We will highlight tools based on single molecule visualization of transcripts that can be used to measure promoter states, transcription rates and mRNA degradation rates in intact tissues. These approaches revealed surprising recurrent patterns in mammalian tissues, that include transcriptional bursting, nuclear retention of mRNA, and coordination of mRNA lifetimes to facilitate rapid adaptation to changing metabolic inputs. The ability to characterize gene expression at the single molecule level can uncover the design principles of gene regulation in metabolic tissues such as the liver and the pancreas.

Argininosuccinate synthase 1 and periportal gene expression in sonic hedgehog hepatocellular adenomas

Genetic alterations define different molecular subclasses of hepatocellular adenoma (HCA) linked with risk factors, histology and clinical behavior. Recently, Argininosuccinate Synthase 1 (ASS1), a major periportal protein, was proposed as a marker of HCA with a high risk of hemorrhage. We aimed to assess the significance of ASS1 expression through the scope of the HCA molecular classification. ASS1 expression was evaluated using RNAseq, quantitative reverse transcriptase polymerase chain reaction (RT-PCR) and Immunohistochemistry. ASS1 and glioma-associated oncogene 1 (GLI1) expression were analyzed in vitro after modulation of GLI1 expression. Using RNAseq in 27 HCA and five nontumor liver samples, ASS1 expression was highly correlated with GLI1 expression (P<0.0001, R=0.75). In the overall series of 408 HCA, ASS1 overexpression was significantly associated with sonic hedgehog HCA (shHCA) compared to other molecular subgroups (P<0.0001), suggesting that sonic hedgehog signaling controls ASS1 expression. GLI1 expression silencing by siRNA induced a downregulation of ASS1 in PLC/PFR5 and SNU878 cell lines. In 390 HCA, we showed that ASS1 expression belonged to the periportal expression program that was maintained in shHCA but down-regulated in all the other HCA subtypes. In contrast, HCA with β-catenin activation showed an activation of a perivenous program. Despite the significant association between GLI1 and ASS1 expression, ASS1 mRNA expression was not associated with specific clinical features. At the protein level using immunohistochemistry, prostaglandin D synthase (PTGDS) was strongly and specifically overexpressed in shHCA.

CONCLUSION

ASS1 is associated with sonic hedgehog activation as part of a periportal program expressed in shHCA, a molecular subgroup defined by INHBE-GLI1 gene fusion. (Hepatology 2018).

Single-Cell Multi-omics: An Engine for New Quantitative Models of Gene Regulation

Cells in a multicellular organism fulfill specific functions by enacting cell-type-specific programs of gene regulation. Single-cell RNA sequencing technologies have provided a transformative view of cell-type-specific gene expression, the output of cell-type-specific gene regulatory programs. This review discusses new single-cell genomic technologies that complement single-cell RNA sequencing by providing additional readouts of cellular state beyond the transcriptome. We highlight regression models as a simple yet powerful approach to relate gene expression to other aspects of cellular state, and in doing so, gain insights into the biochemical mechanisms that are necessary to produce a given gene expression output.

Single-cell genomics to guide human stem cell and tissue engineering

To understand human development and disease, as well as to regenerate damaged tissues, scientists are working to engineer certain cell types in vitro and to create 3D microenvironments in which cells behave physiologically. Single-cell genomics (SCG) technologies are being applied to primary human organs and to engineered cells and tissues to generate atlases of cell diversity in these systems at unparalleled resolution. Moving beyond atlases, SCG methods are powerful tools for gaining insight into the engineering and disease process. Here we discuss how scientists can use single-cell sequencing to optimize human cell and tissue engineering by measuring precision, detecting inefficiencies, and assessing accuracy. We also provide a perspective on how emerging SCG methods can be used to reverse-engineer human cells and tissues and unravel disease mechanisms.

A glass-based, continuously zonated and vascularized human liver acinus microphysiological system (vLAMPS) designed for experimental modeling of diseases and ADME/TOX

The vLAMPS is a human, biomimetic liver MPS, in which the ECM and cell seeding of the intermediate layer prior to assembly, simplifies construction of the model and makes the platform user-friendly. This primarily glass microfluidic device is optimal for real-time imaging, while minimizing the binding of hydrophobic drugs/biologics to the materials that constitute the device. The assembly of the three layered device with primary human hepatocytes and liver sinusoidal endothelial cells (LSECs), and human cell lines for stellate and Kupffer cells, creates a vascular channel separated from the hepatic channel (chamber) by a porous membrane that allows communication between channels, recapitulating the 3D structure of the liver acinus. The vascular channel can be used to deliver drugs, immune cells, as well as various circulating cells and other factors to a stand-alone liver MPS and/or to couple the liver MPS to other organ MPS. We have successfully created continuous oxygen zonation by controlling the flow rates of media in the distinct vascular and hepatic channels and validated the computational modeling of zonation with oxygen sensitive and insensitive beads. This allows the direct investigation of the role of zonation in physiology, toxicology and disease progression. The vascular channel is lined with human LSECs, recapitulating partial immunologic functions within the liver sinusoid, including the activation of LSECs, promoting the binding of polymorphonuclear leukocytes (PMNs) followed by transmigration into the hepatic chamber. The vLAMPS is a valuable platform to investigate the functions of the healthy and diseased human liver using all primary human cell types and/or iPSC-derived cells.

Angiocrine Wnt signaling controls liver growth and metabolic maturation in mice

Postnatal liver development is characterized by hepatocyte growth, proliferation, and functional maturation. Notably, canonical Wnt signaling in hepatocytes has been identified as an important regulator of final adult liver size and metabolic liver zonation. The cellular origin of Wnt ligands responsible for homeostatic liver/body weight ratio (LW/BW) remained unclear, which was also attributable to a lack of suitable endothelial Cre driver mice. To comprehensively analyze the effects of hepatic angiocrine Wnt signaling on liver development and metabolic functions, we used endothelial subtype-specific Stab2-Cre driver mice to delete Wls from hepatic endothelial cells (HECs). The resultant Stab2-Cre^tg/wt ;Wls^fl/fl (Wls-HECKO) mice were viable, but showed a significantly reduced LW/BW. Specifically, ablation of angiocrine Wnt signaling impaired metabolic zonation in the liver, as shown by loss of pericentral, β-catenin-dependent target genes such as glutamine synthase (Glul), RhBg, Axin2, and cytochrome P450 2E1, as well as by extended expression of periportal genes such as arginase 1. Furthermore, endothelial subtype-specific expression of a c-terminally YFP-tagged Wls fusion protein in Wls-HECKO mice (Stab2-Cre^tg/wt ;Wls^fl/fl ;Rosa26:Wls-YFP^fl/wt [Wls-rescue]) restored metabolic liver zonation. Interestingly, lipid metabolism was altered in Wls-HECKO mice exhibiting significantly reduced plasma cholesterol levels, while maintaining normal plasma triglyceride and blood glucose concentrations. On the contrary, zonal expression of Endomucin, LYVE1, and other markers of HEC heterogeneity were not altered in Wls-HECKO livers.

CONCLUSION

Angiocrine Wnt signaling controls liver growth as well as development of metabolic liver zonation in mice, whereas intrahepatic HEC zonation is not affected. (Hepatology 2017).

Shotgun Lipidomics Combined with Laser Capture Microdissection: A Tool To Analyze Histological Zones in Cryosections of Tissues

Shotgun analysis provides a quantitative snapshot of the lipidome composition of cells, tissues, or model organisms; however, it does not elucidate the spatial distribution of lipids. Here we demonstrate that shotgun analysis could quantify low-picomole amounts of lipids isolated by laser capture microdissection (LCM) of hundred micrometer-sized histological zones visualized at the cryosections of tissues. We identified metabolically distinct periportal (pp) and pericentral (pc) zones by immunostaining of 20 μm thick cryosections of a healthy mouse liver. LCM was used to ablate, catapult, and collect the tissue material from 10 to 20 individual zones covering a total area of 0.3-0.5 mm² and containing ca. 500 cells. Top-down shotgun profiling relying upon computational stitching of 61 targeted selective ion monitoring ( t-SIM) spectra quantified more than 200 lipid species from 17 lipid classes including glycero- and glycerophospholipids, sphingolipids, cholesterol esters, and cholesterol. Shotgun LCM revealed the overall commonality of the full lipidome composition of pp and pc zones along with significant ( p < 0.001) difference in the relative abundance of 13 lipid species. Follow-up proteomics analyses of pellets recovered from an aqueous phase saved after the lipid extraction identified 13 known and 7 new protein markers exclusively present in pp or in pc zones and independently validated the specificity of their visualization, isolation, and histological assignment.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

Expansion microscopy: enabling single cell analysis in intact biological systems

There is a need for single cell analysis methods that enable the identification and localization of different kinds of biomolecules throughout cells and intact tissues, thereby allowing characterization and classification of individual cells and their relationships to each other within intact systems. Expansion microscopy (ExM) is a technology that physically magnifies tissues in an isotropic way, thereby achieving super-resolution microscopy on diffraction-limited microscopes, enabling rapid image acquisition and large field of view. As a result, ExM is well-positioned to integrate molecular content and cellular morphology, with the spatial precision sufficient to resolve individual biological building blocks, and the scale and accessibility required to deploy over extended 3-D objects like tissues and organs.

The Human Cell Atlas: Technical approaches and challenges

The Human Cell Atlas is a large, international consortium that aims to identify and describe every cell type in the human body. The comprehensive cellular maps that arise from this ambitious effort have the potential to transform many aspects of fundamental biology and clinical practice. Here, we discuss the technical approaches that could be used today to generate such a resource and also the technical challenges that will be encountered.

Transcriptional burst fraction and size dynamics during lens fiber cell differentiation and detailed insights into the denucleation process

Genes are transcribed in irregular pulses of activity termed transcriptional bursts. Cellular differentiation requires coordinated gene expression; however, it is unknown whether the burst fraction (i.e. the number of active phases of transcription) or size/intensity (the number of RNA molecules produced within a burst) changes during cell differentiation. In the ocular lens, the positions of lens fiber cells correlate precisely with their differentiation status, and the most advanced cells degrade their nuclei. Here, we examined the transcriptional parameters of the β-actin and lens differentiation-specific α-, β-, and γ-crystallin genes by RNA fluorescent in situ hybridization (FISH) in the lenses of embryonic day (E) E12.5, E14.5, and E16.5 mouse embryos and newborns. We found that cellular differentiation dramatically alters the burst fraction in synchronized waves across the lens fiber cell compartment with less dramatic changes in burst intensity. Surprisingly, we observed nascent transcription of multiple genes in nuclei just before nuclear destruction. Nuclear condensation was accompanied by transfer of nuclear proteins, including histone and nonhistone proteins, to the cytoplasm. Although lens-specific deletion of the chromatin remodeler SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A member 5 (Smarca5/Snf2h) interfered with denucleation, persisting nuclei remained transcriptionally competent and exhibited changes in both burst intensity and fraction depending on the gene examined. Our results uncover the mechanisms of nascent transcriptional control during differentiation and chromatin remodeling, confirm the burst fraction as the major factor adjusting gene expression levels, and reveal transcriptional competence of fiber cell nuclei even as they approach disintegration.

Single cell transcriptome sequencing: A new approach for the study of mammalian sex determination

Mammalian sex determination is a highly complex developmental process that is particularly difficult to study due to the limited number of gonadal cells present at the bipotential stage, the large cellular heterogeneity in both testis and ovaries and the rapid sex-dependent differentiation processes. Single-cell RNA-sequencing (scRNA-seq) circumvents the averaging artifacts associated with methods traditionally used to profile bulk populations of cells. It is a powerful tool that allows the identification and classification of cell populations in a comprehensive and unbiased manner. In particular, scRNA-seq enables the tracing of cells along developmental trajectories and characterization of the transcriptional dynamics controlling their differentiation. In this review, we describe the current state-of-the-art experimental methods used for scRNA-seq and discuss their strengths and limitations. Additionally, we summarize the multiple key insights that scRNA-seq has provided to the understanding of mammalian sex determination. Finally, we briefly discuss the future of this technology, as well as complementary applications in single cell -omics in the context of mammalian sex determination.