ASU Scholarship Showcase

We develop a general framework to analyze the controllability of multiplex networks using multiple-relation networks and multiple-layer networks with interlayer couplings as two classes of prototypical systems. In the former, networks associated with different physical variables share the same set of nodes and in the latter, diffusion processes take place. We find that, for a multiple-relation network, a layer exists that dominantly determines the controllability of the whole network and, for a multiple-layer network, a small fraction of the interconnections can enhance the controllability remarkably. Our theory is generally applicable to other types of multiplex networks as well, leading to significant insights into the control of complex network systems with diverse structures and interacting patterns.

Background: Healthy individuals on the lower end of the insulin sensitivity spectrum also have a reduced gene expression response to exercise for specific genes. The goal of this study was to determine the relationship between insulin sensitivity and exercise-induced gene expression in an unbiased, global manner.

Methods and Findings: Euglycemic clamps were used to measure insulin sensitivity and muscle biopsies were done at rest and 30 minutes after a single acute exercise bout in 14 healthy participants. Changes in mRNA expression were assessed using microarrays, and miRNA analysis was performed in a subset of 6 of the participants using sequencing techniques. Following exercise, 215 mRNAs were changed at the probe level (Bonferroni-corrected P<0.00000115). Pathway and Gene Ontology analysis showed enrichment in MAP kinase signaling, transcriptional regulation and DNA binding. Changes in several transcription factor mRNAs were correlated with insulin sensitivity, including MYC, r=0.71; SNF1LK, r=0.69; and ATF3, r= 0.61 (5 corrected for false discovery rate). Enrichment in the 5’-UTRs of exercise-responsive genes suggested regulation by common transcription factors, especially EGR1. miRNA species of interest that changed after exercise included miR-378, which is located in an intron of the PPARGC1B gene.

Conclusions: These results indicate that transcription factor gene expression responses to exercise depend highly on insulin sensitivity in healthy people. The overall pattern suggests a coordinated cycle by which exercise and insulin sensitivity regulate gene expression in muscle.

Although insulin resistance in skeletal muscle is well-characterized, the role of circulating whole blood in the metabolic syndrome phenotype is not well understood. We set out to test the hypothesis that genes involved in inflammation, insulin signaling and mitochondrial function would be altered in expression in the whole blood of individuals with metabolic syndrome. We further wanted to examine whether similar relationships that we have found previously in skeletal muscle exist in peripheral whole blood cells. All subjects (n=184) were Latino descent from the Arizona Insulin Resistance registry. Subjects were classified based on the metabolic syndrome phenotype according to the National Cholesterol Education Program’s Adult Treatment Panel III. Of the 184 Latino subjects in the study, 74 were classified with the metabolic syndrome and 110 were without. Whole blood gene expression profiling was performed using the Agilent 4x44K Whole Human Genome Microarray. Whole blood microarray analysis identified 1,432 probes that were altered in expression ≥1.2 fold and P<0.05 after Benjamini-Hochberg in the metabolic syndrome subjects. KEGG pathway analysis revealed significant enrichment for pathways including ribosome, oxidative phosphorylation and MAPK signaling (all Benjamini-Hochberg P<0.05). Whole blood mRNA expression changes observed in the microarray data were confirmed by quantitative RT-PCR. Transcription factor binding motif enrichment analysis revealed E2F1, ELK1, NF-kappaB, STAT1 and STAT3 significantly enriched after Bonferroni correction (all P<0.05). The results of the present study demonstrate that whole blood is a useful tissue for studying the metabolic syndrome and its underlying insulin resistance although the relationship between blood and skeletal muscle differs.

Background: Although the effect of the fat mass and obesity-associated (FTO) gene on adiposity is well established, there is a lack of evidence whether physical activity (PA) modifies the effect of FTO variants on obesity in Latino populations. Therefore, the purpose of this study was to examine PA influences and interactive effects between FTO variants and PA on measures of adiposity in Latinos.

Results: After controlling for age and sex, participants who did not engage in regular PA exhibited higher BMI, fat mass, HC, and WC with statistical significance (P < 0.001). Although significant associations between the three FTO genotypes and adiposity measures were found, none of the FTO genotype by PA interaction assessments revealed nominally significant associations. However, several of such interactive influences exhibited considerable trend towards association.

Conclusions: These data suggest that adiposity measures are associated with PA and FTO variants in Latinos, but the impact of their interactive influences on these obesity measures appear to be minimal. Future studies with large sample sizes may help to determine whether individuals with specific FTO variants exhibit differential responses to PA interventions.

Background: Obesity is a metabolic disease caused by environmental and genetic factors. However, the epigenetic mechanisms of obesity are incompletely understood. The aim of our study was to investigate the role of skeletal muscle DNA methylation in combination with transcriptomic changes in obesity.

Results: Muscle biopsies were obtained basally from lean (n = 12; BMI = 23.4 ± 0.7 kg/m[superscript 2]) and obese (n = 10; BMI = 32.9 ± 0.7 kg/m[superscript 2]) participants in combination with euglycemic-hyperinsulinemic clamps to assess insulin sensitivity. We performed reduced representation bisulfite sequencing (RRBS) next-generation methylation and microarray analyses on DNA and RNA isolated from vastus lateralis muscle biopsies. There were 13,130 differentially methylated cytosines (DMC; uncorrected P < 0.05) that were altered in the promoter and untranslated (5' and 3'UTR) regions in the obese versus lean analysis. Microarray analysis revealed 99 probes that were significantly (corrected P < 0.05) altered. Of these, 12 genes (encompassing 22 methylation sites) demonstrated a negative relationship between gene expression and DNA methylation. Specifically, sorbin and SH3 domain containing 3 (SORBS3) which codes for the adapter protein vinexin was significantly decreased in gene expression (fold change −1.9) and had nine DMCs that were significantly increased in methylation in obesity (methylation differences ranged from 5.0 to 24.4 %). Moreover, differentially methylated region (DMR) analysis identified a region in the 5'UTR (Chr.8:22,423,530–22,423,569) of SORBS3 that was increased in methylation by 11.2 % in the obese group. The negative relationship observed between DNA methylation and gene expression for SORBS3 was validated by a site-specific sequencing approach, pyrosequencing, and qRT-PCR. Additionally, we performed transcription factor binding analysis and identified a number of transcription factors whose binding to the differentially methylated sites or region may contribute to obesity.

Conclusions: These results demonstrate that obesity alters the epigenome through DNA methylation and highlights novel transcriptomic changes in SORBS3 in skeletal muscle.

Introduction: Decreased insulin sensitivity blunts the normal increase in gene expression from skeletal muscle after exercise. In addition, chronic inflammation decreases insulin sensitivity. Chronic kidney disease (CKD) is an inflammatory state. How CKD and, subsequently, kidney transplantation affects skeletal muscle gene expression after exercise are unknown.

Methods: Study cohort: non-diabetic male/female 4/1, age 52±2 years, with end-stage CKD who underwent successful kidney transplantation. The following were measured both pre-transplant and post-transplant and compared to normals: Inflammatory markers, euglycemic insulin clamp studies determine insulin sensitivity, and skeletal muscle biopsies performed before and within 30 minutes after an acute exercise protocol. Microarray analyses were performed on the skeletal muscle using the 4x44K Whole Human Genome Microarrays. Since nuclear factor of activated T cells (NFAT) plays an important role in T cell activation and calcineurin inhibitors are mainstay immunosuppression, calcineurin/NFAT pathway gene expression was compared at rest and after exercise. Log transformation was performed to prevent skewing of data and regression analyses comparing measures pre- and post-transplant performed.

Result: Markers of inflammation significantly improved post-transplantation. Insulin infusion raised glucose disposal slightly lower post-transplant compared to pre-transplant, but not significantly, thus concluding differences in insulin sensitivity were similar. The overall pattern of gene expression in response to exercise was reduced both pre-and post-transplant compared to healthy volunteers. Although significant changes were observed among NFAT/Calcineurin gene at rest and after exercise in normal cohort, there were no significant differences comparing NFAT/calcineurin pathway gene expression pre- and post-transplant.

Conclusions: Despite an improvement in serum inflammatory markers, no significant differences in glucose disposal were observed post-transplant. The reduced skeletal muscle gene expression, including NFAT/calcineurin gene expression, in response to a single bout of exercise was not improved post-transplant. This study suggests that the improvements in inflammatory mediators post-transplant are unrelated to changes of NFAT/calcineurin gene expression.

Dynamical processes occurring on the edges in complex networks are relevant to a variety of real-world situations. Despite recent advances, a framework for edge controllability is still required for complex networks of arbitrary structure and interaction strength. Generalizing a previously introduced class of processes for edge dynamics, the switchboard dynamics, and exploit- ing the exact controllability theory, we develop a universal framework in which the controllability of any node is exclusively determined by its local weighted structure. This framework enables us to identify a unique set of critical nodes for control, to derive analytic formulas and articulate efficient algorithms to determine the exact upper and lower controllability bounds, and to evaluate strongly structural controllability of any given network. Applying our framework to a large number of model and real-world networks, we find that the interaction strength plays a more significant role in edge controllability than the network structure does, due to a vast range between the bounds determined mainly by the interaction strength. Moreover, transcriptional regulatory networks and electronic circuits are much more strongly structurally controllable (SSC) than other types of real-world networks, directed networks are more SSC than undirected networks, and sparse networks are typically more SSC than dense networks.

Recent works revealed that the energy required to control a complex network depends on the number of driving signals and the energy distribution follows an algebraic scaling law. If one implements control using a small number of drivers, e.g. as determined by the structural controllability theory, there is a high probability that the energy will diverge. We develop a physical theory to explain the scaling behaviour through identification of the fundamental structural elements, the longest control chains (LCCs), that dominate the control energy. Based on the LCCs, we articulate a strategy to drastically reduce the control energy (e.g. in a large number of real-world networks). Owing to their structural nature, the LCCs may shed light on energy issues associated with control of nonlinear dynamical networks.

Given a complex geospatial network with nodes distributed in a two-dimensional region of physical space, can the locations of the nodes be determined and their connection patterns be uncovered based solely on data? We consider the realistic situation where time series/signals can be collected from a single location. A key challenge is that the signals collected are necessarily time delayed, due to the varying physical distances from the nodes to the data collection centre. To meet this challenge, we develop a compressive-sensing-based approach enabling reconstruction of the full topology of the underlying geospatial network and more importantly, accurate estimate of the time delays. A standard triangularization algorithm can then be employed to find the physical locations of the nodes in the network. We further demonstrate successful detection of a hidden node (or a hidden source or threat), from which no signal can be obtained, through accurate detection of all its neighbouring nodes. As a geospatial network has the feature that a node tends to connect with geophysically nearby nodes, the localized region that contains the hidden node can be identified.

Locating sources of diffusion and spreading from minimum data is a significant problem in network science with great applied values to the society. However, a general theoretical framework dealing with optimal source localization is lacking. Combining the controllability theory for complex networks and compressive sensing, we develop a framework with high efficiency and robustness for optimal source localization in arbitrary weighted networks with arbitrary distribution of sources. We offer a minimum output analysis to quantify the source locatability through a minimal number of messenger nodes that produce sufficient measurement for fully locating the sources. When the minimum messenger nodes are discerned, the problem of optimal source localization becomes one of sparse signal reconstruction, which can be solved using compressive sensing. Application of our framework to model and empirical networks demonstrates that sources in homogeneous and denser networks are more readily to be located. A surprising finding is that, for a connected undirected network with random link weights and weak noise, a single messenger node is sufficient for locating any number of sources. The framework deepens our understanding of the network source localization problem and offers efficient tools with broad applications.