Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Lignin biosynthesis in Norway spruce: from a model system to the tree

Lignin is a hydrophobic polymer that is synthesised in the secondary cell walls of all vascular plants. It enables water conduction through the stem, supports the upright growth habit and protects against invading pathogens. In addition, lignin hinders the utilisation of the cellulosic cell walls of plants in pulp and paper industry and as forage. Lignin precursors are synthesised in the cytoplasm through the phenylpropanoid pathway, transported into the cell wall and oxidised by peroxidases or laccases to phenoxy radicals that couple to form the lignin polymer. This study was conducted to characterise the lignin biosynthetic pathway in Norway spruce (Picea abies (L.) Karst.). We focused on the less well-known polymerisation stage, to identify the enzymes and the regulatory mechanisms that are involved. Available data for lignin biosynthesis in gymnosperms is scarce and, for example, the latest improvements in precursor biosynthesis have only been verified in herbaceous plants. Therefore, we also wanted to study in detail the roles of individual gene family members during developmental and stress-induced lignification, using EST sequencing and real-time RT-PCR. We used, as a model, a Norway spruce tissue culture line that produces extracellular lignin into the culture medium, and showed that lignin polymerisation in the tissue culture depends on peroxidase activity. We identified in the culture medium a significant NADH oxidase activity that could generate H2O2 for peroxidases. Two basic culture medium peroxidases were shown to have high affinity to coniferyl alcohol. Conservation of the putative substrate-binding amino acids was observed when the spruce peroxidase sequences were compared with other peroxidases with high affinity to coniferyl alcohol. We also used different peroxidase fractions to produce synthetic in vitro lignins from coniferyl alcohol; however, the linkage pattern of the suspension culture lignin could not be reproduced in vitro with the purified peroxidases, nor with the full complement of culture medium proteins. This emphasised the importance of the precursor radical concentration in the reaction zone, which is controlled by the cells through the secretion of both the lignin precursors and the oxidative enzymes to the apoplast. In addition, we identified basic peroxidases that were reversibly bound to the lignin precipitate. They could be involved, for example, in the oxidation of polymeric lignin, which is required for polymer growth. The dibenzodioxocin substructure was used as a marker for polymer oxidation in the in vitro polymerisation studies, as it is a typical substructure in wood lignin and in the suspension culture lignin. Using immunolocalisation, we found the structure mainly in the S2+S3 layers of the secondary cell walls of Norway spruce tracheids. The structure was primarily formed during the late phases of lignification. Contrary to the earlier assumptions, it appears to be a terminal structure in the lignin macromolecule. Most lignin biosynthetic enzymes are encoded for by several genes, all of which may not participate in lignin biosynthesis. In order to identify the gene family members that are responsible for developmental lignification, ESTs were sequenced from the lignin-forming tissue culture and developing xylem of spruce. Expression of the identified lignin biosynthetic genes was studied using real-time RT-PCR. Candidate genes for developmental lignification were identified by a coordinated, high expression of certain genes within the gene families in all lignin-forming tissues. However, such coordinated expression was not found for peroxidase genes. We also studied stress-induced lignification either during compression wood formation by bending the stems or after Heterobasidion annosum infection. Based on gene expression profiles, stress-induced monolignol biosynthesis appeared similar to the developmental process, and only single PAL and C3H genes were specifically up-regulated by stress. On the contrary, the up-regulated peroxidase genes differed between developmental and stress-induced lignification, indicating specific responses.

Inhibition of Thrombin in Cardiac Surgery : experiments in a porcine model

Cardiac surgery involving cardiopulmonary bypass (CPB) induces activation of inflammation and coagulation systems and is associated with ischemia-reperfusion injury (I/R injury)in various organs including the myocardium, lungs, and intestine. I/R injury is manifested as organ dysfunction. Thrombin, the key enzyme of coagulation , plays a cenral role also in inflammation and contributes to regulation of apoptosis as well. The general aim of this thesis was to evaluate the potential of thrombin inhibition in reducing the adverse effects of I/R injury in myocardium, lungs, and intestine associated with the use of CPB and cardiac surgery. Forty five pigs were used for the studies. Two randomized blinded studies were performed. Animals underwent 75 min of normothermic CPB, 60 min of aortic clamping, and 120 min of reperfusion period. Twenty animals received iv. recombinant hirudin, a selective and effective inbitor of thrombin, or placebo. In a similar setting, twenty animals received an iv-bolus (250 IU/kg) of antithrombin (AT) or placebo. An additional group of 5 animals received 500 IU/kg in an open label setting to test dose response. Generation of thrombin (TAT), coagulation status (ACT), and hemodynamics were measured. Intramucosal pH and pCO2 were measured from the luminal surface of ileum using tonometry simultaneusly with arterial gas analysis. In addition, myocardial, lung, and intestinal biopsies were taken to quantitate leukocyte infiltration (MPO), for histological evaluation, and detection of apoptosis (TUNEL, caspase 3). In conclusion, our data suggest that r-hirudin may be an effective inhibitor of reperfusion induced thrombin generation in addition to being a direct inhibitor of preformed thrombin. Overall, the results suggest that inhibition of thrombin, beyond what is needed for efficient anticoagulation by heparin, has beneficial effects on myocardial I/R injury and hemodynamics during cardiac surgery and CPB. We showed that infusion of the thrombin inhibitor r-hirudin during reperfusion was associated with attenuated post ischemia left ventricular dysfunction and decreased systemic vascular resistance. Consequently microvascular flow was improved during ischemia-reperfusion injury. Improved recovery of myocardium during the post-ischemic reperfusion period was associated with significantly less cardiomyocyte apoptosis and with a trend in anti-inflammatory effects. Thus, inhibition of reperfusion induced thrombin may offer beneficial effects by mechanisms other than direct anticoagulant effects. AT, in doses with a significant anticoagulant effect, did not alleviate myocardial I/R injury in terms of myocardial recovery, histological inflammatory changes or post-ischemic troponin T release. Instead, AT attenuated reperfusion induced increase in pulmonary pressure after CPB. Taken the clinical significance of postoperative pulmonary hemodynamics in patients undergoing cardiopulmonary bypass, the potential positive regulatory role of AT and clinical implications needs to be studied further. Inflammatory response in the gut wall proved to be poorly associated with perturbed mucosal perfusion and the animals with the least neutrophil tissue sequestration and I/R related histological alterations tended to have the most progressive mucosal hypoperfusion. Thus, mechanisms of low-flow reperfusion injury during CPB can differ from the mechanisms seen in total ischemia reperfusion injury.

Curcumin Increases the Pathogenicity of Salmonella enterica Serovar Typhimurium in Murine Model

Curcumin has gained immense importance for its vast therapeutic and prophylactic applications. Contrary to this, our study reveals that it regulates the defense pathways of Salmonella enterica serovar Typhimurium ( S. Typhimurium) to enhance its pathogenicity. In a murine model of typhoid fever, we observed higher bacterial load in Peyer's,patches, mesenteric lymph node, spleen and liver, when infected with curcumin-treated Salmonella. Curcumin increased the resistance of S. Typhimurium against antimicrobial agents like antimicrobial peptides, reactive oxygen and nitrogen species. This increased tolerance might be attributed to the up-regulation of genes involved in resistance against antimicrobial peptides - pmrD and pmrHFIJKLM and genes with antioxidant function - mntH, sodA and sitA. We implicate that iron chelation property of curcumin have a role in regulating mntH and sitA. Interestingly, we see that the curcumin-mediated modulation of pmr genes is through the PhoPQ regulatory system. Curcumin downregulates SPI1 genes, required for entry into epithelial cells and upregulates SPI2 genes required to intracellular survival. Since it is known that the SPI1 and SPI2 system can be regulated by the PhoPQ system, this common regulator could explain curcumin's mode of action. This data urges us to rethink the indiscriminate use of curcumin especially during Salmonella outbreaks.

Intravenous immunoglobulin expands regulatory T cells via induction of cyclooxygenase-2-dependent prostaglandin E2 in human dendritic cells

CD4(+)CD25(+)FoxP3(+) regulatory T cells (Tregs) play a critical role in the maintenance of immune tolerance. Intravenous immunoglobulin (IVIg), a therapeutic preparation of normal pooled human IgG, expands Tregs in various experimental models and in patients. However, the cellular and molecular mechanisms by which IVIg expands Tregs are relatively unknown. As Treg expansion in the periphery requires signaling by antigen-presenting cells such as dendritic cells (DCs) and IVIg has been demonstrated to modulate DC functions, we hypothesized that IVIg induces distinct signaling events in DCs that subsequently mediate Treg expansion. We demonstrate that IVIg expands Tregs via induction of cyclooxygenase (COX)-2-dependent prostaglandin E2 (PGE(2)) in human DCs. However, costimulatory molecules of DCs such as programmed death ligands, OX40 ligand, and inducible T-cell costimulator ligands were not implicated. Inhibition of PGE(2) synthesis by COX-2 inhibitors prevented IVIg-mediated Treg expansion in vitro and significantly diminished IVIg-mediated Treg expansion in vivo and protection from disease in experimental autoimmune encephalomyelitis model. IVIg-mediated COX-2 expression, PGE(2) production, and Treg expansion were mediated in part via interaction of IVIg and F(ab('))(2) fragments of IVIg with DC-specific intercellular adhesion molecule-3-grabbing nonintegrin. Our results thus uncover novel cellular and molecular mechanism by which IVIg expands Tregs.

Endogenous fisheries management in a stochastic model: Why do fishery agencies use TAC

The aim of this paper is to explain under which circumstances using TACs as instrument to manage a fishery along with fishing periods may be interesting from a regulatory point of view. In order to do this, the deterministic analysis of Homans and Wilen (1997)and Anderson (2000) is extended to a stochastic scenario where the resource cannot be measured accurately. The resulting endogenous stochastic model is numerically solved for finding the optimal control rules in the Iberian sardine stock. Three relevant conclusions can be highligted from simulations. First, the higher the uncertainty about the state of the stock is, the lower the probability of closing the fishery is. Second, the use of TACs as management instrument in fisheries already regulated with fishing periods leads to: i) An increase of the optimal season length and harvests, especially for medium and high number of licences, ii) An improvement of the biological and economic variables when the size of the fleet is large; and iii) Eliminate the extinction risk for the resource. And third, the regulator would rather select the number of licences and do not restrict the season length.

USP1 deubiquitinase: cellular functions, regulatory mechanisms and emerging potential as target in cancer therapy

12 p.

Multirating decision model validation: the relevance of the quality of the securitization issues

Rating enables the information asymmetry existing in the issuer-investor relationship to be reduced, particularly for issues with a high degree of complexity, as is the case of securitizations. However, there may be a serious conflict of interest between the issuer’s choice and remuneration of the agency and the credit rating awarded, resulting in lower quality and information power of the published rating. In this paper, we propose an explicative model of the number of ratings requested, by analyzing the relevance of the number of ratings to measure the reliability, where multirating is shown to be associated to the quality, size, liquidity and the degree of information asymmetry relating to the issue. Thus, we consider that the regulatory changes that foster the widespread publication of simultaneous ratings could help to alleviate the problem of rating model arbitrage and the crisis of confidence in credit ratings in general and in the securitization issues, in particular.

Phenotype consequences of myophosphorylase dysfunction: insights from the McArdle mouse model

McArdle disease, caused by inherited deficiency of the enzyme muscle glycogen phosphorylase (GP-MM), is arguably the paradigm of exercise intolerance. The recent knock-in (p.R50X/p.R50X) mouse disease model allows an investigation of the phenotypic consequences of muscle glycogen unavailability and the physiopathology of exercise intolerance. We analysed, in 2-month-old mice [wild-type (wt/wt), heterozygous (p.R50X/wt) and p.R50X/p.R50X)], maximal endurance exercise capacity and the molecular consequences of an absence of GP-MM in the main glycogen metabolism regulatory enzymes: glycogen synthase, glycogen branching enzyme and glycogen debranching enzyme, as well as glycogen content in slow-twitch (soleus), intermediate (gastrocnemius) and glycolytic/fast-twitch (extensor digitorum longus; EDL) muscles.

Targeting the EP1 receptor reduces Fas ligand expression and increases the antitumor immune response in an in vivo model of colon cancer

Despite studies demonstrating that inhibition of cyclooxygenase-2 (COX-2)-derived prostaglandin E2 (PGE2) has significant chemotherapeutic benefits in vitro and in vivo, inhibition of COX enzymes is associated with serious gastrointestinal and cardiovascular side effects, limiting the clinical utility of these drugs. PGE2 signals through four different receptors (EP1–EP4) and targeting individual receptor(s) may avoid these side effects, while retaining significant anticancer benefits. Here, we show that targeted inhibition of the EP1 receptor in the tumor cells and the tumor microenvironment resulted in the significant inhibition of tumor growth in vivo. Both dietary administration and direct injection of the EP1 receptor-specific antagonist, ONO-8713, effectively reduced the growth of established CT26 tumors in BALB/c mice, with suppression of the EP1 receptor in the tumor cells alone less effective in reducing tumor growth. This antitumor effect was associated with reduced Fas ligand expression and attenuated tumor-induced immune suppression. In particular, tumor infiltration by CD4+CD25+Foxp3+ regulatory T cells was decreased, whereas the cytotoxic activity of isolated splenocytes against CT26 cells was increased. F4/80+ macrophage infiltration was also decreased; however, there was no change in macrophage phenotype. These findings suggest that the EP1 receptor represents a potential target for the treatment of colon cancer.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.