In this paper, a robust sequential dictionary learning (DL) algorithm is presented. The proposed algorithm is motivated from the maximum likelihood perspective on dictionary learning and its link to the minimization of the Kullback-Leibler divergence. It is obtained by using a robust loss function in the data fidelity term of the DL objective instead of the usual quadratic loss. The proposed robust loss function is derived from the α-divergence as an alternative to the Kullback-Leibler divergence, which leads to a quadratic loss. Compared to other robust approaches, the proposed loss has the advantage of belonging to class of redescending M-estimators, guaranteeing inference stability from large deviations from the Gaussian nominal noise model. The algorithm is obtained by solving a sequence of penalized rank-1 matrix approximation problems, where the ℓ 1 -norm is introduced as a penalty promoting sparsity and then using a block coordinate descent approach to estimate the unknowns. Performance comparison with similar robust DL algorithms on digit recognition, background removal, and gray-scale image denoising is performed highlighting efficacy of the proposed algorithm.

Oligodendrocytes are the myelin-producing cells of the central nervous system that are responsible for electrically insulating axons to speed the propagation of electrical impulses. A striking feature of oligodendrocyte development within white matter is that the cell bodies of many oligodendrocyte progenitor cells become organised into discrete linear arrays of three or more cells before they differentiate into myelin-producing oligodendrocytes. These linear arrays align parallel to the direction of the axons within white matter tracts and are believed to play an important role in the co-ordination of myelination. Guided by experimental data on the abundance and composition of linear arrays in the corpus callosum of the postnatal mouse brain, we construct discrete and continuous models of linear array generation to specifically investigate the relative influence of cell migration, proliferation, differentiation and death of oligodendroglia upon the genesis of linear arrays during early postnatal development. We demonstrate that only models that incorporate significant cell migration can replicate all of the experimental observations on number of arrays, number of cells in arrays and total cell count of oligodendroglia within a given area of the corpus callosum. These models are also necessary to accurately reflect experimental data on the abundance of linear arrays composed of oligodendrocytes that derive from progenitors of different clonal origins.

The cellular mechanisms that regulate the topographic arrangement of myelin internodes along axons remain largely uncharacterized. Recent clonal analysis of oligodendrocyte morphologies in the mouse optic nerve revealed that adjacent oligodendrocytes frequently formed adjacent internodes on one or more axons in common, whereas oligodendrocytes in the optic nerve were never observed to myelinate the same axon more than once. By modelling the process of axonal selection at the single cell level, we demonstrate that internode length and primary process length constrain the capacity of oligodendrocytes to myelinate the same axon more than once. On the other hand, probabilistic analysis reveals that the observed juxtaposition of myelin internodes among common sets of axons by adjacent oligodendrocytes is highly unlikely to occur by chance. Our analysis may reveal a hitherto unknown level of communication between adjacent oligodendrocytes in the selection of axons for myelination. Together, our analyses provide novel insights into the mechanisms that define the spatial organization of myelin internodes within white matter at the single cell level.

We develop an off-lattice, agent-based model to describe vasculogenesis, the de novo formation of blood vessels from endothelial progenitor cells during development. The endothelial cells that comprise our vessel network are viewed as linearly elastic spheres that move in response to the forces they experience. We distinguish two types of endothelial cells: vessel elements are contained within the network and tip cells are located at the ends of vessels. Tip cells move in response to mechanical forces caused by interactions with neighbouring vessel elements and the local tissue environment, chemotactic forces and a persistence force which accounts for their tendency to continue moving in the same direction. Vessel elements are subject to similar mechanical forces but are insensitive to chemotaxis. An angular persistence force representing interactions with the local tissue is introduced to stabilise buckling instabilities caused by cell proliferation. Only vessel elements proliferate, at rates which depend on their degree of stretch: elongated elements have increased rates of proliferation, and compressed elements have reduced rates. Following division, the fate of the new cell depends on the local mechanical environment: the probability of forming a new sprout is increased if the parent vessel is highly compressed and the probability of being incorporated into the parent vessel increased if the parent is stretched. Simulation results reveal that our hybrid model can reproduce the key qualitative features of vasculogenesis. Extensive parameter sensitivity analyses show that significant changes in network size and morphology are induced by varying the chemotactic sensitivity of tip cells, and the sensitivities of the proliferation rate and the sprouting probability to mechanical stretch. Varying the chemotactic sensitivity directly influences the directionality of the networks. The degree of branching, and thereby the density of the networks, is influenced by the sprouting probability. Glyphs that simultaneously depict several network properties are introduced to show how these and other network quantities change over time and also as model parameters vary. We also show how equivalent glyphs constructed from in vivo data could be used to discriminate between normal and tumour vasculature and, in the longer term, for model validation. We conclude that our biomechanical hybrid model can generate vascular networks that are qualitatively similar to those generated from in vitro and in vivo experiments.

The myelin sheath that insulates some axons in the central nervous system allows for faster signal conduction. Previously, axons were thought to be either unmyelinated or fully myelinated. Recent experimental work has discovered a new pattern of myelination (intermittent myelination) along axons in the mouse brain, in which long unmyelinated axon segments are followed by myelinated segments of comparable length. We use a computational model to explore how myelin distribution (in particular intermittent myelination) affects conduction velocity. We find that although fully myelinated axons minimize conduction velocity, varying the spatial distribution of a fixed amount of myelin along a partially myelinated axon leads to considerable variation in the conduction velocity for action potentials. Whether sodium ion channel number or sodium ion channel density is held constant as the area of the unmyelinated segments increases has a strong influence on the optimal pattern of myelin and the conduction velocity.

Motivated by in vitro time-lapse images of ovarian cancer spheroids inducing mesothelial cell clearance, the traditional agent-based model of cell migration, based on simple volume exclusion, was extended to include the possibility that a cell seeking to move into an occupied location may push the resident cell, and any cells neighbouring it, out of the way to occupy that location. In traditional discrete models of motile cells with volume exclusion such a move would be aborted. We introduce a new shoving mechanism which allows cells to choose the direction to shove cells that expends the least amount of shoving effort (to account for the likely resistance of cells to being pushed). We call this motility rule 'smart shoving'. We examine whether agent-based simulations of different shoving mechanisms can be distinguished on the basis of single realisations and averages over many realisations. We emphasise the difficulty in distinguishing cell mechanisms from cellular automata simulations based on snap-shots of cell distributions, site-occupancy averages and the evolution of the number of cells of each species averaged over many realisations. This difficulty suggests the need for higher resolution cell tracking.

We consider a phenotype-structured model of evolutionary dynamics in a population of cancer cells exposed to the action of a cytotoxic drug. The model consists of a nonlocal parabolic equation governing the evolution of the cell population density function. We develop a novel method for constructing exact solutions to the model equation, which allows for a systematic investigation of the way in which the size and the phenotypic composition of the cell population change in response to variations of the drug dose and other evolutionary parameters. Moreover, we address numerical optimal control for a calibrated version of the model based on biological data from the existing literature, in order to identify the drug delivery schedule that makes it possible to minimise either the population size at the end of the treatment or the average population size during the course of treatment. The results obtained challenge the notion that traditional high-dose therapy represents a “one-fits-all solution” in anticancer therapy by showing that the continuous administration of a relatively low dose of the cytotoxic drug performs more closely to i.e. the optimal dosing regimen to minimise the average size of the cancer cell population during the course of treatment.

In many biological systems, motile agents exhibit random motion with short-term directional persistence, together with crowding effects arising from spatial exclusion. We formulate and study a class of lattice-based models for multiple walkers with motion persistence and spatial exclusion in one and two dimensions, and use a mean-field approximation to investigate relevant population-level partial differential equations in the continuum limit. We show that this model of a persistent exclusion process is in general well described by a nonlinear diffusion equation. With reference to results presented in the current literature, our results reveal that the nonlinearity arises from the combination of motion persistence and volume exclusion, with linearity in terms of the canonical diffusion or heat equation being recovered in either the case of persistence without spatial exclusion, or spatial exclusion without persistence. We generalize our results to include systems of multiple species of interacting, motion-persistent walkers, as well as to incorporate a global drift in addition to persistence. These models are shown to be governed approximately by systems of nonlinear advection-diffusion equations. By comparing the prediction of the mean-field approximation to stochastic simulation results, we assess the performance of our results. Finally, we also address the problem of inferring the presence of persistence from simulation results, with a view to application to experimental cell-imaging data.

The epilepsies affect around 65 million people worldwide and have a substantial missing heritability component. We report a genome-wide mega-analysis involving 15,212 individuals with epilepsy and 29,677 controls, which reveals 16 genome-wide significant loci, of which 11 are novel. Using various prioritization criteria, we pinpoint the 21 most likely epilepsy genes at these loci, with the majority in genetic generalized epilepsies. These genes have diverse biological functions, including coding for ion-channel subunits, transcription factors and a vitamin-B6 metabolism enzyme. Converging evidence shows that the common variants associated with epilepsy play a role in epigenetic regulation of gene expression in the brain. The results show an enrichment for monogenic epilepsy genes as well as known targets of antiepileptic drugs. Using SNP-based heritability analyses we disentangle both the unique and overlapping genetic basis to seven different epilepsy subtypes. Together, these findings provide leads for epilepsy therapies based on underlying pathophysiology.

We report a genome-wide association scan in >6,000 Latin Americans for pigmentation of skin and eyes. We found eighteen signals of association at twelve genomic regions. These include one novel locus for skin pigmentation (in 10q26) and three novel loci for eye pigmentation (in 1q32, 20q13 and 22q12). We demonstrate the presence of multiple independent signals of association in the 11q14 and 15q13 regions (comprising the GRM5/TYR and HERC2/OCA2 genes, respectively) and several epistatic interactions among independently associated alleles. Strongest association with skin pigmentation at 19p13 was observed for an Y182H missense variant (common only in East Asians and Native Americans) in MFSD12, a gene recently associated with skin pigmentation in Africans. We show that the frequency of the derived allele at Y182H is significantly correlated with lower solar radiation intensity in East Asia and infer that MFSD12 was under selection in East Asians, probably after their split from Europeans.