ASU Scholarship Showcase

Linnorm is a novel normalization and transformation method for the analysis of single cell RNA sequencing (scRNA-seq) data. Linnorm is developed to remove technical noises and simultaneously preserve biological variations in scRNA-seq data, such that existing statistical methods can be improved. Using real scRNA-seq data, we compared Linnorm with existing normalization methods, including NODES, SAMstrt, SCnorm, scran, DESeq and TMM. Linnorm shows advantages in speed, technical noise removal and preservation of cell heterogeneity, which can improve existing methods in the discovery of novel subtypes, pseudo-temporal ordering of cells, clustering analysis, etc. Linnorm also performs better than existing DEG analysis methods, including BASiCS, NODES, SAMstrt, Seurat and DESeq2, in false positive rate control and accuracy.

Two distinct monocyte (Mo)/macrophage (Mp) subsets (Ly6C^low and Ly6C^hi) orchestrate cardiac recovery process following myocardial infarction (MI). Prostaglandin (PG) E₂ is involved in the Mo/Mp-mediated inflammatory response, however, the role of its receptors in Mos/Mps in cardiac healing remains to be determined. Here we show that pharmacological inhibition or gene ablation of the Ep3 receptor in mice suppresses accumulation of Ly6C^low Mos/Mps in infarcted hearts. Ep3 deletion in Mos/Mps markedly attenuates healing after MI by reducing neovascularization in peri-infarct zones. Ep3 deficiency diminishes CX3C chemokine receptor 1 (CX3CR1) expression and vascular endothelial growth factor (VEGF) secretion in Mos/Mps by suppressing TGFβ1 signaling and subsequently inhibits Ly6C^low Mos/Mps migration and angiogenesis. Targeted overexpression of Ep3 receptors in Mos/Mps improves wound healing by enhancing angiogenesis. Thus, the PGE₂/Ep3 axis promotes cardiac healing after MI by activating reparative Ly6C^low Mos/Mps, indicating that Ep3 receptor activation may be a promising therapeutic target for acute MI.

Modeling of transcriptional regulatory networks (TRNs) has been increasingly used to dissect the nature of gene regulation. Inference of regulatory relationships among transcription factors (TFs) and genes, especially among multiple TFs, is still challenging. In this study, we introduced an integrative method, LogicTRN, to decode TF–TF interactions that form TF logics in regulating target genes. By combining cis-regulatory logics and transcriptional kinetics into one single model framework, LogicTRN can naturally integrate dynamic gene expression data and TF-DNA-binding signals in order to identify the TF logics and to reconstruct the underlying TRNs. We evaluated the newly developed methodology using simulation, comparison and application studies, and the results not only show their consistence with existing knowledge, but also demonstrate its ability to accurately reconstruct TRNs in biological complex systems.

Precise spatial positioning and isolation of mammalian cells is a critical component of many single cell experimental methods and biological engineering applications. Although a variety of cell patterning methods have been demonstrated, many of these methods subject cells to high stress environments, discriminate against certain phenotypes, or are a challenge to implement. Here, we demonstrate a rapid, simple, indiscriminate, and minimally perturbing cell patterning method using a laser fabricated polymer stencil. The stencil fabrication process requires no stencil-substrate alignment, and is readily adaptable to various substrate geometries and experiments.

Quantitative three-dimensional (3D) computed tomography (CT) imaging of living single cells enables orientation-independent morphometric analysis of the intricacies of cellular physiology. Since its invention, x-ray CT has become indispensable in the clinic for diagnostic and prognostic purposes due to its quantitative absorption-based imaging in true 3D that allows objects of interest to be viewed and measured from any orientation. However, x-ray CT has not been useful at the level of single cells because there is insufficient contrast to form an image. Recently, optical CT has been developed successfully for fixed cells, but this technology called Cell-CT is incompatible with live-cell imaging due to the use of stains, such as hematoxylin, that are not compatible with cell viability. We present a novel development of optical CT for quantitative, multispectral functional 4D (three spatial + one spectral dimension) imaging of living single cells. The method applied to immune system cells offers truly isotropic 3D spatial resolution and enables time-resolved imaging studies of cells suspended in aqueous medium. Using live-cell optical CT, we found a heterogeneous response to mitochondrial fission inhibition in mouse macrophages and differential basal remodeling of small (0.1 to 1 fl) and large (1 to 20 fl) nuclear and mitochondrial structures on a 20- to 30-s time scale in human myelogenous leukemia cells. Because of its robust 3D measurement capabilities, live-cell optical CT represents a powerful new tool in the biomedical research field.

Functional and molecular cell-to-cell variability is pivotal at the cellular, tissue and whole-organism levels. Yet, the ultimate goal of directly correlating the function of the individual cell with its biomolecular profile remains elusive. We present a platform for integrated analysis of functional and transcriptional phenotypes in the same single cells. We investigated changes in the cellular respiration and gene expression diversity resulting from adaptation to repeated episodes of acute hypoxia in a premalignant progression model. We find differential, progression stage-specific alterations in phenotypic heterogeneity and identify cells with aberrant phenotypes. To our knowledge, this study is the first demonstration of an integrated approach to elucidate how heterogeneity at the transcriptional level manifests in the physiologic profile of individual cells in the context of disease progression.

Driven by an increasing number of studies demonstrating its relevance to a broad variety of disease states, the bioenergy production phenotype has been widely characterized at the bulk sample level. Its cell-to-cell variability, a key player associated with cancer cell survival and recurrence, however, remains poorly understood due to ensemble averaging of the current approaches. We present a technology platform for performing oxygen consumption and extracellular acidification measurements of several hundreds to 1,000 individual cells per assay, while offering simultaneous analysis of cellular communication effects on the energy production phenotype. The platform comprises two major components: a tandem optical sensor for combined oxygen and pH detection, and a microwell device for isolation and analysis of single and few cells in hermetically sealed sub-nanoliter chambers. Our approach revealed subpopulations of cells with aberrant energy production profiles and enables determination of cellular response variability to electron transfer chain inhibitors and ion uncouplers.

Migration is a fundamental cellular behavior that plays an indispensable role in development and homeostasis, but can also contribute to pathology such as cancer metastasis. Due to its relevance to many aspects of human health, the ability to accurately measure cell migration is of broad interest, and numerous approaches have been developed. One of the most commonly employed approaches, because of its simplicity and throughput, is the exclusion zone assay in which cells are allowed to migrate into an initially cell-free region. A major drawback of this assay is that it relies on simply counting cells in the exclusion zone and therefore cannot distinguish the effects of proliferation from migration. We report here a simple modification to the exclusion zone migration assay that exclusively measures cell migration and is not affected by proliferation. This approach makes use of a lineage-tracing vital stain that is retained through cell generations and effectively reads out migration relative to the original, parental cell population. This modification is simple, robust, non-perturbing, and inexpensive. We validate the method in a panel of cell lines under conditions that inhibit or promote migration and demonstrate its use in normal and cancer cell lines as well as primary cells.

In carcinogenesis, intercellular interactions within and between cell types are critical but remain poorly understood. We present a study on intercellular interactions between normal and premalignant epithelial cells and their functional relevance in the context of premalignant to malignant progression in Barrett’s esophagus. Using whole transcriptome profiling we found that in the presence of normal epithelial cells, dysplastic cells but not normal cells, exhibit marked down-regulation of a number of key signaling pathways, including the transforming growth factor beta (TGFβ) and epithelial growth factor (EGF). Functional assays revealed both cell types showed repressed proliferation and significant changes in motility (speed, displacement and directionality) as a result of interactions between the two cell types. Cellular interactions appear to be mediated through both direct cell-cell contact and secreted ligands. The findings of this study are important in that they reveal, for the first time, the effects of cellular communication on gene expression and cellular function between premalignant (dysplastic) epithelial cells and their normal counterparts.