ASU Scholarship Showcase

Cellular heterogeneity plays a pivotal role in a variety of functional processes in vivo including carcinogenesis. However, our knowledge about cell-to-cell diversity and how differences in individual cells manifest in alterations at the population level remains very limited mainly due to the lack of appropriate tools enabling studies at the single-cell level. We present a study on changes in cellular heterogeneity in the context of pre-malignant progression in response to hypoxic stress. Utilizing pre-malignant progression of Barrett’s esophagus (BE) as a disease model system we studied molecular mechanisms underlying the progression from metaplastic to dysplastic (pre-cancerous) stage. We used newly developed methods enabling measurements of cell-to-cell differences in copy numbers of mitochondrial DNA, expression levels of a set of mitochondrial and nuclear genes involved in hypoxia response pathways, and mitochondrial membrane potential. In contrast to bulk cell studies reported earlier, our study shows significant differences between metaplastic and dysplastic BE cells in both average values and single-cell parameter distributions of mtDNA copy numbers, mitochondrial function, and mRNA expression levels of studied genes. Based on single-cell data analysis, we propose that mitochondria may be one of the key factors in pre-malignant progression in BE.

We develop a framework to uncover and analyse dynamical anomalies from massive, nonlinear and non-stationary time series data. The framework consists of three steps: preprocessing of massive datasets to eliminate erroneous data segments, application of the empirical mode decomposition and Hilbert transform paradigm to obtain the fundamental components embedded in the time series at distinct time scales, and statistical/scaling analysis of the components. As a case study, we apply our framework to detecting and characterizing high-frequency oscillations (HFOs) from a big database of rat electroencephalogram recordings. We find a striking phenomenon: HFOs exhibit on–off intermittency that can be quantified by algebraic scaling laws. Our framework can be generalized to big data-related problems in other fields such as large-scale sensor data and seismic data analysis.

Bismuth drugs, despite being clinically used for decades, surprisingly remain in use and effective for the treatment of Helicobacter pylori infection, even for resistant strains when co-administrated with antibiotics. However, the molecular mechanisms underlying the clinically sustained susceptibility of H. pylori to bismuth drugs remain elusive. Herein, we report that integration of in-house metalloproteomics and quantitative proteomics allows comprehensive uncovering of the bismuth-associated proteomes, including 63 bismuth-binding and 119 bismuth-regulated proteins from Helicobacter pylori, with over 60% being annotated with catalytic functions. Through bioinformatics analysis in combination with bioassays, we demonstrated that bismuth drugs disrupted multiple essential pathways in the pathogen, including ROS defence and pH buffering, by binding and functional perturbation of a number of key enzymes. Moreover, we discovered that HpDnaK may serve as a new target of bismuth drugs to inhibit bacterium-host cell adhesion. The integrative approach we report, herein, provides a novel strategy to unveil the molecular mechanisms of antimicrobial metals against pathogens in general. This study sheds light on the design of new types of antimicrobial agents with multiple targets to tackle the current crisis of antimicrobial resistance.

Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or ‘cerQTL’, and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.

Infection after renal transplantation remains a major cause of morbidity and death, especially infection from the extensively drug-resistant bacteria, A. baumannii. A total of fourteen A. baumannii isolates were isolated from the donors’ preserved fluid from DCD (donation after cardiac death) renal transplantation and four isolates in the recipients’ draining liquid at the Kidney Disease Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, from March 2013 to November 2014. An outbreak of A. baumannii emerging after DCD renal transplantation was tracked to understand the transmission of the pathogen. PFGE displayed similar DNA patterns between isolates from the same hospital. Antimicrobial susceptibility tests against thirteen antimicrobial agents were determined using the K-B diffusion method and eTest. Whole-genome sequencing was applied to investigate the genetic relationship of the isolates. With the clinical data and research results, we concluded that the A. baumannii isolates 3R1 and 3R2 was probably transmitted from the donor who acquired the bacteria during his stay in the ICU, while isolate 4R1 was transmitted from 3R1 and 3R2 via medical manipulation. This study demonstrated the value of integration of clinical profiles with molecular methods in outbreak investigation and their importance in controlling infection and preventing serious complications after DCD transplantation.

Intercellular interactions play a central role at the tissue and whole organism level modulating key cellular functions in normal and disease states. Studies of cell-cell communications are challenging due to ensemble averaging effects brought about by intrinsic heterogeneity in cellular function which requires such studies to be conducted with small populations of cells. Most of the current methods for producing and studying such small cell populations are complex to implement and require skilled personnel limiting their widespread utility in biomedical research labs. We present a simple and rapid method to produce small populations with varying size of epithelial cells (10–50 cells/population) with high-throughput (~1 population/second) on flat surfaces via patterning of extracellular matrix (ECM) proteins and random seeding of cells. We demonstrate that despite inherent limitations of non-contact, drop-on-demand piezoelectric inkjet printing for protein patterning, varying mixtures of ECM proteins can be deposited with high reproducibility and level of control on glass substrates using a set of dynamically adjustable optimized deposition parameters. We demonstrate high consistency for the number of cells per population (~1 cell standard error of mean), the population’s size (~0.2 coefficient of variation) and shape, as well as accurate spatial placement of and distance between colonies of a panel of metaplastic and dysplastic esophageal epithelial cells with differing adhesion and motility characteristics. The number of cells per colony, colony size and shape can be varied by dynamically varying the amount of ECM proteins deposited per spatial location and the number of spatial locations on the substrate. The method is applicable to a broad range of biological and biomedical studies including cell-cell communications, cellular microenvironment, migration, and stimulus response.

Whole genome sequencing (WGS) is a promising strategy to unravel variants or genes responsible for human diseases and traits. However, there is a lack of robust platforms for a comprehensive downstream analysis. In the present study, we first proposed three novel algorithms, sequence gap-filled gene feature annotation, bit-block encoded genotypes and sectional fast access to text lines to address three fundamental problems. The three algorithms then formed the infrastructure of a robust parallel computing framework, KGGSeq, for integrating downstream analysis functions for whole genome sequencing data. KGGSeq has been equipped with a comprehensive set of analysis functions for quality control, filtration, annotation, pathogenic prediction and statistical tests. In the tests with whole genome sequencing data from 1000 Genomes Project, KGGSeq annotated several thousand more reliable non-synonymous variants than other widely used tools (e.g. ANNOVAR and SNPEff). It took only around half an hour on a small server with 10 CPUs to access genotypes of ∼60 million variants of 2504 subjects, while a popular alternative tool required around one day. KGGSeq's bit-block genotype format used 1.5% or less space to flexibly represent phased or unphased genotypes with multiple alleles and achieved a speed of over 1000 times faster to calculate genotypic correlation.

Nonhyperbolicity, as characterized by the coexistence of Kolmogorov-Arnold-Moser (KAM) tori and chaos in the phase space, is generic in classical Hamiltonian systems. An open but fundamental question in physics concerns the relativistic quantum manifestations of nonhyperbolic dynamics. We choose the mushroom billiard that has been mathematically proven to be nonhyperbolic, and study the resonant tunneling dynamics of a massless Dirac fermion. We find that the tunneling rate as a function of the energy exhibits a striking "clustering" phenomenon, where the majority of the values of the rate concentrate on a narrow region, as a result of the chaos component in the classical phase space. Relatively few values of the tunneling rate, however, spread outside the clustering region due to the integrable component. Resonant tunneling of electrons in nonhyperbolic chaotic graphene systems exhibits a similar behavior. To understand these numerical results, we develop a theoretical framework by combining analytic solutions of the Dirac equation in certain integrable domains and physical intuitions gained from current understanding of the quantum manifestations of chaos. In particular, we employ a theoretical formalism based on the concept of self-energies to calculate the tunneling rate and analytically solve the Dirac equation in one dimension as well as in two dimensions for a circular-ring-type of tunneling systems exhibiting integrable dynamics in the classical limit. Because relatively few and distinct classical periodic orbits are present in the integrable component, the corresponding relativistic quantum states can have drastically different behaviors, leading to a wide spread in the values of the tunneling rate in the energy-rate plane. In contrast, the chaotic component has embedded within itself an infinite number of unstable periodic orbits, which provide far more quantum states for tunneling. Due to the nature of chaos, these states are characteristically similar, leading to clustering of the values of the tunneling rate in a narrow band. The appealing characteristic of our work is a demonstration and physical understanding of the "mixed" role played by chaos and regular dynamics in shaping relativistic quantum tunneling dynamics.

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.