Search

Adenosine is known to modulate dopamine responses in several brain areas. Here, we show that tonic activation of adenosine receptors is able to impede desensitization of D[sub 1] dopamine receptors. As measured by cAMP accumulation in transfected COS-7 cells, long-term exposure to dopamine agonists promoted desensitization of D[sub 1B] receptor but not that of D[sub 1A] receptor. The inability of D[sub 1A] receptor to desensitize was a result of the adenosine present in culture medium acting through activation of adenosine A1 receptors. Cell incubation with either adenosine deaminase, CGS-15943, a generic adenosine receptor antagonist, or the A[sub 1] antagonist DPCPX restored the long-term desensitization time-course of D[sub 1A] receptors. In Ltk cells stably expressing A1 adenosine receptors and D[sub 1A] dopamine receptors, pre-treatment of cells with R(-)-PIA, a full A1 receptor agonist, did not significantly inhibit the acute increase in cAMP levels induced by D1 receptor agonists, but blocked desensitization of D[sub 1A] receptors. However, simultaneous activation of A[sub 1] and D[sub 1A] receptors promoted a delayed D[sub 1A] receptor desensitization. This suggests that functional interaction between A[sub 1] and D[sub 1A] receptors may depend on the activation kinetics of components regulating D1 receptor responses, acting differentially on D[sub 1A] and D[sub 1B] receptors. [ABSTRACT FROM AUTHOR]

Abstract Background Drug susceptible clinical isolates of Candida albicans frequently become highly tolerant to drugs during chemotherapy, with dreadful consequences to patient health. We used RNA sequencing (RNA-seq) to analyze the transcriptomes of a CDR (Candida Drug Resistance) strain and its isogenic drug sensitive counterpart. Results RNA-seq unveiled differential expression of 228 genes including a) genes previously identified as involved in CDR, b) genes not previously associated to the CDR phenotype, and c) novel transcripts whose function as a gene is uncharacterized. In particular, we show for the first time that CDR acquisition is correlated with an overexpression of the transcription factor encoding gene CZF1. CZF1 null mutants were susceptible to many drugs, independently of known multidrug resistance mechanisms. We show that CZF1 acts as a repressor of β-glucan synthesis, thus negatively regulating cell wall integrity. Finally, our RNA-seq data allowed us to identify a new transcribed region, upstream of the TAC1 gene, which encodes the major CDR transcriptional regulator. Conclusion Our results open new perspectives of the role of Czf1 and of our understanding of the transcriptional and post-transcriptional mechanisms that lead to the acquisition of drug resistance in C. albicans, with potential for future improvements of therapeutic strategies.

Abstract Background The identification and characterization of the transcriptional regulatory networks governing the physiology and adaptation of microbial cells is a key step in understanding their behaviour. One such wide-domain regulatory circuit, essential to all cells, is carbon catabolite repression (CCR): it allows the cell to prefer some carbon sources, whose assimilation is of high nutritional value, over less profitable ones. In lower multicellular fungi, the C2H2 zinc finger CreA/CRE1 protein has been shown to act as the transcriptional repressor in this process. However, the complete list of its gene targets is not known. Results Here, we deciphered the CRE1 regulatory range in the model cellulose and hemicellulose-degrading fungus Trichoderma reesei (anamorph of Hypocrea jecorina) by profiling transcription in a wild-type and a delta-cre1 mutant strain on glucose at constant growth rates known to repress and de-repress CCR-affected genes. Analysis of genome-wide microarrays reveals 2.8% of transcripts whose expression was regulated in at least one of the four experimental conditions: 47.3% of which were repressed by CRE1, whereas 29.0% were actually induced by CRE1, and 17.2% only affected by the growth rate but CRE1 independent. Among CRE1 repressed transcripts, genes encoding unknown proteins and transport proteins were overrepresented. In addition, we found CRE1-repression of nitrogenous substances uptake, components of chromatin remodeling and the transcriptional mediator complex, as well as developmental processes. Conclusions Our study provides the first global insight into the molecular physiological response of a multicellular fungus to carbon catabolite regulation and identifies several not yet known targets in a growth-controlled environment.

Abstract Background Branchiomotor neurons comprise an important class of cranial motor neurons that innervate the branchial-arch-derived muscles of the face, jaw and neck. They arise in the ventralmost progenitor domain of the rhombencephalon characterized by expression of the homeodomain transcription factors Nkx2.2 and Phox2b. Phox2b in particular plays a key role in the specification of branchiomotor neurons. In its absence, generic neuronal differentiation is defective in the progenitor domain and no branchiomotor neurons are produced. Conversely, ectopic expression of Phox2b in spinal regions of the neural tube promotes cell cycle exit and neuronal differentiation and, at the same time, induces genes and an axonal phenotype characteristic for branchiomotor neurons. How Phox2b exerts its pleiotropic functions, both as a proneural gene and a neuronal subtype determinant, has remained unknown. Results To gain further insights into the genetic program downstream of Phox2b, we searched for novel Phox2b-regulated genes by cDNA microarray analysis of facial branchiomotor neuron precursors from heterozygous and homozygous Phox2b mutant embryos. We selected for functional studies the genes encoding the axonal growth promoter Gap43, the Wnt antagonist Sfrp1 and the transcriptional regulator Sox13, which were not previously suspected to play roles downstream of Phox2b and whose expression was affected by Phox2b misexpression in the spinal cord. While Gap43 did not produce an obvious phenotype when overexpressed in the neural tube, Sfrp1 induced the interneuron marker Lhx1,5 and Sox13 inhibited neuronal differentiation. We then tested whether Sfrp1 and Sox13, which are down-regulated by Phox2b in the facial neuron precursors, would antagonize some aspects of Phox2b activity. Co-expression of Sfrp1 prevented Phox2b from repressing Lhx1,5 and alleviated the commissural axonal phenotype. When expressed together with Sox13, Phox2b was still able to promote cell cycle exit and neuronal differentiation, but the cells failed to relocate to the mantle layer and to extinguish the neural stem cell marker Sox2. Conclusion Our results suggest novel roles for Sfrp1 and Sox13 in neuronal subtype specification and generic neuronal differentiation, respectively, and indicate that down-regulation of Sfrp1 and Sox13 are essential aspects of the genetic program controlled by Phox2b in cranial motoneurons.

Abstract Background Raw data normalization is a critical step in microarray data analysis because it directly affects data interpretation. Most of the normalization methods currently used are included in the R/BioConductor packages but it is often difficult to identify the most appropriate method. Furthermore, the use of R commands for functions and graphics can introduce mistakes that are difficult to trace. We present here a script written in R that provides a flexible means of access to and monitoring of data normalization for two-color microarrays. This script combines the power of BioConductor and R analysis functions and reduces the amount of R programming required. Results Goulphar was developed in and runs using the R language and environment. It combines and extends functions found in BioConductor packages (limma and marray) to correct for dye biases and spatial artifacts. Goulphar provides a wide range of optional and customizable filters for excluding incorrect signals during the pre-processing step. It displays informative output plots, enabling the user to monitor the normalization process, and helps adapt the normalization method appropriately to the data. All these analyses and graphical outputs are presented in a single PDF report. Conclusion Goulphar provides simple, rapid access to the power of the R/BioConductor statistical analysis packages, with precise control and visualization of the results obtained. Complete documentation, examples and online forms for setting script parameters are available from http://transcriptome.ens.fr/goulphar/.

N-terminal ends of polypeptides are critical for the selective co-translational recruitment of N-terminal modification enzymes. However, it is unknown whether specific N-terminal signatures differentially regulate protein fate according to their cellular functions. In this work, we developed an in-silico approach to detect functional preferences in cellular N-terminomes, and identified in S. cerevisiae more than 200 Gene Ontology terms with specific N-terminal signatures. In particular, we discovered that Mitochondrial Targeting Sequences (MTS) show a strong and specific over-representation at position 2 of hydrophobic residues known to define potential substrates of the N-terminal acetyltransferase NatC. We validated mitochondrial precursors as co-translational targets of NatC by selective purification of translating ribosomes, and found that their N-terminal signature is conserved in Saccharomycotina yeasts. Finally, systematic mutagenesis of the position 2 in a prototypal yeast mitochondrial protein confirmed its critical role in mitochondrial protein import. Our work highlights the hydrophobicity of MTS N-terminal residues and their targeting by NatC as important features for the definition of the mitochondrial proteome, providing a molecular explanation for mitochondrial defects observed in yeast or human NatC-depleted cells. Functional mapping of N-terminal residues thus has the potential to support the discovery of novel mechanisms of protein regulation or targeting., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Nashed et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

This study presents a comprehensive investigation of the aggregation behavior of a malonamide extractant molecule ( N , N '-dimethyl, N , N '-dioctylhexylethoxymalonamide (DMDOHEMA)) in three different solvents, including two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids (1-ethyl-1-butylpiperidinium bis(trifluoromethylsulfonyl)imide ([EBPip + ][NTf 2 - ]) and 1-ethyl-1-octylpiperidinium bis(trifluoromethylsulfonyl)imide ([EOPip + ][NTf 2 - ])) and n -dodecane. By combining polarizable molecular dynamics simulations and small-angle X-ray scattering experiments, we extensively investigated the arrangement of supramolecular assemblies of the extractant molecules. Our results showed that the insertion of the alkyl chains of the extractant molecules into the apolar domain of [EOPip + ][NTf 2 - ] has a significant impact on the aggregation behavior of the extractant molecules, leading to the formation of smaller aggregates having a higher dispersion compared to other solvents. These findings provide new insights into the physicochemical properties of this type of system and are crucial in designing more effective solvents for rare earth metal extraction.

The nanoscale organization of some classes of ionic liquids is responsible for their singular properties. In this paper, we use polarizable molecular dynamics simulations and small-angle X-ray scattering to probe the structure of two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids ([EBPip + ][NTf 2 - ] and [EOPip + ][NTf 2 - ]) that differ in the alkyl chain length of their cation. The X-ray scattering intensities calculated numerically, from the radial distribution functions, are in excellent agreement with the experimental data. The analysis of the different contributions of the X-ray scattering data allowed us to highlight the correlations responsible for the low q peak observed for the long-chain alkyl cations. New angular analyses showed that anions were more likely to align with alkyl chains as their size increased, inducing angular correlation between anions at larger distances. They also showed that the long alkyl chains of the cations aligned more with each other than the short ones. These more aligned alkyl chains induce a smaller volume of the apolar microdomains compared to the well-studied imidazolium-based ionic liquids, leading to the smaller correlation distance for piperidinium-based ionic liquids.

Local translation is a conserved mechanism conferring cells the ability to quickly respond to local stimuli. In the brain, it has been recently reported in astrocytes, whose fine processes contact blood vessels and synapses. Yet the specificity and regulation of astrocyte local translation remain unknown. We study hippocampal perisynaptic astrocytic processes (PAPs) and show that they contain the machinery for translation. Using a refined immunoprecipitation technique, we characterize the entire pool of ribosome-bound mRNAs in PAPs and compare it with the one expressed in the whole astrocyte. We find that a specific pool of mRNAs is highly polarized at the synaptic interface. These transcripts encode an unexpected molecular repertoire, composed of proteins involved in iron homeostasis, translation, cell cycle, and cytoskeleton. Remarkably, we observe alterations in global RNA distribution and ribosome-bound status of some PAP-enriched transcripts after fear conditioning, indicating the role of astrocytic local translation in memory and learning., Competing Interests: Declaration of Interests The authors declare no conflict of interest., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Adsorption of prebiotic building blocks is proposed to have played a role in the emergence of life on Earth. The experimental and theoretical study of this phenomenon should be guided by our knowledge of the geochemistry of the habitable early Earth environments, which could have spanned a large range of settings. Adsorption being an interfacial phenomenon, experiments can be built around the minerals that probably exhibited the largest specific surface areas and were the most abundant, i.e., phyllosilicates. Our current work aims at understanding how nucleotides, the building blocks of RNA and DNA, might have interacted with phyllosilicates under various physico-chemical conditions. We carried out and refined batch adsorption studies to explore parameters such as temperature, pH, salinity, etc. We built a comprehensive, generalized model of the adsorption mechanisms of nucleotides onto phyllosilicate particles, mainly governed by phosphate reactivity. More recently, we used surface chemistry and geochemistry techniques, such as vibrational spectroscopy, low pressure gas adsorption, X-ray microscopy, and theoretical simulations, in order to acquire direct data on the adsorption configurations and localization of nucleotides on mineral surfaces. Although some of these techniques proved to be challenging, questioning our ability to easily detect biosignatures, they confirmed and complemented our pre-established model.

In this work, we used comparative transcriptomics to identify regulatory outliers (ROs) in the human pathogen Candida glabrata . ROs are genes that have very different expression patterns compared to their orthologs in other species. From comparative transcriptome analyses of the response of eight yeast species to toxic doses of selenite, a pleiotropic stress inducer, we identified 38 ROs in C. glabrata . Using transcriptome analyses of C. glabrata response to five different stresses, we pointed out five ROs which were more particularly responsive to iron starvation, a process which is very important for C. glabrata virulence. Global chromatin Immunoprecipitation and gene profiling analyses showed that four of these genes are actually new targets of the iron starvation responsive Aft2 transcription factor in C. glabrata . Two of them ( HBS1 and DOM34b ) are required for C. glabrata optimal growth in iron limited conditions. In S. cerevisiae , the orthologs of these two genes are involved in ribosome rescue by the NO GO decay (NGD) pathway. Hence, our results suggest a specific contribution of NGD co-factors to the C. glabrata adaptation to iron starvation.

Brain mural cells form a heterogeneous family which significantly contributes to the maintenance of the blood-brain barrier and regulation of the cerebral blood flow. Current procedures to isolate them cannot specifically separate their distinct subtypes, in particular vascular smooth muscle cells (VSMCs) and mid-capillary pericytes (mcPCs), which differ among others by their expression of smooth muscle actin (SMA). We herein describe an innovative method allowing SMA + VSMCs and SMA - mcPCs to be freshly isolated from the rat cerebral cortex. Using differential RNA-Seq analysis, we then reveal the specific gene expression profile of each subtype. Our results refine the current description of the role of VSMCs in parenchymal cortical arterioles at the molecular level and provide a unique platform to identify the molecular mechanisms underlying the specific functions of mcPCs in the brain vasculature.

In Drosophila, ubiquitous expression of a short Cyclin G isoform generates extreme developmental noise estimated by fluctuating asymmetry (FA), providing a model to tackle developmental stability. This transcriptional cyclin interacts with chromatin regulators of the Enhancer of Trithorax and Polycomb (ETP) and Polycomb families. This led us to investigate the importance of these interactions in developmental stability. Deregulation of Cyclin G highlights an organ intrinsic control of developmental noise, linked to the ETP-interacting domain, and enhanced by mutations in genes encoding members of the Polycomb Repressive complexes PRC1 and PR-DUB. Deep-sequencing of wing imaginal discs deregulating CycG reveals that high developmental noise correlates with up-regulation of genes involved in translation and down-regulation of genes involved in energy production. Most Cyclin G direct transcriptional targets are also direct targets of PRC1 and RNAPolII in the developing wing. Altogether, our results suggest that Cyclin G, PRC1 and PR-DUB cooperate for developmental stability., Competing Interests: The authors have declared that no competing interests exist.

Background: Study of meta-transcriptomic datasets involving non-model organisms represents bioinformatic challenges. The production of chimeric sequences and our inability to distinguish the taxonomic origins of the sequences produced are inherent and recurrent difficulties in de novo assembly analyses. As the study of holobiont meta-transcriptomes is affected by challenges invoked above, we propose an innovative bioinformatic approach to tackle such difficulties and tested it on marine models as a proof of concept., Results: We considered three holobiont models, of which two transcriptomes were previously published and a yet unpublished transcriptome, to analyze and sort their raw reads using Short Read Connector, a k-mer based similarity method. Before assembly, we thus defined four distinct categories for each holobiont meta-transcriptome: host reads, symbiont reads, shared reads, and unassigned reads. Afterwards, we observed that independent de novo assemblies for each category led to a diminution of the number of chimeras compared to classical assembly methods. Moreover, the separation of each partner's transcriptome offered the independent and comparative exploration of their functional diversity in the holobiont. Finally, our strategy allowed to propose new functional annotations for two well-studied holobionts (a Cnidaria-Dinophyta, a Porifera-Bacteria) and a first meta-transcriptome from a planktonic Radiolaria-Dinophyta system forming widespread symbiotic association for which our knowledge is considerably limited., Conclusions: In contrast to classical assembly approaches, our bioinformatic strategy generates less de novo assembled chimera and allows biologists to study separately host and symbiont data from a holobiont mixture. The pre-assembly separation of reads using an efficient tool as Short Read Connector is an effective way to tackle meta-transcriptomic challenges and offers bright perpectives to study holobiont systems composed of either well-studied or poorly characterized symbiotic lineages and ultimately expand our knowledge about these associations.

Dinoflagellates are one of the most abundant and functionally diverse groups of eukaryotes. Despite an overall scarcity of genomic information for dinoflagellates, constantly emerging high-throughput sequencing resources can be used to characterize and compare these organisms. We assembled de novo and processed 46 dinoflagellate transcriptomes and used a sequence similarity network (SSN) to compare the underlying genomic basis of functional features within the group. This approach constitutes the most comprehensive picture to date of the genomic potential of dinoflagellates. A core-predicted proteome composed of 252 connected components (CCs) of putative conserved protein domains (pCDs) was identified. Of these, 206 were novel and 16 lacked any functional annotation in public databases. Integration of functional information in our network analyses allowed investigation of pCDs specifically associated with functional traits. With respect to toxicity, sequences homologous to those of proteins found in species with toxicity potential (e.g., sxtA4 and sxtG) were not specific to known toxin-producing species. Although not fully specific to symbiosis, the most represented functions associated with proteins involved in the symbiotic trait were related to membrane processes and ion transport. Overall, our SSN approach led to identification of 45,207 and 90,794 specific and constitutive pCDs of, respectively, the toxic and symbiotic species represented in our analyses. Of these, 56% and 57%, respectively (i.e., 25,393 and 52,193 pCDs), completely lacked annotation in public databases. This stresses the extent of our lack of knowledge, while emphasizing the potential of SSNs to identify candidate pCDs for further functional genomic characterization., (© 2018 John Wiley & Sons Ltd.)

Astrocytes are the most abundant glial cells of the central nervous system and have recently been recognized as crucial in the regulation of brain immunity. In most neuropathological conditions, astrocytes are prone to a radical phenotypical change called reactivity, which plays a key role in astrocyte contribution to neuroinflammation. However, how astrocytes regulate brain immunity in healthy conditions is an understudied question. One of the astroglial molecule involved in these regulations might be Connexin 43 (Cx43), a gap junction protein highly enriched in astrocyte perivascular endfeet-terminated processes forming the glia limitans. Indeed, Cx43 deletion in astrocytes (Cx43KO) promotes a continuous immune recruitment and an autoimmune response against an astrocyte protein, without inducing any brain lesion. To investigate the molecular basis of this unique immune response, we characterized the polysomal transcriptome of hippocampal astrocytes deleted for Cx43. Our results demonstrate that, in the absence of Cx43, astrocytes adopt an atypical reactive status with no change in most canonical astrogliosis markers, but with an upregulation of molecules promoting immune recruitment, complement activation as well as anti-inflammatory processes. Intriguingly, while several of these upregulated transcriptional events suggested an activation of the γ-interferon pathway, no increase in this cytokine or activation of related signaling pathways were found in Cx43KO. Finally, deletion of astroglial Cx43 was associated with the upregulation of several angiogenic factors, consistent with an increase in microvascular density in Cx43KO brains. Collectively, these results strongly suggest that Cx43 controls immunoregulatory and angiogenic properties of astrocytes., Competing Interests: The authors declare no conflict of interest.

We present a new educational initiative called Meet-U that aims to train students for collaborative work in computational biology and to bridge the gap between education and research. Meet-U mimics the setup of collaborative research projects and takes advantage of the most popular tools for collaborative work and of cloud computing. Students are grouped in teams of 4-5 people and have to realize a project from A to Z that answers a challenging question in biology. Meet-U promotes "coopetition," as the students collaborate within and across the teams and are also in competition with each other to develop the best final product. Meet-U fosters interactions between different actors of education and research through the organization of a meeting day, open to everyone, where the students present their work to a jury of researchers and jury members give research seminars. This very unique combination of education and research is strongly motivating for the students and provides a formidable opportunity for a scientific community to unite and increase its visibility. We report on our experience with Meet-U in two French universities with master's students in bioinformatics and modeling, with protein-protein docking as the subject of the course. Meet-U is easy to implement and can be straightforwardly transferred to other fields and/or universities. All the information and data are available at www.meet-u.org.

Background/aims: Fatty acid oxidation (FAO), the main source of energy produced by tubular epithelial cells in the kidney, was found to be defective in tubulo-interstitial samples dissected out in kidney biopsies from patients with chronic kidney disease (CKD). Experimental data indicated that this decrease was a strong determinant of renal fibrogenesis, hence a focus for therapeutic interventions. Nevertheless, whether persistently differentiated renal tubules, surviving in a pro-fibrotic environment, also suffer from a decrease in FAO, is currently unknown., Methods: To address this question, we isolated proximal tubules captured ex vivo on the basis of the expression of an intact brush border antigen (Prominin-1) in C57BL6/J mice subjected to a controlled, two-hit model of renal fibrosis (reversible ischemic acute kidney injury (AKI) or sham surgery, followed by angiotensin 2 administration). A transcriptomic high throughput sequencing was performed on total mRNA from these cells, and on whole kidneys., Results: In contrast to mice subjected to sham surgery, mice with a history of AKI displayed histologically more renal fibrosis when exposed to angiotensin 2. High throughput RNA sequencing, principal component analysis and clustering showed marked consistency within experimental groups. As expected, FAO transcripts were decreased in whole fibrotic kidneys. Surprisingly, however, up- rather than down-regulation of metabolic pathways (oxidative phosphorylation, fatty acid metabolism, glycolysis, and PPAR signalling pathway) was a hallmark of the differentiated tubules captured from fibrotic kidneys. Immunofluorescence co-staining analysis confirmed that the expression of FAO enzymes was dependent of tubular trophicity., Conclusions: These data suggest that in differentiated proximal tubules energetic hyperactivity is promoted concurrently with organ fibrogenesis., (© 2018 The Author(s). Published by S. Karger AG, Basel.)

Background: The hydrolysis of biomass to simple sugars used for the production of biofuels in biorefineries requires the action of cellulolytic enzyme mixtures. During the last 50 years, the ascomycete Trichoderma reesei , the main source of industrial cellulase and hemicellulase cocktails, has been subjected to several rounds of classical mutagenesis with the aim to obtain higher production levels. During these random genetic events, strains unable to produce cellulases were generated. Here, whole genome sequencing and transcriptomic analyses of the cellulase-negative strain QM9978 were used for the identification of mutations underlying this cellulase-negative phenotype., Results: Sequence comparison of the cellulase-negative strain QM9978 to the reference strain QM6a identified a total of 43 mutations, of which 33 were located either close to or in coding regions. From those, we identified 23 single-nucleotide variants, nine InDels, and one translocation. The translocation occurred between chromosomes V and VII, is located upstream of the putative transcription factor vib1 , and abolishes its expression in QM9978 as detected during the transcriptomic analyses. Ectopic expression of vib1 under the control of its native promoter as well as overexpression of vib1 under the control of a strong constitutive promoter restored cellulase expression in QM9978, thus confirming that the translocation event is the reason for the cellulase-negative phenotype. Gene deletion of vib1 in the moderate producer strain QM9414 and in the high producer strain Rut-C30 reduced cellulase expression in both cases. Overexpression of vib1 in QM9414 and Rut-C30 had no effect on cellulase production, most likely because vib1 is already expressed at an optimal level under normal conditions., Conclusion: We were able to establish a link between a chromosomal translocation in QM9978 and the cellulase-negative phenotype of the strain. We identified the transcription factor vib1 as a key regulator of cellulases in T. reesei whose expression is absent in QM9978. We propose that in T. reesei , as in Neurospora crassa , vib1 is involved in cellulase induction, although the exact mechanism remains to be elucidated. The data presented here show an example of a combined genome sequencing and transcriptomic approach to explain a specific trait, in this case the QM9978 cellulase-negative phenotype, and how it helps to better understand the mechanisms during cellulase gene regulation. When focusing on mutations on the single base-pair level, changes on the chromosome level can be easily overlooked and through this work we provide an example that stresses the importance of the big picture of the genomic landscape during analysis of sequencing data.

Motivation: Data management and quality control of output from Illumina sequencers is a disk space- and time-consuming task. Thus, we developed Aozan to automatically handle data transfer, demultiplexing, conversion and quality control once a run has finished. This software greatly improves run data management and the monitoring of run statistics via automatic emails and HTML web reports., Availability and Implementation: Aozan is implemented in Java and Python, supported on Linux systems, and distributed under the GPLv3 License at: http://www.outils.genomique.biologie.ens.fr/aozan/ . Aozan source code is available on GitHub: https://github.com/GenomicParisCentre/aozan ., Contact: aozan@biologie.ens.fr., (© The Author (2017). Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com)

Astrocytes send out long processes that are terminated by endfeet at the vascular surface and regulate vascular functions as well as homeostasis at the vascular interface. To date, the astroglial mechanisms underlying these functions have been poorly addressed. Here we demonstrate that a subset of messenger RNAs is distributed in astrocyte endfeet. We identified, among this transcriptome, a pool of messenger RNAs bound to ribosomes, the endfeetome, that primarily encodes for secreted and membrane proteins. We detected nascent protein synthesis in astrocyte endfeet. Finally, we determined the presence of smooth and rough endoplasmic reticulum and the Golgi apparatus in astrocyte perivascular processes and endfeet, suggesting for local maturation of membrane and secreted proteins. These results demonstrate for the first time that protein synthesis occurs in astrocyte perivascular distal processes that may sustain their structural and functional polarization at the vascular interface., Competing Interests: The authors declare no conflict of interest.

The yeast Candida glabrata has become the second cause of systemic candidemia in humans. However, relatively few genome-wide studies have been conducted in this organism and our knowledge of its transcriptional regulatory network is quite limited. In the present work, we combined genome-wide chromatin immunoprecipitation (ChIP-seq), transcriptome analyses, and DNA binding motif predictions to describe the regulatory interactions of the seven Yap (Yeast AP1) transcription factors of C. glabrata. We described a transcriptional network containing 255 regulatory interactions and 309 potential target genes. We predicted with high confidence the preferred DNA binding sites for 5 of the 7 CgYaps and showed a strong conservation of the Yap DNA binding properties between S. cerevisiae and C. glabrata. We provided reliable functional annotation for 3 of the 7 Yaps and identified for Yap1 and Yap5 a core regulon which is conserved in S. cerevisiae, C. glabrata, and C. albicans. We uncovered new roles for CgYap7 in the regulation of iron-sulfur cluster biogenesis, for CgYap1 in the regulation of heme biosynthesis and for CgYap5 in the repression of GRX4 in response to iron starvation. These transcription factors define an interconnected transcriptional network at the cross-roads between redox homeostasis, oxygen consumption, and iron metabolism.

Background: Trichoderma reesei is the main industrial source of cellulases and hemicellulases required for the hydrolysis of biomass to simple sugars, which can then be used in the production of biofuels and biorefineries. The highly productive strains in use today were generated by classical mutagenesis. As byproducts of this procedure, mutants were generated that turned out to be unable to produce cellulases. In order to identify the mutations responsible for this inability, we sequenced the genome of one of these strains, QM9136, and compared it to that of its progenitor T. reesei QM6a., Results: In QM9136, we detected a surprisingly low number of mutagenic events in the promoter and coding regions of genes, i.e. only eight indels and six single nucleotide variants. One of these indels led to a frame-shift in the Zn₂Cys₆ transcription factor XYR1, the general regulator of cellulase and xylanase expression, and resulted in its C-terminal truncation by 140 amino acids. Retransformation of strain QM9136 with the wild-type xyr1 allele fully recovered the ability to produce cellulases, and is thus the reason for the cellulase-negative phenotype. Introduction of an engineered xyr1 allele containing the truncating point mutation into the moderate producer T. reesei QM9414 rendered this strain also cellulase-negative. The correspondingly truncated XYR1 protein was still able to enter the nucleus, but failed to be expressed over the basal constitutive level., Conclusion: The missing 140 C-terminal amino acids of XYR1 are therefore responsible for its previously observed auto-regulation which is essential for cellulases to be expressed. Our data present a working example of the use of genome sequencing leading to a functional explanation of the QM9136 cellulase-negative phenotype.

Background: The filamentous fungus Trichoderma reesei is the main industrial cellulolytic enzyme producer. Several strains have been developed in the past using random mutagenesis, and despite impressive performance enhancements, the pressure for low-cost cellulases has stimulated continuous research in the field. In this context, comparative study of the lower and higher producer strains obtained through random mutagenesis using systems biology tools (genome and transcriptome sequencing) can shed light on the mechanisms of cellulase production and help identify genes linked to performance. Previously, our group published comparative genome sequencing of the lower and higher producer strains NG 14 and RUT C30. In this follow-up work, we examine how these mutations affect phenotype as regards the transcriptome and cultivation behaviour., Results: We performed kinetic transcriptome analysis of the NG 14 and RUT C30 strains of early enzyme production induced by lactose using bioreactor cultivations close to an industrial cultivation regime. RUT C30 exhibited both earlier onset of protein production (3 h) and higher steady-state productivity. A rather small number of genes compared to previous studies were regulated (568), most of them being specific to the NG 14 strain (319). Clustering analysis highlighted similar behaviour for some functional categories and allowed us to distinguish between induction-related genes and productivity-related genes. Cross-comparison of our transcriptome data with previously identified mutations revealed that most genes from our dataset have not been mutated. Interestingly, the few mutated genes belong to the same clusters, suggesting that these clusters contain genes playing a role in strain performance., Conclusions: This is the first kinetic analysis of a transcriptomic study carried out under conditions approaching industrial ones with two related strains of T. reesei showing distinctive cultivation behaviour. Our study sheds some light on some of the events occurring in these strains following induction by lactose. The fact that few regulated genes have been affected by mutagenesis suggests that the induction mechanism is essentially intact compared to that for the wild-type isolate QM6a and might be engineered for further improvement of T. reesei. Genes from two specific clusters might be potential targets for such genetic engineering.

Peak calling is a critical step in ChIPseq data analysis. Choosing the correct algorithm as well as optimized parameters for a specific biological system is an essential task. In this article, we present an original peak-calling method (bPeaks) specifically designed to detect transcription factor (TF) binding sites in small eukaryotic genomes, such as in yeasts. As TF interactions with DNA are strong and generate high binding signals, bPeaks uses simple parameters to compare the sequences (reads) obtained from the immunoprecipitation (IP) with those from the control DNA (input). Because yeasts have small genomes (<20 Mb), our program has the advantage of using ChIPseq information at the single nucleotide level and can explore, in a reasonable computational time, results obtained with different sets of parameter values. Graphical outputs and text files are provided to rapidly assess the relevance of the detected peaks. Taking advantage of the simple promoter structure in yeasts, additional functions were implemented in bPeaks to automatically assign the peaks to promoter regions and retrieve peak coordinates on the DNA sequence for further predictions of regulatory motifs, enriched in the list of peaks. Applications of the bPeaks program to three different ChIPseq datasets from Saccharomyces cerevisiae, Candida albicans and Candida glabrata are presented. Each time, bPeaks allowed us to correctly predict the DNA binding sequence of the studied TF and provided relevant lists of peaks. The bioinformatics tool bPeaks is freely distributed to academic users. Supplementary data, together with detailed tutorials, are available online: http://bpeaks.gene-networks.net., (Copyright © 2014 John Wiley & Sons, Ltd.)

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

Background: Trichoderma is a genus of mycotrophic filamentous fungi (teleomorph Hypocrea) which possess a bright variety of biotrophic and saprotrophic lifestyles. The ability to parasitize and/or kill other fungi (mycoparasitism) is used in plant protection against soil-borne fungal diseases (biological control, or biocontrol). To investigate mechanisms of mycoparasitism, we compared the transcriptional responses of cosmopolitan opportunistic species and powerful biocontrol agents Trichoderma atroviride and T. virens with tropical ecologically restricted species T. reesei during confrontations with a plant pathogenic fungus Rhizoctonia solani., Results: The three Trichoderma spp. exhibited a strikingly different transcriptomic response already before physical contact with alien hyphae. T. atroviride expressed an array of genes involved in production of secondary metabolites, GH16 ß-glucanases, various proteases and small secreted cysteine rich proteins. T. virens, on the other hand, expressed mainly the genes for biosynthesis of gliotoxin, respective precursors and also glutathione, which is necessary for gliotoxin biosynthesis. In contrast, T. reesei increased the expression of genes encoding cellulases and hemicellulases, and of the genes involved in solute transport. The majority of differentially regulated genes were orthologues present in all three species or both in T. atroviride and T. virens, indicating that the regulation of expression of these genes is different in the three Trichoderma spp. The genes expressed in all three fungi exhibited a nonrandom genomic distribution, indicating a possibility for their regulation via chromatin modification., Conclusion: This genome-wide expression study demonstrates that the initial Trichoderma mycotrophy has differentiated into several alternative ecological strategies ranging from parasitism to predation and saprotrophy. It provides first insights into the mechanisms of interactions between Trichoderma and other fungi that may be exploited for further development of biofungicides.

DNA cytosine methylation is a widely conserved epigenetic mark in eukaryotes that appears to have critical roles in the regulation of genome structure and transcription. Genome-wide methylation maps have so far only been established from the supergroups Archaeplastida and Unikont. Here we report the first whole-genome methylome from a stramenopile, the marine model diatom Phaeodactylum tricornutum. Around 6% of the genome is intermittently methylated in a mosaic pattern. We find extensive methylation in transposable elements. We also detect methylation in over 320 genes. Extensive gene methylation correlates strongly with transcriptional silencing and differential expression under specific conditions. By contrast, we find that genes with partial methylation tend to be constitutively expressed. These patterns contrast with those found previously in other eukaryotes. By going beyond plants, animals and fungi, this stramenopile methylome adds significantly to our understanding of the evolution of DNA methylation in eukaryotes.

Unlabelled: We developed a modular and scalable framework called Eoulsan, based on the Hadoop implementation of the MapReduce algorithm dedicated to high-throughput sequencing data analysis. Eoulsan allows users to easily set up a cloud computing cluster and automate the analysis of several samples at once using various software solutions available. Our tests with Amazon Web Services demonstrated that the computation cost is linear with the number of instances booked as is the running time with the increasing amounts of data., Availability and Implementation: Eoulsan is implemented in Java, supported on Linux systems and distributed under the LGPL License at: http://transcriptome.ens.fr/eoulsan/

Chromodomains are found in many regulators of chromatin structure, and most of them recognize methylated lysines on histones. Here, we investigate the role of the Drosophila melanogaster protein Corto's chromodomain. The Enhancer of Trithorax and Polycomb Corto is involved in both silencing and activation of gene expression. Over-expression of the Corto chromodomain (CortoCD) in transgenic flies shows that it is a chromatin-targeting module, critical for Corto function. Unexpectedly, mass spectrometry analysis reveals that polypeptides pulled down by CortoCD from nuclear extracts correspond to ribosomal proteins. Furthermore, real-time interaction analyses demonstrate that CortoCD binds with high affinity RPL12 tri-methylated on lysine 3. Corto and RPL12 co-localize with active epigenetic marks on polytene chromosomes, suggesting that both are involved in fine-tuning transcription of genes in open chromatin. RNA-seq based transcriptomes of wing imaginal discs over-expressing either CortoCD or RPL12 reveal that both factors deregulate large sets of common genes, which are enriched in heat-response and ribosomal protein genes, suggesting that they could be implicated in dynamic coordination of ribosome biogenesis. Chromatin immunoprecipitation experiments show that Corto and RPL12 bind hsp70 and are similarly recruited on gene body after heat shock. Hence, Corto and RPL12 could be involved together in regulation of gene transcription. We discuss whether pseudo-ribosomal complexes composed of various ribosomal proteins might participate in regulation of gene expression in connection with chromatin regulators., Competing Interests: The authors have declared that no competing interests exist.

Background: The identification and characterization of the transcriptional regulatory networks governing the physiology and adaptation of microbial cells is a key step in understanding their behaviour. One such wide-domain regulatory circuit, essential to all cells, is carbon catabolite repression (CCR): it allows the cell to prefer some carbon sources, whose assimilation is of high nutritional value, over less profitable ones. In lower multicellular fungi, the C2H2 zinc finger CreA/CRE1 protein has been shown to act as the transcriptional repressor in this process. However, the complete list of its gene targets is not known., Results: Here, we deciphered the CRE1 regulatory range in the model cellulose and hemicellulose-degrading fungus Trichoderma reesei (anamorph of Hypocrea jecorina) by profiling transcription in a wild-type and a delta-cre1 mutant strain on glucose at constant growth rates known to repress and de-repress CCR-affected genes. Analysis of genome-wide microarrays reveals 2.8% of transcripts whose expression was regulated in at least one of the four experimental conditions: 47.3% of which were repressed by CRE1, whereas 29.0% were actually induced by CRE1, and 17.2% only affected by the growth rate but CRE1 independent. Among CRE1 repressed transcripts, genes encoding unknown proteins and transport proteins were overrepresented. In addition, we found CRE1-repression of nitrogenous substances uptake, components of chromatin remodeling and the transcriptional mediator complex, as well as developmental processes., Conclusions: Our study provides the first global insight into the molecular physiological response of a multicellular fungus to carbon catabolite regulation and identifies several not yet known targets in a growth-controlled environment.

The initial interaction between the Schwann cell and the axon is a complex and poorly understood aspect of the myelination process. To investigate the molecular mechanisms involved in this interaction and to identify novel genes required for myelination, we performed an RNA profiling analysis, comparing Schwann cells cultured alone or in the presence of neurons. This led to the selection of 30 genes, mostly upregulated on Schwann cell-axon interaction. Most of the identified proteins are associated with the extracellular space or signal transduction systems, consistent with possible roles in Schwann cell-axon interaction. We performed a functional analysis of one of these genes, Dok4 (downstream of kinase-4), which encodes a membrane-associated tyrosine kinase substrate. Silencing RNA-mediated knock-down of Dok4 severely affected in vitro myelination. Moreover, Dok4 is required at early stages in the myelination process, including the initial interaction with the axon, and is also involved in Schwann cell migration and proliferation. Finally, this analysis establishes the interest of our gene collection in further understanding of the molecular mechanisms involved in Schwann cell-axon interaction., (Copyright © 2010 Wiley-Liss, Inc.)

Due to its capacity to produce large amounts of cellulases, Trichoderma reesei is increasingly being investigated for second-generation biofuel production from lignocellulosic biomass. The induction mechanisms of T. reesei cellulases have been described recently, but the regulation of the genes involved in their transcription has not been studied thoroughly. Here we report the regulation of expression of the two activator genes xyr1 and ace2, and the corepressor gene ace1, during the induction of cellulase biosynthesis by the inducer lactose in T. reesei QM 9414, a strain producing low levels of cellulase (low producer). We show that all three genes are induced by lactose. xyr1 was also induced by d-galactose, but this induction was independent of d-galactose metabolism. Moreover, ace1 was carbon catabolite repressed, whereas full induction of xyr1 and ace2 in fact required CRE1. Significant differences in these regulatory patterns were observed in the high-producer strain RUT C30 and the hyperproducer strain T. reesei CL847. These observations suggest that a strongly elevated basal transcription level of xyr1 and reduced upregulation of ace1 by lactose may have been important for generating the hyperproducer strain and that thus, these genes are major control elements of cellulase production.

Background: Mycoparasitism, a lifestyle where one fungus is parasitic on another fungus, has special relevance when the prey is a plant pathogen, providing a strategy for biological control of pests for plant protection. Probably, the most studied biocontrol agents are species of the genus Hypocrea/Trichoderma., Results: Here we report an analysis of the genome sequences of the two biocontrol species Trichoderma atroviride (teleomorph Hypocrea atroviridis) and Trichoderma virens (formerly Gliocladium virens, teleomorph Hypocrea virens), and a comparison with Trichoderma reesei (teleomorph Hypocrea jecorina). These three Trichoderma species display a remarkable conservation of gene order (78 to 96%), and a lack of active mobile elements probably due to repeat-induced point mutation. Several gene families are expanded in the two mycoparasitic species relative to T. reesei or other ascomycetes, and are overrepresented in non-syntenic genome regions. A phylogenetic analysis shows that T. reesei and T. virens are derived relative to T. atroviride. The mycoparasitism-specific genes thus arose in a common Trichoderma ancestor but were subsequently lost in T. reesei., Conclusions: The data offer a better understanding of mycoparasitism, and thus enforce the development of improved biocontrol strains for efficient and environmentally friendly protection of plants., (© 2011 Kubicek et al.; licensee BioMed Central Ltd.)

Despite the development of new high-throughput sequencing techniques, microarrays are still attractive tools to study small genome organisms, thanks to sample multiplexing and high-feature densities. However, the oligonucleotide design remains a delicate step for most users. A vast array of software is available to deal with this problem, but each program is developed with its own strategy, which makes the choice of the best solution difficult. Here we describe Teolenn, a universal probe design workflow developed with a flexible and customizable module organization allowing fixed or variable length oligonucleotide generation. In addition, our software is able to supply quality scores for each of the designed probes. In order to assess the relevance of these scores, we performed a real hybridization using a tiling array designed against the Trichoderma reesei fungus genome. We show that our scoring pipeline correlates with signal quality for 97.2% of all the designed probes, allowing for a posteriori comparisons between quality scores and signal intensities. This result is useful in discarding any bad scoring probes during the design step in order to get high-quality microarrays. Teolenn is available at http://transcriptome.ens.fr/teolenn/.

Background: Diatoms are largely responsible for production of biogenic silica in the global ocean. However, in surface seawater, Si(OH)(4) can be a major limiting factor for diatom productivity. Analyzing at the global scale the genes networks involved in Si transport and metabolism is critical in order to elucidate Si biomineralization, and to understand diatoms contribution to biogeochemical cycles., Methodology/principal Findings: Using whole genome expression analyses we evaluated the transcriptional response to Si availability for the model species Phaeodactylum tricornutum. Among the differentially regulated genes we found genes involved in glutamine-nitrogen pathways, encoding putative extracellular matrix components, or involved in iron regulation. Some of these compounds may be good candidates for intracellular intermediates involved in silicic acid storage and/or intracellular transport, which are very important processes that remain mysterious in diatoms. Expression analyses and localization studies gave the first picture of the spatial distribution of a silicic acid transporter in a diatom model species, and support the existence of transcriptional and post-transcriptional regulations., Conclusions/significance: Our global analyses revealed that about one fourth of the differentially expressed genes are organized in clusters, underlying a possible evolution of P. tricornutum genome, and perhaps other pennate diatoms, toward a better optimization of its response to variable environmental stimuli. High fitness and adaptation of diatoms to various Si levels in marine environments might arise in part by global regulations from gene (expression level) to genomic (organization in clusters, dosage compensation by gene duplication), and by post-transcriptional regulation and spatial distribution of SIT proteins., Competing Interests: The authors have declared that no competing interests exist.

Neural crest (NC) cells are a multipotent, highly migratory cell population that generates most of the components of the peripheral nervous system (PNS), including the glial Schwann cells (SC) and boundary cap (BC) cells. These latter cells are located at the interface between the central nervous system and PNS, at the exit/entry points of ventral motor/dorsal sensory axons and give rise to all SC in the nerve roots and to a subset of nociceptive neurons and satellite cells in the dorsal root ganglia. In the present study we have compared BC cells with two closely related cell types, NC and Schwann cell precursors (SCpr), by RNA profiling. This led to the definition of a set of 10 genes that show specific expression in BC cells and/or in their derivatives along the nerve roots. Analysis of the expression of these genes during mouse development revealed novel features, of those most important are: (i) dorsal and ventral nerve root BC cell derivatives express different sets of genes, suggesting that they have distinct properties; (ii) these cells undergo major modifications in their gene expression pattern between embryonic days 14.5 and 17.5, possibly linked to the SCpr-immature Schwann cell transition; (iii) nerve roots SC differ from more distal SC not only in their origins and locations, but also in their gene expression patterns. In conclusion, the identification of these novel makers opens the way for a detailed characterization of BC cells in both mouse and man., ((c) 2009 Wiley-Liss, Inc.)

Trichoderma reesei (teleomorph Hypocrea jecorina) is the main industrial source of cellulases and hemicellulases harnessed for the hydrolysis of biomass to simple sugars, which can then be converted to biofuels such as ethanol and other chemicals. The highly productive strains in use today were generated by classical mutagenesis. To learn how cellulase production was improved by these techniques, we performed massively parallel sequencing to identify mutations in the genomes of two hyperproducing strains (NG14, and its direct improved descendant, RUT C30). We detected a surprisingly high number of mutagenic events: 223 single nucleotides variants, 15 small deletions or insertions, and 18 larger deletions, leading to the loss of more than 100 kb of genomic DNA. From these events, we report previously undocumented non-synonymous mutations in 43 genes that are mainly involved in nuclear transport, mRNA stability, transcription, secretion/vacuolar targeting, and metabolism. This homogeneity of functional categories suggests that multiple changes are necessary to improve cellulase production and not simply a few clear-cut mutagenic events. Phenotype microarrays show that some of these mutations result in strong changes in the carbon assimilation pattern of the two mutants with respect to the wild-type strain QM6a. Our analysis provides genome-wide insights into the changes induced by classical mutagenesis in a filamentous fungus and suggests areas for the generation of enhanced T. reesei strains for industrial applications such as biofuel production.

The increase in feature resolution and the availability of multipack formats from microarray providers has opened the way to various custom genomic applications. However, oligonucleotide design and selection remains a bottleneck of the microarray workflow. Several tools are available to perform this work, and choosing the best one is not an easy task, nor are the choices obvious. Here we review the oligonucleotide design field to help users make their choice. We have first performed a comparative evaluation of the available solutions based on a set of criteria including: ease of installation, user-friendly access, the number of parameters and settings available. In a second step, we chose to submit two real cases to a selection of programs. Finally, we used a set of tests for the in silico benchmark of the oligo sets obtained from each type of software. We show that the design software must be selected according to the goal of the scientist, depending on factors such as the organism used, the number of probes required and their localization on the target sequence. The present work provides keys to the choice of the most relevant software, according to the various parameters we tested.

Background: Branchiomotor neurons comprise an important class of cranial motor neurons that innervate the branchial-arch-derived muscles of the face, jaw and neck. They arise in the ventralmost progenitor domain of the rhombencephalon characterized by expression of the homeodomain transcription factors Nkx2.2 and Phox2b. Phox2b in particular plays a key role in the specification of branchiomotor neurons. In its absence, generic neuronal differentiation is defective in the progenitor domain and no branchiomotor neurons are produced. Conversely, ectopic expression of Phox2b in spinal regions of the neural tube promotes cell cycle exit and neuronal differentiation and, at the same time, induces genes and an axonal phenotype characteristic for branchiomotor neurons. How Phox2b exerts its pleiotropic functions, both as a proneural gene and a neuronal subtype determinant, has remained unknown., Results: To gain further insights into the genetic program downstream of Phox2b, we searched for novel Phox2b-regulated genes by cDNA microarray analysis of facial branchiomotor neuron precursors from heterozygous and homozygous Phox2b mutant embryos. We selected for functional studies the genes encoding the axonal growth promoter Gap43, the Wnt antagonist Sfrp1 and the transcriptional regulator Sox13, which were not previously suspected to play roles downstream of Phox2b and whose expression was affected by Phox2b misexpression in the spinal cord. While Gap43 did not produce an obvious phenotype when overexpressed in the neural tube, Sfrp1 induced the interneuron marker Lhx1,5 and Sox13 inhibited neuronal differentiation. We then tested whether Sfrp1 and Sox13, which are down-regulated by Phox2b in the facial neuron precursors, would antagonize some aspects of Phox2b activity. Co-expression of Sfrp1 prevented Phox2b from repressing Lhx1,5 and alleviated the commissural axonal phenotype. When expressed together with Sox13, Phox2b was still able to promote cell cycle exit and neuronal differentiation, but the cells failed to relocate to the mantle layer and to extinguish the neural stem cell marker Sox2., Conclusion: Our results suggest novel roles for Sfrp1 and Sox13 in neuronal subtype specification and generic neuronal differentiation, respectively, and indicate that down-regulation of Sfrp1 and Sox13 are essential aspects of the genetic program controlled by Phox2b in cranial motoneurons.

The asymmetric localization of mRNA plays an important role in coordinating posttranscriptional events in eukaryotic cells. We investigated the peripheral mitochondrial localization of nuclear-encoded mRNAs (MLR) in various conditions in which the mRNA binding protein context and the translation efficiency were altered. We identified Puf3p, a Pumilio family RNA-binding protein, as the first trans-acting factor controlling the MLR phenomenon. This allowed the characterization of two classes of genes whose mRNAs are translated to the vicinity of mitochondria. Class I mRNAs (256 genes) have a Puf3p binding motif in their 3'UTR region and many of them have their MLR properties deeply affected by PUF3 deletion. Conversely, mutations in the Puf3p binding motif alter the mitochondrial localization of BCS1 mRNA. Class II mRNAs (224 genes) have no Puf3p binding site and their asymmetric localization is not affected by the absence of PUF3. In agreement with a co-translational import process, we observed that the presence of puromycin loosens the interactions between most of the MLR-mRNAs and mitochondria. Unexpectedly, cycloheximide, supposed to solidify translational complexes, turned out to destabilize a class of mRNA-mitochondria interactions. Classes I and II mRNAs, which are therefore transported to the mitochondria through different pathways, correlated with different functional modules. Indeed, Class I genes code principally for the assembly factors of respiratory chain complexes and the mitochondrial translation machinery (ribosomes and translation regulators). Class II genes encode proteins of the respiratory chain or proteins involved in metabolic pathways. Thus, MLR, which is intimately linked to translation control, and the activity of mRNA-binding proteins like Puf3p, may provide the conditions for a fine spatiotemporal control of mitochondrial protein import and mitochondrial protein complex assembly. This work therefore provides new openings for the global study of mitochondria biogenesis.