Growth-rate-dependent and nutrient-specific gene expression resource allocation in fission yeast.

Cellular resources are limited and their relative allocation to gene expression programmes determines physiological states and global properties such as the growth rate. Here, we determined the importance of the growth rate in explaining relative changes in protein and mRNA levels in the simple eukaryote Schizosaccharomyces pombe grown on non-limiting nitrogen sources. Although expression of half of fission yeast genes was significantly correlated with the growth rate, this came alongside wide-spread nutrient-specific regulation. Proteome and transcriptome often showed coordinated regulation but with notable exceptions, such as metabolic enzymes. Genes positively correlated with growth rate participated in every level of protein production apart from RNA polymerase II-dependent transcription. Negatively correlated genes belonged mainly to the environmental stress response programme. Critically, metabolic enzymes, which represent ∼55-70% of the proteome by mass, showed mostly condition-specific regulation. In summary, we provide a rich account of resource allocation to gene expression in a simple eukaryote, advancing our basic understanding of the interplay between growth-rate-dependent and nutrient-specific gene expression.

bayNorm: Bayesian gene expression recovery, imputation and normalization for single-cell RNA-sequencing data.

MOTIVATION: Normalization of single-cell RNA-sequencing (scRNA-seq) data is a prerequisite to their interpretation. The marked technical variability, high amounts of missing observations and batch effect typical of scRNA-seq datasets make this task particularly challenging. There is a need for an efficient and unified approach for normalization, imputation and batch effect correction. RESULTS: Here, we introduce bayNorm, a novel Bayesian approach for scaling and inference of scRNA-seq counts. The method's likelihood function follows a binomial model of mRNA capture, while priors are estimated from expression values across cells using an empirical Bayes approach. We first validate our assumptions by showing this model can reproduce different statistics observed in real scRNA-seq data. We demonstrate using publicly available scRNA-seq datasets and simulated expression data that bayNorm allows robust imputation of missing values generating realistic transcript distributions that match single molecule fluorescence in situ hybridization measurements. Moreover, by using priors informed by dataset structures, bayNorm improves accuracy and sensitivity of differential expression analysis and reduces batch effect compared with other existing methods. Altogether, bayNorm provides an efficient, integrated solution for global scaling normalization, imputation and true count recovery of gene expression measurements from scRNA-seq data. AVAILABILITY AND IMPLEMENTATION: The R package 'bayNorm' is publishd on bioconductor at https://bioconductor.org/packages/release/bioc/html/bayNorm.html. The code for analyzing data in this article is available at https://github.com/WT215/bayNorm_papercode. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single-cell imaging and RNA sequencing reveal patterns of gene expression heterogeneity during fission yeast growth and adaptation.

Phenotypic cell-to-cell variability is a fundamental determinant of microbial fitness that contributes to stress adaptation and drug resistance. Gene expression heterogeneity underpins this variability but is challenging to study genome-wide. Here we examine the transcriptomes of >2,000 single fission yeast cells exposed to various environmental conditions by combining imaging, single-cell RNA sequencing and Bayesian true count recovery. We identify sets of highly variable genes during rapid proliferation in constant culture conditions. By integrating single-cell RNA sequencing and cell-size data, we provide insights into genes that are regulated during cell growth and division, including genes whose expression does not scale with cell size. We further analyse the heterogeneity of gene expression during adaptive and acute responses to changing environments. Entry into the stationary phase is preceded by a gradual, synchronized adaptation in gene regulation that is followed by highly variable gene expression when growth decreases. Conversely, sudden and acute heat shock leads to a stronger, coordinated response and adaptation across cells. This analysis reveals that the magnitude of global gene expression heterogeneity is regulated in response to different physiological conditions within populations of a unicellular eukaryote.

The TAF9 C-terminal conserved region domain is required for SAGA and TFIID promoter occupancy to promote transcriptional activation.

A common function of the TFIID and SAGA complexes, which are recruited by transcriptional activators, is to deliver TBP to promoters to stimulate transcription. Neither the relative contributions of the five shared TBP-associated factor (TAF) subunits in TFIID and SAGA nor the requirement for different domains in shared TAFs for transcriptional activation is well understood. In this study, we uncovered the essential requirement for the highly conserved C-terminal region (CRD) of Taf9, a shared TAF, for transcriptional activation in yeast. Transcriptome profiling performed under Gcn4-activating conditions showed that the Taf9 CRD is required for induced expression of ∼9% of the yeast genome. The CRD was not essential for the Taf9-Taf6 interaction, TFIID or SAGA integrity, or Gcn4 interaction with SAGA in cell extracts. Microarray profiling of a SAGA mutant (spt20Δ) yielded a common set of genes induced by Spt20 and the Taf9 CRD. Chromatin immunoprecipitation (ChIP) assays showed that, although the Taf9 CRD mutation did not impair Gcn4 occupancy, the occupancies of TFIID, SAGA, and the preinitiation complex were severely impaired at several promoters. These results suggest a crucial role for the Taf9 CRD in genome-wide transcription and highlight the importance of conserved domains, other than histone fold domains, as a common determinant for TFIID and SAGA functions.

Set2 methylation of histone H3 lysine 36 suppresses histone exchange on transcribed genes.

Set2-mediated methylation of histone H3 at Lys 36 (H3K36me) is a co-transcriptional event that is necessary for the activation of the Rpd3S histone deacetylase complex, thereby maintaining the coding region of genes in a hypoacetylated state. In the absence of Set2, H3K36 or Rpd3S acetylated histones accumulate on open reading frames (ORFs), leading to transcription initiation from cryptic promoters within ORFs. Although the co-transcriptional deacetylation pathway is well characterized, the factors responsible for acetylation are as yet unknown. Here we show that, in yeast, co-transcriptional acetylation is achieved in part by histone exchange over ORFs. In addition to its function of targeting and activating the Rpd3S complex, H3K36 methylation suppresses the interaction of H3 with histone chaperones, histone exchange over coding regions and the incorporation of new acetylated histones. Thus, Set2 functions both to suppress the incorporation of acetylated histones and to signal for the deacetylation of these histones in transcribed genes. By suppressing spurious cryptic transcripts from initiating within ORFs, this pathway is essential to maintain the accuracy of transcription by RNA polymerase II.

Characterization of Leishmania donovani aquaporins shows presence of subcellular aquaporins similar to tonoplast intrinsic proteins of plants.

Leishmania donovani, a protozoan parasite, resides in the macrophages of the mammalian host. The aquaporin family of proteins form important components of the parasite-host interface. The parasite-host interface could be a potential target for chemotherapy. Analysis of L. major and L. infantum genomes showed the presence of five aquaporins (AQPs) annotated as AQP9 (230aa), AQP putative (294aa), AQP-like protein (279aa), AQP1 (314aa) and AQP-like protein (596aa). We report here the structural modeling, localization and functional characterization of the AQPs from L. donovani. LdAQP1, LdAQP9, LdAQP2860 and LdAQP2870 have the canonical NPA-NPA motifs, whereas LdAQP putative has a non-canonical NPM-NPA motif. In the carboxyl terminal to the second NPA box of all AQPs except AQP1, a valine/alanine residue was found instead of the arginine. In that respect these four AQPs are similar to tonoplast intrinsic proteins in plants, which are localized to intracellular organelles. Confocal microscopy of L. donovani expressing GFP-tagged AQPs showed an intracellular localization of LdAQP9 and LdAQP2870. Real-time PCR assays showed expression of all aquaporins except LdAQP2860, whose level was undetectable. Three-dimensional homology modeling of the AQPs showed that LdAQP1 structure bears greater topological similarity to the aquaglyceroporin than to aquaporin of E. coli. The pore of LdAQP1 was very different from the rest in shape and size. The cavity of LdAQP2860 was highly irregular and undefined in geometry. For functional characterization, four AQP proteins were heterologously expressed in yeast. In the fps1Δ yeast cells, which lacked the key aquaglyceroporin, LdAQP1 alone displayed an osmosensitive phenotype indicating glycerol transport activity. However, expression of LdAQP1 and LdAQP putative in a yeast gpd1Δ strain, deleted for glycerol production, conferred osmosensitive phenotype indicating water transport activity or aquaporin function. Our analysis for the first time shows the presence of subcellular aquaporins and provides structural and functional characterization of aquaporins in Leishmania donovani.