Implicit methods for qualitative modeling of gene regulatory networks.

Advancements in high-throughput technologies to measure increasingly complex biological phenomena at the genomic level are rapidly changing the face of biological research from the single-gene single-protein experimental approach to studying the behavior of a gene in the context of the entire genome (and proteome). This shift in research methodologies has resulted in a new field of network biology that deals with modeling cellular behavior in terms of network structures such as signaling pathways and gene regulatory networks. In these networks, different biological entities such as genes, proteins, and metabolites interact with each other, giving rise to a dynamical system. Even though there exists a mature field of dynamical systems theory to model such network structures, some technical challenges are unique to biology such as the inability to measure precise kinetic information on gene-gene or gene-protein interactions and the need to model increasingly large networks comprising thousands of nodes. These challenges have renewed interest in developing new computational techniques for modeling complex biological systems. This chapter presents a modeling framework based on Boolean algebra and finite-state machines that are reminiscent of the approach used for digital circuit synthesis and simulation in the field of very-large-scale integration (VLSI). The proposed formalism enables a common mathematical framework to develop computational techniques for modeling different aspects of the regulatory networks such as steady-state behavior, stochasticity, and gene perturbation experiments.

Gene transfer of a soluble IL-1 type 2 receptor-Ig fusion protein improves cardiac allograft survival in rats.

OBJECTIVE: Interleukin-1 (IL-1) mediates ischemia-reperfusion injury and graft inflammation after heart transplantation. IL-1 affects target cells through two distinct types of transmembrane receptors, type-1 receptor (IL-1R1), which transduces the signal, and the non-signaling type-2 receptor (IL-1R2), which acts as a ligand sink that subtracts IL-1beta from IL-1R1. We analyzed the efficacy of adenovirus (Ad)-mediated gene transfer of a soluble IL-1R2-Ig fusion protein in delaying cardiac allograft rejection and the mechanisms underlying the protective effect. METHODS: IL-1 inhibition by IL-1R2-Ig was tested using an in vitro functional assay whereby endothelial cells preincubated with AdIL-1R2-Ig or control virus were stimulated with recombinant IL-1beta or tumor necrosis factor-alpha (TNF-alpha), and urokinase-type plasminogen activator (u-PA) induction was measured by zymography. AdIL-1R2-Ig was delivered to F344 rat donor hearts ex vivo, which were placed in the abdominal position in LEW hosts. Intragraft inflammatory cell infiltrates and proinflammatory cytokine expression were analyzed by immunohistochemistry and real-time reverse transcriptase-polymerase chain reaction (RT-PCR), respectively. RESULTS: IL-1R2-Ig specifically inhibited IL-1beta-induced u-PA responses in vitro. IL-1R2-Ig gene transfer reduced intragraft monocytes/macrophages and CD4(+) cell infiltrates (p<0.05), TNF-alpha and transforming growth factor-beta (TGF-beta) expression (p<0.05), and prolonged graft survival (15.6+/-5.7 vs 10.3+/-2.5 days with control vector and 10.1+/-2.1 days with buffer alone; p<0.01). AdIL-1R2-Ig combined with a subtherapeutic regimen of cyclosporin A (CsA) was superior to CsA alone (19.4+/-3.0 vs 15.9+/-1.8 days; p<0.05). CONCLUSIONS: Soluble IL-1 type-2 receptor gene transfer attenuates cardiac allograft rejection in a rat model. IL-1 inhibition may be useful as an adjuvant therapy in heart transplantation.

Rhythmic growth explained by coincidence between internal and external cues.

Most organisms use circadian oscillators to coordinate physiological and developmental processes such as growth with predictable daily environmental changes like sunrise and sunset. The importance of such coordination is highlighted by studies showing that circadian dysfunction causes reduced fitness in bacteria and plants, as well as sleep and psychological disorders in humans. Plant cell growth requires energy and water-factors that oscillate owing to diurnal environmental changes. Indeed, two important factors controlling stem growth are the internal circadian oscillator and external light levels. However, most circadian studies have been performed in constant conditions, precluding mechanistic study of interactions between the clock and diurnal variation in the environment. Studies of stem elongation in diurnal conditions have revealed complex growth patterns, but no mechanism has been described. Here we show that the growth phase of Arabidopsis seedlings in diurnal light conditions is shifted 8-12 h relative to plants in continuous light, and we describe a mechanism underlying this environmental response. We find that the clock regulates transcript levels of two basic helix-loop-helix genes, phytochrome-interacting factor 4 (PIF4) and PIF5, whereas light regulates their protein abundance. These genes function as positive growth regulators; the coincidence of high transcript levels (by the clock) and protein accumulation (in the dark) allows them to promote plant growth at the end of the night. Thus, these two genes integrate clock and light signalling, and their coordinated regulation explains the observed diurnal growth rhythms. This interaction may serve as a paradigm for understanding how endogenous and environmental signals cooperate to control other processes.

A stroma-related gene signature predicts resistance to neoadjuvant chemotherapy in breast cancer.

To better understand the relationship between tumor-host interactions and the efficacy of chemotherapy, we have developed an analytical approach to quantify several biological processes observed in gene expression data sets. We tested the approach on tumor biopsies from individuals with estrogen receptor-negative breast cancer treated with chemotherapy. We report that increased stromal gene expression predicts resistance to preoperative chemotherapy with 5-fluorouracil, epirubicin and cyclophosphamide (FEC) in subjects in the EORTC 10994/BIG 00-01 trial. The predictive value of the stromal signature was successfully validated in two independent cohorts of subjects who received chemotherapy but not in an untreated control group, indicating that the signature is predictive rather than prognostic. The genes in the signature are expressed in reactive stroma, according to reanalysis of data from microdissected breast tumor samples. These findings identify a previously undescribed resistance mechanism to FEC treatment and suggest that antistromal agents may offer new ways to overcome resistance to chemotherapy.

Expression of the mammalian Xenotropic Polytropic Virus Receptor 1 (XPR1) in tobacco leaves leads to phosphate export.

Phosphate homeostasis in multicellular eukaryotes depends on both phosphate influx and efflux. The mammalian Xenotropic Polytropic Virus Receptor 1 (XPR1) shares homology to the Arabidopsis PHO1, a phosphate exporter expressed in roots. However, phosphate export activity of XPR1 has not yet been demonstrated in a heterologous system. Here, wedemonstrate that transient expression in tobacco leaves of XPR1-GFP leads to specific phosphate export. Like PHO1-GFP, XPR1-GFP is localized predominantly to the endomembrane system in tobacco cells. These results show that tobacco leaves are a good heterologous system to study the transport activity of members of the PHO1/XPR1 family.

Expression and secretion of the atrial natriuretic peptide in human adipose tissue and preadipocytes.

OBJECTIVE: Atrial natriuretic peptide (ANP) is a secretory hormone displaying diuretic, natriuretic, and vasorelaxant activities. Recently, its lipolytic activity has been reported. Since the expression of ANP in adipose tissue has not been documented, we used real-time reverse transcriptase polymerase chain reaction (RT-PCR) to investigate the expression of ANP in human adipose tissue and preadipocytes. RESEARCH METHODS AND PROCEDURES: RNA was extracted from the human adipose tissue of severely obese premenopausal women as well as from human preadipocytes. For human preadipocytes, two cell systems were investigated: the human preadipose immortalized (Chub-S7) cells, a well-characterized human preadipose cell line, and primary preadipocytes derived from the stromal vascular fraction of the human adipose tissue. We measured the mRNA of ANP, of corin (a transmembrane serine protease involved in the conversion of pro-ANP to ANP) and of uncoupling protein 2 (UCP2; a control gene known to be ubiquitously expressed). The expression of ANP was also investigated using immunofluorescence and radioimmunoassay in Chub-S7 cells and human primary preadipocytes in culture. RESULTS: Our results indicate that ANP and corin are expressed at the mRNA level in human adipose tissue and preadipocytes. Immunofluorescence experiments demonstrated that pro-ANP was expressed in Chub-S7 cells. In addition, ANP secretion could be measured in Chub-S7 cells and human primary preadipocytes in culture. Rosiglitazone, a selective peroxisome proliferator-activated receptor type gamma (PPAR-gamma) agonist promoting adipocyte differentiation, was found to modulate both ANP expression and secretion in preadipocytes. DISCUSSION: Our findings suggest the existence of an autocrine/paracrine system for ANP in the human adipose tissue whose implications in lipolysis and cardiovascular function need to be further explored.

Influence of HLA-C expression level on HIV control.

A variant upstream of human leukocyte antigen C (HLA-C) shows the most significant genome-wide effect on HIV control in European Americans and is also associated with the level of HLA-C expression. We characterized the differential cell surface expression levels of all common HLA-C allotypes and tested directly for effects of HLA-C expression on outcomes of HIV infection in 5243 individuals. Increasing HLA-C expression was associated with protection against multiple outcomes independently of individual HLA allelic effects in both African and European Americans, regardless of their distinct HLA-C frequencies and linkage relationships with HLA-B and HLA-A. Higher HLA-C expression was correlated with increased likelihood of cytotoxic T lymphocyte responses and frequency of viral escape mutation. In contrast, high HLA-C expression had a deleterious effect in Crohn's disease, suggesting a broader influence of HLA expression levels in human disease.

Pob1 participates in the Cdc42 regulation of fission yeast actin cytoskeleton.

Rho GTPases regulate the actin cytoskeleton in all eukaryotes. Fission yeast Cdc42 is involved in actin cable assembly and formin For3 regulation. We isolated cdc42-879 as a thermosensitive strain with actin cable and For3 localization defects. In a multicopy suppressor screening, we identified pob1(+) as suppressor of cdc42-879 thermosensitivity. Pob1 overexpression also partially restores actin cables and localization of For3 in the mutant strain. Pob1 interacts with Cdc42 and this GTPase regulates Pob1 localization and/or stability. The C-terminal pleckstrin homology (PH) domain of Pob1 is required for Cdc42 binding. Pob1 also binds to For3 through its N-terminal sterile alpha motif (SAM) domain and contributes to the formin localization at the cell tips. The previously described pob1-664 mutant strain (Mol. Biol. Cell. 10, 2745-2757, 1999), which carries a mutation in the PH domain, as well as pob1 mutant strains in which Pob1 lacks the N-terminal region (pob1DeltaN) or the SAM domain (pob1DeltaSAM), have cytoskeletal defects similar to that of cdc42-879 cells. Expression of constitutively active For3DAD* partially restores actin organization in cdc42-879, pob1-664, pob1DeltaN, and pob1DeltaSAM. Therefore, we propose that Pob1 is required for For3 localization to the tips and facilitates Cdc42-mediated relief of For3 autoinhibition to stimulate actin cable formation.

Evolution of gene expression changes in newborn rats after mechanical ventilation with reversible intubation

Résumant mon travail de thèse, l'article qui suit décrit un nouveau modèle animal servant à étudier l'impact combiné d'une ventilation mécanique (VM), d'une oxygénothérapie et d'une inflammation sur des poumons immatures. Cette étude permet, pour la première fois, de mesurer l'expression de gènes à distance d'une VM pour en analyser la cinétique. La VM représente un traitement intégral dans la prise en charge de prématurés. Sauvant des vies, elle est cependant non-physiologique et décrite comme nocive à court et à long terme, empêchant le bon développement pulmonaire. Nombreuses études se sont intéressées à l'impact immédiat de la VM sur les poumons, mais il n'existe à ce jour aucun modèle de rongeur pour en analyser les effets tardifs. Par analogie avec la clinique, nous avons créé un modèle avec un animal dont le stade développemental pulmonaire est comparable aux prématurés humains et consistant en une oxygénothérapie, une VM modérée avec intubation non chirurgicale, similaire à la pratique quotidienne, et un contexte inflammatoire mimant celui de chorioamnionite dans lequel bien des prématurés naissent. Nous avons ensuite réalisé une extubation pour permettre une période de rétablissement, puis fait des analyses et sur le plan structurel par histologie conventionnelle et en 3D, et sur le plan biologique, par analyse de l'expression de gènes et de protéines. Ce travail a permis de valider ce nouveau modèle comme outil de recherche pour réaliser des mesures à distance d'une VM chez des rats nouveau-nés. Comparant ces mesures à celles prises à la fin de la VM, nous observons: une augmentation initiale et transitoire des médiateurs impliqués dans la cascade inflammatoire dont le corrélat histologique est une maladie inflammatoire pulmonaire et, tardivement, une altération plus développementale de la structure pulmonaire avec diminution de l'alvéolarisation. Ceci pourrait être en partie dû à une expression asynchrone de gènes décrits comme importants pour la formation des alvéoles (matrix metalloproteinase 9, elastine). Offrant une nouvelle approche pour la recherche pulmonaire chez les rongeurs, ce modèle servira comme futur outil pour approfondir nos connaissances de la physiopathologie conduisant aux altérations structurelles retrouvées dans les poumons d'anciens prématurés soumis à une VM (dysplasie broncho-pulmonaire), pour tester l'influence de certains traitements (p.ex. surfactant) et pour étudier les effets de la VM en l'appliquant à des modèles transgéniques.

Enhanced cerebral expression of MCT1 and MCT2 in a rat ischemia model occurs in activated microglial cells.

Monocarboxylate transporters (MCTs) are essential for the use of lactate, an energy substrate known to be overproduced in brain during an ischemic episode. The expression of MCT1 and MCT2 was investigated at 48 h of reperfusion from focal ischemia induced by unilateral extradural compression in Wistar rats. Increased MCT1 mRNA expression was detected in the injured cortex and hippocampus of compressed animals compared to sham controls. In the contralateral, uncompressed hemisphere, increases in MCT1 mRNA level in the cortex and MCT2 mRNA level in the hippocampus were noted. Interestingly, strong MCT1 and MCT2 protein expression was found in peri-lesional macrophages/microglia and in an isolectin B4+/S100beta+ cell population in the corpus callosum. In vitro, MCT1 and MCT2 protein expression was observed in the N11 microglial cell line, whereas an enhancement of MCT1 expression by tumor necrosis factor-alpha (TNF-alpha) was shown in these cells. Modulation of MCT expression in microglia suggests that these transporters may help sustain microglial functions during recovery from focal brain ischemia. Overall, our study indicates that changes in MCT expression around and also away from the ischemic area, both at the mRNA and protein levels, are a part of the metabolic adaptations taking place in the brain after ischemia.

Genome-wide prediction of matrix attachment regions that increase gene expression in mammalian cells.

Gene transfer in eukaryotic cells and organisms suffers from epigenetic effects that result in low or unstable transgene expression and high clonal variability. Use of epigenetic regulators such as matrix attachment regions (MARs) is a promising approach to alleviate such unwanted effects. Dissection of a known MAR allowed the identification of sequence motifs that mediate elevated transgene expression. Bioinformatics analysis implied that these motifs adopt a curved DNA structure that positions nucleosomes and binds specific transcription factors. From these observations, we computed putative MARs from the human genome. Cloning of several predicted MARs indicated that they are much more potent than the previously known element, boosting the expression of recombinant proteins from cultured cells as well as mediating high and sustained expression in mice. Thus we computationally identified potent epigenetic regulators, opening new strategies toward high and stable transgene expression for research, therapeutic production or gene-based therapies.

Connexin 43 expression in human hypertrophied heart due to pressure and volume overload.

Left ventricular hypertrophy (LVH) is due to pressure overload or mechanical stretch and is thought to be associated with remodeling of gap-junctions. We investigated whether the expression of connexin 43 (Cx43) is altered in humans in response to different degrees of LVH. The expression of Cx43 was analyzed by quantitative polymerase chain reaction, Western blot analysis and immunohistochemistry on left ventricular biopsies from patients undergoing aortic or mitral valve replacement. Three groups were analyzed: patients with aortic stenosis with severe LVH (n=9) versus only mild LVH (n=7), and patients with LVH caused by mitral regurgitation (n=5). Cx43 mRNA expression and protein expression were similar in the three groups studied. Furthermore, immunohistochemistry revealed no change in Cx43 distribution. We can conclude that when compared with mild LVH or with LVH due to volume overload, severe LVH due to chronic pressure overload is not accompanied by detectable changes of Cx43 expression or spatial distribution.

The direct effects of tacrolimus and cyclosporin A on isolated human islets: A functional, survival and gene expression study.

The use of immunosuppressive drugs in transplanted patients is associated with the development of diabetes, possibly due to β-cell toxicity. To better understand the mechanisms leading to post-transplant diabetes, we investigated the actions of prolonged exposure of isolated human islets to therapeutical levels of tacrolimus (Tac) or cyclosporin A (CsA). Islets were isolated from the pancreas of multiorgan donors by enzymatic digestion and density gradient centrifugation. Functional, survival and molecular studies were then performed after 4 days of incubation with therapeutical concentrations of Tac or  CsA. Glucose-induced insulin secretion was significantly decreased in Tac, but not in CsA exposed islets, which was associated with a reduction of the amount of insulin granules as shown by electron microscopy. The percentage of apoptotic β-cells was higher in Tac than CsA exposed islets. Microarray experiments followed by Gene Set Enrichment Analysis revealed that gene expression was more markedly affected upon Tac treatment. In conclusion, Tac and CsA affect features of beta-cell differently, with several changes occurring at the molecular level.

PAX3/FOXO1 fusion gene status is the key prognostic molecular marker in rhabdomyosarcoma and significantly improves current risk stratification.

PURPOSE: To improve the risk stratification of patients with rhabdomyosarcoma (RMS) through the use of clinical and molecular biologic data. PATIENTS AND METHODS: Two independent data sets of gene-expression profiling for 124 and 101 patients with RMS were used to derive prognostic gene signatures by using a meta-analysis. These and a previously published metagene signature were evaluated by using cross validation analyses. A combined clinical and molecular risk-stratification scheme that incorporated the PAX3/FOXO1 fusion gene status was derived from 287 patients with RMS and evaluated. RESULTS: We showed that our prognostic gene-expression signature and the one previously published performed well with reproducible and significant effects. However, their effect was reduced when cross validated or tested in independent data and did not add new prognostic information over the fusion gene status, which is simpler to assay. Among nonmetastatic patients, patients who were PAX3/FOXO1 positive had a significantly poorer outcome compared with both alveolar-negative and PAX7/FOXO1-positive patients. Furthermore, a new clinicomolecular risk score that incorporated fusion gene status (negative and PAX3/FOXO1 and PAX7/FOXO1 positive), Intergroup Rhabdomyosarcoma Study TNM stage, and age showed a significant increase in performance over the current risk-stratification scheme. CONCLUSION: Gene signatures can improve current stratification of patients with RMS but will require complex assays to be developed and extensive validation before clinical application. A significant majority of their prognostic value was encapsulated by the fusion gene status. A continuous risk score derived from the combination of clinical parameters with the presence or absence of PAX3/FOXO1 represents a robust approach to improving current risk-adapted therapy for RMS.