A Direct in vivo RNAi screen identifies MKK4 as a key regulator of liver regeneration.

The liver harbors a distinct capacity for endogenous regeneration; however, liver regeneration is often impaired in disease and therefore insufficient to compensate for the loss of hepatocytes and organ function. Here we describe a functional genetic approach for the identification of gene targets that can be exploited to increase the regenerative capacity of hepatocytes. Pools of small hairpin RNAs (shRNAs) were directly and stably delivered into mouse livers to screen for genes modulating liver regeneration. Our studies identify the dual-specific kinase MKK4 as a master regulator of liver regeneration. MKK4 silencing robustly increased the regenerative capacity of hepatocytes in mouse models of liver regeneration and acute and chronic liver failure. Mechanistically, induction of MKK7 and a JNK1-dependent activation of the AP1 transcription factor ATF2 and the Ets factor ELK1 are crucial for increased regeneration of hepatocytes with MKK4 silencing.

Cyclin-dependent kinase 6 is a chromatin-bound cofactor for NF-κB-dependent gene expression.

Given the intimate link between inflammation and dysregulated cell proliferation in cancer, we investigated cytokine-triggered gene expression in different cell cycle stages. Transcriptome analysis revealed that G1 release through cyclin-dependent kinase 6 (CDK6) and CDK4 primes and cooperates with the cytokine-driven gene response. CDK6 physically and functionally interacts with the NF-κB subunit p65 in the nucleus and is found at promoters of many transcriptionally active NF-κB target genes. CDK6 recruitment to distinct chromatin regions of inflammatory genes was essential for proper loading of p65 to its cognate binding sites and for the function of p65 coactivators, such as TRIP6. Furthermore, cytokine-inducible nuclear translocation and chromatin association of CDK6 depends on the kinase activity of TAK1 and p38. These results have widespread biological implications, as aberrant CDK6 expression or activation that is frequently observed in human tumors modulates NF-κB to shape the cytokine and chemokine repertoires in chronic inflammation and cancer.

Birth, evolution, and transmission of satellite-free mammalian centromeric domains.

Mammalian centromeres are associated with highly repetitive DNA (satellite DNA), which has so far hindered molecular analysis of this chromatin domain. Centromeres are epigenetically specified, and binding of the CENPA protein is their main determinant. In previous work, we described the first example of a natural satellite-free centromere on Equus caballus Chromosome 11. Here, we investigated the satellite-free centromeres of Equus asinus by using ChIP-seq with anti-CENPA antibodies. We identified an extraordinarily high number of centromeres lacking satellite DNA (16 of 31). All of them lay in LINE- and AT-rich regions. A subset of these centromeres is associated with DNA amplification. The location of CENPA binding domains can vary in different individuals, giving rise to epialleles. The analysis of epiallele transmission in hybrids (three mules and one hinny) showed that centromeric domains are inherited as Mendelian traits, but their position can slide in one generation. Conversely, centromere location is stable during mitotic propagation of cultured cells. Our results demonstrate that the presence of more than half of centromeres void of satellite DNA is compatible with genome stability and species survival. The presence of amplified DNA at some centromeres suggests that these arrays may represent an intermediate stage toward satellite DNA formation during evolution. The fact that CENPA binding domains can move within relatively restricted regions (a few hundred kilobases) suggests that the centromeric function is physically limited by epigenetic boundaries.

A Novel Mechanism of Host-Pathogen Interaction through sRNA in Bacterial Outer Membrane Vesicles.

Bacterial outer membrane vesicle (OMV)-mediated delivery of proteins to host cells is an important mechanism of host-pathogen communication. Emerging evidence suggests that OMVs contain differentially packaged short RNAs (sRNAs) with the potential to target host mRNA function and/or stability. In this study, we used RNA-Seq to characterize differentially packaged sRNAs in Pseudomonas aeruginosa OMVs, and to show transfer of OMV sRNAs to human airway cells. We selected one sRNA for further study based on its stable secondary structure and predicted mRNA targets. Our candidate sRNA (sRNA52320), a fragment of a P. aeruginosa methionine tRNA, was abundant in OMVs and reduced LPS-induced as well as OMV-induced IL-8 secretion by cultured primary human airway epithelial cells. We also showed that sRNA52320 attenuated OMV-induced KC cytokine secretion and neutrophil infiltration in mouse lung. Collectively, these findings are consistent with the hypothesis that sRNA52320 in OMVs is a novel mechanism of host-pathogen interaction whereby P. aeruginosa reduces the host immune response.

SUMOylation of the polycomb group protein L3MBTL2 facilitates repression of its target genes.

Lethal(3) malignant brain tumour like 2 (L3MBTL2) is an integral component of the polycomb repressive complex 1.6 (PRC1.6) and has been implicated in transcriptional repression and chromatin compaction. Here, we show that L3MBTL2 is modified by SUMO2/3 at lysine residues 675 and 700 close to the C-terminus. SUMOylation of L3MBTL2 neither affected its repressive activity in reporter gene assays nor it's binding to histone tails in vitro. In order to analyse whether SUMOylation affects binding of L3MBTL2 to chromatin, we performed ChIP-Seq analysis with chromatin of wild-type HEK293 cells and with chromatin of HEK293 cells stably expressing either FLAG-tagged SUMOylation-competent or SUMOylation-defective L3MBTL2. Wild-type FLAG-L3MBTL2 and the SUMOylation-defective FLAG-L3MBTL2 K675/700R mutant essentially occupied the same sites as endogenous L3MBTL2 suggesting that SUMOylation of L3MBTL2 does not affect chromatin binding. However, a subset of L3MBTL2-target genes, particularly those with low L3MBTL2 occupancy including pro-inflammatory genes, was de-repressed in cells expressing the FLAG-L3MBTL2 K675/700R mutant. Finally, we provide evidence that SUMOylation of L3MBTL2 facilitates repression of these PRC1.6-target genes by balancing the local H2Aub1 levels established by the ubiquitinating enzyme RING2 and the de-ubiquitinating PR-DUB complex.

Characterization of the p53 cistrome--DNA binding cooperativity dissects p53's tumor suppressor functions.

p53 protects us from cancer by transcriptionally regulating tumor suppressive programs designed to either prevent the development or clonal expansion of malignant cells. How p53 selects target genes in the genome in a context- and tissue-specific manner remains largely obscure. There is growing evidence that the ability of p53 to bind DNA in a cooperative manner prominently influences target gene selection with activation of the apoptosis program being completely dependent on DNA binding cooperativity. Here, we used ChIP-seq to comprehensively profile the cistrome of p53 mutants with reduced or increased cooperativity. The analysis highlighted a particular relevance of cooperativity for extending the p53 cistrome to non-canonical binding sequences characterized by deletions, spacer insertions and base mismatches. Furthermore, it revealed a striking functional separation of the cistrome on the basis of cooperativity; with low cooperativity genes being significantly enriched for cell cycle and high cooperativity genes for apoptotic functions. Importantly, expression of high but not low cooperativity genes was correlated with superior survival in breast cancer patients. Interestingly, in contrast to most p53-activated genes, p53-repressed genes did not commonly contain p53 binding elements. Nevertheless, both the degree of gene activation and repression were cooperativity-dependent, suggesting that p53-mediated gene repression is largely indirect and mediated by cooperativity-dependently transactivated gene products such as CDKN1A, E2F7 and non-coding RNAs. Since both activation of apoptosis genes with non-canonical response elements and repression of pro-survival genes are crucial for p53's apoptotic activity, the cistrome analysis comprehensively explains why p53-induced apoptosis, but not cell cycle arrest, strongly depends on the intermolecular cooperation of p53 molecules as a possible safeguard mechanism protecting from accidental cell killing.

Differential roles for MBD2 and MBD3 at methylated CpG islands, active promoters and binding to exon sequences.

The heterogeneous collection of nucleosome remodelling and deacetylation (NuRD) complexes can be grouped into the MBD2- or MBD3-containing complexes MBD2-NuRD and MBD3-NuRD. MBD2 is known to bind to methylated CpG sequences in vitro in contrast to MBD3. Although functional differences have been described, a direct comparison of MBD2 and MBD3 in respect to genome-wide binding and function has been lacking. Here, we show that MBD2-NuRD, in contrast to MBD3-NuRD, converts open chromatin with euchromatic histone modifications into tightly compacted chromatin with repressive histone marks. Genome-wide, a strong enrichment for MBD2 at methylated CpG sequences is found, whereas CpGs bound by MBD3 are devoid of methylation. MBD2-bound genes are generally lower expressed as compared with MBD3-bound genes. When depleting cells for MBD2, the MBD2-bound genes increase their activity, whereas MBD2 plus MBD3-bound genes reduce their activity. Most strikingly, MBD3 is enriched at active promoters, whereas MBD2 is bound at methylated promoters and enriched at exon sequences of active genes.

LINT, a novel dL(3)mbt-containing complex, represses malignant brain tumour signature genes.

Mutations in the l(3)mbt tumour suppressor result in overproliferation of Drosophila larval brains. Recently, the derepression of different gene classes in l(3)mbt mutants was shown to be causal for transformation. However, the molecular mechanisms of dL(3)mbt-mediated gene repression are not understood. Here, we identify LINT, the major dL(3)mbt complex of Drosophila. LINT has three core subunits-dL(3)mbt, dCoREST, and dLint-1-and is expressed in cell lines, embryos, and larval brain. Using genome-wide ChIP-Seq analysis, we show that dLint-1 binds close to the TSS of tumour-relevant target genes. Depletion of the LINT core subunits results in derepression of these genes. By contrast, histone deacetylase, histone methylase, and histone demethylase activities are not required to maintain repression. Our results support a direct role of LINT in the repression of brain tumour-relevant target genes by restricting promoter access.

Recruitment of the ATP-dependent chromatin remodeler dMi-2 to the transcribed region of active heat shock genes.

The ATP-dependent chromatin remodeler dMi-2 can play both positive and negative roles in gene transcription. Recently, we have shown that dMi-2 is recruited to the hsp70 gene in a heat shock-dependent manner, and is required to achieve high transcript levels. Here, we use chromatin immunoprecipitation sequencing (ChIP-Seq) to identify other chromatin regions displaying increased dMi-2 binding upon heat shock and to characterize the distribution of dMi-2 over heat shock genes. We show that dMi-2 is recruited to the body of at least seven heat shock genes. Interestingly, dMi-2 binding extends several hundred base pairs beyond the polyadenylation site into the region where transcriptional termination occurs. We find that dMi-2 does not associate with the entire nucleosome-depleted hsp70 locus 87A. Rather, dMi-2 binding is restricted to transcribed regions. Our results suggest that dMi-2 distribution over active heat shock genes are determined by transcriptional activity.

Genome-wide localization and expression profiling establish Sp2 as a sequence-specific transcription factor regulating vitally important genes.

The transcription factor Sp2 is essential for early mouse development and for proliferation of mouse embryonic fibroblasts in culture. Yet its mechanisms of action and its target genes are largely unknown. In this study, we have combined RNA interference, in vitro DNA binding, chromatin immunoprecipitation sequencing and global gene-expression profiling to investigate the role of Sp2 for cellular functions, to define target sites and to identify genes regulated by Sp2. We show that Sp2 is important for cellular proliferation that it binds to GC-boxes and occupies proximal promoters of genes essential for vital cellular processes including gene expression, replication, metabolism and signalling. Moreover, we identified important key target genes and cellular pathways that are directly regulated by Sp2. Most significantly, Sp2 binds and activates numerous sequence-specific transcription factor and co-activator genes, and represses the whole battery of cholesterol synthesis genes. Our results establish Sp2 as a sequence-specific regulator of vitally important genes.

Transcriptome Dynamics of Pseudomonas aeruginosa during Transition from Overlapping To Non-Overlapping Cell Cycles.

Bacteria either duplicate their chromosome once per cell division or a new round of replication is initiated before the cells divide, thus cell cycles overlap. Here, we show that the opportunistic pathogen Pseudomonas aeruginosa switches from fast growth with overlapping cell cycles to sustained slow growth with only one replication round per cell division when cultivated under standard laboratory conditions. The transition was characterized by fast-paced, sequential changes in transcriptional activity along the ori-ter axis of the chromosome reflecting adaptation to the metabolic needs during both growth phases. Quorum sensing (QS) activity was highest at the onset of the slow growth phase with non-overlapping cell cycles. RNA sequencing of subpopulations of these cultures sorted based on their DNA content, revealed a strong gene dosage effect as well as specific expression patterns for replicating and nonreplicating cells. Expression of flagella and mexE, involved in multidrug efflux was restricted to cells that did not replicate, while those that did showed a high activity of the cell division locus and recombination genes. A possible role of QS in the formation of these subpopulations upon switching to non-overlapping cell cycles could be a subject of further research. IMPORTANCE The coordination of gene expression with the cell cycle has so far been studied only in a few bacteria, the bottleneck being the need for synchronized cultures. Here, we determined replication-associated effects on transcription by comparing Pseudomonas aeruginosa cultures that differ in their growth mode and number of replicating chromosomes. We further show that cell cycle-specific gene regulation can be principally identified by RNA sequencing of subpopulations from cultures that replicate only once per cell division and that are sorted according to their DNA content. Our approach opens the possibility to study asynchronously growing bacteria from a wide phylogenetic range and thereby enhance our understanding of the evolution of cell cycle control on the transcriptional level.

Genome sequence of the thermophilic strain Bacillus coagulans 2-6, an efficient producer of high-optical-purity L-lactic acid.

Bacillus coagulans 2-6 is an efficient producer of lactic acid. The genome of B. coagulans 2-6 has the smallest genome among the members of the genus Bacillus known to date. The frameshift mutation at the start of the d-lactate dehydrogenase sequence might be responsible for the production of high-optical-purity l-lactic acid.

Evolutionary conservation of essential and highly expressed genes in Pseudomonas aeruginosa.

BACKGROUND: The constant increase in development and spread of bacterial resistance to antibiotics poses a serious threat to human health. New sequencing technologies are now on the horizon that will yield massive increases in our capacity for DNA sequencing and will revolutionize the drug discovery process. Since essential genes are promising novel antibiotic targets, the prediction of gene essentiality based on genomic information has become a major focus. RESULTS: In this study we demonstrate that pooled sequencing is applicable for the analysis of sequence variations of strain collections with more than 10 individual isolates. Pooled sequencing of 36 clinical Pseudomonas aeruginosa isolates revealed that essential and highly expressed proteins evolve at lower rates, whereas extracellular proteins evolve at higher rates. We furthermore refined the list of experimentally essential P. aeruginosa genes, and identified 980 genes that show no sequence variation at all. Among the conserved nonessential genes we found several that are involved in regulation, motility and virulence, indicating that they represent factors of evolutionary importance for the lifestyle of a successful environmental bacterium and opportunistic pathogen. CONCLUSION: The detailed analysis of a comprehensive set of P. aeruginosa genomes in this study clearly disclosed detailed information of the genomic makeup and revealed a large set of highly conserved genes that play an important role for the lifestyle of this microorganism. Sequencing strain collections enables for a detailed and extensive identification of sequence variations as potential bacterial adaptation processes, e.g., during the development of antibiotic resistance in the clinical setting and thus may be the basis to uncover putative targets for novel treatment strategies.

Complete genome sequence of the myxobacterium Sorangium cellulosum.

The genus Sorangium synthesizes approximately half of the secondary metabolites isolated from myxobacteria, including the anti-cancer metabolite epothilone. We report the complete genome sequence of the model Sorangium strain S. cellulosum So ce56, which produces several natural products and has morphological and physiological properties typical of the genus. The circular genome, comprising 13,033,779 base pairs, is the largest bacterial genome sequenced to date. No global synteny with the genome of Myxococcus xanthus is apparent, revealing an unanticipated level of divergence between these myxobacteria. A large percentage of the genome is devoted to regulation, particularly post-translational phosphorylation, which probably supports the strain's complex, social lifestyle. This regulatory network includes the highest number of eukaryotic protein kinase-like kinases discovered in any organism. Seventeen secondary metabolite loci are encoded in the genome, as well as many enzymes with potential utility in industry.

Haplotypes of the porcine peroxisome proliferator-activated receptor delta gene are associated with backfat thickness.

BACKGROUND: Peroxisome proliferator-activated receptor delta belongs to the nuclear receptor superfamily of ligand-inducible transcription factors. It is a key regulator of lipid metabolism. The peroxisome proliferator-activated receptor delta gene (PPARD) has been assigned to a region on porcine chromosome 7, which harbours a quantitative trait locus for backfat. Thus, PPARD is considered a functional and positional candidate gene for backfat thickness. The purpose of this study was to test this candidate gene hypothesis in a cross of breeds that were highly divergent in lipid deposition characteristics. RESULTS: Screening for genetic variation in porcine PPARD revealed only silent mutations. Nevertheless, significant associations between PPARD haplotypes and backfat thickness were observed in the F2 generation of the Mangalitsa x Piétrain cross as well as a commercial German Landrace population. Haplotype 5 is associated with increased backfat in F2 Mangalitsa x Piétrain pigs, whereas haplotype 4 is associated with lower backfat thickness in the German Landrace population. Haplotype 4 and 5 carry the same alleles at all but one SNP. Interestingly, the opposite effects of PPARD haplotypes 4 and 5 on backfat thickness are reflected by opposite effects of these two haplotypes on PPAR-delta mRNA levels. Haplotype 4 significantly increases PPAR-delta mRNA levels, whereas haplotype 5 decreases mRNA levels of PPAR-delta. CONCLUSION: This study provides evidence for an association between PPARD and backfat thickness. The association is substantiated by mRNA quantification. Further studies are required to clarify, whether the observed associations are caused by PPARD or are the result of linkage disequilibrium with a causal variant in a neighbouring gene.

Integrated Transcriptional Regulatory Network of Quorum Sensing, Replication Control, and SOS Response in Dinoroseobacter shibae.

Quorum sensing (QS) coordinates population wide gene expression of bacterial species. Highly adaptive traits like gene transfer agents (GTA), morphological heterogeneity, type 4 secretion systems (T4SS), and flagella are QS controlled in Dinoroseobacter shibae, a Roseobacter model organism. Its QS regulatory network is integrated with the CtrA phosphorelay that controls cell division in alphaproteobacteria. To elucidate the network topology, we analyzed the transcriptional response of the QS-negative D. shibae strain ΔluxI1 toward externally added autoinducer (AI) over a time period of 3 h. The signaling cascade is initiated by the CtrA phosphorelay, followed by the QS genes and other target genes, including the second messenger c-di-GMP, competence, flagella and pili. Identification of transcription factor binding sites in promoters of QS induced genes revealed the integration of QS, CtrA phosphorelay and the SOS stress response mediated by LexA. The concentration of regulatory genes located close to the origin or terminus of replication suggests that gene regulation and replication are tightly coupled. Indeed, addition of AI first stimulates and then represses replication. The restart of replication comes along with increased c-di-GMP levels. We propose a model in which QS induces replication followed by differentiation into GTA producing and non-producing cells. CtrA-activity is controlled by the c-di-GMP level, allowing some of the daughter cells to replicate again. The size of the GTA producing subpopulation is tightly controlled by QS via the AI Synthase LuxI2. Finally, induction of the SOS response allows for integration of GTA DNA into the host chromosome.

Spiroketal polyketide formation in Sorangium: identification and analysis of the biosynthetic gene cluster for the highly cytotoxic spirangienes.

Natural products constitute important lead structures in drug discovery. In bacteria, they are often synthesized by large, modular multienzyme complexes. Detailed analysis of the biosynthetic machinery should enable its directed engineering and production of desirable analogs. The myxobacterium Sorangium cellulosum So ce90 produces the cytotoxic spiroketal polyketide spirangien, for which we describe the identification and functional analysis of the biosynthetic pathway. The gene cluster spans 88 kb and encodes 7 type I polyketide synthases and additional enzymes such as a stand-alone thioesterase and 2 methyltransferases. Inactivation of two cytochrome P(450) monooxygenase genes resulted in the production of acyclic spirangien derivatives, providing direct evidence for the involvement of these enzymes in spiroketal formation. The presence of large DNA repeats is consistent with multiple rounds of gene duplication during the evolution of the biosynthetic gene locus.

Molecular and biochemical studies of chondramide formation-highly cytotoxic natural products from Chondromyces crocatus Cm c5.

The jaspamide/chondramide family of depsipeptides are mixed PKS/NRPS natural products isolated from marine sponges and a terrestrial myxobacterium that potently affect the function of the actin cytoskeleton. As a first step to improve production in heterologous host cells and permit genetic approaches to novel analogs, we have cloned and characterized the chondramide biosynthetic genes from the myxobacterium Chondromyces crocatus Cm c5. In addition to the expected PKS and NRPS genes, the cluster encodes a rare tyrosine aminomutase for beta-tyrosine formation and a previously unknown tryptophan-2-halogenase. Conditions for gene transfer into C. crocatus Cm c5 were developed, and inactivation of several genes corroborated their proposed function and served to define the boundaries of the cluster. Biochemical characterization of the final NRPS adenylation domain confirmed the direct activation of beta-tyrosine, and fluorinated chondramides were produced through precursor-directed biosynthesis.

Divergent co-transcriptomes of different host cells infected with Toxoplasma gondii reveal cell type-specific host-parasite interactions.

The apicomplexan parasite Toxoplasma gondii infects various cell types in avian and mammalian hosts including humans. Infection of immunocompetent hosts is mostly asymptomatic or benign, but leads to development of largely dormant bradyzoites that persist predominantly within neurons and muscle cells. Here we have analyzed the impact of the host cell type on the co-transcriptomes of host and parasite using high-throughput RNA sequencing. Murine cortical neurons and astrocytes, skeletal muscle cells (SkMCs) and fibroblasts differed by more than 16,200 differentially expressed genes (DEGs) before and after infection with T. gondii. However, only a few hundred of them were regulated by infection and these largely diverged in neurons, SkMCs, astrocytes and fibroblasts indicating host cell type-specific transcriptional responses after infection. The heterogeneous transcriptomes of host cells before and during infection coincided with ~5,400 DEGs in T. gondii residing in different cell types. Finally, we identified gene clusters in both T. gondii and its host, which correlated with the predominant parasite persistence in neurons or SkMCs as compared to astrocytes or fibroblasts. Thus, heterogeneous expression profiles of different host cell types and the parasites' ability to adapting to them may govern the parasite-host cell interaction during toxoplasmosis.

Genomic peculiarity of coding sequences and metabolic potential of probiotic Escherichia coli strain Nissle 1917 inferred from raw genome data.

Probiotic Escherichia coli strain Nissle 1917 (O6:K5:H1) is a commensal E. coli isolate that has a long tradition in medicine for the treatment of various intestinal disorders in humans. To elucidate the molecular basis of its probiotic nature, we started sequencing the genome of this organism with a whole-genome shotgun approach. A 7.8-fold coverage of the genomic sequence has been generated and is now in the finishing stage. To exploit the genome data as early as possible and to generate hypotheses for functional studies, the unfinished sequencing data were analyzed in this work using a new method [Sun, J., Zeng, A.P., 2004. IdentiCS--identification of coding sequence and in silico reconstruction of the metabolic network directly from unannotated low-coverage bacterial genome sequence. BMC Bioinformatics 5, 112] which is particularly suitable for the prediction of coding sequences (CDSs) from unannotated genome sequence. The CDSs predicted for E. coli Nissle 1917 were compared with those of all five other sequenced E. coli strains (E. coli K-12 MG1655, E. coli K-12 W3110, E. coli CFT073, EHEC O157:H7 EDL933 and EHEC O157:H7 Sakai) published to date. Five thousand one hundred and ninety-two CDSs were predicted for E. coli Nissle 1917, of which 1065 were assigned with enzyme EC numbers. The comparison of all predicted CDSs of E. coli Nissle 1917 to the other E. coli strains revealed 108 CDSs specific for this isolate. They are organized as four big genome islands and many other smaller gene clusters. Based on CDSs with EC numbers for enzymes, the potential metabolic network of Nissle 1917 was reconstructed and compared to those of the other five E. coli strains. Overall, the comparative genomic analysis sheds light on the genomic peculiarity of the probiotic E. coli strain Nissle 1917 and is helpful for designing further functional studies long before the sequencing project is completely finished.