Local control: a hub-based model for the c-di-GMP network

The genome of Pseudomonas fluorescens encodes >50 proteins predicted to play a role in bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP)-mediated biofilm formation. We built a network representation of protein-protein interactions and extracted key information via multidimensional scaling (i.e., principal component analysis) of node centrality measures, which measure features of proteins in a network. Proteins of different domain types (diguanylate cyclase, dual domain, phosphodiesterase, PilZ) exhibit unique network behavior and can be accurately classified by their network centrality values (i.e., roles in the network). The predictive power of protein-protein interactions in biofilm formation indicates the possibility of localized pools of c-di-GMP. A regression model showed a statistically significant impact of protein-protein interactions on the extent of biofilm formation in various environments. These results highlight the importance of a localized c-di-GMP signaling, extend our understanding of signaling by this second messenger beyond the current "Bow-tie Model," support a newly proposed "Hub Model," and suggest future avenues of investigation.

Sequential stacking link prediction algorithms for temporal networks

Link prediction algorithms are indispensable tools in many scientific applications by speeding up network data collection and imputing missing connections. However, in many systems, links change over time and it remains unclear how to optimally exploit such temporal information for link predictions in such networks. Here, we show that many temporal topological features, in addition to having high computational cost, are less accurate in temporal link prediction than sequentially stacked static network features. This sequential stacking link prediction method uses 41 static network features that avoid detailed feature engineering choices and is capable of learning a highly accurate predictive distribution of future connections from historical data. We demonstrate that this algorithm works well for both partially observed and completely unobserved target layers, and on two temporal stochastic block models achieves near-oracle-level performance when combined with other single predictor methods as an ensemble learning method. Finally, we empirically illustrate that stacking multiple predictive methods together further improves performance on 19 real-world temporal networks from different domains.

Frontoparietal functional connectivity moderates the link between time spent on social media and subsequent negative affect in daily life

Evidence on the harms and benefits of social media use is mixed, in part because the effects of social media on well-being depend on a variety of individual difference moderators. Here, we explored potential neural moderators of the link between time spent on social media and subsequent negative affect. We specifically focused on the strength of correlation among brain regions within the frontoparietal system, previously associated with the top-down cognitive control of attention and emotion. Participants (N = 54) underwent a resting state functional magnetic resonance imaging scan. Participants then completed 28 days of ecological momentary assessment and answered questions about social media use and negative affect, twice a day. Participants who spent more than their typical amount of time on social media since the previous time point reported feeling more negative at the present moment. This within-person temporal association between social media use and negative affect was mainly driven by individuals with lower resting state functional connectivity within the frontoparietal system. By contrast, time spent on social media did not predict subsequent affect for individuals with higher frontoparietal functional connectivity. Our results highlight the moderating role of individual functional neural connectivity in the relationship between social media and affect.

Psychological distance intervention reminders reduce alcohol consumption frequency in daily life

Modifying behaviors, such as alcohol consumption, is difficult. Creating psychological distance between unhealthy triggers and one's present experience can encourage change. Using two multisite, randomized experiments, we examine whether theory-driven strategies to create psychological distance-mindfulness and perspective-taking-can change drinking behaviors among young adults without alcohol dependence via a 28-day smartphone intervention (Study 1, N = 108 participants, 5492 observations; Study 2, N = 218 participants, 9994 observations). Study 2 presents a close replication with a fully remote delivery during the COVID-19 pandemic. During weeks when they received twice-a-day intervention reminders, individuals in the distancing interventions reported drinking less frequently than on control weeks-directionally in Study 1, and significantly in Study 2. Intervention reminders reduced drinking frequency but did not impact amount. We find that smartphone-based mindfulness and perspective-taking interventions, aimed to create psychological distance, can change behavior. This approach requires repeated reminders, which can be delivered via smartphones.

Dense Sampling Approaches for Psychiatry Research: Combining Scanners and Smartphones

Together, data from brain scanners and smartphones have sufficient coverage of biology, psychology, and environment to articulate between-person differences in the interplay within and across biological, psychological, and environmental systems thought to underlie psychopathology. An important next step is to develop frameworks that combine these two modalities in ways that leverage their coverage across layers of human experience to have maximum impact on our understanding and treatment of psychopathology. We review literature published in the last 3 years highlighting how scanners and smartphones have been combined to date, outline and discuss the strengths and weaknesses of existing approaches, and sketch a network science framework heretofore underrepresented in work combining scanners and smartphones that can push forward our understanding of health and disease.

Age-related changes in network controllability are mitigated by redundancy in large-scale brain networks

The aging brain undergoes major changes in its topology. The mechanisms by which the brain mitigates age-associated changes in topology to maintain robust control of brain networks are unknown. Here we used diffusion MRI data from cognitively intact participants (n=480, ages 40-90) to study age-associated changes in the controllability of structural brain networks, features that could mitigate these changes, and the overall effect on cognitive function. We found age-associated declines in controllability in control hubs and large-scale networks, particularly within the and frontoparietal control and default mode networks. Redundancy, quantified via the assessment of multi-step paths within networks, mitigated the effects of changes in topology on network controllability. Lastly, network controllability, redundancy, and grey matter volume each played important complementary roles in cognitive function. In sum, our results highlight the importance of redundancy for robust control of brain networks and in cognitive function in healthy-aging.

Within-Person Temporal Associations Among Self-Reported Physical Activity, Sleep, and Well-Being in College Students

OBJECTIVE

A holistic understanding of the naturalistic dynamics among physical activity, sleep, emotions, and purpose in life as part of a system reflecting wellness is key to promoting well-being. The main aim of this study was to examine the day-to-day dynamics within this wellness system.

METHODS

Using self-reported emotions (happiness, sadness, anger, anxiousness) and physical activity periods collected twice per day, and daily reports of sleep and purpose in life via smartphone experience sampling, more than 28 days as college students ( n = 226 young adults; mean [standard deviation] = 20.2 [1.7] years) went about their daily lives, we examined day-to-day temporal and contemporaneous dynamics using multilevel vector autoregressive models that consider the network of wellness together.

RESULTS

Network analyses revealed that higher physical activity on a given day predicted an increase of happiness the next day. Higher sleep quality on a given night predicted a decrease in negative emotions the next day, and higher purpose in life predicted decreased negative emotions up to 2 days later. Nodes with the highest centrality were sadness, anxiety, and happiness in the temporal network and purpose in life, anxiety, and anger in the contemporaneous network.

CONCLUSIONS

Although the effects of sleep and physical activity on emotions and purpose in life may be shorter term, a sense of purpose in life is a critical component of wellness that can have slightly longer effects, bleeding into the next few days. High-arousal emotions and purpose in life are central to motivating people into action, which can lead to behavior change.

E. coli catheter-associated urinary tract infections are associated with distinctive virulence and biofilm gene determinants

Urinary catheterization facilitates urinary tract colonization by E. coli and increases infection risk. Here, we aimed to identify strain-specific characteristics associated with the transition from colonization to infection in catheterized patients. In a single-site study population, we compared E. coli isolates from patients with catheter-associated asymptomatic bacteriuria (CAASB) to those with catheter-associated urinary tract infection (CAUTI). CAUTI isolates were dominated by a phylotype B2 subclade containing the multidrug-resistant ST131 lineage relative to CAASB isolates, which were phylogenetically more diverse. A distinctive combination of virulence-associated genes was present in the CAUTI-associated B2 subclade. Catheter-associated biofilm formation was widespread among isolates and did not distinguish CAUTI from CAASB strains. Preincubation with CAASB strains could inhibit catheter colonization by multiple ST131 CAUTI isolates. Comparative genomic analysis identified a group of variable genes associated with high catheter biofilm formation present in both CAUTI and CAASB strains. Among these, ferric citrate transport (Fec) system genes were experimentally associated with enhanced catheter biofilm formation using reporter and fecA deletion strains. These results are consistent with a variable role for catheter biofilm formation in promoting CAUTI by ST131-like strains or resisting CAUTI by lower-risk strains that engage in niche exclusion.

Purpose in life, neural alcohol cue reactivity and daily alcohol use in social drinkers

BACKGROUND AND AIM

Alcohol craving is an urge to consume alcohol that commonly precedes drinking; however, craving does not lead to drinking for all people under all circumstances. The current study measured the correlation between neural reactivity and alcohol cues as a risk, and purpose in daily life as a protective factor that may influence the link between alcohol craving and the subsequent amount of consumption.

DESIGN

Observational study that correlated functional magnetic resonance imaging (fMRI) data on neural cue reactivity and ecological momentary assessments (EMA) on purpose in life and alcohol use.

SETTING

Two college campuses in the United States.

PARTICIPANTS

A total of 54 college students (37 women, 16 men, and 1 other) recruited via campus-based groups from January 2019 to October 2020.

MEASUREMENTS

Participants underwent fMRI while viewing images of alcohol; we examined activity within the ventral striatum, a key region of interest implicated in reward and craving. Participants then completed 28 days of EMA and answered questions about daily levels of purpose in life and alcohol use, including how much they craved and consumed alcohol.

FINDINGS

A significant three-way interaction indicated that greater alcohol cue reactivity within the ventral striatum was associated with heavier alcohol use following craving in daily life only when people were previously feeling a lower than usual sense of purpose. By contrast, individuals with heightened neural alcohol cue reactivity drank less in response to craving if they were feeling a stronger than their usual sense of purpose in the preceding moments (b_interaction  = -0.086, P < 0.001, 95% CI = -0.137, -0.035).

CONCLUSIONS

Neural sensitivity to alcohol cues within the ventral striatum appears to be a potential risk for increased alcohol use in social drinkers, when people feel less purposeful. Enhancing daily levels of purpose in life may promote alcohol moderation among social drinkers who show relatively higher reactivity to alcohol cues.

Altered neural flexibility in children with attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is one of the most common neurodevelopmental disorders of childhood, and is often characterized by altered executive functioning. Executive function has been found to be supported by flexibility in dynamic brain reconfiguration. Thus, we applied multilayer community detection to resting-state fMRI data in 180 children with ADHD and 180 typically developing children (TDC) to identify alterations in dynamic brain reconfiguration in children with ADHD. We specifically evaluated MR derived neural flexibility, which is thought to underlie cognitive flexibility, or the ability to selectively switch between mental processes. Significantly decreased neural flexibility was observed in the ADHD group at both the whole brain (raw p = 0.0005) and sub-network levels (p < 0.05, FDR corrected), particularly for the default mode network, attention-related networks, executive function-related networks, and primary networks. Furthermore, the subjects with ADHD who received medication exhibited significantly increased neural flexibility (p = 0.025, FDR corrected) when compared to subjects with ADHD who were medication naïve, and their neural flexibility was not statistically different from the TDC group (p = 0.74, FDR corrected). Finally, regional neural flexibility was capable of differentiating ADHD from TDC (Accuracy: 77% for tenfold cross-validation, 74.46% for independent test) and of predicting ADHD severity using clinical measures of symptom severity (R2: 0.2794 for tenfold cross-validation, 0.156 for independent test). In conclusion, the present study found that neural flexibility is altered in children with ADHD and demonstrated the potential clinical utility of neural flexibility to identify children with ADHD, as well as to monitor treatment responses and disease severity.

A bipartite graph-based expected networks approach identifies DDR genes not associated with TMB yet predictive of immune checkpoint blockade response

Immune checkpoint blockade (ICB) has had remarkable success for treatment of solid tumors. However, as only a subset of patients exhibit responses, there is a continued need for biomarker development. Numerous reports have shown a link between tumor mutational burden (TMB) and ICB response, while others have identified a link between ICB response and mutation in DNA damage repair (DDR) genes. However, it remains unclear to what extent mutations in DDR genes hold predictive value above and beyond their association with TMB. Herein, we present a networks-based test and bipartite graph-based expected TMB score (BiG-BETS) with higher specificity for discriminating DDR genes and pathways that are associated with elevated TMB. Moreover, we find that mutations in certain DDR genes that are not associated with elevated TMB (low BiG-BETS) are nevertheless predictive of ICB benefit in high TMB patients, demonstrating that their inactivation contributes to ICB response in a TMB-independent manner.

A Simultaneous Feature Selection and Compositional Association Test for Detecting Sparse Associations in High-Dimensional Metagenomic Data

Numerous metagenomic studies aim to discover associations between the microbial composition of an environment (e.g., gut, skin, oral) and a phenotype of interest. Multivariate analysis is often performed in these studies without critical a priori knowledge of which taxa are associated with the phenotype being studied. This approach typically reduces statistical power in settings where the true associations among only a few taxa are obscured by high dimensionality (i.e., sparse association signals). At the same time, low sample size and compositional sample space constraints may reduce beyond-study generalizability if not properly accounted for. To address these difficulties, we developed the Selection-Energy-Permutation (SelEnergyPerm) method, a nonparametric group association test with embedded feature selection that directly accounts for compositional constraints using parsimonious logratio signatures between taxonomic features, for characterizing and understanding alterations in microbial community structure. Simulation results show SelEnergyPerm selects small independent sets of logratios that capture strong associations in a range of scenarios. Additionally, our simulation results demonstrate SelEnergyPerm consistently detects/rejects associations in synthetic data with sparse, dense, or no association signals. We demonstrate the novel benefits of our method in four case studies utilizing publicly available 16S amplicon and whole-genome sequencing datasets. Our R implementation of Selection-Energy-Permutation, including an example demonstration and the code to generate all of the scenarios used here, is available at https://www.github.com/andrew84830813/selEnergyPermR.

Comparing transmission potential networks based on social network surveys, close contacts and environmental overlap in rural Madagascar

Social and spatial network analysis is an important approach for investigating infectious disease transmission, especially for pathogens transmitted directly between individuals or via environmental reservoirs. Given the diversity of ways to construct networks, however, it remains unclear how well networks constructed from different data types effectively capture transmission potential. We used empirical networks from a population in rural Madagascar to compare social network survey and spatial data-based networks of the same individuals. Close contact and environmental pathogen transmission pathways were modelled with the spatial data. We found that naming social partners during the surveys predicted higher close-contact rates and the proportion of environmental overlap on the spatial data-based networks. The spatial networks captured many strong and weak connections that were missed using social network surveys alone. Across networks, we found weak correlations among centrality measures (a proxy for superspreading potential). We conclude that social network surveys provide important scaffolding for understanding disease transmission pathways but miss contact-specific heterogeneities revealed by spatial data. Our analyses also highlight that the superspreading potential of individuals may vary across transmission modes. We provide detailed methods to construct networks for close-contact transmission pathogens when not all individuals simultaneously wear GPS trackers.

The association between hippocampal volume and memory in pathological aging is mediated by functional redundancy

Hippocampal neurodegeneration, a primary component of Alzheimer's disease pathology, relates to poor cognition; however, the mechanisms underlying this relationship are not well understood. Using a sample of cognitively normal older adults and individuals with mild cognitive impairment, this study aims to determine the topological properties of functional networks accompanying hippocampal atrophy in aging, along with their association to cognition and clinical progression. We considered two conceptually differing topological properties: redundancy (the existence of alternative channels of functional commutation) and local efficiency (the efficiency of local information exchange). Hippocampal redundancy, but not local efficiency, mediated the association between low hippocampal volume and low memory in both the whole sample and in ß-amyloid positive participants. Additionally, participants with high hippocampal volume, redundancy, and memory clustered separately from those with low values on all three measures, with the latter group showing higher conversion rates to dementia within three years. Together, these results demonstrate that reduced hippocampal redundancy is one mechanism through which hippocampal atrophy associates with memory impairment in healthy and pathological aging.

Synchronization of coupled Kuramoto oscillators under resource constraints

A fundamental understanding of synchronized behavior in multiagent systems can be acquired by studying analytically tractable Kuramoto models. However, such models typically diverge from many real systems whose dynamics evolve under nonnegligible resource constraints. Here we construct a system of coupled Kuramoto oscillators that consume or produce resources as a function of their oscillation frequency. At high coupling, we observe strongly synchronized dynamics, whereas at low coupling, we observe independent oscillator dynamics as expected from standard Kuramoto models. For intermediate coupling, which typically induces a partially synchronized state, we empirically observe that (and theoretically explain why) the system can exist in either: (i) a state in which the order parameter oscillates in time, or (ii) a state in which multiple synchronization states are simultaneously stable. Whether (i) or (ii) occurs depends upon whether the oscillators consume or produce resources, respectively. Relevant for systems as varied as coupled neurons and social groups, our paper lays important groundwork for future efforts to develop quantitative predictions of synchronized dynamics for systems embedded in environments marked by sparse resources.

Optimizing Emergency Stroke Transport Strategies Using Physiological Models

Supplemental Digital Content is available in the text.

The criteria for choosing between drip and ship and mothership transport strategies in emergency stroke care is widely debated. Although existing data-driven probability models can inform transport decision-making at an epidemiological level, we propose a novel mathematical, physiologically derived framework that provides insight into how patient characteristics underlying infarct core growth influence these decisions.

Accrual of functional redundancy along the lifespan and its effects on cognition

Despite the necessity to understand how the brain endures the initial stages of age-associated cognitive decline, no brain mechanism has been quantitatively specified to date. The brain may withstand the effects of cognitive aging through redundancy, a design feature in engineered and biological systems, which entails the presence of substitute elements to protect it against failure. Here, we investigated the relationship between functional network redundancy and age over the human lifespan and their interaction with cognition, analyzing resting-state functional MRI images and cognitive measures from 579 subjects. Network-wide redundancy was significantly associated with age, showing a stronger link with age than other major topological measures, presenting a pattern of accumulation followed by old-age decline. Critically, redundancy significantly mediated the association between age and executive function, with lower anti-correlation between age and cognition in subjects with high redundancy. The results suggest that functional redundancy accrues throughout the lifespan, mitigating the effects of age on cognition.

Fecal IgA, Antigen Absorption, and Gut Microbiome Composition Are Associated With Food Antigen Sensitization in Genetically Susceptible Mice

Food allergy is a potentially fatal disease affecting 8% of children and has become increasingly common in the past two decades. Despite the prevalence and severe nature of the disease, the mechanisms underlying sensitization remain to be further elucidated. The Collaborative Cross is a genetically diverse panel of inbred mice that were specifically developed to study the influence of genetics on complex diseases. Using this panel of mouse strains, we previously demonstrated CC027/GeniUnc mice, but not C3H/HeJ mice, develop peanut allergy after oral exposure to peanut in the absence of a Th2-skewing adjuvant. Here, we investigated factors associated with sensitization in CC027/GeniUnc mice following oral exposure to peanut, walnut, milk, or egg. CC027/GeniUnc mice mounted antigen-specific IgE responses to peanut, walnut and egg, but not milk, while C3H/HeJ mice were not sensitized to any antigen. Naïve CC027/GeniUnc mice had markedly lower total fecal IgA compared to C3H/HeJ, which was accompanied by stark differences in gut microbiome composition. Sensitized CC027/GeniUnc mice had significantly fewer CD3+ T cells but higher numbers of CXCR5+ B cells and T follicular helper cells in the mesenteric lymph nodes compared to C3H/HeJ mice, which is consistent with their relative immunoglobulin production. After oral challenge to the corresponding food, peanut- and walnut-sensitized CC027/GeniUnc mice experienced anaphylaxis, whereas mice exposed to milk and egg did not. Ara h 2 was detected in serum collected post-challenge from peanut-sensitized mice, indicating increased absorption of this allergen, while Bos d 5 and Gal d 2 were not detected in mice exposed to milk and egg, respectively. Machine learning on the change in gut microbiome composition as a result of food protein exposure identified a unique signature in CC027/GeniUnc mice that experienced anaphylaxis, including the depletion of Akkermansia. Overall, these results demonstrate several factors associated with enteral sensitization in CC027/GeniUnc mice, including diminished total fecal IgA, increased allergen absorption and altered gut microbiome composition. Furthermore, peanuts and tree nuts may have inherent properties distinct from milk and eggs that contribute to allergy.

Lower functional hippocampal redundancy in mild cognitive impairment

With an increasing prevalence of mild cognitive impairment (MCI) and Alzheimer's disease (AD) in response to an aging population, it is critical to identify and understand neuroprotective mechanisms against cognitive decline. One potential mechanism is redundancy: the existence of duplicate elements within a system that provide alternative functionality in case of failure. As the hippocampus is one of the earliest sites affected by AD pathology, we hypothesized that functional hippocampal redundancy is protective against cognitive decline. We compared hippocampal functional redundancy derived from resting-state functional MRI networks in cognitively normal older adults, with individuals with early and late MCI, as well as the relationship between redundancy and cognition. Posterior hippocampal redundancy was reduced between cognitively normal and MCI groups, plateauing across early and late MCI. Higher hippocampal redundancy was related to better memory performance only for cognitively normal individuals. Critically, functional hippocampal redundancy did not come at the expense of network efficiency. Our results provide support that hippocampal redundancy protects against cognitive decline in aging.

Positive psychological states and stress responses in caregivers of adults receiving an allogeneic bone marrow transplant: A study protocol

AIMS

This protocol directs a study that aims to: (a) describe the caregiver's experience over 8-12 weeks after an index adult patient's allogeneic bone marrow transplant (BMT) for advanced cancer using a case-oriented approach and mixed methods, with qualitative methods in the foreground; and (b) explore networks of relationships among psycho-neurological symptoms, positive psychological states and caregiver health.

DESIGN

Case-oriented longitudinal design using multiple data types and analytic approaches.

METHODS

Data will be collected from 10-12 caregivers. The sample will be recruited from a large public hospital in the southeastern United States using maximum variation sampling (e.g., caregiver race/ethnicity, relationship to patient, age, education, and number of caregiving roles). Participants will be asked to complete weekly surveys, have their blood drawn bi-weekly and participate in an interview each month during the study period (~100 days). Aim 1 analysis will include directed content analysis and case-oriented visual analysis. Aim 2 analysis will include symptom network estimation of psycho-neurological symptoms, positive psychological states, and caregiver health. Institutional review board approval was obtained August 2018.

DISCUSSION

Results will provide an in-depth description of caregivers' experiences in the 100 days after BMT. Findings will inform generation of hypotheses and identification of targets for interventions to improve caregiver's experiences after BMT.

IMPACT

This in-depth multi-method longitudinal study to describe caregivers of adult patients receiving an allogeneic BMT is an essential step in understanding caregivers' complex responses to chronic stress and the role of positive psychological states. The results from this study will inform future research on chronic stress processes, intense caregiving, and intervention development.

Dynamics of social network emergence explain network evolution

Networked systems emerge and subsequently evolve. Although several models describe the process of network evolution, researchers know far less about the initial process of network emergence. Here, we report temporal survey results of a real-world social network starting from its point of inception. We find that individuals' ties undergo an initial cycle of rapid expansion and contraction. This process helps to explain the eventual interactions and working structure in the network (in this case, scientific collaboration). We propose a stylized concept and model of "churn" to describe the process of network emergence and stabilization. Our empirical and simulation results suggest that these network emergence dynamics may be instrumental for explaining network details, as well as behavioral outcomes at later time periods.

Abstract WP291: Stochastic Methods Can Resolve the Dilemma of Emergency Stroke Transport

Introduction: Drip and Ship (DNS) and Mothership (MS) are well-known emergency transport strategies in acute stroke care, but the criteria for choosing between the two is widely debated. Existing models define time-dependent outcomes but cannot resolve this debate with statistical significance because the independent variables are deterministic. We propose a novel stochastic framework that quantifies statistical significance between DNS and MS in a network of primary and comprehensive stroke centers.

Methods: We represented the physiology of ischemic core growth as a stochastic first-order differential equation, enabling infarct volume at time of reperfusion to be calculated and mapped to 90-day mRS. Using Texas as a case study, we configured the state’s stroke network within 15,811 geographic blocks as defined by census data. For each block, we ran Monte Carlo simulations to generate Beta distributions of large- and small-vessel infarct volumes, which were then translated into cumulative distribution functions of mRS. A two-sample Kolmogorov-Smirnov test for significance, and Cohen’s d effect size statistic for practical significance were computed between each DNS and MS pair. Stable effect sizes were assured by sampling > 5,000 total infarct volumes for each block. All model parameters were established from large cohort studies or trials.

Results: Of the 13,113 blocks where the primary stroke center is the closest hospital from origin, DNS produces significantly better stroke outcomes than MS in 79.0% (0.3% SEM; P < 0.05; 0.2 < d < 0.5). For the subset of patients with large-vessel strokes, MS produces significantly better outcomes in 44.6% of blocks (1.3% SEM; P < 0.05; 0.4 < d < 0.85).

Conclusion: Stochastic methods enable the use of clinically relevant metrics for comparative significance of DNS and MS in a geographic region. This formalism, which has not been incorporated in previous models, can be further generalized beyond stochastic infarct volumes if sufficiently large datasets become available. For example, the kinetic growth model can integrate the statistical distributions of times (pre-hospital and hospital) leading up to intervention, and patient attributes that affect outcomes, such as the degree of collateral flow and comorbidities.

Emotion semantics show both cultural variation and universal structure

Many human languages have words for emotions such as “anger” and “fear,” yet it is not clear whether these emotions have similar meanings across languages, or why their meanings might vary. We estimate emotion semantics across a sample of 2474 spoken languages using “colexification”—a phenomenon in which languages name semantically related concepts with the same word. Analyses show significant variation in networks of emotion concept colexification, which is predicted by the geographic proximity of language families. We also find evidence of universal structure in emotion colexification networks, with all families differentiating emotions primarily on the basis of hedonic valence and physiological activation. Our findings contribute to debates about universality and diversity in how humans understand and experience emotion.

Concurrency and reachability in treelike temporal networks

Network properties govern the rate and extent of various spreading processes, from simple contagions to complex cascades. Recently, the analysis of spreading processes has been extended from static networks to temporal networks, where nodes and links appear and disappear. We focus on the effects of accessibility, whether there is a temporally consistent path from one node to another, and reachability, the density of the corresponding accessibility graph representation of the temporal network. The level of reachability thus inherently limits the possible extent of any spreading process on the temporal network. We study reachability in terms of the overall levels of temporal concurrency between edges and the structural cohesion of the network agglomerating over all edges. We use simulation results and develop heterogeneous mean-field model predictions for random networks to better quantify how the properties of the underlying temporal network regulate reachability.

Map equation with metadata: Varying the role of attributes in community detection

Much of the community detection literature studies structural communities, communities defined solely by the connectivity patterns of the network. Often networks contain additional metadata which can inform community detection such as the grade and gender of students in a high school social network. In this work, we introduce a tuning parameter to the content map equation that allows users of the Infomap community detection algorithm to control the metadata's relative importance for identifying network structure. On synthetic networks, we show that our algorithm can overcome the structural detectability limit when the metadata are well aligned with community structure. On real-world networks, we show how our algorithm can achieve greater mutual information with the metadata at a cost in the traditional map equation. Our tuning parameter, like the focusing knob of a microscope, allows users to "zoom in" and "zoom out" on communities with varying levels of focus on the metadata.

Comparing Shared Patient Networks Across Payers

BACKGROUND

Measuring care coordination in administrative data facilitates important research to improve care quality.

OBJECTIVE

To compare shared patient networks constructed from administrative claims data across multiple payers.

DESIGN

Social network analysis of pooled cross sections of physicians treating prevalent colorectal cancer patients between 2003 and 2013.

PARTICIPANTS

Surgeons, medical oncologists, and radiation oncologists identified from North Carolina Central Cancer Registry data linked to Medicare claims (N = 1735) and private insurance claims (N = 1321).

MAIN MEASURES

Provider-level measures included the number of patients treated, the number of providers with whom they share patients (by specialty), the extent of patient sharing with each specialty, and network centrality. Network-level measures included the number of providers and shared patients, the density of shared-patient relationships among providers, and the size and composition of clusters of providers with a high level of patient sharing.

RESULTS

For 24.5% of providers, total patient volume rank differed by at least one quintile group between payers. Medicare claims missed 14.6% of all shared patient relationships between providers, but captured a greater number of patient-sharing relationships per provider compared with the private insurance database, even after controlling for the total number of patients (27.242 vs 26.044, p < 0.001). Providers in the private network shared a higher fraction of patients with other providers (0.226 vs 0.127, p < 0.001) compared to the Medicare network. Clustering coefficients for providers, weighted betweenness, and eigenvector centrality varied greatly across payers. Network differences led to some clusters of providers that existed in the combined network not being detected in Medicare alone.

CONCLUSION

Many features of shared patient networks constructed from a single-payer database differed from similar networks constructed from other payers' data. Depending on a study's goals, shortcomings of single-payer networks should be considered when using claims data to draw conclusions about provider behavior.

EndNote: Feature-based classification of networks

Network representations of systems from various scientific and societal domains are neither completely random nor fully regular, but instead appear to contain recurring structural features. These features tend to be shared by networks belonging to the same broad class, such as the class of social networks or the class of biological networks. Within each such class, networks describing similar systems tend to have similar features. This occurs presumably because networks representing similar systems would be expected to be generated by a shared set of domain specific mechanisms, and it should therefore be possible to classify networks based on their features at various structural levels. Here we describe and demonstrate a new hybrid approach that combines manual selection of network features of potential interest with existing automated classification methods. In particular, selecting well-known network features that have been studied extensively in social network analysis and network science literature, and then classifying networks on the basis of these features using methods such as random forest, which is known to handle the type of feature collinearity that arises in this setting, we find that our approach is able to achieve both higher accuracy and greater interpretability in shorter computation time than other methods.

Metabolomic networks connect host-microbiome processes to human Clostridioides difficile infections

Clostridioides difficile infection (CDI) accounts for a substantial proportion of deaths attributable to antibiotic-resistant bacteria in the United States. Although C. difficile can be an asymptomatic colonizer, its pathogenic potential is most commonly manifested in patients with antibiotic-modified intestinal microbiomes. In a cohort of 186 hospitalized patients, we showed that host and microbe-associated shifts in fecal metabolomes had the potential to distinguish patients with CDI from those with non-C. difficile diarrhea and C. difficile colonization. Patients with CDI exhibited a chemical signature of Stickland amino acid fermentation that was distinct from those of uncolonized controls. This signature suggested that C. difficile preferentially catabolizes branched chain amino acids during CDI. Unexpectedly, we also identified a series of noncanonical, unsaturated bile acids that were depleted in patients with CDI. These bile acids may derive from an extended host-microbiome dehydroxylation network in uninfected patients. Bile acid composition and leucine fermentation defined a prototype metabolomic model with potential to distinguish clinical CDI from asymptomatic C. difficile colonization.

Social clustering in epidemic spread on coevolving networks

Even though transitivity is a central structural feature of social networks, its influence on epidemic spread on coevolving networks has remained relatively unexplored. Here we introduce and study an adaptive susceptible-infected-susceptible (SIS) epidemic model wherein the infection and network coevolve with nontrivial probability to close triangles during edge rewiring, leading to substantial reinforcement of network transitivity. This model provides an opportunity to study the role of transitivity in altering the SIS dynamics on a coevolving network. Using numerical simulations and approximate master equations (AMEs), we identify and examine a rich set of dynamical features in the model. In many cases, AMEs including transitivity reinforcement provide accurate predictions of stationary-state disease prevalence and network degree distributions. Furthermore, for some parameter settings, the AMEs accurately trace the temporal evolution of the system. We show that higher transitivity reinforcement in the model leads to lower levels of infective individuals in the population, when closing a triangle is the dominant rewiring mechanism. These methods and results may be useful in developing ideas and modeling strategies for controlling SIS-type epidemics.

Assessing the robustness of cluster solutions obtained from sparse count matrices

Psychological researchers often seek to obtain cluster solutions from sparse count matrices (e.g., social networks; counts of symptoms that are in common for 2 given individuals; structural brain imaging). Increasingly, community detection methods are being used to subset the data in a data-driven manner. While many of these approaches perform well in simulation studies and thus offer some improvement upon traditional clustering approaches, there is no readily available approach for evaluating the robustness of these solutions in empirical data. Researchers have no way of knowing if their results are due to noise. We describe here 2 approaches novel to the field of psychology that enable evaluation of cluster solution robustness. This tutorial also explains the use of an associated R package, perturbR, which provides researchers with the ability to use the methods described herein. In the first approach, the cluster assignment from the original matrix is compared against cluster assignments obtained by randomly perturbing the edges in the matrix. Stable cluster solutions should not demonstrate large changes in the presence of small perturbations. For the second approach, Monte Carlo simulations of random matrices that have the same properties as the original matrix are generated. The distribution of quality scores ("modularity") obtained from the cluster solutions from these matrices are then compared with the score obtained from the original matrix results. From this, one can assess if the results are better than what would be expected by chance. perturbR automates these 2 methods, providing an easy-to-use resource for psychological researchers. We demonstrate the utility of this package using benchmark simulated data generated from a previous study and then apply the methods to publicly available empirical data obtained from social networks and structural neuroimaging. (PsycINFO Database Record (c) 2019 APA, all rights reserved).

Epidemic spreading in localized environments with recurrent mobility patterns

The spreading of epidemics is very much determined by the structure of the contact network, which may be impacted by the mobility dynamics of the individuals themselves. In confined scenarios where a small, closed population spends most of its time in localized environments and has easily identifiable mobility patterns—such as workplaces, university campuses, or schools—it is of critical importance to identify the factors controlling the rate of disease spread. Here, we present a discrete-time, metapopulation-based model to describe the transmission of susceptible-infected-susceptible-like diseases that take place in confined scenarios where the mobilities of the individuals are not random but, rather, follow clear recurrent travel patterns. This model allows analytical determination of the onset of epidemics, as well as the ability to discern which contact structures are most suited to prevent the infection to spread. It thereby determines whether common prevention mechanisms, as isolation, are worth implementing in such a scenario and their expected impact.

Functional classification of long non-coding RNAs by k-mer content

The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called kmers. We found that lncRNAs of related function often had similar kmer profiles despite lacking linear homology, and that kmer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. Kmer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.

RIGID GRAPH COMPRESSION: MOTIF-BASED RIGIDITY ANALYSIS FOR DISORDERED FIBER NETWORKS

Using particle-scale models to accurately describe property enhancements and phase transitions in macroscopic behavior is a major engineering challenge in composite materials science. To address some of these challenges, we use the graph theoretic property of rigidity to model mechanical reinforcement in composites with stiff rod-like particles. We develop an efficient algorithmic approach called rigid graph compression (RGC) to describe the transition from floppy to rigid in disordered fiber networks ("rod-hinge systems"), which form the reinforcing phase in many composite systems. To establish RGC on a firm theoretical foundation, we adapt rigidity matroid theory to identify primitive topological network motifs that serve as rules for composing interacting rigid particles into larger rigid components. This approach is computationally efficient and stable, because RGC requires only topological information about rod interactions (encoded by a sparse unweighted network) rather than geometrical details such as rod locations or pairwise distances (as required in rigidity matroid theory). We conduct numerical experiments on simulated two-dimensional rod-hinge systems to demonstrate that RGC closely approximates the rigidity percolation threshold for such systems, through comparison with the pebble game algorithm (which is exact in two dimensions). Importantly, whereas the pebble game is derived from Laman's condition and is only valid in two dimensions, the RGC approach naturally extends to higher dimensions.

NETWORK-ENSEMBLE COMPARISONS WITH STOCHASTIC REWIRING AND VON NEUMANN ENTROPY

Assessing whether a given network is typical or atypical for a random-network ensemble (i.e., network-ensemble comparison) has widespread applications ranging from null-model selection and hypothesis testing to clustering and classifying networks. We develop a framework for network-ensemble comparison by subjecting the network to stochastic rewiring. We study two rewiring processes-uniform and degree-preserved rewiring-which yield random-network ensembles that converge to the Erdős-Rényi and configuration-model ensembles, respectively. We study convergence through von Neumann entropy (VNE)-a network summary statistic measuring information content based on the spectra of a Laplacian matrix-and develop a perturbation analysis for the expected effect of rewiring on VNE. Our analysis yields an estimate for how many rewires are required for a given network to resemble a typical network from an ensemble, offering a computationally efficient quantity for network-ensemble comparison that does not require simulation of the corresponding rewiring process.

Evolutionary prisoner's dilemma games coevolving on adaptive networks

We study a model for switching strategies in the Prisoner's Dilemma game on adaptive networks of player pairings that coevolve as players attempt to maximize their return. We use a node-based strategy model wherein each player follows one strategy at a time (cooperate or defect) across all of its neighbors, changing that strategy and possibly changing partners in response to local changes in the network of player pairing and in the strategies used by connected partners. We compare and contrast numerical simulations with existing pair approximation differential equations for describing this system, as well as more accurate equations developed here using the framework of approximate master equations. We explore the parameter space of the model, demonstrating the relatively high accuracy of the approximate master equations for describing the system observations made from simulations. We study two variations of this partner-switching model to investigate the system evolution, predict stationary states, and compare the total utilities and other qualitative differences between these two model variants.

The scaling structure of the global road network

Because of increasing global urbanization and its immediate consequences, including changes in patterns of food demand, circulation and land use, the next century will witness a major increase in the extent of paved roads built worldwide. To model the effects of this increase, it is crucial to understand whether possible self-organized patterns are inherent in the global road network structure. Here, we use the largest updated database comprising all major roads on the Earth, together with global urban and cropland inventories, to suggest that road length distributions within croplands are indistinguishable from urban ones, once rescaled to account for the difference in mean road length. Such similarity extends to road length distributions within urban or agricultural domains of a given area. We find two distinct regimes for the scaling of the mean road length with the associated area, holding in general at small and at large values of the latter. In suitably large urban and cropland domains, we find that mean and total road lengths increase linearly with their domain area, differently from earlier suggestions. Scaling regimes suggest that simple and universal mechanisms regulate urban and cropland road expansion at the global scale. As such, our findings bear implications for global road infrastructure growth based on land-use change and for planning policies sustaining urban expansions.

Post-Processing Partitions to Identify Domains of Modularity Optimization

We introduce the Convex Hull of Admissible Modularity Partitions (CHAMP) algorithm to prune and prioritize different network community structures identified across multiple runs of possibly various computational heuristics. Given a set of partitions, CHAMP identifies the domain of modularity optimization for each partition—i.e., the parameter-space domain where it has the largest modularity relative to the input set—discarding partitions with empty domains to obtain the subset of partitions that are “admissible” candidate community structures that remain potentially optimal over indicated parameter domains. Importantly, CHAMP can be used for multi-dimensional parameter spaces, such as those for multilayer networks where one includes a resolution parameter and interlayer coupling. Using the results from CHAMP, a user can more appropriately select robust community structures by observing the sizes of domains of optimization and the pairwise comparisons between partitions in the admissible subset. We demonstrate the utility of CHAMP with several example networks. In these examples, CHAMP focuses attention onto pruned subsets of admissible partitions that are 20-to-1785 times smaller than the sets of unique partitions obtained by community detection heuristics that were input into CHAMP.

Super-Resolution Community Detection for Layer-Aggregated Multilayer Networks

Applied network science often involves preprocessing network data before applying a network-analysis method, and there is typically a theoretical disconnect between these steps. For example, it is common to aggregate time-varying network data into windows prior to analysis, and the trade-offs of this preprocessing are not well understood. Focusing on the problem of detecting small communities in multilayer networks, we study the effects of layer aggregation by developing random-matrix theory for modularity matrices associated with layer-aggregated networks with N nodes and L layers, which are drawn from an ensemble of Erdős-Rényi networks with communities planted in subsets of layers. We study phase transitions in which eigenvectors localize onto communities (allowing their detection) and which occur for a given community provided its size surpasses a detectability limit K* . When layers are aggregated via a summation, we obtain [Formula: see text], where T is the number of layers across which the community persists. Interestingly, if T is allowed to vary with L, then summation-based layer aggregation enhances small-community detection even if the community persists across a vanishing fraction of layers, provided that T/L decays more slowly than 𝒪(L-1/2). Moreover, we find that thresholding the summation can, in some cases, cause K* to decay exponentially, decreasing by orders of magnitude in a phenomenon we call super-resolution community detection. In other words, layer aggregation with thresholding is a nonlinear data filter enabling detection of communities that are otherwise too small to detect. Importantly, different thresholds generally enhance the detectability of communities having different properties, illustrating that community detection can be obscured if one analyzes network data using a single threshold.

EIGENVECTOR-BASED CENTRALITY MEASURES FOR TEMPORAL NETWORKS. Multiscale Model Simul 2017;15:537-574. [PMID: 29046619 PMCID: PMC5643020 DOI: 10.1137/16m1066142]

Numerous centrality measures have been developed to quantify the importances of nodes in time-independent networks, and many of them can be expressed as the leading eigenvector of some matrix. With the increasing availability of network data that changes in time, it is important to extend such eigenvector-based centrality measures to time-dependent networks. In this paper, we introduce a principled generalization of network centrality measures that is valid for any eigenvector-based centrality. We consider a temporal network with N nodes as a sequence of T layers that describe the network during different time windows, and we couple centrality matrices for the layers into a supra-centrality matrix of size NT × NT whose dominant eigenvector gives the centrality of each node i at each time t. We refer to this eigenvector and its components as a joint centrality, as it reflects the importances of both the node i and the time layer t. We also introduce the concepts of marginal and conditional centralities, which facilitate the study of centrality trajectories over time. We find that the strength of coupling between layers is important for determining multiscale properties of centrality, such as localization phenomena and the time scale of centrality changes. In the strong-coupling regime, we derive expressions for time-averaged centralities, which are given by the zeroth-order terms of a singular perturbation expansion. We also study first-order terms to obtain first-order-mover scores, which concisely describe the magnitude of nodes' centrality changes over time. As examples, we apply our method to three empirical temporal networks: the United States Ph.D. exchange in mathematics, costarring relationships among top-billed actors during the Golden Age of Hollywood, and citations of decisions from the United States Supreme Court.

Scientific Training in the Era of Big Data: A New Pedagogy for Graduate Education

The era of "big data" has radically altered the way scientific research is conducted and new knowledge is discovered. Indeed, the scientific method is rapidly being complemented and even replaced in some fields by data-driven approaches to knowledge discovery. This paradigm shift is sometimes referred to as the "fourth paradigm" of data-intensive and data-enabled scientific discovery. Interdisciplinary research with a hard emphasis on translational outcomes is becoming the norm in all large-scale scientific endeavors. Yet, graduate education remains largely focused on individual achievement within a single scientific domain, with little training in team-based, interdisciplinary data-oriented approaches designed to translate scientific data into new solutions to today's critical challenges. In this article, we propose a new pedagogy for graduate education: data-centered learning for the domain-data scientist. Our approach is based on four tenets: (1) Graduate training must incorporate interdisciplinary training that couples the domain sciences with data science. (2) Graduate training must prepare students for work in data-enabled research teams. (3) Graduate training must include education in teaming and leadership skills for the data scientist. (4) Graduate training must provide experiential training through academic/industry practicums and internships. We emphasize that this approach is distinct from today's graduate training, which offers training in either data science or a domain science (e.g., biology, sociology, political science, economics, and medicine), but does not integrate the two within a single curriculum designed to prepare the next generation of domain-data scientists. We are in the process of implementing the proposed pedagogy through the development of a new graduate curriculum based on the above four tenets, and we describe herein our strategy, progress, and lessons learned. While our pedagogy was developed in the context of graduate education, the general approach of data-centered learning can and should be applied to students and professionals at any stage of their education, including at the K-12, undergraduate, graduate, and professional levels. We believe that the time is right to embed data-centered learning within our educational system and, thus, generate the talent required to fully harness the potential of big data.

Transitivity reinforcement in the coevolving voter model

One of the fundamental structural properties of many networks is triangle closure. Whereas the influence of this transitivity on a variety of contagion dynamics has been previously explored, existing models of coevolving or adaptive network systems typically use rewiring rules that randomize away this important property, raising questions about their applicability. In contrast, we study here a modified coevolving voter model dynamics that explicitly reinforces and maintains such clustering. Carrying out numerical simulations for a variety of parameter settings, we establish that the transitions and dynamical states observed in coevolving voter model networks without clustering are altered by reinforcing transitivity in the model. We then use a semi-analytical framework in terms of approximate master equations to predict the dynamical behaviors of the model for a variety of parameter settings.

Climate Shocks and Migration: An Agent-Based Modeling Approach

This is a study of migration responses to climate shocks. We construct an agent-based model that incorporates dynamic linkages between demographic behaviors, such as migration, marriage, and births, and agriculture and land use, which depend on rainfall patterns. The rules and parameterization of our model are empirically derived from qualitative and quantitative analyses of a well-studied demographic field site, Nang Rong district, Northeast Thailand. With this model, we simulate patterns of migration under four weather regimes in a rice economy: 1) a reference, 'normal' scenario; 2) seven years of unusually wet weather; 3) seven years of unusually dry weather; and 4) seven years of extremely variable weather. Results show relatively small impacts on migration. Experiments with the model show that existing high migration rates and strong selection factors, which are unaffected by climate change, are likely responsible for the weak migration response.

Enhanced Detectability of Community Structure in Multilayer Networks through Layer Aggregation

Many systems are naturally represented by a multilayer network in which edges exist in multiple layers that encode different, but potentially related, types of interactions, and it is important to understand limitations on the detectability of community structure in these networks. Using random matrix theory, we analyze detectability limitations for multilayer (specifically, multiplex) stochastic block models (SBMs) in which L layers are derived from a common SBM. We study the effect of layer aggregation on detectability for several aggregation methods, including summation of the layers' adjacency matrices for which we show the detectability limit vanishes as O(L^{-1/2}) with increasing number of layers, L. Importantly, we find a similar scaling behavior when the summation is thresholded at an optimal value, providing insight into the common-but not well understood-practice of thresholding pairwise-interaction data to obtain sparse network representations.

Clustering network layers with the strata multilayer stochastic block model

Multilayer networks are a useful data structure for simultaneously capturing multiple types of relationships between a set of nodes. In such networks, each relational definition gives rise to a layer. While each layer provides its own set of information, community structure across layers can be collectively utilized to discover and quantify underlying relational patterns between nodes. To concisely extract information from a multilayer network, we propose to identify and combine sets of layers with meaningful similarities in community structure. In this paper, we describe the "strata multilayer stochastic block model" (sMLSBM), a probabilistic model for multilayer community structure. The central extension of the model is that there exist groups of layers, called "strata", which are defined such that all layers in a given stratum have community structure described by a common stochastic block model (SBM). That is, layers in a stratum exhibit similar node-to-community assignments and SBM probability parameters. Fitting the sMLSBM to a multilayer network provides a joint clustering that yields node-to-community and layer-to-stratum assignments, which cooperatively aid one another during inference. We describe an algorithm for separating layers into their appropriate strata and an inference technique for estimating the SBM parameters for each stratum. We demonstrate our method using synthetic networks and a multilayer network inferred from data collected in the Human Microbiome Project.

Spatiotemporal patterns and trends of Indian monsoonal rainfall extremes

In this study, we provide a comprehensive analysis of trends in the extremes during the Indian summer monsoon (ISM) months (June to September) at different temporal and spatial scales. Our goal is to identify and quantify spatiotemporal patterns and trends that have emerged during the recent decades and may be associated with changing climatic conditions. Our analysis primarily relies on quantile regression that avoids making any subjective choices on spatial, temporal, or intensity pattern of extreme rainfall events. Our analysis divides the Indian monsoon region into climatic compartments that show different and partly opposing trends. These include strong trends towards intensified droughts in Northwest India, parts of Peninsular India, and Myanmar; in contrast, parts of Pakistan, Northwest Himalaya, and Central India show increased extreme daily rain intensity leading to higher flood vulnerability. Our analysis helps explain previously contradicting results of trends in average ISM rainfall.

Network Structure and Biased Variance Estimation in Respondent Driven Sampling

This paper explores bias in the estimation of sampling variance in Respondent Driven Sampling (RDS). Prior methodological work on RDS has focused on its problematic assumptions and the biases and inefficiencies of its estimators of the population mean. Nonetheless, researchers have given only slight attention to the topic of estimating sampling variance in RDS, despite the importance of variance estimation for the construction of confidence intervals and hypothesis tests. In this paper, we show that the estimators of RDS sampling variance rely on a critical assumption that the network is First Order Markov (FOM) with respect to the dependent variable of interest. We demonstrate, through intuitive examples, mathematical generalizations, and computational experiments that current RDS variance estimators will always underestimate the population sampling variance of RDS in empirical networks that do not conform to the FOM assumption. Analysis of 215 observed university and school networks from Facebook and Add Health indicates that the FOM assumption is violated in every empirical network we analyze, and that these violations lead to substantially biased RDS estimators of sampling variance. We propose and test two alternative variance estimators that show some promise for reducing biases, but which also illustrate the limits of estimating sampling variance with only partial information on the underlying population social network.

Network Analysis Reveals Sex- and Antibiotic Resistance-Associated Antivirulence Targets in Clinical Uropathogens

Increasing antibiotic resistance among uropathogenic Escherichia coli (UPEC) is driving interest in therapeutic targeting of nonconserved virulence factor (VF) genes. The ability to formulate efficacious combinations of antivirulence agents requires an improved understanding of how UPEC deploy these genes. To identify clinically relevant VF combinations, we applied contemporary network analysis and biclustering algorithms to VF profiles from a large, previously characterized inpatient clinical cohort. These mathematical approaches identified four stereotypical VF combinations with distinctive relationships to antibiotic resistance and patient sex that are independent of traditional phylogenetic grouping. Targeting resistance- or sex-associated VFs based upon these contemporary mathematical approaches may facilitate individualized anti-infective therapies and identify synergistic VF combinations in bacterial pathogens.

Think locally, act locally: detection of small, medium-sized, and large communities in large networks

It is common in the study of networks to investigate intermediate-sized (or "meso-scale") features to try to gain an understanding of network structure and function. For example, numerous algorithms have been developed to try to identify "communities," which are typically construed as sets of nodes with denser connections internally than with the remainder of a network. In this paper, we adopt a complementary perspective that communities are associated with bottlenecks of locally biased dynamical processes that begin at seed sets of nodes, and we employ several different community-identification procedures (using diffusion-based and geodesic-based dynamics) to investigate community quality as a function of community size. Using several empirical and synthetic networks, we identify several distinct scenarios for "size-resolved community structure" that can arise in real (and realistic) networks: (1) the best small groups of nodes can be better than the best large groups (for a given formulation of the idea of a good community); (2) the best small groups can have a quality that is comparable to the best medium-sized and large groups; and (3) the best small groups of nodes can be worse than the best large groups. As we discuss in detail, which of these three cases holds for a given network can make an enormous difference when investigating and making claims about network community structure, and it is important to take this into account to obtain reliable downstream conclusions. Depending on which scenario holds, one may or may not be able to successfully identify "good" communities in a given network (and good communities might not even exist for a given community quality measure), the manner in which different small communities fit together to form meso-scale network structures can be very different, and processes such as viral propagation and information diffusion can exhibit very different dynamics. In addition, our results suggest that, for many large realistic networks, the output of locally biased methods that focus on communities that are centered around a given seed node (or set of seed nodes) might have better conceptual grounding and greater practical utility than the output of global community-detection methods. They also illustrate structural properties that are important to consider in the development of better benchmark networks to test methods for community detection.