ViLoN-a multi-layer network approach to data integration demonstrated for patient stratification.

With more and more data being collected, modern network representations exploit the complementary nature of different data sources as well as similarities across patients. We here introduce the Variation of information fused Layers of Networks algorithm (ViLoN), a novel network-based approach for the integration of multiple molecular profiles. As a key innovation, it directly incorporates prior functional knowledge (KEGG, GO). In the constructed network of patients, patients are represented by networks of pathways, comprising genes that are linked by common functions and joint regulation in the disease. Patient stratification remains a key challenge both in the clinic and for research on disease mechanisms and treatments. We thus validated ViLoN for patient stratification on multiple data type combinations (gene expression, methylation, copy number), showing substantial improvements and consistently competitive performance for all. Notably, the incorporation of prior functional knowledge was critical for good results in the smaller cohorts (rectum adenocarcinoma: 90, esophageal carcinoma: 180), where alternative methods failed.

Distinguishing between different percolation regimes in noisy dynamic networks with an application to epileptic seizures.

In clinical neuroscience, epileptic seizures have been associated with the sudden emergence of coupled activity across the brain. The resulting functional networks-in which edges indicate strong enough coupling between brain regions-are consistent with the notion of percolation, which is a phenomenon in complex networks corresponding to the sudden emergence of a giant connected component. Traditionally, work has concentrated on noise-free percolation with a monotonic process of network growth, but real-world networks are more complex. We develop a class of random graph hidden Markov models (RG-HMMs) for characterizing percolation regimes in noisy, dynamically evolving networks in the presence of edge birth and edge death. This class is used to understand the type of phase transitions undergone in a seizure, and in particular, distinguishing between different percolation regimes in epileptic seizures. We develop a hypothesis testing framework for inferring putative percolation mechanisms. As a necessary precursor, we present an EM algorithm for estimating parameters from a sequence of noisy networks only observed at a longitudinal subsampling of time points. Our results suggest that different types of percolation can occur in human seizures. The type inferred may suggest tailored treatment strategies and provide new insights into the fundamental science of epilepsy.

Predicting Emission Wavelengths in Benzobisoxazole-Based OLEDs with Gradient Boosted Ensemble Models.

We demonstrate the use of gradient-boosted ensemble models that accurately predict emission wavelengths in benzobis[1,2-d:4,5-d']oxazole (BBO) based fluorescent emitters. We have curated a database of 50 molecules from previously published data by the Jeffries-EL group using density functional theory (DFT) computed ground and excited state features. We consider two machine learning (ML) models based on (i) whole cruciform molecules and (ii) their constituent fragment molecules. Both ML models provide accurate predictions with root-mean-square errors between 30 and 36 nm, competitive with state-of-the-art deep learning models trained on orders of magnitude more molecules, and this accuracy holds even when tested on four new BBO emitters unseen by the models. We also provide an interpretable feature importance analysis and discuss the relevant relationships between DFT and changes in predicted emission wavelength.

Early Introduction and Rise of the Omicron Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variant in Highly Vaccinated University Populations.

BACKGROUND: The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly transmissible in vaccinated and unvaccinated populations. The dynamics that govern its establishment and propensity toward fixation (reaching 100% frequency in the SARS-CoV-2 population) in communities remain unknown. Here, we describe the dynamics of Omicron at 3 institutions of higher education (IHEs) in the greater Boston area. METHODS: We use diagnostic and variant-specifying molecular assays and epidemiological analytical approaches to describe the rapid dominance of Omicron following its introduction into 3 IHEs with asymptomatic surveillance programs. RESULTS: We show that the establishment of Omicron at IHEs precedes that of the state and region and that the time to fixation is shorter at IHEs (9.5-12.5 days) than in the state (14.8 days) or region. We show that the trajectory of Omicron fixation among university employees resembles that of students, with a 2- to 3-day delay. Finally, we compare cycle threshold values in Omicron vs Delta variant cases on college campuses and identify lower viral loads among college affiliates who harbor Omicron infections. CONCLUSIONS: We document the rapid takeover of the Omicron variant at IHEs, reaching near-fixation within the span of 9.5-12.5 days despite lower viral loads, on average, than the previously dominant Delta variant. These findings highlight the transmissibility of Omicron, its propensity to rapidly dominate small populations, and the ability of robust asymptomatic surveillance programs to offer early insights into the dynamics of pathogen arrival and spread.

Estimation of heterogeneous instantaneous reproduction numbers with application to characterize SARS-CoV-2 transmission in Massachusetts counties.

The reproductive number is an important metric that has been widely used to quantify the infectiousness of communicable diseases. The time-varying instantaneous reproductive number is useful for monitoring the real-time dynamics of a disease to inform policy making for disease control. Local estimation of this metric, for instance at a county or city level, allows for more targeted interventions to curb transmission. However, simultaneous estimation of local reproductive numbers must account for potential sources of heterogeneity in these time-varying quantities-a key element of which is human mobility. We develop a statistical method that incorporates human mobility between multiple regions for estimating region-specific instantaneous reproductive numbers. The model also can account for exogenous cases imported from outside of the regions of interest. We propose two approaches to estimate the reproductive numbers, with mobility data used to adjust incidence in the first approach and to inform a formal priori distribution in the second (Bayesian) approach. Through a simulation study, we show that region-specific reproductive numbers can be well estimated if human mobility is reasonably well approximated by available data. We use this approach to estimate the instantaneous reproductive numbers of COVID-19 for 14 counties in Massachusetts using CDC case report data and the human mobility data collected by SafeGraph. We found that, accounting for mobility, our method produces estimates of reproductive numbers that are distinct across counties. In contrast, independent estimation of county-level reproductive numbers tends to produce similar values, as trends in county case-counts for the state are fairly concordant. These approaches can also be used to estimate any heterogeneity in transmission, for instance, age-dependent instantaneous reproductive number estimates. As people are more mobile and interact frequently in ways that permit transmission, it is important to account for this in the estimation of the reproductive number.

Assessing the impact of aircraft arrival on ambient ultrafine particle number concentrations in near-airport communities in Boston, Massachusetts.

Aircraft emissions contribute to overall ambient air pollution, including ultrafine particle (UFP) concentrations. However, accurately ascertaining aviation contributions to UFP is challenging due to high spatiotemporal variability along with intermittent aviation emissions. The objective of this study was to evaluate the impact of arrival aircraft on particle number concentration (PNC), a proxy for UFP, across six study sites 3-17 km from a major arrival aircraft flight path into Boston Logan International Airport by utilizing real-time aircraft activity and meteorological data. Ambient PNC at all monitoring sites was similar at the median but had greater variation at the 95th and 99th percentiles with more than two-fold increases in PNC observed at sites closer to the airport. PNC was elevated during the hours with high aircraft activity with sites closest to the airport exhibiting stronger signals when downwind from the airport. Regression models indicated that the number of arrival aircraft per hour was associated with measured PNC at all six sites, with a maximum contribution of 50% of total PNC at a monitor 3 km from the airport during hours with arrival activity on the flight path of interest (26% across all hours). Our findings suggest strong but intermittent contributions from arrival aircraft to ambient PNC in communities near airports.

Projecting Quarantine Utilization During a Pandemic.

Objectives. To develop an approach to project quarantine needs during an outbreak, particularly for communally housed individuals who interact with outside individuals. Methods. We developed a method that uses basic surveillance data to do short-term projections of future quarantine needs. The development of this method was rigorous, but it is conceptually simple and easy to implement and allows one to anticipate potential superspreading events. We demonstrate how this method can be used with data from the fall 2020 semester of a large urban university in Boston, Massachusetts, that provided quarantine housing for students living on campus in response to the COVID-19 pandemic. Our approach accounted for potentially infectious interactions between individuals living in university housing and those who did not. Results. Our approach was able to accurately project 10-day-ahead quarantine utilization for on-campus students in a large urban university. Our projections were most accurate when we anticipated weekend superspreading events around holidays. Conclusions. We provide an easy-to-use software tool to project quarantine utilization for institutions that can account for mixing with outside populations. This software tool has potential application for universities, corrections facilities, and the military. (Am J Public Health. 2022;112(2):277-283. https://doi.org/10.2105/AJPH.2021.306573).

Sensor-based localization of epidemic sources on human mobility networks.

We investigate the source detection problem in epidemiology, which is one of the most important issues for control of epidemics. Mathematically, we reformulate the problem as one of identifying the relevant component in a multivariate Gaussian mixture model. Focusing on the study of cholera and diseases with similar modes of transmission, we calibrate the parameters of our mixture model using human mobility networks within a stochastic, spatially explicit epidemiological model for waterborne disease. Furthermore, we adopt a Bayesian perspective, so that prior information on source location can be incorporated (e.g., reflecting the impact of local conditions). Posterior-based inference is performed, which permits estimates in the form of either individual locations or regions. Importantly, our estimator only requires first-arrival times of the epidemic by putative observers, typically located only at a small proportion of nodes. The proposed method is demonstrated within the context of the 2000-2002 cholera outbreak in the KwaZulu-Natal province of South Africa.

Estimation of local time-varying reproduction numbers in noisy surveillance data.

A valuable metric in understanding local infectious disease dynamics is the local time-varying reproduction number, i.e. the expected number of secondary local cases caused by each infected individual. Accurate estimation of this quantity requires distinguishing cases arising from local transmission from those imported from elsewhere. Realistically, we can expect identification of cases as local or imported to be imperfect. We study the propagation of such errors in estimation of the local time-varying reproduction number. In addition, we propose a Bayesian framework for estimation of the true local time-varying reproduction number when identification errors exist. And we illustrate the practical performance of our estimator through simulation studies and with outbreaks of COVID-19 in Hong Kong and Victoria, Australia. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.

Convex combination sequence kernel association test for rare-variant studies.

We propose a novel variant set test for rare-variant association studies, which leverages multiple single-nucleotide variant (SNV) annotations. Our approach optimizes a convex combination of different sequence kernel association test (SKAT) statistics, where each statistic is constructed from a different annotation and combination weights are optimized through a multiple kernel learning algorithm. The combination test statistic is evaluated empirically through data splitting. In simulations, we find our method preserves type I error at α=2.5×10-6 and has greater power than SKAT(-O) when SNV weights are not misspecified and sample sizes are large ( N≥5,000 ). We utilize our method in the Framingham Heart Study (FHS) to identify SNV sets associated with fasting glucose. While we are unable to detect any genome-wide significant associations between fasting glucose and 4-kb windows of rare variants ( p<10-7 ) in 6,419 FHS participants, our method identifies suggestive associations between fasting glucose and rare variants near ROCK2 ( p=2.1×10-5 ) and within CPLX1 ( p=5.3×10-5 ). These two genes were previously reported to be involved in obesity-mediated insulin resistance and glucose-induced insulin secretion by pancreatic beta-cells, respectively. These findings will need to be replicated in other cohorts and validated by functional genomic studies.

Detection of multiple perturbations in multi-omics biological networks.

Cellular mechanism-of-action is of fundamental concern in many biological studies. It is of particular interest for identifying the cause of disease and learning the way in which treatments act against disease. However, pinpointing such mechanisms is difficult, due to the fact that small perturbations to the cell can have wide-ranging downstream effects. Given a snapshot of cellular activity, it can be challenging to tell where a disturbance originated. The presence of an ever-greater variety of high-throughput biological data offers an opportunity to examine cellular behavior from multiple angles, but also presents the statistical challenge of how to effectively analyze data from multiple sources. In this setting, we propose a method for mechanism-of-action inference by extending network filtering to multi-attribute data. We first estimate a joint Gaussian graphical model across multiple data types using penalized regression and filter for network effects. We then apply a set of likelihood ratio tests to identify the most likely site of the original perturbation. In addition, we propose a conditional testing procedure to allow for detection of multiple perturbations. We demonstrate this methodology on paired gene expression and methylation data from The Cancer Genome Atlas (TCGA).

Initial evaluation of the effects of an environmental-focused problem-solving intervention for transition-age young people with developmental disabilities: Project TEAM.

AIM: Project TEAM (Teens making Environment and Activity Modifications) teaches transition-age young people with developmental disabilities, including those with co-occurring intellectual or cognitive disabilities, to identify and resolve environmental barriers to participation. We examined its effects on young people's attainment of participation goals, knowledge, problem-solving, self-determination, and self-efficacy. METHOD: We used a quasi-experimental, repeated measures design (initial, outcome, 6-week follow-up) with two groups: (1) Project TEAM (28 males, 19 females; mean age 17y 6mo); and (2) goal-setting comparison (21 males, 14 females; mean age 17y 6mo). A matched convenience sample was recruited in two US states. Attainment of participation goals and goal attainment scaling (GAS) T scores were compared at outcome. Differences between groups for all other outcomes were analyzed using linear mixed effects models. RESULTS: At outcome, Project TEAM participants demonstrated greater knowledge (estimated mean difference: 1.82; confidence interval [CI]: 0.90, 2.74) and ability to apply knowledge during participation (GAS: t[75]=4.21; CI: 5.21, 14.57) compared to goal-setting. While both groups achieved significant improvements in knowledge, problem-solving, and self-determination, increases in parent reported self-determination remained at 6-week follow-up only for Project TEAM (estimated mean difference: 4.65; CI: 1.32, 7.98). Significantly more Project TEAM participants attained their participation goals by follow-up (Project TEAM=97.6%, goal-setting=77.1%, p=0.009). INTERPRETATION: Both approaches support attainment of participation goals. Although inconclusive, Project TEAM may uniquely support young people with developmental disabilities to act in a self-determined manner and apply an environmental problem-solving approach over time. WHAT THIS PAPER ADDS: Individualized goal-setting, alone or during Project TEAM (Teens making Environment and Activity Modifications) appears to support attainment of participation goals. Project TEAM appears to support young people with developmental disabilities to apply an environmental problem-solving approach to participation barriers. Parents of young people with developmental disabilities report sustained changes in self-determination 6 weeks after Project TEAM.

Molecular adaptations of striatal spiny projection neurons during levodopa-induced dyskinesia.

Levodopa treatment is the major pharmacotherapy for Parkinson's disease. However, almost all patients receiving levodopa eventually develop debilitating involuntary movements (dyskinesia). Although it is known that striatal spiny projection neurons (SPNs) are involved in the genesis of this movement disorder, the molecular basis of dyskinesia is not understood. In this study, we identify distinct cell-type-specific gene-expression changes that occur in subclasses of SPNs upon induction of a parkinsonian lesion followed by chronic levodopa treatment. We identify several hundred genes, the expression of which is correlated with levodopa dose, many of which are under the control of activator protein-1 and ERK signaling. Despite homeostatic adaptations involving several signaling modulators, activator protein-1-dependent gene expression remains highly dysregulated in direct pathway SPNs upon chronic levodopa treatment. We also discuss which molecular pathways are most likely to dampen abnormal dopaminoceptive signaling in spiny projection neurons, hence providing potential targets for antidyskinetic treatments in Parkinson's disease.

Meta-Analysis for Penalized Regression Methods with Multi-Cohort Genome-Wide Association Studies.

OBJECTIVE: Penalized regression has been successfully applied in genome-wide association studies. While meta-analysis is often conducted to increase power and protect patients' confidentiality, methods for meta-analyzing results of penalized regression in multi-cohort setting are still under development. METHODS: We propose to use a data-splitting method to obtain valid p values (or equivalently, coefficient estimates and standard errors) for meta-analysis across multiple cohorts. We examine two ways of splitting data in multi-cohort setting and propose three methods to conduct meta-analysis based on p values. We compare the three meta-analysis methods to mega-analysis, which consists of pooling individual level data. We also apply our proposed meta-analysis approaches to the Framingham Heart Study data, where we divide the original dataset into four parts to create a multi-cohort scenario. RESULTS: The simulations suggest that splitting cohorts has better performance than splitting data within each cohort. The real data application also shows that this method provides results that are similar to the mega-analysis. CONCLUSION: After comparing the three methods that we proposed to conduct meta-analysis, we recommend splitting cohorts rather than datasets to obtain valid p values for meta-analysis of results from penalized regression in multi-cohort setting.

Modeling gene-covariate interactions in sparse regression with group structure for genome-wide association studies.

In genome-wide association studies (GWAS), it is of interest to identify genetic variants associated with phenotypes. For a given phenotype, the associated genetic variants are usually a sparse subset of all possible variants. Traditional Lasso-type estimation methods can therefore be used to detect important genes. But the relationship between genotypes at one variant and a phenotype may be influenced by other variables, such as sex and life style. Hence it is important to be able to incorporate gene-covariate interactions into the sparse regression model. In addition, because there is biological knowledge on the manner in which genes work together in structured groups, it is desirable to incorporate this information as well. In this paper, we present a novel sparse regression methodology for gene-covariate models in association studies that not only allows such interactions but also considers biological group structure. Simulation results show that our method substantially outperforms another method, in which interaction is considered, but group structure is ignored. Application to data on total plasma immunoglobulin E (IgE) concentrations in the Framingham Heart Study (FHS), using sex and smoking status as covariates, yields several potentially interesting gene-covariate interactions.

On the Question of Effective Sample Size in Network Modeling: An Asymptotic Inquiry.

The modeling and analysis of networks and network data has seen an explosion of interest in recent years and represents an exciting direction for potential growth in statistics. Despite the already substantial amount of work done in this area to date by researchers from various disciplines, however, there remain many questions of a decidedly foundational nature - natural analogues of standard questions already posed and addressed in more classical areas of statistics - that have yet to even be posed, much less addressed. Here we raise and consider one such question in connection with network modeling. Specifically, we ask, "Given an observed network, what is the sample size?" Using simple, illustrative examples from the class of exponential random graph models, we show that the answer to this question can very much depend on basic properties of the networks expected under the model, as the number of vertices nV in the network grows. In particular, adopting the (asymptotic) scaling of the variance of the maximum likelihood parameter estimates as a notion of effective sample size, say neff, we show that whether the networks are sparse or not under our model (i.e., having relatively few or many edges between vertices, respectively) is sufficient to yield an order of magnitude difference in neff, from O(nV ) to [Formula: see text]. We then explore some practical implications of this result, using both simulation and data on food-sharing from Lamalera, Indonesia.

Network-guided sparse regression modeling for detection of gene-by-gene interactions.

MOTIVATION: Genetic variants identified by genome-wide association studies to date explain only a small fraction of total heritability. Gene-by-gene interaction is one important potential source of unexplained total heritability. We propose a novel approach to detect such interactions that uses penalized regression and sparse estimation principles, and incorporates outside biological knowledge through a network-based penalty. RESULTS: We tested our new method on simulated and real data. Simulation showed that with reasonable outside biological knowledge, our method performs noticeably better than stage-wise strategies (i.e. selecting main effects first, and interactions second, from those main effects selected) in finding true interactions, especially when the marginal strength of main effects is weak. We applied our method to Framingham Heart Study data on total plasma immunoglobulin E (IgE) concentrations and found a number of interactions among different classes of human leukocyte antigen genes that may interact to influence the risk of developing IgE dysregulation and allergy. AVAILABILITY: The proposed method is implemented in R and available at http://math.bu.edu/people/kolaczyk/software.html. CONTACT: chenlu@bu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Network-based prediction for sources of transcriptional dysregulation using latent pathway identification analysis.

Understanding the systemic biological pathways and the key cellular mechanisms that dictate disease states, drug response, and altered cellular function poses a significant challenge. Although high-throughput measurement techniques, such as transcriptional profiling, give some insight into the altered state of a cell, they fall far short of providing by themselves a complete picture. Some improvement can be made by using enrichment-based methods to, for example, organize biological data of this sort into collections of dysregulated pathways. However, such methods arguably are still limited to primarily a transcriptional view of the cell. Augmenting these methods still further with networks and additional -omics data has been found to yield pathways that play more fundamental roles. We propose a previously undescribed method for identification of such pathways that takes a more direct approach to the problem than any published to date. Our method, called latent pathway identification analysis (LPIA), looks for statistically significant evidence of dysregulation in a network of pathways constructed in a manner that implicitly links pathways through their common function in the cell. We describe the LPIA methodology and illustrate its effectiveness through analysis of data on (i) metastatic cancer progression, (ii) drug treatment in human lung carcinoma cells, and (iii) diagnosis of type 2 diabetes. With these analyses, we show that LPIA can successfully identify pathways whose perturbations have latent influences on the transcriptionally altered genes.

Emergence of persistent networks in long-term intracranial EEG recordings.

Over the past two decades, the increased ability to analyze network relationships among neural structures has provided novel insights into brain function. Most network approaches, however, focus on static representations of the brain's physical or statistical connectivity. Few studies have examined how brain functional networks evolve spontaneously over long epochs of continuous time. To address this, we examine functional connectivity networks deduced from continuous long-term electrocorticogram recordings. For a population of six human patients, we identify a persistent pattern of connections that form a frequency-band-dependent network template, and a set of core connections that appear frequently and together. These structures are robust, emerging from brief time intervals (~100 s) regardless of cognitive state. These results suggest that a metastable, frequency-band-dependent scaffold of brain connectivity exists from which transient activity emerges and recedes.

Binge drinking trajectories from adolescence to young adulthood: the effects of peer social network.

This study investigates an association between social network characteristics and binge drinking from adolescence to young adulthood, utilizing National Longitudinal Study of Adolescent Health (n = 7,966) and employing social network and longitudinal analysis. Lower integration and socialization with alcohol-using peers had immediate risks of binge drinking during adolescence; however, over time, the effects of socialization with alcohol-using peers had the most dramatic reduction. The most prestigious adolescents had the highest longitudinal risks of binge drinking, although they had no immediate risk. Alcohol consumption-related interventions overlooking longitudinal dynamics of social networks may not effectively prevent adolescents from binge drinking in young adulthood.