Studies of the Human Transcriptome

Gene expression is one of the most critical factors influencing the phenotype of a cell. As a result of several technological advances, measuring gene expression levels has become one of the most common molecular biological measurements to study the behaviour of cells. The scientific community has produced enormous and constantly increasing collection of gene expression data from various human cells both from healthy and pathological conditions. However, while each of these studies is informative and enlighting in its own context and research setup, diverging methods and terminologies make it very challenging to integrate existing gene expression data to a more comprehensive view of human transcriptome function. On the other hand, bioinformatic science advances only through data integration and synthesis. The aim of this study was to develop biological and mathematical methods to overcome these challenges and to construct an integrated database of human transcriptome as well as to demonstrate its usage. Methods developed in this study can be divided in two distinct parts. First, the biological and medical annotation of the existing gene expression measurements needed to be encoded by systematic vocabularies. There was no single existing biomedical ontology or vocabulary suitable for this purpose. Thus, new annotation terminology was developed as a part of this work. Second part was to develop mathematical methods correcting the noise and systematic differences/errors in the data caused by various array generations. Additionally, there was a need to develop suitable computational methods for sample collection and archiving, unique sample identification, database structures, data retrieval and visualization. Bioinformatic methods were developed to analyze gene expression levels and putative functional associations of human genes by using the integrated gene expression data. Also a method to interpret individual gene expression profiles across all the healthy and pathological tissues of the reference database was developed. As a result of this work 9783 human gene expression samples measured by Affymetrix microarrays were integrated to form a unique human transcriptome resource GeneSapiens. This makes it possible to analyse expression levels of 17330 genes across 175 types of healthy and pathological human tissues. Application of this resource to interpret individual gene expression measurements allowed identification of tissue of origin with 92.0% accuracy among 44 healthy tissue types. Systematic analysis of transcriptional activity levels of 459 kinase genes was performed across 44 healthy and 55 pathological tissue types and a genome wide analysis of kinase gene co-expression networks was done. This analysis revealed biologically and medically interesting data on putative kinase gene functions in health and disease. Finally, we developed a method for alignment of gene expression profiles (AGEP) to perform analysis for individual patient samples to pinpoint gene- and pathway-specific changes in the test sample in relation to the reference transcriptome database. We also showed how large-scale gene expression data resources can be used to quantitatively characterize changes in the transcriptomic program of differentiating stem cells. Taken together, these studies indicate the power of systematic bioinformatic analyses to infer biological and medical insights from existing published datasets as well as to facilitate the interpretation of new molecular profiling data from individual patients.

Diel variation in fig volatiles across syconium development: making sense of scents

Plants produce volatile organic compounds (VOCs) in a variety of contexts that include response to abiotic and biotic stresses, attraction of pollinators and parasitoids, and repulsion of herbivores. Some of these VOCs may also exhibit diel variation in emission. In Ficus racemosa, we examined variation in VOCs released by fig syconia throughout syconium development and between day and night. Syconia are globular enclosed inflorescences that serve as developing nurseries for pollinating and parasitic fig wasps. Syconia are attacked by gallers early in their development, serviced by pollinators in mid phase, and are attractive to parasitoids in response to the development of gallers at later stages. VOC bouquets of the different development phases of the syconium were distinctive, as were their day and night VOC profiles. VOCs such as alpha-muurolene were characteristic of the pollen-receptive diurnal phase, and may serve to attract the diurnally-active pollinating wasps. Diel patterns of release of volatiles could not be correlated with their predicted volatility as determined by Henry's law constants at ambient temperatures. Therefore, factors other than Henry's law constant such as stomatal conductance or VOC synthesis must explain diel variation in VOC emission. A novel use of weighted gene co-expression network analysis (WGCNA) on the volatilome resulted in seven distinct modules of co-emitted VOCs that could be interpreted on the basis of syconium ecology. Some modules were characterized by the response of fig syconia to early galling by parasitic wasps and consisted largely of green leaf volatiles (GLVs). Other modules, that could be characterized by a combination of syconia response to oviposition and tissue feeding by larvae of herbivorous galler pollinators as well as of parasitized wasps, consisted largely of putative herbivore-induced plant volatiles (HIPVs). We demonstrated the usefulness of WGCNA analysis of the volatilome in making sense of the scents produced by the syconia at different stages and diel phases of their development.

A naturally occurring amino acid substitution in the voltage-dependent sodium channel selectivity filter affects channel gating

The pore of sodium channels contains a selectivity filter made of 4 amino acids, D/E/K/A. In voltage sensitive sodium channel (Nav) channels from jellyfish to human the fourth amino acid is Ala. This Ala, when mutated to Asp, promotes slow inactivation. In some Nav channels of pufferfishes, the Ala is replaced with Gly. We studied the biophysical properties of an Ala-to-Gly substitution (A1529G) in rat Nav1.4 channel expressed in Xenopus oocytes alone or with a beta 1 subunit. The Ala-to-Gly substitution does not affect monovalent cation selectivity and positively shifts the voltage-dependent inactivation curve, although co-expression with a beta 1 subunit eliminates the difference between A1529G and WT. There is almost no difference in channel fast inactivation, but the beta 1 subunit accelerates WT current inactivation significantly more than it does the A1529G channels. The Ala-to-Gly substitution mainly influences the rate of recovery from slow inactivation. Again, the beta 1 subunit is less effective on speeding recovery of A1529G than the WT. We searched Nav channels in numerous databases and noted at least four other independent Ala-to-Gly substitutions in Nav channels in teleost fishes. Thus, the Ala-to-Gly substitution occurs more frequently than previously realized, possibly under selection for alterations of channel gating.

GBNV encoded movement protein (NSm) remodels ER network via C-terminal coiled coil domain

Plant viruses exploit the host machinery for targeting the viral genome-movement protein complex to plasmodesmata (PD). The mechanism by which the non-structural protein m (NSm) of Groundnut bud necrosis virus (GBNV) is targeted to PD was investigated using Agrobacterium mediated transient expression of NSm and its fusion proteins in Nicotiana benthamiana. GFP:NSm formed punctuate structures that colocalized with mCherry:plasmodesmata localized protein la (PDLP la) confirming that GBNV NSm localizes to PD. Unlike in other movement proteins, the C-terminal coiled coil domain of GBNV NSm was shown to be involved in the localization of NSm to PD, as deletion of this domain resulted in the cytoplasmic localization of NSm. Treatment with Brefeldin A demonstrated the role of ER in targeting GFP NSm to PD. Furthermore, mCherry:NSm co-localized with ER-GFP (endoplasmic reticulum targeting peptide (HDEL peptide fused with GFP). Co-expression of NSm with ER-GFP showed that the ER-network was transformed into vesicles indicating that NSm interacts with ER and remodels it. Mutations in the conserved hydrophobic region of NSm (residues 130-138) did not abolish the formation of vesicles. Additionally, the conserved prolines at positions 140 and 142 were found to be essential for targeting the vesicles to the cell membrane. Further, systematic deletion of amino acid residues from N- and C-terminus demonstrated that N-terminal 203 amino acids are dispensable for the vesicle formation. On the other hand, the C-terminal coiled coil domain when expressed alone could also form vesicles. These results suggest that GBNV NSm remodels the ER network by forming vesicles via its interaction through the C-terminal coiled coil domain. Interestingly, NSm interacts with NP in vitro and coexpression of these two proteins in planta resulted in the relocalization of NP to PD and this relocalization was abolished when the N-terminal unfolded region of NSm was deleted. Thus, the NSm interacts with NP via its N-terminal unfolded region and the NSm-NP complex could in turn interact with the ER membrane via the C-terminal coiled coil domain of NSm to form vesicles that are targeted to PD and there by assist the cell to cell movement of the viral genome complex. (C) 2015 Elsevier Inc. All rights reserved.

Estudo fenotípico da região Aorta-Gônada-Mesonefros, de embriões de galinhas Gallus gallus domesticus, com ênfase na aorta dorsal

Acredita-se que os primeiros progenitores da hematopoese definitiva surjam da diferenciação do endotélio da aorta dorsal, na altura da região da Aorta-Gônada-Mesonefros (AGM). Com o intuito de estudar esta região e o fenótipo das células do endotélio da aorta dorsal nesta posição topográfica, ovos galados de Gallus gallus domesticus L. foram incubados em chocadeira, classificados em estádios de E16 a E25 e processados histotecnologicamente para obtenção de secções seriadas na altura da região AGM. Estas passaram por coloração por Hematoxilina-Eosina, histoquímica para PAS, PAS-diastase e Alcian Blue pH 1.0 e pH 2.5, histoquímica por lectinas fluoresceinadas e imunofluorescência para moléculas de superfície, citoesqueleto e matriz extracelular. Foi observada hipertrofia endotelial no assoalho da aorta nos estádios observados, o qual se apresentava positivo ao PAS, com ocorrência frequente de vacuolizações basais PAS negativas, e o surgimento ocasional de grupamentos celulares intravasculares. Nestes, as células que se destacavam da membrana basal do endotélio expressavam progressivamente mais material PAS positivo, o qual, no entanto, em nenhum momento pareceu se tratar de glicogênio. Em relação às glicosaminoglicanas, notamos a presença predominante de ácido hialurônico por todo o mesênquima da região e em outras estruturas como periferia da notocorda, tubo neural e mesoderma lateral. Ocorreu co-expressão de fibronectina e α-actina de músculo liso em células circunjacentes à aorta, na face ventral do vaso. GFAP e BMP-4 são expressas entre as células do tubo neural e em sua periferia, assim como na notocorda do embrião. As lectinas Abrus precatorius, Lens culinarise Ricinus communis mostraram-se positivas principalmente na região subedotelial do assoalho da aorta nos estádios observados neste trabalho. Bandeiraea simplicifolia exibiu pouca marcação na aorta dorsal e a Arachis hypogeae foi negativa. Outras estruturas da região AGM também expressaram resíduos de açúcares revelados por estas lectinas, tais como: notocorda, tubo neural, mesênquima, intestino primitivo e saco vitelínico. Estes resultados acrescentam elementos morfológicos e bioquímicos ao conhecimento sobre a região AGM de embriões de galinha e sobre o endotélio, possivelmente hemogênico, da aorta dorsal.

Structure-function analysis of USP1: insights into the role of Ser313 phosphorylation site and the effect of cancer-associated mutations on autocleavage

Background: In complex with its cofactor UAF1, the USP1 deubiquitinase plays an important role in cellular processes related to cancer, including the response to DNA damage. The USP1/UAF1 complex is emerging as a novel target in cancer therapy, but several aspects of its function and regulation remain to be further clarified. These include the role of the serine 313 phosphorylation site, the relative contribution of different USP1 sequence motifs to UAF1 binding, and the potential effect of cancer-associated mutations on USP1 regulation by autocleavage. Methods: We have generated a large set of USP1 structural variants, including a catalytically inactive form (C90S), non-phosphorylatable (S313A) and phosphomimetic (S313D) mutants, deletion mutants lacking potential UAF1 binding sites, a mutant (GG/AA) unable to undergo autocleavage at the well-characterized G670/G671 diglycine motif, and four USP1 mutants identified in tumor samples that cluster around this cleavage site (G667A, L669P, K673T and A676T). Using cell-based assays, we have determined the ability of these mutants to bind UAF1, to reverse DNA damage-induced monoubiquitination of PCNA, and to undergo autocleavage. Results: A non-phosphorylatable S313A mutant of USP1 retained the ability to bind UAF1 and to reverse PCNA ubiquitination in cell-based assays. Regardless of the presence of a phosphomimetic S313D mutation, deletion of USP1 fragment 420-520 disrupted UAF1 binding, as determined using a nuclear relocation assay. The UAF1 binding site in a second UAF1-interacting DUB, USP46, was mapped to a region homologous to USP1(420-520). Regarding USP1 autocleavage, co-expression of the C90S and GG/AA mutants did not result in cleavage, while the cancer-associated mutation L669P was found to reduce cleavage efficiency. Conclusions: USP1 phosphorylation at S313 is not critical for PCNA deubiquitination, neither for binding to UAF1 in a cellular environment. In this context, USP1 amino acid motif 420-520 is necessary and sufficient for UAF1 binding. This motif, and a homologous amino acid segment that mediates USP46 binding to UAF1, map to the Fingers sub-domain of these DUBs. On the other hand, our results support the view that USP1 autocleavage may occur in cis, and can be altered by a cancer-associated mutation.

Discovering transcriptional modules by Bayesian data integration.

MOTIVATION: We present a method for directly inferring transcriptional modules (TMs) by integrating gene expression and transcription factor binding (ChIP-chip) data. Our model extends a hierarchical Dirichlet process mixture model to allow data fusion on a gene-by-gene basis. This encodes the intuition that co-expression and co-regulation are not necessarily equivalent and hence we do not expect all genes to group similarly in both datasets. In particular, it allows us to identify the subset of genes that share the same structure of transcriptional modules in both datasets. RESULTS: We find that by working on a gene-by-gene basis, our model is able to extract clusters with greater functional coherence than existing methods. By combining gene expression and transcription factor binding (ChIP-chip) data in this way, we are better able to determine the groups of genes that are most likely to represent underlying TMs. AVAILABILITY: If interested in the code for the work presented in this article, please contact the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNF220, an E3 ubiquitin ligase that targets Sin3B for ubiquitination

Modification of proteins by ubiquitination plays important roles in various cellular processes. During this process, the target specificity is determined by ubiquitin ligases. Here we identify RNF220 (RING finger protein 220) as a novel ubiquitin ligase for Sin3B. As a conserved RING protein, RNF220 can bind E2 and mediate auto-ubiquitination of itself. Through a yeast two-hybrid screen, we isolated Sin3B as one of its targets, which is a scaffold protein of the Sin3/HDAC (histone deacetylase) corepressor complex. RNF220 specifically interacts with Sin3B both in vitro and in vivo. Sin3B can be regulated by the ubiquitin-proteasome system. Co-expression of RNF220 promotes the ubiquitination and proteasomal degradation of Sin3B. Taken together, these results reveal a new mechanism for regulating the Sin3/HDAC complex. (C) 2010 Elsevier Inc. All rights reserved.

Characterization of TRAF6 mediated ubiquitination of presenilins and γ-secretase substrates

Post-translational modification of the γ-secretase protease complexes and their substrates has an important role in controlling receptor-initiated signalling events, which are critically important in the pathogenesis of cancer, inflammatory and Alzheimer’s disease. Our lab has previously characterised an interaction between TRAF6 and presenilin-1, which lead to the identification of interleukin-1 (IL-1) receptor type 1 (IL-1R1) and Toll-like receptor-4 (TLR4) as novel γ-secretase substrates. Subsequently our group showed that TRAF6 promoted ubiquitination and γ-secretase cleavage of IL-1R1. The aim of this project is to study the association between TRAF6 and the presenilins, the critical γ-secretase complex components, and to determine the functional importance of TRAF6-mediated ubiquitination of γ-secretase substrates. Firstly, we show that the full-length presenilins are novel substrates of TRAF6-mediated Lysine-63-linked polyubiquitination. Secondly, we show that co-expression of TRAF6 and the presenilins increases the stability and alters the turnover of the presenilins. Thirdly, we reveal that TRAF6-mediated ubiquitination of presenilin does not affect γ-secretase enzyme activity, but may regulate the full-length presenilin functions such as ER Ca2+ signalling. Previously, we have reported IL-1R1 as a novel substrate of TRAF6-mediated ubiquitination. In this study, we identified five lysine residues in the IL-1R1 intracellular domain targeted by TRAF6-mediated polyubiquitination. Furthermore, mutagenesis of these five lysine residues led to decreased IL-1R1 cell surface expression, precluded the ectodomain shedding and attenuated the responsiveness to IL-1β stimulation, demonstrating the critical role of TRAF6 in IL-1R1 trafficking.

Investigation of p75 neurotrophin receptor and the role of its post-translational modifications in tumour cell motility and invasiveness

The p75 neurotrophin receptor (p75NTR) is a member of the tumour necrosis factor superfamily, which relies on the recruitment of cytosolic protein partners - including the TNF receptor associated factor 6 (TRAF6) E3 ubiquitin ligase - to produce cellular responses such as apoptosis, survival, and inhibition of neurite outgrowth. Recently,p75NTR was also shown to undergo γ-secretase-mediated regulated intramembrane proteolysis, and the receptor ICD was found to migrate to the nucleus where it regulates gene transcription. Moreover, γ-secretase-mediated proteolysis was shown to be involved in glioblastoma cell migration and invasion. In this study we report that TRAF6-mediated K63-linked polyubiquitination at multiple or alternative lysine residues influences p75NTR-ICD stability in vitro. In addition, we found that TRAF6-mediated ubiquitination of p75NTR is not influenced by inhibition of dynamin. Moreover, we report beta-transducin repeats-containing protein (β-TrCP) as a novel E3- ligase that ubiquitinates p75NTR, which is independent of serine phosphorylation of the p75NTR destruction motif. In contrast to its influence on other substrates, co-expression of β-TrCP did not reduce p75NTR stability. We created U87-MG glioblastoma cell lines stably expressing wild type, γ-secretaseresistant and constitutively cleaved receptor, as well as the ICD-stabilized mutant K301R. Interestingly, only wild-type p75NTR induces increased glioblastoma cell migration, which could be reversed by application of γ-secretase inhibitor. Microarray and qRT-PCR analysis of mRNA transcripts in these cell lines yielded several promising genes that might be involved in glioblastoma cell migration and invasion, such as cadherin 11 and matrix metalloproteinase 12. Analysis of potential transcription factor binding sites revealed that transcription of these genes might be regulated by well known p75NTR signalling cascades such as NF-κB or JNK signalling, which are independent of γ-secretase-mediated cleavage of the receptor. In contrast, while p75NTR overexpression was confirmed in melanoma cell lines and a patient sample of melanoma metastasis to the brain, inhibition of γ-secretase did not influence melanoma cell migration. Collectively, this study provides several avenues to better understand the physiological importance of posttranslational modifications of p75NTR and the significance of the receptor in glioblastoma cell migration and invasion.

A sensitive bioassay for the mycotoxin aflatoxin B1, which also responds to the mycotoxins aflatoxin G1 and T-2 toxin, using engineered baker's yeast

A novel aflatoxin B(1) bioassay was created by introducing a Lipomyces kononenkoae alpha-amylase gene into a strain of S. cerevisiae capable of expressing the human cytochrome P450 3A4 (CYP3A4), and the cognate human CYP450 reductase. This strain and a dextranase-expressing strain were used in the development of a microtitre plate mycotoxin bioassay, which employed methanol as the solvent and polymyxin B nonapeptide as a permeation enhancer. Stable co-expression of the CYP3A4 gene system and of the dextranase and amylase genes in the two bioassay strains was demonstrated. The bioassay signalled toxicity as inhibition of secreted carbohydrase activity, using sensitive fluorimetric assays. The amylase-expressing strain could detect aflatoxin B(1) at 2 ng/ml, and was more sensitive than the dextranase-expressing strain. Aflatoxin G(1) could be detected at 2 microg/ml, and the trichothecene mycotoxin T-2 toxin was detectable at 100 ng/ml.

Kv1.5 is a major component underlying the A-type potassium current in retinal arteriolar smooth muscle.

Little is known about the molecular characteristics of the voltage-activated K(+) (K(v)) channels that underlie the A-type K(+) current in vascular smooth muscle cells of the systemic circulation. We investigated the molecular identity of the A-type K(+) current in retinal arteriolar myocytes using patch-clamp techniques, RT-PCR, immunohistochemistry, and neutralizing antibody studies. The A-type K(+) current was resistant to the actions of specific inhibitors for K(v)3 and K(v)4 channels but was blocked by the K(v)1 antagonist correolide. No effects were observed with pharmacological agents against K(v)1.1/2/3/6 and 7 channels, but the current was partially blocked by riluzole, a K(v)1.4 and K(v)1.5 inhibitor. The current was not altered by the removal of extracellular K(+) but was abolished by flecainide, indicative of K(v)1.5 rather than K(v)1.4 channels. Transcripts encoding K(v)1.5 and not K(v)1.4 were identified in freshly isolated retinal arterioles. Immunofluorescence labeling confirmed a lack of K(v)1.4 expression and revealed K(v)1.5 to be localized to the plasma membrane of the arteriolar smooth muscle cells. Anti-K(v)1.5 antibody applied intracellularly inhibited the A-type K(+) current, whereas anti-K(v)1.4 antibody had no effect. Co-expression of K(v)1.5 with K(v)beta1 or K(v)beta3 accessory subunits is known to transform K(v)1.5 currents from delayed rectifers into A-type currents. K(v)beta1 mRNA expression was detected in retinal arterioles, but K(v)beta3 was not observed. K(v)beta1 immunofluorescence was detected on the plasma membrane of retinal arteriolar myocytes. The findings of this study suggest that K(v)1.5, most likely co-assembled with K(v)beta1 subunits, comprises a major component underlying the A-type K(+) current in retinal arteriolar smooth muscle cells

Regulation of 3-phosphoinositide-dependent protein kinase-1 (PDK1) by src involves tyrosine phosphorylation of PDK1 and Src homology 2 domain binding

3-Phosphoinositide-dependent protein kinase-1 (PDK1) appears to play a central regulatory role in many cell signalings between phosphoinositide-3 kinase and various intracellular serine/threonine kinases. In resting cells, PDK1 is known to be constitutively active and is further activated by tyrosine phosphorylation (Tyr(9) and Tyr(373/376)) following the treatment of the cell with insulin or pervanadate. However, little is known about the mechanisms for this additional activation of PDK1. Here, we report that the SH2 domain of Src, Crk, and GAP recognized tyrosine-phosphorylated PDK1 in vitro. Destabilization of PDK1 induced by geldanamycin (a Hsp90 inhibitor) was partially blocked in HEK 293 cells expressing PDK1- Y9F. Co-expression of Hsp90 enhanced PDK1-Src complex formation and led to further increased PDK1 activity toward PKB and SGK. Immunohistochemical analysis with anti- phospho-Tyr9 antibodies showed that the level of Tyr9 phosphorylation was markedly increased in tumor samples compared with normal. Taken together, these data suggest that phosphorylation of PDK1 on Tyr9, distinct from Tyr373/376, is important for PDK1/Src complex formation, leading to PDK1 activation. Furthermore, Tyr9 phosphorylation is critical for the stabilization of both PDK1 and the PDK1/Src complex via Hsp90-mediated protection of PDK1 degradation.

Identification of tyrosine phosphorylation sites on 3-phosphoinositide-dependent protein kinase-1 and their role in regulating kinase activity

3-Phosphoinositide-dependent protein kinase-1 (PDK1) plays a central role in signal transduction pathways that activate phosphoinositide 3-kinase. Despite its key role as an upstream activator of enzymes such as protein kinase B and p70 ribosomal protein S6 kinase, the regulatory mechanisms controlling PDK1 activity are poorly understood. PDK1 has been reported to be constitutively active in resting cells and not further activated by growth factor stimulation (Casamayor, A., Morrice, N. A., and Alessi, D. R. (1999) Biochem. J. 342, 287-292). Here, we report that PDK1 becomes tyrosine-phosphorylated and translocates to the plasma membrane in response to pervanadate and insulin. Following pervanadate treatment, PDK1 kinase activity increased 1.5- to 3-fold whereas the activity of PDK1 associated with the plasma membrane increased similar to6-fold. The activity of PDK1 localized to the plasma membrane was also increased by insulin treatment. Three tyrosine phosphorylation sites of PDK1 (Tyr-9 and Tyr-373/376) were identified using in vivo labeling and mass spectrometry. Using site-directed mutants, we show that, although phosphorylation on Tyr-373/376 is important for PDK1 activity, phosphorylation on Tyr-9 has no effect on the activity of the kinase. Both of these residues can be phosphorylated by v-Src tyrosine kinase in vitro, and co-expression of v-Src leads to tyrosine phosphorylation and activation of PDK1. Thus, these data suggest that PDK1 activity is regulated by reversible phosphorylation, possibly by a member of the Src kinase family.

The GalF protein of Escherichia coli is not a UDP-glucose pyrophosphorylase but interacts with the GalU protein possibly to regulate cellular levels of UDP-glucose

We report the functional characterization of the galF gene of strain VW187 (Escherichia coli O7:K1), which encodes a polypeptide displaying structural features common to bacterial UDP-glucose pyrophosphorylases, including the E. coli GalU protein. These enzymes catalyse a reversible reaction converting UTP and glucose-1-phosphate into UDP-glucose and PPi. We show that, although the GalF protein is expressed in vivo, GalF-expressing plasmids cannot complement the phenotype of a galU mutant and extracts from this mutant which only produces GalF are enzymatically inactive. In contrast, the presence of GalU and GalF proteins in the same cell-free extract caused a significant reduction in the rate of pyrophosphorolysis (conversion of UDP-glucose into glucose-1-phosphate) but no significant effect on the kinetics of synthesis of UDP-glucose. The presence of GalF also increased the thermal stability of the enzyme in vitro. The effect of GalF in the biochemical properties of the UDP-glucose pyrophosphorylase required the co-synthesis of GalF and GalU, suggesting that they could interact as components of the oligomeric enzyme. The physical interaction of GalU and GalF was demonstrated in vivo by the co-expression of both proteins as fusion products using a yeast two-hybrid system. Furthermore, using a pair of galF-/galU+ and galF/galU+ isogenic strains, we demonstrated that the presence of GalF is associated with an increased concentration of intracellular UDP-glucose as well as with an enhancement of the thermal stability of the UDP-glucose pyrophosphorylase in vivo. We propose that GalF is a non-catalytic subunit of the UDP-glucose pyrophosphorylase modulating the enzyme activity to increase the formation of UDP-glucose, and this function is important for bacterial adaptation to conditions of stress.