Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. Results: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.

Identification of novel functional TBP-binding sites and general factor repertoires.

Our current knowledge of the general factor requirement in transcription by the three mammalian RNA polymerases is based on a small number of model promoters. Here, we present a comprehensive chromatin immunoprecipitation (ChIP)-on-chip analysis for 28 transcription factors on a large set of known and novel TATA-binding protein (TBP)-binding sites experimentally identified via ChIP cloning. A large fraction of identified TBP-binding sites is located in introns or lacks a gene/mRNA annotation and is found to direct transcription. Integrated analysis of the ChIP-on-chip data and functional studies revealed that TAF12 hitherto regarded as RNA polymerase II (RNAP II)-specific was found to be also involved in RNAP I transcription. Distinct profiles for general transcription factors and TAF-containing complexes were uncovered for RNAP II promoters located in CpG and non-CpG islands suggesting distinct transcription initiation pathways. Our study broadens the spectrum of general transcription factor function and uncovers a plethora of novel, functional TBP-binding sites in the human genome.

Functional analysis of the 5' flanking sequence of a vaccinia virus late gene.

A series of mutations, including 5' and 3' deletions, as well as insertions were introduced into the 5' flanking nucleotide sequence of a vaccinia virus late gene. This DNA has been shown previously to contain all the necessary elements for correct regulation of the gene most probably transcribed by the viral RNA polymerase. To facilitate the assays, the mutated DNA was fused to the chloramphenicol acetyltransferase gene and inserted into the genome of live vaccinia virus. The effects of the mutations on expression of the chimeric gene were studied by both enzyme assays and nuclease S1 analysis. The results showed that 5' deletions up to about 15 bp from the putative initiation site of transcription still yielded high levels of gene expression. All mutations, however, that deleted the authentic late mRNA start site, abolished promoter activity.

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.

The GENCODE pseudogene resource.

BACKGROUND: Pseudogenes have long been considered as nonfunctional genomic sequences. However, recent evidence suggests that many of them might have some form of biological activity, and the possibility of functionality has increased interest in their accurate annotation and integration with functional genomics data. RESULTS: As part of the GENCODE annotation of the human genome, we present the first genome-wide pseudogene assignment for protein-coding genes, based on both large-scale manual annotation and in silico pipelines. A key aspect of this coupled approach is that it allows us to identify pseudogenes in an unbiased fashion as well as untangle complex events through manual evaluation. We integrate the pseudogene annotations with the extensive ENCODE functional genomics information. In particular, we determine the expression level, transcription-factor and RNA polymerase II binding, and chromatin marks associated with each pseudogene. Based on their distribution, we develop simple statistical models for each type of activity, which we validate with large-scale RT-PCR-Seq experiments. Finally, we compare our pseudogenes with conservation and variation data from primate alignments and the 1000 Genomes project, producing lists of pseudogenes potentially under selection. CONCLUSIONS: At one extreme, some pseudogenes possess conventional characteristics of functionality; these may represent genes that have recently died. On the other hand, we find interesting patterns of partial activity, which may suggest that dead genes are being resurrected as functioning non-coding RNAs. The activity data of each pseudogene are stored in an associated resource, psiDR, which will be useful for the initial identification of potentially functional pseudogenes.

RpoN (sigma54) controls production of antifungal compounds and biocontrol activity in Pseudomonas fluorescens CHA0.

Pseudomonas fluorescens CHA0 is an effective biocontrol agent of root diseases caused by fungal pathogens. The strain produces the antibiotics 2,4-diacetylphloroglucinol (DAPG) and pyoluteorin (PLT) that make essential contributions to pathogen suppression. This study focused on the role of the sigma factor RpoN (sigma54) in regulation of antibiotic production and biocontrol activity in P. fluorescens. An rpoN in-frame-deletion mutant of CHAO had a delayed growth, was impaired in the utilization of several carbon and nitrogen sources, and was more sensitive to salt stress. The rpoN mutant was defective for flagella and displayed drastically reduced swimming and swarming motilities. Interestingly, the rpoN mutant showed a severalfold enhanced production of DAPG and expression of the biosynthetic gene phlA compared with the wild type and the mutant complemented with monocopy rpoN+. By contrast, loss of RpoN function resulted in markedly lowered PLT production and plt gene expression, suggesting that RpoN controls the balance of the two antibiotics in strain CHA0. In natural soil microcosms, the rpoN mutant was less effective in protecting cucumber from a root rot caused by Pythium ultimum. Remarkably, the mutant was not significantly impaired in its root colonization capacity, even at early stages of root infection by Pythium spp. Taken together, our results establish RpoN for the first time as a major regulator of biocontrol activity in Pseudomonas fluorescens.

mTORC1 Directly Phosphorylates and Regulates Human MAF1.

mTORC1 is a central regulator of growth in response to nutrient availability, but few direct targets have been identified. RNA polymerase (pol) III produces a number of essential RNA molecules involved in protein synthesis, RNA maturation, and other processes. Its activity is highly regulated, and deregulation can lead to cell transformation. The human phosphoprotein MAF1 becomes dephosphorylated and represses pol III transcription after various stresses, but neither the significance of the phosphorylations nor the kinase involved is known. We find that human MAF1 is absolutely required for pol III repression in response to serum starvation or TORC1 inhibition by rapamycin or Torin1. The protein is phosphorylated mainly on residues S60, S68, and S75, and this inhibits its pol III repression function. The responsible kinase is mTORC1, which phosphorylates MAF1 directly. Our results describe molecular mechanisms by which mTORC1 controls human MAF1, a key repressor of RNA polymerase III transcription, and add a new branch to the signal transduction cascade immediately downstream of TORC1.

Molecular characterization of a human matrix attachment region epigenetic regulator.

Matrix attachment regions (MAR) generally act as epigenetic regulatory sequences that increase gene expression, and they were proposed to partition chromosomes into loop-forming domains. However, their molecular mode of action remains poorly understood. Here, we assessed the possible contribution of the AT-rich core and adjacent transcription factor binding motifs to the transcription augmenting and anti-silencing effects of human MAR 1-68. Either flanking sequences together with the AT-rich core were required to obtain the full MAR effects. Shortened MAR derivatives retaining full MAR activity were constructed from combinations of the AT-rich sequence and multimerized transcription factor binding motifs, implying that both transcription factors and the AT-rich microsatellite sequence are required to mediate the MAR effect. Genomic analysis indicated that MAR AT-rich cores may be depleted of histones and enriched in RNA polymerase II, providing a molecular interpretation of their chromatin domain insulator and transcriptional augmentation activities.

Identification of genes potentially involved in solute stress response in Sphingomonas wittichii RW1 by transposon mutant recovery.

The term water stress refers to the effects of low water availability on microbial growth and physiology. Water availability has been proposed as a major constraint for the use of microorganisms in contaminated sites with the purpose of bioremediation. Sphingomonas wittichii RW1 is a bacterium capable of degrading the xenobiotic compounds dibenzofuran and dibenzo-p-dioxin, and has potential to be used for targeted bioremediation. The aim of the current work was to identify genes implicated in water stress in RW1 by means of transposon mutagenesis and mutant growth experiments. Conditions of low water potential were mimicked by adding NaCl to the growth media. Three different mutant selection or separation method were tested which, however recovered different mutants. Recovered transposon mutants with poorer growth under salt-induced water stress carried insertions in genes involved in proline and glutamate biosynthesis, and further in a gene putatively involved in aromatic compound catabolism. Transposon mutants growing poorer on medium with lowered water potential also included ones that had insertions in genes involved in more general functions such as transcriptional regulation, elongation factor, cell division protein, RNA polymerase β or an aconitase.

Functional importance of cardiac enhancer-associated noncoding RNAs in heart development and disease.

The key information processing units within gene regulatory networks are enhancers. Enhancer activity is associated with the production of tissue-specific noncoding RNAs, yet the existence of such transcripts during cardiac development has not been established. Using an integrated genomic approach, we demonstrate that fetal cardiac enhancers generate long noncoding RNAs (lncRNAs) during cardiac differentiation and morphogenesis. Enhancer expression correlates with the emergence of active enhancer chromatin states, the initiation of RNA polymerase II at enhancer loci and expression of target genes. Orthologous human sequences are also transcribed in fetal human hearts and cardiac progenitor cells. Through a systematic bioinformatic analysis, we identified and characterized, for the first time, a catalog of lncRNAs that are expressed during embryonic stem cell differentiation into cardiomyocytes and associated with active cardiac enhancer sequences. RNA-sequencing demonstrates that many of these transcripts are polyadenylated, multi-exonic long noncoding RNAs. Moreover, knockdown of two enhancer-associated lncRNAs resulted in the specific downregulation of their predicted target genes. Interestingly, the reactivation of the fetal gene program, a hallmark of the stress response in the adult heart, is accompanied by increased expression of fetal cardiac enhancer transcripts. Altogether, these findings demonstrate that the activity of cardiac enhancers and expression of their target genes are associated with the production of enhancer-derived lncRNAs.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

Characterization of a MAF1 knockout mouse model

L'ARN polymérase 3 transcrit un petit groupe de gènes fortement exprimés et impliqués dans plusieurs mécanismes moléculaires. Les ARNs de transfert ou ARNt représentent plus ou moins la moitié du transcriptome de l'ARN polymérase 3. Ils sont directement impliqués dans la traduction des protéines en agissant comme transporteurs d'acides aminés qui sont incorporés à la chaîne naissante de polypeptides. Chez des levures cultivées dans un milieu jusqu'à épuisement des nutriments, Maf1 réprime la transcription par l'ARN polymérase 3, favorisant ainsi l'économie énergétique cellulaire. Dans un modèle de cellules de mammifères, MAF1 réprime aussi la transcription de l'ARN polymérase 3 dans des conditions de stress, cependant il n'existe aucune donnée quant à son rôle chez un mammifère vivant. Pendant mon doctorat, j'ai utilisé une souris délétée pour le gène Maf1 afin de connaître les effets de ce gène chez un mammifère. Etonnamment, la souris Maf1-­‐/-­‐ est résistante à l'obésité même si celle-­‐ci est nourrie avec une nourriture riche en matières grasses. Des études moléculaires et de métabolomiques ont montré qu'il existe des cycles futiles de production et dégradation des lipides et des ARNt, ce qui entraîne une augmentation de la dépense énergique et favorise la résistance à l'obésité. En plus de la caractérisation de la souris Maf1-­‐/-­‐, pendant ma thèse j'ai également développé une méthode afin de normaliser les données de ChIP-­‐sequencing. Cette méthode est fondée sur l'utilisation d'un contrôle interne, représenté ici par l'ajout d'une quantité fixe de chromatine provenant d'un organisme différent de celui étudié. La méthode a amélioré considérablement la reproductibilité des valeurs entre réplicas biologiques. Elle a aussi révélé des différences entre échantillons issus de conditions différentes. Une occupation supérieure de l'ARN polymérase 3 sur les gènes Pol 3 chez les souris Maf1 KO entraîne une augmentation du niveau de précurseurs d'ARNt, ayant pour effet probable la saturation de la machinerie de maturation des ARNt. En effet, chez les souris Maf1 KO, le pourcentage d'ARNt modifiés est plus faible que chez les souris type sauvage. Ce déséquilibre entre le niveau de précurseurs et d'ARNt matures entraîne une diminution de la traduction protéique. Ces résultats ont permis d'identifier de nouvelles fonctions pour la protéine MAF1, comme étant une protéine régulatrice à la fois de la transcription mais aussi de la traduction et en étant un cible potentielle au traitement à l'obésité. -- RNA polymerase III (Pol 3) transcribes a small set of highly expressed genes involved in different molecular mechanisms. tRNAs account for almost half of the Pol 3 transcriptome and are involved in translation, bringing a new amino into the nascent polypeptide chain. In yeast, under nutrient deprivation, Maf1 acts for cell energetic economy by repressing Pol 3 transcription. In mammalian cells, MAF1 also represses Pol 3 activity under conditions of serum deprivation or DNA damages but nothing is known about its role in a mammalian organism. During my thesis studies, I used a Maf1 KO mouse model to characterize the effects of Maf1 deletion in a living animal. Surprisingly, the MAF1 KO mouse developed an unexpected phenotype, being resistant to high fat diet-­‐induced obesity and displaying an extended lifespan. Molecular and metabolomics characterizations revealed futile cycles of lipids and tRNAs, which are produced and immediately degraded, which increases energy consumption in the Maf1 KO mouse and probably explains in part the protection to obesity. Additionally to the mouse characterization, I also developed a method to normalize ChIP-­‐seq data, based on the addition of a foreign chromatin to be used as an internal control. The method improved reproducibility between replicates and revealed differences of Pol 3 occupancy between WT and Maf1 KO samples that were not seen without normalization to the internal control. I then established that increased Pol 3 occupancy in the Maf1 KO mouse liver was associated with increased levels of tRNA precursor but not of mature tRNAs, the effective molecules involved in translation. The overproduction of precursor tRNAs associated with the deletion of Maf1 apparently overwhelms the tRNA processing machinery as the Maf1 KO mice have lower levels of fully modified tRNAs. This maturation defect directly impacts on translation efficiency as polysomic fractions and newly synthetized protein levels were reduced in the liver of the Maf1 KO mouse. Altogether, these results indicate new functions for MAF1, a regulator of both transcription and translation as well as a potential target for obesity treatment.

RNA-driven cyclin-dependent kinase regulation: when CDK9/cyclin T subunits of P-TEFb meet their ribonucleoprotein partners.

The positive transcription elongation factor (P-TEFb) consists of CDK9, a cyclin-dependent kinase and its cyclin T partner. It is required for transcription of most class II genes. Its activity is regulated by non-coding RNAs. The 7SK cellular RNA turns the HEXIM cellular protein into a P-TEFb inhibitor that binds its cyclin T subunit. Thus, P-TEFb activity responds to variations in global cellular transcriptional activity and to physiological conditions linked to cell differentiation, proliferation or cardiac hypertrophy. In contrast, the Tat activation region RNA plays an activating role. This feature at the 5' end of the human immunodeficiency (HIV) viral transcript associates with the viral protein Tat that in turn binds cyclin T1 and recruits active P-TEFb to the HIV promoter. This results in enhanced P-TEFb activity, which is critical for an efficient production of viral transcripts. Although discovered recently, the regulation of P-TEFb becomes a paradigm for non-coding RNAs that regulate transcription factors. It is also a unique example of RNA-driven regulation of a cyclindependent kinase.

NrsZ: a novel, processed, nitrogen-dependent, small non-coding RNA that regulates Pseudomonas aeruginosa PAO1 virulence.

The opportunistic pathogen Pseudomonas aeruginosa PAO1 has a remarkable capacity to adapt to various environments and to survive with limited nutrients. Here, we report the discovery and characterization of a novel small non-coding RNA: NrsZ (nitrogen-regulated sRNA). We show that under nitrogen limitation, NrsZ is induced by the NtrB/C two component system, an important regulator of nitrogen assimilation and P. aeruginosa's swarming motility, in concert with the alternative sigma factor RpoN. Furthermore, we demonstrate that NrsZ modulates P. aeruginosa motility by controlling the production of rhamnolipid surfactants, virulence factors notably needed for swarming motility. This regulation takes place through the post-transcriptional control of rhlA, a gene essential for rhamnolipids synthesis. Interestingly, we also observed that NrsZ is processed in three similar short modules, and that the first short module encompassing the first 60 nucleotides is sufficient for NrsZ regulatory functions.

A TARBP2-dependent miRNA expression profile underlies cancer stem cell properties and provides candidate therapeutic reagents in Ewing sarcoma.

We have recently demonstrated that human pediatric mesenchymal stem cells can be reprogrammed toward a Ewing sarcoma family tumor (ESFT) cancer stem cell (CSC) phenotype by mechanisms that implicate microRNAs (miRNAs). Here, we show that the miRNA profile of ESFT CSCs is shared by embryonic stem cells and CSCs from divergent tumor types. We also provide evidence that the miRNA profile of ESFT CSCs is the result of reversible disruption of TARBP2-dependent miRNA maturation. Restoration of TARBP2 activity and systemic delivery of synthetic forms of either of two of its targets, miRNA-143 or miRNA-145, inhibited ESFT CSC clonogenicity and tumor growth in vivo. Our observations suggest that CSC self-renewal and tumor maintenance may depend on deregulation of TARBP2-dependent miRNA expression.