Background: Publicly available RNA-seq datasets are often underused although being helpful to improve functional annotation of eukaryotic genomes. This is especially true for filamentous fungi genomes which structure differs from most well annotated yeast genomes. Podospora anserina is a filamentous fungal model, which genome has been sequenced and annotated in 2008. Still, the current annotation lacks information about cis-regulatory elements, including promoters, transcription starting sites and terminators, which are instrumental to integrate epigenomic features into global gene regulation strategies. Results: Here we took advantage of 37 RNA-seq experiments that were obtained in contrasted developmental and physiological conditions, to complete the functional annotation of P. anserina genome. Out of the 10,800 previously annotated genes, 5’UTR and 3’UTR were defined for 7554, among which, 3328 showed differential transcriptional signal starts and/or transcriptional end sites. In addition, alternative splicing events were detected for 2350 genes, mostly due alternative 3’splice sites and 1732 novel transcriptionally active regions (nTARs) in unannotated regions were identified. Conclusions: Our study provides a comprehensive genome-wide functional annotation of P. anserina genome, including chromatin features, cis-acting elements such as UTRs, alternative splicing events and transcription of non-coding regions. These new findings will likely improve our understanding of gene regulation strategies in compact genomes, such as those of filamentous fungi. Characterization of alternative transcripts and nTARs paves the way to the discovery of putative new genes, alternative peptides or regulatory non-coding RNAs.

Multinucleate cells are found in many eukaryotes, but how multiple nuclei coordinate their functions is still poorly understood. In the cytoplasm of the ciliate Paramecium tetraurelia, two micronuclei (MIC) serving sexual reproduction coexist with a somatic macronucleus (MAC) dedicated to gene expression. During sexual processes, the MAC is progressively destroyed while still ensuring transcription, and new MACs develop from copies of the zygotic MIC. Several gene clusters are successively induced and switched off before vegetative growth resumes. Concomitantly, programmed genome rearrangement (PGR) removes transposons and their relics from the new MACs. Development of the new MACs is controlled by the old MAC, since the latter expresses genes involved in PGR, including the PGM gene encoding the essential PiggyMac endonuclease that cleaves the ends of eliminated sequences. Using RNA deep sequencing and transcriptome analysis, we show that impairing PGR upregulates key known PGR genes, together with ∼600 other genes possibly also involved in PGR. Among these genes, 42% are no longer induced when no new MACs are formed, including 180 genes that are co-expressed with PGM under all tested conditions. We propose that bi-directional crosstalk between the two coexisting generations of MACs links gene expression to the progression of MAC development.

Candida glabrata is a cause of life-threatening invasive infections especially in elderly and immunocompromised patients. Part of human digestive and urogenital microbiota, C. glabrata faces varying iron availability, low during infection or high in digestive and urogenital tracts. To maintain its homeostasis, C. glabrata must get enough iron for essential cellular processes and resist toxic iron excess. The response of this pathogen to both depletion and lethal excess of iron at 30°C have been described in the literature using different strains and iron sources. However, adaptation to iron variations at 37°C, the human body temperature and to gentle overload, is poorly known. In this study, we performed transcriptomic experiments at 30°C and 37°C with low and high but sub-lethal ferrous concentrations. We identified iron responsive genes and clarified the potential effect of temperature on iron homeostasis. Our exploration of the datasets was facilitated by the inference of functional networks of co-expressed genes, which can be accessed through a web interface. Relying on stringent selection and independently of existing knowledge, we characterized a list of 214 genes as key elements of C. glabrata iron homeostasis and interesting candidates for medical applications.

The functions of eukaryotic chromosomes and their spatial architecture in the nucleus are reciprocally dependent. Hi-C experiments are routinely used to study chromosome 3D organization by probing chromatin interactions. Standard representation of the data has relied on contact maps that show the frequency of interactions between parts of the genome. In parallel, it has become easier to build 3D models of the entire genome based on the same Hi-C data, and thus benefit from the methodology and visualization tools developed for structural biology. 3D modeling of entire genomes leverages the understanding of their spatial organization. However, this opportunity for original and insightful modeling is underexploited. In this paper, we show how seeing the spatial organization of chromosomes can bring new perspectives to omics data integration. We assembled state-of-the-art tools into a workflow that goes from Hi-C raw data to fully annotated 3D models and we re-analysed public omics datasets available for three fungal species. Besides the well-described properties of the spatial organization of their chromosomes (Rabl conformation, hypercoiling and chromosome territories), our results highlighted (i) in Saccharomyces cerevisiae, the backbones of the cohesin anchor regions, which were aligned all along the chromosomes, (ii) in Schizosaccharomyces pombe, the oscillations of the coiling of chromosome arms throughout the cell cycle and (iii) in Neurospora crassa, the massive relocalization of histone marks in mutants of heterochromatin regulators. 3D modeling of the chromosomes brings new opportunities for visual integration of omics data. This holistic perspective supports intuition and lays the foundation for building new concepts.

Label-free quantitative proteomics has emerged as a powerful strategy to obtain high quality quantitative measures of the proteome with only a very small quantity of total protein extract. Because our research projects were requiring the application of bottom-up shotgun mass spectrometry proteomics in the pathogenic yeasts Candida glabrata and Candida albicans, we performed preliminary experiments to (i) obtain a precise list of all the proteins for which measures of abundance could be obtained and (ii) assess the reproducibility of the results arising respectively from biological and technical replicates. Three time-courses were performed in each Candida species, and an alkaline pH stress was induced for two of them. Cells were collected 10 and 60 min after stress induction and proteins were extracted. Samples were analysed two times by mass spectrometry. Our final dataset thus comprises label-free quantitative proteomics results for 24 samples (two species, three time-courses, two time points and two runs of mass spectrometry). Statistical procedures were applied to identify proteins with differential abundances between stressed and unstressed situations. Considering that C. glabrata and C. albicans are human pathogens, which face important pH fluctuations during a human host infection, this dataset has a potential value to other researchers in the field.

Numerous genes are overexpressed in the radioresistant bacterium Deinococcus radiodurans after exposure to radiation or prolonged desiccation. The DdrO and IrrE proteins play a major role in regulating the expression of approximately predicted twenty of these genes. The transcriptional repressor DdrO blocks the expression of these genes under normal growth conditions. After exposure to genotoxic agents, the IrrE metalloprotease cleaves DdrO and relieves gene repression. Bioinformatic analyzes showed that this mechanism seems to be conserved in several species of Deinococcus , but many questions remain as such the number of genes regulated by DdrO. Here, by RNA-seq and CHiP-seq assays performed at a genome-wide scale coupled with bioinformatic analyses, we show that, the DdrO regulon in D. radiodurans includes many other genes than those previously described. These results thus pave the way to better understand the radioresistance mechanisms encoded by this bacterium. Author Summary: The main response pathway to genotoxic conditions in the radioresistant bacterium Deinococcus radiodurans is regulated by the constitutively expressed metalloprotease IrrE that cleaves the transcriptional repressor DdrO, leading to the expression of the genes repressed by DdrO. One of the major goals to better understand how pathways involved in radioresistance are coordinated into this fascinating bacterium is to highlight genes regulated by DdrO. In this study, we mapped in vivo the DdrO regulon in D. radiodurans by using two genome-scale approaches, ChIP-seq and RNA-seq analyses, coupled with bioinformatic analyses. As homologs of these two proteins are also found in many other bacteria, these results also pave the way to compare the stress-induced responses mediated by this couple of proteins in diverse bacteria.