EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

Allele variants in functional MicroRNA target sites of the neurotrophin-3 receptor gene (NTRK3) as susceptibility factors for anxiety disorders

Genetic and functional data indicate that variation in the expression of the neurotrophin-3 receptor gene (NTRK3) may have an impact on neuronal plasticity, suggesting a role for NTRK3 in the pathophysiology of anxiety disorders. MicroRNA (miRNA) posttranscriptional gene regulators act by base-pairing to specific sequence sites, usually at the 3'UTR of the target mRNA. Variants at these sites might result in gene expression changes contributing to disease susceptibility. We investigated genetic variation in two different isoforms of NTRK3 as candidate susceptibility factors for anxiety by resequencing their 3'UTRs in patients with panic disorder (PD), obsessive-compulsive disorder (OCD), and in controls. We have found the C allele of rs28521337, located in a functional target site for miR-485-3p in the truncated isoform of NTRK3, to be significantly associated with the hoarding phenotype of OCD. We have also identified two new rare variants in the 3'UTR of NTRK3, ss102661458 and ss102661460, each present only in one chromosome of a patient with PD. The ss102661458 variant is located in a functional target site for miR-765, and the ss102661460 in functional target sites for two miRNAs, miR-509 and miR-128, the latter being a brain-enriched miRNA involved in neuronal differentiation and synaptic processing. Interestingly, these two variants significantly alter the miRNA-mediated regulation of NTRK3, resulting in recovery of gene expression. These data implicate miRNAs as key posttranscriptional regulators of NTRK3 and provide a framework for allele-specific miRNA regulation of NTRK3 in anxiety disorders.

Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility

Background: Germline genetic variation is associated with the differential expression of many human genes. The phenotypic effects of this type of variation may be important when considering susceptibility to common genetic diseases. Three regions at 8q24 have recently been identified to independently confer risk of prostate cancer. Variation at 8q24 has also recently been associated with risk of breast and colorectal cancer. However, none of the risk variants map at or relatively close to known genes, with c-MYC mapping a few hundred kilobases distally. Results: This study identifies cis-regulators of germline c-MYC expression in immortalized lymphocytes of HapMap individuals. Quantitative analysis of c-MYC expression in normal prostate tissues suggests an association between overexpression and variants in Region 1 of prostate cancer risk. Somatic c-MYC overexpression correlates with prostate cancer progression and more aggressive tumor forms, which was also a pathological variable associated with Region 1. Expression profiling analysis and modeling of transcriptional regulatory networks predicts a functional association between MYC and the prostate tumor suppressor KLF6. Analysis of MYC/Myc-driven cell transformation and tumorigenesis substantiates a model in which MYC overexpression promotes transformation by down-regulating KLF6. In this model, a feedback loop through E-cadherin down-regulation causes further transactivation of c-MYC.Conclusion: This study proposes that variation at putative 8q24 cis-regulator(s) of transcription can significantly alter germline c-MYC expression levels and, thus, contribute to prostate cancer susceptibility by down-regulating the prostate tumor suppressor KLF6 gene.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Tumor necrosis factor and transforming growth factor β regulate clock genes by controlling the expression of the cold inducible RNA-binding protein (CIRBP).

The circadian clock drives the rhythmic expression of a broad array of genes that orchestrate metabolism, sleep wake behavior, and the immune response. Clock genes are transcriptional regulators engaged in the generation of circadian rhythms. The cold inducible RNA-binding protein (CIRBP) guarantees high amplitude expression of clock. The cytokines TNF and TGFβ impair the expression of clock genes, namely the period genes and the proline- and acidic amino acid-rich basic leucine zipper (PAR-bZip) clock-controlled genes. Here, we show that TNF and TGFβ impair the expression of Cirbp in fibroblasts and neuronal cells. IL-1β, IL-6, IFNα, and IFNγ do not exert such effects. Depletion of Cirbp is found to increase the susceptibility of cells to the TNF-mediated inhibition of high amplitude expression of clock genes and modulates the TNF-induced cytokine response. Our findings reveal a new mechanism of cytokine-regulated expression of clock genes.

Characterization of the role of AKAP-Lbc/p38 complex in the hypertrophic response of the heart

In response to chronic stress the heart undergoes an adverse remodeling process associated with cardiomyocyte hypertrophy, increased cellular apoptosis and fibrosis, which ultimately causes cardiac dysfunction and heart failure. Increasing evidence suggest the role of scaffolding and anchoring proteins in coordinating different signaling pathways that mediate the hypertrophic response of the heart. In this context, the family of Α-kinase anchoring proteins (AKAPs) emerged as important regulators of the cardiac function. During my thesis work I have conducted two independent projects, both of them aiming at elucidating the role of AKAPs in the heart. It has been shown that AKAP-Lbc, an anchoring protein that possesses an intrinsic Rho- specific exchange factor activity, organizes a signaling complex that links AKAP-Lbc- dependent activation of RhoA with the mitogen activated protein kinase (MAPK) p38. The first aim of my thesis was to study the role of this novel transduction pathway in the context of cardiac hypertrophy. Here we show that transgenic mice overexpressing in cardiomyocytes a competitor fragment of AKAP-Lbc, which specifically disrupts endogenous AKAP-Lbc / p38 complexes, developed early dilated cardiomyopathy in response to two weeks of transverse aortic constriction (TAC) as compared to controls. Interestingly, inhibition of the AKAP-Lbc / p38 transduction pathway significantly reduced the hypertrophic growth of single cardiomyocytes induced by pressure overload. Therefore, it appears that the AKAP- Lbc / p38 complex is crucially involved in the regulation of stress-induced cardiomyocyte hypertrophy and that disruption of this signaling pathway is detrimental for the heart under conditions of sustained hemodynamic stress. Secondly, in order to identify new AKAPs involved in the regulation of cardiac function, we followed a proteomic approach which allowed us to characterize AKAP2 as a major AKAP in the heart. Importantly, here we show that AKAP2 interacts with several proteins known to be involved in the control of gene transcription, such as the nuclear receptor coactivator 3 (NCoA3) or the ATP-dependent SWI/SNF chromatin remodeling complex. Thus, we propose AKAP2 as a novel mediator of cardiac gene expression through its interaction with these transcriptional regulators.

Effects of epigenetic regulatory elements on telomeric gene expression

Transgene expression in eukaryotic cells strongly depends on the locus of integration in the host genome and often results in limited transcription level because of unfavorable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which are believed to act on the chromatin structure and may protect transgenes from this so-called position effect. Despite being extensively used to increase transgene expression, the mechanism of action of many of these elements remains largely unknown. Here we evaluated different epigenetic regulatory DNA elements for their ability to protect transgene transcription at telomeres, a defined chromatin environment associated to low or inconsistent expression caused by the Telomere Position Effect (TPE). For the assessment of the effects of epigenetic regulators at telomeres, a novel dual reporter system had to be designed. Telomeric integration of the newly-developed dual reporter system carrying different epigenetic regulators showed that MARs (Matrix Attachment Regions), a UCOE (Ubiquitous Chromatin-Opening Element) or the chicken cHS4 insulator have strong barrier activity which prevented TPE from spreading toward the centromere, resulting in stable and in some cases increased expression of a telomeric-distal reporter gene. In addition, MARs and STAR element 40 resulted in an increase of cells expressing the telomeric-proximal reporter gene, suggesting also an anti-silencing effect. Chromatin immunoprecipitation assays revealed that at telomeres MARs actively promote the deposition of euchromatic histone marks, especially acetylation of both histone H3 and H4, which might be involved in MARs' barrier and transcriptional activator activities. Differently, the chromatin in proximity of the UCOE element was depleted of several repressive chromatin marks, such as trimethylated lysine 9 and lysine 27 on histone H3 and trimethylated lysine 20 of histone H4, possibly favoring the preservation of an open chromatin structure at the integration site. We conclude that epigenetic regulatory elements that may be used to enhance and sustain transgene expression have all a specific epigenetic signature which might be at the basis of their mechanism of action, and that a combination of different classes of epigenetic regulators might be advantageous when high levels of protein expression are required. - L'expression des transgènes dans les cellules eucaryotes est fortement influencée par leur site d'intégration dans le génome. En effet, une structure chromatinienne défavorable au niveau du site d'intégration peut fortement limiter le degré d'expression d'un transgène. Il existe toutefois des séquences d'ADN qui, en agissant sur la structure de la chromatine, permettent de limiter cet effet de position et, par conséquent, de promouvoir l'expression soutenue d'un transgène. Ces éléments génomiques, connus comme régulateurs épigénétiques, sont largement utilisés dans plusieurs domaines où une expression élevée et soutenue est requise, malgré un mode de fonctionnement parfois méconnu. Dans cette étude, j'ai évalué la capacité de différents régulateurs épigénétiques à protéger la transcription de transgènes au niveau des télomères, régions chromatiniennes bien définies qui ont été associées à un fort effet de silençage, connu comme «effet de position télomérique». Pour cela, un nouveau système à deux gènes rapporteurs a été développé. Lorsque des MARs (Matrix Attachment Regions, séquences d'ADN pouvant s'associer à la matrice nucléaire), un UCOE (Ubiquitous Chromatin-Opening Element, élément pouvant ouvrir la chromatine) ou l'isolateur génétique cHS4 (dérivé du locus de la β-globine de poulet) sont placés entre les deux gènes rapporteurs, une forte activité barrière bloquant la propagation de la chromatine répressive télomérique est observée, résultant en un plus grand nombre de cellules exprimant le gène télomérique-distal. D'autre part, une augmentation du nombre de cellules exprimant le gène télomérique-proximal, observée en présence des éléments MAR et STAR 40 (Stabilizing Anti-Repressor element 40, un élément pouvant prévenir la répression génique), suggère aussi un faible effet anti-silençage pour ces éléments. Des expériences d'immunoprécipitation de la chromatine démontrent qu'au télomère, les MARs favorisent l'assemblage de marqueurs de la chromatine active, surtout l'acétylation des histones H3 et H4, qui pourraient être à la base de l'activité barrière et de celle d'activateur transcriptionel. Différemment, la chromatine à proximité de l'élément UCOE est particulièrement pauvre en marqueurs de la chromatine silencieuse, comme la trimethylation des lysines 9 et 27 de l'histone H3, ainsi que la trimethylation de la lysine 20 de l'histone H4. Cela suggère que UCOE pourrait préserver une structure chromatinienne ouverte au site d'intégration, favorisant l'expression des gènes à sa proximité. En conclusion, les régulateurs épigénétiques analysés lors de cette étude ont tous montré une signature épigénétique spécifique qui pourrait être à la base de leurs mécanismes de fonctionnement, suggérant aussi qu'une utilisation d'éléments épigénétiques de classe différente dans un même vecteur d'expression pourrait être avantageuse lorsque de hauts et soutenus niveaux d'expression sont nécessaires.

Sex and clock to discover new PPARα functions

Summary : PPARα is a ligand-activated transcription factor that is a member of the nuclear receptor superfamily. In rodents, PPARα is highly expressed in liver, especially in parenchymal cells, where it has an impact on several hepatic functions such as nutrient metabolism, inflammation and metabolic stress. Ligands for PPARα comprise long chain unsaturated fatty acids, eicosanoids and lipid lowering fibrate drugs. In liver, many metabolic processes are orchestrated by the hepatic circadian clock. The aim of the hepatic clock is to synchronize cellular pathways allowing animals to adapt their metabolism to predictable daily changes in the environment. Indeed, similar to PPARα, the hepatic clock influences nutrient metabolism and detoxification through circadian output regulators :the PAR-domain basic leucine zipper proteins called PAR blip proteins. In this report, we showed that through a positive feedback loop mechanism, PAR. blip, proteins participate to the availability of PPARα endogenous ligands that contribute to the circadian expression and functions of PPARα. Interestingly, we also discovered some unexpected hepatic sexual dimorphic functions of PPARα. These functions are determined b PPARα sumoylation, interaction with DNA methylation mechanism and with unexpected proteins with gender specificity. The connection between circadian clock and hepatic sexual dimorphism opens new perspectives regarding the chronobiology of PPARα activity and the beneficial effects of PPARα agonist in the treatment of diseases related to steroid hormones metabolism characterized by inflammation and hepatotoxicity. Résumé : PPARα est un facteur de transcription activé par un ligand, membre de la superfamille des récepteurs nucléaires. Chez les rongeurs, PPARα est fortement exprimé dans le foie, spécialement dans les cellules du parenchyme dans lesquelles il joue un role important dans les fonctions hépatiques tels que le métabolisme des nutriments, l'inflammation et les stress métaboliques. Les ligands pour PPARα comprennent les acides gras à longues chaînes, les eicosanoides et les médicaments hypolipidémiques (fibrates). Dans le foie, beaucoup de processus métaboliques sont orchestrés par l'horloge circadienne hépatique. Le but de cette horloge est de synchroniser les voies métaboliqués permettant aux animaux d'adapter leurs métabolismes aux changements journaliers. Ainsi, l'horloge hépatique influence le métabolisme des nutriments tels que l'utilisation des lipides à travers certains régulateurs circadians appelés facteurs de transcription PAR bZips. Dans ce mémoire, nous avons montré qu'à travers une boucle de régulation, les protéines PAR bZip contrôlent la production des ligands endogènes à PPARα, jouant un rôle dans l'expression circadienne et les fonctions de PPARα. Nous avons également découvert des aspects méconnus des fonctions liées au dimorphisme sexuel de PPARα. Nous avons montré que PPARα est différemment sumoylisé entre les sexes et interagit avec la méthylation de l'ADN ainsi qu'avec des protéines insoupçonnées comme partenaires de PPARα. De part leur lien avec l'horloge circadienne et le dimorphisme sexuel, nos découvertes ouvrent de nouvelles perspectives concernant la chronobiologie de l'activité de PPARα et les effets bénéfiques des ses activateurs dans le traitement des maladies liées au métabolisme des hormones stéroides.

Stereospecific recognition of pyochelin and enantio-pyochelin by the PchR proteins in fluorescent pseudomonads.

The siderophore pyochelin of Pseudomonas aeruginosa promotes growth under iron limitation and induces the expression of its biosynthesis genes via the transcriptional AraC/XylS-type regulator PchR. Pseudomonas fluorescens strain CHA0 makes the optical antipode of pyochelin termed enantio-pyochelin, which also promotes growth and induces the expression of its biosynthesis genes when iron is scarce. Growth promotion and signalling by pyochelin and enantio-pyochelin are highly stereospecific and are known to involve the pyochelin and enantio-pyochelin outer-membrane receptors FptA and FetA, respectively. Here we show that stereospecificity in signalling is also based on the stereospecificity of the homologous PchR proteins of P. aeruginosa and P. fluorescens towards their respective siderophore effectors. We found that PchR functioned in the heterologous species only if supplied with its native ligand and that the FptA and FetA receptors enhanced the efficiency of signalling. By constructing and expressing hybrid and truncated PchR regulators we showed that the weakly conserved N-terminal domain of PchR is responsible for siderophore specificity. Thus, both uptake and transcriptional regulation confer stereospecificity to pyochelin and enantio-pyochelin biosynthesis.

Patterns of plant subcellular responses to successful oomycete infections reveal differences in host cell reprogramming and endocytic trafficking.

Adapted filamentous pathogens such as the oomycetes Hyaloperonospora arabidopsidis (Hpa) and Phytophthora infestans (Pi) project specialized hyphae, the haustoria, inside living host cells for the suppression of host defence and acquisition of nutrients. Accommodation of haustoria requires reorganization of the host cell and the biogenesis of a novel host cell membrane, the extrahaustorial membrane (EHM), which envelops the haustorium separating the host cell from the pathogen. Here, we applied live-cell imaging of fluorescent-tagged proteins labelling a variety of membrane compartments and investigated the subcellular changes associated with accommodating oomycete haustoria in Arabidopsis and N. benthamiana. Plasma membrane-resident proteins differentially localized to the EHM. Likewise, secretory vesicles and endosomal compartments surrounded Hpa and Pi haustoria revealing differences between these two oomycetes, and suggesting a role for vesicle trafficking pathways for the pathogen-controlled biogenesis of the EHM. The latter is supported by enhanced susceptibility of mutants in endosome-mediated trafficking regulators. These observations point at host subcellular defences and specialization of the EHM in a pathogen-specific manner. Defence-associated haustorial encasements, a double-layered membrane that grows around mature haustoria, were frequently observed in Hpa interactions. Intriguingly, all tested plant proteins accumulated at Hpa haustorial encasements suggesting the general recruitment of default vesicle trafficking pathways to defend pathogen access. Altogether, our results show common requirements of subcellular changes associated with oomycete biotrophy, and highlight differences between two oomycete pathogens in reprogramming host cell vesicle trafficking for haustoria accommodation. This provides a framework for further dissection of the pathogen-triggered reprogramming of host subcellular changes.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Role of islet microRNAs in diabetes: which model for which question?

MicroRNAs are important regulators of gene expression. The vast majority of the cells in our body rely on hundreds of these tiny non-coding RNA molecules to precisely adjust their protein repertoire and faithfully accomplish their tasks. Indeed, alterations in the microRNA profile can lead to cellular dysfunction that favours the appearance of several diseases. A specific set of microRNAs plays a crucial role in pancreatic beta cell differentiation and is essential for the fine-tuning of insulin secretion and for compensatory beta cell mass expansion in response to insulin resistance. Recently, several independent studies reported alterations in microRNA levels in the islets of animal models of diabetes and in islets isolated from diabetic patients. Surprisingly, many of the changes in microRNA expression observed in animal models of diabetes were not detected in the islets of diabetic patients and vice versa. These findings are unlikely to merely reflect species differences because microRNAs are highly conserved in mammals. These puzzling results are most probably explained by fundamental differences in the experimental approaches which selectively highlight the microRNAs directly contributing to diabetes development, the microRNAs predisposing individuals to the disease or the microRNAs displaying expression changes subsequent to the development of diabetes. In this review we will highlight the suitability of the different models for addressing each of these questions and propose future strategies that should allow us to obtain a better understanding of the contribution of microRNAs to the development of diabetes mellitus in humans.

The endoplasmic reticulum: a sensor of cellular stress that modulates immune responses.

Many inflammatory and infectious diseases are characterized by the activation of signaling pathways steaming from the endoplasmic reticulum (ER). These pathways, primarily associated with loss of ER homeostasis, are emerging as key regulators of inflammation and infection. Recent advances shed light on the mechanisms linking ER-stress and immune responses.

Study of protein complexes

Summary : A lot of information can be obtained on proteins when proteomics methods are used. In our study, we aimed to characterize complexes containing pro-apoptotic proteins by different proteomics methods and finally focused on PIDD (p53-induced protein with a death domain), for which the most interesting results were obtained. PIDD has been shown to function as a molecular switch between genotoxic stress-induced apoptotis and genotoxic stress-induced cell survival through NF-κB activation. To exert these two functions, PIDD forms alternate complexes respectively with caspase2 and CRADD on one hand and RIP 1 and NEMO on the other hand. The first part of our study focuses on the processing of PIDD. PIDD full length (FL) is constitutively cleaved into three fragments, an N-terminal one (PIDD-N) and two fragments containing the C-terminus (PIDD-C and PIDD-CC). Localization of the two PIDD cleavage sites by mass spectrometry (MS) allowed to understand that PIDD is probably not cleaved by proteases but is subject to protein (self-)splicing and also to map the PIDD-N, PIDD-C and PIDD-CC fragments exactly. Further characterization of these three fragments by Tinel et al. (Tinel et al., 2007) showed that PIDD-C is involved in activation of an apoptotic pathway while PIDD-CC is involved in NF-κB activation. We also found that PIDD is subject to proline-directed phosphorylation at two serine residues in PIDD-N, the regulatory fragment of PIDD. The second part of the study aimed at identifying by proteomics techniques proteins that co-purify with PIDD and therefore are putative cellular interaction partners. In this respect we analyzed samples obtained in different conditions or with different PIDD constructs corresponding to processed fragments. This allowed us to identify a large number of potential interactors for PIDD. For example, by comparing data obtained from PIDD-C and PIDD-FL affinity purifications, we found that the Hsp90 chaperone system interacts strongly with PIDD-N. In the third part of this study, we developed methods to selectively and rapidly quantify by MS proteins of interest in PIDD affinity purifications or negative controls. Using these tools we detected significant changes in PIDD-FL-copurifying proteins treated by heat shock. Overall, our studies provide informative data on the processing of PIDD and its possible involvement in several molecular pathways.

Promises of apoptosis-inducing peptides in cancer therapeutics.

Until recently, most research efforts aimed at developing anti-cancer tools were focusing on small molecules. Alternative compounds are now being increasingly assessed for their potential anti-cancer properties, including peptides and their derivatives. One earlier limitation to the use of peptides was their limited capacity to cross membranes but this limitation was alleviated with the characterization of cell-permeable sequences. Additionally, means are designed to target peptides to their malignant targets. Most anti-cancer peptidic compounds induce apoptosis of tumor cells by modulating the activity of Bcl-2 family members that control the release of death factors from the mitochondria or by inhibiting negative regulators of caspases, the proteases that mediate the apoptotic response in cells. Some of these peptides have been shown to inhibit the growth of tumors in mouse models. Hopefully, pro-apoptotic anti-tumor peptides will soon be tested for their efficacy in patients with cancers.