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.

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.

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.

We investigate the emergence of extreme events in interdependent networks. We introduce an inter-layer traffic resource competing mechanism to account for the limited capacity associated with distinct network layers. A striking finding is that, when the number of network layers and/or the overlap among the layers are increased, extreme events can emerge in a cascading manner on a global scale. Asymptotically, there are two stable absorption states: a state free of extreme events and a state of full of extreme events, and the transition between them is abrupt. Our results indicate that internal interactions in the multiplex system can yield qualitatively distinct phenomena associated with extreme events that do not occur for independent network layers. An implication is that, e.g., public resource competitions among different service providers can lead to a higher resource requirement than naively expected. We derive an analytical theory to understand the emergence of global-scale extreme events based on the concept of effective betweenness. We also articulate a cost-effective control scheme through increasing the capacity of very few hubs to suppress the cascading process of extreme events so as to protect the entire multi-layer infrastructure against global-scale breakdown.

A remarkable phenomenon in spatiotemporal dynamical systems is chimera state, where the structurally and dynamically identical oscillators in a coupled networked system spontaneously break into two groups, one exhibiting coherent motion and another incoherent. This phenomenon was typically studied in the setting of non-local coupling configurations. We ask what can happen to chimera states under systematic changes to the network structure when links are removed from the network in an orderly fashion but the local coupling topology remains invariant with respect to an index shift. We find the emergence of multicluster chimera states. Remarkably, as a parameter characterizing the amount of link removal is increased, chimera states of distinct numbers of clusters emerge and persist in different parameter regions. We develop a phenomenological theory, based on enhanced or reduced interactions among oscillators in different spatial groups, to explain why chimera states of certain numbers of clusters occur in certain parameter regions. The theoretical prediction agrees well with numerics.

Successful identification of directed dynamical influence in complex systems is relevant to significant problems of current interest. Traditional methods based on Granger causality and transfer entropy have issues such as difficulty with nonlinearity and large data requirement. Recently a framework based on nonlinear dynamical analysis was proposed to overcome these difficulties. We find, surprisingly, that noise can counterintuitively enhance the detectability of directed dynamical influence. In fact, intentionally injecting a proper amount of asymmetric noise into the available time series has the unexpected benefit of dramatically increasing confidence in ascertaining the directed dynamical influence in the underlying system. This result is established based on both real data and model time series from nonlinear ecosystems. We develop a physical understanding of the beneficial role of noise in enhancing detection of directed dynamical influence.

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.