MARIS: method for analyzing RNA following intracellular sorting

Conserved and Divergent Features of Mesenchymal Progenitor Cell Types within the Cortical Nephrogenic Niche of the Human and Mouse Kidney

Cellular interactions among nephron, interstitial, and collecting duct progenitors drive mammalian kidney development. In mice, Six2⁺ nephron progenitor cells (NPCs) and Foxd1⁺ interstitial progenitor cells (IPCs) form largely distinct lineage compartments at the onset of metanephric kidney development. Here, we used the method for analyzing RNA following intracellular sorting (MARIS) approach, single-cell transcriptional profiling, in situ hybridization, and immunolabeling to characterize the presumptive NPC and IPC compartments of the developing human kidney. As in mice, each progenitor population adopts a stereotypical arrangement in the human nephron-forming niche: NPCs capped outgrowing ureteric branch tips, whereas IPCs were sandwiched between the NPCs and the renal capsule. Unlike mouse NPCs, human NPCs displayed a transcriptional profile that overlapped substantially with the IPC transcriptional profile, and key IPC determinants, including FOXD1, were readily detected within SIX2⁺ NPCs. Comparative gene expression profiling in human and mouse Six2/SIX2⁺ NPCs showed broad agreement between the species but also identified species-biased expression of some genes. Notably, some human NPC-enriched genes, including DAPL1 and COL9A2, are linked to human renal disease. We further explored the cellular diversity of mesenchymal cell types in the human nephrogenic niche through single-cell transcriptional profiling. Data analysis stratified NPCs into two main subpopulations and identified a third group of differentiating cells. These findings were confirmed by section in situ hybridization with novel human NPC markers predicted through the single-cell studies. This study provides a benchmark for the mesenchymal progenitors in the human nephrogenic niche and highlights species-variability in kidney developmental programs.

Isolation of intact RNA from murine CD4+ T cells after intracellular cytokine staining and fluorescence-activated cell sorting

Intracellular cytokine staining (ICS) is a powerful method for identifying functionally distinct lymphocyte subsets, and for isolating these by fluorescence activated cell sorting (FACS). Although transcriptomic analysis of cells sorted on the basis of ICS has many potential applications, this is rarely performed because of the difficulty in isolating intact RNA from cells processed using standard fixation and permeabilization buffers for ICS. To address this issue, we compared three buffers shown previously to preserve RNA in nonhematopoietic cells subjected to intracellular staining for their effects on RNA isolated from T lymphocytes processed for ICS. Our results showed that buffers containing the recombinant ribonuclease inhibitor RNasin or high molar concentrations of salt yielded intact RNA from fixed and permeabilized T cells. As proof of principle, we successfully used the buffer containing RNasin to isolate intact RNA from CD4⁺ T cells that were sorted by FACS on the basis of specific cytokine production, thus demonstrating the potential of this approach for coupling ICS with transcriptomic analysis.

Conserved and Divergent Features of Mesenchymal Progenitor Cell Types within the Cortical Nephrogenic Niche of the Human and Mouse Kidney

Cellular interactions among nephron, interstitial, and collecting duct progenitors drive mammalian kidney development. In mice, Six2⁺ nephron progenitor cells (NPCs) and Foxd1⁺ interstitial progenitor cells (IPCs) form largely distinct lineage compartments at the onset of metanephric kidney development. Here, we used the method for analyzing RNA following intracellular sorting (MARIS) approach, single-cell transcriptional profiling, in situ hybridization, and immunolabeling to characterize the presumptive NPC and IPC compartments of the developing human kidney. As in mice, each progenitor population adopts a stereotypical arrangement in the human nephron-forming niche: NPCs capped outgrowing ureteric branch tips, whereas IPCs were sandwiched between the NPCs and the renal capsule. Unlike mouse NPCs, human NPCs displayed a transcriptional profile that overlapped substantially with the IPC transcriptional profile, and key IPC determinants, including FOXD1, were readily detected within SIX2⁺ NPCs. Comparative gene expression profiling in human and mouse Six2/SIX2⁺ NPCs showed broad agreement between the species but also identified species-biased expression of some genes. Notably, some human NPC-enriched genes, including DAPL1 and COL9A2, are linked to human renal disease. We further explored the cellular diversity of mesenchymal cell types in the human nephrogenic niche through single-cell transcriptional profiling. Data analysis stratified NPCs into two main subpopulations and identified a third group of differentiating cells. These findings were confirmed by section in situ hybridization with novel human NPC markers predicted through the single-cell studies. This study provides a benchmark for the mesenchymal progenitors in the human nephrogenic niche and highlights species-variability in kidney developmental programs.

Using intracellular markers to identify a novel set of surface markers for live cell purification from a heterogeneous hIPSC culture

Human embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs) can provide sources for midbrain dopaminergic (mDA) neural progenitors (NPCs) for cell therapy to treat Parkinson's disease (PD) patients. However, the well-known line-to-cell line variability in the differentiation capacity of individual cell lines needs to be improved for the success of this therapy. To address this issue, we sought to identify mDA NPC specific cell surface markers for fluorescence activated cell sorting (FACS). Through RNA isolation after sorting for NPCs based on staining for cell-specific transcription factors followed by microarray, we identified two positive cell surface markers (CORIN and CD166) and one negative cell surface marker (CXCR4) for mDA NPC sorting. These three markers can enrich floor plate NPCs to 90% purity, and the sorted NPCs more efficiently differentiate to mature dopaminergic neurons compared to unsorted or CORIN⁺ alone mDA NPCs. This surface marker identification strategy can be used broadly to facilitate isolation of cell subtypes of interest from heterogeneous cultures.

Precision Medicine for Acute Kidney Injury (AKI): Redefining AKI by Agnostic Kidney Tissue Interrogation and Genetics

Acute kidney injury (AKI) currently is diagnosed by a temporal trend of a single blood analyte: serum creatinine. This measurement is neither sensitive nor specific to kidney injury or its protean forms. Newer biomarkers, neutrophil gelatinase-associated lipocalin (NGAL, Lipocalin 2, Siderocalin), or kidney injury molecule-1 (KIM-1, Hepatitis A Virus Cellular Receptor 1), accelerate the diagnosis of AKI as well as prospectively distinguish rapidly reversible from prolonged causes of serum creatinine increase. Nonetheless, these biomarkers lack the capacity to subfractionate AKI further (eg, sepsis versus ischemia versus nephrotoxicity from medications, enzymes, or metals) or inform us about the primary and secondary sites of injury. It also is unknown whether all nephrons are injured in AKI, whether all cells in a nephron are affected, and whether injury responses can be stimulus-specific or cell type-specific or both. In this review, we summarize fully agnostic tissue interrogation approaches that may help to redefine AKI in cellular and molecular terms, including single-cell and single-nuclei RNA sequencing technology. These approaches will empower a shift in the current paradigm of AKI diagnosis, classification, and staging, and provide the renal community with a significant advance toward precision medicine in the analysis AKI.

Maintaining RNA Integrity for Transcriptomic Profiling of Ex Vivo Cultured Limbal Epithelial Stem Cells after Fluorescence-Activated Cell Sorting (FACS)

Background

Transcriptomic profiling of ex vivo cultured human limbal epithelial stem cells (hLESCs) will foster better understanding of corneal physiology and novel treatment paradigms to limbal stem cell deficiency (LSCD). However, currently such profiling studies are hampered due to difficulties with producing sufficient amounts of intact mRNA for deep RNA sequencing (RNA-seq) from subpopulations sorted on the basis of co-expression of membrane and intracellular antigens by fluorescence-activated cell sorting (FACS).

Methods

To address this problem, we systematically analyzed the critical steps, and found that ethanol fixation together with optimized downstream procedures provided a pipeline that yielded high quality total RNA in amounts to readily support the RNA-seq procedure, while still preserving good discrimination between the individual hLESC immunophenotypes.

Results

The average RNA integrity number (RIN) was 7.7 ± 0.4, and the average yield was 4.6 ± 1.7 pg of RNA per cell. The sequencing analysis of the isolated RNA produced high quality data with more than 70% of read pairs mapping uniformly to the reference genome and 80% of bases with a Phred score of at least 30.

Conclusion

In this study, we developed a reliable FACS-based procedure using ethanol as a fixative that would support accurate isolation of limbal epithelial progenitor subpopulations along with RNA yield and quality sufficient to enable deep transcriptomic profiling.

Proteomic analysis of cell cycle progression in asynchronous cultures, including mitotic subphases, using PRIMMUS

The temporal regulation of protein abundance and post-translational modifications is a key feature of cell division. Recently, we analysed gene expression and protein abundance changes during interphase under minimally perturbed conditions (Ly et al., 2014, 2015). Here, we show that by using specific intracellular immunolabelling protocols, FACS separation of interphase and mitotic cells, including mitotic subphases, can be combined with proteomic analysis by mass spectrometry. Using this PRIMMUS (PRoteomic analysis of Intracellular iMMUnolabelled cell Subsets) approach, we now compare protein abundance and phosphorylation changes in interphase and mitotic fractions from asynchronously growing human cells. We identify a set of 115 phosphorylation sites increased during G2, termed 'early risers'. This set includes phosphorylation of S738 on TPX2, which we show is important for TPX2 function and mitotic progression. Further, we use PRIMMUS to provide the first a proteome-wide analysis of protein abundance remodeling between prophase, prometaphase and anaphase.

Biobanking strategy and sample preprocessing for integrative research in monoclonal gammopathies

AIMS

Some types of monoclonal gammopathies are typified by a very limited availability of aberrant cells. Modern research use high throughput technologies and an integrated approach for detailed characterisation of abnormal cells. This strategy requires relatively high amounts of starting material which cannot be obtained from every diagnosis without causing inconvenience to the patient. The aim of this methodological paper is to reflect our long experience with laboratory work and describe the best protocols for sample collection, sorting and further preprocessing in terms of the available number of cells and intended downstream application in monoclonal gammopathies research. Potential pitfalls are also discussed.

METHODS

Comparison and optimisation of freezing and sorting protocols for plasma cells in monoclonal gammopathies, followed by testing of various nucleic acid isolation and amplification techniques to establish a guideline for sample processing in haemato-oncology research.

RESULTS

We show the average numbers of aberrant cells that can be obtained from various monoclonal gammopathies (monoclonal gammopathy of undetermined significance/light chain amyloidosis/multiple myeloma (MM)/MM circulating plasma cells/ minimal residual disease MM-10 123/22 846/305 501/68 641/4000 aberrant plasma cells of 48/30/10/16/37×106 bone marrow mononuclear cells) and the expected yield of nucleic acids provided from multiple isolation kits (DNA/RNA yield from 1 to 200×103 cells was 2.14-427/0.12-123 ng).

CONCLUSIONS

Tested kits for parallel isolation deliver outputs comparable with kits specialised for just one type of molecule. We also present our positive experience with the whole genome amplification method, which can serve as a very powerful tool to gain complex information from a very small cell population.

Unravelling biology and shifting paradigms in cancer with single-cell sequencing

The fundamental operative unit of a cancer is the genetically and epigenetically innovative single cell. Whether proliferating or quiescent, in the primary tumour mass or disseminated elsewhere, single cells govern the parameters that dictate all facets of the biology of cancer. Thus, single-cell analyses provide the ultimate level of resolution in our quest for a fundamental understanding of this disease. Historically, this quest has been hampered by technological shortcomings. In this Opinion article, we argue that the rapidly evolving field of single-cell sequencing has unshackled the cancer research community of these shortcomings. From furthering an elemental understanding of intra-tumoural genetic heterogeneity and cancer genome evolution to illuminating the governing principles of disease relapse and metastasis, we posit that single-cell sequencing promises to unravel the biology of all facets of this disease.

Insulin-producing cells derived from 'induced pluripotent stem cells' of patients with fulminant type 1 diabetes: Vulnerability to cytokine insults and increased expression of apoptosis-related genes

AIMS/INTRODUCTION

The present study was carried out to generate induced pluripotent stem cells (iPSCs) from patients with fulminant type 1 diabetes, and evaluate the cytokine-induced apoptotic reactions of β-like insulin-producing cells differentiated from the iPSCs.

MATERIALS AND METHODS

iPSCs were generated from fibroblasts of patients with fulminant type 1 diabetes by inducing six reprogramming factors. Insulin-producing cells were differentiated from the iPSCs in vitro. The proportion of cleaved caspase-3-positive or terminal deoxynucleotidyl transferase 2'-deoxyuridine, 5'-triphosphate nick end labeling-positive cells among insulin (INS)-positive cells derived from fulminant type 1 diabetes iPSC and control human iPSC lines was evaluated under treatment with tumor necrosis factor-α, interleukin-1β and interferon-γ. Ribonucleic acid sequencing was carried out to compare gene expressions in INS-positive cells derived from fulminant type 1 diabetes iPSC and control human iPSC lines.

RESULTS

Two iPSC clones were established from each of three patients with fulminant type 1 diabetes. The differentiation of insulin-producing cells from fulminant type 1 diabetes iPSC was confirmed by immunofluorescence analysis and KCl-induced C-peptide secretion. After treatment with pro-inflammatory cytokines, these INS-positive cells showed higher expression of cleaved caspase-3 than those derived from control human iPSCs. Altered expression levels of several apoptosis-related genes were observed in INS-positive cells derived from the fulminant type 1 diabetes iPSCs by ribonucleic acid sequencing.

CONCLUSIONS

We generated iPSCs from patients with fulminant type 1 diabetes and differentiated them into insulin-producing cells. This in vitro disease model can be used to elucidate the disease mechanisms of fulminant type 1 diabetes.

Sox5 regulates beta-cell phenotype and is reduced in type 2 diabetes

Type 2 diabetes (T2D) is characterized by insulin resistance and impaired insulin secretion, but the mechanisms underlying insulin secretion failure are not completely understood. Here, we show that a set of co-expressed genes, which is enriched for genes with islet-selective open chromatin, is associated with T2D. These genes are perturbed in T2D and have a similar expression pattern to that of dedifferentiated islets. We identify Sox5 as a regulator of the module. Sox5 knockdown induces gene expression changes similar to those observed in T2D and diabetic animals and has profound effects on insulin secretion, including reduced depolarization-evoked Ca²⁺-influx and β-cell exocytosis. SOX5 overexpression reverses the expression perturbations observed in a mouse model of T2D, increases the expression of key β-cell genes and improves glucose-stimulated insulin secretion in human islets from donors with T2D. We suggest that human islets in T2D display changes reminiscent of dedifferentiation and highlight SOX5 as a regulator of β-cell phenotype and function.

Interrogating islets in health and disease with single-cell technologies

BACKGROUND

Blood glucose levels are tightly controlled by the coordinated actions of hormone-producing endocrine cells that reside in pancreatic islets. Islet cell malfunction underlies diabetes development and progression. Due to the cellular heterogeneity within islets, it has been challenging to uncover how specific islet cells contribute to glucose homeostasis and diabetes pathogenesis. Recent advances in single-cell technologies and computational methods have opened up new avenues to resolve islet heterogeneity and study islet cell states in health and disease.

SCOPE OF REVIEW

In the past year, a multitude of studies have been published that used single-cell approaches to interrogate the transcriptome and proteome of the different islet cell types. Here, we summarize the conclusions of these studies, as well as discuss the technologies used and the challenges faced with computational analysis of single-cell data from islet studies.

MAJOR CONCLUSIONS

By analyzing single islet cells from rodents and humans at different ages and disease states, the studies reviewed here have provided new insight into endocrine cell function and facilitated a high resolution molecular characterization of poorly understood processes, including regeneration, maturation, and diabetes pathogenesis. Gene expression programs and pathways identified in these studies pave the way for the discovery of new targets and approaches to prevent, monitor, and treat diabetes.

Single-Cell Dynamic Analysis of Mitosis in Haploid Embryonic Stem Cells Shows the Prolonged Metaphase and Its Association with Self-diploidization

The recent establishment of mammalian haploid embryonic stem cells (ESCs) provides new possibilities for genetic screening and for understanding genome evolution and function. However, the dynamics of mitosis in haploid ESCs is still unclear. Here, we report that the duration of mitosis in haploid ESCs, especially the metaphase, is significantly longer than that in diploid ESCs. Delaying mitosis by chemicals increased self-diploidization of haploid ESCs, while shortening mitosis stabilized haploid ESCs. Taken together, our study suggests that the delayed mitosis of haploid ESCs is associated with self-diploidization.

•

Mitosis is prolonged in haploid ESCs, especially in metaphase

•

Self-diploidization might be associated with the prolonged mitosis of haploid ESCs

Resolving host-pathogen interactions by dual RNA-seq

The transcriptome is a powerful proxy for the physiological state of a cell, healthy or diseased. As a result, transcriptome analysis has become a key tool in understanding the molecular changes that accompany bacterial infections of eukaryotic cells. Until recently, such transcriptomic studies have been technically limited to analyzing mRNA expression changes in either the bacterial pathogen or the infected eukaryotic host cell. However, the increasing sensitivity of high-throughput RNA sequencing now enables "dual RNA-seq" studies, simultaneously capturing all classes of coding and noncoding transcripts in both the pathogen and the host. In the five years since the concept of dual RNA-seq was introduced, the technique has been applied to a range of infection models. This has not only led to a better understanding of the physiological changes in pathogen and host during the course of an infection but has also revealed hidden molecular phenotypes of virulence-associated small noncoding RNAs that were not visible in standard infection assays. Here, we use the knowledge gained from these recent studies to suggest experimental and computational guidelines for the design of future dual RNA-seq studies. We conclude this review by discussing prospective applications of the technique.

Intracellular Staining and Flow Cytometry to Identify Lymphocyte Subsets within Murine Aorta, Kidney and Lymph Nodes in a Model of Hypertension

It is now well known that T lymphocytes play a critical role in the development of several cardiovascular diseases1,2,3,4,5. For example, studies from our group have shown that hypertension is associated with an excessive accumulation of T cells in the vessels and kidney during the development of experimental hypertension6. Once in these tissues, T cells produce several cytokines that affect both vascular and renal function leading to vasoconstriction and sodium and water retention1,2. To fully understand how T cells cause cardiovascular and renal diseases, it is important to be able to identify and quantify the specific T cell subsets present in these tissues. T cell subsets are defined by a combination of surface markers, the cytokines they secrete, and the transcription factors they express. The complexity of the T cell population makes flow cytometry and intracellular staining an invaluable technique to dissect the phenotypes of the lymphocytes present in tissues. Here, we provide a detailed protocol to identify the surface and intracellular markers (cytokines and transcription factors) in T cells isolated from murine kidney, aorta and aortic draining lymph nodes in a model of angiotensin II induced hypertension. The following steps are described in detail: isolation of the tissues, generation of the single cell suspensions, ex vivo stimulation, fixation, permeabilization and staining. In addition, several fundamental principles of flow cytometric analyses including choosing the proper controls and appropriate gating strategies are discussed.

(Re)Building a Kidney

(Re)Building a Kidney is a National Institute of Diabetes and Digestive and Kidney Diseases-led consortium to optimize approaches for the isolation, expansion, and differentiation of appropriate kidney cell types and the integration of these cells into complex structures that replicate human kidney function. The ultimate goals of the consortium are two-fold: to develop and implement strategies for in vitro engineering of replacement kidney tissue, and to devise strategies to stimulate regeneration of nephrons in situ to restore failing kidney function. Projects within the consortium will answer fundamental questions regarding human gene expression in the developing kidney, essential signaling crosstalk between distinct cell types of the developing kidney, how to derive the many cell types of the kidney through directed differentiation of human pluripotent stem cells, which bioengineering or scaffolding strategies have the most potential for kidney tissue formation, and basic parameters of the regenerative response to injury. As these projects progress, the consortium will incorporate systematic investigations in physiologic function of in vitro and in vivo differentiated kidney tissue, strategies for engraftment in experimental animals, and development of therapeutic approaches to activate innate reparative responses.

Metabolic heterogeneity of activated beige/brite adipocytes in inguinal adipose tissue

Sustained β3 adrenergic receptor (ADRB3) activation simultaneously upregulates fatty acid synthesis and oxidation in mouse brown, beige, and white adipose tissues; however, the cellular basis of this dual regulation is not known. Treatment of mice with the ADRB3 agonist CL316,243 (CL) increased expression of fatty acid synthase (FASN) and medium chain acyl-CoA dehydrogenase (MCAD) protein within the same cells in classic brown and white adipose tissues. Surprisingly, in inguinal adipose tissue, CL-upregulated FASN and MCAD in distinct cell populations: high MCAD expression occurred in multilocular adipocytes that co-expressed UCP1+, whereas high FASN expression occurred in paucilocular adipocytes lacking detectable UCP1. Genetic tracing with UCP1-cre, however, indicated nearly half of adipocytes with a history of UCP1 expression expressed high levels of FASN without current expression of UCP1. Global transcriptomic analysis of FACS-isolated adipocytes confirmed the presence of distinct anabolic and catabolic phenotypes, and identified differential expression of transcriptional pathways known to regulate lipid synthesis and oxidation. Surprisingly, paternally-expressed genes of the non-classical gene imprinted network were strikingly enriched in anabolic phenotypes, suggesting possible involvement in maintaining the balance of metabolic phenotypes. The results indicate that metabolic heterogeneity is a distinct property of activated beige/brite adipocytes that might be under epigenetic control.

Editing the Neuronal Genome: a CRISPR View of Chromatin Regulation in Neuronal Development, Function, and Plasticity

The dynamic orchestration of gene expression is crucial for the proper differentiation, function, and adaptation of cells. In the brain, transcriptional regulation underlies the incredible diversity of neuronal cell types and contributes to the ability of neurons to adapt their function to the environment. Recently, novel methods for genome and epigenome editing have begun to revolutionize our understanding of gene regulatory mechanisms. In particular, the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system has proven to be a particularly accessible and adaptable technique for genome engineering. Here, we review the use of CRISPR/Cas9 in neurobiology and discuss how these studies have advanced understanding of nervous system development and plasticity. We cover four especially salient applications of CRISPR/Cas9: testing the consequences of enhancer mutations, tagging genes and gene products for visualization in live cells, directly activating or repressing enhancers in vivo, and manipulating the epigenome. In each case, we summarize findings from recent studies and discuss evolving adaptations of the method.

Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes

Blood glucose levels are tightly controlled by the coordinated action of at least four cell types constituting pancreatic islets. Changes in the proportion and/or function of these cells are associated with genetic and molecular pathophysiology of monogenic, type 1, and type 2 (T2D) diabetes. Cellular heterogeneity impedes precise understanding of the molecular components of each islet cell type that govern islet (dys)function, particularly the less abundant delta and gamma/pancreatic polypeptide (PP) cells. Here, we report single-cell transcriptomes for 638 cells from nondiabetic (ND) and T2D human islet samples. Analyses of ND single-cell transcriptomes identified distinct alpha, beta, delta, and PP/gamma cell-type signatures. Genes linked to rare and common forms of islet dysfunction and diabetes were expressed in the delta and PP/gamma cell types. Moreover, this study revealed that delta cells specifically express receptors that receive and coordinate systemic cues from the leptin, ghrelin, and dopamine signaling pathways implicating them as integrators of central and peripheral metabolic signals into the pancreatic islet. Finally, single-cell transcriptome profiling revealed genes differentially regulated between T2D and ND alpha, beta, and delta cells that were undetectable in paired whole islet analyses. This study thus identifies fundamental cell-type–specific features of pancreatic islet (dys)function and provides a critical resource for comprehensive understanding of islet biology and diabetes pathogenesis.

An Abbreviated Protocol for In Vitro Generation of Functional Human Embryonic Stem Cell-Derived Beta-Like Cells

The ability to yield glucose-responsive pancreatic beta-cells from human pluripotent stem cells in vitro will facilitate the development of the cell replacement therapies for the treatment of Type 1 Diabetes. Here, through the sequential in vitro targeting of selected signaling pathways, we have developed an abbreviated five-stage protocol (25–30 days) to generate human Embryonic Stem Cell-Derived Beta-like Cells (ES-DBCs). We showed that Geltrex, as an extracellular matrix, could support the generation of ES-DBCs more efficiently than that of the previously described culture systems. The activation of FGF and Retinoic Acid along with the inhibition of BMP, SHH and TGF-beta led to the generation of 75% NKX6.1⁺/NGN3⁺ Endocrine Progenitors. The inhibition of Notch and tyrosine kinase receptor AXL, and the treatment with Exendin-4 and T3 in the final stage resulted in 35% mono-hormonal insulin positive cells, 1% insulin and glucagon positive cells and 30% insulin and NKX6.1 co-expressing cells. Functionally, ES-DBCs were responsive to high glucose in static incubation and perifusion studies, and could secrete insulin in response to successive glucose stimulations. Mitochondrial metabolic flux analyses using Seahorse demonstrated that the ES-DBCs could efficiently metabolize glucose and generate intracellular signals to trigger insulin secretion. In conclusion, targeting selected signaling pathways for 25–30 days was sufficient to generate ES-DBCs in vitro. The ability of ES-DBCs to secrete insulin in response to glucose renders them a promising model for the in vitro screening of drugs, small molecules or genes that may have potential to influence beta-cell function.

End Sequence Analysis Toolkit (ESAT) expands the extractable information from single-cell RNA-seq data

RNA-seq protocols that focus on transcript termini are well suited for applications in which template quantity is limiting. Here we show that, when applied to end-sequencing data, analytical methods designed for global RNA-seq produce computational artifacts. To remedy this, we created the End Sequence Analysis Toolkit (ESAT). As a test, we first compared end-sequencing and bulk RNA-seq using RNA from dendritic cells stimulated with lipopolysaccharide (LPS). As predicted by the telescripting model for transcriptional bursts, ESAT detected an LPS-stimulated shift to shorter 3′-isoforms that was not evident by conventional computational methods. Then, droplet-based microfluidics was used to generate 1000 cDNA libraries, each from an individual pancreatic islet cell. ESAT identified nine distinct cell types, three distinct β-cell types, and a complex interplay between hormone secretion and vascularization. ESAT, then, offers a much-needed and generally applicable computational pipeline for either bulk or single-cell RNA end-sequencing.

Age-Dependent Pancreatic Gene Regulation Reveals Mechanisms Governing Human β Cell Function

Intensive efforts are focused on identifying regulators of human pancreatic islet cell growth and maturation to accelerate development of therapies for diabetes. After birth, islet cell growth and function are dynamically regulated; however, establishing these age-dependent changes in humans has been challenging. Here, we describe a multimodal strategy for isolating pancreatic endocrine and exocrine cells from children and adults to identify age-dependent gene expression and chromatin changes on a genomic scale. These profiles revealed distinct proliferative and functional states of islet α cells or β cells and histone modifications underlying age-dependent gene expression changes. Expression of SIX2 and SIX3, transcription factors without prior known functions in the pancreas and linked to fasting hyperglycemia risk, increased with age specifically in human islet β cells. SIX2 and SIX3 were sufficient to enhance insulin content or secretion in immature β cells. Our work provides a unique resource to study human-specific regulators of islet cell maturation and function.

Single-cell transcriptomes reveal characteristic features of human pancreatic islet cell types

Pancreatic islets of Langerhans contain several specialized endocrine cell types, which are commonly identified by the expression of single marker genes. However, the established marker genes cannot capture the complete spectrum of cellular heterogeneity in human pancreatic islets, and existing bulk transcriptome datasets provide averages across several cell populations. To dissect the cellular composition of the human pancreatic islet and to establish transcriptomes for all major cell types, we performed single‐cell RNA sequencing on 70 cells sorted from human primary tissue. We used this dataset to validate previously described marker genes at the single‐cell level and to identify specifically expressed transcription factors for all islet cell subtypes. All data are available for browsing and download, thus establishing a useful resource of single‐cell expression profiles for endocrine cells in human pancreatic islets.

Fixed single-cell transcriptomic characterization of human radial glial diversity

The human neocortex is created from diverse intermixed progenitors in the prenatal germinal zones. These progenitors have been difficult to characterize since progenitors—particularly radial glia (RG)—are rare, and are defined by a combination of intracellular markers, position and morphology. To circumvent these problems we developed a method called FRISCR for transcriptome profiling of individual fixed, stained and sorted cells. After validation of FRISCR using human embryonic stem cells, we profiled primary human RG that constitute only 1% of the mid-gestation cortex. These RG could be classified into ventricular zone-enriched RG (vRG) that express ANXA1 and CRYAB, and outer subventricular zone-localized RG (oRG) that express HOPX. Our study identifies the first markers and molecular profiles of vRG and oRG cells, and provides an essential step for understanding molecular networks driving the lineage of human neocortical progenitors. Furthermore, FRISCR allows targeted single-cell transcriptomic profiling of tissues that lack live-cell markers.

Novel Observations From Next-Generation RNA Sequencing of Highly Purified Human Adult and Fetal Islet Cell Subsets

Understanding distinct gene expression patterns of normal adult and developing fetal human pancreatic α- and β-cells is crucial for developing stem cell therapies, islet regeneration strategies, and therapies designed to increase β-cell function in patients with diabetes (type 1 or 2). Toward that end, we have developed methods to highly purify α-, β-, and δ-cells from human fetal and adult pancreata by intracellular staining for the cell-specific hormone content, sorting the subpopulations by flow cytometry, and, using next-generation RNA sequencing, we report the detailed transcriptomes of fetal and adult α- and β-cells. We observed that human islet composition was not influenced by age, sex, or BMI, and transcripts for inflammatory gene products were noted in fetal β-cells. In addition, within highly purified adult glucagon-expressing α-cells, we observed surprisingly high insulin mRNA expression, but not insulin protein expression. This transcriptome analysis from highly purified islet α- and β-cell subsets from fetal and adult pancreata offers clear implications for strategies that seek to increase insulin expression in type 1 and type 2 diabetes.

Genome-wide RNA-Seq of Human Motor Neurons Implicates Selective ER Stress Activation in Spinal Muscular Atrophy

Spinal muscular atrophy (SMA) is caused by mutations in the SMN1 gene. Because this gene is expressed ubiquitously, it remains poorly understood why motor neurons (MNs) are one of the most affected cell types. To address this question, we carried out RNA sequencing studies using fixed, antibody-labeled, and purified MNs produced from control and SMA patient-derived induced pluripotent stem cells (iPSCs). We found SMA-specific changes in MNs, including hyper-activation of the ER stress pathway. Functional studies demonstrated that inhibition of ER stress improves MN survival in vitro even in MNs expressing low SMN. In SMA mice, systemic delivery of an ER stress inhibitor that crosses the blood-brain barrier led to the preservation of spinal cord MNs. Therefore, our study implies that selective activation of ER stress underlies MN death in SMA. Moreover, the approach we have taken would be broadly applicable to the study of disease-prone human cells in heterogeneous cultures.

Testing pancreatic islet function at the single cell level by calcium influx with associated marker expression

Studying the response of islet cells to glucose stimulation is important for understanding cell function in healthy and disease states. Most functional assays are performed on whole islets or cell populations, resulting in averaged observations and loss of information at the single cell level. We demonstrate methods to examine calcium fluxing in individual cells of intact islets in response to multiple glucose challenges. Wild-type mouse islets predominantly contained cells that responded to three (out of three) sequential high glucose challenges, whereas cells of diabetic islets (db/db or NOD) responded less frequently or not at all. Imaged islets were also immunostained for endocrine markers to associate the calcium flux profile of individual cells with gene expression. Wild-type mouse islet cells that robustly fluxed calcium expressed β cell markers (INS/NKX6.1), whereas islet cells that inversely fluxed at low glucose expressed α cell markers (GCG). Diabetic mouse islets showed a higher proportion of dysfunctional β cells that responded poorly to glucose challenges. Most of the failed calcium influx responses in β cells were observed in the second and third high glucose challenges, emphasizing the importance of multiple sequential glucose challenges for assessing the full function of islet cells. Human islet cells were also assessed and showed functional α and β cells. This approach to analyze islet responses to multiple glucose challenges in correlation with gene expression assays expands the understanding of β cell function and the diseased state.

DeCoN: genome-wide analysis of in vivo transcriptional dynamics during pyramidal neuron fate selection in neocortex

Neuronal development requires a complex choreography of transcriptional decisions to obtain specific cellular identities. Realizing the ultimate goal of identifying genome-wide signatures that define and drive specific neuronal fates has been hampered by enormous complexity in both time and space during development. Here, we have paired high-throughput purification of pyramidal neuron subclasses with deep profiling of spatiotemporal transcriptional dynamics during corticogenesis to resolve lineage choice decisions. We identified numerous features ranging from spatial and temporal usage of alternative mRNA isoforms and promoters to a host of mRNA genes modulated during fate specification. Notably, we uncovered numerous long noncoding RNAs with restricted temporal and cell-type-specific expression. To facilitate future exploration, we provide an interactive online database to enable multidimensional data mining and dissemination. This multifaceted study generates a powerful resource and informs understanding of the transcriptional regulation underlying pyramidal neuron diversity in the neocortex.

VIDEO ABSTRACT

Differentiated human stem cells resemble fetal, not adult, β cells

Human pluripotent stem cells (hPSCs) have the potential to generate any human cell type, and one widely recognized goal is to make pancreatic β cells. To this end, comparisons between differentiated cell types produced in vitro and their in vivo counterparts are essential to validate hPSC-derived cells. Genome-wide transcriptional analysis of sorted insulin-expressing (INS(+)) cells derived from three independent hPSC lines, human fetal pancreata, and adult human islets points to two major conclusions: (i) Different hPSC lines produce highly similar INS(+) cells and (ii) hPSC-derived INS(+) (hPSC-INS(+)) cells more closely resemble human fetal β cells than adult β cells. This study provides a direct comparison of transcriptional programs between pure hPSC-INS(+) cells and true β cells and provides a catalog of genes whose manipulation may convert hPSC-INS(+) cells into functional β cells.