Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Glucose release from GLUT2-null hepatocytes: characterization of a major and a minor pathway.

We previously reported that glucose can be released from GLUT2-null hepatocytes through a membrane traffic-based pathway issued from the endoplasmic reticulum. Here, we further characterized this glucose release mechanism using biosynthetic labeling protocols. In continuous pulse-labeling experiments, we determined that glucose secretion proceeded linearly and with the same kinetics in control and GLUT2-null hepatocytes. In GLUT2-deficient hepatocytes, however, a fraction of newly synthesized glucose accumulated intracellularly. The linear accumulation of glucose in the medium was inhibited in mutant, but not in control, hepatocytes by progesterone and low temperature, as previously reported, but, importantly, also by microtubule disruption. The intracellular pool of glucose was shown to be present in the cytosol, and, in pulse-chase experiments, it was shown to be released at a relatively slow rate. Release was not inhibited by S-4048 (an inhibitor of glucose-6-phosphate translocase), cytochalasin B, or progesterone. It was inhibited by phloretin, carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone, and low temperature. We conclude that the major release pathway segregates glucose away from the cytosol by use of a membrane traffic-based, microtubule-dependent mechanism and that the release of the cytosolic pool of newly synthesized glucose, through an as yet unidentified plasma membrane transport system, cannot account for the bulk of glucose release.

Arabidopsis jasmonate signaling pathway.

Jasmonates control defense gene expression, growth, and fertility throughout the plant kingdom and have been studied extensively in Arabidopsis thaliana. The prohormone jasmonic acid (JA) is conjugated to amino acids such as isoleucine to form the active hormone jasmonoyl-isoleucine (JA-Ile). A series of breakthroughs has identified the SCF [SCF consists of four subunits: a cullin, SKP1 (S-phase kinase-associated protein 1), a RING finger protein (RBX1/HRT1/ROC1), and an F-box protein] CORONATINE INSENSITIVE1 (COI1) E3 ubiquitin ligase complex and the JASMONATE ZIM-DOMAIN (JAZ) proteins as central components in the perception of and transcriptional response to JA-Ile. JAZ proteins (most probably as dimers) bind transcription factors such as MYC2 before JA-Ile production. JA-Ile binds to COI1 to facilitate the formation of COI1-JAZ complexes, leading to ubiquitination and subsequent degradation of JAZ proteins. The degradation of JAZ proteins liberates transcription factors that function in the presence of the RNA polymerase II coregulatory complex Mediator to permit the expression of a number of jasmonate-regulated genes. Recent developments include the identification of COI1 as a receptor for jasmonates. Upstream of the signaling events, microRNA319 (miR319) negatively regulates the production of JA and JA-derived signals.

Plants and tortoises: mutations in the Arabidopsis jasmonate pathway increase feeding in a vertebrate herbivore.

Photosynthetic tissues, the major food source of many invertebrates and vertebrates, are well defended. Many defence traits in leaves are controlled via the jasmonate signalling pathway in which jasmonate acts as a hormone by binding to a receptor to activate responses that lead to increased resistance to invertebrate folivores. We predicted that mutations in jasmonate synthesis might also increase the vulnerability of leaves to vertebrate folivores and tested this hypothesis using the Eastern Hermann's tortoise (Eurotestudo boettgeri) and an Arabidopsis thaliana (Brassicaceae) allene oxide synthase (aos) mutant unable to synthesize jasmonate. Tortoises preferred the aos mutant over the wild type (WT). Based on these results, we then investigated the effect of mutating jasmonate perception using a segregating population of the recessive A. thaliana jasmonate receptor mutant coronatine insensitive1-1 (coi1-1). Genotyping of these plants after tortoise feeding revealed that the homozygous coi1-1 receptor mutant was consumed more readily than the heterozygous mutant or the WT. Therefore, the plant's ability to synthesize or perceive jasmonate reduces feeding by a vertebrate herbivore. We also tested whether or not tortoise feeding behaviour was influenced by glucosinolates, the principal defence chemicals in Arabidopsis leaves with known roles in defence against many generalist insects. However, in contrast to what has been observed with such insects, leaves in which the levels of these compounds were reduced genetically were consumed at a similar rate to those of the WT.

Integrated Ladinian bio-chronostratigraphy and geochrononology of Monte San Giorgio (Southern Alps, Switzerland)

New biostratigraphic data significantly improve the age assignment of the Ladinian succession of Monte San Giorgio (UNESCO World Heritage List site, Southern Alps, Switzerland), whose world-famous fossil marine vertebrate faunas are now dated to the substage and zone levels. High-resolution single-zircon U-Pb dating was performed using ID-TIMS and chemical abrasion (CA) pre-treatment technique on volcanic ash layers intercalated in the biostratigraphically-defined intervals of the Meride Limestone. It yielded ages of 241.07 +/- 0.13 Ma (Cava superiore beds, P. gredleri Zone), 240.63 +/- 0.13 Ma (Cassina beds, P gredleri/P. archelaus transition Zone) and 239.51 +/- 0.15 Ma (Lower Kalkschieferzone, P. archelaus Zone). Our results suggest that the time interval including the vertebrate-bearing Middle Triassic section spans around 4 Myr and is thus significantly shorter than so far assumed. The San Giorgio Dolomite and the Meride Limestone correlate with intervals of the Buchenstein Formation and the Wengen Formation in the reference section at Bagolino, where the Global boundary Stratotype Section and Point (GSSP) for the base of the Ladinian was defined. The new radio-isotopic ages of the Meride Limestone are up to 2 Myr older than those published for the biostratigraphically-equivalent intervals at Bagolino but they are consistent with the recent re-dating of the underlying Besano Formation, also performed using the CA technique. Average sedimentation rates at Monte San Giorgio are by more than an order of magnitude higher compared to those assumed for the Buchenstein Formation, which formed under sediment-starved pelagic conditions, and reflect prevailing high subsidence and high carbonate mud supply from the adjoining Salvatore/Esino platforms. Finally, the high-resolution U-Pb ages allow a correlation of the vertebrate faunas of the Cava superiore/Cava inferiore beds with the marine vertebrate record of the Prosanto Formation (Upper Austroalpine), so far precluded by the poor biostratigraphic control of the latter.

The puzzles of the prokineticin 2 pathway in human reproduction.

Prokineticin, 1 (PROK1) and prokineticin 2 (PROK2), are two closely related proteins that were identified as the mammalian homologs of their two amphibian homologs, mamba intestinal toxin (MIT-1) and Bv8. MIT-1 was initially identified as a non-toxic constituent in the venom of the black mamba snake (Dendroaspis polylepis) (Joubert and Strydom, 1980) while Bv8 was identified in the skin secretion of the toad, Bombina variegate (Mollay et al., 1999). All three homologs stimulate gastrointestinal motility thus accounting for their family name "prokineticins" (Schweitz et al., 1990, 1999). However, since its initial description, both PROK1 and PROK2 have been found to regulate a dazzling array of biological functions throughout the body. In particular, PROK1 acts as a potent angiogenic mitogen on endocrine vascular epithelium, thus earning its other name, Endocrine gland-vascular endothelial factor (EG-VEGF) (LeCouter et al., 2002). In contrast, the PROK2 signaling pathway is a critical regulator of olfactory bulb morphogenesis and sexual maturation in mammals and this function is the focus of this review.

Pharmacological inhibition of nicotinamide phosphoribosyltransferase/visfatin enzymatic activity identifies a new inflammatory pathway linked to NAD.

Nicotinamide phosphoribosyltransferase (NAMPT), also known as visfatin, is the rate-limiting enzyme in the salvage pathway of NAD biosynthesis from nicotinamide. Since its expression is upregulated during inflammation, NAMPT represents a novel clinical biomarker in acute lung injury, rheumatoid arthritis, and Crohn's disease. However, its role in disease progression remains unknown. We report here that NAMPT is a key player in inflammatory arthritis. Increased expression of NAMPT was confirmed in mice with collagen-induced arthritis, both in serum and in the arthritic paw. Importantly, a specific competitive inhibitor of NAMPT effectively reduced arthritis severity with comparable activity to etanercept, and decreased pro-inflammatory cytokine secretion in affected joints. Moreover, NAMPT inhibition reduced intracellular NAD concentration in inflammatory cells and circulating TNFalpha levels during endotoxemia in mice. In vitro pharmacological inhibition of NAMPT reduced the intracellular concentration of NAD and pro-inflammatory cytokine secretion by inflammatory cells. Thus, NAMPT links NAD metabolism to inflammatory cytokine secretion by leukocytes, and its inhibition might therefore have therapeutic efficacy in immune-mediated inflammatory disorders.

Novel acid phosphatase in Candida glabrata suggests selective pressure and niche specialization in the phosphate signal transduction pathway.

Evolution through natural selection suggests unnecessary genes are lost. We observed that the yeast Candida glabrata lost the gene encoding a phosphate-repressible acid phosphatase (PHO5) present in many yeasts including Saccharomyces cerevisiae. However, C. glabrata still had phosphate starvation-inducible phosphatase activity. Screening a C. glabrata genomic library, we identified CgPMU2, a member of a three-gene family that contains a phosphomutase-like domain. This small-scale gene duplication event could allow for sub- or neofunctionalization. On the basis of phylogenetic and biochemical characterizations, CgPMU2 has neofunctionalized to become a broad range, phosphate starvation-regulated acid phosphatase, which functionally replaces PHO5 in this pathogenic yeast. We determined that CgPmu2, unlike ScPho5, is not able to hydrolyze phytic acid (inositol hexakisphosphate). Phytic acid is present in fruits and seeds where S. cerevisiae grows, but is not abundant in mammalian tissues where C. glabrata grows. We demonstrated that C. glabrata is limited from an environment where phytic acid is the only source of phosphate. Our work suggests that during evolutionary time, the selection for the ancestral PHO5 was lost and that C. glabrata neofunctionalized a weak phosphatase to replace PHO5. Convergent evolution of a phosphate starvation-inducible acid phosphatase in C. glabrata relative to most yeast species provides an example of how small changes in signal transduction pathways can mediate genetic isolation and uncovers a potential speciation gene.

Blocking the B7-H4 pathway with novel recombinant antibodies enhances T cell-mediated antitumor responses.

B7-H4 inhibits T-cell activation and is widely expressed by solid neoplasms. We have recently demonstrated that the expression of B7-H4 on the surface of malignant cells in vivo is inducible, and that novel anti-B7-H4 recombinant antibodies can reverse the inhibition of tumor-specific T cells. Thus, antibodies targeting the B7-H4 pathways may extend the survival of cancer patients by restoring T cell-mediated antitumor responses.

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.

ER stress activates the NLRP3 inflammasome via an UPR-independent pathway.

Uncontrolled endoplasmic reticulum (ER) stress responses are proposed to contribute to the pathology of chronic inflammatory diseases such as type 2 diabetes or atherosclerosis. However, the connection between ER stress and inflammation remains largely unexplored. Here, we show that ER stress causes activation of the NLRP3 inflammasome, with subsequent release of the pro-inflammatory cytokine interleukin-1β. This ER-triggered proinflammatory signal shares the same requirement for reactive oxygen species production and potassium efflux compared with other known NLRP3 inflammasome activators, but is independent of the classical unfolded protein response (UPR). We thus propose that the NLRP3 inflammasome senses and responds to ER stress downstream of a previously uncharacterized ER stress response signaling pathway distinct from the UPR, thus providing mechanistic insight to the link between ER stress and chronic inflammatory diseases.

Hypoxia-induced inhibition of epithelial Na(+) channels in the lung. Role of Nedd4-2 and the ubiquitin-proteasome pathway.

Transepithelial sodium transport via alveolar epithelial Na(+) channels (ENaC) and Na(+),K(+)-ATPase constitutes the driving force for removal of alveolar edema fluid. Alveolar hypoxia associated with pulmonary edema may impair ENaC activity and alveolar Na(+) absorption through a decrease of ENaC subunit expression at the apical membrane of alveolar epithelial cells (AECs). Here, we investigated the mechanism(s) involved in this process in vivo in the β-Liddle mouse strain mice carrying a truncation of β-ENaC C-terminus abolishing the interaction between β-ENaC and the ubiquitin protein-ligase Nedd4-2 that targets the channel for endocytosis and degradation and in vitro in rat AECs. Hypoxia (8% O2 for 24 h) reduced amiloride-sensitive alveolar fluid clearance by 69% in wild-type mice but had no effect in homozygous mutated β-Liddle littermates. In vitro, acute exposure of AECs to hypoxia (0.5-3% O2 for 1-6 h) rapidly decreased transepithelial Na(+) transport as assessed by equivalent short-circuit current Ieq and the amiloride-sensitive component of Na(+) current across the apical membrane, reflecting ENaC activity. Hypoxia induced a decrease of ENaC subunit expression in the apical membrane of AECs with no change in intracellular expression and induced a 2-fold increase in α-ENaC polyubiquitination. Hypoxic inhibition of amiloride-sensitive Ieq was fully prevented by preincubation with the proteasome inhibitors MG132 and lactacystin or with the antioxidant N-acetyl-cysteine. Our data strongly suggest that Nedd4-2-mediated ubiquitination of ENaC leading to endocytosis and degradation of apical Na(+) channels is a key feature of hypoxia-induced inhibition of transepithelial alveolar Na(+) transport.

Evidence for different expression profiles for c-Met, EGFR, PTEN and the mTOR pathway in low and high grade endometrial carcinomas in a cohort of consecutive women. Occurrence of PIK3CA and K-Ras mutations and microsatellite instability.

Molecular and genetic investigations in endometrial carcinogenesis may have prognostic and therapeutic implications. We studied the expression of EGFR, c-Met, PTEN and the mTOR signalling pathway (phospho-AKT/phospho-mTOR/phospho-RPS6) in 69 consecutive tumours and 16 tissue microarrays. We also analysed PIK3CA, K-Ras mutations and microsatellite instability (MSI). We distinguished two groups: group 1 (grade 1 and 2 endometrioid cancers) and group 2 (grade 3 endometrioid and type II clear and serous cell cancers). We hypothesised that these histological groups might have different features. We found that a) survival was higher in group 1 with less aggressive tumours (P⟨0.03); b) EGFR (P=0.01), PTEN and the AKT/mTOR/RPS6 signalling pathway were increased in group 1 versus group 2 (P=0.05 for phospho-mTOR); c) conversely, c-Met was higher (P⟨0.03) in group 2 than in group 1; d) In group 1, EGFR was correlated with c-Met, phospho-mTOR, phospho-RPS6 and the global activity of the phospho-AKT/phospho-mTOR/phospho-RPS6 pathway. In group 2, EGFR was correlated only with the phospho-AKT/phospho-mTOR/phospho-RPS6 pathway, whereas c-Met was correlated with PTEN; e) survival was higher for tumours with more than 50% PTEN-positive cells; f) K-RAS and PIK3CA mutations occurred in 10-12% of the available tumours and MSI in 40.4%, with a loss of MLH1 and PMS2 expression. Our results for endometrial cancers provide the first evidence for a difference in status between groups 1 and 2. The patients may benefit from different targeted treatments, anti-EGFR agents and rapamycin derivatives (anti-mTOR) for group 1 and an anti c-MET/ligand complex for group 2.

An integrated encyclopedia of DNA elements in the human genome.

The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.

Nedd4-2 modulates renal Na+-Cl- cotransporter via the aldosterone-SGK1-Nedd4-2 pathway.

Regulation of renal Na(+) transport is essential for controlling blood pressure, as well as Na(+) and K(+) homeostasis. Aldosterone stimulates Na(+) reabsorption by the Na(+)-Cl(-) cotransporter (NCC) in the distal convoluted tubule (DCT) and by the epithelial Na(+) channel (ENaC) in the late DCT, connecting tubule, and collecting duct. Aldosterone increases ENaC expression by inhibiting the channel's ubiquitylation and degradation; aldosterone promotes serum-glucocorticoid-regulated kinase SGK1-mediated phosphorylation of the ubiquitin-protein ligase Nedd4-2 on serine 328, which prevents the Nedd4-2/ENaC interaction. It is important to note that aldosterone increases NCC protein expression by an unknown post-translational mechanism. Here, we present evidence that Nedd4-2 coimmunoprecipitated with NCC and stimulated NCC ubiquitylation at the surface of transfected HEK293 cells. In Xenopus laevis oocytes, coexpression of NCC with wild-type Nedd4-2, but not its catalytically inactive mutant, strongly decreased NCC activity and surface expression. SGK1 prevented this inhibition in a kinase-dependent manner. Furthermore, deficiency of Nedd4-2 in the renal tubules of mice and in cultured mDCT(15) cells upregulated NCC. In contrast to ENaC, Nedd4-2-mediated inhibition of NCC did not require the PY-like motif of NCC. Moreover, the mutation of Nedd4-2 at either serine 328 or 222 did not affect SGK1 action, and mutation at both sites enhanced Nedd4-2 activity and abolished SGK1-dependent inhibition. Taken together, these results suggest that aldosterone modulates NCC protein expression via a pathway involving SGK1 and Nedd4-2 and provides an explanation for the well-known aldosterone-induced increase in NCC protein expression.