Marginal stability and traveling fronts in two-phase mixtures

Toward Predicting the Spatio-Temporal Dynamics of Alopecia Areata Lesions Using Partial Differential Equation Analysis

Hair loss in the autoimmune disease, alopecia areata (AA), is characterized by the appearance of circularly spreading alopecic lesions in seemingly healthy skin. The distinct spatial patterns of AA lesions form because the immune system attacks hair follicle cells that are in the process of producing hair shaft, catapults the mini-organs that produce hair from a state of growth (anagen) into an apoptosis-driven regression state (catagen), and causes major hair follicle dystrophy along with rapid hair shaft shedding. In this paper, we develop a model of partial differential equations (PDEs) to describe the spatio-temporal dynamics of immune system components that clinical and experimental studies show are primarily involved in the disease development. Global linear stability analysis reveals there is a most unstable mode giving rise to a pattern. The most unstable mode indicates a spatial scale consistent with results of the humanized AA mouse model of Gilhar et al. (Autoimmun Rev 15(7):726-735, 2016) for experimentally induced AA lesions. Numerical simulations of the PDE system confirm our analytic findings and illustrate the formation of a pattern that is characteristic of the spatio-temporal AA dynamics. We apply marginal linear stability analysis to examine and predict the pattern propagation.

Pattern formation in multiphase models of chemotactic cell aggregation

We develop a continuum model for the aggregation of cells cultured in a nutrient-rich medium in a culture well. We consider a 2D geometry, representing a vertical slice through the culture well, and assume that the cell layer depth is small compared with the typical lengthscale of the culture well. We adopt a continuum mechanics approach, treating the cells and culture medium as a two-phase mixture. Specifically, the cells and culture medium are treated as fluids. Additionally, the cell phase can generate forces in response to environmental cues, which include the concentration of a chemoattractant that is produced by the cells within the culture medium. The model leads to a system of coupled nonlinear partial differential equations for the volume fraction and velocity of the cell phase, the culture medium pressure and the chemoattractant concentration, which must be solved subject to appropriate boundary and initial conditions. To gain insight into the system, we consider two model reductions, appropriate when the cell layer depth is thin compared to the typical length scale of the culture well: a (simple) 1D and a (more involved) thin-film extensional flow reduction. By investigating the resulting systems of equations analytically and numerically, we identify conditions under which small amplitude perturbations to a homogeneous steady state (corresponding to a spatially uniform cell distribution) can lead to a spatially varying steady state (pattern formation). Our analysis reveals that the simpler 1D reduction has the same qualitative features as the thin-film extensional flow reduction in the linear and weakly nonlinear regimes, motivating the use of the simpler 1D modelling approach when a qualitative understanding of the system is required. However, the thin-film extensional flow reduction may be more appropriate when detailed quantitative agreement between modelling predictions and experimental data is desired. Furthermore, full numerical simulations of the two model reductions in regions of parameter space when the system is not close to marginal stability reveal significant differences in the evolution of the volume fraction and velocity of the cell phase, and chemoattractant concentration.

Differential Activity-Driven Instabilities in Biphasic Active Matter

Active stresses can cause instabilities in contractile gels and living tissues. Here we provide a generic hydrodynamic theory that treats these systems as a mixture of two phases of varying activity and different mechanical properties. We find that differential activity between the phases causes a uniform mixture to undergo a demixing instability. We follow the nonlinear evolution of the instability and characterize a phase diagram of the resulting patterns. Our study complements other instability mechanisms in mixtures driven by differential adhesion, differential diffusion, differential growth, and differential motion.