Search

We consider a run-and-tumble particle on a half-line with an absorbing target at the origin. The particle has an internal velocity state that switches between two fixed values at Poisson-distributed times. The particles evolves according to an overdamped Langevin dynamics, with a periodic force field, such that every in a given period interval are accessible to the particle. The Laplace transform of the backward Fokker--Planck equation satisfied by the survival probability of the particle yields systems of equation satisfied by the moments of the first-passage time of the particle at the origin. The mean first passage time has already been calculated assuming the particle reaches the origin in finite time almost surely. We calculate the probability that the particle reaches the origin at finite time, given its initial position and velocity. We obtain an integral condition on the force, under which the exit probability is less than one. If the force field is a constant drift, the condition is satisfied by positive values of the drift. The mean first-passage time of the particle (averaged over trajectories that reach the origin) is obtained in closed form as a consequence in both regimes., Comment: 35 pages, 2 figures

A non-conserving zero-range process with extensive creation, annihilation and hopping rates is subjected to local resetting. The model is formulated on a large, fully-connected network of states. The states are equipped with a (bounded) fitness level: particles are added to each state at a rate proportional to the fitness level of the state. Moreover, particles are annihilated at a constant rate, and hop at a fixed rate to a uniformly-drawn state in the network. This model has been interpreted in terms of population dynamics: the fitness is the reproductive fitness in a haploid population, and the hopping process models mutation. It has also been interpreted as a model of network growth with a fixed set of nodes (in which particles occupying a state are interpreted as links pointing to this state). In the absence of resetting, the model is known to reach a steady state, which in a certain limit may exhibit a condensate at maximum fitness. If the model is subjected to global resetting by annihilating all particles at Poisson-distributed times, there is no condensation in the steady state. If the system is subjected to local resetting, the occupation numbers of each state are reset to zero at independent random times. These times are distributed according to a Poisson process whose rate (the resetting rate) depends on the fitness. We derive the evolution equation satisfied by the probability law of the occupation numbers. We calculate the average occupation numbers in the steady state. The existence of a condensate is found to depend on the local behavior of the resetting rate at maximum fitness: if the resetting rate vanishes at least linearly at high fitness, a condensate appears at maximum fitness in the limit where the sum of the annihilation and hopping rates is equal to the maximum fitness., Comment: 26 pages, 1 figure; v2: 27 pages, some clarifications, typos corrected, references added; v3: 28 pages, some more clarifications and corrections, references added

Renewal processes are zero-dimensional processes defined by independent intervals of time between zero crossings of a random walker. We subject renewal processes them to stochastic resetting by setting the position of the random walker to the origin at Poisson-distributed time with rate $r$. We introduce an additional parameter, the probability $\beta$ of keeping the sign state of the system at resetting time. Moreover, we introduce a trap at the origin, which absorbs the process with a fixed probability at each zero crossing. We obtain the mean lifetime of the process in closed form. For time intervals drawn from a L\'evy stable distribution of parameter $\theta$, the mean lifetime is finite for every positive value of the resetting rate, but goes to infinity when $r$ goes to zero. If the sign-keeping probability $\beta$ is higher than a critical level $\beta_c(\theta)$ (and strictly lower than $1$), the mean lifetime exhibits two extrema as a function of the resetting rate. Moreover, it goes to zero as $r^{-1}$ when $r$ goes to infinity. On the other hand, there is a single minimum if $\beta$ is set to one., Comment: 20 pages, LaTeX; v2: misprints corrected, references added; v3: more misprints corrected, appendices added

The voter model is a toy model of consensus formation based on nearest-neighbor interactions. A voter sits at each vertex in a hypercubic lattice (of dimension $d$) and is in one of two possible opinion states. The opinion state of each voter flips randomly, at a rate proportional to the fraction of the nearest neighbors that disagree with the voter. If the voters are initially independent and undecided, the model is known to lead to a consensus if and only if $d\leq 2$. In this paper the model is subjected to stochastic resetting: the voters revert independently to their initial opinion according to a Poisson process of fixed intensity (the resetting rate). This resetting prescription induces kinetic equations for the average opinion state and for the two-point function of the model. For initial conditions consisting of undecided voters except for one decided voter at the origin, the one-point function evolves as the probability of presence of a diffusive random walker on the lattice, whose position is stochastically reset to the origin. The resetting prescription leads to a non-equilibrium steady state. For an initial state consisting of independent undecided voters, the density of domain walls in the steady state is expressed in closed form as a function of the resetting rate. This function is differentiable at zero if and only if $d\geq 5$., Comment: 17 pages, 1 figure; V2: correction to Section 4.3, typos corrected, references added; V3: 25 pages, two-point function corrected, figures and numerical simulations added, generalization to higher dimensions added; V4: 38 pages, appendices added, references added; V5: 48 pages, typos; corrected, clarifications added; V6: references added, Fig. 1 added

We consider a random walker on a ring, subjected to resetting at Poisson-distributed times to the initial position (the walker takes the shortest path along the ring to the initial position at resetting times). In the case of a Brownian random walker the mean first-completion time of a turn is expressed in closed form as a function of the resetting rate. The value is shorter than in the ordinary process if the resetting rate is low enough. Moreover, the mean first-completion time of a turn can be minimised in the resetting rate. At large time the distribution of winding numbers does not reach a steady state, which is in contrast with the non-compact case of a Brownian particle under resetting on the real line. The mean total number of turns (and the variance of the net number of turns) grow linearly with time, with a proportionality constant equal to the inverse of the mean first-completion time of a turn., Comment: 16 pages, 2 figures; V2: typos corrected; V3: more typos corrected, references added

The model of binary aggregation with constant kernel is subjected to stochastic resetting: aggregates of any size explode into monomers at independent stochastic times. These resetting times are Poisson distributed, and the rate of the process is called the resetting rate. The master equation yields a Bernoulli-type equation in the generating function of the concentration of aggregates of any size, which can be solved exactly. This resetting prescription leads to a non-equilibrium steady state for the densities of aggregates, which is a function of the size of the aggregate, rescaled by a function of the resetting rate. The steady-state density of aggregates of a given size is maximised if the resetting rate is set to the quotient of the aggregation rate by the size of the aggregate (minus one)., Comment: 11 pages, LaTeX; v2: typos corrected, clarifications and references added; v3: more typos corrected; v4: 19 pages, 3 figures added, discussion expanded, more typos corrected, accepted version

We propose a model of run-and-tumble particles (RTPs) on a line with a fertile site at the origin. After going through the fertile site, a run-and-tumble particle gives rise to new particles until it flips direction. The process of creation of new particles is modelled by a fertility function (of the distance to the fertile site), multiplied by a fertility rate. If the initial conditions correspond to a single RTP with even probability density, the system is parity-invariant. The equations of motion can be solved in the Laplace domain, in terms of the density of right-movers at the origin. At large time, this density is shown to grow exponentially, at a rate that depends only on the fertility function and fertility rate. Moreover, the total density of RTPs (divided by the density of right-movers at the origin), reaches a stationary state that does not depend on the initial conditions, and presents a local minimum at the fertile site., Comment: 28 pages, LaTeX; v2: typos corrected; v3: Eq. 81 corrected, more typos corrected; v4: more typos corrected, references added; v5: more typos corrected; v6: more typos corrected, reference added, discussion expanded, accepted version

We consider the Larkin model of a directed polymer with Gaussian-distributed random forces, with the addition of a resetting process whereby the transverse position of the end-point of the polymer is reset to zero with constant rate $r$. We express the average over disorder of the mean time to absorption by an absorbing target at a fixed value of the transverse position. Thanks to the independence properties of the distribution of the random forces, this expression is analogous to the mean time to absorption for a diffusive particle under resetting, which possesses a single minimum at an optimal value $r^\ast$ of the resetting rate . Moreover, the mean time to absorption can be expanded as a power series of the amplitude of the disorder, around the value $r^\ast$ of the resetting rate. We obtain the susceptibility of the optimal resetting rate to disorder in closed form, and find it to be positive., Comment: 23 pages, LaTeX; v2: typos corrected; v3: discussion added, published version

The zero-temperature limit of the backgammon model under resetting is studied. The model is a balls-in-boxes model whose relaxation dynamics is governed by the density of boxes containing just one particle. As these boxes become rare at large times, the model presents an entropy barrier. As a preliminary step, a related model with faster relaxation, known to be mapped to a symmetric random walk, is studied by mapping recent results on diffusion with resetting onto the balls-in-boxes problem. Diffusion with an absorbing target at the origin (and diffusion constant equal to one), stochastically reset to the unit position, is a continuum approximation to the dynamics of the balls-in-boxes model, with resetting to a configuration maximising the number of boxes containing just one ball. In the limit of a large system, the relaxation time of the balls-in-boxes model under resetting is finite. The backgammon model subject to a constant resetting rate is then studied using an adiabatic approximation., Comment: 16 pages, LaTeX; V2, typos corrected; V3, clarifications added, references added, more typos crrected, accepted version

We consider a non-conserving zero-range process with hopping rate proportional to the number of particles at each site. Particles are added to the system with a site-dependent creation rate, and vanish with a uniform annihilation rate. On a fully-connected lattice with a large number of sites, the mean-field geometry leads to a negative binomial law for the number of particles at each site, with parameters depending on the hopping, creation and annihilation rates. This model can be mapped to population dynamics (if the creation rates are reproductive fitnesses in a haploid population, and the hopping rate is the mutation rate). It can also be mapped to a Bianconi--Barab\'asi model of a growing network with random rewiring of links (if creation rates are the rates of acquisition of links by nodes, and the hopping rate is the rewiring rate). The steady state has recently been worked out and gives rise to occupation numbers that reproduce Kingman's house-of-cards model of selection and mutation. In this paper we solve the master equation using a functional method, which yields integral equations satisfied by the occupation numbers. The occupation numbers are shown to forget initial conditions at an exponential rate that decreases linearly with the fitness level. Moreover, they can be computed exactly in the Laplace domain, which allows to obtain the steady state of the system under resetting. The result modifies the house-of-cards result by simply adding a skewed version of the initial conditions, and by adding the resetting rate to the hopping rate., Comment: 21 pages; v2: Section 2 clarified, published version

We consider a non-conserving zero-range process with hopping rate proportional to the number of particles at each site. Particles are added to the system with a site-dependent creation rate, and removed from the system with a uniform annihilation rate. On a fully-connected lattice with a large number of sites, the mean-field geometry leads to a negative binomial law for the number of particles at each site, with parameters depending on the hopping, creation and annihilation rates. This model of particles is mapped to a model of population dynamics: the site label is interpreted as a level of fitness, the site-dependent creation rate is interpreted as a selection function, and the hopping process is interpreted as the introduction of mutants. In the limit of large density, the fraction of the total population occupying each site approaches the limiting distribution in the house-of-cards model of selection-mutation, introduced by Kingman. A single site can be occupied by a macroscopic fraction of the particles if the mutation rate is below a critical value (which matches the critical value worked out in the house-of-cards model). This feature generalises to classes of selection functions that increase sufficiently fast at high fitness. The process can be mapped to a model of evolving networks, inspired by the Bianconi--Barab\'asi model, but involving a large and fixed set of nodes. Each node forms links at a rate biased by its fitness, moreover links are destroyed at a uniform rate, and redirected at a certain rate. If this redirection rate matches the mutation rate, the number of links pointing to nodes of a given fitness level is distributed as the numbers of particles in the non-conserving zero-range process. There is a finite critical redirection rate if the density of quenched fitnesses goes to zero sufficiently fast at high fitness., Comment: 16 pages, 2 figures; V2, typos corrected; V3, more typos corrected; V4, new Section 6 with application to a model evolving networks (accepted version)

The wiring diagram of the mouse brain has recently been mapped at a mesoscopic scale in the Allen Mouse Brain Connectivity Atlas. Axonal projections from brain regions were traced using green fluoresent proteins. The resulting data were registered to a common three-dimensional reference space. They yielded a matrix of connection strengths between 213 brain regions. Global features such as closed loops formed by connections of similar intensity can be inferred using tools from persistent homology. We map the wiring diagram of the mouse brain to a simplicial complex (filtered by connection strengths). We work out generators of the first homology group. Some regions, including nucleus accumbens, are connected to the entire brain by loops, whereas no region has non-zero connection strength to all brain regions. Thousands of loops go through the isocortex, the striatum and the thalamus. On the other hand, medulla is the only major brain compartment that contains more than 100 loops., Comment: 19 pages, 6 figures, 4 tables; V2: typos corrected, references added; V3: more typos corrected

The coordinate Bethe Ansatz solution of the log-gamma polymer is extended to boundary conditions with one fixed end and the other attached to one half of a one-dimensional lattice. The large-time limit is studied using a saddle-point approximation,and the cumulative distribution function of the rescaled free energy of a long polymer is expressed as a Fredholm determinant. Scaling limits of the kernel are identified, leading to a crossover from the GUE to the GOE Tracy--Widom distributions. The continuum limit reproduces the crossover from droplet to flat initial conditions of the Kardar--Parisi--Zhang equation., Comment: 20 pages

The discrete polymer model with random Boltzmann weights with homogeneous inverse gamma distribution, introduced by Sepp\"al\"ainen, is studied in the case of a polymer with one fixed and one free end. The model with two fixed ends has been integrated by Thiery and Le Doussal, using coordinate Bethe Ansatz techniques and an analytic-continuation prescription. The probability distribution of the free energy has been obtained through the replica method, even though the moments of the partition sum do not exist at all orders due to the fat tail in the distribution of Boltzmann weights. To extend this approach to the polymer with one free end, we argue that the contribution to the partition sums in the thermodynamic limit is localised on parity-invariant string states. This situation is analogous to the case of the continuum polymer with one free end, related to the Kardar--Parisi--Zhang equation with flat boundary conditions and solved by Le Doussal and Calabrese. The expansion of the generating function of the partition sum in terms of numbers of strings can also be transposed to the log-gamma polymer model, with the induced Fredholm determinant structure. We derive the large-time limit of the rescaled cumulative distribution function, and relate it to the GOE Tracy--Widom distribution. The derivation is conjectural in the sense that it assumes completeness of a family of string states (and expressions of their norms already used in the fixed-end problem) and extends heuristically the order of moments of the partition sum to the complex plane., Comment: 29 pages, no figures; v2: 36 pages, 4 figures, section 7 added (numerical simulations), version accepted in JSTAT

A centipede made of $N$ quantum walkers on a one-dimensional lattice is considered. The distance between two consecutive legs is either one or two lattice spacings, and a global constraint is imposed: the maximal distance between the first and last leg is $N+1$. This is the strongest global constraint compatible with walking. For an initial value of the wave function corresponding to a localized configuration at the origin, the probability law of the first leg of the centipede can be expressed in closed form in terms of Bessel functions. The dispersion relation and the group velocities are worked out exactly. Their maximal group velocity goes to zero when $N$ goes to infinity, which is in contrast with the behaviour of group velocities of quantum centipedes without global constraint, which were recently shown by Krapivsky, Luck and Mallick to give rise to ballistic spreading of extremal wave-front at non-zero velocity in the large-$N$ limit. The corresponding Hamiltonians are implemented numerically, based on a block structure of the space of configurations corresponding to compositions of the integer $N$. The growth of the maximal group velocity when the strong constraint is gradually relaxed is explored, and observed to be linear in the density of gaps allowed in the configurations. Heuristic arguments are presented to infer that the large-$N$ limit of the globally constrained model can yield finite group velocities provided the allowed number of gaps is a finite fraction of $N$., Comment: 8 pages, LaTeX; v2: 18 pages, 4 figures, extension to more general configurations, numerical checks added

The Allen Brain Atlas (ABA) of the adult mouse consists of digitized expression profiles of thousands of genes in the mouse brain, co-registered to a common three-dimensional template (the Allen Reference Atlas). This brain-wide, genome-wide data set has triggered a renaissance in neuroanatomy. Its voxelized version (with cubic voxels of side 200 microns) can be analyzed on a desktop computer using MATLAB. On the other hand, brain cells exhibit a great phenotypic diversity (in terms of size, shape and electrophysiological activity), which has inspired the names of some well-studied cell types, such as granule cells and medium spiny neurons. However, no exhaustive taxonomy of brain cells is available. A genetic classification of brain cells is under way, and some cell types have been characterized by their transcriptome profiles. However, given a cell type characterized by its transcriptome, it is not clear where else in the brain similar cells can be found. The ABA can been used to solve this region-specificity problem in a data-driven way: rewriting the brain-wide expression profiles of all genes in the atlas as a sum of cell-type-specific transcriptome profiles is equivalent to solving a quadratic optimization problem at each voxel in the brain. However, the estimated brain-wide densities of 64 cell types published recently were based on one series of co-registered coronal in situ hybridization (ISH) images per gene, whereas the online ABA contains several image series per gene, including sagittal ones. In the presented work, we simulate the variability of cell-type densities in a Monte Carlo way by repeatedly drawing a random image series for each gene and solving optimization problems. This yields error bars on the region-specificity of cell types., Comment: Prepared for the International Conference on Mathematical Modeling in Physical Sciences, 5th-8th June 2015, Mykonos Island, Greece

The Allen Atlas of the adult mouse brain is a brain-wide, genome-wide data set that has been made available online, triggering a renaissance in neuroanatomy. In particular, it has been used to define brain regions in a computational, data-driven way, and to estimate the region-specificity of cell types characterized independently by their transcriptional activity. However, these results were based on one series of co-registered (coronal) ISH image series per gene, whereas the online ABA contains several image series per genes, including sagittal ones. Since the sagittal series cover mostly the left hemisphere, we can simulate the variability of results by repeatedly drawing a random image series for each gene and restricting the computation to the left hemisphere. This gives rise to an estimate of error bars on the results of computational neuroanatomy., Comment: prepared for the International Conference on Mathematical Modeling in Physical Sciences June 5-8, 2015, Mykonos Island, Greece

Stereotactic targeting is a commonly used technique for performing injections in the brains of mice and other animals. The most common method for targeting stereoscopic injections uses the skull indentations bregma and lambda as reference points and is limited in its precision by factors such as skull curvature and individual variation, as well as an incomplete correspondence between skull landmarks and brain locations. In this software tool, a 3D laser scan of the mouse skull is taken in vitro and registered onto a reference skull using a point cloud matching algorithm, and the parameters of the transformation are used to position a glass pipette to place tracer injections. The software was capable of registering sample skulls with less than 100 micron error, and was able to target an injection in a mouse with error of roughly 500 microns. These results indicate that using skull scan registration has the potential to be widely applicable in automating stereotactic targeting of tracer injections., Comment: 8 pages

Two cliques of genes identified computationally for their high co-expression in the mouse brain according to the Allen Brain Atlas, and for their enrichment in genes related to autism spectrum disorder, have recently been shown to be highly co-expressed in the cerebellar cortex, compared to what could be expected by chance. Moreover, the expression of these cliques of genes is not homogeneous across the cerebellum, and it has been noted that their gene expression pattern seems to highlight the granular layer. However, this observation was only made by eye, and recent advances in computational neuroanatomy allow to rank cell types in the mouse brain (characterized by their transcriptome profiles) according to the similarity between their density profiles and the expression profiles of the cliques. We establish by Monte Carlo simulation that with probability at least 99 percent, the expression profiles of the two cliques are more similar to the density profile of granule cells than 99 percent of the expression of cliques containing the same number of genes (Purkinje cells also score above 99 percent in one of the cliques). Thresholding the expression profiles shows that the signal is more intense in the granular layer., Comment: 13 pages, 3 figures

The voxelized Allen Atlas of the adult mouse brain (at a resolution of 200 microns) has been used in [arXiv:1303.0013] to estimate the region-specificity of 64 cell types whose transcriptional profile in the mouse brain has been measured in microarray experiments. In particular, the model yields estimates for the brain-wide density of each of these cell types. We conduct numerical experiments to estimate the errors in the estimated density profiles. First of all, we check that a simulated thalamic profile based on 200 well-chosen genes can transfer signal from cerebellar Purkinje cells to the thalamus. This inspires us to sub-sample the atlas of genes by repeatedly drawing random sets of 200 genes and refitting the model. This results in a random distribution of density profiles, that can be compared to the predictions of the model. This results in a ranking of cell types by the overlap between the original and sub-sampled density profiles. Cell types with high rank include medium spiny neurons, several samples of cortical pyramidal neurons, hippocampal pyramidal neurons, granule cells and cholinergic neurons from the brain stem. In some cases with lower rank, the average sub-sample can have better contrast properties than the original model (this is the case for amygdalar neurons and dopaminergic neurons from the ventral midbrain). Finally, we add some noise to the cell-type-specific transcriptomes by mixing them using a scalar parameter weighing a random matrix. After refitting the model, we observe than a mixing parameter of $5\%$ leads to modifications of density profiles that span the same interval as the ones resulting from sub-sampling., Comment: 178 pages, 207 figures, 7 tables; v2: typos corrected; v3: more typos corrected, missing pseudo-code in section 3.5 written, image attachment changed in Figure 6, misuse of notation $r^{signal}$ corrected, conclusions unchanged

The Allen Atlas of the adult mouse brain is used to estimate the region-specificity of 64 cell types whose transcriptional profile in the mouse brain has been measured in microarray experiments. We systematically analyze the preliminary results presented in [arXiv:1111.6217], using the techniques implemented in the Brain Gene Expression Analysis toolbox. In particular, for each cell-type-specific sample in the study, we compute a brain-wide correlation profile to the Allen Atlas, and estimate a brain-wide density profile by solving a quadratic optimization problem at each voxel in the mouse brain. We characterize the neuroanatomical properties of the correlation and density profiles by ranking the regions of the left hemisphere delineated in the Allen Reference Atlas. We compare these rankings to prior biological knowledge of the brain region from which the cell-type-specific sample was extracted., Comment: V1: 88 pages, 37 figures, 43 tables; V2: 137 pages, new section discussing different fitting panels

We review quantitative methods and software developed to analyze genome-scale, brain-wide spatially-mapped gene-expression data. We expose new methods based on the underlying high-dimensional geometry of voxel space and gene space, and on simulations of the distribution of co-expression networks of a given size. We apply them to the Allen Atlas of the adult mouse brain, and to the co-expression network of a set of genes related to nicotine addiction retrieved from the NicSNP database. The computational methods are implemented in {\ttfamily{BrainGeneExpressionAnalysis}}, a Matlab toolbox available for download., Comment: 25 pages, 8 figures, accepted in Quantitative Biology (2012) 0002

The Allen Brain Atlas project (ABA) generated a genome-scale collection of gene-expression profiles using in-situ hybridization. These profiles were co-registered to the three-dimensional Allen Reference Atlas (ARA) of the adult mouse brain. A set of more than 4,000 such volumetric data are available for the full brain, at a resolution of 200 microns. These data are presented in a voxel-by-gene matrix. The ARA comes with several systems of annotation, hierarchical (40 cortical regions, 209 sub-cortical regions in the whole brain), or non-hierarchical (12 regions in the left hemisphere, with refinement into 94 regions, and cortical layers). The high-dimensional nature of this dataset and the possible connection between anatomy and gene expression pose challenges to data analysis. We developed the Brain Gene Expression Analysis Toolbox, whose functionalities include: determination of marker genes for brain regions, statistical analysis of brain-wide co-expression patterns, and the computation of brain-wide correlation maps with cell-type specific microarray data., Comment: 59 pages; v2: bug fixed on page 5, preface added; v3: more bugs fixed, cell-type-specific data released (chapter 5); v4: references added, chapter 2 expanded; v5: download link fixed, more typos corrected; v6: download link fixed

We use the Allen Gene Expression Atlas (AGEA) and the OMA ortholog dataset to investigate the evolution of mouse-brain neuroanatomy from the standpoint of the molecular evolution of brain-specific genes. For each such gene, using the phylogenetic tree for all fully sequenced species and the presence of orthologs of the gene in these species, we construct and assign a discrete measure of evolutionary age. The gene expression profile of all gene of similar age, relative to the average gene expression profile, distinguish regions of the brain that are over-represented in the corresponding evolutionary timescale. We argue that the conclusions one can draw on evolution of twelve major brain regions from such a molecular level analysis supplements existing knowledge of mouse brain evolution and introduces new quantitative tools, especially for comparative studies, when AGEA-like data sets for other species become available. Using the functional role of the genes representational of a certain evolutionary timescale and brain region we compare and contrast, wherever possible, our observations with existing knowledge in evolutionary neuroanatomy., Comment: 14 pages

The three-dimensional data-driven Allen Gene Expression Atlas of the adult mouse brain consists of numerized in-situ hybridization data for thousands of genes, co-registered to the Allen Reference Atlas. We propose quantitative criteria to rank genes as markers of a brain region, based on the localization of the gene expression and on its functional fitting to the shape of the region. These criteria lead to natural generalizations to sets of genes. We find sets of genes weighted with coefficients of both signs with almost perfect localization in all major regions of the left hemisphere of the brain, except the pallidum. Generalization of the fitting criterion with positivity constraint provides a lesser improvement of the markers, but requires sparser sets of genes., Comment: 6 pages, 5 figures

Much of the genome is expressed in the vertebrate brain, with individual genes exhibiting different spatially-varying patterns of expression. These variations are not independent, with pairs of genes exhibiting complex patterns of co-expression, such that two genes may be similarly expressed in one region, but differentially expressed in other regions. These correlations have been previously studied quantitatively, particularly for the gene expression atlas of the mouse brain, but the biological meaning of the co-expression patterns remains obscure. We propose a simple model of the co-expression patterns in terms of spatial distributions of underlying cell types. We establish the plausibility of the model in terms of a test set of cell types for which both the gene expression profiles and the spatial distributions are known., Comment: 41 pages; v2: typos corrected

We propose a quantitative method to estimate the statistical properties of sets of genes for which expression data are available and co-registered to a reference atlas of the brain. It is based on graph-theoretic properties of co-expression coefficients between pairs of genes. We apply this method to mouse genes from the Allen Gene Expression Atlas. Co-expression patterns of a list of several hundreds of genes related to addiction are analyzed, using ISH data produced for the mouse brain at the Allen Institute. It appears that large subsets of this set of genes are much more highly co-expressed than expected by chance., Comment: 11 pages, 3 figures; v2: 2 figures added, references added

Quantitative criteria are proposed to identify genes (and sets of genes) whose expression marks a specific brain region (or a set of brain regions). Gene-expression energies, obtained for thousands of mouse genes by numerization of in-situ hybridization images in the Allen Gene Expression Atlas, are used to test these methods in the mouse brain. Individual genes are ranked using integrals of their expression energies across brain regions. The ranking is generalized to sets of genes and the problem of optimal markers of a classical region receives a linear-algebraic solution. Moreover, the goodness of the fitting of the expression profile of a gene to the profile of a brain region is closely related to the co-expression of genes. The geometric interpretation of this fact leads to a quantitative criterion to detect markers of pairs of brain regions. Local properties of the gene-expression profiles are also used to detect genes that separate a given grain region from its environment., Comment: 26 pages, LaTeX

Given an atlas of the brain and a number of injections to be performed in order to map out the connections between parts of the brain, we propose an algorithm to compute the coordinates of the injections. The algorithm is designed to sample the brain in the most homogeneous way compatible with the separation of brain regions. It can be applied to other species for which a neuroanatomical atlas is available. The computation is tested on the annotation at a resolution of 25 microns corresponding to the Allen Reference Atlas, which is hierarchical and consists of 209 regions. The resulting injection coordinates are being used for the injection protocol of the Mouse Brain Architecture project. Due to its large size and layered structure, the cerebral cortex is treated in a separate algorithm, which is more adapted to its geometry., Comment: 13 pages, LaTeX

Fibrations of flux backgrounds by supersymmetric cycles are investigated. For an internal six-manifold M with static SU(2) structure and mirror \hat{M}, it is argued that the product M x \hat{M} is doubly fibered by supersymmetric three-tori, with both sets of fibers transverse to M and \hat{M}. The mirror map is then realized by T-dualizing the fibers. Mirror-symmetric properties of the fluxes, both geometric and non-geometric, are shown to agree with previous conjectures based on the requirement of mirror symmetry for Killing prepotentials. The fibers are conjectured to be destabilized by fluxes on generic SU(3)xSU(3) backgrounds, though they may survive at type-jumping points. T-dualizing the surviving fibers ensures the exchange of pure spinors under mirror symmetry., Comment: 30 pages, 3 figures, LaTeX; v2: references added

Various approaches to T-duality with NSNS three-form flux are reconciled. Non-commutative torus fibrations are shown to be the open-string version of T-folds. The non-geometric T-dual of a three-torus with uniform flux is embedded into a generalized complex six-torus, and the non-geometry is probed by D0-branes regarded as generalized complex submanifolds. The non-commutativity scale, which is present in these compactifications, is given by a holomorphic Poisson bivector that also encodes the variation of the dimension of the world-volume of D-branes under monodromy. This bivector is shown to exist in SU(3) x SU(3) structure compactifications, which have been proposed as mirrors to NSNS-flux backgrounds. The two SU(3)-invariant spinors are generically not parallel, thereby giving rise to a non-trivial Poisson bivector. Furthermore we show that for non-geometric T-duals, the Poisson bivector may not be decomposable into the tensor product of vectors., Comment: 25 pages, LaTeX; v2: typos corrected, references added

In generalized complex geometry, D-branes can be seen as maximally isotropic spaces and are thus in one-to-one correspondence with pure spinors. When considered on the sum of the tangent and cotangent bundles to the ambient space, all the branes are of the same dimension and the transverse scalars enter on par with the gauge fields; the split between the longitudinal and transverse directions is done in accordance with the type of the pure spinor corresponding to the given D-brane. We elaborate on the relation of this picture to the T-duality transformations and stability of D-branes. A discussion of tachyon condensation in the context of the generalized complex geometry is given, linking the description of D-branes as generalized complex submanifolds to their K-theoretic classification., Comment: 19 pages, LaTeX

This thesis is devoted to derivative corrections to the effective action of D-branes, and to mirror symmetry with D-branes. Series of derivative corrections first predicted by non-commutative gauge theory are completed by couplings between the metric and the gauge field. The result is interpreted as a deformation of the non-commutative gauge theory, whose structure survives. The derivation is applied to the tachyon field, whose potential is shown to be deformed by the very same corrections. Moreover, a prescription is given for the coupling of p-adic strings to a magnetic field, thus allowing to study p-adic solitons using non-commutative field-theory techniques. The link with topological D-branes is provided by the non-commutative description of D-branes in the B-model. The fibre bundles supported by the D-branes are still holomorphic in this description. Establishing this property involves the realization of D-branes as boundary conditions, within the framework of generalized complex geometry. This geometric framework is then used to describe mirror symmetry with D-branes on a Calabi--Yau manifold admitting a $T^3$-fibration. Two pure spinors, involved in the stability equations for topological D-branes, and modified by gauge fields, are exchanged, thus unifying Lagrangian and non-Lagrangian D-branes of the A-model as mirrors of stable D-branes of the B-model., Comment: PhD thesis, defended 27/05/2005; text in French, 98 pages

It has been argued recently that mirror symmetry exchanges two pure spinors characterizing a generic manifold with SU(3)-structure. We show how pure spinors are modified in the presence of topological D-branes, so that they are still exchanged by mirror symmetry. This exchange emerges from the fact that the modified pure spinors come out as moment maps for the symmetries of A and B-models. The modification by the gauge field is argued to ensure the inclusion into the mirror exchange of the A-model non-Lagrangian branes endowed with a non-flat connection. Treating the connection as a distribution on an ambient six-manifold, assumed to be T^3-fibered, the proposed mirror formula is established by fiberwise T-duality., Comment: 13 pages, LaTeX

Products defined in the context of noncommutative gauge theory allow for an interpolation between exact results on tachyon potentials at zero and large background B-fields. Techniques for computations of effective actions are transposed from the framework of gauge theory to the framework of boundary string field theory, resulting in deformations of the noncommutative tachyon potential by anomalous dimensions., Comment: 10 pages, LaTeX

A new term in the p-adic world-sheet action is proposed, which couples a constant B-field to the boundary of the world-sheet at disk level. The induced deformation of tachyon scattering amplitudes by star-products is derived. This is in agreement with the deformation of effective action postulated in recent investigations of noncommutative solitons in p-adic string theory., Comment: 6 pages, LaTeX; note and references added

It was recently suggested by A. Kapustin that turning on a $B$-field, and allowing some discrepancy between the left and and right-moving complex structures, must induce an identification of B-branes with holomorphic line bundles on a non-commutative complex torus. We translate the stability condition for the branes into this language and identify the stable topological branes with previously proposed non-commutative instanton equations. This involves certain topological identities whose derivation has become familiar in non-commutative field theory. It is crucial for these identities that the instantons are localized. We therefore explore the case of non-constant field strength, whose non-linearities are dealt with thanks to the rank-one Seiberg--Witten map., Comment: 12 pages, LaTeX

Derivative corrections to the Wess--Zumino couplings of open-string effective actions are computed at all orders in derivatives, taking the open-string metric into account. This leads to a set of deformed star-products beyond the Seiberg--Witten limit, and allows to reinterpret the couplings in terms of a deformed integration prescription along a Wilson line in the non-commutative set-up. Moreover, the recursive definition of the star-products induces deformations of U(1) non-commutative Yang--Mills theory., Comment: 9 pages; references added; a few clarifications and a consistency check added, version to appear in Phys. Lett. B

The boundary state formalism is used to confirm predictions from non-commutativity for the derivative corrections to the Dirac--Born--Infeld and Chern--Simons actions, at all orders in derivatives. As anticipated by S. Mukhi, the method applies by induction to every coupling in the Chern--Simons action. It is also used to derive the corrections to the Dirac--Born--Infeld action at quadratic order in the field strength., Comment: 16 pages, references added, a few points clarified. Completed list of references