Search

The human immune system depends on a highly diverse collection of antibody-making B cells. B cell receptor sequence diversity is generated by a random recombination process called “rearrangement” forming progenitor B cells, then a Darwinian process of lineage diversification and selection called “affinity maturation.” The resulting receptors can be sequenced in high throughput for research and diagnostics. Such a collection of sequences contains a mixture of various lineages, each of which may be quite numerous, or may consist of only a single member. As a step to understanding the process and result of this diversification, one may wish to reconstruct lineage membership, i.e. to cluster sampled sequences according to which came from the same rearrangement events. We call this clustering problem “clonal family inference.” In this paper we describe and validate a likelihood-based framework for clonal family inference based on a multi-hidden Markov Model (multi-HMM) framework for B cell receptor sequences. We describe an agglomerative algorithm to find a maximum likelihood clustering, two approximate algorithms with various trade-offs of speed versus accuracy, and a third, fast algorithm for finding specific lineages. We show that under simulation these algorithms greatly improve upon existing clonal family inference methods, and that they also give significantly different clusters than previous methods when applied to two real data sets., Author Summary Antibodies must recognize a great diversity of antigens to protect us from infectious disease. The binding properties of antibodies are determined by the DNA sequences of their corresponding B cell receptors (BCRs). These BCR sequences are created in naive form by VDJ recombination, which randomly selects and trims the ends of V, D, and J genes, then joins the resulting segments together with additional random nucleotides. If they pass initial screening and bind an antigen, these sequences then undergo an evolutionary process of reproduction, mutation, and selection, revising the BCR to improve binding to its cognate antigen. It has recently become possible to determine the BCR sequences resulting from this process in high throughput. Although these sequences implicitly contain a wealth of information about both antigen exposure and the process by which we learn to resist pathogens, this information can only be extracted using computer algorithms. In this paper we describe a likelihood-based statistical method to determine, given a collection of BCR sequences, which of them are derived from the same recombination events. It is based on a hidden Markov model (HMM) of VDJ rearrangement which is able to calculate likelihoods for many sequences at once.

Constraining the many biological parameters that govern cortical dynamics is computationally and conceptually difficult because of the curse of dimensionality. This paper addresses these challenges by proposing (1) a novel data-informed mean-field (MF) approach to efficiently map the parameter space of network models; and (2) an organizing principle for studying parameter space that enables the extraction biologically meaningful relations from this high-dimensional data. We illustrate these ideas using a large-scale network model of the Macaque primary visual cortex. Of the 10-20 model parameters, we identify 7 that are especially poorly constrained, and use the MF algorithm in (1) to discover the firing rate contours in this 7D parameter cube. Defining a “biologically plausible” region to consist of parameters that exhibit spontaneous Excitatory and Inhibitory firing rates compatible with experimental values, we find that this region is a slightly thickened codimension-1 submanifold. An implication of this finding is that while plausible regimes depend sensitively on parameters, they are also robust and flexible provided one compensates appropriately when parameters are varied. Our organizing principle for conceptualizing parameter dependence is to focus on certain 2D parameter planes that govern lateral inhibition: Intersecting these planes with the biologically plausible region leads to very simple geometric structures which, when suitably scaled, have a universal character independent of where the intersections are taken. In addition to elucidating the geometry of the plausible region, this invariance suggests useful approximate scaling relations. Our study offers, for the first time, a complete characterization of the set of all biologically plausible parameters for a detailed cortical model, which has been out of reach due to the high dimensionality of parameter space., Author summary Cortical circuits are characterized by a high degree of structural and dynamical complexity, and this biological reality is reflected in the large number of parameters in even semi-realistic cortical models. A fundamental task of computational neuroscience is to understand how these parameters govern network dynamics. While some neuronal parameters can be measured in vivo, many remain poorly constrained due to limitations of available experimental techniques. Computational models can address this problem by relating difficult-to-measure parameters to observable quantities, but to do so one must overcome two challenges: (1) the computational expense of mapping a high dimensional parameter space, and (2) extracting biological insights from such a map. This study aims to address these challenges in the following ways: First, we propose a parsimonious data-informed algorithm that efficiently predicts spontaneous cortical activity, thereby speeding up the mapping of parameter landscapes. Second, we show that lateral inhibition provides a basis for conceptualizing cortical parameter space, enabling us to begin to make sense of its geometric structure and attendant scaling relations. We illustrate our approach on a biologically realistic model of the monkey primary visual cortex.

Since the outbreak of the COVID-19 pandemic, many healthcare facilities have suffered from shortages in medical resources, particularly in Personal Protective Equipment (PPE). In this paper, we propose a game-theoretic approach to schedule PPE orders among healthcare facilities. In this PPE game, each independent healthcare facility optimises its own storage utilisation in order to keep its PPE cost at a minimum. Such a model can reduce peak demand considerably when applied to a variable PPE consumption profile. Experiments conducted for NHS England regions using actual data confirm that the challenge of securing PPE supply during disasters such as COVID-19 can be eased if proper stock management procedures are adopted. These procedures can include early stockpiling, increasing storage capacities and implementing measures that can prolong the time period between successive infection waves, such as social distancing measures. Simulation results suggest that the provision of PPE dedicated storage space can be a viable solution to avoid straining PPE supply chains in case a second wave of COVID-19 infections occurs., Comment: 22 pages, 7 figures, published in PLOS ONE https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0246110

The sex ratio at birth (SRB) in India has been reported imbalanced since the 1970s. Previous studies have shown a great variation in the SRB across geographic locations in India till 2016. As one of the most populous countries and in view of its great regional heterogeneity, it is crucial to produce probabilistic projections for the SRB in India at state level for the purpose of population projection and policy planning. In this paper, we implement a Bayesian hierarchical time series model to project SRB in India by state. We generate SRB probabilistic projections from 2017 to 2030 for 29 States and Union Territories (UTs) in India, and present results in 21 States/UTs with data from the Sample Registration System. Our analysis takes into account two state-specific factors that contribute to sex-selective abortion and resulting sex imbalances at birth: intensity of son preference and fertility squeeze. We project that the largest contribution to female births deficits is in Uttar Pradesh, with cumulative number of missing female births projected to be 2.0 (95% credible interval [1.9; 2.2]) million from 2017 to 2030. The total female birth deficits during 2017-2030 for the whole India is projected to be 6.8 [6.6; 7.0] million.

Occupationally-induced back pain is a leading cause of reduced productivity in industry. Detecting when a worker is lifting incorrectly and at increased risk of back injury presents significant possible benefits. These include increased quality of life for the worker due to lower rates of back injury and fewer workers' compensation claims and missed time for the employer. However, recognizing lifting risk provides a challenge due to typically small datasets and subtle underlying features in accelerometer and gyroscope data. A novel method to classify a lifting dataset using a 2D convolutional neural network (CNN) and no manual feature extraction is proposed in this paper; the dataset consisted of 10 subjects lifting at various relative distances from the body with 720 total trials. The proposed deep CNN displayed greater accuracy (90.6%) compared to an alternative CNN and multilayer perceptron (MLP). A deep CNN could be adapted to classify many other activities that traditionally pose greater challenges in industrial environments due to their size and complexity., Comment: 21 pages, 10 figures

Despite a very effective vaccine, active conflict and community distrust during the ongoing DRC Ebola epidemic are undermining control efforts, including a ring vaccination strategy that requires the prompt immunization of close contacts of infected individuals. However, in April 2019, it was reported 20% or more of close contacts cannot be reached or refuse vaccination, and it is predicted that the ring vaccination strategy would not be effective with such a high level of inaccessibility. The vaccination strategy is now incorporating a “third ring” community-level vaccination that targets members of communities even if they are not known contacts of Ebola cases. To assess the impact of vaccination strategies for controlling Ebola epidemics in the context of variable levels of community accessibility, we employed an individual-level stochastic transmission model that incorporates four sources of heterogeneity: a proportion of the population is inaccessible for contact tracing and vaccination due to lack of confidence in interventions or geographic inaccessibility, two levels of population mixing resembling household and community transmission, two types of vaccine doses with different time periods until immunity, and transmission rates that depend on spatial distance. Our results indicate that a ring vaccination strategy alone would not be effective for containing the epidemic in the context of significant delays to vaccinating contacts even for low levels of household inaccessibility and affirm the positive impact of a supplemental community vaccination strategy. Our key results are that as levels of inaccessibility increase, there is a qualitative change in the effectiveness of the vaccination strategy. For higher levels of vaccine access, the probability that the epidemic will end steadily increases over time, even if probabilities are lower than they would be otherwise with full community participation. For levels of vaccine access that are too low, however, the vaccination strategies are not expected to be successful in ending the epidemic even though they help lower incidence levels, which saves lives, and makes the epidemic easier to contain and reduces spread to other communities. This qualitative change occurs for both types of vaccination strategies: ring vaccination is effective for containing an outbreak until the levels of inaccessibility exceeds approximately 10% in the context of significant delays to vaccinating contacts, a combined ring and community vaccination strategy is effective until the levels of inaccessibility exceeds approximately 50%. More broadly, our results underscore the need to enhance community engagement to public health interventions in order to enhance the effectiveness of control interventions to ensure outbreak containment., Author summary In the context of the ongoing Ebola epidemic in DRC, active conflict and community distrust are undermining control efforts, including vaccination strategies. In this paper, we employed an individual-level stochastic structured transmission model to assess the impact of vaccination strategies on epidemic control in the context of variable levels of household inaccessibility. We found that a ring vaccination strategy of close contacts would not be effective for containing the epidemic in the context of significant delays to vaccinating contacts even for low levels of household inaccessibility and evaluate the impact of a supplemental community vaccination strategy. For lower levels of inaccessibility, the probability of epidemic containment increases over time. For higher levels of inaccessibility, even the combined ring and community vaccination strategies are not expected to contain the epidemic even though they help lower incidence levels, which saves lives, makes the epidemic easier to contain and reduces spread to other communities. We found that ring vaccination is effective for containing an outbreak until the levels of inaccessibility exceeds approximately 10%, a combined ring and community vaccination strategy is effective until the levels of inaccessibility exceeds approximately 50%. Our findings underscore the need to enhance community engagement to public health interventions.

A central problem of neuroscience involves uncovering the principles governing the organization of nervous systems which ensure robustness in brain development. The nematode Caenorhabditis elegans provides us with a model organism for studying this question. In this paper, we focus on the invariant connection structure and spatial arrangement of the neurons comprising the somatic neuronal network of this organism to understand the key developmental constraints underlying its design. We observe that neurons with certain shared characteristics—such as, neural process lengths, birth time cohort, lineage and bilateral symmetry—exhibit a preference for connecting to each other. Recognizing the existence of such homophily and their relative degree of importance in determining connection probability within neurons (for example, in synapses, symmetric pairing is the most dominant factor followed by birth time cohort, process length and lineage) helps in connecting specific neuronal attributes to the topological organization of the network. Further, the functional identities of neurons appear to dictate the temporal hierarchy of their appearance during the course of development. Providing crucial insights into principles that may be common across many organisms, our study shows how the trajectory in the developmental landscape constrains the structural organization of a nervous system., Author summary Nervous systems of even the simplest of organisms involve a large array of neurons communicating with each other through a bewildering pattern of interconnections. What are the forces that govern the formation of this complex network? Using the somatic nervous system of the nematode Caenorhabditis elegans, we show how developmental constraints reflected by the symmetric pairing, birth times, process lengths and cell lineages of individual neurons are related to the observed spatial arrangement of neurons, as well as, the structural and functional organization of the network in the adult nematode. Our work suggests that focusing on the developmental trajectory, which complements approaches involving optimality considerations, may reveal the strategies underlying the development of a nervous system that may be robust across organisms.

Electroencephalogram (EEG) is a common base signal used to monitor brain activity and diagnose sleep disorders. Manual sleep stage scoring is a time-consuming task for sleep experts and is limited by inter-rater reliability. In this paper, we propose an automatic sleep stage annotation method called SleepEEGNet using a single-channel EEG signal. The SleepEEGNet is composed of deep convolutional neural networks (CNNs) to extract time-invariant features, frequency information, and a sequence to sequence model to capture the complex and long short-term context dependencies between sleep epochs and scores. In addition, to reduce the effect of the class imbalance problem presented in the available sleep datasets, we applied novel loss functions to have an equal misclassified error for each sleep stage while training the network. We evaluated the proposed method on different single-EEG channels (i.e., Fpz-Cz and Pz-Oz EEG channels) from the Physionet Sleep-EDF datasets published in 2013 and 2018. The evaluation results demonstrate that the proposed method achieved the best annotation performance compared to current literature, with an overall accuracy of 84.26%, a macro F1-score of 79.66% and Cohen's Kappa coefficient = 0.79. Our developed model is ready to test with more sleep EEG signals and aid the sleep specialists to arrive at an accurate diagnosis. The source code is available at https://github.com/SajadMo/SleepEEGNet.

During the last years, the new science of municipalities has been established as a fertile quantitative approach to systematically understand the urban phenomena. One of its main pillars is the proposition that urban systems display universal scaling behavior regarding socioeconomic, infrastructural and individual basic services variables. This paper discusses the extension of the universality proposition by testing it against a broad range of urban metrics in a developing country urban system. We present an exploration of the scaling exponents for over 6$ variables for the Brazilian urban system. As Brazilian municipalities can deviate significantly from urban settlements, urban-like municipalities were selected based on a systematic density cut-off procedure and the scaling exponents were estimated for this new subset of municipalities. To validate our findings we compared the results for overlaying variables with other studies based on alternative methods. It was found that the analyzed socioeconomic variables follow a superlinear scaling relationship with the population size, and most of the infrastructure and individual basic services variables follow expected sublinear and linear scaling, respectively. However, some infrastructural and individual basic services variables deviated from their expected regimes, challenging the universality hypothesis of urban scaling. We propose that these deviations are a product of top-down decisions/policies. Our analysis spreads over a time-range of 10 years, what is not enough to draw conclusive observations, nevertheless we found hints that the scaling exponent of these variables are evolving towards the expected scaling regime, indicating that the deviations might be temporally constrained and that the urban systems might eventually reach the expected scaling regime., Comment: 14 pages, 6 figures

Cognitive brain imaging is accumulating datasets about the neural substrate of many different mental processes. Yet, most studies are based on few subjects and have low statistical power. Analyzing data across studies could bring more statistical power; yet the current brain-imaging analytic framework cannot be used at scale as it requires casting all cognitive tasks in a unified theoretical framework. We introduce a new methodology to analyze brain responses across tasks without a joint model of the psychological processes. The method boosts statistical power in small studies with specific cognitive focus by analyzing them jointly with large studies that probe less focal mental processes. Our approach improves decoding performance for 80% of 35 widely-different functional-imaging studies. It finds commonalities across tasks in a data-driven way, via common brain representations that predict mental processes. These are brain networks tuned to psychological manipulations. They outline interpretable and plausible brain structures. The extracted networks have been made available; they can be readily reused in new neuro-imaging studies. We provide a multi-study decoding tool to adapt to new data., Author summary Brain-imaging findings in cognitive neuroscience often have low statistical power, despite the availability of functional imaging data across hundreds of studies. Yet, with current analytic frameworks, combining data across studies that map responses to different tasks discards the nuances of the cognitive questions they ask. In this paper, we propose a new approach for fMRI analysis, where a predictive model is used to extract the shared information from many studies together, while respecting their original paradigms. Our method extracts cognitive representations that associate a wide variety of functions to specific brain structures. This provides quantitative improvements and cognitive insights when analyzing together 35 task-fMRI studies; the breadth of the functional data we consider is much higher than in previous work. Reusing the representations learned by our approach also improves statistical power in studies outside the training corpus.

Injuries have a great impact on professional soccer, due to their large influence on team performance and the considerable costs of rehabilitation for players. Existing studies in the literature provide just a preliminary understanding of which factors mostly affect injury risk, while an evaluation of the potential of statistical models in forecasting injuries is still missing. In this paper, we propose a multi-dimensional approach to injury forecasting in professional soccer that is based on GPS measurements and machine learning. By using GPS tracking technology, we collect data describing the training workload of players in a professional soccer club during a season. We then construct an injury forecaster and show that it is both accurate and interpretable by providing a set of case studies of interest to soccer practitioners. Our approach opens a novel perspective on injury prevention, providing a set of simple and practical rules for evaluating and interpreting the complex relations between injury risk and training performance in professional soccer.

In many social complex systems, in which agents are linked by non-linear interactions, the history of events strongly influences the whole network dynamics. However, a class of "commonly accepted beliefs" seems rarely studied. In this paper, we examine how the growth process of a (social) network is influenced by past circumstances. In order to tackle this cause, we simply modify the well known preferential attachment mechanism by imposing a time dependent kernel function in the network evolution equation. This approach leads to a fractional order Barabasi-Albert (BA) differential equation, generalizing the BA model. Our results show that, with passing time, an aging process is observed for the network dynamics. The aging process leads to a decay for the node degree values, thereby creating an opposing process to the preferential attachment mechanism. On one hand, based on the preferential attachment mechanism, nodes with a high degree are more likely to absorb links; but, on the other hand, a node's age has a reduced chance for new connections. This competitive scenario allows an increased chance for younger members to become a hub. Simulations of such a network growth with aging constraint confirm the results found from solving the fractional BA equation. We also report, as an exemplary application, an investigation of the collaboration network between Hollywood movie actors. It is undubiously shown that a decay in the dynamics of their collaboration rate is found, - even including a sex difference. Such findings suggest a widely universal application of the so generalized BA model., Comment: 13 pages; 5 figures; 71 references; as prepared for submission to PLOS ONE

The presence of bacteria or fungi in the bloodstream of patients is abnormal and can lead to life-threatening conditions. A computational model based on a bidirectional long short-term memory artificial neural network, is explored to assist doctors in the intensive care unit to predict whether examination of blood cultures of patients will return positive. As input it uses nine monitored clinical parameters, presented as time series data, collected from 2177 ICU admissions at the Ghent University Hospital. Our main goal is to determine if general machine learning methods and more specific, temporal models, can be used to create an early detection system. This preliminary research obtains an area of 71.95% under the precision recall curve, proving the potential of temporal neural networks in this context.

The friendship paradox is the phenomenon that in social networks, people on average have fewer friends than their friends do. The generalized friendship paradox is an extension to attributes other than the number of friends. The friendship paradox and its generalized version have gathered recent attention due to the information they provide about network structure and local inequalities. In this paper, we propose several measures of nodal qualities which capture different aspects of their activities and influence in online social networks. Using these measures we analyse the prevalence of the generalized friendship paradox over Twitter and we report high levels of prevalence (up to over 90\% of nodes). We contend that this prevalence of the friendship paradox and its generalized version arise because of the hierarchical nature of the connections in the network. This hierarchy is nested as opposed to being star-like. We conclude that these paradoxes are collective phenomena not created merely by a minority of well-connected or high-attribute nodes. Moreover, our results show that a large fraction of individuals can experience the generalized friendship paradox even in the absence of a significant correlation between degrees and attributes.

Fingerprint recognition plays an important role in many commercial applications and is used by millions of people every day, e.g. for unlocking mobile phones. Fingerprint image segmentation is typically the first processing step of most fingerprint algorithms and it divides an image into foreground, the region of interest, and background. Two types of error can occur during this step which both have a negative impact on the recognition performance: 'true' foreground can be labeled as background and features like minutiae can be lost, or conversely 'true' background can be misclassified as foreground and spurious features can be introduced. The contribution of this paper is threefold: firstly, we propose a novel factorized directional bandpass (FDB) segmentation method for texture extraction based on the directional Hilbert transform of a Butterworth bandpass (DHBB) filter interwoven with soft-thresholding. Secondly, we provide a manually marked ground truth segmentation for 10560 images as an evaluation benchmark. Thirdly, we conduct a systematic performance comparison between the FDB method and four of the most often cited fingerprint segmentation algorithms showing that the FDB segmentation method clearly outperforms these four widely used methods. The benchmark and the implementation of the FDB method are made publicly available. Open-Access Publikationsfonds 2016 peerReviewed

Purpose: Respiratory gated radiation therapy (RGRT) is used to minimize the radiation dose to normal tissue in lung cancer patients. Determinations of the gating window in the respiratory phase of patients are important in RGRT but it is not easy. The objective of this study was to determine the optimal gating window with a visible guiding system in RGRT. Materials and Methods: Between April and October in 2014 the breathing signals of 23 lung cancer patients were recorded with a Real-time Position Management (RPM) respiratory gating system (Varian, USA). We performed statistical analysis with breathing signals to find the optimal gating window for the guided breathing for RGRT. Results: 19 of the 23 patients showed statistically significant differences (p < 0.05) when the breathing signals obtained before and after breathing training were compared, The standard deviation of the respiration signals after breathing training was the lowest in the phase of 30 % - 70 % (p < 0.05). Conclusions: RGRT with RPM system was confirmed the usefulness of the visible guiding system. Respiratory regularity was significantly improved with the use of the RPM with our visible guidance system, therefore, it would help improve the accuracy and efficiency of RGRT., This paper has been withdrawn by the author due to a fatal errors in tables and figures