ASU Scholarship Showcase

Two classes of scaling behaviours, namely the super-linear scaling of links or activities, and the sub-linear scaling of area, diversity, or time elapsed with respect to size have been found to prevail in the growth of complex networked systems. Despite some pioneering modelling approaches proposed for specific systems, whether there exists some general mechanisms that account for the origins of such scaling behaviours in different contexts, especially in socioeconomic systems, remains an open question. We address this problem by introducing a geometric network model without free parameter, finding that both super-linear and sub-linear scaling behaviours can be simultaneously reproduced and that the scaling exponents are exclusively determined by the dimension of the Euclidean space in which the network is embedded. We implement some realistic extensions to the basic model to offer more accurate predictions for cities of various scaling behaviours and the Zipf distribution reported in the literature and observed in our empirical studies. All of the empirical results can be precisely recovered by our model with analytical predictions of all major properties. By virtue of these general findings concerning scaling behaviour, our models with simple mechanisms gain new insights into the evolution and development of complex networked systems.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

The heterocyclic indole-alkaloid scytonemin is a sunscreen found exclusively among cyanobacteria. An 18-gene cluster is responsible for scytonemin production in Nostoc punctiforme ATCC 29133. The upstream genes scyABCDEF in the cluster are proposed to be responsible for scytonemin biosynthesis from aromatic amino acid substrates. In vitro studies of ScyA, ScyB, and ScyC proved that these enzymes indeed catalyze initial pathway reactions. Here we characterize the role of ScyD, ScyE, and ScyF, which were logically predicted to be responsible for late biosynthetic steps, in the biological context of N. punctiforme. In-frame deletion mutants of each were constructed (ΔscyD, ΔscyE, and ΔscyF) and their phenotypes studied. Expectedly, ΔscyE presents a scytoneminless phenotype, but no accumulation of the predicted intermediaries. Surprisingly, ΔscyD retains scytonemin production, implying that it is not required for biosynthesis. Indeed, scyD presents an interesting evolutionary paradox: it likely originated in a duplication event from scyE, and unlike other genes in the operon, it has not been subjected to purifying selection. This would suggest that it is a pseudogene, and yet scyD is highly conserved in the scytonemin operon of cyanobacteria. ΔscyF also retains scytonemin production, albeit exhibiting a reduction of the production yield compared with the wild-type. This indicates that ScyF is not essential but may play an adjuvant role for scytonemin synthesis. Altogether, our findings suggest that these downstream genes are not responsible, as expected, for the late steps of scytonemin synthesis and we must look for those functions elsewhere. These findings are particularly important for biotechnological production of this sunscreen through heterologous expression of its genes in more tractable organisms.

N₂ fixation and ammonia oxidation (AO) are the two most important processes in the nitrogen (N) cycle of biological soil crusts (BSCs). We studied the short-term response of acetylene reduction assay (ARA) rates, an indicator of potential N₂ fixation, and AO rates to temperature (T, -5°C to 35°C) in BSC of different successional stages along the BSC ecological succession and geographic origin (hot Chihuahuan and cooler Great Basin deserts). ARA in all BSCs increased with T until saturation occurred between 15 and 20°C, and declined at 30–35°C. Culture studies using cyanobacteria isolated from these crusts indicated that the saturating effect was traceable to their inability to grow well diazotrophically within the high temperature range. Below saturation, temperature response was exponential, with Q₁₀ significantly different in the two areas (~ 5 for Great Basin BSCs; 2–3 for Chihuahuan BSCs), but similar between the two successional stages. However, in contrast to ARA, AO showed a steady increase to 30–35°C in Great Basin, and Chihuhuan BSCs showed no inhibition at any tested temperature. The T response of AO also differed significantly between Great Basin (Q₁₀ of 4.5–4.8) and Chihuahuan (Q₁₀ of 2.4–2.6) BSCs, but not between successional stages. Response of ARA rates to T did not differ from that of AO in either desert. Thus, while both processes scaled to T in unison until 20°C, they separated to an increasing degree at higher temperature. As future warming is likely to occur in the regions where BSCs are often the dominant living cover, this predicted decoupling is expected to result in higher proportion of nitrates in soil relative to ammonium. As nitrate is more easily lost as leachate or to be reduced to gaseous forms, this could mean a depletion of soil N over large landscapes globally.

Background: Medical and public health scientists are using evolution to devise new strategies to solve major health problems. But based on a 2003 survey, medical curricula may not adequately prepare physicians to evaluate and extend these advances. This study assessed the change in coverage of evolution in North American medical schools since 2003 and identified opportunities for enriching medical education.

Methods: In 2013, curriculum deans for all North American medical schools were invited to rate curricular coverage and perceived importance of 12 core principles, the extent of anticipated controversy from adding evolution, and the usefulness of 13 teaching resources. Differences between schools were assessed by Pearson’s chi-square test, Student’s t-test, and Spearman’s correlation. Open-ended questions sought insight into perceived barriers and benefits.

Results: Despite repeated follow-up, 60 schools (39%) responded to the survey. There was no evidence of sample bias. The three evolutionary principles rated most important were antibiotic resistance, environmental mismatch, and somatic selection in cancer. While importance and coverage of principles were correlated (r = 0.76, P < 0.01), coverage (at least moderate) lagged behind importance (at least moderate) by an average of 21% (SD = 6%). Compared to 2003, a range of evolutionary principles were covered by 4 to 74% more schools. Nearly half (48%) of responders anticipated igniting controversy at their medical school if they added evolution to their curriculum. The teaching resources ranked most useful were model test questions and answers, case studies, and model curricula for existing courses/rotations. Limited resources (faculty expertise) were cited as the major barrier to adding more evolution, but benefits included a deeper understanding and improved patient care.

Conclusion: North American medical schools have increased the evolution content in their curricula over the past decade. However, coverage is not commensurate with importance. At a few medical schools, anticipated controversy impedes teaching more evolution. Efforts to improve evolution education in medical schools should be directed toward boosting faculty expertise and crafting resources that can be easily integrated into existing curricula.

Cyanobacteria are considered good models for biohydrogen production because they are relatively simple organisms with a demonstrable ability to generate H₂ under certain physiological conditions. However, most produce only little H₂, revert readily to H₂ consumption, and suffer from hydrogenase sensitivity to O₂. Strains of the cyanobacteria Lyngbya aestuarii and Microcoleus chthonoplastes obtained from marine intertidal cyanobacterial mats were recently found to display much better H₂ production potential. Because of their ecological origin in environments that become quickly anoxic in the dark, we hypothesized that this differential ability may have evolved to serve a role in the fermentation of the photosynthate. Here we show that, when forced to ferment internal substrate, these cyanobacteria display desirable characteristics of physiological H₂ production. Among them, the strain L. aestuarii BL J had the fastest specific rates and attained the highest H₂ concentrations during fermentation of photosynthate, which proceeded via a mixed acid fermentation pathway to yield acetate, ethanol, lactate, H₂, CO₂, and pyruvate. Contrary to expectations, the H₂ yield per mole of glucose was only average compared to that of other cyanobacteria. Thermodynamic analyses point to the use of electron donors more electronegative than NAD(P)H in Lyngbya hydrogenases as the basis for its strong H₂ production ability. In any event, the high specific rates and H₂ concentrations coupled with the lack of reversibility of the enzyme, at the expense of internal, photosynthetically generated reductants, makes L. aestuarii BL J and/or its enzymes, a potentially feasible platform for large-scale H₂ production.

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.

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.

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.

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.