ChIPping away at gene regulation

The coordinated regulation of gene expression in higher eukaryotes is complex and poorly understood. Recent technological advances have allowed the first insights into these networks on a genome-wide scale. These investigations have identified transcription factor target sites in the genome and successfully predicted cooperative interactions with other factors. However, a detailed understanding of the processes that coordinate gene expression remains elusive. Here, we highlight the advances that have been made using current methods, and the need for new technologies to address the gaps in our knowledge and to map these complex pathways further.

Evolution, gene regulation and functional analysis of BMP2 in fish

Bone morphogenetic proteins (BMPs) are multifunctional growth factors belonging to the transforming growth factor β (TGFβ) superfamily with a central role in bone formation and mineralization. BMP2, a founding member of this family, has demonstrated remarkable osteogenic properties and is clinically used to promote bone repair and fracture healing. Lack of basic data on factors regulating BMP2 expression and activity have hampered a better understanding of its role in bone formation and bone-related diseases. The objective of this work was to collect new functional data and determine spatiotemporal expression patterns in a fish system aiming towards a better understanding of BMP2 function and regulation. Transcriptional and post-transcriptional regulation of gilthead seabream BMP2 gene was inferred from luciferase reporter systems. Several bone- and cartilage-related transcription factors (e.g. RUNX3, MEF2c, SOX9 and ETS1) were found to regulate BMP2 transcription, while microRNA 20a was shown to affect stability of the BMP2 transcript and thus the mineralogenic capacity of fish bone-derived host cells. The regulation of BMP2 activity through an interaction with the matrix Gla protein (MGP) was investigated in vitro using BMP responsive elements (BRE) coupled to luciferase reporter gene. Although we demonstrated the functionality of the experimental system in a fish cell line and the activation of BMP signaling pathway by seabream BMP2, no conclusive evidence could be collected on a possible interaction beween MGP and BMP2. The evolutionary relationship among the members of BMP2/4/16 subfamily was inferred from taxonomic and phylogenetic analyses. BMP16 diverged prior to BMP2 and BMP4 and should be the result of an ancient genome duplication that occurred early in vertebrate evolution. Structural and functional data suggested that all three proteins are effectors of the BMP signaling pathway, but expression data revealed different spatiotemporal patterns in teleost fish suggesting distinct mechanisms of regulation. In this work, through the collection of novel data, we provide additional insight into the regulation, the structure and the phylogenetic relationship of BMP2 and its closely related family members.

Role of the Protein Tyrosine Kinase 7 gene in human neural tube defects

Les anomalies du tube neural (ATN) sont des anomalies développementales où le tube neural reste ouvert (1-2/1000 naissances). Afin de prévenir cette maladie, une connaissance accrue des processus moléculaires est nécessaire. L’étiologie des ATN est complexe et implique des facteurs génétiques et environnementaux. La supplémentation en acide folique est reconnue pour diminuer les risques de développer une ATN de 50-70% et cette diminution varie en fonction du début de la supplémentation et de l’origine démographique. Les gènes impliqués dans les ATN sont largement inconnus. Les études génétiques sur les ATN chez l’humain se sont concentrées sur les gènes de la voie métabolique des folates du à leur rôle protecteur dans les ATN et les gènes candidats inférés des souris modèles. Ces derniers ont montré une forte association entre la voie non-canonique Wnt/polarité cellulaire planaire (PCP) et les ATN. Le gène Protein Tyrosine Kinase 7 est un membre de cette voie qui cause l’ATN sévère de la craniorachischisis chez les souris mutantes. Ptk7 interagit génétiquement avec Vangl2 (un autre gène de la voie PCP), où les doubles hétérozygotes montrent une spina bifida. Ces données font de PTK7 comme un excellent candidat pour les ATN chez l’humain. Nous avons re-séquencé la région codante et les jonctions intron-exon de ce gène dans une cohorte de 473 patients atteints de plusieurs types d’ATN. Nous avons identifié 6 mutations rares (fréquence allélique <1%) faux-sens présentes chez 1.1% de notre cohorte, dont 3 sont absentes dans les bases de données publiques. Une variante, p.Gly348Ser, a agi comme un allèle hypermorphique lorsqu'elle est surexprimée dans le modèle de poisson zèbre. Nos résultats impliquent la mutation de PTK7 comme un facteur de risque pour les ATN et supporte l'idée d'un rôle pathogène de la signalisation PCP dans ces malformations.

On the function of the Dictyostelium Argonaute A protein (AgnA) in epigenetic gene regulation

With molecular biology methods and bioinformatics, the Argonaute proteins in Dictyostelium discoideum were characterized, and the function of the AgnA protein in RNAi and DNA methylation was investigated, as well as cellular features. Also interaction partners of the PAZ-Piwi domain of AgnA (PAZ-PiwiAgnA) were discovered. The Dictyostelium genome encodes five Argonaute proteins, termed AgnA/B/C/D/E. The expression level of Argonaute proteins was AgnB/D/E > AgnA > AgnC. All these proteins contain the characteristic conserved of PAZ and Piwi domains. Fluorescence microscopy revealed that the overexpressed C-terminal GFP-fusion of PAZ-PiwiAgnA (PPWa-GFP) localized to the cytoplasm. Overexpression of PPWa-GFP leaded to an increased gene silencing efficiency mediated by RNAi but not by antisense RNA. This indicated that PAZ-PiwiAgnA is involved in the RNAi pathway, but not in the antisense pathway. An analysis of protein-protein interactions by a yeast-two-hybrid screen on a cDNA library from vegetatively grown Dictyostelium revealed that several proteins, such as EF2, EF1-I, IfdA, SahA, SamS, RANBP1, UAE1, CapA, and GpdA could interact with PAZ-PiwiAgnA. There was no interaction between PAZ-PiwiAgnA and HP1, HelF and DnmA detected by direct yeast-two-hybrid analysis. The fluorescence microscopy images showed that the overexpressed GFP-SahA or IfdA fusion proteins localized to both cytoplasm and nuclei, while the overexpressed GFP-SamS localized to the cytoplasm. The expression of SamS in AgnA knock down mutants was strongly down regulated on cDNA and mRNA level in, while the expression of SahA was only slightly down regulated. AgnA knock down mutants displayed defects in growth and phagocytosis, which suggested that AgnA affects also cell biological features. The inhibition of DNA methylation on DIRS-1 and Skipper retroelements, as well as the endogenous mvpB and telA gene, observed for the same strains, revealed that AgnA is involved in the DNA methylation pathway. Northern blot analysis showed that Skipper and DIRS-1 were rarely expressed in Ax2, but the expression of Skipper was upregulated in AgnA knock down mutants, while the expression of DIRS-1 was not changed. A knock out of the agnA gene failed even though the homologous recombination of the disruption construct occurred at the correct site, which indicated that there was a duplication of the agnA gene in the genome. The same phenomenon was also observed in ifdA knock out experiments.

Global iron-dependent gene regulation in Escherichia coli - A new mechanism for iron homeostasis

Organisms generally respond to iron deficiency by increasing their capacity to take up iron and by consuming intracellular iron stores. Escherichia coli, in which iron metabolism is particularly well understood, contains at least 7 iron-acquisition systems encoded by 35 iron-repressed genes. This Fe-dependent repression is mediated by a transcriptional repressor, Fur ( ferric uptake regulation), which also controls genes involved in other processes such as iron storage, the Tricarboxylic Acid Cycle, pathogenicity, and redox-stress resistance. Our macroarray-based global analysis of iron- and Fur-dependent gene expression in E. coli has revealed several novel Fur-repressed genes likely to specify at least three additional iron- transport pathways. Interestingly, a large group of energy metabolism genes was found to be iron and Fur induced. Many of these genes encode iron- rich respiratory complexes. This iron- and Fur-dependent regulation appears to represent a novel iron-homeostatic mechanism whereby the synthesis of many iron- containing proteins is repressed under iron- restricted conditions. This mechanism thus accounts for the low iron contents of fur mutants and explains how E. coli can modulate its iron requirements. Analysis of Fe-55-labeled E. coli proteins revealed a marked decrease in iron- protein composition for the fur mutant, and visible and EPR spectroscopy showed major reductions in cytochrome b and d levels, and in iron- sulfur cluster contents for the chelator-treated wild-type and/or fur mutant, correlating well with the array and quantitative RT-PCR data. In combination, the results provide compelling evidence for the regulation of intracellular iron consumption by the Fe2+-Fur complex.

Chronic hypoxia in the human neuroblastoma SH-SY5Y causes reduced expression of the putative alpha-secretases, ADAM10 and TACE, without altering their mRNA levels

Alzheimer's disease is more frequent following an ischemic or hypoxic episode, with levels of beta-amyloid peptides elevated in brains from patients. Similar increases are found after experimental ischemia in animals. It is possible that increased beta-amyloid deposition arises from alterations in amyloid precursor protein (APP) metabolism, indeed, we have shown that exposing cells of neuronal origin to chronic hypoxia decreased the secretion of soluble APP (sAPPalpha) derived by action of alpha-secretase on APP, coinciding with a decrease in protein levels of ADAM10, a disintegrin metalloprotease which is thought to be the major alpha-secretase. In the current study, we extended those observations to determine whether the expression of ADAM10 and another putative alpha-secretase, TACE, as well as the beta-secretase, BACE1 were regulated by chronic hypoxia at the level of protein and mRNA. Using Western blotting and RT-PCR, we demonstrate that after 48 h chronic hypoxia, such that sAPPalpha secretion is decreased by over 50%, protein levels of ADAM10 and TACE and by approximately 60% and 40% respectively with no significant decrease in BACE1 levels. In contrast, no change in the expression of the mRNA for these proteins could be detected. Thus, we conclude that under CH the level of the putative alpha-secretases, ADAM10 and TACE are regulated by post-translational mechanisms, most probably proteolysis, rather than at the level of transcription.

p90 ribosomal S6 kinases play a significant role in early gene regulation in the cardiomyocyte response to Gq protein-coupled receptor stimuli, endothelin-1 and α1-adrenergic receptor agonists

Extracellular signal-regulated kinases 1/2 (ERK1/2) and their substrates, p90 ribosomal S6 kinases (RSKs), phosphorylate different transcription factors, contributing differentially to transcriptomic profiles. In cardiomyocytes, ERK1/2 are required for >70% of the transcriptomic response to endothelin-1. Here, we investigated the role of RSKs in the transcriptomic responses to Gq protein-coupled receptor agonists, endothelin-1, phenylephrine (generic α1-adrenergic receptor agonist) and A61603 (α1A-adrenergic receptor selective). Phospho-ERK1/2 and phospho-RSKs appeared in cardiomyocyte nuclei within 2-3 min of stimulation (endothelin-1>a61603≈phenylephrine). All agonists increased nuclear RSK2, but only endothelin-1 increased nuclear RSK1 content. PD184352 (inhibits ERK1/2 activation) and BI-D1870 (inhibits RSKs) were used to dissect the contribution of RSKs to the endothelin-1-responsive transcriptome. Of 213 RNAs upregulated at 1 h, 51% required RSKs for upregulation whereas 29% required ERK1/2 but not RSKs. The transcriptomic response to phenylephrine overlapped with, but was not identical to, endothelin-1. As with endothelin-1, PD184352 inhibited upregulation of most phenylephrine-responsive transcripts, but the greater variation in effects of BI-D1870 suggests that differential RSK signalling influences global gene expression. A61603 induced similar changes in RNA expression in cardiomyocytes as phenylephrine, indicating that the signal was mediated largely through α1A-adrenergic receptors. A61603 also increased expression of immediate early genes in perfused adult rat hearts and, as in cardiomyocytes, upregulation of the majority of genes was inhibited by PD184352. PD184352 or BI-D1870 prevented the increased surface area induced by endothelin-1 in cardiomyocytes. Thus, RSKs play a significant role in regulating cardiomyocyte gene expression and hypertrophy in response to Gq protein-coupled receptor stimulation.

Modelling negative feedback networks for activating transcription factor 3 predicts a dominant role for miRNAs in immediate early gene regulation

Activating transcription factor 3 (Atf3) is rapidly and transiently upregulated in numerous systems, and is associated with various disease states. Atf3 is required for negative feedback regulation of other genes, but is itself subject to negative feedback regulation possibly by autorepression. In cardiomyocytes, Atf3 and Egr1 mRNAs are upregulated via ERK1/2 signalling and Atf3 suppresses Egr1 expression. We previously developed a mathematical model for the Atf3-Egr1 system. Here, we adjusted and extended the model to explore mechanisms of Atf3 feedback regulation. Introduction of an autorepressive loop for Atf3 tuned down its expression and inhibition of Egr1 was lost, demonstrating that negative feedback regulation of Atf3 by Atf3 itself is implausible in this context. Experimentally, signals downstream from ERK1/2 suppress Atf3 expression. Mathematical modelling indicated that this cannot occur by phosphorylation of pre-existing inhibitory transcriptional regulators because the time delay is too short. De novo synthesis of an inhibitory transcription factor (ITF) with a high affinity for the Atf3 promoter could suppress Atf3 expression, but (as with the Atf3 autorepression loop) inhibition of Egr1 was lost. Developing the model to include newly-synthesised miRNAs very efficiently terminated Atf3 protein expression and, with a 4-fold increase in the rate of degradation of mRNA from the mRNA/miRNA complex, profiles for Atf3 mRNA, Atf3 protein and Egr1 mRNA approximated to the experimental data. Combining the ITF model with that of the miRNA did not improve the profiles suggesting that miRNAs are likely to play a dominant role in switching off Atf3 expression post-induction.

Symmetry in biology: from genetic code to stochastic gene regulation

Mathematical models, as instruments for understanding the workings of nature, are a traditional tool of physics, but they also play an ever increasing role in biology - in the description of fundamental processes as well as that of complex systems. In this review, the authors discuss two examples of the application of group theoretical methods, which constitute the mathematical discipline for a quantitative description of the idea of symmetry, to genetics. The first one appears, in the form of a pseudo-orthogonal (Lorentz like) symmetry, in the stochastic modelling of what may be regarded as the simplest possible example of a genetic network and, hopefully, a building block for more complicated ones: a single self-interacting or externally regulated gene with only two possible states: ` on` and ` off`. The second is the algebraic approach to the evolution of the genetic code, according to which the current code results from a dynamical symmetry breaking process, starting out from an initial state of complete symmetry and ending in the presently observed final state of low symmetry. In both cases, symmetry plays a decisive role: in the first, it is a characteristic feature of the dynamics of the gene switch and its decay to equilibrium, whereas in the second, it provides the guidelines for the evolution of the coding rules.

GAD2 on chromosome 10p12 is a candidate gene for human obesity

The gene GAD2 encoding the glutamic acid decarboxylase enzyme (GAD65) is a positional candidate gene for obesity on Chromosome 10p11–12, a susceptibility locus for morbid obesity in four independent ethnic populations. GAD65 catalyzes the formation of γ-aminobutyric acid (GABA), which interacts with neuropeptide Y in the paraventricular nucleus to contribute to stimulate food intake. A case-control study (575 morbidly obese and 646 control subjects) analyzing GAD2 variants identified both a protective haplotype, including the most frequent alleles of single nucleotide polymorphisms (SNPs) +61450 C>A and +83897 T>A (OR = 0.81, 95% CI [0.681–0.972], p = 0.0049) and an at-risk SNP (−243 A>G) for morbid obesity (OR = 1.3, 95% CI [1.053–1.585], p = 0.014). Furthermore, familial-based analyses confirmed the association with the obesity of SNP +61450 C>A and +83897 T>A haplotype (χ2 = 7.637, p = 0.02). In the murine insulinoma cell line βTC3, the G at-risk allele of SNP −243 A>G increased six times GAD2 promoter activity (p < 0.0001) and induced a 6-fold higher affinity for nuclear extracts. The −243 A>G SNP was associated with higher hunger scores (p = 0.007) and disinhibition scores (p = 0.028), as assessed by the Stunkard Three-Factor Eating Questionnaire. As GAD2 is highly expressed in pancreatic β cells, we analyzed GAD65 antibody level as a marker of β-cell activity and of insulin secretion. In the control group, −243 A>G, +61450 C>A, and +83897 T>A SNPs were associated with lower GAD65 autoantibody levels (p values of 0.003, 0.047, and 0.006, respectively). SNP +83897 T>A was associated with lower fasting insulin and insulin secretion, as assessed by the HOMA-B% homeostasis model of β-cell function (p = 0.009 and 0.01, respectively). These data support the hypothesis of the orexigenic effect of GABA in humans and of a contribution of genes involved in GABA metabolism in the modulation of food intake and in the development of morbid obesity.

Chronic disease burden and human capital investment: evidence from the chagas disease campaign in Brazil

We investigate the effects of augmented life expectancy and health improvements on human capital investment, labor supply and fertility decisions. Our main motivation is the prediction of human capital theory that a longer and healthier life encourages educational investment and female labor force participation, while discouraging fertility. To assess the magnitude of these effects, we explore a national campaign against Chagas disease in Brazil as an exogenous source of adult mortality decline and improvement in health conditions. We show that, relative to non-endemic areas, previously endemic regions saw higher increases in educational investment, measured by literacy, school attendance and years of schooling, following the campaign. Additionally, we find that labor force participation increased in high prevalence areas relative to low prevalence ones. Furthermore, we estimate a substantially higher effect on female labor force participation relative to male, suggesting that longevity gains and health improvements affected women's incentives to work, encouraging women to join the labor force. We do not find significant effects on fertility decisions.

Aspects of gene regulation in the diploid and tetraploid Odontophrynus americanus (Amphibia, Anura, Leptodactylidae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)