QT-GWAS: A novel method for unveiling biosynthetic loci affecting qualitative metabolic traits.

Although the plant kingdom provides an enormous diversity of metabolites with potentially beneficial applications for humankind, a large fraction of these metabolites and their biosynthetic pathways remain unknown. Resolving metabolite structures and their biosynthetic pathways is key to gaining biological understanding and to allow metabolic engineering. In order to retrieve novel biosynthetic genes involved in specialized metabolism, we developed a novel untargeted method designated as qualitative trait GWAS (QT-GWAS) that subjects qualitative metabolic traits to a genome-wide association study, while the conventional metabolite GWAS (mGWAS) mainly considers the quantitative variation of metabolites. As a proof of the validity of QT-GWAS, 23 and 15 of the retrieved associations identified in Arabidopsis thaliana by QT-GWAS and mGWAS, respectively, were supported by previous research. Furthermore, seven gene-metabolite associations retrieved by QT-GWAS were confirmed in this study through reverse genetics combined with metabolomics and/or in vitro enzyme assays. As such, we established that CYTOCHROME P450 706A5 (CYP706A5) is involved in the biosynthesis of chroman derivatives, UDP-GLYCOSYLTRANSFERASE 76C3 (UGT76C3) is able to hexosylate guanine in vitro and in planta, and SULFOTRANSFERASE 202B1 (SULT202B1) catalyzes the sulfation of neolignans in vitro. Collectively, our study demonstrates that the untargeted QT-GWAS method can retrieve valid gene-metabolite associations at the level of enzyme-encoding genes, even new associations that cannot be found by the conventional mGWAS, providing a new approach for dissecting qualitative metabolic traits.

Opposing effects of trans- and cis-cinnamic acid during rice coleoptile elongation.

The phenylpropanoid cinnamic acid (CA) is a plant metabolite that can occur under a trans- or cis-form. In contrast to the proven bioactivity of the cis-form (c-CA), the activity of trans-CA (t-CA) is still a matter of debate. We tested both compounds using a submerged rice coleoptile assay and demonstrated that they have opposite effects on cell elongation. Notably, in the tip of rice coleoptile t-CA showed an inhibiting and c-CA a stimulating activity. By combining transcriptomics and (untargeted) metabolomics with activity assays and genetic and pharmacological experiments, we aimed to explain the underlying mechanistic processes. We propose a model in which c-CA treatment activates proton pumps and stimulates acidification of the apoplast, which in turn leads to the loosening of the cell wall, necessary for elongation. We hypothesize that c-CA also inactivates auxin efflux transporters, which might cause a local auxin accumulation in the tip of the coleoptile. For t-CA, the phenotype can partially be explained by a stimulation of cell wall polysaccharide feruloylation, leading to a more rigid cell wall. Metabolite profiling also demonstrated that salicylic acid (SA) derivatives are increased upon t-CA treatment. As SA is a known antagonist of auxin, the shift in SA homeostasis provides an additional explanation of the observed t-CA-mediated restriction on cell growth.

Nanobody-Dependent Delocalization of Endocytic Machinery in Arabidopsis Root Cells Dampens Their Internalization Capacity.

Plant cells perceive and adapt to an ever-changing environment by modifying their plasma membrane (PM) proteome. Whereas secretion deposits new integral membrane proteins, internalization by endocytosis removes membrane proteins and associated ligands, largely with the aid of adaptor protein (AP) complexes and the scaffolding molecule clathrin. Two AP complexes function in clathrin-mediated endocytosis at the PM in plant cells, the heterotetrameric AP-2 complex and the hetero-octameric TPLATE complex (TPC). Whereas single subunit mutants in AP-2 develop into viable plants, genetic mutation of a single TPC subunit causes fully penetrant male sterility and silencing single subunits leads to seedling lethality. To address TPC function in somatic root cells, while minimizing indirect effects on plant growth, we employed nanobody-dependent delocalization of a functional, GFP-tagged TPC subunit, TML, in its respective homozygous genetic mutant background. In order to decrease the amount of functional TPC at the PM, we targeted our nanobody construct to the mitochondria and fused it to TagBFP2 to visualize it independently of its bait. We furthermore limited the effect of our delocalization to those tissues that are easily accessible for live-cell imaging by expressing it from the PIN2 promoter, which is active in root epidermal and cortex cells. With this approach, we successfully delocalized TML from the PM. Moreover, we also show co-recruitment of TML-GFP and AP2A1-TagRFP to the mitochondria, suggesting that our approach delocalized complexes, rather than individual adaptor complex subunits. In line with the specific expression domain, we only observed minor effects on root growth, yet realized a clear reduction of endocytic flux in epidermal root cells. Nanobody-dependent delocalization in plants, here exemplified using a TPC subunit, has the potential to be widely applicable to achieve specific loss-of-function analysis of otherwise lethal mutants.

Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators.

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

Seedling developmental defects upon blocking CINNAMATE-4-HYDROXYLASE are caused by perturbations in auxin transport.

The phenylpropanoid pathway serves a central role in plant metabolism, providing numerous compounds involved in diverse physiological processes. Most carbon entering the pathway is incorporated into lignin. Although several phenylpropanoid pathway mutants show seedling growth arrest, the role for lignin in seedling growth and development is unexplored. We use complementary pharmacological and genetic approaches to block CINNAMATE-4-HYDROXYLASE (C4H) functionality in Arabidopsis seedlings and a set of molecular and biochemical techniques to investigate the underlying phenotypes. Blocking C4H resulted in reduced lateral rooting and increased adventitious rooting apically in the hypocotyl. These phenotypes coincided with an inhibition in AUX transport. The upstream accumulation in cis-cinnamic acid was found to be likely to cause polar AUX transport inhibition. Conversely, a downstream depletion in lignin perturbed phloem-mediated AUX transport. Restoring lignin deposition effectively reestablished phloem transport and, accordingly, AUX homeostasis. Our results show that the accumulation of bioactive intermediates and depletion in lignin jointly cause the aberrant phenotypes upon blocking C4H, and demonstrate that proper deposition of lignin is essential for the establishment of AUX distribution in seedlings. Our data position the phenylpropanoid pathway and lignin in a new physiological framework, consolidating their importance in plant growth and development.

Genetic Variability of Arabidopsis thaliana Mature Root System Architecture and Genome-Wide Association Study.

Root system architecture (RSA) has a direct influence on the efficiency of nutrient uptake and plant growth, but the genetics of RSA are often studied only at the seedling stage. To get an insight into the genetic blueprint of a more mature RSA, we exploited natural variation and performed a detailed in vitro study of 241 Arabidopsis thaliana accessions using large petri dishes. A comprehensive analysis of 17 RSA traits showed high variability among the different accessions, unveiling correlations between traits and conditions of the natural habitat of the plants. A sub-selection of these accessions was grown in water-limiting conditions in a rhizotron set-up, which revealed that especially the spatial distribution showed a high consistency between in vitro and ex vitro conditions, while in particular, a large root area in the lower zone favored drought tolerance. The collected RSA phenotype data were used to perform genome-wide association studies (GWAS), which stands out from the previous studies by its exhaustive measurements of RSA traits on more mature Arabidopsis accessions used for GWAS. As a result, we found not only several genes involved in the lateral root (LR) development or auxin signaling pathways to be associated with RSA traits but also new candidate genes that are potentially involved in the adaptation to the natural habitats.

Comparative transcriptomics enables the identification of functional orthologous genes involved in early leaf growth.

Leaf growth is a complex trait for which many similarities exist in different plant species, suggesting functional conservation of the underlying pathways. However, a global view of orthologous genes involved in leaf growth showing conserved expression in dicots and monocots is currently missing. Here, we present a genome-wide comparative transcriptome analysis between Arabidopsis and maize, identifying conserved biological processes and gene functions active during leaf growth. Despite the orthology complexity between these distantly related plants, 926 orthologous gene groups including 2829 Arabidopsis and 2974 maize genes with similar expression during leaf growth were found, indicating conservation of the underlying molecular networks. We found 65% of these genes to be involved in one-to-one orthology, whereas only 28.7% of the groups with divergent expression had one-to-one orthology. Within the pool of genes with conserved expression, 19 transcription factor families were identified, demonstrating expression conservation of regulators active during leaf growth. Additionally, 25 Arabidopsis and 25 maize putative targets of the TCP transcription factors with conserved expression were determined based on the presence of enriched transcription factor binding sites. Based on large-scale phenotypic data, we observed that genes with conserved expression have a higher probability to be involved in leaf growth and that leaf-related phenotypes are more frequently present for genes having orthologues between dicots and monocots than clade-specific genes. This study shows the power of integrating transcriptomic with orthology data to identify or select candidates for functional studies during leaf development in flowering plants.

Capturing the phosphorylation and protein interaction landscape of the plant TOR kinase.

The target of rapamycin (TOR) kinase is a conserved regulatory hub that translates environmental and nutritional information into permissive or restrictive growth decisions. Despite the increased appreciation of the essential role of the TOR complex in plants, no large-scale phosphoproteomics or interactomics studies have been performed to map TOR signalling events in plants. To fill this gap, we combined a systematic phosphoproteomics screen with a targeted protein complex analysis in the model plant Arabidopsis thaliana. Integration of the phosphoproteome and protein complex data on the one hand shows that both methods reveal complementary subspaces of the plant TOR signalling network, enabling proteome-wide discovery of both upstream and downstream network components. On the other hand, the overlap between both data sets reveals a set of candidate direct TOR substrates. The integrated network embeds both evolutionarily-conserved and plant-specific TOR signalling components, uncovering an intriguing complex interplay with protein synthesis. Overall, the network provides a rich data set to start addressing fundamental questions about how TOR controls key processes in plants, such as autophagy, auxin signalling, chloroplast development, lipid metabolism, nucleotide biosynthesis, protein translation or senescence.

Functional analysis of Arabidopsis and maize transgenic lines overexpressing the ADP-ribose/NADH pyrophosphohydrolase, AtNUDX7.

The conserved poly(ADP-ribosyl)ation (PAR) pathway consists of three genetic components that are potential targets to modulate the plant's energy homeostasis upon stress with the aim to improve yield stability in crops and help secure food supply. We studied the role of the PAR pathway component ADP-ribose/NADH pyrophosphohydrolase (AtNUDX7) in yield and mild drought stress by using a transgenic approach in Arabidopsis thaliana and maize (Zea mays). Arabidopsis AtNUDX7 cDNA was overexpressed in Arabidopsis and maize by means of the constitutive Cauliflower Mosaic Virus 35S promoter and the strong constitutive Brachypodium distachyon pBdEF1α promoter, respectively. Overexpression of AtNUDX7 in Arabidopsis improved seed parameters that were measured by a novel, automated method, accelerated flowering and reduced inflorescence height. This combination of beneficial traits suggested that AtNUDX7 overexpression in Arabidopsis might enhance the ADP-ribose recycling step and maintain energy levels by supplying an ATP source in the poly(ADP-ribosyl)ation energy homeostasis pathway. Arabidopsis and maize lines with high, medium and low overexpression levels of the AtNUDX7 gene were analysed in automated platforms and the inhibition of several growth parameters was determined under mild drought stress conditions. The data showed that the constitutive overexpression of the Arabidopsis AtNUDX7 gene in Arabidopsis and maize at varying levels did not improve tolerance to mild drought stress, but knocking down AtNUDX7 expression did, however at the expense of general growth under normal conditions.

A user-friendly platform for yeast two-hybrid library screening using next generation sequencing.

Yeast two-hybrid (Y2H) is a well-established genetics-based system that uses yeast to selectively display binary protein-protein interactions (PPIs). To meet the current need to unravel complex PPI networks, several adaptations have been made to establish medium- to high-throughput Y2H screening platforms, with several having successfully incorporated the use of the next-generation sequencing (NGS) technology to increase the scale and sensitivity of the method. However, these have been to date mainly restricted to the use of fully annotated custom-made open reading frame (ORF) libraries and subject to complex downstream data processing. Here, a streamlined Y2H library screening strategy, based on integration of Y2H with NGS, called Y2H-seq, was developed, which allows efficient and reliable screening of Y2H cDNA libraries. To generate proof of concept, the method was applied to screen for interaction partners of two key components of the jasmonate signaling machinery in the model plant Arabidopsis thaliana, resulting in the identification of several previously reported as well as hitherto unknown interactors. Our Y2H-seq method offers a user-friendly, specific and sensitive screening method that allows identification of PPIs without prior knowledge of the organism's ORFs, thereby extending the method to organisms of which the genome has not entirely been annotated yet. The quantitative NGS readout allows to increase genome coverage, thereby overcoming some of the bottlenecks of current Y2H technologies, which will further strengthen the value of the Y2H technology as a discovery platform.

A Spatiotemporal DNA Endoploidy Map of the Arabidopsis Root Reveals Roles for the Endocycle in Root Development and Stress Adaptation.

Somatic polyploidy caused by endoreplication is observed in arthropods, molluscs, and vertebrates but is especially prominent in higher plants, where it has been postulated to be essential for cell growth and fate maintenance. However, a comprehensive understanding of the physiological significance of plant endopolyploidy has remained elusive. Here, we modeled and experimentally verified a high-resolution DNA endoploidy map of the developing Arabidopsis thaliana root, revealing a remarkable spatiotemporal control of DNA endoploidy levels across tissues. Fitting of a simplified model to publicly available data sets profiling root gene expression under various environmental stress conditions suggested that this root endoploidy patterning may be stress-responsive. Furthermore, cellular and transcriptomic analyses revealed that inhibition of endoreplication onset alters the nuclear-to-cellular volume ratio and the expression of cell wall-modifying genes, in correlation with the appearance of cell structural changes. Our data indicate that endopolyploidy might serve to coordinate cell expansion with structural stability and that spatiotemporal endoreplication pattern changes may buffer for stress conditions, which may explain the widespread occurrence of the endocycle in plant species growing in extreme or variable environments.

The 'TranSeq' 3'-end sequencing method for high-throughput transcriptomics and gene space refinement in plant genomes.

High-throughput RNA sequencing has proven invaluable not only to explore gene expression but also for both gene prediction and genome annotation. However, RNA sequencing, carried out on tens or even hundreds of samples, requires easy and cost-effective sample preparation methods using minute RNA amounts. Here, we present TranSeq, a high-throughput 3'-end sequencing procedure that requires 10- to 20-fold fewer sequence reads than the current transcriptomics procedures. TranSeq significantly reduces costs and allows a greater increase in size of sample sets analyzed in a single experiment. Moreover, in comparison with other 3'-end sequencing methods reported to date, we demonstrate here the reliability and immediate applicability of TranSeq and show that it not only provides accurate transcriptome profiles but also produces precise expression measurements of specific gene family members possessing high sequence similarity. This is difficult to achieve in standard RNA-seq methods, in which sequence reads cover the entire transcript. Furthermore, mapping TranSeq reads to the reference tomato genome facilitated the annotation of new transcripts improving >45% of the existing gene models. Hence, we anticipate that using TranSeq will boost large-scale transcriptome assays and increase the spatial and temporal resolution of gene expression data, in both model and non-model plant species. Moreover, as already performed for tomato (ITAG3.0; www.solgenomics.net), we strongly advocate its integration into current and future genome annotations.

Nonselective Chemical Inhibition of Sec7 Domain-Containing ARF GTPase Exchange Factors.

Small GTP-binding proteins from the ADP-ribosylation factor (ARF) family are important regulators of vesicle formation and cellular trafficking in all eukaryotes. ARF activation is accomplished by a protein family of guanine nucleotide exchange factors (GEFs) that contain a conserved catalytic Sec7 domain. Here, we identified and characterized Secdin, a small-molecule inhibitor of Arabidopsis thaliana ARF-GEFs. Secdin application caused aberrant retention of plasma membrane (PM) proteins in late endosomal compartments, enhanced vacuolar degradation, impaired protein recycling, and delayed secretion and endocytosis. Combined treatments with Secdin and the known ARF-GEF inhibitor Brefeldin A (BFA) prevented the BFA-induced PM stabilization of the ARF-GEF GNOM, impaired its translocation from the Golgi to the trans-Golgi network/early endosomes, and led to the formation of hybrid endomembrane compartments reminiscent of those in ARF-GEF-deficient mutants. Drug affinity-responsive target stability assays revealed that Secdin, unlike BFA, targeted all examined Arabidopsis ARF-GEFs, but that the interaction was probably not mediated by the Sec7 domain because Secdin did not interfere with the Sec7 domain-mediated ARF activation. These results show that Secdin and BFA affect their protein targets through distinct mechanisms, in turn showing the usefulness of Secdin in studies in which ARF-GEF-dependent endomembrane transport cannot be manipulated with BFA.

The reduction in maize leaf growth under mild drought affects the transition between cell division and cell expansion and cannot be restored by elevated gibberellic acid levels.

Growth is characterized by the interplay between cell division and cell expansion, two processes that occur separated along the growth zone at the maize leaf. To gain further insight into the transition between cell division and cell expansion, conditions were investigated in which the position of this transition zone was positively or negatively affected. High levels of gibberellic acid (GA) in plants overexpressing the GA biosynthesis gene GA20-OXIDASE (GA20OX-1OE ) shifted the transition zone more distally, whereas mild drought, which is associated with lowered GA biosynthesis, resulted in a more basal positioning. However, the increased levels of GA in the GA20OX-1OE line were insufficient to convey tolerance to the mild drought treatment, indicating that another mechanism in addition to lowered GA levels is restricting growth during drought. Transcriptome analysis with high spatial resolution indicated that mild drought specifically induces a reprogramming of transcriptional regulation in the division zone. 'Leaf Growth Viewer' was developed as an online searchable tool containing the high-resolution data.

From network to phenotype: the dynamic wiring of an Arabidopsis transcriptional network induced by osmotic stress.

Plants have established different mechanisms to cope with environmental fluctuations and accordingly fine-tune their growth and development through the regulation of complex molecular networks. It is largely unknown how the network architectures change and what the key regulators in stress responses and plant growth are. Here, we investigated a complex, highly interconnected network of 20 Arabidopsis transcription factors (TFs) at the basis of leaf growth inhibition upon mild osmotic stress. We tracked the dynamic behavior of the stress-responsive TFs over time, showing the rapid induction following stress treatment, specifically in growing leaves. The connections between the TFs were uncovered using inducible overexpression lines and were validated with transient expression assays. This study resulted in the identification of a core network, composed of ERF6, ERF8, ERF9, ERF59, and ERF98, which is responsible for most transcriptional connections. The analyses highlight the biological function of this core network in environmental adaptation and its redundancy. Finally, a phenotypic analysis of loss-of-function and gain-of-function lines of the transcription factors established multiple connections between the stress-responsive network and leaf growth.

Altered expression of maize PLASTOCHRON1 enhances biomass and seed yield by extending cell division duration.

Maize is the highest yielding cereal crop grown worldwide for grain or silage. Here, we show that modulating the expression of the maize PLASTOCHRON1 (ZmPLA1) gene, encoding a cytochrome P450 (CYP78A1), results in increased organ growth, seedling vigour, stover biomass and seed yield. The engineered trait is robust as it improves yield in an inbred as well as in a panel of hybrids, at several locations and over multiple seasons in the field. Transcriptome studies, hormone measurements and the expression of the auxin responsive DR5rev:mRFPer marker suggest that PLA1 may function through an increase in auxin. Detailed analysis of growth over time demonstrates that PLA1 stimulates the duration of leaf elongation by maintaining dividing cells in a proliferative, undifferentiated state for a longer period of time. The prolonged duration of growth also compensates for growth rate reduction caused by abiotic stresses.

The Transcription Factor MYB29 Is a Regulator of ALTERNATIVE OXIDASE1a.

Plants sense and integrate a variety of signals from the environment through different interacting signal transduction pathways that involve hormones and signaling molecules. Using ALTERNATIVE OXIDASE1a (AOX1a) gene expression as a model system of retrograde or stress signaling between mitochondria and the nucleus, MYB DOMAIN PROTEIN29 (MYB29) was identified as a negative regulator (regulator of alternative oxidase1a 7 [rao7] mutant) in a genetic screen of Arabidopsis (Arabidopsis thaliana). rao7/myb29 mutants have increased levels of AOX1a transcript and protein compared to wild type after induction with antimycin A. A variety of genes previously associated with the mitochondrial stress response also display enhanced transcript abundance, indicating that RAO7/MYB29 negatively regulates mitochondrial stress responses in general. Meta-analysis of hormone-responsive marker genes and identification of downstream transcription factor networks revealed that MYB29 functions in the complex interplay of ethylene, jasmonic acid, salicylic acid, and reactive oxygen species signaling by regulating the expression of various ETHYLENE RESPONSE FACTOR and WRKY transcription factors. Despite an enhanced induction of mitochondrial stress response genes, rao7/myb29 mutants displayed an increased sensitivity to combined moderate light and drought stress. These results uncover interactions between mitochondrial retrograde signaling and the regulation of glucosinolate biosynthesis, both regulated by RAO7/MYB29. This common regulator can explain why perturbation of the mitochondrial function leads to transcriptomic responses overlapping with responses to biotic stress.

Forever Young: The Role of Ubiquitin Receptor DA1 and E3 Ligase BIG BROTHER in Controlling Leaf Growth and Development.

The final size of plant organs is determined by a combination of cell proliferation and cell expansion. Leaves account for a large part of above-ground biomass and provide energy to complete the plant's life cycle. Although the final size of leaves is remarkably constant under fixed environmental conditions, several genes have been described to enhance leaf growth when their expression is modulated. In Arabidopsis (Arabidopsis thaliana), mutations in DA1 and BB increase leaf size, an effect that is synergistically enhanced in the double mutant. Here, we show that overexpression of a dominant-negative version of DA1 enhances leaf size in a broad range of natural accessions of this species, indicating a highly conserved role of this protein in controlling organ size. We also found that during early stages of development, leaves of da1-1 and bb/eod1-2 mutants were already larger than the isogenic Col-0 wild type, but this phenotype was triggered by different cellular mechanisms. Later during development, da1-1 and bb/eod1-2 leaves showed a prolonged longevity, which was enhanced in the double mutant. Conversely, ectopic expression of DA1 or BB restricted growth and promoted leaf senescence. In concert, shortly upon induction of DA1 and BB expression, several marker genes for the transition from proliferation to expansion were highly up-regulated. Additionally, multiple genes involved in maintaining the mitotic cell cycle were rapidly down-regulated and senescence genes were strongly up-regulated, particularly upon BB induction. With these results, we demonstrate that DA1 and BB restrict leaf size and promote senescence through converging and different mechanisms.

Natural Variation of Molecular and Morphological Gibberellin Responses.

Although phytohormones such as gibberellins are essential for many conserved aspects of plant physiology and development, plants vary greatly in their responses to these regulatory compounds. Here, we use genetic perturbation of endogenous gibberellin levels to probe the extent of intraspecific variation in gibberellin responses in natural accessions of Arabidopsis (Arabidopsis thaliana). We find that these accessions vary greatly in their ability to buffer the effects of overexpression of GA20ox1, encoding a rate-limiting enzyme for gibberellin biosynthesis, with substantial differences in bioactive gibberellin concentrations as well as transcriptomes and growth trajectories. These findings demonstrate a surprising level of flexibility in the wiring of regulatory networks underlying hormone metabolism and signaling.

Emergent adaptive behaviour of GRN-controlled simulated robots in a changing environment.

We developed a bio-inspired robot controller combining an artificial genome with an agent-based control system. The genome encodes a gene regulatory network (GRN) that is switched on by environmental cues and, following the rules of transcriptional regulation, provides output signals to actuators. Whereas the genome represents the full encoding of the transcriptional network, the agent-based system mimics the active regulatory network and signal transduction system also present in naturally occurring biological systems. Using such a design that separates the static from the conditionally active part of the gene regulatory network contributes to a better general adaptive behaviour. Here, we have explored the potential of our platform with respect to the evolution of adaptive behaviour, such as preying when food becomes scarce, in a complex and changing environment and show through simulations of swarm robots in an A-life environment that evolution of collective behaviour likely can be attributed to bio-inspired evolutionary processes acting at different levels, from the gene and the genome to the individual robot and robot population.