Atypical brainstem modulation of pain might contribute to changes in sensory processing typical of migraine. The study objective was to investigate whether migraine is associated with brainstem structural alterations that correlate with this altered pain processing. MRI T1-weighted images of 55 migraine patients and 58 healthy controls were used to: (1) create deformable mesh models of the brainstem that allow for shape analyses; (2) calculate volumes of the midbrain, pons, medulla and the superior cerebellar peduncles; (3) interrogate correlations between regional brainstem volumes, cutaneous heat pain thresholds, and allodynia symptoms. Migraineurs had smaller midbrain volumes (healthy controls = 61.28 mm³, SD = 5.89; migraineurs = 58.80 mm³, SD = 6.64; p = 0.038), and significant (p < 0.05) inward deformations in the ventral midbrain and pons, and outward deformations in the lateral medulla and dorsolateral pons relative to healthy controls. Migraineurs had a negative correlation between ASC-12 allodynia symptom severity with midbrain volume (r = − 0.32; p = 0.019) and a positive correlation between cutaneous heat pain thresholds with medulla (r = 0.337; p = 0.012) and cerebellar peduncle volumes (r = 0.435; p = 0.001). Migraineurs with greater symptoms of allodynia have smaller midbrain volumes and migraineurs with lower heat pain thresholds have smaller medulla and cerebellar peduncles. The brainstem likely plays a role in altered sensory processing in migraine and brainstem structure might reflect severity of allodynia and hypersensitivity to pain in migraine.

MicroRNAs (miRNAs) regulate gene output by targeting degenerate elements in mRNAs and have undergone drastic expansions in higher metazoan genomes. The evolutionary advantage of maintaining copies of highly similar miRNAs is not well understood, nor is it clear what unique functions, if any, miRNA family members possess. Here, we study evolutionary patterns of metazoan miRNAs, focusing on the targeting preferences of the let-7 and miR-10 families. These studies reveal hotspots for sequence evolution with implications for targeting and secondary structure. High-throughput screening for functional targets reveals that each miRNA represses sites with distinct features and regulates a large number of genes with cooperative function in regulatory networks. Unexpectedly, given the high degree of similarity, single-nucleotide changes grant miRNA family members with distinct targeting preferences. Together, our data suggest complex functional relationships among miRNA duplications, novel expression patterns, sequence change, and the acquisition of new targets.

Network reconstruction is a fundamental problem for understanding many complex systems with unknown interaction structures. In many complex systems, there are indirect interactions between two individuals without immediate connection but with common neighbors. Despite recent advances in network reconstruction, we continue to lack an approach for reconstructing complex networks with indirect interactions. Here we introduce a two-step strategy to resolve the reconstruction problem, where in the first step, we recover both direct and indirect interactions by employing the Lasso to solve a sparse signal reconstruction problem, and in the second step, we use matrix transformation and optimization to distinguish between direct and indirect interactions. The network structure corresponding to direct interactions can be fully uncovered. We exploit the public goods game occurring on complex networks as a paradigm for characterizing indirect interactions and test our reconstruction approach. We find that high reconstruction accuracy can be achieved for both homogeneous and heterogeneous networks, and a number of empirical networks in spite of insufficient data measurement contaminated by noise. Although a general framework for reconstructing complex networks with arbitrary types of indirect interactions is yet lacking, our approach opens new routes to separate direct and indirect interactions in a representative complex system.

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

Improved tools for providing specific intraoperative diagnoses could improve patient care. In neurosurgery, intraoperatively differentiating non-operative lesions such as CNS B-cell lymphoma from operative lesions can be challenging, often necessitating immunohistochemical (IHC) procedures which require up to 24-48 hours. Here, we evaluate the feasibility of generating rapid ex vivo specific labeling using a novel lymphoma-specific fluorescent switchable aptamer. Our B-cell lymphoma-specific switchable aptamer produced only low-level fluorescence in its unbound conformation and generated an 8-fold increase in fluorescence once bound to its target on CD20-positive lymphoma cells. The aptamer demonstrated strong binding to B-cell lymphoma cells within 15 minutes of incubation as observed by flow cytometry. We applied the switchable aptamer to ex vivo xenograft tissue harboring B-cell lymphoma and astrocytoma, and within one hour specific visual identification of lymphoma was routinely possible. In this proof-of-concept study in human cell culture and orthotopic xenografts, we conclude that a fluorescent switchable aptamer can provide rapid and specific labeling of B-cell lymphoma, and that developing aptamer-based labeling approaches could simplify tissue staining and drastically reduce time to histopathological diagnoses compared with IHC-based methods. We propose that switchable aptamers could enhance expeditious, accurate intraoperative decision-making.

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

Background: Genetic profiling represents the future of neuro-oncology but suffers from inadequate biopsies in heterogeneous tumors like Glioblastoma (GBM). Contrast-enhanced MRI (CE-MRI) targets enhancing core (ENH) but yields adequate tumor in only ~60% of cases. Further, CE-MRI poorly localizes infiltrative tumor within surrounding non-enhancing parenchyma, or brain-around-tumor (BAT), despite the importance of characterizing this tumor segment, which universally recurs. In this study, we use multiple texture analysis and machine learning (ML) algorithms to analyze multi-parametric MRI, and produce new images indicating tumor-rich targets in GBM.

Methods: We recruited primary GBM patients undergoing image-guided biopsies and acquired pre-operative MRI: CE-MRI, Dynamic-Susceptibility-weighted-Contrast-enhanced-MRI, and Diffusion Tensor Imaging. Following image coregistration and region of interest placement at biopsy locations, we compared MRI metrics and regional texture with histologic diagnoses of high- vs low-tumor content (≥80% vs <80% tumor nuclei) for corresponding samples. In a training set, we used three texture analysis algorithms and three ML methods to identify MRI-texture features that optimized model accuracy to distinguish tumor content. We confirmed model accuracy in a separate validation set.

Results: We collected 82 biopsies from 18 GBMs throughout ENH and BAT. The MRI-based model achieved 85% cross-validated accuracy to diagnose high- vs low-tumor in the training set (60 biopsies, 11 patients). The model achieved 81.8% accuracy in the validation set (22 biopsies, 7 patients).

Conclusion: Multi-parametric MRI and texture analysis can help characterize and visualize GBM’s spatial histologic heterogeneity to identify regional tumor-rich biopsy targets.

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.