Carlos Villarreal

Single-cell and coupled GRN models of cell patterning in the Arabidopsis thaliana root stem cell niche

BackgroundRecent experimental work has uncovered some of the genetic components required to maintain the Arabidopsis thaliana root stem cell niche (SCN) and its structure. Two main pathways are involved. One pathway depends on the genes SHORTROOT and SCARECROW and the other depends on the PLETHORA genes, which have been proposed to constitute the auxin readouts. Recent evidence suggests that a regulatory circuit, composed of WOX5 and CLE40, also contributes to the SCN maintenance. Yet, we still do not understand how the niche is dynamically maintained and patterned or if the uncovered molecular components are sufficient to recover the observed gene expression configurations that characterize the cell types within the root SCN. Mathematical and computational tools have proven useful in understanding the dynamics of cell differentiation. Hence, to further explore root SCN patterning, we integrated available experimental data into dynamic Gene Regulatory Network (GRN) models and addressed if these are sufficient to attain observed gene expression configurations in the root SCN in a robust and autonomous manner.ResultsWe found that an SCN GRN model based only on experimental data did not reproduce the configurations observed within the root SCN. We developed several alternative GRN models that recover these expected stable gene configurations. Such models incorporate a few additional components and interactions in addition to those that have been uncovered. The recovered configurations are stable to perturbations, and the models are able to recover the observed gene expression profiles of almost all the mutants described so far. However, the robustness of the postulated GRNs is not as high as that of other previously studied networks.ConclusionsThese models are the first published approximations for a dynamic mechanism of the A. thaliana root SCN cellular pattering. Our model is useful to formally show that the data now available are not sufficient to fully reproduce root SCN organization and genetic profiles. We then highlight some experimental holes that remain to be studied and postulate some novel gene interactions. Finally, we suggest the existence of a generic dynamical motif that can be involved in both plant and animal SCN maintenance.

A Minimal Regulatory Network of Extrinsic and Intrinsic Factors Recovers Observed Patterns of CD4+ T Cell Differentiation and Plasticity.

CD4+ T cells orchestrate the adaptive immune response in vertebrates. While both experimental and modeling work has been conducted to understand the molecular genetic mechanisms involved in CD4+ T cell responses and fate attainment, the dynamic role of intrinsic (produced by CD4+ T lymphocytes) versus extrinsic (produced by other cells) components remains unclear, and the mechanistic and dynamic understanding of the plastic responses of these cells remains incomplete. In this work, we studied a regulatory network for the core transcription factors involved in CD4+ T cell-fate attainment. We first show that this core is not sufficient to recover common CD4+ T phenotypes. We thus postulate a minimal Boolean regulatory network model derived from a larger and more comprehensive network that is based on experimental data. The minimal network integrates transcriptional regulation, signaling pathways and the micro-environment. This network model recovers reported configurations of most of the characterized cell types (Th0, Th1, Th2, Th17, Tfh, Th9, iTreg, and Foxp3-independent T regulatory cells). This transcriptional-signaling regulatory network is robust and recovers mutant configurations that have been reported experimentally. Additionally, this model recovers many of the plasticity patterns documented for different T CD4+ cell types, as summarized in a cell-fate map. We tested the effects of various micro-environments and transient perturbations on such transitions among CD4+ T cell types. Interestingly, most cell-fate transitions were induced by transient activations, with the opposite behavior associated with transient inhibitions. Finally, we used a novel methodology was used to establish that T-bet, TGF-β and suppressors of cytokine signaling proteins are keys to recovering observed CD4+ T cell plastic responses. In conclusion, the observed CD4+ T cell-types and transition patterns emerge from the feedback between the intrinsic or intracellular regulatory core and the micro-environment. We discuss the broader use of this approach for other plastic systems and possible therapeutic interventions.

Exact surface impedance formulation of the Casimir force: application to spatially dispersive metals

We show that exact results are obtained for the calculation of Casimir forces between arbitrary materials using the concept of surface impedances, obtaining in a trivial way the force in the limit of perfect conductors and also Lifshitz formula in the limit of semi-infinite media. As an example we present a full and rigorous calculation of the Casimir force between two metallic half-spaces described by a hydrodynamic nonlocal dielectric response.

Gene Regulatory Network Models for Floral Organ Determination

Understanding how genotypes map unto phenotypes implies an integrative understanding of the processes regulating cell differentiation and morphogenesis, which comprise development. Such a task requires the use of theoretical and computational approaches to integrate and follow the concerted action of multiple genetic and nongenetic components that hold highly nonlinear interactions. Gene regulatory network (GRN) models have been proposed to approach such task. GRN models have become very useful to understand how such types of interactions restrict the multi-gene expression patterns that characterize different cell-fates. More recently, such temporal single-cell models have been extended to recover the temporal and spatial components of morphogenesis. Since the complete genomic GRN is still unknown and intractable for any organism, and some clear developmental modules have been identified, we focus here on the analysis of well-curated and experimentally grounded small GRN modules. One of the first experimentally grounded GRN that was proposed and validated corresponds to the regulatory module involved in floral organ determination. In this chapter we use this GRN as an example of the methodologies involved in: (1) formalizing and integrating molecular genetic data into the logical functions (Boolean functions) that rule gene interactions and dynamics in a Boolean GRN; (2) the algorithms and computational approaches used to recover the steady-states that correspond to each cell type, as well as the set of initial GRN configurations that lead to each one of such states (i.e., basins of attraction); (3) the approaches used to validate a GRN model using wild type and mutant or overexpression data, or to test the robustness of the GRN being proposed; (4) some of the methods that have been used to incorporate random fluctuations in the GRN Boolean functions and enable stochastic GRN models to address the temporal sequence with which gene configurations and cell fates are attained; (5) the methodologies used to approximate discrete Boolean GRN to continuous systems and their use in further dynamic analyses. The methodologies explained for the GRN of floral organ determination developed here in detail can be applied to any other functional developmental module.

HIV-Envelope–Dependent Cell-Cell Fusion: Quantitative Studies

Interaction in vitro between cells infected with human immunodeficiency virus (HIV) and surrounding, uninfected, target cells often leads to cell fusion and the formation of multinucleated cells, called syncytia. The presence in HIV-infected individuals of virus strains able to induce syncytia in cultures of T cells is associated with disease progression and AIDS. Even in the asymptomatic stage of infection, multinucleated cells have been observed in different organs, indicating that fused cells may be generated and remain viable in the tissues of patients. We used lymphocytic cells transfected for the expression of the HIV-envelope (Env) glycoproteins to develop a method for the direct quantification of fusion events by flow cytometry (Huerta et al., 2006, J. Virol. Methods 138, 17–23; López-Balderas et al., 2007, Virus Res. 123, 138–146). The method involves the staining of fusion partners with lipophilic probes and the use of fluorescence resonance energy transfer (FRET) to distinguish between fused and aggregated cells. We have shown that such a flow-cytometry assay is appropriate for the screening of compounds that have the potential to modulate HIV-Env–mediated cell fusion. Even those syncytia that are small or few in numbers can be detected. Quantitative analysis of the fusion products was performed with this technique; the results indicated that the time of reaction and initial proportion of fusion partners determine the number, relative size, and average cellular composition of syncytia. Heterogeneity of syncytia generated by HIV-Env–mediated cell-cell fusion may result in a variety of possible outcomes that, in turn, may influence the biological properties of the syncytia and surrounding cells, as well as replication of virus. Given the myriad immune abnormalities leading to AIDS, the full understanding of the extent, diverse composition, and role of fused cells in the pathogenesis of, and immune response to, HIV infection is an important, pending issue.

Casimir Forces in Nanostructures

We present a theoretical calculation of Casimir forces in structures made of parallel slabs that can be made of dispersive and absorptive materials. The materials are characterized by the reflectivity amplitude coefficients of the vacuum modes between the slabs. In particular, we present results for an antisymmetric configuration in which one plate is a metal and the other one a dielectric material. As a reference, we also calculate the force for the symmetric case (two metallic or two dielectric slabs). Our results show that the Casimir force could have a relevant contribution to the interaction between the tip and sample in atomic force microscopy experiments.

High-multipolar effects on the Casimir force: The non-retarded limit

We show that the dispersive force between a spherical nanoparticle (with a radius

Reshaping the epigenetic landscape during early flower development: induction of attractor transitions by relative differences in gene decay rates

\le

Spatial dispersion in Casimir forces: a brief review

100 nm) and a substrate is enhanced by several orders of magnitude when the sphere is near to the substrate. We calculate exactly the dispersive force in the non-retarded limit by incorporating the contributions to the interaction from of all the multipolar electromagnetic modes. We show that as the sphere approaches the substrate, the fluctuations of the electromagnetic field, induced by the vacuum and the presence of the substrate, the dispersive force is enhanced by orders of magnitude. We discuss this effect as a function of the size of the sphere.We implement a spectral representation formalism to calculate exactly the non-retarded Casimir force or dispersive van der Waals force, between a spherical nanoparticle and a planar substrate beyond the dipolar approximation. For a sphere of radius R separated a distance z from the substrate, we find that high-order multipole interactions induce extra factors R/z in the force, as compared to the dipole power law of ~ R3/z3. As a consequence, at small separations the non-retarded Casimir force increases by several orders of magnitude with respect to that estimated in the usual dipole approximation.

Spectral representation of the nonretarded dispersive force between a sphere and a substrate

BackgroundGene regulatory network (GRN) dynamical models are standard systems biology tools for the mechanistic understanding of developmental processes and are enabling the formalization of the epigenetic landscape (EL) model.MethodsIn this work we propose a modeling framework which integrates standard mathematical analyses to extend the simple GRN Boolean model in order to address questions regarding the impact of gene specific perturbations in cell-fate decisions during development.ResultsWe systematically tested the propensity of individual genes to produce qualitative changes to the EL induced by modification of gene characteristic decay rates reflecting the temporal dynamics of differentiation stimuli. By applying this approach to the flower specification GRN (FOS-GRN) we uncovered differences in the functional (dynamical) role of their genes. The observed dynamical behavior correlates with biological observables. We found a relationship between the propensity of undergoing attractor transitions between attraction basins in the EL and the direction of differentiation during early flower development - being less likely to induce up-stream attractor transitions as the course of development progresses. Our model also uncovered a potential mechanism at play during the transition from EL basins defining inflorescence meristem to those associated to flower organs meristem. Additionally, our analysis provided a mechanistic interpretation of the homeotic property of the ABC genes, being more likely to produce both an induced inter-attractor transition and to specify a novel attractor. Finally, we found that there is a close relationship between a gene’s topological features and its propensity to produce attractor transitions.ConclusionsThe study of how the state-space associated with a dynamical model of a GRN can be restructured by modulation of genes’ characteristic expression times is an important aid for understanding underlying mechanisms occurring during development. Our contribution offers a simple framework to approach such problem, as exemplified here by the case of flower development. Different GRN models and the effect of diverse inductive signals can be explored within the same framework. We speculate that the dynamical role of specific genes within a GRN, as uncovered here, might give information about which genes are more likely to link a module to other regulatory circuits and signaling transduction pathways.