Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Light-regulated interactions with SPA proteins underlie cryptochrome-mediated gene expression.

Cryptochromes are a class of photosensory receptors that control important processes in animals and plants primarily by regulating gene expression. How photon absorption by cryptochromes leads to changes in gene expression has remained largely elusive. Three recent studies, including Lian and colleagues (pp. 1023-1028) and Liu and colleagues (pp. 1029-1034) in this issue of Genes & Development, demonstrate that the interaction of light-activated Arabidopsis cryptochromes with a class of regulatory components of E3 ubiquitin ligase complexes leads to environmentally controlled abundance of transcriptional regulators.

Transcriptional repression of peroxisome proliferator-activated receptor beta/delta in murine keratinocytes by CCAAT/enhancer-binding proteins.

The roles of peroxisome proliferator-activated receptors (PPARs) and CCAAT/enhancer-binding proteins (C/EBPs) in keratinocyte and sebocyte differentiation suggest that both families of transcription factors closely interact in the skin. Initial characterization of the mouse PPARbeta promoter revealed an AP-1 site that is crucial for the regulation of PPARbeta expression in response to inflammatory cytokines in the skin. We now present evidence for a novel regulatory mechanism of the expression of the PPARbeta gene by which two members of the C/EBP family of transcription factors inhibit its basal promoter activity in mouse keratinocytes. We first demonstrate that C/EBPalpha and C/EBPbeta, but not C/EBPdelta, inhibit the expression of PPARbeta through the recruitment of a transcriptional repressor complex containing HDAC-1 to a specific C/EBP binding site on the PPARbeta promoter. Consistent with this repression, the expression patterns of PPARbeta and C/EBPs are mutually exclusive in keratinocytes of the interfollicular epidermis and hair follicles in mouse developing skin. This work reveals the importance of the regulatory interplay between PPARbeta and C/EBP transcription factors in the control of proliferation and differentiation in this organ. Such insights are crucial for the understanding of the molecular control regulating the balance between proliferation and differentiation in many cell types including keratinocytes.

Spatio-temporal sequence of cross-regulatory events in root meristem growth.

A central question in developmental biology is how multicellular organisms coordinate cell division and differentiation to determine organ size. In Arabidopsis roots, this balance is controlled by cytokinin-induced expression of SHORT HYPOCOTYL 2 (SHY2) in the so-called transition zone of the meristem, where SHY2 negatively regulates auxin response factors (ARFs) by protein-protein interaction. The resulting down-regulation of PIN-FORMED (PIN) auxin efflux carriers is considered the key event in promoting differentiation of meristematic cells. Here we show that this regulation involves additional, intermediary factors and is spatio-temporally constrained. We found that the described cytokinin-auxin crosstalk antagonizes BREVIS RADIX (BRX) activity in the developing protophloem. BRX is an auxin-responsive target of the prototypical ARF MONOPTEROS (MP), a key promoter of vascular development, and transiently enhances PIN3 expression to promote meristem growth in young roots. At later stages, cytokinin induction of SHY2 in the vascular transition zone restricts BRX expression to down-regulate PIN3 and thus limit meristem growth. Interestingly, proper SHY2 expression requires BRX, which could reflect feedback on the auxin responsiveness of SHY2 because BRX protein can directly interact with MP, likely acting as a cofactor. Thus, cross-regulatory antagonism between BRX and SHY2 could determine ARF activity in the protophloem. Our data suggest a model in which the regulatory interactions favor BRX expression in the early proximal meristem and SHY2 prevails because of supplementary cytokinin induction in the later distal meristem. The complex equilibrium of this regulatory module might represent a universal switch in the transition toward differentiation in various developmental contexts.

Biosynthesis of pyochelin and dihydroaeruginoic acid requires the iron-regulated pchDCBA operon in Pseudomonas aeruginosa.

The high-affinity siderophore salicylate is an intermediate in the biosynthetic pathway of pyochelin, another siderophore and chelator of transition metal ions, in Pseudomonas aeruginosa. The 2.5-kb region upstream of the salicylate biosynthetic genes pchBA was sequenced and found to contain two additional, contiguous genes, pchD and pchC, having the same orientation. The deduced amino acid sequence of the 60-kDa PchD protein was similar to those of the EntE protein (2,3-dihydroxybenzoate-AMP ligase) of Escherichia coli and other adenylate-forming enzymes, suggesting that salicylate might be adenylated at the carboxyl group by PchD. The 28-kDa PchC protein showed similarities to thioesterases of prokaryotic and eukaryotic origin and might participate in the release of the product(s) formed from activated salicylate. One potential product, dihydroaeruginoate (Dha), was identified in culture supernatants of iron-limited P. aeruginosa cells. The antifungal antibiotic Dha is thought to arise from the reaction of salicylate with cysteine, followed by cyclization of cysteine. Inactivation of the chromosomal pchD gene by insertion of the transcription and translation stop element omega Sm/Sp abolished the production of Dha and pyochelin, implying that PchD-mediated activation of salicylate may be a common first step in the synthesis of both metabolites. Furthermore, the pchD::omega Sm/Sp mutation had a strong polar effect on the expression of the pchBA genes, i.e., on salicylate synthesis, indicating that the pchDCBA genes constitute a transcriptional unit. A full-length pchDCBA transcript of ca. 4.4 kb could be detected in iron-deprived, growing cells of P. aeruginosa. Transcription of pchD started at tandemly arranged promoters, which overlapped with two Fur boxes (binding sites for the ferric uptake regulator) and the promoter of the divergently transcribed pchR gene encoding an activator of pyochelin biosynthesis. This promoter arrangement allows tight iron-mediated repression of the pchDCBA operon.

Characterisation of the putative effector interaction site of the regulatory HbpR protein from Pseudomonas azelaica by site-directed mutagenesis.

Bacterial transcription activators of the XylR/DmpR subfamily exert their expression control via σ(54)-dependent RNA polymerase upon stimulation by a chemical effector, typically an aromatic compound. Where the chemical effector interacts with the transcription regulator protein to achieve activation is still largely unknown. Here we focus on the HbpR protein from Pseudomonas azelaica, which is a member of the XylR/DmpR subfamily and responds to biaromatic effectors such as 2-hydroxybiphenyl. We use protein structure modeling to predict folding of the effector recognition domain of HbpR and molecular docking to identify the region where 2-hydroxybiphenyl may interact with HbpR. A large number of site-directed HbpR mutants of residues in- and outside the predicted interaction area was created and their potential to induce reporter gene expression in Escherichia coli from the cognate P(C) promoter upon activation with 2-hydroxybiphenyl was studied. Mutant proteins were purified to study their conformation. Critical residues for effector stimulation indeed grouped near the predicted area, some of which are conserved among XylR/DmpR subfamily members in spite of displaying different effector specificities. This suggests that they are important for the process of effector activation, but not necessarily for effector specificity recognition.

Transcriptional regulatory patterns of the myelin basic protein and malic enzyme genes by the thyroid hormone receptors alpha1 and beta1.

While there is evidence that the two ubiquitously expressed thyroid hormone (T3) receptors, TRalpha1 and TRbeta1, have distinct functional specificities, the mechanism by which they discriminate potential target genes remains largely unexplained. In this study, we demonstrate that the thyroid hormone response elements (TRE) from the malic enzyme and myelin basic protein genes (METRE and MBPTRE) respectively, are not functionally equivalent. The METRE, which is a direct repeat motif with a 4-base pair gap between the two half-site hexamers binds thyroid hormone receptor as a heterodimer with 9-cis-retinoic acid receptor (RXR) and mediates a high T3-dependent activation in response to TRalpha1 or TRbeta1 in NIH3T3 cells. In contrast, the MBPTRE, which consists of an inverted palindrome formed by two hexamers spaced by 6 base pairs, confers an efficient transactivation by TRbeta1 but a poor transactivation by TRalpha1. While both receptors form heterodimers with RXR on MBPTRE, the poor transactivation by TRalpha1 correlates also with its ability to bind efficiently as a monomer. This monomer, which is only observed with TRalpha1 bound to MBPTRE, interacts neither with N-CoR nor with SRC-1, explaining its functional inefficacy. However, in Xenopus oocytes, in which RXR proteins are not detectable, the transactivation mediated by TRalpha1 and TRbeta1 is equivalent and independent of a RXR supply, raising the question of the identity of the thyroid hormone receptor partner in these cells. Thus, in mammalian cells, the binding characteristics of TRalpha1 to MBPTRE (i.e. high monomer binding efficiency and low transactivation activity) might explain the particular pattern of T3 responsiveness of MBP gene expression during central nervous system development.

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.

Regulatory Mechanisms involved in Th2 Cell Differentiation. A Proteomics Approach.

Th2-solujen erilaistumista ohjaavat säätelyverkostot ja niiden tutkiminen proteomiikan avulla Astma ja allergiat ovat laajalle levinneitä ja vakavia sairauksia, joista kärsivät miljoonat ihmiset ympäri maailmaa. Koe-eläimillä tehdyt tutkimukset osoittavat, että interleukiini-4 (IL-4) on tärkeä allergisen astman ja allergioiden kehittymiselle ja kroonistumiselle. Se ohjaa T-auttajasolujen (Th-solujen) kehittymistä Th2-tyypin soluiksi, joilla on merkittävä rooli näiden tautien puhkeamisessa. Th2-solut tuottavat myös itse IL-4:ä, joka edesauttaa taudin seuraavien vaiheiden kehittymistä. Erityisesti STAT6-proteiini, joka aktivoituu IL-4-stimulaation seurauksena, on tarpeen Th2- vasteen syntymiselle ja kroonistumiselle antigeenin aiheuttamassa keuhkoputkien astmaattisessa tulehduksessa. Väitöskirjatyöni tarkoituksena oli käyttää kaksidimensionaaliseen elektroforeesiin (2- DE) perustuvaa proteomiikkaa ja massaspektrometriaa uusien Th2-solujen erilaistumista säätelevien proteiinien tunnistamiseksi. Erilaistumattomat Th-solut eristettiin vastasyntyneen napaverestä tai hiiren pernasta. Solut aktivoitiin Tsolureseptorin ja ns. ko-stimulatoristen reseptorien kautta ja erilaistettiin joko Th1- tai Th2-suuntaan vastaavasti erilaistavien IL-12- ja IL-4-sytokiinien avulla. Ensimmäisessä tutkimuksessa in vitro -erilaistettujen Th1- ja Th2-solujen proteomeja verrattiin keskenään proteiinien ilmenemisessä tai proteiinimodifikaatioissa olevien erojen tunnistamiseksi. Kaksi muuta päätutkimusta keskittyivät IL-4:n aiheuttamaan proteiinitason säätelyyn ensimmäisen vuorokauden aikana T-soluaktivaation jälkeen. Näistä ensimmäisessä IL-4:n aiheuttamia eroja tunnistettiin aktivoiduista ihmisen Thsoluista. IL-4:n todettiin säätelevän useita proteiineja kaspaasien välittämissä signalointiteissä sekä lisäävän T-solujen elävyyttä ja aktivoitumista. Toisessa tutkimuksessa STAT6-poistogeenisten hiirien lymfosyyttien proteomia verrattiin villityypin kontrollisoluihin T-soluaktivaation ja IL-4-stimulaation jälkeen. Näissä tutkimuksissa karakterisoitiin useita uusia IL-4:n ja STAT6:n kohdeproteiineja ja löydettiin uusia säätelyverkostoja. Tutkimustulokset ovat johtaneet uusiin Th2-erilaistumismekanismeja koskeviin hypoteeseihin.

Essential PchG-dependent reduction in pyochelin biosynthesis of Pseudomonas aeruginosa.

The biosynthetic genes pchDCBA and pchEF, which are known to be required for the formation of the siderophore pyochelin and its precursors salicylate and dihydroaeruginoate (Dha), are clustered with the pchR regulatory gene on the chromosome of Pseudomonas aeruginosa. The 4.6-kb region located downstream of the pchEF genes was found to contain three additional, contiguous genes, pchG, pchH, and pchI, probably forming a pchEFGHI operon. The deduced amino acid sequences of PchH and PchI are similar to those of ATP binding cassette transport proteins with an export function. PchG is a homolog of the Yersinia pestis and Y. enterocolitica proteins YbtU and Irp3, which are involved in the biosynthesis of yersiniabactin. A null mutation in pchG abolished pyochelin formation, whereas mutations in pchH and pchI did not affect the amounts of salicylate, Dha, and pyochelin produced. The pyochelin biosynthetic genes were expressed from a vector promoter, uncoupling them from Fur-mediated repression by iron and PchR-dependent induction by pyochelin. In a P. aeruginosa mutant lacking the entire pyochelin biosynthetic gene cluster, the expressed pchDCBA and pchEFG genes were sufficient for salicylate, Dha, and pyochelin production. Pyochelin formation was also obtained in the heterologous host Escherichia coli expressing pchDCBA and pchEFG together with the E. coli entD gene, which provides a phosphopantetheinyl transferase necessary for PchE and PchF activation. The PchG protein was purified and used in combination with PchD and phosphopantetheinylated PchE and PchF in vitro to produce pyochelin from salicylate, L-cysteine, ATP, NADPH, and S-adenosylmethionine. Based on this assay, a reductase function was attributed to PchG. In summary, this study completes the identification of the biosynthetic genes required for pyochelin formation from chorismate in P. aeruginosa.

Grx5 is a mitochondrial glutaredoxin required for the activity of iron/sulfur enzymes

Yeast cells contain a family of three monothiol glutaredoxins: Grx3, 4, and 5. Absence of Grx5 leads to constitutive oxidative damage, exacerbating that caused by external oxidants. Phenotypic defects associated with the absence of Grx5 are suppressed by overexpression ofSSQ1 and ISA2, two genes involved in the synthesis and assembly of iron/sulfur clusters into proteins. Grx5 localizes at the mitochondrial matrix, like other proteins involved in the synthesis of these clusters, and the mature form lacks the first 29 amino acids of the translation product. Absence of Grx5 causes: 1) iron accumulation in the cell, which in turn could promote oxidative damage, and 2) inactivation of enzymes requiring iron/sulfur clusters for their activity. Reduction of iron levels in grx5 null mutants does not restore the activity of iron/sulfur enzymes, and cell growth defects are not suppressed in anaerobiosis or in the presence of disulfide reductants. Hence, Grx5 forms part of the mitochondrial machinery involved in the synthesis and assembly of iron/sulfur centers.

Evolution and cellular function of monothiol glutaredoxins: involvement in iron-sulphur cluster assembly

A number of bacterial species, mostly proteobacteria, possess monothiol glutaredoxins homologous to the Saccharomyces cerevisiae mitochondrial protein Grx5, which is involved in iron–sulphur cluster synthesis. Phylogenetic profiling is used to predict that bacterial monothiol glutaredoxins also participate in the iron–sulphur cluster (ISC) assembly machinery, because their phylogenetic profiles are similar to the profiles of the bacterial homologues of yeast ISC proteins. High evolutionary cooccurrence is observed between the Grx5 homologues and the homologues of the Yah1 ferredoxin, the scaffold proteins Isa1 and Isa2, the frataxin protein Yfh1 and the Nfu1 protein. This suggests that a specific functional interaction exists between these ISC machinery proteins. Physical interaction analyses using low-definition protein docking predict the formation of strong and specific complexes between Grx5 and several components of the yeast ISC machinery. Two-hybrid analysis has confirmed the in vivo interaction between Grx5 and Isa1. Sequence comparison techniques and cladistics indicate that the other two monothiol glutaredoxins of S. cerevisiae, Grx3 and Grx4, have evolved from the fusion of a thioredoxin gene with a monothiol glutaredoxin gene early in the eukaryotic lineage, leading to differential functional specialization. While bacteria do not contain these chimaeric glutaredoxins, in many eukaryotic species Grx5 and Grx3/4-type monothiol glutaredoxins coexist in the cell.

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

Establishment of a Developmental Compartment Requires Interactions between Three Synergistic Cis-regulatory Modules.

The subdivision of cell populations in compartments is a key event during animal development. In Drosophila, the gene apterous (ap) divides the wing imaginal disc in dorsal vs ventral cell lineages and is required for wing formation. ap function as a dorsal selector gene has been extensively studied. However, the regulation of its expression during wing development is poorly understood. In this study, we analyzed ap transcriptional regulation at the endogenous locus and identified three cis-regulatory modules (CRMs) essential for wing development. Only when the three CRMs are combined, robust ap expression is obtained. In addition, we genetically and molecularly analyzed the trans-factors that regulate these CRMs. Our results propose a three-step mechanism for the cell lineage compartment expression of ap that includes initial activation, positive autoregulation and Trithorax-mediated maintenance through separable CRMs.

Ribosome profiling reveals the rhythmic liver translatome and circadian clock regulation by upstream open reading frames.

Mammalian gene expression displays widespread circadian oscillations. Rhythmic transcription underlies the core clock mechanism, but it cannot explain numerous observations made at the level of protein rhythmicity. We have used ribosome profiling in mouse liver to measure the translation of mRNAs into protein around the clock and at high temporal and nucleotide resolution. We discovered, transcriptome-wide, extensive rhythms in ribosome occupancy and identified a core set of approximately 150 mRNAs subject to particularly robust daily changes in translation efficiency. Cycling proteins produced from nonoscillating transcripts revealed thus-far-unknown rhythmic regulation associated with specific pathways (notably in iron metabolism, through the rhythmic translation of transcripts containing iron responsive elements), and indicated feedback to the rhythmic transcriptome through novel rhythmic transcription factors. Moreover, estimates of relative levels of core clock protein biosynthesis that we deduced from the data explained known features of the circadian clock better than did mRNA expression alone. Finally, we identified uORF translation as a novel regulatory mechanism within the clock circuitry. Consistent with the occurrence of translated uORFs in several core clock transcripts, loss-of-function of Denr, a known regulator of reinitiation after uORF usage and of ribosome recycling, led to circadian period shortening in cells. In summary, our data offer a framework for understanding the dynamics of translational regulation, circadian gene expression, and metabolic control in a solid mammalian organ.