Modeling microbes: New methods for integrated metabolic and regulatory network reconstruction

PhD Thesis in Bioengineering

Evolution of the gene regulatory network controlling flower dorsoventral asymmetry

Tese de Doutoramento em Biologia de Plantas.

Reconstruction of the regulatory network for Bacillus subtilis and reconciliation with gene expression data

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb. 2016.00275

Unraveling the complex network of cuticular structure and function.

A hydrophobic cuticle is deposited at the outermost extracellular matrix of the epidermis in primary tissues of terrestrial plants. Besides forming a protective shield against the environment, the cuticle is potentially involved in several developmental processes during plant growth. A high degree of variation in cuticle composition and structure exists between different plant species and tissues. Lots of progress has been made recently in understanding the different steps of biosynthesis, transport, and deposition of cuticular components. However, the molecular mechanisms that underlie cuticular function remain largely elusive.

GLUT2 surface expression and intracellular transport via the constitutive pathway in pancreatic beta cells and insulinoma: evidence for a block in trans-Golgi network exit by brefeldin A.

The biosynthesis, intracellular transport, and surface expression of the beta cell glucose transporter GLUT2 was investigated in isolated islets and insulinoma cells. Using a trypsin sensitivity assay to measure cell surface expression, we determined that: (a) greater than 95% of GLUT2 was expressed on the plasma membrane; (b) GLUT2 did not recycle in intracellular vesicles; and (c) after trypsin treatment, reexpression of the intact transporter occurred with a t1/2 of approximately 7 h. Kinetics of intracellular transport of GLUT2 was investigated in pulse-labeling experiments combined with glycosidase treatment and the trypsin sensitivity assay. We determined that transport from the endoplasmic reticulum to the trans-Golgi network (TGN) occurred with a t1/2 of 15 min and that transport from the TGN to the plasma membrane required a similar half-time. When added at the start of a pulse-labeling experiment, brefeldin A prevented exit of GLUT2 from the endoplasmic reticulum. When the transporter was first accumulated in the TGN during a 15-min period of chase, but not following a low temperature (22 degrees C) incubation, addition of brefeldin A (BFA) prevented subsequent surface expression of the transporter. This indicated that brefeldin A prevented GLUT2 exit from the TGN by acting at a site proximal to the 22 degrees C block. Together, these data demonstrate that GLUT2 surface expression in beta cells is via the constitutive pathway, that transport can be blocked by BFA at two distinct steps and that once on the surface, GLUT2 does not recycle in intracellular vesicles.

Testing of Diagnostic Methods for Detection of Influenza Virusfor Optimal Performance in the Context of an Influenza Surveillance Network

Influenza surveillance networks must detect early the viruses that will cause the forthcoming annual epidemics and isolate the strains for further characterization. We obtained the highest sensitivity (95.4%) with a diagnostic tool that combined a shell-vial assay and reverse transcription-PCR on cell culture supernatants at 48 h, and indeed, recovered the strain

Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.

Small cell (neuroendocrine) carcinoma of the prostate: etiology, diagnosis, prognosis, and therapeutic implications--a retrospective study of 30 patients from the rare cancer network.

Within the framework of the Rare Cancer Network Study, we examined 30 patients suffering from small cell neuroendocrine prostate cancer, either in an early/localized or an advanced/metastatic stage. Patients were treated with cisplatin-based chemotherapy, with or without pelvic radiotherapy. Two patients with early disease achieved complete remission for a duration of 19 and 22 months. Three patients with advanced disease achieved complete remission for 6, 7, and 54 months, respectively. Twenty-five patients succumbed to massive local and/or distant failure. No patient presented with brain metastases as the initial site of relapse. Small cell neuroendocrine prostate carcinoma is a very aggressive disease with a poor prognosis, even in its localized form. Despite initial response, the common cisplatin-based chemotherapy plus radiotherapy failed to improve outcome markedly. Improvement will come from understanding the biology of the disease and integrating new targeted therapies into the treatment of this rare and aggressive tumor.

The results of surgery, with or without radiotherapy, for primary spinal myxopapillary ependymoma: a retrospective study from the rare cancer network.

PURPOSE: The aim of this study was to assess the outcome of patients with primary spinal myxopapillary ependymoma (MPE). MATERIALS AND METHODS: Data from a series of 85 (35 females, 50 males) patients with spinal MPE were collected in this retrospective multicenter study. Thirty-eight (45%) underwent surgery only and 47 (55%) received postoperative radiotherapy (RT). Median administered radiation dose was 50.4 Gy (range, 22.2-59.4). Median follow-up of the surviving patients was 60.0 months (range, 0.2-316.6). RESULTS: The 5-year progression-free survival (PFS) was 50.4% and 74.8% for surgery only and surgery with postoperative low- (<50.4 Gy) or high-dose (>or=50.4 Gy) RT, respectively. Treatment failure was observed in 24 (28%) patients. Fifteen patients presented treatment failure at the primary site only, whereas 2 and 1 patients presented with brain and distant spinal failure only. Three and 2 patients with local failure presented with concomitant spinal distant seeding and brain failure, respectively. One patient failed simultaneously in the brain and spine. Age greater than 36 years (p = 0.01), absence of neurologic symptoms at diagnosis (p = 0.01), tumor size >or=25 mm (p = 0.04), and postoperative high-dose RT (p = 0.05) were variables predictive of improved PFS on univariate analysis. In multivariate analysis, only postoperative high-dose RT was independent predictors of PFS (p = 0.04). CONCLUSIONS: The observed pattern of failure was mainly local, but one fifth of the patients presented with a concomitant spinal or brain component. Postoperative high-dose RT appears to significantly reduce the rate of tumor progression.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Outcome and patterns of failure in testicular lymphoma: a multicenter Rare Cancer Network study.

PURPOSE: To assess the outcome and patterns of failure in patients with testicular lymphoma treated by chemotherapy (CT) and/or radiation therapy (RT). METHODS AND MATERIALS: Data from a series of 36 adult patients with Ann Arbor Stage I (n = 21), II (n = 9), III (n = 3), or IV (n = 3) primary testicular lymphoma, consecutively treated between 1980 and 1999, were collected in a retrospective multicenter study by the Rare Cancer Network. Median age was 64 years (range: 21-91 years). Full staging workup (chest X-ray, testicular ultrasound, abdominal ultrasound, and/or thoracoabdominal computer tomography, bone marrow assessment, full blood count, lactate dehydrogenase, and cerebrospinal fluid evaluation) was completed in 18 (50%) patients. All but one patient underwent orchidectomy, and spermatic cord infiltration was found in 9 patients. Most patients (n = 29) had CT, consisting in most cases of cyclophosphamide, doxorubicin, vincristine, and prednisolone (CHOP) with (n = 17) or without intrathecal CT. External RT was delivered to scrotum alone (n = 12) or testicular, iliac, and para-aortic regions (n = 8). The median RT dose was 31 Gy (range: 20-44 Gy) in a median of 17 fractions (10-24), using a median of 1.8 Gy (range: 1.5-2.5 Gy) per fraction. The median follow-up period was 42 months (range: 6-138 months). RESULTS: After a median period of 11 months (range: 1-76 months), 14 patients presented lymphoma progression, mostly in the central nervous system (CNS) (n = 8). Among the 17 patients who received intrathecal CT, 4 had a CNS relapse (p = NS). No testicular, iliac, or para-aortic relapse was observed in patients receiving RT to these regions. The 5-year overall, lymphoma-specific, and disease-free survival was 47%, 66%, and 43%, respectively. In univariate analyses, statistically significant factors favorably influencing the outcome were early-stage and combined modality treatment. Neither RT technique nor total dose influenced the outcome. Multivariate analysis revealed that the most favorable independent factors predicting the outcome were younger age, early-stage disease, and combined modality treatment. CONCLUSIONS: In this multicenter retrospective study, CNS was found to be the principal site of relapse, and no extra-CNS lymphoma progression was observed in the irradiated volumes. More effective CNS prophylaxis, including combined modalities, should be prospectively explored in this uncommon site of extranodal lymphoma.

Primary mucosa-associated lymphoid tissue lymphoma of the salivary glands: a multicenter Rare Cancer Network study.

Purpose: Involvement of salivary glands with mucosa-associated lymphoid tissue (MALT) lymphoma is rare. This retrospective study was performed to assess the clinical profile, treatment outcome, and prognostic factors of MALT lymphoma of the salivary glands.Methods and Materials: Thirteen member centers of the Rare Cancer Network from 10 countries participated, providing data on 63 patients. The median age was 58 years; 47 patients were female and 16 were male. The parotid glands were involved in 49 cases, submandibular in 15, and minor glands in 3. Multiple glands were involved in 9 patients. Staging was as follows: IE in 34, IIE in 12, IIIE in 2, and IV in 15 patients.Results: Surgery (S) alone was performed in 9, radiotherapy (RI) alone in 8, and chemotherapy (CT) alone in 4 patients. Forty-one patients received combined modality treatment (S + RT in 23, S + CT in 8, RT + CT in 4, and all three modalities in 6 patients). No active treatment was given in one case. After initial treatment there was no tumor in 57 patients and residual tumor in 5. Tumor progression was observed in 23 (36.5%) (local in 1, other salivary glands in 10, lymph nodes in 11, and elsewhere in 6). Five patients died of disease progression and the other 5 of other causes. The 5-year disease-free survival, disease-specific survival, and overall survival were 54.4%, 93.2%, and 81.7%, respectively. Factors influencing disease-free survival were use of RI, stage, and residual tumor (p < 0.01). Factors influencing disease-specific survival were stage, recurrence, and residual tumor (p < 0.01).Conclusions: To our knowledge, this report represents the largest series of MALT lymphomas of the salivary glands published to date. This disease may involve all salivary glands either initially or subsequently in 30% of patients. Recurrences may occur in up to 35% of patients at 5 years; however, survival is not affected. Radiotherapy is the only treatment modality that improves disease-free survival. (C) 2012 Elsevier Inc.

A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network.

BACKGROUND: The nuclear receptors are a large family of eukaryotic transcription factors that constitute major pharmacological targets. They exert their combinatorial control through homotypic heterodimerisation. Elucidation of this dimerisation network is vital in order to understand the complex dynamics and potential cross-talk involved. RESULTS: Phylogeny, protein-protein interactions, protein-DNA interactions and gene expression data have been integrated to provide a comprehensive and up-to-date description of the topology and properties of the nuclear receptor interaction network in humans. We discriminate between DNA-binding and non-DNA-binding dimers, and provide a comprehensive interaction map, that identifies potential cross-talk between the various pathways of nuclear receptors. CONCLUSION: We infer that the topology of this network is hub-based, and much more connected than previously thought. The hub-based topology of the network and the wide tissue expression pattern of NRs create a highly competitive environment for the common heterodimerising partners. Furthermore, a significant number of negative feedback loops is present, with the hub protein SHP [NR0B2] playing a major role. We also compare the evolution, topology and properties of the nuclear receptor network with the hub-based dimerisation network of the bHLH transcription factors in order to identify both unique themes and ubiquitous properties in gene regulation. In terms of methodology, we conclude that such a comprehensive picture can only be assembled by semi-automated text-mining, manual curation and integration of data from various sources.

The essential role of radiotherapy in the treatment of Merkel cell carcinoma: a study from the Rare Cancer Network.

PURPOSE: To evaluate the role of postoperative radiotherapy (RT) in Merkel cell carcinoma (MCC). METHODS AND MATERIALS: A retrospective multicenter study was performed in 180 patients with MCC treated between February 1988 and September 2009. Patients who had had surgery alone were compared with patients who received surgery and postoperative RT or radical RT. Local relapse-free survival (LRFS), regional relapse-free survival (RRFS), and distant metastasis-free survival (DMFS) rates were assessed together with disease-free survival (DFS), cancer-specific survival (CSS), and overall survival (OS) rates. RESULTS: Seventy-nine patients were male and 101 patients were female, and the median age was 73 years old (range, 38-93 years). The majority of patients had localized disease (n = 146), and the remaining patients had regional lymph node metastasis (n = 34). Forty-nine patients underwent surgery for the primary tumor without postoperative RT to the primary site; the other 131 patients received surgery for the primary tumor, followed by postoperative RT (n = 118) or a biopsy of the primary tumor followed by radical RT (n = 13). Median follow-up was 5 years (range, 0.2-16.5 years). Patients in the RT group had improved LRFS (93% vs. 64%; p < 0.001), RRFS (76% vs. 27%; p < 0.001), DMFS (70% vs. 42%; p = 0.01), DFS (59% vs. 4%; p < 0.001), and CSS (65% vs. 49%; p = 0.03) rates compared to patients who underwent surgery for the primary tumor alone; LRFS, RRFS, DMFS, and DFS rates remained significant with multivariable Cox regression analysis. However OS was not significantly improved by postoperative RT (56% vs. 46%; p = 0.2). CONCLUSIONS: After multivariable analysis, postoperative RT was associated with improved outcome and seems to be an important component in the multimodality treatment of MCC.

Outcome and prognostic factors in olfactory neuroblastoma: a rare cancer network study.

Purpose: To assess the outcome in patients with olfactory neuroblastoma (ONB). Methods and Materials: Seventy-seven patients treated for nonmetastatic ONB between 1971 and 2004 were included. According to Kadish classification, there were 11 patients with Stage A, 29 with Stage B, and 37 with Stage C. T-classification included 9 patients with T1, 26 with T2, 16 with T3, 15 with T4a, and 11 with T4b tumors. Sixty-eight patients presented with N0 (88%) disease. Results: Most of the patients (n = 56, 73 %) benefited from surgery (S), and total excision was possible in 44 patients (R0 in 32, R1 in 13, R2 in 11). All but five patients benefited from RT, and chemotherapy was given in 21(27%). Median follow-up period was 72 months (range, 6-315). The 5-year overall survival (OS), disease-free survival (DES), locoregional control, and local control were 64%, 57%, 62%, and 70%, respectively. In univariate analyses, favorable factors were Kadish A or B disease, T1 T3 tumors, no nodal involvement, curative surgery, R0/R1 resection, and RT-dose 54 Gy or higher. Multivariate analysis revealed that the best independent factors predicting the outcome were T1 T3, N0, R0/R1 resection, and total RT dose (54 Gy or higher). Conclusion: In this multicenter retrospective study, patients with ONB treated with R0 or R1 surgical resection followed by at least 54-Gy postoperative RT had the best outcome. Novel strategies including concomitant chemotherapy and/or higher dose RT should be prospectively investigated in this rare disease for which local failure remains a problem.