PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Pathway modelling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48Linked to original sources

Nuclear and mitochondrial DNA repair: similar pathways?

Mitochondrial DNA (mtDNA) alterations are implicated in a broad range of human diseases and alterations of the mitochondrial genome are assumed to be a result of its high susceptibility to oxidative damage and its limited DNA repair compared to nuclear DNA (nDNA). Characterization of DNA repair mechanisms has generally focused on these processes in nDNA but increasing interest and research effort have contributed to our knowledge of the mechanisms underlying DNA repair in mitochondria. In this review, we make comparisons between nDNA and mtDNA repair pathways and propose a model for how these pathways interact in mitochondria.

Animals↗

In vitro analysis of processing at the 3'-end of precursors of M1 RNA, the catalytic subunit of Escherichia coli RNase P: multiple pathways and steps for the processing.

M1 RNA of 377 nucleotides, the catalytic subunit of Escherichia coli RNase P, is produced by a 3' processing reaction from precursor M1 RNA, a major transcript from the rnpB gene. We analyzed products and intermediates generated by the in vitro processing reaction using a 40% ammonium sulfate precipitate of the S30 fraction (ASP-40) and determined their involvement in the processing. From the results we proposed a model of two pathways for 3' processing of M1 RNA. In this model, one pathway (pathway I) involves +385/+386 intermediates and the other pathway (pathway II) does not. The position of the 3'-end of the precursor molecule determined the choice of the pathways. The precursor having the 3'-end of +413 was processed by both pathways while that having the +415 end was processed only by pathway II. The ASP-40 fraction generated processing products (termed +378/+379 RNA) containing one or two more nucleotides at the 3'-end than M1 RNA, regardless of which pathway was used. Therefore, both pathways require the final 3' trimming for complete processing. The endonucleolytic generation of +378/+379 RNA by pathway II was blocked by the rne-3071 mutation, suggesting that this step is carried out by RNase E.

Catalytic Domain↗

InsP3 receptor is essential for growth and differentiation but not for vision in Drosophila.

Phospholipase C (PLC) is the focal point for two major signal transduction pathways: one initiated by G protein-coupled receptors and the other by tyrosine kinase receptors. Active PLC hydrolyzes phosphatidylinositol bisphosphate (PIP2) into the two second messengers inositol 1,4,5-trisphosphate (InsP3) and diacyl glycerol (DAG). DAG activates protein kinase C, and InsP3 mobilizes calcium from intracellular stores via the InsP3 receptor. Changes in [Ca2+]i regulate the function of a wide range of target proteins, including ion channels, kinases, phosphatases, proteases, and transcription factors (Berridge, 1993). In the mouse, there are three InsP3R genes, and type 1 InsP3R mutants display ataxia and epileptic seizures (Matsumoto et al., 1996). In Drosophila, only one InsP3 receptor (InsP3R) gene is known, and it is expressed ubiquitously throughout development (Hasan and Rosbash, 1992; Yoshikawa et al., 1992; Raghu and Hasan, 1995). Here, we characterize Drosophila InsP3R mutants and demonstrate that the InsP3R is essential for embryonic and larval development. Interestingly, maternal InsP3R mRNA is sufficient for progression through the embryonic stages, but larval organs show asynchronous and defective cell divisions, and imaginal discs arrest early and fail to differentiate. We also generated adult mosaic animals and demonstrate that phototransduction, a model PLC pathway thought to require InsP3R, does not require InsP3R for signaling.

Animals↗

Studies of acetylcholine receptor subunit gene expression: chromatin structural changes during myogenesis.

Myogenesis proceeds stepwise from pluripotential stem cell to differentiated myotube. The precise number of transitions that occur along the developmental pathway remains to be determined. We examined the myogenic pathway as modelled by mouse mesodermal stem cell and muscle cell lines for stage-specific alterations in the chromatin structure of the acetylcholine receptor delta and gamma subunit genes. We reasoned that such an analysis would allow us to observe either the primary events in the activation of these muscle-specific genes or processes secondary to the binding of muscle-specific regulatory proteins. We probed chromatin structure with DNase I (deoxyribonuclease I) and precisely mapped to the 5' ends of the delta and gamma genes DNase I hypersensitive (DH) sites whose induction is unique to each myogenic stage. Putative mesodermal stem cells have the simplest pattern of DH sites with no sites near the 5' ends of the delta and gamma genes, whereas differentiated myotubes express the most complex pattern; the myoblast pattern is intermediate and of two types. In muscle cell lines where differentiation must be induced the myoblasts have a simple pattern (one more site than stem cells); in muscle lines where differentiation is spontaneous the myoblasts express a complex pattern of DH sites (one less site than myotubes). Inducible myoblasts seem to be arrested in an earlier step in the myogenic pathway than spontaneously differentiating myoblasts. Thus, myogenic activation of acetylcholine receptor subunit genes appears to be a stepwise process that can be detected by chromatin structural changes specific to four distinct stages of muscle development: stem cell, early myoblast, late myoblast, and differentiated myotube.

Animals↗

A biologically based dynamic model for predicting the disposition of methanol and its metabolites in animals and humans.

A multicompartment biologically based dynamic model was developed to describe the time evolution of methanol and its metabolites in the whole body and in accessible biological matrices of rats, monkeys, and humans following different exposure scenarios. The dynamic of intercompartment exchanges was described mathematically by a mass balance differential equation system. The model's conceptual and functional representation was the same for rats, monkeys, and humans, but relevant published data specific to the species of interest served to determine the critical parameters of the kinetics. Simulations provided a close approximation to kinetic data available in the published literature. The average pulmonary absorption fraction of methanol was estimated to be 0.60 in rats, 0.69 in monkeys, and 0.58-0.82 in human volunteers. The corresponding average elimination half-life of absorbed methanol through metabolism to formaldehyde was estimated to be 1.3, 0.7-3.2, and 1.7 h. Saturation of methanol metabolism appeared to occur at a lower exposure in rats than in monkeys and humans. Also, the main species difference in the kinetics was attributed to a metabolism rate constant of whole body formaldehyde to formate estimated to be twice as high in rats as in monkeys. Inversely, in monkeys and in humans, a larger fraction of body burden of formaldehyde is rapidly transferred to a long-term component. The latter represents the formaldehyde that (directly or after oxidation to formate) binds to various endogenous molecules or is taken up by the tetrahydrofolic-acid-dependent one-carbon pathway to become the building block of synthetic pathways. This model can be used to quantitatively relate methanol or its metabolites in biological matrices to the absorbed dose and tissue burden at any point in time in rats, monkeys, and humans for different exposures, thus reducing uncertainties in the dose-response relationship, and animal-to-human and exposure scenario comparisons. The model, adapted to kinetic data in human volunteers exposed acutely to methanol vapors, predicts that 8-h inhalation exposures ranging from 500 to 2000 ppm, without physical activities, are needed to increase concentrations of blood formate and urinary formic acid above mean background values reported by various authors (4.9-10.3 and 6.3-13 mg/liter, respectively). This leaves blood and urinary methanol concentrations as the most sensitive biomarkers of absorbed methanol.

Air Pollutants, Occupational↗

Severe hypercholesterolemia, hypertriglyceridemia, and atherosclerosis in mice lacking both leptin and the low density lipoprotein receptor.

Leptin-deficient mice (ob/ob) are an excellent murine model for obesity, insulin resistance, and diabetes, all of which are components of a multiple risk factor syndrome that, along with hypercholesterolemia, precipitates a potential high risk for atherosclerosis. In the current study, we show an unexpectedly severe hyperlipidemia in ob/ob mice on a background of low density lipoprotein receptor (LDLR) deficiency (-/-). Doubly mutant mice (LDLR-/-;ob/ob) exhibited striking elevations in both total plasma cholesterol (TC) and triglyceride (TG) levels (1715 +/- 87 and 1016 +/- 172 mg/dl, respectively), at age 3-4 months, resulting in extensive atherosclerotic lesions throughout the aorta by 6 months. Lipoprotein analyses revealed the elevated TC and TG levels to be due to a large increase in an apoB-containing broad-beta remnant lipoprotein fraction. While fasting, diet restriction, and low level leptin treatment significantly lowered TG levels, they caused only slight changes in TC levels. Hepatic cholesterol and triglyceride contents as well as mRNA levels of cholesterologenic and lipogenic enzymes suggest that leptin deficiency increased hepatic triglyceride production but did not change cholesterol production in ob/ob mice regardless of their LDLR genotype. These data provide evidence that the hypertriglyceridemia and hypercholesterolemia in the doubly mutant mice are caused by distinct mechanisms and point to the possibility that leptin might have some impact on plasma cholesterol metabolism, possibly through an LDLR-independent pathway. This model will be an excellent tool for future studies on the relationship between impaired fuel metabolism, increased plasma remnant lipoproteins, diabetes, and atherosclerosis.

Animals↗

A model for the selective loss of major histocompatibility complex self-restricted T cell immune responses during the development of acquired immune deficiency syndrome (AIDS).

Functional analyses of peripheral blood leukocytes (PBL) from high risk HTLV-III antibody-negative and antibody-positive donors, as well as from patients with acquired immune deficiency syndrome (AIDS), lymphadenopathy syndrome (LAS), and AIDS-related complex (ARC) were performed by the in vitro generation of cytotoxic T lymphocyte (CTL) and proliferative responses to the HLA self-restricted antigens of influenza virus (S + X) and to nonself restricted HLA alloantigens (ALLO). All 40 antibody-negative donors tested responded to both S + X and ALLO in the CTL response, whereas six of 14 antibody-positive, two of three LAS, four of five ARC, and seven of 17 AIDS patients exhibited a selective absence of CTL to S + X, but generated normal or elevated CTL responses to ALLO. Of the remaining 10 AIDS patients, nine did not respond to either S + X or ALLO, and one responded to both S + X and ALLO. A similar selective loss of the proliferative response to S + X was found. We also observed antibody-positive donors who initially generated CTL responses to S + X and ALLO, but lost the S + X response as a function of time. We were able to restore the selective loss of S + X CTL activity in vitro by the addition of IL 2 and, to some extent, by co-stimulation with S + X plus ALLO. Depletion of CD4+ T helper cells and removal of autologous antigen-presenting cells from the PBL of healthy antibody-negative donors indicated that distinct T helper cell subsets exist in human PBL, and that S + X responses must use a CD4+ T helper population, whereas ALLO responses can utilize an alternate CD4- T helper pathway. A model is presented indicating the selective depletion of CD4+ T helper function in the developmental stages of AIDS. The functional test for T helper activity to self restricted antigens may be the earliest indicator of immune functional loss in the development of AIDS, and may precede a reduction in the absolute number of CD4+ cells.

AIDS-Related Complex↗

An integrated modeling-experimental strategy for the analysis of metabolic pathways.

An inherent problem in studying the behavior of a metabolic pathway is the impossibility of developing a complete, detailed model that includes all the cellular processes that have an impact on the set of fluxes in such a pathway. Lacking this, one requires some means of modeling the interactions between a metabolic pathway and other cellular processes for the purpose of analyzing pathway characteristics within the cell (e.g., determining sensitivity coefficients for various steps in the pathway) with a minimal amount of time and effort. A general framework is developed for studying these issues in a rigorous manner. Using this framework, detailed knowledge about a metabolic pathway (i.e., a set of rate expressions for steps in the pathway) can be combined with the results from a relatively simple set of experiments in order to obtain estimates for the sensitivity of the pathway to enzyme activities, inhibition constants, and other parameters that determine the pathway's behavior, while accounting for the pathway's interaction with the rest of the cellular metabolism. A model system representing amino acid production is used to illustrate the problem and to provide results based on computational experiments. The modeling strategy described here should be useful in genetic design to improve pathway fluxes and metabolic network selectivity.

Amino Acids↗

Allelic loss of p53 gene is associated with genesis and maintenance, but not invasion, of mouse carcinoma in situ of the bladder.

Carcinoma in situ (CIS) of the bladder has recently been proposed to be a heterogeneous group of diseases with varied histogenesis and biological behavior. In this study, we describe the sequential steps of CIS development and progression in a transgenic mouse model expressing low levels of the SV40 large T antigen. We found that CIS in transgenic mice arose from urothelial dysplasia, that CIS could persist for an extended period of time without invasion, and that the majority of CIS eventually evolved into high-grade, superficial, papillary tumors before a small fraction of them advanced to invasion/metastasis. A genome-wide search of chromosomal imbalances by comparative genomic hybridization revealed that 9 of 11 (82%) of CIS had losses on chromosome 11. Southern blotting demonstrated the allelic loss of the p53 gene, which resides on mouse chromosome 11, in four comparative genomic hybridization-tested tumors and 10 of 11 (91%) additional CIS examined. Consistent with the reduced p53 gene dosage because of the allelic loss and the functional inactivation of p53 protein of the remaining allele by SV40T antigen, there was a dramatic decrease in CIS of Mdm-2, a major p53 target. In contrast, the level of p21, another p53 target, was largely unaltered, suggesting that p21 expression can be regulated by p53-independent mechanisms. These results delineate the early stages of bladder tumorigenesis and suggest that the loss of a p53-bearing chromosome is an early event in bladder tumorigenesis and is crucial for the genesis and the maintenance, but not the progression, of bladder CIS. On the basis of our current and previous transgenic studies, we have proposed an integrated pathway progression model of bladder cancer.

Alleles↗

Impact of correlated synaptic input on output firing rate and variability in simple neuronal models.

Cortical neurons are typically driven by thousands of synaptic inputs. The arrival of a spike from one input may or may not be correlated with the arrival of other spikes from different inputs. How does this interdependence alter the probability that the postsynaptic neuron will fire? We constructed a simple random walk model in which the membrane potential of a target neuron fluctuates stochastically, driven by excitatory and inhibitory spikes arriving at random times. An analytic expression was derived for the mean output firing rate as a function of the firing rates and pairwise correlations of the inputs. This stochastic model made three quantitative predictions. (1) Correlations between pairs of excitatory or inhibitory inputs increase the fluctuations in synaptic drive, whereas correlations between excitatory-inhibitory pairs decrease them. (2) When excitation and inhibition are fully balanced (the mean net synaptic drive is zero), firing is caused by the fluctuations only. (3) In the balanced case, firing is irregular. These theoretical predictions were in excellent agreement with simulations of an integrate-and-fire neuron that included multiple conductances and received hundreds of synaptic inputs. The results show that, in the balanced regime, weak correlations caused by signals shared among inputs may have a multiplicative effect on the input-output rate curve of a postsynaptic neuron, i.e. they may regulate its gain; in the unbalanced regime, correlations may increase firing probability mainly around threshold, when output rate is low; and in all cases correlations are expected to increase the variability of the output spike train.

Action Potentials↗

Coupling of transcription termination to RNAi.

In metazoans, the mechanisms of transcriptional termination by RNA polymerase II (Pol II) and accelerated decay of messenger RNA (mRNA) following transcription shutdown are linked by sharing the same sequence elements and mRNA elongation, processing and termination factors. This begs the question, how could one process have two opposite outcomes, making or degrading mRNA? An integrated "allosteric-GENEi-torpedo" model that could explain this paradox predicts participation of two novel factors: (1) An allosteric factor, regulated by a physiological repressor, binds to a unique sequence element of a gene near the site of cleavage and polyadenylation, poly(A) site, and acts on the homologous site on the nascent transcript to cause its cleavage. The conformational changes of this factor determine the fate of nascent RNA, either to get cleaved and processed to mature mRNA for directing protein synthesis, or not to get cleaved and become template for double-stranded (ds) RNA synthesis. (2) A general transcription termination factor, recruited by transcribing Pol II at the poly(A) site, allostrically alters and induces Pol II to switch template from DNA to nascent RNA several hundred nucleotides downstream of the poly(A) site. The template switch disengages Pol II from DNA and effectively terminates transcription. The Pol II with newly acquired RNA-dependent RNA polymerase activity retraces its path, back along the nascent RNA, so generating dsRNA. The extent to which it can retrace this path is determined by the factors influencing the cleavage of the pre-mRNA at the site of polyA addition. If cleavage and polyadenylation occur, the retracing is cut short, the 3' RNA is degraded by an exonuclease and the polymerase is liberated to reinitiate transcription. If the cleavage is inhibited, then a full-length dsRNA can be produced. This can then be subject to cleavage by "Dicer", which generates fragments of approximately 22bp that guide degradation of the cognate mRNA via the RNA interference (RNAi) pathway. This model complements the current "allosteric-torpedo" model of transcription termination, and could explain the apparent paradox of the divergent results of a common biological process.

DEAD-box RNA Helicases↗

Estimating resting motor thresholds in transcranial magnetic stimulation research and practice: a computer simulation evaluation of best methods.

OBJECTIVES: Resting motor threshold is the basic unit of dosing in transcranial magnetic stimulation (TMS) research and practice. There is little consensus on how best to estimate resting motor threshold with TMS, and only a few tools and resources are readily available to TMS researchers. The current study investigates the accuracy and efficiency of 5 different approaches to motor threshold assessment for TMS research and practice applications. METHODS: Computer simulation models are used to test the efficiency and accuracy of 5 different adaptive parameter estimation by sequential testing (PEST) procedures. For each approach, data are presented with respect to the mean number of TMS trials necessary to reach the motor threshold estimate as well as the mean accuracy of the estimates. RESULTS: A simple nonparametric PEST procedure appears to provide the most accurate motor threshold estimates, but takes slightly longer (on average, 3.48 trials) to complete than a popular parametric alternative (maximum likelihood PEST). Recommendations are made for the best starting values for each of the approaches to maximize both efficiency and accuracy. CONCLUSIONS: In light of the computer simulation data provided in this article, the authors review and suggest which techniques might best fit different TMS research and clinical situations. Lastly, a free user-friendly software package is described and made available on the world wide web that allows users to run all of the motor threshold estimation procedures discussed in this article for clinical and research applications.

Computer Simulation↗

Estimation of metabolic flux rates in liver purine catabolism of tumour-bearing mice by computer simulation of radioactive tracer experiments.

Mouse hepatocytes from healthy control mice and from Ehrlich ascites tumour-bearing mice were used for tracer-kinetic studies of purine catabolism of liver cells during different periods of tumour growth. The dynamics of the radioactive tracers were modelled mathematically by a system of differential equations. Computer simulations, i.e. direct fitting of numerical solutions of these equations to the observed time-courses of metabolites and specific radioactivities, enables one to estimate unknown kinetic parameters of a simplified model of pathways of hepatic purine catabolism in tumour-bearing mice. There occurred great differences of metabolic flux rates between control hepatocytes, hepatocytes of mice during the proliferating period of tumour growth (6th day after inoculation of the tumour) and hepatocytes of mice during the resting period of tumour growth (12th day after inoculation of the tumour). The final purine degradation of hepatocytes prepared during the proliferating period was lower in comparison with that of control hepatocytes, but it was markedly higher in hepatocytes prepared during the resting period of tumour growth. The changes in hepatocyte purine catabolism during the proliferating period of tumour growth argue for transitions which aim at the maintenance of high purine nucleotide levels in the liver itself rather than for an increased nucleoside and nucleobase supply for the tumour. This suggestion is in accordance with the increased ATP level of the liver during the proliferating phase of tumour growth.(ABSTRACT TRUNCATED AT 250 WORDS)

Animals↗

Modeling the MHC class I pathway by combining predictions of proteasomal cleavage, TAP transport and MHC class I binding.

Epitopes presented by major histocompatibility complex (MHC) class I molecules are selected by a multi-step process. Here we present the first computational prediction of this process based on in vitro experiments characterizing proteasomal cleavage, transport by the transporter associated with antigen processing (TAP) and MHC class I binding. Our novel prediction method for proteasomal cleavages outperforms existing methods when tested on in vitro cleavage data. The analysis of our predictions for a new dataset consisting of 390 endogenously processed MHC class I ligands from cells with known proteasome composition shows that the immunological advantage of switching from constitutive to immunoproteasomes is mainly to suppress the creation of peptides in the cytosol that TAP cannot transport. Furthermore, we show that proteasomes are unlikely to generate MHC class I ligands with a C-terminal lysine residue, suggesting processing of these ligands by a different protease that may be tripeptidyl-peptidase II (TPPII).

ATP-Binding Cassette Transporters↗

Advanced computing for systems biology.

Systems biology is based on computational modelling and simulation of large networks of interacting components. Models may be intended to capture processes, mechanisms, components and interactions at different levels of fidelity. Input data are often large and geographically disperse, and may require the computation to be moved to the data, not vice versa. In addition, complex system-level problems require collaboration across institutions and disciplines. Grid computing can offer robust, scaleable solutions for distributed data, compute and expertise. We illustrate some of the range of computational and data requirements in systems biology with three case studies: one requiring large computation but small data (orthologue mapping in comparative genomics), a second involving complex terabyte data (the Visible Cell project) and a third that is both computationally and data-intensive (simulations at multiple temporal and spatial scales). Authentication, authorisation and audit systems are currently not well scalable and may present bottlenecks for distributed collaboration particularly where outcomes may be commercialised. Challenges remain in providing lightweight standards to facilitate the penetration of robust, scalable grid-type computing into diverse user communities to meet the evolving demands of systems biology.

Animals↗

Neuroanatomical targets of anxiogenic drugs in the hindbrain as revealed by Fos immunocytochemistry.

It is speculated that specific hindbrain transmitter pathways centred on the periaqueductal gray and locus coeruleus are an important integrative neural substrate for the expression of anxiety and the somatic symptoms and cardiovascular changes that accompany severe anxiety states, such as in panic disorder. Here we investigated the effects of various drugs, known to induce panic in humans and to be anxiogenic in animals, on Fos expression in the periaqueductal gray, locus coeruleus and other parts of the rat hindbrain. The drugs tested were the benozodiazepine inverse agonist FG-7142, the alpha(2)-adrenoceptor antagonist yohimbine, the non-selective 5-hydroxytryptamine(2C) receptor agonist m-chlorophenyl piperazine, the adenosine antagonist caffeine and the cholecystokinin analogue BOC-CCK(4). A clear-cut finding was that administration of each anxiogenic drug caused a striking region-specific pattern of Fos expression within the hindbrain. In particular, the drugs commonly increased Fos-like immunoreactivity in the periaqueductal gray and locus coeruleus. Increased Fos expression in the periaqueductal gray was specific to the rostral dorsolateral and caudal ventrolateral regions. All the anxiogenic drugs also increased Fos-like immunoreactivity in the lateral parabrachial nucleus and nucleus of the solitary tract and all but one (BOC-CCK(4)) increased Fos in the dorsal raphe nucleus. Rats habituated to the test environment and injected with saline vehicle displayed little or no Fos-like immunoreactivity in the hindbrain areas investigated. In summary, each of the anxiogenic drugs tested (FG-7142, yohimbine, m-chlorophenyl piperazine, caffeine and BOC-CCK(4)) increased Fos expression in a restricted number of hindbrain regions, including the periaqueductal gray and locus coeruleus. Previous Fos studies have found that these same regions are activated by various fearful environmental stimuli. Therefore, a specific set of hindbrain circuits may be commonly involved in the processing of anxiety-related information evoked by pharmacological and environmental manipulation. The present findings also raise the possibility that measurement of the effect of anxiogenic drugs on Fos expression might be a useful way to model hindbrain pathways activated by anxiety and possibly panic.

Adrenergic alpha-Antagonists↗

Programmed cell death in mature erythrocytes: a model for investigating death effector pathways operating in the absence of mitochondria.

Human mature erythrocytes have been considered as unable to undergo programmed cell death (PCD), due to their lack of mitochondria, nucleus and other organelles, and to the finding that they survive two conditions that induce PCD in vitro in all human nucleated cells, treatment with staurosporine and serum deprivation. Here we report that mature erythrocytes can undergo a rapid self-destruction process sharing several features with apoptosis, including cell shrinkage, plasma membrane microvesiculation, phosphatidylserine externalization, and leading to erythrocyte disintegration, or, in the presence of macrophages, to macrophage ingestion of dying erythrocytes. This regulated form of PCD was induced by Ca(2+) influx, and prevented by cysteine protease inhibitors that allowed erythrocyte survival in vitro and in vivo. The cysteine proteinases involved seem not to be caspases, since (i) proforms of caspase 3, while present in erythrocytes, were not activated during erythrocyte death; (ii) cytochrome c, a critical component of the apoptosome, was lacking; and (iii) cell-free assays did not detect activated effectors of nuclear apoptosis in dying erythrocytes. Our findings provide the first identification that a death program can operate in the absence of mitochondria. They indicate that mature erythrocytes share with all other mammalian cell types the capacity to self-destruct in response to environmental signals, and imply that erythrocyte survival may be modulated by therapeutic intervention.

Animals↗

Consensus by democracy. Using meta-analyses of microarray and genomic data to model the cold acclimation signaling pathway in Arabidopsis.

The whole-genome response of Arabidopsis (Arabidopsis thaliana) exposed to different types and durations of abiotic stress has now been described by a wealth of publicly available microarray data. When combined with studies of how gene expression is affected in mutant and transgenic Arabidopsis with altered ability to transduce the low temperature signal, these data can be used to test the interactions between various low temperature-associated transcription factors and their regulons. We quantized a collection of Affymetrix microarray data so that each gene in a particular regulon could vote on whether a cis-element found in its promoter conferred induction (+1), repression (-1), or no transcriptional change (0) during cold stress. By statistically comparing these election results with the voting behavior of all genes on the same gene chip, we verified the bioactivity of novel cis-elements and defined whether they were inductive or repressive. Using in silico mutagenesis we identified functional binding consensus variants for the transcription factors studied. Our results suggest that the previously identified ICEr1 (induction of CBF expression region 1) consensus does not correlate with cold gene induction, while the ICEr3/ICEr4 consensuses identified using our algorithms are present in regulons of genes that were induced coordinate with observed ICE1 transcript accumulation and temporally preceding genes containing the dehydration response element. Statistical analysis of overlap and cis-element enrichment in the ICE1, CBF2, ZAT12, HOS9, and PHYA regulons enabled us to construct a regulatory network supported by multiple lines of evidence that can be used for future hypothesis testing.

Acclimatization↗