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.

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.

The ongoing Zika virus (ZIKV) epidemic in the Americas poses a major global public health emergency. While ZIKV is transmitted from human to human by bites of Aedes mosquitoes, recent evidence indicates that ZIKV can also be transmitted via sexual contact with cases of sexually transmitted ZIKV reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA. Yet, the role of sexual transmission on the spread and control of ZIKV infection is not well-understood. We introduce a mathematical model to investigate the impact of mosquito-borne and sexual transmission on the spread and control of ZIKV and calibrate the model to ZIKV epidemic data from Brazil, Colombia, and El Salvador. Parameter estimates yielded a basic reproduction number R_0 = 2.055 (95% CI: 0.523–6.300), in which the percentage contribution of sexual transmission is 3.044% (95% CI: 0.123–45.73). Our sensitivity analyses indicate that R_0 is most sensitive to the biting rate and mortality rate of mosquitoes while sexual transmission increases the risk of infection and epidemic size and prolongs the outbreak. Prevention and control efforts against ZIKV should target both the mosquito-borne and sexual transmission routes.

Our ability to uncover complex network structure and dynamics from data is fundamental to understanding and controlling collective dynamics in complex systems. Despite recent progress in this area, reconstructing networks with stochastic dynamical processes from limited time series remains to be an outstanding problem. Here we develop a framework based on compressed sensing to reconstruct complex networks on which stochastic spreading dynamics take place. We apply the methodology to a large number of model and real networks, finding that a full reconstruction of inhomogeneous interactions can be achieved from small amounts of polarized (binary) data, a virtue of compressed sensing. Further, we demonstrate that a hidden source that triggers the spreading process but is externally inaccessible can be ascertained and located with high confidence in the absence of direct routes of propagation from it. Our approach thus establishes a paradigm for tracing and controlling epidemic invasion and information diffusion in complex networked systems.

Background: Hemorrhagic fever with renal syndrome (HFRS), a rodent-borne infectious disease, is one of the most serious public health threats in China. Increasing our understanding of the spatial and temporal patterns of HFRS infections could guide local prevention and control strategies.

Methodology/Principal Findings: We employed statistical models to analyze HFRS case data together with environmental data from the Dongting Lake district during 2005–2010. Specifically, time-specific ecologic niche models (ENMs) were used to quantify and identify risk factors associated with HFRS transmission as well as forecast seasonal variation in risk across geographic areas. Results showed that the Maximum Entropy model provided the best predictive ability (AUC = 0.755). Time-specific Maximum Entropy models showed that the potential risk areas of HFRS significantly varied across seasons. High-risk areas were mainly found in the southeastern and southwestern areas of the Dongting Lake district. Our findings based on models focused on the spring and winter seasons showed particularly good performance. The potential risk areas were smaller in March, May and August compared with those identified for June, July and October to December. Both normalized difference vegetation index (NDVI) and land use types were found to be the dominant risk factors.

Conclusions/Significance: Our findings indicate that time-specific ENMs provide a useful tool to forecast the spatial and temporal risk of HFRS.

Background:
Pandemic influenza is said to 'shift mortality' to younger age groups; but also to spare a subpopulation of the elderly population. Does one of these effects dominate? Might this have important ramifications?

Methods: We estimated age-specific excess mortality rates for all-years for which data were available in the 20th century for Australia, Canada, France, Japan, the UK, and the USA for people older than 44 years of age. We modeled variation with age, and standardized estimates to allow direct comparison across age groups and countries. Attack rate data for four pandemics were assembled.

Results: For nearly all seasons, an exponential model characterized mortality data extremely well. For seasons of emergence and a variable number of seasons following, however, a subpopulation above a threshold age invariably enjoyed reduced mortality. 'Immune escape', a stepwise increase in mortality among the oldest elderly, was observed a number of seasons after both the A(H2N2) and A(H3N2) pandemics. The number of seasons from emergence to escape varied by country. For the latter pandemic, mortality rates in four countries increased for younger age groups but only in the season following that of emergence. Adaptation to both emergent viruses was apparent as a progressive decrease in mortality rates, which, with two exceptions, was seen only in younger age groups. Pandemic attack rate variation with age was estimated to be similar across four pandemics with very different mortality impact.

Conclusions: In all influenza pandemics of the 20th century, emergent viruses resembled those that had circulated previously within the lifespan of then-living people. Such individuals were relatively immune to the emergent strain, but this immunity waned with mutation of the emergent virus. An immune subpopulation complicates and may invalidate vaccine trials. Pandemic influenza does not 'shift' mortality to younger age groups; rather, the mortality level is reset by the virulence of the emerging virus and is moderated by immunity of past experience. In this study, we found that after immune escape, older age groups showed no further mortality reduction, despite their being the principal target of conventional influenza vaccines. Vaccines incorporating variants of pandemic viruses seem to provide little benefit to those previously immune. If attack rates truly are similar across pandemics, it must be the case that immunity to the pandemic virus does not prevent infection, but only mitigates the consequences.

Methicillin resistant Staphylococcus aureus (MRSA) is currently a major cause of skin and soft tissue infections (SSTI) in the United States. Seasonal variation of MRSA infections in hospital settings has been widely observed. However, systematic time-series analysis of incidence data is desirable to understand the seasonality of community acquired (CA)-MRSA infections at the population level. In this paper, using data on monthly SSTI incidence in children aged 0–19 years and enrolled in Medicaid in Maricopa County, Arizona, from January 2005 to December 2008, we carried out time-series and nonlinear regression analysis to determine the periodicity, trend, and peak timing in SSTI incidence in children at different age: 0-4 years, 5-9 years, 10-14 years, and 15-19 years. We also assessed the temporal correlation between SSTI incidence and meteorological variables including average temperature and humidity. Our analysis revealed a strong annual seasonal pattern of SSTI incidence with peak occurring in early September. This pattern was consistent across age groups. Moreover, SSTIs followed a significantly increasing trend over the 4-year study period with annual incidence increasing from 3.36% to 5.55% in our pediatric population of approximately 290,000. We also found a significant correlation between the temporal variation in SSTI incidence and mean temperature and specific humidity. Our findings could have potential implications on prevention and control efforts against CA-MRSA.

Background: Prior research shows that work in agriculture and construction/extraction occupations increases the risk of environmental heat-associated death.

Purpose: To assess the risk of environmental heat-associated death by occupation.

Methods: This was a case-control study. Cases were heat-caused and heat-related deaths occurring from May-October during the period 2002–2009 in Maricopa County, Arizona. Controls were selected at random from non-heat-associated deaths during the same period in Maricopa County. Information on occupation, age, sex, and race-ethnicity was obtained from death certificates. Logistic regression analysis was used to estimate odds ratios for heat-associated death.

Results: There were 444 cases of heat-associated deaths in adults (18+ years) and 925 adult controls. Of heat-associated deaths, 332 (75%) occurred in men; a construction/extraction or agriculture occupation was described on the death certificate in 115 (35%) of these men. In men, the age-adjusted odds ratios for heat-associated death were 2.32 (95% confidence interval 1.55, 3.48) in association with construction/extraction and 3.50 (95% confidence interval 1.94, 6.32) in association with agriculture occupations. The odds ratio for heat-associated death was 10.17 (95% confidence interval 5.38, 19.23) in men with unknown occupation. In women, the age-adjusted odds ratio for heat-associated death was 6.32 (95% confidence interval 1.48, 27.08) in association with unknown occupation. Men age 65 years and older in agriculture occupations were at especially high risk of heat-associated death.

Conclusion: The occurrence of environmental heat-associated death in men in agriculture and construction/extraction occupations in a setting with predictable periods of high summer temperatures presents opportunities for prevention.

In this study we characterized the relationship between temperature and mortality in central Arizona desert cities that have an extremely hot climate. Relationships between daily maximum apparent temperature (AT_max) and mortality for eight condition-specific causes and all-cause deaths were modeled for all residents and separately for males and females ages <65 and ≥65 during the months May-October for years 2000-2008. The most robust relationship was between ATmax on day of death and mortality from direct exposure to high environmental heat. For this condition-specific cause of death, the heat thresholds in all gender and age groups (AT_max = 90–97 °F; 32.2‒36.1 °C) were below local median seasonal temperatures in the study period (AT_max = 99.5 °F; 37.5 °C). Heat threshold was defined as AT_max at which the mortality ratio begins an exponential upward trend. Thresholds were identified in younger and older females for cardiac disease/stroke mortality (AT_max = 106 and 108 °F; 41.1 and 42.2 °C) with a one-day lag. Thresholds were also identified for mortality from respiratory diseases in older people (AT_max = 109 °F; 42.8 °C) and for all-cause mortality in females (AT_max = 107 °F; 41.7 °C) and males <65 years (AT_max = 102 °F; 38.9 °C). Heat-related mortality in a region that has already made some adaptations to predictable periods of extremely high temperatures suggests that more extensive and targeted heat-adaptation plans for climate change are needed in cities worldwide.