Nonlinear DNA methylation trajectories in aging male mice.

Although DNA methylation data yields highly accurate age predictors, little is known about the dynamics of this quintessential epigenomic biomarker during lifespan. To narrow the gap, we investigate the methylation trajectories of male mouse colon at five different time points of aging. Our study indicates the existence of sudden hypermethylation events at specific stages of life. Precisely, we identify two epigenomic switches during early-to-midlife (3-9 months) and mid-to-late-life (15-24 months) transitions, separating the rodents' life into three stages. These nonlinear methylation dynamics predominantly affect genes associated with the nervous system and enrich in bivalently marked chromatin regions. Based on groups of nonlinearly modified loci, we construct a clock-like classifier STageR (STage of aging estimatoR) that accurately predicts murine epigenetic stage. We demonstrate the universality of our clock in an independent mouse cohort and with publicly available datasets.

Bashing irreproducibility with shournal.

Arguably, the most important tool for many computational scientists is the Linux shell. Processing steps carried out there are critical for a large number of analyses. While the manual documentation of the work is time-consuming and error-prone, existing tools do not integrate well into the shell or suffer from a large overhead. Here, we present shournal, which integrates tightly into the shell and automatically records all shell commands along with their associated file events. Thus, for all files, it can later be told how they were generated and processed. Additionally, it allows the creation of detailed reports for whole project folders. shournal retrieves its data directly from the Linux kernel and allows the monitoring of whole process trees with low overhead.

Multi-omics analysis identifies RFX7 targets involved in tumor suppression and neuronal processes.

Recurrently mutated in lymphoid neoplasms, the transcription factor RFX7 is emerging as a tumor suppressor. Previous reports suggested that RFX7 may also have a role in neurological and metabolic disorders. We recently reported that RFX7 responds to p53 signaling and cellular stress. Furthermore, we found RFX7 target genes to be dysregulated in numerous cancer types also beyond the hematological system. However, our understanding of RFX7's target gene network and its role in health and disease remains limited. Here, we generated RFX7 knock-out cells and employed a multi-omics approach integrating transcriptome, cistrome, and proteome data to obtain a more comprehensive picture of RFX7 targets. We identify novel target genes linked to RFX7's tumor suppressor function and underscoring its potential role in neurological disorders. Importantly, our data reveal RFX7 as a mechanistic link that enables the activation of these genes in response to p53 signaling.

The landscape of human p53-regulated long non-coding RNAs reveals critical host gene co-regulation.

The role of long non-coding RNAs (lncRNAs) in p53-mediated tumor suppression has become increasingly appreciated in the past decade. Thus, the identification of p53-regulated lncRNAs can be a promising starting point to select and prioritize lncRNAs for functional analyses. By integrating transcriptome and transcription factor-binding data, we identified 379 lncRNAs that are recurrently differentially regulated by p53. Dissecting the mechanisms by which p53 regulates many of them, we identified sets of lncRNAs regulated either directly by p53 or indirectly through the p53-RFX7 and p53-p21-DREAM/RB:E2F pathways. Importantly, we identified multiple p53-responsive lncRNAs that are co-regulated with their protein-coding host genes, revealing an important mechanism by which p53 may regulate lncRNAs. Further analysis of transcriptome data and clinical data from cancer patients showed that recurrently p53-regulated lncRNAs are associated with patient survival. Together, the integrative analysis of the landscape of p53-regulated lncRNAs provides a powerful resource facilitating the identification of lncRNA function and displays the mechanisms of p53-dependent regulation that could be exploited for developing anticancer approaches.

TargetGeneReg 2.0: a comprehensive web-atlas for p53, p63, and cell cycle-dependent gene regulation.

In recent years, our web-atlas at www.TargetGeneReg.org has enabled many researchers to uncover new biological insights and to identify novel regulatory mechanisms that affect p53 and the cell cycle - signaling pathways that are frequently dysregulated in diseases like cancer. Here, we provide a substantial upgrade of the database that comprises an extension to include non-coding genes and the transcription factors ΔNp63 and RFX7. TargetGeneReg 2.0 combines gene expression profiling and transcription factor DNA binding data to determine, for each gene, the response to p53, ΔNp63, and cell cycle signaling. It can be used to dissect common, cell type and treatment-specific effects, identify the most promising candidates, and validate findings. We demonstrate the increased power and more intuitive layout of the resource using realistic examples.

p53-mediated AKT and mTOR inhibition requires RFX7 and DDIT4 and depends on nutrient abundance.

In recent years the tumor suppressor p53 has been increasingly recognized as a potent regulator of the cell metabolism and for its ability to inhibit the critical pro-survival kinases AKT and mTOR. The mechanisms through which p53 controls AKT and mTOR, however, are largely unclear. Here, we demonstrate that p53 activates the metabolic regulator DDIT4 indirectly through the regulatory factor X 7 (RFX7). We provide evidence that DDIT4 is required for p53 to inhibit mTOR complex 2 (mTORC2)-dependent AKT activation. Most strikingly, we also find that the DDIT4 regulator RFX7 is required for p53-mediated inhibition of mTORC1 and AKT. Our results suggest that AMPK activation plays no role and p53-mediated AKT inhibition is not critical for p53-mediated mTORC1 inhibition. Moreover, using recently developed physiological cell culture media we uncover that basal p53 and RFX7 activity can play a critical role in restricting mTORC1 activity under physiological nutrient conditions, and we propose a nutrient-dependent model for p53-RFX7-mediated mTORC1 inhibition. These results establish RFX7 and its downstream target DDIT4 as essential effectors in metabolic control elicited by p53.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress.

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7's role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

Dissecting the DNA binding landscape and gene regulatory network of p63 and p53.

The transcription factor p53 is the best-known tumor suppressor, but its sibling p63 is a master regulator of epidermis development and a key oncogenic driver in squamous cell carcinomas (SCC). Despite multiple gene expression studies becoming available, the limited overlap of reported p63-dependent genes has made it difficult to decipher the p63 gene regulatory network. Particularly, analyses of p63 response elements differed substantially among the studies. To address this intricate data situation, we provide an integrated resource that enables assessing the p63-dependent regulation of any human gene of interest. We use a novel iterative de novo motif search approach in conjunction with extensive ChIP-seq data to achieve a precise global distinction between p53-and p63-binding sites, recognition motifs, and potential co-factors. We integrate these data with enhancer:gene associations to predict p63 target genes and identify those that are commonly de-regulated in SCC representing candidates for prognosis and therapeutic interventions.

Multiple Roots of Fruiting Body Formation in Amoebozoa.

Establishment of multicellularity represents a major transition in eukaryote evolution. A subgroup of Amoebozoa, the dictyosteliids, has evolved a relatively simple aggregative multicellular stage resulting in a fruiting body supported by a stalk. Protosteloid amoeba, which are scattered throughout the amoebozoan tree, differ by producing only one or few single stalked spores. Thus, one obvious difference in the developmental cycle of protosteliids and dictyosteliids seems to be the establishment of multicellularity. To separate spore development from multicellular interactions, we compared the genome and transcriptome of a Protostelium species (Protostelium aurantium var. fungivorum) with those of social and solitary members of the Amoebozoa. During fruiting body formation nearly 4,000 genes, corresponding to specific pathways required for differentiation processes, are upregulated. A comparison with genes involved in the development of dictyosteliids revealed conservation of >500 genes, but most of them are also present in Acanthamoeba castellanii for which fruiting bodies have not been documented. Moreover, expression regulation of those genes differs between P. aurantium and Dictyostelium discoideum. Within Amoebozoa differentiation to fruiting bodies is common, but our current genome analysis suggests that protosteliids and dictyosteliids used different routes to achieve this. Most remarkable is both the large repertoire and diversity between species in genes that mediate environmental sensing and signal processing. This likely reflects an immense adaptability of the single cell stage to varying environmental conditions. We surmise that this signaling repertoire provided sufficient building blocks to accommodate the relatively simple demands for cell-cell communication in the early multicellular forms.

Customized workflow development and data modularization concepts for RNA-Sequencing and metatranscriptome experiments.

RNA-Sequencing (RNA-Seq) has become a widely used approach to study quantitative and qualitative aspects of transcriptome data. The variety of RNA-Seq protocols, experimental study designs and the characteristic properties of the organisms under investigation greatly affect downstream and comparative analyses. In this review, we aim to explain the impact of structured pre-selection, classification and integration of best-performing tools within modularized data analysis workflows and ready-to-use computing infrastructures towards experimental data analyses. We highlight examples for workflows and use cases that are presented for pro-, eukaryotic and mixed dual RNA-Seq (meta-transcriptomics) experiments. In addition, we are summarizing the expertise of the laboratories participating in the project consortium "Structured Analysis and Integration of RNA-Seq experiments" (de.STAIR) and its integration with the Galaxy-workbench of the RNA Bioinformatics Center (RBC).

Differential Effects of Vitamins A and D on the Transcriptional Landscape of Human Monocytes during Infection.

Vitamin A and vitamin D are essential nutrients with a wide range of pleiotropic effects in humans. Beyond their well-documented roles in cellular differentiation, embryogenesis, tissue maintenance and bone/calcium homeostasis, both vitamins have attracted considerable attention due to their association with-immunological traits. Nevertheless, our knowledge of their immunomodulatory potential during infection is restricted to single gene-centric studies, which do not reflect the complexity of immune processes. In the present study, we performed a comprehensive RNA-seq-based approach to define the whole immunomodulatory role of vitamins A and D during infection. Using human monocytes as host cells, we characterized the differential role of both vitamins upon infection with three different pathogens: Aspergillus fumigatus, Candida albicans and Escherichia coli. Both vitamins showed an unexpected ability to counteract the pathogen-induced transcriptional responses. Upon infection, we identified 346 and 176 immune-relevant genes that were regulated by atRA and vitD, respectively. This immunomodulatory activity was dependent on the inflammatory stimulus, allowing us to distinguish regulatory patterns which were specific for each stimulatory setting. Moreover, we explored possible direct and indirect mechanisms of vitamin-mediated regulation of the immune response. Our findings highlight the importance of vitamin-monitoring in critically ill patients. Moreover, our results underpin the potential of atRA and vitD as therapeutic options for anti-inflammatory treatment.

Massive Effect on LncRNAs in Human Monocytes During Fungal and Bacterial Infections and in Response to Vitamins A and D.

Mycoses induced by C.albicans or A.fumigatus can cause important host damage either by deficient or exaggerated immune response. Regulation of chemokine and cytokine signaling plays a crucial role for an adequate inflammation, which can be modulated by vitamins A and D. Non-coding RNAs (ncRNAs) as transcription factors or cis-acting antisense RNAs are known to be involved in gene regulation. However, the processes during fungal infections and treatment with vitamins in terms of therapeutic impact are unknown. We show that in monocytes both vitamins regulate ncRNAs involved in amino acid metabolism and immune system processes using comprehensive RNA-Seq analyses. Compared to protein-coding genes, fungi and bacteria induced an expression change in relatively few ncRNAs, but with massive fold changes of up to 4000. We defined the landscape of long-ncRNAs (lncRNAs) in response to pathogens and observed variation in the isoforms composition for several lncRNA following infection and vitamin treatment. Most of the involved antisense RNAs are regulated and positively correlated with their sense protein-coding genes. We investigated lncRNAs with stimulus specific immunomodulatory activity as potential marker genes: LINC00595, SBF2-AS1 (A.fumigatus) and RP11-588G21.2, RP11-394l13.1 (C.albicans) might be detectable in the early phase of infection and serve as therapeutic targets in the future.

Correction: Gene Expansion Shapes Genome Architecture in the Human Pathogen Lichtheimia corymbifera: An Evolutionary Genomics Analysis in the Ancient Terrestrial Mucorales (Mucoromycotina).

[This corrects the article DOI: 10.1371/journal.pgen.1004496.].

Differential transcriptional responses to Ebola and Marburg virus infection in bat and human cells.

The unprecedented outbreak of Ebola in West Africa resulted in over 28,000 cases and 11,000 deaths, underlining the need for a better understanding of the biology of this highly pathogenic virus to develop specific counter strategies. Two filoviruses, the Ebola and Marburg viruses, result in a severe and often fatal infection in humans. However, bats are natural hosts and survive filovirus infections without obvious symptoms. The molecular basis of this striking difference in the response to filovirus infections is not well understood. We report a systematic overview of differentially expressed genes, activity motifs and pathways in human and bat cells infected with the Ebola and Marburg viruses, and we demonstrate that the replication of filoviruses is more rapid in human cells than in bat cells. We also found that the most strongly regulated genes upon filovirus infection are chemokine ligands and transcription factors. We observed a strong induction of the JAK/STAT pathway, of several genes encoding inhibitors of MAP kinases (DUSP genes) and of PPP1R15A, which is involved in ER stress-induced cell death. We used comparative transcriptomics to provide a data resource that can be used to identify cellular responses that might allow bats to survive filovirus infections.

Defining the transcriptomic landscape of Candida glabrata by RNA-Seq.

Candida glabrata is the second most common pathogenic Candida species and has emerged as a leading cause of nosocomial fungal infections. Its reduced susceptibility to antifungal drugs and its close relationship to Saccharomyces cerevisiae make it an interesting research focus. Although its genome sequence was published in 2004, little is known about its transcriptional dynamics. Here, we provide a detailed RNA-Seq-based analysis of the transcriptomic landscape of C. glabrata in nutrient-rich media, as well as under nitrosative stress and during pH shift. Using RNA-Seq data together with state-of-the-art gene prediction tools, we refined the annotation of the C. glabrata genome and predicted 49 novel protein-coding genes. Of these novel genes, 14 have homologs in S. cerevisiae and six are shared with other Candida species. We experimentally validated four novel protein-coding genes of which two are differentially regulated during pH shift and interaction with human neutrophils, indicating a potential role in host-pathogen interaction. Furthermore, we identified 58 novel non-protein-coding genes, 38 new introns and condition-specific alternative splicing. Finally, our data suggest different patterns of adaptation to pH shift and nitrosative stress in C. glabrata, Candida albicans and S. cerevisiae and thus further underline a distinct evolution of virulence in yeast.

Gene expansion shapes genome architecture in the human pathogen Lichtheimia corymbifera: an evolutionary genomics analysis in the ancient terrestrial mucorales (Mucoromycotina).

Lichtheimia species are the second most important cause of mucormycosis in Europe. To provide broader insights into the molecular basis of the pathogenicity-associated traits of the basal Mucorales, we report the full genome sequence of L. corymbifera and compared it to the genome of Rhizopus oryzae, the most common cause of mucormycosis worldwide. The genome assembly encompasses 33.6 MB and 12,379 protein-coding genes. This study reveals four major differences of the L. corymbifera genome to R. oryzae: (i) the presence of an highly elevated number of gene duplications which are unlike R. oryzae not due to whole genome duplication (WGD), (ii) despite the relatively high incidence of introns, alternative splicing (AS) is not frequently observed for the generation of paralogs and in response to stress, (iii) the content of repetitive elements is strikingly low (<5%), (iv) L. corymbifera is typically haploid. Novel virulence factors were identified which may be involved in the regulation of the adaptation to iron-limitation, e.g. LCor01340.1 encoding a putative siderophore transporter and LCor00410.1 involved in the siderophore metabolism. Genes encoding the transcription factors LCor08192.1 and LCor01236.1, which are similar to GATA type regulators and to calcineurin regulated CRZ1, respectively, indicating an involvement of the calcineurin pathway in the adaption to iron limitation. Genes encoding MADS-box transcription factors are elevated up to 11 copies compared to the 1-4 copies usually found in other fungi. More findings are: (i) lower content of tRNAs, but unique codons in L. corymbifera, (ii) Over 25% of the proteins are apparently specific for L. corymbifera. (iii) L. corymbifera contains only 2/3 of the proteases (known to be essential virulence factors) in comparison to R. oryzae. On the other hand, the number of secreted proteases, however, is roughly twice as high as in R. oryzae.

Genomewide comparison and novel ncRNAs of Aquificales.

BACKGROUND: The Aquificales are a diverse group of thermophilic bacteria that thrive in terrestrial and marine hydrothermal environments. They can be divided into the families Aquificaceae, Desulfurobacteriaceae and Hydrogenothermaceae. Although eleven fully sequenced and assembled genomes are available, only little is known about this taxonomic order in terms of RNA metabolism. RESULTS: In this work, we compare the available genomes, extend their protein annotation, identify regulatory sequences, annotate non-coding RNAs (ncRNAs) of known function, predict novel ncRNA candidates, show idiosyncrasies of the genetic decoding machinery, present two different types of transfer-messenger RNAs and variations of the CRISPR systems. Furthermore, we performed a phylogenetic analysis of the Aquificales based on entire genome sequences, and extended this by a classification among all bacteria using 16S rRNA sequences and a set of orthologous proteins.Combining several in silico features (e.g. conserved and stable secondary structures, GC-content, comparison based on multiple genome alignments) with an in vivo dRNA-seq transcriptome analysis of Aquifex aeolicus, we predict roughly 100 novel ncRNA candidates in this bacterium. CONCLUSIONS: We have here re-analyzed the Aquificales, a group of bacteria thriving in extreme environments, sharing the feature of a small, compact genome with a reduced number of protein and ncRNA genes. We present several classical ncRNAs and riboswitch candidates. By combining in silico analysis with dRNA-seq data of A. aeolicus we predict nearly 100 novel ncRNA candidates.