ASU Scholarship Showcase

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.

Recently, the phenomenon of quantum-classical correspondence breakdown was uncovered in optomechanics, where in the classical regime the system exhibits chaos but in the corresponding quantum regime the motion is regular - there appears to be no signature of classical chaos whatsoever in the corresponding quantum system, generating a paradox. We find that transient chaos, besides being a physically meaningful phenomenon by itself, provides a resolution. Using the method of quantum state diffusion to simulate the system dynamics subject to continuous homodyne detection, we uncover transient chaos associated with quantum trajectories. The transient behavior is consistent with chaos in the classical limit, while the long term evolution of the quantum system is regular. Transient chaos thus serves as a bridge for the quantum-classical transition (QCT). Strikingly, as the system transitions from the quantum to the classical regime, the average chaotic transient lifetime increases dramatically (faster than the Ehrenfest time characterizing the QCT for isolated quantum systems). We develop a physical theory to explain the scaling law.

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.

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 role of ambiguity tolerance in career decision making was examined in a sample of college students (n = 275). Three hypotheses were proposed regarding the direct prediction of ambiguity tolerance on career indecision, the indirect prediction of ambiguity tolerance on career indecision through environmental and self explorations, and the moderation effect of ambiguity tolerance on the link of environmental and self explorations with career indecision. Results supported the significance of ambiguity tolerance with respect to career indecision, finding that it directly predicted general indecisiveness, dysfunctional beliefs, lack of information, and inconsistent information, and moderated the prediction of environmental exploration on inconsistent information. The implications of this study are discussed and suggestions for future research are provided.

Type 2 diabetes (T2D) is a complex metabolic disease that is more prevalent in ethnic groups such as Mexican Americans, and is strongly associated with the risk factors obesity and insulin resistance. The goal of this study was to perform whole genome gene expression profiling in adipose tissue to detect common patterns of gene regulation associated with obesity and insulin resistance. We used phenotypic and genotypic data from 308 Mexican American participants from the Veterans Administration Genetic Epidemiology Study (VAGES). Basal fasting RNA was extracted from adipose tissue biopsies from a subset of 75 unrelated individuals, and gene expression data generated on the Illumina BeadArray platform. The number of gene probes with significant expression above baseline was approximately 31,000. We performed multiple regression analysis of all probes with 15 metabolic traits. Adipose tissue had 3,012 genes significantly associated with the traits of interest (false discovery rate, FDR ≤ 0.05). The significance of gene expression changes was used to select 52 genes with significant (FDR ≤ 10^-4) gene expression changes across multiple traits. Gene sets/Pathways analysis identified one gene, alcohol dehydrogenase 1B (ADH1B) that was significantly enriched (P < 10^-60) as a prime candidate for involvement in multiple relevant metabolic pathways. Illumina BeadChip derived ADH1B expression data was consistent with quantitative real time PCR data. We observed significant inverse correlations with waist circumference (2.8 x 10[superscript -9]), BMI (5.4 x 10^-6), and fasting plasma insulin (P < 0.001). These findings are consistent with a central role for ADH1B in obesity and insulin resistance and provide evidence for a novel genetic regulatory mechanism for human metabolic diseases related to these traits.

Evolutionary games model a common type of interactions in a variety of complex, networked, natural systems and social systems. Given such a system, uncovering the interacting structure of the underlying network is key to understanding its collective dynamics. Based on compressive sensing, we develop an efficient approach to reconstructing complex networks under game-based interactions from small amounts of data. The method is validated by using a variety of model networks and by conducting an actual experiment to reconstruct a social network. While most existing methods in this area assume oscillator networks that generate continuous-time data, our work successfully demonstrates that the extremely challenging problem of reverse engineering of complex networks can also be addressed even when the underlying dynamical processes are governed by realistic, evolutionary-game type of interactions in discrete time.

Anthropogenic water sources (AWS) are developed water sources used as a management tool for desert wildlife species. Studies documenting the effects of AWS are often focused on game species; whereas, the effects on non-target wildlife are less understood. We used live trapping techniques to investigate rodent abundance, biomass, and diversity metrics near AWS and paired control sites; we sampled vegetation to determine rodent-habitat associations in the Sauceda Mountains of the Sonoran Desert in Arizona. A total of 370 individual mammals representing three genera and eight species were captured in 4,800 trap nights from winter 2011 to spring 2012.

A multi-response permutation procedure was used to identify differences in small mammal community abundance and biomass by season and treatment. Rodent abundance, biomass, and richness were greater at AWS compared to control sites. Patterns of abundance and biomass were driven by the desert pocket mouse (Chaetodipus penicillatus) which was the most common capture and two times more numerous at AWS compared to controls. Vegetation characteristics, explored using principal components analysis, were similar between AWS and controls. Two species that prefer vegetation structure, Bailey’s pocket mouse (C. baileyi) and white-throated woodrat (Neotoma albigula), had greater abundances and biomass near AWS and were associated with habitat having high cactus density. Although small mammals do not drink free-water, perhaps higher abundances of some species of desert rodents at AWS could be related to artificial structure associated with construction or other resources. Compared to the 30-year average of precipitation for the area, the period of our study occurred during a dry winter. During dry periods, perhaps AWS provide resources to rodents related to moisture.

Since nitrogen (N) is often limiting in permafrost soils, we investigated the N[subscript 2]-fixing genetic potential and the inferred taxa harboring those genes by sequencing nifH gene fragments in samples taken along a permafrost thaw gradient in an Alaskan boreal soil. Samples from minimally, moderately and extensively thawed sites were taken to a depth of 79 cm to encompass zones above and below the depth of the water table. NifH reads were translated with frameshift correction and 112,476 sequences were clustered at 5% amino acid dissimilarity resulting in 1,631 OTUs. Sample depth in relation to water table depth was correlated to differences in the NifH sequence classes with those most closely related to group I nifH-harboring Alpha- and Beta-Proteobacteria in higher abundance above water table depth while those related to group III nifH-harboring Delta Proteobacteria more abundant below. The most dominant below water table depth NifH sequences, comprising 1/3 of the total, were distantly related to Verrucomicrobia-Opitutaceae. Overall, these results suggest that permafrost thaw alters the class-level composition of N[subscript 2]-fixing communities in the thawed soil layers and that this distinction corresponds to the depth of the water table. These nifH data were also compared to nifH sequences obtained from a study at an Alaskan taiga site, and to those of other geographically distant, non-permafrost sites. The two Alaska sites were differentiated largely by changes in relative abundances of the same OTUs, whereas the non-Alaska sites were differentiated by the lack of many Alaskan OTUs, and the presence of unique halophilic, sulfate- and iron-reducing taxa in the Alaska sites.

We examined the effect of different soil sample sizes obtained from an agricultural field, under a single cropping system uniform in soil properties and aboveground crop responses, on bacterial and fungal community structure and microbial diversity indices. DNA extracted from soil sample sizes of 0.25, 1, 5, and 10 g using MoBIO kits and from 10 and 100 g sizes using a bead-beating method (SARDI) were used as templates for high-throughput sequencing of 16S and 28S rRNA gene amplicons for bacteria and fungi, respectively, on the Illumina MiSeq and Roche 454 platforms. Sample size significantly affected overall bacterial and fungal community structure, replicate dispersion and the number of operational taxonomic units (OTUs) retrieved. Richness, evenness and diversity were also significantly affected. The largest diversity estimates were always associated with the 10 g MoBIO extractions with a corresponding reduction in replicate dispersion. For the fungal data, smaller MoBIO extractions identified more unclassified Eukaryota incertae sedis and unclassified glomeromycota while the SARDI method retrieved more abundant OTUs containing unclassified Pleosporales and the fungal genera Alternaria and Cercophora. Overall, these findings indicate that a 10 g soil DNA extraction is most suitable for both soil bacterial and fungal communities for retrieving optimal diversity while still capturing rarer taxa in concert with decreasing replicate variation.