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.

Controlling complex networks has become a forefront research area in network science and engineering. Recent efforts have led to theoretical frameworks of controllability to fully control a network through steering a minimum set of driver nodes. However, in realistic situations not every node is accessible or can be externally driven, raising the fundamental issue of control efficacy: if driving signals are applied to an arbitrary subset of nodes, how many other nodes can be controlled? We develop a framework to determine the control efficacy for undirected networks of arbitrary topology. Mathematically, based on non-singular transformation, we prove a theorem to determine rigorously the control efficacy of the network and to identify the nodes that can be controlled for any given driver nodes. Physically, we develop the picture of diffusion that views the control process as a signal diffused from input signals to the set of controllable nodes. The combination of mathematical theory and physical reasoning allows us not only to determine the control efficacy for model complex networks and a large number of empirical networks, but also to uncover phenomena in network control, e.g., hub nodes in general possess lower control centrality than an average node in undirected networks.

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.

A challenging problem in network science is to control complex networks. In existing frameworks of structural or exact controllability, the ability to steer a complex network toward any desired state is measured by the minimum number of required driver nodes. However, if we implement actual control by imposing input signals on the minimum set of driver nodes, an unexpected phenomenon arises: due to computational or experimental error there is a great probability that convergence to the final state cannot be achieved. In fact, the associated control cost can become unbearably large, effectively preventing actual control from being realized physically. The difficulty is particularly severe when the network is deemed controllable with a small number of drivers. Here we develop a physical controllability framework based on the probability of achieving actual control. Using a recently identified fundamental chain structure underlying the control energy, we offer strategies to turn physically uncontrollable networks into physically controllable ones by imposing slightly augmented set of input signals on properly chosen nodes. Our findings indicate that, although full control can be theoretically guaranteed by the prevailing structural controllability theory, it is necessary to balance the number of driver nodes and control cost to achieve physical control.

In spite of the recent interest and advances in linear controllability of complex networks, controlling nonlinear network dynamics remains an outstanding problem. Here we develop an experimentally feasible control framework for nonlinear dynamical networks that exhibit multistability. The control objective is to apply parameter perturbation to drive the system from one attractor to another, assuming that the former is undesired and the latter is desired. To make our framework practically meaningful, we consider restricted parameter perturbation by imposing two constraints: it must be experimentally realizable and applied only temporarily. We introduce the concept of attractor network, which allows us to formulate a quantifiable controllability framework for nonlinear dynamical networks: a network is more controllable if the attractor network is more strongly connected. We test our control framework using examples from various models of experimental gene regulatory networks and demonstrate the beneficial role of noise in facilitating control.

The histone deacetylase (HDAC) inhibitor vorinostat has received significant attention in recent years as an ‘epigenetic’ drug used to treat solid tumors. However, its mechanisms of action are not entirely understood, particularly with regard to its interaction with the aberrations in 3D nuclear structure that accompany neoplastic progression. We investigated the impact of vorinostat on human esophageal epithelial cell lines derived from normal, metaplastic (pre-cancerous), and malignant tissue. Using a combination of novel optical computed tomography (CT)-based quantitative 3D absorption microscopy and conventional confocal fluorescence microscopy, we show that subjecting malignant cells to vorinostat preferentially alters their 3D nuclear architecture relative to non-cancerous cells. Optical CT (cell CT) imaging of fixed single cells showed that drug-treated cancer cells exhibit significant alterations in nuclear morphometry. Confocal microscopy revealed that vorinostat caused changes in the distribution of H3K9ac-marked euchromatin and H3K9me3-marked constitutive heterochromatin. Additionally, 3D immuno-FISH showed that drug-induced expression of the DNA repair gene MGMT was accompanied by spatial relocation toward the center of the nucleus in the nuclei of metaplastic but not in non-neoplastic cells. Our data suggest that vorinostat’s differential modulation of 3D nuclear architecture in normal and abnormal cells could play a functional role in its anti-cancer action.

Nocturnal cooling of urban areas governs the evolution of thermal state and many thermal-driven environmental issues in cities, especially those suffer strong urban heat island (UHI) effect. Advances in the fundamental understanding of the underlying physics of nighttime UHI involve disentangling complex contributing effects and remains an open challenge. In this study, we develop new numerical algorithms to characterize the thermodynamics of urban nocturnal cooling based on solving the energy balance equations for both the landscape surface and the overlying atmosphere. Further, a scaling law is proposed to relate the UHI intensity to a range of governing mechanisms, including the vertical and horizontal transport of heat in the surface layer, the urban-rural breeze, and the possible urban expansion. The accuracy of proposed methods is evaluated against in-situ urban measurements collected in cities with different geographic and climatic conditions. It is found that the vertical and horizontal contributors modulate the nocturnal UHI at distinct elevation in the atmospheric boundary layer.

It has long been accepted that modern reproductive patterns are likely contributors to breast cancer susceptibility because of their influence on hormones such as estrogen and the importance of these hormones in breast cancer. We conducted a meta-analysis to assess whether this ‘evolutionary mismatch hypothesis’ can explain susceptibility to both estrogen receptor positive (ER-positive) and estrogen receptor negative (ER-negative) cancer. Our meta-analysis includes a total of 33 studies and examines parity, age of first birth and age of menarche broken down by estrogen receptor status. We found that modern reproductive patterns are more closely linked to ER-positive than ER-negative breast cancer. Thus, the evolutionary mismatch hypothesis for breast cancer can account for ER-positive breast cancer susceptibility but not ER-negative breast cancer.

In a meta-analysis published by myself and co-authors, we report differences in the life history risk factors for estrogen receptor negative (ER−) and estrogen receptor positive (ER+) breast cancers. Our meta-analysis did not find the association of ER− breast cancer risk with fast life history characteristics that Hidaka and Boddy suggest in their response to our article. There are a number of possible explanations for the differences between their conclusions and the conclusions we drew from our meta-analysis, including limitations of our meta-analysis and methodological challenges in measuring and categorizing estrogen receptor status. These challenges, along with the association of ER+ breast cancer with slow life history characteristics, may make it challenging to find a clear signal of ER− breast cancer with fast life history characteristics, even if that relationship does exist. The contradictory results regarding breast cancer risk and life history characteristics illustrate a more general challenge in evolutionary medicine: often different sub-theories in evolutionary biology make contradictory predictions about disease risk. In this case, life history models predict that breast cancer risk should increase with faster life history characteristics, while the evolutionary mismatch hypothesis predicts that breast cancer risk should increase with delayed reproduction. Whether life history tradeoffs contribute to ER− breast cancer is still an open question, but current models and several lines of evidence suggest that it is a possibility.

Network reconstruction is a fundamental problem for understanding many complex systems with unknown interaction structures. In many complex systems, there are indirect interactions between two individuals without immediate connection but with common neighbors. Despite recent advances in network reconstruction, we continue to lack an approach for reconstructing complex networks with indirect interactions. Here we introduce a two-step strategy to resolve the reconstruction problem, where in the first step, we recover both direct and indirect interactions by employing the Lasso to solve a sparse signal reconstruction problem, and in the second step, we use matrix transformation and optimization to distinguish between direct and indirect interactions. The network structure corresponding to direct interactions can be fully uncovered. We exploit the public goods game occurring on complex networks as a paradigm for characterizing indirect interactions and test our reconstruction approach. We find that high reconstruction accuracy can be achieved for both homogeneous and heterogeneous networks, and a number of empirical networks in spite of insufficient data measurement contaminated by noise. Although a general framework for reconstructing complex networks with arbitrary types of indirect interactions is yet lacking, our approach opens new routes to separate direct and indirect interactions in a representative complex system.