Source author record

Boris I. Shraiman

Boris I. Shraiman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2020arXiv

Multi-defect Dynamics in Active Nematics

Recent experiments and numerical studies have drawn attention to the dynamics of active nematics. Two-dimensional active nematics flow spontaneously and exhibit spatiotemporal chaotic flows with proliferation of topological defects in the nematic texture. It has been proposed that the dynamics of active nematics can be understood in terms of the dynamics of interacting defects, propelled by active stress. Previous work has derived effective equations of motion for individual defects as quasi-particles moving in the mean field generated by other defects, but an effective theory governing multi-defect dynamics has remained out of reach. In this paper, we examine the dynamics of 2D active nematics in the limit of strong order and overdamped compressible flow. The activity-induced defect dynamics is formulated as a perturbation of the manifold of quasi-static nematic textures explicitly parameterized by defect positions. This makes it possible to derive a set of coupled ordinary differential equations governing defect (and therefore texture) dynamics. Interestingly, because of the non-orthogonality of textures associated with individual defects, their motion is coupled through a position dependent ``collective mobility" matrix. In addition to the familiar active self-propulsion of the $+1/2$ defect, we obtain new collective effects of activity that can be interpreted in terms of non-central and non-reciprocal interactions between defects.

preprint2018arXiv

Geometry of epithelial cells provides a robust method for image based inference of stress within tissues

Cellular mechanics plays an important role in epithelial morphogenesis, a process wherein cells reshape and rearrange to produce tissue-scale deformations. However, the study of tissue-scale mechanics is impaired by the difficulty of direct measurement of stress in-vivo. Alternative, image-based inference schemes aim to estimate stress from snapshots of cellular geometry but are challenged by sensitivity to fluctuations and measurement noise as well as the dependence on boundary conditions. Here we overcome these difficulties by introducing a new variational approach - the Geometrical Variation Method (GVM) - which exploits the fundamental duality between stress and cellular geometry that exists in the state of mechanical equilibrium of discrete mechanical networks that approximate cellular tissues. In the Geometrical Variation Method, the two dimensional apical geometry of an epithelial tissue is approximated by a 2D tiling with Circular Arc Polygons (CAP) in which the arcs represent intercellular interfaces defined by the balance of local line tension and pressure differentials between adjacent cells. We take advantage of local constraints that mechanical equilibrium imposes on CAP geometry to define a variational procedure that extracts the best fitting equilibrium configuration from images of epithelial monolayers. The GVM-based stress inference algorithm has been validated by the comparison of the predicted cellular and mesoscopic scale stress and measured myosin II patterns in the epithelial tissue during Drosophila embryogenesis. GVM prediction of mesoscopic stress tensor correlates at the 80% level with the measured myosin distribution and reveals that most of the myosin II activity is involved in a static internal force balance within the epithelial layer. Lastly, this study provides a practical method for non-destructive estimation of stress in live epithelial tissues.

preprint2016arXiv

Leaf growth is conformal

Growth pattern dynamics lie at the heart of morphogenesis. Here, we investigate the growth of plant leaves. We compute the conformal transformation that maps the contour of a leaf at a given stage onto the contour of the same leaf at a later stage. Based on the mapping we predict the local displacement field in the leaf blade and find it to agree with the experimentally measured displacement field to 92%. This approach is applicable to any two-dimensional system with locally isotropic growth, enabling the deduction of the whole growth field just from observation of the tissue contour.

preprint2015arXiv

Inferring Epigenetic Dynamics from Kin Correlations

Populations of isogenic embryonic stem cells or clonal bacteria often exhibit extensive phenotypic heterogeneity which arises from stochastic intrinsic dynamics of cells. The internal state of the cell can be transmitted epigenetically in cell division, leading to correlations in the phenotypic states of cells related by descent. Therefore, a phenotypic snapshot of a collection of cells with known genealogical structure, contains information on phenotypic dynamics. Here we use a model of phenotypic dynamics on a genealogical tree to define an inference method which allows to extract an approximate probabilistic description of phenotypic dynamics based on measured correlations as a function of the degree of kinship. The approach is tested and validated on the example of Pyoverdine dynamics in P. aeruginosa colonies. Interestingly, we find that correlations among pairs and triples of distant relatives have a simple but non-trivial structure indicating that observed phenotypic dynamics on the genealogical tree is approximately conformal - a symmetry characteristic of critical behavior in physical systems. Proposed inference method is sufficiently general to be applied in any system where lineage information is available.

preprint2015arXiv

Prediction, dynamics, and visualization of antigenic phenotypes of seasonal influenza viruses

Human seasonal influenza viruses evolve rapidly, enabling the virus population to evade immunity and re-infect previously infected individuals. Antigenic properties are largely determined by the surface glycoprotein hemagglutinin (HA) and amino acid substitutions at exposed epitope sites in HA mediate loss of recognition by antibodies. Here, we show that antigenic differences measured through serological assay data are well described by a sum of antigenic changes along the path connecting viruses in a phylogenetic tree. This mapping onto the tree allows prediction of antigenicity from HA sequence data alone. The mapping can further be used to make predictions about the makeup of the future seasonal influenza virus population, and we compare predictions between models with serological and sequence data. To make timely model output readily available, we developed a web browser based application that visualizes antigenic data on a continuously updated phylogeny.

preprint2014arXiv

Predicting evolution from the shape of genealogical trees

Given a sample of genome sequences from an asexual population, can one predict its evolutionary future? Here we demonstrate that the branching patterns of reconstructed genealogical trees contains information about the relative fitness of the sampled sequences and that this information can be used to predict successful strains. Our approach is based on the assumption that evolution proceeds by accumulation of small effect mutations, does not require species specific input and can be applied to any asexual population under persistent selection pressure. We demonstrate its performance using historical data on seasonal influenza A/H3N2 virus. We predict the progenitor lineage of the upcoming influenza season with near optimal performance in 30% of cases and make informative predictions in 16 out of 19 years. Beyond providing a tool for prediction, our ability to make informative predictions implies persistent fitness variation among circulating influenza A/H3N2 viruses.

preprint2013arXiv

Coalescence, genetic diversity in sexual populations under selection

In sexual populations, selection operates neither on the whole genome, which is repeatedly taken apart and reassembled by recombination, nor on individual alleles that are tightly linked to the chromosomal neighborhood. The resulting interference between linked alleles reduces the efficiency of selection and distorts patterns of genetic diversity. Inference of evolutionary history from diversity shaped by linked selection requires an understanding of these patterns. Here, we present a simple but powerful scaling analysis identifying the unit of selection as the genomic "linkage block" with a characteristic length determined in a self-consistent manner by the condition that the rate of recombination within the block is comparable to the fitness differences between different alleles of the block. We find that an asexual model with the strength of selection tuned to that of the linkage block provides an excellent description of genetic diversity and the site frequency spectra when compared to computer simulations. This linkage block approximation is accurate for the entire spectrum of strength of selection and is particularly powerful in scenarios with many weakly selected loci. The latter limit allows us to characterize coalescence, genetic diversity, and the speed of adaptation in the infinitesimal model of quantitative genetics.

preprint2012arXiv

Emergence of clones in sexual populations

In sexual population, recombination reshuffles genetic variation and produces novel combinations of existing alleles, while selection amplifies the fittest genotypes in the population. If recombination is more rapid than selection, populations consist of a diverse mixture of many genotypes, as is observed in many populations. In the opposite regime, which is realized for example in the facultatively sexual populations that outcross in only a fraction of reproductive cycles, selection can amplify individual genotypes into large clones. Such clones emerge when the fitness advantage of some of the genotypes is large enough that they grow to a significant fraction of the population despite being broken down by recombination. The occurrence of this "clonal condensation" depends, in addition to the outcrossing rate, on the heritability of fitness. Clonal condensation leads to a strong genetic heterogeneity of the population which is not adequately described by traditional population genetics measures, such as Linkage Disequilibrium. Here we point out the similarity between clonal condensation and the freezing transition in the Random Energy Model of spin glasses. Guided by this analogy we explicitly calculate the probability, Y, that two individuals are genetically identical as a function of the key parameters of the model. While Y is the analog of the spin-glass order parameter, it is also closely related to rate of coalescence in population genetics: Two individuals that are part of the same clone have a recent common ancestor.

preprint2012arXiv

Fluctuations of fitness distributions and the rate of Muller's ratchet

The accumulation of deleterious mutations is driven by rare fluctuations which lead to the loss of all mutation free individuals, a process known as Muller's ratchet. Even though Muller's ratchet is a paradigmatic process in population genetics, a quantitative understanding of its rate is still lacking. The difficulty lies in the nontrivial nature of fluctuations in the fitness distribution which control the rate of extinction of the fittest genotype. We address this problem using the simple but classic model of mutation selection balance with deleterious mutations all having the same effect on fitness. We show analytically how fluctuations among the fittest individuals propagate to individuals of lower fitness and have a dramatically amplified effects on the bulk of the population at a later time. If a reduction in the size of the fittest class reduces the mean fitness only after a delay, selection opposing this reduction is also delayed. This delayed restoring force speeds up Muller's ratchet. We show how the delayed response can be accounted for using a path integral formulation of the stochastic dynamics and provide an expression for the rate of the ratchet that is accurate across a broad range of parameters.

preprint2011arXiv

Collective and single cell behavior in epithelial contact inhibition

Control of cell proliferation is a fundamental aspect of tissue physiology central to morphogenesis, wound healing and cancer. Although many of the molecular genetic factors are now known, the system level regulation of growth is still poorly understood. A simple form of inhibition of cell proliferation is encountered in vitro in normally differentiating epithelial cell cultures and is known as "contact inhibition". The study presented here provides a quantitative characterization of contact inhibition dynamics on tissue-wide and single cell levels. Using long-term tracking of cultured MDCK cells we demonstrate that inhibition of cell division in a confluent monolayer follows inhibition of cell motility and sets in when mechanical constraint on local expansion causes divisions to reduce cell area. We quantify cell motility and cell cycle statistics in the low density confluent regime and their change across the transition to epithelial morphology which occurs with increasing cell density. We then study the dynamics of cell area distribution arising through reductive division, determine the average mitotic rate as a function of cell size and demonstrate that complete arrest of mitosis occurs when cell area falls below a critical value. We also present a simple computational model of growth mechanics which captures all aspects of the observed behavior. Our measurements and analysis show that contact inhibition is a consequence of mechanical interaction and constraint rather than interfacial contact alone, and define quantitative phenotypes that can guide future studies of molecular mechanisms underlying contact inhibition.

preprint2011arXiv

Genetic Draft and Quasi-Neutrality in Large Facultatively Sexual Populations

Large populations may contain numerous simultaneously segregating polymorphisms subject to natural selection. Since selection acts on individuals whose fitness depends on many loci, different loci affect each other's dynamics. This leads to stochastic fluctuations of allele frequencies above and beyond genetic drift - an effect known as genetic draft. Since recombination disrupts associations between alleles, draft is strong when recombination is rare. Here, we study a facultatively outcrossing population in a regime where the frequency of out-crossing and recombination, r, is small compared to the characteristic scale of fitness differences σ. In this regime, fit genotypes expand clonally, leading to large fluctuations in the number of recombinant offspring genotypes. The power law tail in the distribution of the latter makes it impossible to capture the dynamics of draft by an effective neutral model. Instead, we find that the fixation time of a neutral allele increases only slowly with the population size but depends sensitively on the ratio r/σ. The efficacy of selection is reduced dramatically and alleles behave "quasi-neutrally" even for Ns>> 1, provided that |s|< s_c, where s_c depends strongly on r/σ, but only weakly on population size N. In addition, the anomalous fluctuations due to draft change the spectrum of (quasi)-neutral alleles from f(ν)\sim 1/ν, corresponding to drift, to \sim1/ν^2. Finally, draft accelerates the rate of two step adaptations through deleterious intermediates.

preprint2011arXiv

Mechanical Stress Inference for Two Dimensional Cell Arrays

Many morphogenetic processes involve mechanical rearrangement of epithelial tissues that is driven by precisely regulated cytoskeletal forces and cell adhesion. The mechanical state of the cell and intercellular adhesion are not only the targets of regulation, but are themselves likely signals that coordinate developmental process. Yet, because it is difficult to directly measure mechanical stress {\it in vivo} on sub-cellular scale, little is understood about the role of mechanics of development. Here we present an alternative approach which takes advantage of the recent progress in live imaging of morphogenetic processes and uses computational analysis of high resolution images of epithelial tissues to infer relative magnitude of forces acting within and between cells. We model intracellular stress in terms of bulk pressure and interfacial tension, allowing these parameters to vary from cell to cell and from interface to interface. Assuming that epithelial cell layers are close to mechanical equilibrium, we use the observed geometry of the two dimensional cell array to infer interfacial tensions and intracellular pressures. Here we present the mathematical formulation of the proposed Mechanical Inverse method and apply it to the analysis of epithelial cell layers observed at the onset of ventral furrow formation in the {\it Drosophila} embryo and in the process of hair-cell determination in the avian cochlea. The analysis reveals mechanical anisotropy in the former process and mechanical heterogeneity, correlated with cell differentiation, in the latter process. The method opens a way for quantitative and detailed experimental tests of models of cell and tissue mechanics.

preprint2011arXiv

Rare beneficial mutations can halt Muller's ratchet

The vast majority of mutations are deleterious, and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet turns into a mutational "meltdown" leading to a rapid degradation of population fitness. Evidently, the long term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Here we propose that the stable evolutionary state of a population in a static environment is a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that this state exists for any population size N and mutation rate $U$. Assuming that beneficial and deleterious mutations have the same fitness effect s, we calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low εsuffices to maintain stability, even in small populations in the face of high mutation rates and weak selection. This may explain the maintenance of mitochondria and other asexual genomes, and has implications for the expected statistics of genetic diversity in these populations.

preprint2011arXiv

Rate of Adaptation in Large Sexual Populations

Adaptation often involves the acquisition of a large number of genomic changes which arise as mutations in single individuals. In asexual populations, combinations of mutations can fix only when they arise in the same lineage, but for populations in which genetic information is exchanged, beneficial mutations can arise in different individuals and be combined later. In large populations, when the product of the population size N and the total beneficial mutation rate U_b is large, many new beneficial alleles can be segregating in the population simultaneously. We calculate the rate of adaptation, v, in several models of such sexual populations and show that v is linear in NU_b only in sufficiently small populations. In large populations, v increases much more slowly as log NU_b. The prefactor of this logarithm, however, increases as the square of the recombination rate. This acceleration of adaptation by recombination implies a strong evolutionary advantage of sex.

preprint2011arXiv

Statistical Genetics and Evolution of Quantitative Traits

The distribution and heritability of many traits depends on numerous loci in the genome. In general, the astronomical number of possible genotypes makes the system with large numbers of loci difficult to describe. Multilocus evolution, however, greatly simplifies in the limit of weak selection and frequent recombination. In this limit, populations rapidly reach Quasi-Linkage Equilibrium (QLE) in which the dynamics of the full genotype distribution, including correlations between alleles at different loci, can be parameterized by the allele frequencies. This review provides a simplified exposition of the concept and mathematics of QLE which is central to the statistical description of genotypes in sexual populations. We show how key results of Quantitative Genetics such as the generalized Fisher's "Fundamental Theorem", along with Wright's Adaptive Landscape, emerge within QLE from the dynamics of the genotype distribution. We then discuss under what circumstances QLE is applicable, and what the breakdown of QLE implies for the population structure and the dynamics of selection. Understanding of the fundamental aspects of multilocus evolution obtained through simplified models may be helpful in providing conceptual and computational tools to address the challenges arising in the studies of complex quantitative phenotypes of practical interest.