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.

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.

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.

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.

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.

Complete Sequence of Probiotic Symbioflor 2 Escherichia coli Strain G3/10 and Draft Sequences of Symbioflor 2 E. coli Strains G1/2, G4/9, G5, G6/7, and G8.

The complete genome of probiotic Escherichia coli strain G3/10 is presented here. In addition, the probiotic E. coli strains G1/2, G4/9, G5, G6/7, and G8 are presented in draft form. These six strains together comprise the probiotic product Symbioflor 2 (DSM 17252).

A functional insulator screen identifies NURF and dREAM components to be required for enhancer-blocking.

Chromatin insulators of higher eukaryotes functionally divide the genome into active and inactive domains. Furthermore, insulators regulate enhancer/promoter communication, which is evident from the Drosophila bithorax locus in which a multitude of regulatory elements control segment specific gene activity. Centrosomal protein 190 (CP190) is targeted to insulators by CTCF or other insulator DNA-binding factors. Chromatin analyses revealed that insulators are characterized by open and nucleosome depleted regions. Here, we wanted to identify chromatin modification and remodelling factors required for an enhancer blocking function. We used the well-studied Fab-8 insulator of the bithorax locus to apply a genome-wide RNAi screen for factors that contribute to the enhancer blocking function of CTCF and CP190. Among 78 genes required for optimal Fab-8 mediated enhancer blocking, all four components of the NURF complex as well as several subunits of the dREAM complex were most evident. Mass spectrometric analyses of CTCF or CP190 bound proteins as well as immune precipitation confirmed NURF and dREAM binding. Both co-localise with most CP190 binding sites in the genome and chromatin immune precipitation showed that CP190 recruits NURF and dREAM. Nucleosome occupancy and histone H3 binding analyses revealed that CP190 mediated NURF binding results in nucleosomal depletion at CP190 binding sites. Thus, we conclude that CP190 binding to CTCF or to other DNA binding insulator factors mediates recruitment of NURF and dREAM. Furthermore, the enhancer blocking function of insulators is associated with nucleosomal depletion and requires NURF and dREAM.

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.

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.

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.

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.

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.

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.

The Pseudomonas aeruginosa transcriptome in planktonic cultures and static biofilms using RNA sequencing.

In this study, we evaluated how gene expression differs in mature Pseudomonas aeruginosa biofilms as opposed to planktonic cells by the use of RNA sequencing technology that gives rise to both quantitative and qualitative information on the transcriptome. Although a large proportion of genes were consistently regulated in both the stationary phase and biofilm cultures as opposed to the late exponential growth phase cultures, the global biofilm gene expression pattern was clearly distinct indicating that biofilms are not just surface attached cells in stationary phase. A large amount of the genes found to be biofilm specific were involved in adaptation to microaerophilic growth conditions, repression of type three secretion and production of extracellular matrix components. Additionally, we found many small RNAs to be differentially regulated most of them similarly in stationary phase cultures and biofilms. A qualitative analysis of the RNA-seq data revealed more than 3000 putative transcriptional start sites (TSS). By the use of rapid amplification of cDNA ends (5'-RACE) we confirmed the presence of three different TSS associated with the pqsABCDE operon, two in the promoter of pqsA and one upstream of the second gene, pqsB. Taken together, this study reports the first transcriptome study on P. aeruginosa that employs RNA sequencing technology and provides insights into the quantitative and qualitative transcriptome including the expression of small RNAs in P. aeruginosa biofilms.

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.

Genomewide analyses define different modes of transcriptional regulation by peroxisome proliferator-activated receptor-ß/δ (PPARß/δ).

Peroxisome proliferator-activated receptors (PPARs) are nuclear receptors with essential functions in lipid, glucose and energy homeostasis, cell differentiation, inflammation and metabolic disorders, and represent important drug targets. PPARs heterodimerize with retinoid X receptors (RXRs) and can form transcriptional activator or repressor complexes at specific DNA elements (PPREs). It is believed that the decision between repression and activation is generally governed by a ligand-mediated switch. We have performed genomewide analyses of agonist-treated and PPARß/δ-depleted human myofibroblasts to test this hypothesis and to identify global principles of PPARß/δ-mediated gene regulation. Chromatin immunoprecipitation sequencing (ChIP-Seq) of PPARß/δ, H3K4me3 and RNA polymerase II enrichment sites combined with transcriptional profiling enabled the definition of 112 bona fide PPARß/δ target genes showing either of three distinct types of transcriptional response: (I) ligand-independent repression by PPARß/δ; (II) ligand-induced activation and/or derepression by PPARß/δ; and (III) ligand-independent activation by PPARß/δ. These data identify PPRE-mediated repression as a major mechanism of transcriptional regulation by PPARß/δ, but, unexpectedly, also show that only a subset of repressed genes are activated by a ligand-mediated switch. Our results also suggest that the type of transcriptional response by a given target gene is connected to the structure of its associated PPRE(s) and the biological function of its encoded protein. These observations have important implications for understanding the regulatory PPAR network and PPARß/δ ligand-based drugs.

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.

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.