A gene regulatory network model that recovers the abaxial-adaxial polarity in Arabidopsis thaliana leaf primordium

Megaphylls, present in the majority of vascular plants, show in many plant lineages an abaxial-adaxial polarity in their dorsoventral axis. This polarity commonly translates into different tissues developing on each side of the leaf blade. This is important because it promotes better photosynthetic efficiency as related to light absorption and gas exchange. Many researchers have studied the molecular bases of the emergence of leaf abaxial-adaxial polarity, showing that it is produced by the interaction and differential expression of particular genes and other molecules. However, until now, it is still unclear if the molecular components documented thus far are sufficient to explain the emergence of leaf polarity. In this work, we integrated the available experimental data to construct a graph of the Gene Regulatory Network (GRN) involved in the formation of abaxial-adaxial polarity in the leaf primordium of Arabidopsis thaliana. This graph consisted of 21 nodes and 47 regulations. We extracted the main components of the graph to obtain a Minimum Network consisting of six genes and 22 possible regulations. Then, we used the Boolean network (BN) formalism to describe the dynamics of this Minimum Network. We identified 1905 distinct BNs that comprised the regulations of the Minimum Network and exclusively generated the two attractors representing the abaxial and adaxial cell types. This highlights the fact that most graphs, including our network, can describe experimentally observed behaviors with many BN dynamics. By performing mutant simulations and robustness analysis, we found that two of the 1905 BNs better reproduce experimentally available information. To produce the expected attractors, both BNs predict the same missing regulations, which we propose should be experimentally analyzed to confirm their existence. Interestingly, these two BNs have low robustness to perturbations compared with previously analyzed GRNs. This was an unexpected result since abaxial-adaxial polarity is a robust biological trait, which suggests more components or regulations of the network are missing.