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.

Background:
Ketogenic diets are high fat and low carbohydrate or very low carbohydrate diets, which render high production of ketones upon consumption known as nutritional ketosis (NK). Ketosis is also produced during fasting periods, which is known as fasting ketosis (FK). Recently, the combinations of NK and FK, as well as NK alone, have been used as resources for weight loss management and treatment of epilepsy.

Methods:
A crossover study design was applied to 11 healthy individuals, who maintained moderately sedentary lifestyle, and consumed three types of diet randomly assigned over a three-week period. All participants completed the diets in a randomized and counterbalanced fashion. Each weekly diet protocol included three phases: Phase 1 - A mixed diet with ratio of fat: (carbohydrate + protein) by mass of 0.18 or the equivalence of 29% energy from fat from Day 1 to Day 5. Phase 2- A mixed or a high-fat diet with ratio of fat: (carbohydrate + protein) by mass of approximately 0.18, 1.63, or 3.80 on Day 6 or the equivalence of 29%, 79%, or 90% energy from fat, respectively. Phase 3 - A fasting diet with no calorie intake on Day 7. Caloric intake from diets on Day 1 to Day 6 was equal to each individual’s energy expenditure. On Day 7, ketone buildup from FK was measured.

Results:
A statistically significant effect of Phase 2 (Day 6) diet was found on FK of Day 7, as indicated by repeated analysis of variance (ANOVA), F(2,20) = 6.73, p < 0.0058. Using a Fisher LDS pair-wise comparison, higher significant levels of acetone buildup were found for diets with 79% fat content and 90% fat content vs. 29% fat content (with p = 0.00159**, and 0.04435**, respectively), with no significant difference between diets with 79% fat content and 90% fat content. In addition, independent of the diet, a significantly higher ketone buildup capability of subjects with higher resting energy expenditure (R[superscript 2] = 0.92), and lower body mass index (R[superscript 2] = 0.71) was observed during FK.

A relatively unexplored issue in cybersecurity science and engineering is whether there exist intrinsic patterns of cyberattacks. Conventional wisdom favors absence of such patterns due to the overwhelming complexity of the modern cyberspace. Surprisingly, through a detailed analysis of an extensive data set that records the time-dependent frequencies of attacks over a relatively wide range of consecutive IP addresses, we successfully uncover intrinsic spatiotemporal patterns underlying cyberattacks, where the term “spatio” refers to the IP address space. In particular, we focus on analyzing macroscopic properties of the attack traffic flows and identify two main patterns with distinct spatiotemporal characteristics: deterministic and stochastic. Strikingly, there are very few sets of major attackers committing almost all the attacks, since their attack “fingerprints” and target selection scheme can be unequivocally identified according to the very limited number of unique spatiotemporal characteristics, each of which only exists on a consecutive IP region and differs significantly from the others. We utilize a number of quantitative measures, including the flux-fluctuation law, the Markov state transition probability matrix, and predictability measures, to characterize the attack patterns in a comprehensive manner. A general finding is that the attack patterns possess high degrees of predictability, potentially paving the way to anticipating and, consequently, mitigating or even preventing large-scale cyberattacks using macroscopic approaches.

Supply-demand processes take place on a large variety of real-world networked systems ranging from power grids and the internet to social networking and urban systems. In a modern infrastructure, supply-demand systems are constantly expanding, leading to constant increase in load requirement for resources and consequently, to problems such as low efficiency, resource scarcity, and partial system failures. Under certain conditions global catastrophe on the scale of the whole system can occur through the dynamical process of cascading failures. We investigate optimization and resilience of time-varying supply-demand systems by constructing network models of such systems, where resources are transported from the supplier sites to users through various links. Here by optimization we mean minimization of the maximum load on links, and system resilience can be characterized using the cascading failure size of users who fail to connect with suppliers.

We consider two representative classes of supply schemes: load driven supply and fix fraction supply. Our findings are: (1) optimized systems are more robust since relatively smaller cascading failures occur when triggered by external perturbation to the links; (2) a large fraction of links can be free of load if resources are directed to transport through the shortest paths; (3) redundant links in the performance of the system can help to reroute the traffic but may undesirably transmit and enlarge the failure size of the system; (4) the patterns of cascading failures depend strongly upon the capacity of links; (5) the specific location of the trigger determines the specific route of cascading failure, but has little effect on the final cascading size; (6) system expansion typically reduces the efficiency; and (7) when the locations of the suppliers are optimized over a long expanding period, fewer suppliers are required. These results hold for heterogeneous networks in general, providing insights into designing optimal and resilient complex supply-demand systems that expand constantly in time.

Many drugs are effective in the early stage of treatment, but patients develop drug resistance after a certain period of treatment, causing failure of the therapy. An important example is Herceptin, a popular monoclonal antibody drug for breast cancer by specifically targeting human epidermal growth factor receptor 2 (Her2). Here we demonstrate a quantitative binding kinetics analysis of drug-target interactions to investigate the molecular scale origin of drug resistance. Using a surface plasmon resonance imaging, we measured the in situ Herceptin-Her2 binding kinetics in single intact cancer cells for the first time, and observed significantly weakened Herceptin-Her2 interactions in Herceptin-resistant cells, compared to those in Herceptin-sensitive cells. We further showed that the steric hindrance of Mucin-4, a membrane protein, was responsible for the altered drug-receptor binding. This effect of a third molecule on drug-receptor interactions cannot be studied using traditional purified protein methods, demonstrating the importance of the present intact cell-based binding kinetics analysis.

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.

Quantifying the interactions of bacteria with external ligands is fundamental to the understanding of pathogenesis, antibiotic resistance, immune evasion, and mechanism of antimicrobial action. Due to inherent cell-to-cell heterogeneity in a microbial population, each bacterium interacts differently with its environment. This large variability is washed out in bulk assays, and there is a need of techniques that can quantify interactions of bacteria with ligands at the single bacterium level. In this work, we present a label-free and real-time plasmonic imaging technique to measure the binding kinetics of ligand interactions with single bacteria, and perform statistical analysis of the heterogeneity. Using the technique, we have studied interactions of antibodies with single Escherichia coli O157:H7 cells and demonstrated a capability of determining the binding kinetic constants of single live bacteria with ligands, and quantify heterogeneity in a microbial population.

Piezoresistivity is a fundamental property of materials that has found many device applications. Here we report piezoresistivity in double helical DNA molecules. By studying the dependence of molecular conductance and piezoresistivity of single DNA molecules with different sequences and lengths, and performing molecular orbital calculations, we show that the piezoresistivity of DNA is caused by force-induced changes in the π–π electronic coupling between neighbouring bases, and in the activation energy of hole hopping. We describe the results in terms of thermal activated hopping model together with the ladder-based mechanical model for DNA proposed by de Gennes.

Studying the thermoelectric effect in DNA is important for unravelling charge transport mechanisms and for developing relevant applications of DNA molecules. Here we report a study of the thermoelectric effect in single DNA molecules. By varying the molecular length and sequence, we tune the charge transport in DNA to either a hopping- or tunnelling-dominated regimes. The thermoelectric effect is small and insensitive to the molecular length in the hopping regime. In contrast, the thermoelectric effect is large and sensitive to the length in the tunnelling regime. These findings indicate that one may control the thermoelectric effect in DNA by varying its sequence and length. We describe the experimental results in terms of hopping and tunnelling charge transport models.