The Translational Landscape of the Human Heart.

Gene expression in human tissue has primarily been studied on the transcriptional level, largely neglecting translational regulation. Here, we analyze the translatomes of 80 human hearts to identify new translation events and quantify the effect of translational regulation. We show extensive translational control of cardiac gene expression, which is orchestrated in a process-specific manner. Translation downstream of predicted disease-causing protein-truncating variants appears to be frequent, suggesting inefficient translation termination. We identify hundreds of previously undetected microproteins, expressed from lncRNAs and circRNAs, for which we validate the protein products in vivo. The translation of microproteins is not restricted to the heart and prominent in the translatomes of human kidney and liver. We associate these microproteins with diverse cellular processes and compartments and find that many locate to the mitochondria. Importantly, dozens of microproteins are translated from lncRNAs with well-characterized noncoding functions, indicating previously unrecognized biology.

Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

Natural variation of histone modification and its impact on gene expression in the rat genome.

Histone modifications are epigenetic marks that play fundamental roles in many biological processes including the control of chromatin-mediated regulation of gene expression. Little is known about interindividual variability of histone modification levels across the genome and to what extent they are influenced by genetic variation. We annotated the rat genome with histone modification maps, identified differences in histone trimethyl-lysine levels among strains, and described their underlying genetic basis at the genome-wide scale using ChIP-seq in heart and liver tissues in a panel of rat recombinant inbred and their progenitor strains. We identified extensive variation of histone methylation levels among individuals and mapped hundreds of underlying cis- and trans-acting loci throughout the genome that regulate histone methylation levels in an allele-specific manner. Interestingly, most histone methylation level variation was trans-linked and the most prominent QTL identified influenced H3K4me3 levels at 899 putative promoters throughout the genome in the heart. Cis- acting variation was enriched in binding sites of distinct transcription factors in heart and liver. The integrated analysis of DNA variation together with histone methylation and gene expression levels showed that histoneQTLs are an important predictor of gene expression and that a joint analysis significantly enhanced the prediction of gene expression traits (eQTLs). Our data suggest that genetic variation has a widespread impact on histone trimethylation marks that may help to uncover novel genotype-phenotype relationships.

A trans-acting locus regulates an anti-viral expression network and type 1 diabetes risk.

Combined analyses of gene networks and DNA sequence variation can provide new insights into the aetiology of common diseases that may not be apparent from genome-wide association studies alone. Recent advances in rat genomics are facilitating systems-genetics approaches. Here we report the use of integrated genome-wide approaches across seven rat tissues to identify gene networks and the loci underlying their regulation. We defined an interferon regulatory factor 7 (IRF7)-driven inflammatory network (IDIN) enriched for viral response genes, which represents a molecular biomarker for macrophages and which was regulated in multiple tissues by a locus on rat chromosome 15q25. We show that Epstein-Barr virus induced gene 2 (Ebi2, also known as Gpr183), which lies at this locus and controls B lymphocyte migration, is expressed in macrophages and regulates the IDIN. The human orthologous locus on chromosome 13q32 controlled the human equivalent of the IDIN, which was conserved in monocytes. IDIN genes were more likely to associate with susceptibility to type 1 diabetes (T1D)-a macrophage-associated autoimmune disease-than randomly selected immune response genes (P = 8.85 × 10(-6)). The human locus controlling the IDIN was associated with the risk of T1D at single nucleotide polymorphism rs9585056 (P = 7.0 × 10(-10); odds ratio, 1.15), which was one of five single nucleotide polymorphisms in this region associated with EBI2 (GPR183) expression. These data implicate IRF7 network genes and their regulatory locus in the pathogenesis of T1D.

Mutations in the gene encoding the 3'-5' DNA exonuclease TREX1 are associated with systemic lupus erythematosus.

TREX1 acts in concert with the SET complex in granzyme A-mediated apoptosis, and mutations in TREX1 cause Aicardi-Goutières syndrome and familial chilblain lupus. Here, we report monoallelic frameshift or missense mutations and one 3' UTR variant of TREX1 present in 9/417 individuals with systemic lupus erythematosus but absent in 1,712 controls (P = 4.1 x 10(-7)). We demonstrate that two mutant TREX1 alleles alter subcellular targeting. Our findings implicate TREX1 in the pathogenesis of SLE.

Protease inhibitor 15, a candidate gene for abdominal aortic internal elastic lamina ruptures in the rat.

The inbred Brown Norway (BN) rat develops spontaneous ruptures of the internal elastic lamina (RIEL) of the abdominal aorta (AA) and iliac arteries. Prior studies with crosses of the BN/Orl RJ (susceptible) and LOU/M (resistant) showed the presence of a significant QTL on chromosome 5 and the production of congenic rats proved the involvement of this locus. In this study, we further dissected the above-mentioned QTL by creating a new panel of LOU.BN(chr5) congenic and subcongenic lines and reduced the locus to 5.2 Mb. Then we studied 1,002 heterogeneous stock (HS) rats, whose phenotyping revealed a low prevalence and high variability for RIEL. High-resolution mapping in the HS panel detected the major locus on chromosome 5 (log P > 35) and refined it to 1.4 Mb. Subsequently, RNA-seq analysis on AA of BN, congenics, and LOU revealed expression differences for only protease inhibitor 15 (Pi15) gene and a putative long intergenic noncoding RNA (lincRNA) within the linkage region. The high abundance of lincRNA with respect to reduced Pi15 expression, in conjunction with exertion of longitudinal strain, may be related to RIEL, indicating the potential importance of proteases in biological processes related to defective aortic internal elastic lamina structure. Similar mechanisms may be involved in aneurysm initiation in the human AA.

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease.

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.

Evolution of translational control and the emergence of genes and open reading frames in human and non-human primate hearts.

Evolutionary innovations can be driven by changes in the rates of RNA translation and the emergence of new genes and small open reading frames (sORFs). In this study, we characterized the transcriptional and translational landscape of the hearts of four primate and two rodent species through integrative ribosome and transcriptomic profiling, including adult left ventricle tissues and induced pluripotent stem cell-derived cardiomyocyte cell cultures. We show here that the translational efficiencies of subunits of the mitochondrial oxidative phosphorylation chain complexes IV and V evolved rapidly across mammalian evolution. Moreover, we discovered hundreds of species-specific and lineage-specific genomic innovations that emerged during primate evolution in the heart, including 551 genes, 504 sORFs and 76 evolutionarily conserved genes displaying human-specific cardiac-enriched expression. Overall, our work describes the evolutionary processes and mechanisms that have shaped cardiac transcription and translation in recent primate evolution and sheds light on how these can contribute to cardiac development and disease.

PITX2 deficiency leads to atrial mitochondrial dysfunction.

AIM: Reduced left atrial PITX2 is associated with atrial cardiomyopathy and atrial fibrillation. PITX2 is restricted to left atrial cardiomyocytes in the adult heart. The links between PITX2 deficiency, atrial cardiomyopathy and atrial fibrillation are not fully understood. METHODS AND RESULTS: To identify mechanisms linking PITX2 deficiency to atrial fibrillation, we generated and characterized PITX2-deficient human atrial cardiomyocytes derived from human induced pluripotent stem cells (hiPSC) and their controls. PITX2-deficient hiPSC-derived atrial cardiomyocytes showed shorter and disorganised sarcomeres and increased mononucleation. Electron microscopy found an increased number of smaller mitochondria compared to the control. Mitochondrial protein expression was altered in PITX2-deficient hiPSC-derived atrial cardiomyocytes. Single-nuclear RNA-sequencing found differences in cellular respiration pathways and differentially expressed mitochondrial and ion channel genes in PITX2-deficient hiPSC-derived atrial cardiomyocytes. PITX2 repression in hiPSC-derived atrial cardiomyocytes replicated dysregulation of cellular respiration. Mitochondrial respiration was shifted to increased glycolysis in PITX2-deficient hiPSC-derived atrial cardiomyocytes. PITX2-deficient human hiPSC-derived atrial cardiomyocytes showed higher spontaneous beating rates. Action potential duration was more variable with an overall prolongation of early repolarization, consistent with metabolic defects. Gene expression analyses confirmed changes in mitochondrial genes in left atria from 42 patients with atrial fibrillation compared to 43 patients in sinus rhythm. Dysregulation of left atrial mitochondrial (COX7C) and metabolic (FOXO1) genes was associated with PITX2 expression in human left atria. CONCLUSIONS: In summary, PITX2 deficiency causes mitochondrial dysfunction and a metabolic shift to glycolysis in human atrial cardiomyocytes. PITX2-dependent metabolic changes can contribute to the structural and functional defects found in PITX2-deficient atria.

The genome sequence of the spontaneously hypertensive rat: Analysis and functional significance.

The spontaneously hypertensive rat (SHR) is the most widely studied animal model of hypertension. Scores of SHR quantitative loci (QTLs) have been mapped for hypertension and other phenotypes. We have sequenced the SHR/OlaIpcv genome at 10.7-fold coverage by paired-end sequencing on the Illumina platform. We identified 3.6 million high-quality single nucleotide polymorphisms (SNPs) between the SHR/OlaIpcv and Brown Norway (BN) reference genome, with a high rate of validation (sensitivity 96.3%-98.0% and specificity 99%-100%). We also identified 343,243 short indels between the SHR/OlaIpcv and reference genomes. These SNPs and indels resulted in 161 gain or loss of stop codons and 629 frameshifts compared with the BN reference sequence. We also identified 13,438 larger deletions that result in complete or partial absence of 107 genes in the SHR/OlaIpcv genome compared with the BN reference and 588 copy number variants (CNVs) that overlap with the gene regions of 688 genes. Genomic regions containing genes whose expression had been previously mapped as cis-regulated expression quantitative trait loci (eQTLs) were significantly enriched with SNPs, short indels, and larger deletions, suggesting that some of these variants have functional effects on gene expression. Genes that were affected by major alterations in their coding sequence were highly enriched for genes related to ion transport, transport, and plasma membrane localization, providing insights into the likely molecular and cellular basis of hypertension and other phenotypes specific to the SHR strain. This near complete catalog of genomic differences between two extensively studied rat strains provides the starting point for complete elucidation, at the molecular level, of the physiological and pathophysiological phenotypic differences between individuals from these strains.

HOXB4 can enhance the differentiation of embryonic stem cells by modulating the hematopoietic niche.

Hematopoietic differentiation of embryonic stem cells (ESCs) in vitro has been used as a model to study early hematopoietic development, and it is well documented that hematopoietic differentiation can be enhanced by overexpression of HOXB4. HOXB4 is expressed in hematopoietic progenitor cells (HPCs) where it promotes self-renewal, but it is also expressed in the primitive streak of the gastrulating embryo. This led us to hypothesize that HOXB4 might modulate gene expression in prehematopoietic mesoderm and that this property might contribute to its prohematopoietic effect in differentiating ESCs. To test our hypothesis, we developed a conditionally activated HOXB4 expression system using the mutant estrogen receptor (ER(T2)) and showed that a pulse of HOXB4 prior to HPC emergence in differentiating ESCs led to an increase in hematopoietic differentiation. Expression profiling revealed an increase in the expression of genes associated with paraxial mesoderm that gives rise to the hematopoietic niche. Therefore, we considered that HOXB4 might modulate the formation of the hematopoietic niche as well as the production of hematopoietic cells per se. Cell mixing experiments supported this hypothesis demonstrating that HOXB4 activation can generate a paracrine as well as a cell autonomous effect on hematopoietic differentiation. We provide evidence to demonstrate that this activity is partly mediated by the secreted protein FRZB.

Natural genetic variation quantitatively regulates heart rate and dimension.

The polygenic contribution to heart development and function along the health-disease continuum remains unresolved. To gain insight into the genetic basis of quantitative cardiac phenotypes, we utilize highly inbred Japanese rice fish models, Oryzias latipes, and Oryzias sakaizumii. Employing automated quantification of embryonic heart rates as core metric, we profiled phenotype variability across five inbred strains. We observed maximal phenotypic contrast between individuals of the HO5 and the HdrR strain. HO5 showed elevated heart rates associated with embryonic ventricular hypoplasia and impaired adult cardiac function. This contrast served as the basis for genome-wide mapping. In a segregation population of 1192 HO5 x HdrR F2 embryos, we mapped 59 loci (173 genes) associated with heart rate. Experimental validation of the top 12 candidate genes in loss-of-function models revealed their causal and distinct impact on heart rate, development, ventricle size, and arrhythmia. Our study uncovers new diagnostic and therapeutic targets for developmental and electrophysiological cardiac diseases and provides a novel scalable approach to investigate the intricate genetic architecture of the vertebrate heart.

Decoding fMRI alcohol cue reactivity and its association with drinking behaviour.

BACKGROUND: Cue reactivity, the enhanced sensitivity to conditioned cues, is associated with habitual and compulsive alcohol consumption. However, most previous studies in alcohol use disorder (AUD) compared brain activity between alcohol and neutral conditions, solely as cue-triggered neural reactivity. OBJECTIVE: This study aims to find the neural subprocesses during the processing of visual alcohol cues in AUD individuals, and how these neural patterns are predictive for relapse. METHODS: Using cue reactivity and rating tasks, we separately modelled the patterns decoding the processes of visual object recognition and reward appraisal of alcohol cues with representational similarity analysis, and compared the decoding involvements (ie, distance between neural responses and hypothesised decoding models) between AUD and healthy individuals. We further explored connectivity between the identified neural systems and the whole brain and predicted relapse within 6 months using decoding involvements of the neural patterns. FINDINGS: AUD individuals, compared with healthy individuals, showed higher involvement of motor-related brain regions in decoding visual features, and their reward, habit and executive networks were more engaged in appraising reward values. Connectivity analyses showed the involved neural systems were widely connected with higher cognitive networks during alcohol cue processing in AUD individuals, and decoding involvements of frontal eye fields and dorsolateral prefrontal cortex could contribute to relapse prediction. CONCLUSIONS: These findings provide insight into how AUD individuals differently decode alcohol cues compared with healthy participants, from the componential processes of visual object recognition and reward appraisal. CLINICAL IMPLICATIONS: The identified patterns are suggested as biomarkers and potential therapeutic targets in AUD.

Characterization of the genomic structure and function of regions influencing renin and angiogenesis in the SS rat.

Impaired regulation of renin in Dahl salt-sensitive rats (SS/JRHsdMcwi, SS) contributes to attenuated angiogenesis in this strain. This study examined angiogenic function and genomic structure of regions surrounding the renin gene using subcongenic strains of the SS and BN/NHsdMcwi (BN) rat to identify important genomic variations between SS and BN involved in angiogenesis. Three candidate regions on Chr 13 were studied: two congenic strains containing 0.89 and 2.62 Mb portions of BN Chr 13 that excluded the BN renin allele and a third strain that contained a 2.02 Mb overlapping region that included the BN renin allele. Angiogenesis induced by electrical stimulation of the tibialis anterior muscle was attenuated in the SS compared with the BN. Congenics carrying the SS renin allele had impaired angiogenesis, while strains carrying the BN renin allele had angiogenesis restored. The exception was a congenic including a region of BN genome 0.4 Mb distal to renin that restored both renin regulation and angiogenesis. This suggests that there is a distant regulatory element in the BN capable of restoring normal regulation of the SS renin allele. The importance of ANG II in the restored angiogenic response was demonstrated by blocking with losartan. Sequencing of the 4.05 Mb candidate region in SS and BN revealed a total of 8,850 SNPs and other sequence variants. An analysis of the genes and their variants in the region suggested a number of pathways that may explain the impaired regulation of renin and angiogenesis in the SS rat.

Prolactin-sensitive olfactory sensory neurons regulate male preference in female mice by modulating responses to chemosensory cues.

Chemosensory cues detected in the nose need to be integrated with the hormonal status to trigger appropriate behaviors, but the neural circuits linking the olfactory and the endocrine system are insufficiently understood. Here, we characterize olfactory sensory neurons in the murine nose that respond to the pituitary hormone prolactin. Deletion of prolactin receptor in these cells results in impaired detection of social odors and blunts male preference in females. The prolactin-responsive olfactory sensory neurons exhibit a distinctive projection pattern to the brain that is similar across different individuals and express a limited subset of chemosensory receptors. Prolactin modulates the responses within these neurons to discrete chemosensory cues contained in male urine, providing a mechanism by which the hormonal status can be directly linked with distinct olfactory cues to generate appropriate behavioral responses.

A trans locus causes a ribosomopathy in hypertrophic hearts that affects mRNA translation in a protein length-dependent fashion.

BACKGROUND: Little is known about the impact of trans-acting genetic variation on the rates with which proteins are synthesized by ribosomes. Here, we investigate the influence of such distant genetic loci on the efficiency of mRNA translation and define their contribution to the development of complex disease phenotypes within a panel of rat recombinant inbred lines. RESULTS: We identify several tissue-specific master regulatory hotspots that each control the translation rates of multiple proteins. One of these loci is restricted to hypertrophic hearts, where it drives a translatome-wide and protein length-dependent change in translational efficiency, altering the stoichiometric translation rates of sarcomere proteins. Mechanistic dissection of this locus across multiple congenic lines points to a translation machinery defect, characterized by marked differences in polysome profiles and misregulation of the small nucleolar RNA SNORA48. Strikingly, from yeast to humans, we observe reproducible protein length-dependent shifts in translational efficiency as a conserved hallmark of translation machinery mutants, including those that cause ribosomopathies. Depending on the factor mutated, a pre-existing negative correlation between protein length and translation rates could either be enhanced or reduced, which we propose to result from mRNA-specific imbalances in canonical translation initiation and reinitiation rates. CONCLUSIONS: We show that distant genetic control of mRNA translation is abundant in mammalian tissues, exemplified by a single genomic locus that triggers a translation-driven molecular mechanism. Our work illustrates the complexity through which genetic variation can drive phenotypic variability between individuals and thereby contribute to complex disease.

Characterization of the phosphoinositide 3-kinase-dependent transcriptome in murine embryonic stem cells: identification of novel regulators of pluripotency.

Phosphoinositide 3-kinase (PI3K)-dependent signaling has been implicated in the regulation of embryonic stem (ES) cell fate. To gain further insight into the mechanisms regulated by PI3Ks in murine ES cells, we have performed expression profiling using Affymetrix GeneChips to characterize the transcriptional changes that arise as a result of inhibition of PI3K-dependent signaling. Using filtering of greater than 1.5-fold change in expression and an analysis of variance significance level of p < .05, we have defined a dataset comprising 646 probe sets that detect changes in transcript expression (469 down and 177 up) on inhibition of PI3Ks. Changes in expression of selected genes have been validated by quantitative reverse transcription polymerase chain reaction. Gene ontology analyses reveal significant over-representation of transcriptional regulators within our dataset. In addition, several known regulators of ES cell pluripotency, for example, Nanog, Esrrb, Tbx3, and Tcl-1, are among the downregulated genes. To evaluate the functional involvement of selected genes in regulation of ES cell self-renewal, we have used short interfering RNA-mediated knockdown. These studies identify genes not previously associated with control of ES cell fate that are involved in regulating ES cell pluripotency, including the protein tyrosine phosphatase Shp-1 and the Zscan4 family of zinc finger proteins. Further gain-of-function analyses demonstrate the importance of Zscan4c in regulation of ES cell pluripotency.

Novel STAT3 target genes exert distinct roles in the inhibition of mesoderm and endoderm differentiation in cooperation with Nanog.

Leukemia inhibitory factor (LIF) activates the transcription factor signal transducer and activator of transcription 3 (STAT3), which results in the maintenance of mouse embryonic stem cells in the pluripotent state by inhibiting both mesodermal and endodermal differentiation. How the LIF/STAT3 pathway inhibits commitment to both mesoderm and endoderm lineages is presently unknown. Using a hormone-dependent STAT3 and with microarray analysis, we identified 58 targets of STAT3 including 20 unknown genes. Functional analysis showed that 22 among the 23 STAT3 target genes analyzed contribute to the maintenance of the undifferentiated state, as evidenced by an increase in the frequency of differentiated colonies in a self-renewal assay and a concomitant elevation of early differentiation markers upon knockdown. Fourteen of them, including Dact1, Klf4, Klf5, Rgs16, Smad7, Ccrn4l, Cnnm1, Ocln, Ier3, Pim1, Cyr61, and Sgk, were also regulated by Nanog. Analysis of lineage-specific markers showed that the STAT3 target genes fell into three distinct categories, depending on their capacity to inhibit either mesoderm or endoderm differentiation or both. The identification of genes that harness self-renewal and are downstream targets of both STAT3 and Nanog shed light on the mechanisms underlying functional redundancy between STAT3 and Nanog in mouse embryonic stem cells.

Phosphatidylinositol 4-kinase ß mutations cause nonsyndromic sensorineural deafness and inner ear malformation.

Congenital hearing loss is a common disorder worldwide. Heterogeneous gene variation accounts for approximately 20-25% of such patients. We investigated a five-generation Chinese family with autosomal-dominant nonsyndromic sensorineural hearing loss (SNHL). No wave was detected in the pure-tone audiometry, and the auditory brainstem response was absent in all patients. Computed tomography of the patients, as well as of two sporadic SNHL cases, showed bilateral inner ear anomaly, cochlear maldevelopment, absence of the osseous spiral lamina, and an enlarged vestibular aqueduct. Such findings were absent in nonaffected persons. We used linkage analysis and exome sequencing and uncovered a heterozygous missense mutation in the PI4KB gene (p.Gln121Arg) encoding phosphatidylinositol 4-kinase ß (PI4KB) from the patients in this family. In addition, 3 missense PI4KB (p.Val434Gly, p.Glu667Lys, and p.Met739Arg) mutations were identified in five patients with nonsyndromic SNHL from 57 sporadic cases. No such mutations were present within 600 Chinese controls, the 1000 genome project, gnomAD, or similar databases. Depleting pi4kb mRNA expression in zebrafish caused inner ear abnormalities and audiosensory impairment, mimicking the patient phenotypes. Moreover, overexpression of 4 human missense PI4KB mutant mRNAs in zebrafish embryos resulted in impaired hearing function, suggesting dominant-negative effects. Taken together, our results reveal that PI4KB mutations can cause SNHL and inner ear malformation. PI4KB should be included in neonatal deafness screening.

Three LIF-dependent signatures and gene clusters with atypical expression profiles, identified by transcriptome studies in mouse ES cells and early derivatives.

BACKGROUND: Mouse embryonic stem (ES) cells remain pluripotent in vitro when grown in the presence of the cytokine Leukaemia Inhibitory Factor (LIF). Identification of LIF targets and of genes regulating the transition between pluripotent and early differentiated cells is a critical step for understanding the control of ES cell pluripotency. RESULTS: By gene profiling studies carried out with mRNAs from ES cells and their early derivatives treated or not with LIF, we have identified i) LIF-dependent genes, highly expressed in pluripotent cells, whose expression level decreases sharply upon LIF withdrawal [Pluri genes], ii) LIF induced genes [Lifind genes] whose expression is differentially regulated depending upon cell context and iii) genes specific to the reversible or irreversible committed states. In addition, by hierarchical gene clustering, we have identified, among eight independent gene clusters, two atypical groups of genes, whose expression level was highly modulated in committed cells only. Computer based analyses led to the characterization of different sub-types of Pluri and Lifind genes, and revealed their differential modulation by Oct4 or Nanog master genes. Individual knock down of a selection of Pluri and Lifind genes leads to weak changes in the expression of early differentiation markers, in cell growth conditions in which these master genes are still expressed. CONCLUSION: We have identified different sets of LIF-regulated genes depending upon the cell state (reversible or irreversible commitment), which allowed us to present a novel global view of LIF responses. We are also reporting on the identification of genes whose expression is strictly regulated during the commitment step. Furthermore, our studies identify sub-networks of genes with a restricted expression in pluripotent ES cells, whose down regulation occurs while the master knot (composed of OCT4, SOX2 and NANOG) is still expressed and which might be down-regulated together for driving cells towards differentiation.