RESUMO

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

Assuntos

RESUMO

Sulfur (S) is an essential macronutrient for plants and its availability in soils is an important determinant for growth and development. Current regulatory policies aimed at reducing industrial S emissions together with changes in agronomical practices have led to a decline in S contents in soils worldwide. Deficiency of sulfate-the primary form of S accessible to plants in soil-has adverse effects on both crop yield and nutritional quality. Hence, recent research has increasingly focused on unraveling the molecular mechanisms through which plants detect and adapt to a limiting supply of sulfate. A significant part of these studies involves the use of omics technologies and has generated comprehensive catalogs of sulfate deficiency-responsive genes and processes, principally in Arabidopsis together with a few studies centering on crop species such as wheat, rice, or members of the Brassica genus. Although we know that sulfate deficiency elicits an important reprogramming of the transcriptome, the transcriptional regulators orchestrating this response are not yet well understood. In this review, we summarize our current knowledge of gene expression responses to sulfate deficiency and recent efforts towards the identification of the transcription factors that are involved in controlling these responses. We further compare the transcriptional response and putative regulators between Arabidopsis and two important crop species, rice and tomato, to gain insights into common mechanisms of the response to sulfate deficiency.

Assuntos

RESUMO

Alternative polyadenylation (APA) increases transcript diversity through the generation of isoforms with varying 3' untranslated region (3' UTR) lengths. As the 3' UTR harbors regulatory element target sites, such as miRNAs or RNA-binding proteins, changes in this region can impact post-transcriptional regulation and translation. Moreover, the APA landscape can change based on the cell type, cell state, or condition. Given that APA events can impact protein expression, investigating translational control is crucial for comprehending the overall cellular regulation process. Revisiting data from polysome profiling followed by RNA sequencing, we investigated the cardiomyogenic differentiation of pluripotent stem cells by identifying the transcripts that show dynamic 3' UTR lengthening or shortening, which are being actively recruited to ribosome complexes. Our findings indicate that dynamic 3' UTR lengthening is not exclusively associated with differential expression during cardiomyogenesis but rather with recruitment to polysomes. We confirm that the differentiated state of cardiomyocytes shows a preference for shorter 3' UTR in comparison to the pluripotent stage although preferences vary during the days of the differentiation process. The most distinct regulatory changes are seen in day 4 of differentiation, which is the mesoderm commitment time point of cardiomyogenesis. After identifying the miRNAs that would target specifically the alternative 3' UTR region of the isoforms, we constructed a gene regulatory network for the cardiomyogenesis process, in which genes related to the cell cycle were identified. Altogether, our work sheds light on the regulation and dynamic 3' UTR changes of polysome-recruited transcripts that take place during the cardiomyogenic differentiation of pluripotent stem cells.

RESUMO

BACKGROUND: Understanding the molecular underpinnings of phenotypic variations is critical for enhancing poultry breeding programs. The Brazilian broiler (TT) and laying hen (CC) lines exhibit striking differences in body weight, growth potential, and muscle mass. Our work aimed to compare the global transcriptome of wing and pectoral tissues during the early development (days 2.5 to 3.5) of these chicken lines, unveiling disparities in gene expression and regulation. RESULTS: Different and bona-fide transcriptomic profiles were identified for the compared lines. A similar number of up- and downregulated differentially expressed genes (DEGs) were identified, considering the broiler line as a reference. Upregulated DEGs displayed an enrichment of protease-encoding genes, whereas downregulated DEGs exhibited a prevalence of receptors and ligands. Gene Ontology analysis revealed that upregulated DEGs were mainly associated with hormone response, mitotic cell cycle, and different metabolic and biosynthetic processes. In contrast, downregulated DEGs were primarily linked to communication, signal transduction, cell differentiation, and nervous system development. Regulatory networks were constructed for the mitotic cell cycle and cell differentiation biological processes, as their contrasting roles may impact the development of distinct postnatal traits. Within the mitotic cell cycle network, key upregulated DEGs included CCND1 and HSP90, with central regulators being NF-κB subunits (RELA and REL) and NFATC2. The cell differentiation network comprises numerous DEGs encoding transcription factors (e.g., HOX genes), receptors, ligands, and histones, while the main regulatory hubs are CREB, AR and epigenetic modifiers. Clustering analyses highlighted PIK3CD as a central player within the differentiation network. CONCLUSIONS: Our study revealed distinct developmental transcriptomes between Brazilian broiler and layer lines. The gene expression profile of broiler embryos seems to favour increased cell proliferation and delayed differentiation, which may contribute to the subsequent enlargement of pectoral tissues during foetal and postnatal development. Our findings pave the way for future functional studies and improvement of targeted traits of economic interest in poultry.

Assuntos

RESUMO

Studying gene regulatory networks associated with cancer provides valuable insights for therapeutic purposes, given that cancer is fundamentally a genetic disease. However, as the number of genes in the system increases, the complexity arising from the interconnections between network components grows exponentially. In this study, using Boolean logic to adjust the existing relationships between network components has facilitated simplifying the modeling process, enabling the generation of attractors that represent cell phenotypes based on breast cancer RNA-seq data. A key therapeutic objective is to guide cells, through targeted interventions, to transition from the current cancer attractor to a physiologically distinct attractor unrelated to cancer. To achieve this, we developed a computational method that identifies network nodes whose inhibition can facilitate the desired transition from one tumor attractor to another associated with apoptosis, leveraging transcriptomic data from cell lines. To validate the model, we utilized previously published in vitro experiments where the downregulation of specific proteins resulted in cell growth arrest and death of a breast cancer cell line. The method proposed in this manuscript combines diverse data sources, conducts structural network analysis, and incorporates relevant biological knowledge on apoptosis in cancer cells. This comprehensive approach aims to identify potential targets of significance for personalized medicine.

Assuntos

RESUMO

The microbiome has shown a correlation with the diet and lifestyle of each population in health and disease, the ability to communicate at the cellular level with the host through innate and adaptative immune receptors, and therefore an important role in modulating inflammatory process related to the establishment and progression of cancer. The oral cavity is one of the most important interaction windows between the human body and the environment, allowing the entry of an important number of microorganisms and their passage across the gastrointestinal tract and lungs. In this review, the contribution of the microbiome network to the establishment of systemic diseases like cancer is analyzed through their synergistic interactions and bidirectional crosstalk in the oral-gut-lung axis as well as its communication with the host cells. Moreover, the impact of the characteristic microbiota of each population in the formation of the multiomics molecular metafirm of the oral-gut-lung axis is also analyzed through state-of-the-art sequencing techniques, which allow a global study of the molecular processes involved of the flow of the microbiota environmental signals through cancer-related cells and its relationship with the establishment of the transcription factor network responsible for the control of regulatory processes involved with tumorigenesis.

Assuntos

RESUMO

Introduction: Pseudomonas aeruginosa infections are one of the leading causes of death in immunocompromised patients with cystic fibrosis, diabetes, and lung diseases such as pneumonia and bronchiectasis. Furthermore, P. aeruginosa is one of the main multidrug-resistant bacteria responsible for nosocomial infections worldwide, including the multidrug-resistant CCBH4851 strain isolated in Brazil. Methods: One way to analyze their dynamic cellular behavior is through computational modeling of the gene regulatory network, which represents interactions between regulatory genes and their targets. For this purpose, Boolean models are important predictive tools to analyze these interactions. They are one of the most commonly used methods for studying complex dynamic behavior in biological systems. Results and discussion: Therefore, this research consists of building a Boolean model of the gene regulatory network of P. aeruginosa CCBH4851 using data from RNA-seq experiments. Next, the basins of attraction are estimated, as these regions and the transitions between them can help identify the attractors, representing long-term behavior in the Boolean model. The essential genes of the basins were associated with the phenotypes of the bacteria for two conditions: biofilm formation and polymyxin B treatment. Overall, the Boolean model and the analysis method proposed in this work can identify promising control actions and indicate potential therapeutic targets, which can help pinpoint new drugs and intervention strategies.

RESUMO

Intestinal gastric cancer (IGC) carcinogenesis results from a complex interplay between environmental and molecular factors, ultimately contributing to disease development. We used integrative bioinformatic analysis to investigate IGC high-throughput molecular data to uncover interactions among differentially expressed genes, microRNAs, and proteins and their roles in IGC. An integrated network was generated based on experimentally validated microRNA-gene/protein interaction data, with three regulatory circuits involved in a complex network contributing to IGC progression. Key regulators were determined, including 23 microRNA and 15 gene/protein hubs. The regulatory circuit networks were associated with hallmarks of cancer, e.g., cell death, apoptosis and the cell cycle, the immune response, and epithelial-to-mesenchymal transition, indicating that different mechanisms of gene regulation impact similar biological functions. Altered expression of hubs was related to the clinicopathological characteristics of IGC patients and showed good performance in discriminating tumors from adjacent nontumor tissues and in relation to T stage and overall survival (OS). Interestingly, expression of upregulated hub hsa-mir-200b and its downregulated target hub gene/protein CFL2 were related not only to pathological T staging and OS but also to changes during IGC carcinogenesis. Our study suggests that regulation of CFL2 by hsa-miR-200b is a dynamic process during tumor progression and that this control plays essential roles in IGC development. Overall, the results indicate that this regulatory interaction is an important component in IGC pathogenesis. Also, we identified a novel molecular interplay between microRNAs, proteins, and genes associated with IGC in a complex biological network and the hubs closely related to IGC carcinogenesis as potential biomarkers.

RESUMO

The Orchidaceae is a mega-diverse plant family with ca. 29,000 species with a large variety of life forms that can colonize transitory habitats. Despite this diversity, little is known about their flowering integrators in response to specific environmental factors. During the reproductive transition in flowering plants a vegetative apical meristem (SAM) transforms into an inflorescence meristem (IM) that forms bracts and flowers. In model grasses, like rice, a flowering genetic regulatory network (FGRN) controlling reproductive transitions has been identified, but little is known in the Orchidaceae. In order to analyze the players of the FRGN in orchids, we performed comprehensive phylogenetic analyses of CONSTANS-like/CONSTANS-like 4 (COL/COL4), FLOWERING LOCUS D (FD), FLOWERING LOCUS C/FRUITFULL (FLC/FUL) and SUPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene lineages. In addition to PEBP and AGL24/SVP genes previously analyzed, here we identify an increase of orchid homologs belonging to COL4, and FUL gene lineages in comparison with other monocots, including grasses, due to orchid-specific gene lineage duplications. Contrariwise, local duplications in Orchidaceae are less frequent in the COL, FD and SOC1 gene lineages, which points to a retention of key functions under strong purifying selection in essential signaling factors. We also identified changes in the protein sequences after such duplications, variation in the evolutionary rates of resulting paralogous clades and targeted expression of isolated homologs in different orchids. Interestingly, vernalization-response genes like VERNALIZATION1 (VRN1) and FLOWERING LOCUS C (FLC) are completely lacking in orchids, or alternatively are reduced in number, as is the case of VERNALIZATION2/GHD7 (VRN2). Our findings point to non-canonical factors sensing temperature changes in orchids during reproductive transition. Expression data of key factors gathered from Elleanthus auratiacus, a terrestrial orchid in high Andean mountains allow us to characterize which copies are actually active during flowering. Altogether, our data lays down a comprehensive framework to assess gene function of a restricted number of homologs identified more likely playing key roles during the flowering transition, and the changes of the FGRN in neotropical orchids in comparison with temperate grasses.

RESUMO

Gene regulatory networks are graph models representing cellular transcription events. Networks are far from complete due to time and resource consumption for experimental validation and curation of the interactions. Previous assessments have shown the modest performance of the available network inference methods based on gene expression data. Here, we study several caveats on the inference of regulatory networks and methods assessment through the quality of the input data and gold standard, and the assessment approach with a focus on the global structure of the network. We used synthetic and biological data for the predictions and experimentally-validated biological networks as the gold standard (ground truth). Standard performance metrics and graph structural properties suggest that methods inferring co-expression networks should no longer be assessed equally with those inferring regulatory interactions. While methods inferring regulatory interactions perform better in global regulatory network inference than co-expression-based methods, the latter is better suited to infer function-specific regulons and co-regulation networks. When merging expression data, the size increase should outweigh the noise inclusion and graph structure should be considered when integrating the inferences. We conclude with guidelines to take advantage of inference methods and their assessment based on the applications and available expression datasets.

RESUMO

Context: Inferring gene regulatory networks (GRN) from high-throughput gene expression data is a challenging task for which different strategies have been developed. Nevertheless, no ever-winning method exists, and each method has its advantages, intrinsic biases, and application domains. Thus, in order to analyze a dataset, users should be able to test different techniques and choose the most appropriate one. This step can be particularly difficult and time consuming, since most methods' implementations are made available independently, possibly in different programming languages. The implementation of an open-source library containing different inference methods within a common framework is expected to be a valuable toolkit for the systems biology community. Results: In this work, we introduce GReNaDIne (Gene Regulatory Network Data-driven Inference), a Python package that implements 18 machine learning data-driven gene regulatory network inference methods. It also includes eight generalist preprocessing techniques, suitable for both RNA-seq and microarray dataset analysis, as well as four normalization techniques dedicated to RNA-seq. In addition, this package implements the possibility to combine the results of different inference tools to form robust and efficient ensembles. This package has been successfully assessed under the DREAM5 challenge benchmark dataset. The open-source GReNaDIne Python package is made freely available in a dedicated GitLab repository, as well as in the official third-party software repository PyPI Python Package Index. The latest documentation on the GReNaDIne library is also available at Read the Docs, an open-source software documentation hosting platform. Contribution: The GReNaDIne tool represents a technological contribution to the field of systems biology. This package can be used to infer gene regulatory networks from high-throughput gene expression data using different algorithms within the same framework. In order to analyze their datasets, users can apply a battery of preprocessing and postprocessing tools and choose the most adapted inference method from the GReNaDIne library and even combine the output of different methods to obtain more robust results. The results format provided by GReNaDIne is compatible with well-known complementary refinement tools such as PYSCENIC.

Assuntos

RESUMO

The transcriptomic analysis of microarray and RNA-Seq datasets followed our own bioinformatic pipeline to identify a transcriptional regulatory network of lung cancer. Twenty-six transcription factors are dysregulated and co-expressed in most of the lung cancer and pulmonary arterial hypertension datasets, which makes them the most frequently dysregulated transcription factors. Co-expression, gene regulatory, coregulatory, and transcriptional regulatory networks, along with fibration symmetries, were constructed to identify common connection patterns, alignments, main regulators, and target genes in order to analyze transcription factor complex formation, as well as its synchronized co-expression patterns in every type of lung cancer. The regulatory function of the most frequently dysregulated transcription factors over lung cancer deregulated genes was validated with ChEA3 enrichment analysis. A Kaplan-Meier plotter analysis linked the dysregulation of the top transcription factors with lung cancer patients' survival. Our results indicate that lung cancer has unique and common deregulated genes and transcription factors with pulmonary arterial hypertension, co-expressed and regulated in a coordinated and cooperative manner by the transcriptional regulatory network that might be associated with critical biological processes and signaling pathways related to the acquisition of the hallmarks of cancer, making them potentially relevant tumor biomarkers for lung cancer early diagnosis and targets for the development of personalized therapies against lung cancer.

RESUMO

Piscirickettsia salmonis is the most important health problem facing Chilean Aquaculture. Previous reports suggest that P. salmonis can survive in salmonid macrophages by interfering with the host immune response. However, the relevant aspects of the molecular pathogenesis of P. salmonis have been poorly characterized. In this work, we evaluated the transcriptomic changes in macrophage-like cell line SHK-1 infected with P. salmonis at 24- and 48-hours post-infection (hpi) and generated network models of the macrophage response to the infection using co-expression analysis and regulatory transcription factor-target gene information. Transcriptomic analysis showed that 635 genes were differentially expressed after 24- and/or 48-hpi. The pattern of expression of these genes was analyzed by weighted co-expression network analysis (WGCNA), which classified genes into 4 modules of expression, comprising early responses to the bacterium. Induced genes included genes involved in metabolism and cell differentiation, intracellular transportation, and cytoskeleton reorganization, while repressed genes included genes involved in extracellular matrix organization and RNA metabolism. To understand how these expression changes are orchestrated and to pinpoint relevant transcription factors (TFs) controlling the response, we established a curated database of TF-target gene regulatory interactions in Salmo salar, SalSaDB. Using this resource, together with co-expression module data, we generated infection context-specific networks that were analyzed to determine highly connected TF nodes. We found that the most connected TF of the 24- and 48-hpi response networks is KLF17, an ortholog of the KLF4 TF involved in the polarization of macrophages to an M2-phenotype in mammals. Interestingly, while KLF17 is induced by P. salmonis infection, other TFs, such as NOTCH3 and NFATC1, whose orthologs in mammals are related to M1-like macrophages, are repressed. In sum, our results suggest the induction of early regulatory events associated with an M2-like phenotype of macrophages that drives effectors related to the lysosome, RNA metabolism, cytoskeleton organization, and extracellular matrix remodeling. Moreover, the M1-like response seems delayed in generating an effective response, suggesting a polarization towards M2-like macrophages that allows the survival of P. salmonis. This work also contributes to SalSaDB, a curated database of TF-target gene interactions that is freely available for the Atlantic salmon community.

Assuntos

RESUMO

Cyclic attractors generated from Boolean models may explain the adaptability of a cell in response to a dynamical complex tumor microenvironment. In contrast to this idea, we postulate that cyclic attractors in certain cases could be a systemic mechanism to face the perturbations coming from the environment. To justify our conjecture, we present a dynamic analysis of a highly curated transcriptional regulatory network of macrophages constrained into a cancer microenvironment. We observed that when M1-associated transcription factors (STAT1 or NF-κB) are perturbed and the microenvironment balances to a hyper-inflammation condition, cycle attractors activate genes whose signals counteract this effect implicated in tissue damage. The same behavior happens when the M2-associated transcription factors are disturbed (STAT3 or STAT6); cycle attractors will prevent a hyper-regulation scenario implicated in providing a suitable environment for tumor growth. Therefore, here we propose that cyclic macrophage phenotypes can serve as a reservoir for balancing the phenotypes when a specific phenotype-based transcription factor is perturbed in the regulatory network of macrophages. We consider that cyclic attractors should not be simply ignored, but it is necessary to carefully evaluate their biological importance. In this work, we suggest one conjecture: the cyclic attractors can serve as a reservoir to balance the inflammatory/regulatory response of the network under external perturbations.

Assuntos

RESUMO

The use of a new bioinformatics pipeline allowed the identification of deregulated transcription factors (TFs) coexpressed in lung cancer that could become biomarkers of tumor establishment and progression. A gene regulatory network (GRN) of lung cancer was created with the normalized gene expression levels of differentially expressed genes (DEGs) from the microarray dataset GSE19804. Moreover, coregulatory and transcriptional regulatory network (TRN) analyses were performed for the main regulators identified in the GRN analysis. The gene targets and binding motifs of all potentially implicated regulators were identified in the TRN and with multiple alignments of the TFs' target gene sequences. Six transcription factors (E2F3, FHL2, ETS1, KAT6B, TWIST1, and RUNX2) were identified in the GRN as essential regulators of gene expression in non-small-cell lung cancer (NSCLC) and related to the lung tumoral process. Our findings indicate that RUNX2 could be an important regulator of the lung cancer GRN through the formation of coregulatory complexes with other TFs related to the establishment and progression of lung cancer. Therefore, RUNX2 could become an essential biomarker for developing diagnostic tools and specific treatments against tumoral diseases in the lung after the experimental validation of its regulatory function.

RESUMO

BACKGROUND: Research on gene duplication is abundant and comes from a wide range of approaches, from high-throughput analyses and experimental evolution to bioinformatics and theoretical models. Notwithstanding, a consensus is still lacking regarding evolutionary mechanisms involved in evolution through gene duplication as well as the conditions that affect them. We argue that a better understanding of evolution through gene duplication requires considering explicitly that genes do not act in isolation. It demands studying how the perturbation that gene duplication implies percolates through the web of gene interactions. Due to evolution's contingent nature, the paths that lead to the final fate of duplicates must depend strongly on the early stages of gene duplication, before gene copies have accumulated distinctive changes. METHODS: Here we use a widely-known model of gene regulatory networks to study how gene duplication affects network behavior in early stages. Such networks comprise sets of genes that cross-regulate. They organize gene activity creating the gene expression patterns that give cells their phenotypic properties. We focus on how duplication affects two evolutionarily relevant properties of gene regulatory networks: mitigation of the effect of new mutations and access to new phenotypic variants through mutation. RESULTS: Among other observations, we find that those networks that are better at maintaining the original phenotype after duplication are usually also better at buffering the effect of single interaction mutations and that duplication tends to enhance further this ability. Moreover, the effect of mutations after duplication depends on both the kind of mutation and genes involved in it. We also found that those phenotypes that had easier access through mutation before duplication had higher chances of remaining accessible through new mutations after duplication. CONCLUSION: Our results support that gene duplication often mitigates the impact of new mutations and that this effect is not merely due to changes in the number of genes. The work that we put forward helps to identify conditions under which gene duplication may enhance evolvability and robustness to mutations.

Assuntos

RESUMO

We reconstructed a transcriptional regulatory network for adrenocortical carcinoma (ACC) using transcriptomic and clinical data from The Cancer Genome Atlas (TCGA)-ACC cohort. We investigated the association of transcriptional regulatory units (regulons) with overall survival, molecular phenotypes, and immune signatures. We annotated the ACC regulons with cancer hallmarks and assessed single sample regulon activities in the European Network for the Study of Adrenal Tumors (ENSAT) cohort. We found 369 regulons associated with overall survival and subdivided them into four clusters: RC1 and RC2, associated with good prognosis, and RC3 and RC4, associated with worse outcomes. The RC1 and RC3 regulons were highly correlated with the 'Steroid Phenotype,' while the RC2 and RC4 regulons were highly correlated with a molecular proliferation signature. We selected two regulons, NR5A1 (steroidogenic factor 1, SF-1) and CENPA (Centromeric Protein A), that were consistently associated with overall survival for further downstream analyses. The CENPA regulon was the primary regulator of MKI-67 (a marker of proliferation KI-67), while the NR5A1 regulon is a well-described transcription factor (TF) in ACC tumorigenesis. We also found that the ZBTB4 (Zinc finger and BTB domain-containing protein 4) regulon, which is negatively associated with CENPA in our transcriptional regulatory network, is also a druggable anti-tumorigenic TF. We anticipate that the ACC regulons may be used as a reference for further investigations concerning the complex molecular interactions in ACC tumors.

Assuntos

RESUMO

A comparative proteomic study at 6 h of growth in minimal medium (MM) and bacteroids at 18 days of symbiosis of Rhizobium etli CFN42 with the Phaseolus vulgaris leguminous plant was performed. A gene ontology classification of proteins in MM and bacteroid, showed 31 and 10 pathways with higher or equal than 30 and 20% of proteins with respect to genome content per pathway, respectively. These pathways were for energy and environmental compound metabolism, contributing to understand how Rhizobium is adapted to the different conditions. Metabolic maps based on orthology of the protein profiles, showed 101 and 74 functional homologous proteins in the MM and bacteroid profiles, respectively, which were grouped in 34 different isoenzymes showing a great impact in metabolism by covering 60 metabolic pathways in MM and symbiosis. Taking advantage of co-expression of transcriptional regulators (TF's) in the profiles, by selection of genes whose matrices were clustered with matrices of TF's, Transcriptional Regulatory networks (TRN´s) were deduced by the first time for these metabolic stages. In these clustered TF-MM and clustered TF-bacteroid networks, containing 654 and 246 proteins, including 93 and 46 TFs, respectively, showing valuable information of the TF's and their regulated genes with high stringency. Isoenzymes were specific for adaptation to the different conditions and a different transcriptional regulation for MM and bacteroid was deduced. The parameters of the TRNs of these expected biological networks and biological networks of E. coli and B. subtilis segregate from the random theoretical networks. These are useful data to design experiments on TF gene-target relationships for bases to construct a TRN.