Implementation of a roadmap for the comprehensive diagnosis, follow-up, and research of childhood leukemias in vulnerable regions of Mexico: results from the PRONAII Strategy

The main objective of the National Project for Research and Incidence of Childhood Leukemias is to reduce early mortality rates for these neoplasms in the vulnerable regions of Mexico. This project was conducted in the states of Oaxaca, Puebla, and Tlaxcala. A key strategy of the project is the implementation of an effective roadmap to ensure that leukemia patients are the target of maximum benefit of interdisciplinary collaboration between researchers, clinicians, surveyors, and laboratories. This strategy guarantees the comprehensive management of diagnosis and follow-up samples of pediatric patients with leukemia, centralizing, managing, and analyzing the information collected. Additionally, it allows for a precise diagnosis and monitoring of the disease through immunophenotype and measurable residual disease (MRD) studies, enhancing research and supporting informed clinical decisions for the first time in these regions through a population-based study. This initiative has significantly improved the diagnostic capacity of leukemia in girls, boys, and adolescents in the regions of Oaxaca, Puebla, and Tlaxcala, providing comprehensive, high-quality care with full coverage in the region. Likewise, it has strengthened collaboration between health institutions, researchers, and professionals in the sector, which contributes to reducing the impact of the disease on the community.

The MADS-box genes SOC1 and AGL24 antagonize XAL2 functions in Arabidopsis thaliana root development

MADS-domain transcription factors play pivotal roles in numerous developmental processes in Arabidopsis thaliana. While their involvement in flowering transition and floral development has been extensively examined, their functions in root development remain relatively unexplored. Here, we explored the function and genetic interaction of three MADS-box genes (XAL2, SOC1 and AGL24) in primary root development. By analyzing loss-of-function and overexpression lines, we found that SOC1 and AGL24, both critical components in flowering transition, redundantly act as repressors of primary root growth as the loss of function of either SOC1 or AGL24 partially recovers the primary root growth, meristem cell number, cell production rate, and the length of fully elongated cells of the short-root mutant xal2-2. Furthermore, we observed that the simultaneous overexpression of AGL24 and SOC1 leads to short-root phenotypes, affecting meristem cell number and fully elongated cell size, whereas SOC1 overexpression is sufficient to affect columella stem cell differentiation. Additionally, qPCR analyses revealed that these genes exhibit distinct modes of transcriptional regulation in roots compared to what has been previously reported for aerial tissues. We identified 100 differentially expressed genes in xal2-2 roots by RNA-seq. Moreover, our findings revealed that the expression of certain genes involved in cell differentiation, as well as stress responses, which are either upregulated or downregulated in the xal2-2 mutant, reverted to WT levels in the absence of SOC1 or AGL24.

Hybrid lineages of CD4+ T cells: a handbook update

CD4+ T lymphocytes have been classified into several lineages, according to their gene expression profiles and their effector responses. Interestingly, recent evidence is showing that many lineages could yield hybrid phenotypes with unique properties and functions. It has been reported that such hybrid lineages might underlie pathologies or may function as effector cells with protection capacities against molecular threats. In this work, we reviewed the characteristics of the hybrid lineages reported in the literature, in order to identify the expression profiles that characterize them and the markers that could be used to identify them. We also review the differentiation cues that elicit their hybrid origin and what is known about their physiological roles.

Subclassification of B-acute lymphoblastic leukemia according to age, immunophenotype and microenvironment, predicts MRD risk in Mexican children from vulnerable regions

Introduction

The decisive key to disease-free survival in B-cell precursor acute lymphoblastic leukemia in children, is the combination of diagnostic timeliness and treatment efficacy, guided by accurate patient risk stratification. Implementation of standardized and high-precision diagnostic/prognostic systems is particularly important in the most marginalized geographic areas in Mexico, where high numbers of the pediatric population resides and the highest relapse and early death rates due to acute leukemias are recorded even in those cases diagnosed as standard risk.

Methods

By using a multidimensional and integrated analysis of the immunophenotype of leukemic cells, the immunological context and the tumor microenvironment, this study aim to capture the snapshot of acute leukemia at disease debut of a cohort of Mexican children from vulnerable regions in Puebla, Oaxaca and Tlaxcala and its potential use in risk stratification.

Results and discussion

Our findings highlight the existence of a distinct profile of ProB-ALL in children older than 10 years, which is associated with a six-fold increase in the risk of developing measurable residual disease (MRD). Along with the absence of CD34+ seminal cells for normal hematopoiesis, this ProB-ALL subtype exhibited several characteristics related to poor prognosis, including the high expression level of myeloid lineage markers such as MPO and CD33, as well as upregulation of CD19, CD34, CD24, CD20 and nuTdT. In contrast, it showed a trend towards decreased expression of CD9, CD81, CD123, CD13, CD15 and CD21. Of note, the mesenchymal stromal cell compartment constituting their leukemic niche in the bone marrow, displayed characteristics of potential suppressive microenvironment, such as the expression of Gal9 and IDO1, and the absence of the chemokine CXCL11. Accordingly, adaptive immunity components were poorly represented. Taken together, our results suggest, for the first time, that a biologically distinct subtype of ProB-ALL emerges in vulnerable adolescents, with a high risk of developing MRD. Rigorous research on potential enhancing factors, environmental or lifestyle, is crucial for its detection and prevention. The use of the reported profile for early risk stratification is suggested.

XAANTAL1 Reveals an Additional Level of Flowering Regulation in the Shoot Apical Meristem in Response to Light and Increased Temperature in Arabidopsis

Light and photoperiod are environmental signals that regulate flowering transition. In plants like Arabidopsis thaliana, this regulation relies on CONSTANS, a transcription factor that is negatively posttranslational regulated by phytochrome B during the morning, while it is stabilized by PHYA and cryptochromes 1/2 at the end of daylight hours. CO induces the expression of FT, whose protein travels from the leaves to the apical meristem, where it binds to FD to regulate some flowering genes. Although PHYB delays flowering, we show that light and PHYB positively regulate XAANTAL1 and other flowering genes in the shoot apices. Also, the genetic data indicate that XAL1 and FD participate in the same signaling pathway in flowering promotion when plants are grown under a long-day photoperiod at 22 °C. By contrast, XAL1 functions independently of FD or PIF4 to induce flowering at higher temperatures (27 °C), even under long days. Furthermore, XAL1 directly binds to FD, SOC1, LFY, and AP1 promoters. Our findings lead us to propose that light and temperature influence the floral network at the meristem level in a partially independent way of the signaling generated from the leaves.

Local and Regional Dynamics of Native Maize Seed Lot Use by Small-Scale Producers and Their Impact on Transgene Presence in Three Mexican States

Mexico harbors over 50% of maize's genetic diversity in the Americas. Native maize varieties are actively managed by small-scale producers within a diverse array of cultivation systems. Seed lot use, exchange and admixture has consequences for the in situ conservation of such varieties. Here we analyze native maize seed management dynamics from 906 small-scale producers surveyed in three Mexican states: Mexico City, Oaxaca and Chiapas. Furthermore, we analyze how their management practices can relate to transgene presence, which was experimentally documented for maize samples associated with the applied surveys. Through a data mining approach, we investigated which practices might be related with a higher probability of transgene presence. The variables found to have a strong spatial association with transgene presence were: for Mexico City, maize producers with larger parcels; for Oaxaca, producer's age (43-46 years) and the sale of seed; for Chiapas, the use of agricultural machinery and younger producers (37-43 years). Additionally, transgene presence and frequency within the socioeconomic regions of Oaxaca and Chiapas was analyzed. In Oaxaca, higher transgene frequencies occurred in regions where transgene presence had been previously reported. In Chiapas, the border regions with Guatemala as well as a region where reproduction of improved seed takes place, the highest proportion of positive samples were found. A detailed mapping of regional seed markets and seed exchange sites together with deployment of national and local biosecurity measures, could help prevent the further spread of transgenes into native maize varieties, as well as improve conservation efforts.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

A Green Light to Switch on Genes: Revisiting Trithorax on Plants

The Trithorax Group (TrxG) is a highly conserved multiprotein activation complex, initially defined by its antagonistic activity with the PcG repressor complex. TrxG regulates transcriptional activation by the deposition of H3K4me3 and H3K36me3 marks. According to the function and evolutionary origin, several proteins have been defined as TrxG in plants; nevertheless, little is known about their interactions and if they can form TrxG complexes. Recent evidence suggests the existence of new TrxG components as well as new interactions of some TrxG complexes that may be acting in specific tissues in plants. In this review, we bring together the latest research on the topic, exploring the interactions and roles of TrxG proteins at different developmental stages, required for the fine-tuned transcriptional activation of genes at the right time and place. Shedding light on the molecular mechanism by which TrxG is recruited and regulates transcription.

Combined Approach of GWAS and Phylogenetic Analyses to Identify New Candidate Genes That Participate in Arabidopsis thaliana Primary Root Development Using Cellular Measurements and Primary Root Length

Genome-wide association studies (GWAS) have allowed the identification of different loci associated with primary root (PR) growth, and Arabidopsis is an excellent model for these studies. The PR length is controlled by cell proliferation, elongation, and differentiation; however, the specific contribution of proliferation and differentiation in the control of PR growth is still poorly studied. To this end, we analyzed 124 accessions and used a GWAS approach to identify potential causal genomic regions related to four traits: PR length, growth rate, cell proliferation and cell differentiation. Twenty-three genes and five statistically significant SNPs were identified. The SNP with the highest score mapped to the fifth exon of NAC048 and this change makes a missense variant in only 33.3% of the accessions with a large PR, compared with the accessions with a short PR length. Moreover, we detected five more SNPs in this gene and in NAC3 that allow us to discover closely related accessions according to the phylogenetic tree analysis. We also found that the association between genetic variants among the 18 genes with the highest scores in our GWAS and the phenotypic classes into which we divided our accessions are not straightforward and likely follow historical patterns.

Live Plant Cell Tracking: Fiji plugin to analyze cell proliferation dynamics and understand morphogenesis

Arabidopsis (Arabidopsis thaliana) primary and lateral roots (LRs) are well suited for 3D and 4D microscopy, and their development provides an ideal system for studying morphogenesis and cell proliferation dynamics. With fast-advancing microscopy techniques used for live-imaging, whole tissue data are increasingly available, yet present the great challenge of analyzing complex interactions within cell populations. We developed a plugin "Live Plant Cell Tracking" (LiPlaCeT) coupled to the publicly available ImageJ image analysis program and generated a pipeline that allows, with the aid of LiPlaCeT, 4D cell tracking and lineage analysis of populations of dividing and growing cells. The LiPlaCeT plugin contains ad hoc ergonomic curating tools, making it very simple to use for manual cell tracking, especially when the signal-to-noise ratio of images is low or variable in time or 3D space and when automated methods may fail. Performing time-lapse experiments and using cell-tracking data extracted with the assistance of LiPlaCeT, we accomplished deep analyses of cell proliferation and clonal relations in the whole developing LR primordia and constructed genealogical trees. We also used cell-tracking data for endodermis cells of the root apical meristem (RAM) and performed automated analyses of cell population dynamics using ParaView software (also publicly available). Using the RAM as an example, we also showed how LiPlaCeT can be used to generate information at the whole-tissue level regarding cell length, cell position, cell growth rate, cell displacement rate, and proliferation activity. The pipeline will be useful in live-imaging studies of roots and other plant organs to understand complex interactions within proliferating and growing cell populations. The plugin includes a step-by-step user manual and a dataset example that are available at https://www.ibt.unam.mx/documentos/diversos/LiPlaCeT.zip.

Unraveling the role of epigenetic regulation in asymmetric cell division during plant development

Asymmetric cell divisions are essential to generate different cellular lineages. In plants, asymmetric cell divisions regulate the correct formation of the embryo, stomatal cells, apical and root meristems, and lateral roots. Current knowledge of regulation of asymmetric cell divisions suggests that, in addition to the function of key transcription factor networks, epigenetic mechanisms play crucial roles. Therefore, we highlight the importance of epigenetic regulation and chromatin dynamics for integration of signals and specification of cells that undergo asymmetric cell divisions, as well as for cell maintenance and cell fate establishment of both progenitor and daughter cells. We also discuss the polarization and segregation of cell components to ensure correct epigenetic memory or resetting of epigenetic marks during asymmetric cell divisions.

The flowering transition pathways converge into a complex gene regulatory network that underlies the phase changes of the shoot apical meristem in Arabidopsis thaliana

Post-embryonic plant development is characterized by a period of vegetative growth during which a combination of intrinsic and extrinsic signals triggers the transition to the reproductive phase. To understand how different flowering inducing and repressing signals are associated with phase transitions of the Shoot Apical Meristem (SAM), we incorporated available data into a dynamic gene regulatory network model for Arabidopsis thaliana. This Flowering Transition Gene Regulatory Network (FT-GRN) formally constitutes a dynamic system-level mechanism based on more than three decades of experimental data on flowering. We provide novel experimental data on the regulatory interactions of one of its twenty-three components: a MADS-box transcription factor XAANTAL2 (XAL2). These data complement the information regarding flowering transition under short days and provides an example of the type of questions that can be addressed by the FT-GRN. The resulting FT-GRN is highly connected and integrates developmental, hormonal, and environmental signals that affect developmental transitions at the SAM. The FT-GRN is a dynamic multi-stable Boolean system, with 2²³ possible initial states, yet it converges into only 32 attractors. The latter are coherent with the expression profiles of the FT-GRN components that have been experimentally described for the developmental stages of the SAM. Furthermore, the attractors are also highly robust to initial states and to simulated perturbations of the interaction functions. The model recovered the meristem phenotypes of previously described single mutants. We also analyzed the attractors landscape that emerges from the postulated FT-GRN, uncovering which set of signals or components are critical for reproductive competence and the time-order transitions observed in the SAM. Finally, in the context of such GRN, the role of XAL2 under short-day conditions could be understood. Therefore, this model constitutes a robust biological module and the first multi-stable, dynamical systems biology mechanism that integrates the genetic flowering pathways to explain SAM phase transitions.

Structural Analysis of SARS-CoV-2 ORF8 Protein: Pathogenic and Therapeutic Implications

Current therapeutic strategies and vaccines against SARS-CoV-2 are mainly focused on the Spike protein despite there are other viral proteins with important roles in COVID-19 pathogenicity. For example, ORF8 restructures vesicular trafficking in the host cell, impacts intracellular immunity through the IFN-I signaling, and growth pathways through the mitogen-activated protein kinases (MAPKs). In this mini-review, we analyze the main structural similarities of ORF8 with immunological molecules such as IL-1, contributing to the immunological deregulation observed in COVID-19. We also propose that the blockage of some effector functions of ORF8 with Rapamycin, such as the mTORC1 activation through MAPKs 40 pathway, with Rapamycin, can be a promising approach to reduce COVID-19 mortality.

The Epigenetic Faces of ULTRAPETALA1

ULTRAPETALA1 (ULT1) is a versatile plant-exclusive protein, initially described as a trithorax group (TrxG) factor that regulates transcriptional activation and counteracts polycomb group (PcG) repressor function. As part of TrxG, ULT1 interacts with ARABIDOPSIS TRITHORAX1 (ATX1) to regulate H3K4me3 activation mark deposition. However, our recent studies indicate that ULT1 can also act independently of ATX1. Moreover, the ULT1 ability to interact with transcription factors (TFs) and PcG proteins indicates that it is a versatile protein with other roles. Therefore, in this work we revised recent information about the function of Arabidopsis ULT1 to understand the roles of ULT1 in plant development. Furthermore, we discuss the molecular mechanisms of ULT1, highlighting its epigenetic role, in which ULT1 seems to have characteristics of an epigenetic molecular switch that regulates repression and activation processes via TrxG and PcG complexes.

Spatio-Temporal Dynamics of the Patterning of Arabidopsis Flower Meristem

The qualitative model presented in this work recovers the onset of the four fields that correspond to those of each floral organ whorl of Arabidopsis flower, suggesting a mechanism for the generation of the positional information required for the differential expression of the A, B, and C identity genes according to the ABC model for organ determination during early stages of flower development. Our model integrates a previous model for the emergence of WUS pattern in the floral meristem, and shows that this pre-pattern is a necessary but not sufficient condition for the posterior information of the four fields predicted by the ABC model. Furthermore, our model predicts that LFY diffusion along the L1 layer of cells is not a necessary condition for the patterning of the floral meristem.

Hormonal Regulation of Stem Cell Proliferation at the Arabidopsis thaliana Root Stem Cell Niche

The root stem cell niche (SCN) of Arabidopsis thaliana consists of the quiescent center (QC) cells and the surrounding initial stem cells that produce progeny to replenish all the tissues of the root. The QC cells divide rather slowly relative to the initials, yet most root tissues can be formed from these cells, depending on the requirements of the plant. Hormones are fundamental cues that link such needs with the cell proliferation and differentiation dynamics at the root SCN. Nonetheless, the crosstalk between hormone signaling and the mechanisms that regulate developmental adjustments is still not fully understood. Developmental transcriptional regulatory networks modulate hormone biosynthesis, metabolism, and signaling, and conversely, hormonal responses can affect the expression of transcription factors involved in the spatiotemporal patterning at the root SCN. Hence, a complex genetic-hormonal regulatory network underlies root patterning, growth, and plasticity in response to changing environmental conditions. In this review, we summarize the scientific literature regarding the role of hormones in the regulation of QC cell proliferation and discuss how hormonal signaling pathways may be integrated with the gene regulatory network that underlies cell fate in the root SCN. The conceptual framework we present aims to contribute to the understanding of the mechanisms by which hormonal pathways act as integrators of environmental cues to impact on SCN activity.

Integrative Roles of Phytohormones on Cell Proliferation, Elongation and Differentiation in the Arabidopsis thaliana Primary Root

The growth of multicellular organisms relies on cell proliferation, elongation and differentiation that are tightly regulated throughout development by internal and external stimuli. The plasticity of a growth response largely depends on the capacity of the organism to adjust the ratio between cell proliferation and cell differentiation. The primary root of Arabidopsis thaliana offers many advantages toward understanding growth homeostasis as root cells are continuously produced and move from cell proliferation to elongation and differentiation that are processes spatially separated and could be studied along the longitudinal axis. Hormones fine tune plant growth responses and a huge amount of information has been recently generated on the role of these compounds in Arabidopsis primary root development. In this review, we summarized the participation of nine hormones in the regulation of the different zones and domains of the Arabidopsis primary root. In some cases, we found synergism between hormones that function either positively or negatively in proliferation, elongation or differentiation. Intriguingly, there are other cases where the interaction between hormones exhibits unexpected results. Future analysis on the molecular mechanisms underlying crosstalk hormone action in specific zones and domains will unravel their coordination over PR development.

Interplay between Hormones and Several Abiotic Stress Conditions on Arabidopsis thaliana Primary Root Development

As sessile organisms, plants must adjust their growth to withstand several environmental conditions. The root is a crucial organ for plant survival as it is responsible for water and nutrient acquisition from the soil and has high phenotypic plasticity in response to a lack or excess of them. How plants sense and transduce their external conditions to achieve development, is still a matter of investigation and hormones play fundamental roles. Hormones are small molecules essential for plant growth and their function is modulated in response to stress environmental conditions and internal cues to adjust plant development. This review was motivated by the need to explore how Arabidopsis thaliana primary root differentially sense and transduce external conditions to modify its development and how hormone-mediated pathways contribute to achieve it. To accomplish this, we discuss available data of primary root growth phenotype under several hormone loss or gain of function mutants or exogenous application of compounds that affect hormone concentration in several abiotic stress conditions. This review shows how different hormones could promote or inhibit primary root development in A. thaliana depending on their growth in several environmental conditions. Interestingly, the only hormone that always acts as a promoter of primary root development is gibberellins.

Multi-level gene regulatory network models to understand complex mechanisms underlying plant development

Patterning in plant development is the emergent outcome of the feedback-based interplay between tissue-coupled intracellular regulatory networks and physicochemical fields. This interplay gives rise to dynamics that evolve on a wide spectrum of spatiotemporal scales. This imposes important challenges for computational approaches to model the dynamics of plant development. These challenges are being tackled in recent times by computational and mathematical advances that have made progress in the modelling of regulatory networks, as well as in approaches to couple the latter to physicochemical fields. Efforts in this direction are fundamental to identify the dynamical constraints that emerge from non-cellular autonomous activity in cell-fate decisions and patterning, and requires an understanding of how multi-level and multi-scale processes are coupled. Here, we discuss the use of multi-level modeling and simulation tools for the study of multicellular systems, with emphasis on plants. As illustrative examples, we discuss recent works elucidating the mechanisms that underlie patterning in the root meristem of Arabidopsis thaliana, and in plant responses to environmental conditions.

Dynamical modeling predicts an inflammation-inducible CXCR7+ B cell precursor with potential implications in lymphoid blockage pathologies

Background

The blockage at the early B lymphoid cell development pathway within the bone marrow is tightly associated with hematopoietic and immune diseases, where the disruption of basal regulatory networks prevents the continuous replenishment of functional B cells. Dynamic computational models may be instrumental for the comprehensive understanding of mechanisms underlying complex differentiation processes and provide novel prediction/intervention platforms to reinvigorate the system.

Methods

By reconstructing a three-module regulatory network including genetic transcription, intracellular transduction, and microenvironment communication, we have investigated the early B lineage cell fate decisions in normal and pathological settings. The early B cell differentiation network was simulated as a Boolean model and then transformed, using fuzzy logic, to a continuous model. We tested null and overexpression mutants to analyze the emergent behavior of the network. Due to its importance in inflammation, we investigated the effect of NFkB induction at different early B cell differentiation stages.

Results

While the exhaustive synchronous and asynchronous simulation of the early B cell regulatory network (eBCRN) reproduced the configurations of the hematopoietic progenitors and early B lymphoid precursors of the pathway, its simulation as a continuous model with fuzzy logics suggested a transient IL-7R⁺ ProB-to-Pre-B subset expressing pre-BCR and a series of dominant B-cell transcriptional factors. This conspicuous differentiating cell population up-regulated CXCR7 and reduced CXCR4 and FoxO1 expression levels. Strikingly, constant but intermediate NFkB signaling at specific B cell differentiation stages allowed stabilization of an aberrant CXCR7⁺ pre-B like phenotype with apparent affinity to proliferative signals, while under constitutive overactivation of NFkB, such cell phenotype was aberrantly exacerbated from the earliest stage of common lymphoid progenitors. Our mutant models revealed an abnormal delay in the BCR assembly upon NFkB activation, concomitant to sustained Flt3 signaling, down-regulation of Ebf1, Irf4 and Pax5 genes transcription, and reduced Ig recombination, pointing to a potential lineage commitment blockage.

Discussion

For the first time, an inducible CXCR7^hi B cell precursor endowed with the potential capability of shifting central lymphoid niches, is inferred by computational modeling. Its phenotype is compatible with that of leukemia-initiating cells and might be the foundation that bridges inflammation with blockage-related malignancies and a wide range of immunological diseases. Besides the predicted differentiation impairment, inflammation-inducible phenotypes open the possibility of newly formed niches colonized by the reported precursor. Thus, emergent bone marrow ecosystems are predicted following a pro-inflammatory induction, that may lead to hematopoietic instability associated to blockage pathologies.

Beyond What Your Retina Can See: Similarities of Retinoblastoma Function between Plants and Animals, from Developmental Processes to Epigenetic Regulation

The Retinoblastoma protein (pRb) is a key cell cycle regulator conserved in a wide variety of organisms. Experimental analysis of pRb's functions in animals and plants has revealed that this protein participates in cell proliferation and differentiation processes. In addition, pRb in animals and its orthologs in plants (RBR), are part of highly conserved protein complexes which suggest the possibility that analogies exist not only between functions carried out by pRb orthologs themselves, but also in the structure and roles of the protein networks where these proteins are involved. Here, we present examples of pRb/RBR participation in cell cycle control, cell differentiation, and in the regulation of epigenetic changes and chromatin remodeling machinery, highlighting the similarities that exist between the composition of such networks in plants and animals.

The COVID-19 Pandemic and Paradigm Change in Global Scientific Research

The current pandemic has rocked the lives of human beings every-where in ways never imagined, forcing us to question where our civilization is headed. In this article, we explore and discuss scien-tifi c evidence that helps explain recent events in the context of the COVID-19 pandemic.COVID-19 is caused by infection with a zoonotic-origin novel virus, SARS-CoV-2, that is genetically close to two coronavirus types iso-lated in bats. The transmission dynamics to humans from the original and intermediary hosts remain poorly understood, but it is highly likely that the SARS-CoV-2 virus infected humans after undergoing an inter-species transfer from bats to an intermediate species, and from there to human beings. Crossing the species barrier is largely fostered by industrial-scale agricultural practices that simplify original ecosystem connections by reducing biodiversity, facilitating the emergence of new infectious diseases. The scientifi c community has played an exemplary role in responding to this global emergency, working to fi nd timely, relevant solutions for governments and society as a whole. We need to take this opportunity to promote a global and open science that delves into the interrelation-ships of the biological, environmental, social and economic dimen-sions of this and other diseases while questioning current modes of production and their impact on the environment, and thus on human health worldwide. Keywords: Coronavirus infections; communicable diseases; zoonoses; ecosystems; technology, industry, and agriculture; pandemics; global health; Mexico.

A Computational Model of the Endothelial to Mesenchymal Transition

Endothelial cells (ECs) form the lining of lymph and blood vessels. Changes in tissue requirements or wounds may cause ECs to behave as tip or stalk cells. Alternatively, they may differentiate into mesenchymal cells (MCs). These processes are known as EC activation and endothelial-to-mesenchymal transition (EndMT), respectively. EndMT, Tip, and Stalk EC behaviors all require SNAI1, SNAI2, and Matrix metallopeptidase (MMP) function. However, only EndMT inhibits the expression of VE-cadherin, PECAM1, and VEGFR2, and also leads to EC detachment. Physiologically, EndMT is involved in heart valve development, while a defective EndMT regulation is involved in the physiopathology of cardiovascular malformations, congenital heart disease, systemic and organ fibrosis, pulmonary arterial hypertension, and atherosclerosis. Therefore, the control of EndMT has many promising potential applications in regenerative medicine. Despite the fact that many molecular components involved in EC activation and EndMT have been characterized, the system-level molecular mechanisms involved in this process have not been elucidated. Toward this end, hereby we present Boolean network model of the molecular involved in the regulation of EC activation and EndMT. The simulated dynamic behavior of our model reaches fixed and cyclic patterns of activation that correspond to the expected EC and MC cell types and behaviors, recovering most of the specific effects of simple gain and loss-of-function mutations as well as the conditions associated with the progression of several diseases. Therefore, our model constitutes a theoretical framework that can be used to generate hypotheses and guide experimental inquiry to comprehend the regulatory mechanisms behind EndMT. Our main findings include that both the extracellular microevironment and the pattern of molecular activity within the cell regulate EndMT. EndMT requires a lack of VEGFA and sufficient oxygen in the extracellular microenvironment as well as no FLI1 and GATA2 activity within the cell. Additionally Tip cells cannot undergo EndMT directly. Furthermore, the specific conditions that are sufficient to trigger EndMT depend on the specific pattern of molecular activation within the cell.

ULTRAPETALA1 maintains Arabidopsis root stem cell niche independently of ARABIDOPSIS TRITHORAX1

During plant development, morphogenetic processes rely on the activity of meristems. Meristem homeostasis depends on a complex regulatory network constituted by different factors and hormone signaling that regulate gene expression to coordinate the correct balance between cell proliferation and differentiation. ULTRAPETALA1, a transcriptional regulatory protein described as an Arabidopsis Trithorax group factor, has been characterized as a regulator of the shoot and floral meristems activity. Here, we highlight the role of ULTRAPETALA1 in root stem cell niche maintenance. We found that ULTRAPETALA1 is required to regulate both the quiescent center cell division rate and auxin signaling at the root tip. Furthermore, ULTRAPETALA1 regulates columella stem cell differentiation. These roles are independent of the ARABIDOPSIS TRITHORAX1, suggesting a different mechanism by which ULTRAPETALA1 can act in the root apical meristem of Arabidopsis. This work introduces a new component of the regulatory network needed for the root stem cell niche maintenance.

MADS-Box Genes Are Key Components of Genetic Regulatory Networks Involved in Abiotic Stress and Plastic Developmental Responses in Plants

Plants, as sessile organisms, adapt to different stressful conditions, such as drought, salinity, extreme temperatures, and nutrient deficiency, via plastic developmental and growth responses. Depending on the intensity and the developmental phase in which it is imposed, a stress condition may lead to a broad range of responses at the morphological, physiological, biochemical, and molecular levels. Transcription factors are key components of regulatory networks that integrate environmental cues and concert responses at the cellular level, including those that imply a stressful condition. Despite the fact that several studies have started to identify various members of the MADS-box gene family as important molecular components involved in different types of stress responses, we still lack an integrated view of their role in these processes. In this review, we analyze the function and regulation of MADS-box gene family members in response to drought, salt, cold, heat, and oxidative stress conditions in different developmental processes of several plants. In addition, we suggest that MADS-box genes are key components of gene regulatory networks involved in plant responses to stress and plant developmental plasticity in response to seasonal changes in environmental conditions.

Identifying the transition to the maturation zone in three ecotypes of Arabidopsis thaliana roots

The Arabidopsis thaliana (hereafter Arabidopsis) root has become a useful model for studying how organ morphogenesis emerge from the coordination and balance of cell proliferation and differentiation, as both processes may be observed and quantified in the root at different stages of development. Hence, being able to objectively identify and delimit the different stages of root development has been very important. Up to now, three different zones along the longitudinal axis of the primary root of Arabidopsis, have been identified: the root apical meristematic zone (RAM) with two domains [the proliferative (PD) and the transition domain (TD)], the elongation zone (EZ) and the maturation zone (MZ). We previously reported a method to quantify the length of the cells of the meristematic and the elongation zone, as well as the boundaries or transitions between the root domains along the growing part of the Arabidopsis root. In this study, we provide a more accurate criterion to identify the MZ. Traditionally, the transition between the EZ to the MZ has been established by the emergence of the first root-hair bulge in the epidermis, because this emergence coincides with cell maturation in this cell type. But we have found here that after the emergence of the first root-hair bulge some cells continue to elongate and we have confirmed this in three different Arabidopsis ecotypes. We established the limit between the EZ and the MZ by looking for the closest cortical cell with a longer length than the average cell length of 10 cells after the cortical cell closest to the epidermal cell with the first root-hair bulge in these three ecotypes. In Col-0 and Ws this cell is four cells above the one with the root hair bulge and, in the Ler ecotype, this cell is five cells above. To unambiguously identifying the site at which cells stop elongating and attain their final length and fate at the MZ, we propose to calculate the length of completely elongated cortical cells counting 10 cells starting from the sixth cell above the cortical cell closest to the epidermal cell with the first root-hair bulge. We validated this proposal in the three ecotypes analyzed and consider that this proposal may aid at having a more objective way to characterize root phenotypes and compare among them.

A dynamic genetic-hormonal regulatory network model explains multiple cellular behaviors of the root apical meristem of Arabidopsis thaliana

The study of the concerted action of hormones and transcription factors is fundamental to understand cell differentiation and pattern formation during organ development. The root apical meristem of Arabidopsis thaliana is a useful model to address this. It has a stem cell niche near its tip conformed of a quiescent organizer and stem or initial cells around it, then a proliferation domain followed by a transition domain, where cells diminish division rate before transiting to the elongation zone; here, cells grow anisotropically prior to their final differentiation towards the plant base. A minimal model of the gene regulatory network that underlies cell-fate specification and patterning at the root stem cell niche was proposed before. In this study, we update and couple such network with both the auxin and cytokinin hormone signaling pathways to address how they collectively give rise to attractors that correspond to the genetic and hormonal activity profiles that are characteristic of different cell types along A. thaliana root apical meristem. We used a Boolean model of the genetic-hormonal regulatory network to integrate known and predicted regulatory interactions into alternative models. Our analyses show that, after adding some putative missing interactions, the model includes the necessary and sufficient components and regulatory interactions to recover attractors characteristic of the root cell types, including the auxin and cytokinin activity profiles that correlate with different cellular behaviors along the root apical meristem. Furthermore, the model predicts the existence of activity configurations that could correspond to the transition domain. The model also provides a possible explanation for apparently paradoxical cellular behaviors in the root meristem. For example, how auxin may induce and at the same time inhibit WOX5 expression. According to the model proposed here the hormonal regulation of WOX5 might depend on the cell type. Our results illustrate how non-linear multi-stable qualitative network models can aid at understanding how transcriptional regulators and hormonal signaling pathways are dynamically coupled and may underlie both the acquisition of cell fate and the emergence of hormonal activity profiles that arise during complex organ development.

The combination of the functionalities of feedback circuits is determinant for the attractors' number and size in pathway-like Boolean networks

Molecular regulation was initially assumed to follow both a unidirectional and a hierarchical organization forming pathways. Regulatory processes, however, form highly interlinked networks with non-hierarchical and non-unidirectional structures that contain statistically overrepresented circuits or motifs. Here, we analyze the behavior of pathways containing non-unidirectional (i.e. bidirectional) and non-hierarchical interactions that create motifs. In comparison with unidirectional and hierarchical pathways, our pathways have a high diversity of behaviors, characterized by the size and number of attractors. Motifs have been studied individually showing that feedback circuit motifs regulate the number and size of attractors. It is less clear what happens in molecular networks that usually contain multiple feedbacks. Here, we find that the way feedback circuits couple to each other (i.e., the combination of the functionalities of feedback circuits) regulate both the number and size of the attractors. We show that the different expected results of epistasis analysis (a method to infer regulatory interactions) are produced by many non-hierarchical and non-unidirectional structures. Thus, these structures cannot be correctly inferred by epistasis analysis. Finally, we show that the combinations of functionalities, combined with other network properties, allow for a better characterization of regulatory structures.

Boolean Dynamic Modeling Approaches to Study Plant Gene Regulatory Networks: Integration, Validation, and Prediction

Mathematical models based on dynamical systems theory are well-suited tools for the integration of available molecular experimental data into coherent frameworks in order to propose hypotheses about the cooperative regulatory mechanisms driving developmental processes. Computational analysis of the proposed models using well-established methods enables testing the hypotheses by contrasting predictions with observations. Within such framework, Boolean gene regulatory network dynamical models have been extensively used in modeling plant development. Boolean models are simple and intuitively appealing, ideal tools for collaborative efforts between theorists and experimentalists. In this chapter we present protocols used in our group for the study of diverse plant developmental processes. We focus on conceptual clarity and practical implementation, providing directions to the corresponding technical literature.

Longitudinal zonation pattern in Arabidopsis root tip defined by a multiple structural change algorithm

Background and Aims The Arabidopsis thaliana root is a key experimental system in developmental biology. Despite its importance, we are still lacking an objective and broadly applicable approach for identification of number and position of developmental domains or zones along the longitudinal axis of the root apex or boundaries between them, which is essential for understanding the mechanisms underlying cell proliferation, elongation and differentiation dynamics during root development. Methods We used a statistics approach, the multiple structural change algorithm (MSC), for estimating the number and position of developmental transitions in the growing portion of the root apex. Once the positions of the transitions between domains and zones were determined, linear models were used to estimate the critical size of dividing cells (LcritD) and other parameters. Key Results The MSC approach enabled identification of three discrete regions in the growing parts of the root that correspond to the proliferation domain (PD), the transition domain (TD) and the elongation zone (EZ). Simultaneous application of the MSC approach and G2-to-M transition (CycB1;1DB:GFP) and endoreduplication (pCCS52A1:GUS) molecular markers confirmed the presence and position of the TD. We also found that the MADS-box gene XAANTAL1 (XAL1) is required for the wild-type (wt) PD increase in length during the first 2 weeks of growth. Contrary to wt, in the xal1 loss-of-function mutant the increase and acceleration of root growth were not detected. We also found alterations in LcritD in xal1 compared with wt, which was associated with longer cell cycle duration in the mutant. Conclusions The MSC approach is a useful, objective and versatile tool for identification of the PD, TD and EZ and boundaries between them in the root apices and can be used for the phenotyping of different genetic backgrounds, experimental treatments or developmental changes within a genotype. The tool is publicly available at www.ibiologia.com.mx/MSC_analysis.

Root Architecture Diversity and Meristem Dynamics in Different Populations of Arabidopsis thaliana

Arabidopsis thaliana has been an excellent model system for molecular genetic approaches to development and physiology. More recently, the potential of studying various accessions collected from diverse habitats has been started to exploit. Col-0 has been the best-studied accession but we now know that several traits show significant divergences among them. In this work, we focused in the root that has become a key system for development. We studied root architecture and growth dynamics of 12 Arabidopsis accessions. Our data reveal a wide variability in root architecture and root length among accessions. We also found variability in the root apical meristem (RAM), explained mainly by cell size at the RAM transition domain and possibly by peculiar forms of organization at the stem cell niche in some accessions. Contrary to Col-0 reports, in some accessions the RAM size not always explains the variations in the root length; indicating that elongated cell size could be more relevant in the determination of root length than the RAM size itself. This study contributes to investigations dealing with understanding the molecular and cellular basis of phenotypic variation, the role of plasticity on adaptation, and the developmental mechanisms that may restrict phenotypic variation in response to contrasting environmental conditions.

The impact of Polycomb group (PcG) and Trithorax group (TrxG) epigenetic factors in plant plasticity

Current advances indicate that epigenetic mechanisms play important roles in the regulatory networks involved in plant developmental responses to environmental conditions. Hence, understanding the role of such components becomes crucial to understanding the mechanisms underlying the plasticity and variability of plant traits, and thus the ecology and evolution of plant development. We now know that important components of phenotypic variation may result from heritable and reversible epigenetic mechanisms without genetic alterations. The epigenetic factors Polycomb group (PcG) and Trithorax group (TrxG) are involved in developmental processes that respond to environmental signals, playing important roles in plant plasticity. In this review, we discuss current knowledge of TrxG and PcG functions in different developmental processes in response to internal and environmental cues and we also integrate the emerging evidence concerning their function in plant plasticity. Many such plastic responses rely on meristematic cell behavior, including stem cell niche maintenance, cellular reprogramming, flowering and dormancy as well as stress memory. This information will help to determine how to integrate the role of epigenetic regulation into models of gene regulatory networks, which have mostly included transcriptional interactions underlying various aspects of plant development and its plastic response to environmental conditions.

A Dynamic Gene Regulatory Network Model That Recovers the Cyclic Behavior of Arabidopsis thaliana Cell Cycle

Cell cycle control is fundamental in eukaryotic development. Several modeling efforts have been used to integrate the complex network of interacting molecular components involved in cell cycle dynamics. In this paper, we aimed at recovering the regulatory logic upstream of previously known components of cell cycle control, with the aim of understanding the mechanisms underlying the emergence of the cyclic behavior of such components. We focus on Arabidopsis thaliana, but given that many components of cell cycle regulation are conserved among eukaryotes, when experimental data for this system was not available, we considered experimental results from yeast and animal systems. We are proposing a Boolean gene regulatory network (GRN) that converges into only one robust limit cycle attractor that closely resembles the cyclic behavior of the key cell-cycle molecular components and other regulators considered here. We validate the model by comparing our in silico configurations with data from loss- and gain-of-function mutants, where the endocyclic behavior also was recovered. Additionally, we approximate a continuous model and recovered the temporal periodic expression profiles of the cell-cycle molecular components involved, thus suggesting that the single limit cycle attractor recovered with the Boolean model is not an artifact of its discrete and synchronous nature, but rather an emergent consequence of the inherent characteristics of the regulatory logic proposed here. This dynamical model, hence provides a novel theoretical framework to address cell cycle regulation in plants, and it can also be used to propose novel predictions regarding cell cycle regulation in other eukaryotes.

In multicellular organisms, cells undergo a cyclic behavior of DNA duplication and delivery of a copy to daughter cells during cell division. In each of the main cell-cycle (CC) stages different sets of proteins are active and genes are expressed. Understanding how such cycling cellular behavior emerges and is robustly maintained in the face of changing developmental and environmental conditions, remains a fundamental challenge of biology. The molecular components that cycle through DNA duplication and citokinesis are interconnected in a complex regulatory network. Several models of such network have been proposed, although the regulatory network that robustly recovers a limit-cycle steady state that resembles the behavior of CC molecular components has been recovered only in a few cases, and no comprehensive model exists for plants. In this paper we used the plant Arabidopsis thaliana, as a study system to propose a core regulatory network to recover a cyclic attractor that mimics the oscillatory behavior of the key CC components. Our analyses show that the proposed GRN model is robust to transient alterations, and is validated with the loss- and gain-of-function mutants of the CC components. The interactions proposed for Arabidopsis thaliana CC can inspire predictions for further uncovering regulatory motifs in the CC of other organisms including human.

XAANTAL2 (AGL14) Is an Important Component of the Complex Gene Regulatory Network that Underlies Arabidopsis Shoot Apical Meristem Transitions

In Arabidopsis thaliana, multiple genes involved in shoot apical meristem (SAM) transitions have been characterized, but the mechanisms required for the dynamic attainment of vegetative, inflorescence, and floral meristem (VM, IM, FM) cell fates during SAM transitions are not well understood. Here we show that a MADS-box gene, XAANTAL2 (XAL2/AGL14), is necessary and sufficient to induce flowering, and its regulation is important in FM maintenance and determinacy. xal2 mutants are late flowering, particularly under short-day (SD) condition, while XAL2 overexpressing plants are early flowering, but their flowers have vegetative traits. Interestingly, inflorescences of the latter plants have higher expression levels of LFY, AP1, and TFL1 than wild-type plants. In addition we found that XAL2 is able to bind the TFL1 regulatory regions. On the other hand, the basipetal carpels of the 35S::XAL2 lines lose determinacy and maintain high levels of WUS expression under SD condition. To provide a mechanistic explanation for the complex roles of XAL2 in SAM transitions and the apparently paradoxical phenotypes of XAL2 and other MADS-box (SOC1, AGL24) overexpressors, we conducted dynamic gene regulatory network (GRN) and epigenetic landscape modeling. We uncovered a GRN module that underlies VM, IM, and FM gene configurations and transition patterns in wild-type plants as well as loss and gain of function lines characterized here and previously. Our approach thus provides a novel mechanistic framework for understanding the complex basis of SAM development.

Role of transcriptional regulation in the evolution of plant phenotype: A dynamic systems approach

A growing body of evidence suggests that alterations in transcriptional regulation of genes involved in modulating development are an important part of phenotypic evolution, and this can be documented among species and within populations. While the effects of differential transcriptional regulation in organismal development have been preferentially studied in animal systems, this phenomenon has also been addressed in plants. In this review, we summarize evidence for cis-regulatory mutations, trans-regulatory changes and epigenetic modifications as molecular events underlying important phenotypic alterations, and thus shaping the evolution of plant development. We postulate that a mechanistic understanding of why such molecular alterations have a key role in development, morphology and evolution will have to rely on dynamic models of complex regulatory networks that consider the concerted action of genetic and nongenetic components, and that also incorporate the restrictions underlying the genotype to phenotype mapping process. Developmental Dynamics 244:1074-1095, 2015. © 2015 Wiley Periodicals, Inc.

A model of the regulatory network involved in the control of the cell cycle and cell differentiation in the Caenorhabditis elegans vulva

BACKGROUND

There are recent experimental reports on the cross-regulation between molecules involved in the control of the cell cycle and the differentiation of the vulval precursor cells (VPCs) of Caenorhabditis elegans. Such discoveries provide novel clues on how the molecular mechanisms involved in the cell cycle and cell differentiation processes are coordinated during vulval development. Dynamic computational models are helpful to understand the integrated regulatory mechanisms affecting these cellular processes.

RESULTS

Here we propose a simplified model of the regulatory network that includes sufficient molecules involved in the control of both the cell cycle and cell differentiation in the C. elegans vulva to recover their dynamic behavior. We first infer both the topology and the update rules of the cell cycle module from an expected time series. Next, we use a symbolic algorithmic approach to find which interactions must be included in the regulatory network. Finally, we use a continuous-time version of the update rules for the cell cycle module to validate the cyclic behavior of the network, as well as to rule out the presence of potential artifacts due to the synchronous updating of the discrete model. We analyze the dynamical behavior of the model for the wild type and several mutants, finding that most of the results are consistent with published experimental results.

CONCLUSIONS

Our model shows that the regulation of Notch signaling by the cell cycle preserves the potential of the VPCs and the three vulval fates to differentiate and de-differentiate, allowing them to remain completely responsive to the concentration of LIN-3 and lateral signal in the extracellular microenvironment.

The MADS transcription factor XAL2/AGL14 modulates auxin transport during Arabidopsis root development by regulating PIN expression

Elucidating molecular links between cell-fate regulatory networks and dynamic patterning modules is a key for understanding development. Auxin is important for plant patterning, particularly in roots, where it establishes positional information for cell-fate decisions. PIN genes encode plasma membrane proteins that serve as auxin efflux transporters; mutations in members of this gene family exhibit smaller roots with altered root meristems and stem-cell patterning. Direct regulators of PIN transcription have remained elusive. Here, we establish that a MADS-box gene (XAANTAL2, XAL2/AGL14) controls auxin transport via PIN transcriptional regulation during Arabidopsis root development; mutations in this gene exhibit altered stem-cell patterning, root meristem size, and root growth. XAL2 is necessary for normal shootward and rootward auxin transport, as well as for maintaining normal auxin distribution within the root. Furthermore, this MADS-domain transcription factor upregulates PIN1 and PIN4 by direct binding to regulatory regions and it is required for PIN4-dependent auxin response. In turn, XAL2 expression is regulated by auxin levels thus establishing a positive feedback loop between auxin levels and PIN regulation that is likely to be important for robust root patterning.

Cell patterns emerge from coupled chemical and physical fields with cell proliferation dynamics: the Arabidopsis thaliana root as a study system

A central issue in developmental biology is to uncover the mechanisms by which stem cells maintain their capacity to regenerate, yet at the same time produce daughter cells that differentiate and attain their ultimate fate as a functional part of a tissue or an organ. In this paper we propose that, during development, cells within growing organs obtain positional information from a macroscopic physical field that is produced in space while cells are proliferating. This dynamical interaction triggers and responds to chemical and genetic processes that are specific to each biological system. We chose the root apical meristem of Arabidopsis thaliana to develop our dynamical model because this system is well studied at the molecular, genetic and cellular levels and has the key traits of multicellular stem-cell niches. We built a dynamical model that couples fundamental molecular mechanisms of the cell cycle to a tension physical field and to auxin dynamics, both of which are known to play a role in root development. We perform extensive numerical calculations that allow for quantitative comparison with experimental measurements that consider the cellular patterns at the root tip. Our model recovers, as an emergent pattern, the transition from proliferative to transition and elongation domains, characteristic of stem-cell niches in multicellular organisms. In addition, we successfully predict altered cellular patterns that are expected under various applied auxin treatments or modified physical growth conditions. Our modeling platform may be extended to explicitly consider gene regulatory networks or to treat other developmental systems.

The emergence of tumors results from altered cell differentiation and proliferation during organ and tissue development. Understanding how such altered or normal patterns are established is still a challenge. Molecular genetic approaches to understanding pattern formation have searched for key central genetic controllers. However, biological patterns emerge as a consequence of coupled complex genetic and non-genetic sub-systems operating at various spatial and temporal scales and levels of organization. We present a two-dimensional model and simulation benchmark that considers the integrated dynamics of physical and chemical fields that result from cell proliferation. We aim at understanding how the cellular patterns of stem-cell niches emerge. In these, organizer cells with very low rates of proliferation are surrounded by stem cells with slightly higher proliferation rates that transit to a domain of active proliferation and then of elongation and differentiation. We quantified such cellular patterns in the Arabidopsis thaliana root to test our theoretical propositions. The results of our simulations closely mimic observed root cellular patterns, thus providing a proof of principle that coupled physical fields and chemical processes under active cell proliferation give rise to stem-cell patterns. Our framework may be extended to other developmental systems and to consider gene regulatory networks.

When ABC becomes ACB

Understanding how the information contained in genes is mapped onto the phenotypes, and deriving formal frameworks to search for generic aspects of developmental constraints and evolution remains one of the main challenges of contemporary biological research. The Mexican endemic triurid Lacandonia schismatica (Lacandoniaceae), a mycoheterotrophic monocotyledonous plant with hermaphroditic reproductive axes is alone among 250,000 species of angiosperms, as it has central stamens surrounded by a peripheral gynoecium, representing a natural instance of a homeotic mutant. Based on the classical ABC model of flower development, it has recently been shown that the B-function gene APETALA3 (AP3), essential for stamen identity, was displaced toward the flower centre in L. schismatica (ABC to ACB) from the early stages of flower development. A functional conservation of B-function genes from L. schismatica through the rescue of B-gene mutants in Arabidopsis thaliana, as well as conserved protein interactions, has also been demonstrated. Thus, it has been shown that relatively simple genetic alterations may underlie large morphological shifts fixed in extant natural populations. Nevertheless, critical questions remain in order to have a full and sufficient explanation of the molecular genetic mechanisms underlying L. schismatica's unique floral arrangement. Evolutionary approaches to developmental mechanisms and systems biology, including high-throughput functional genomic studies and models of complex developmental gene regulatory networks, constitute two main approaches to meet such a challenge. In this review, the aim is to address some of the pending questions with the ultimate goal of investigating further the mechanisms of L. schismatica's unique homeotic flower arrangement and its evolution.

B-function expression in the flower center underlies the homeotic phenotype of Lacandonia schismatica (Triuridaceae)

Spontaneous homeotic transformations have been described in natural populations of both plants and animals, but little is known about the molecular-genetic mechanisms underlying these processes in plants. In the ABC model of floral organ identity in Arabidopsis thaliana, the B- and C-functions are necessary for stamen morphogenesis, and C alone is required for carpel identity. We provide ABC model-based molecular-genetic evidence that explains the unique inside-out homeotic floral organ arrangement of the monocotyledonous mycoheterotroph species Lacandonia schismatica (Triuridaceae) from Mexico. Whereas a quarter million flowering plant species bear central carpels surrounded by stamens, L. schismatica stamens occur in the center of the flower and are surrounded by carpels. The simplest explanation for this is that the B-function is displaced toward the flower center. Our analyses of the spatio-temporal pattern of B- and C-function gene expression are consistent with this hypothesis. The hypothesis is further supported by conservation between the B-function genes of L. schismatica and Arabidopsis, as the former are able to rescue stamens in Arabidopsis transgenic complementation lines, and Ls-AP3 and Ls-PI are able to interact with each other and with the corresponding Arabidopsis B-function proteins in yeast. Thus, relatively simple molecular modifications may underlie important morphological shifts in natural populations of extant plant taxa.

B-function expression in the flower center underlies the homeotic phenotype of Lacandonia schismatica (Triuridaceae)

Spontaneous homeotic transformations have been described in natural populations of both plants and animals, but little is known about the molecular-genetic mechanisms underlying these processes in plants. In the ABC model of floral organ identity in Arabidopsis thaliana, the B- and C-functions are necessary for stamen morphogenesis, and C alone is required for carpel identity. We provide ABC model-based molecular-genetic evidence that explains the unique inside-out homeotic floral organ arrangement of the monocotyledonous mycoheterotroph species Lacandonia schismatica (Triuridaceae) from Mexico. Whereas a quarter million flowering plant species bear central carpels surrounded by stamens, L. schismatica stamens occur in the center of the flower and are surrounded by carpels. The simplest explanation for this is that the B-function is displaced toward the flower center. Our analyses of the spatio-temporal pattern of B- and C-function gene expression are consistent with this hypothesis. The hypothesis is further supported by conservation between the B-function genes of L. schismatica and Arabidopsis, as the former are able to rescue stamens in Arabidopsis transgenic complementation lines, and Ls-AP3 and Ls-PI are able to interact with each other and with the corresponding Arabidopsis B-function proteins in yeast. Thus, relatively simple molecular modifications may underlie important morphological shifts in natural populations of extant plant taxa.

Floral morphogenesis: stochastic explorations of a gene network epigenetic landscape

In contrast to the classical view of development as a preprogrammed and deterministic process, recent studies have demonstrated that stochastic perturbations of highly non-linear systems may underlie the emergence and stability of biological patterns. Herein, we address the question of whether noise contributes to the generation of the stereotypical temporal pattern in gene expression during flower development. We modeled the regulatory network of organ identity genes in the Arabidopsis thaliana flower as a stochastic system. This network has previously been shown to converge to ten fixed-point attractors, each with gene expression arrays that characterize inflorescence cells and primordial cells of sepals, petals, stamens, and carpels. The network used is binary, and the logical rules that govern its dynamics are grounded in experimental evidence. We introduced different levels of uncertainty in the updating rules of the network. Interestingly, for a level of noise of around 0.5–10%, the system exhibited a sequence of transitions among attractors that mimics the sequence of gene activation configurations observed in real flowers. We also implemented the gene regulatory network as a continuous system using the Glass model of differential equations, that can be considered as a first approximation of kinetic-reaction equations, but which are not necessarily equivalent to the Boolean model. Interestingly, the Glass dynamics recover a temporal sequence of attractors, that is qualitatively similar, although not identical, to that obtained using the Boolean model. Thus, time ordering in the emergence of cell-fate patterns is not an artifact of synchronous updating in the Boolean model. Therefore, our model provides a novel explanation for the emergence and robustness of the ubiquitous temporal pattern of floral organ specification. It also constitutes a new approach to understanding morphogenesis, providing predictions on the population dynamics of cells with different genetic configurations during development.