Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

BACKGROUND

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and conditions. How the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear.

RESULTS

We analyze existing bulk and single-cell transcriptomic datasets to uncover the prevalence of transcriptional adaptation in mammalian systems across diverse contexts and cell types. We perform regulon gene expression analyses of transcription factor target sets in both bulk and pooled single-cell genetic perturbation datasets. Our results reveal greater robustness in expression of regulons of transcription factors exhibiting transcriptional adaptation compared to those of transcription factors that do not. Stochastic mathematical modeling of minimal compensatory gene networks qualitatively recapitulates several aspects of transcriptional adaptation, including paralog upregulation and robustness to mutation. Combined with machine learning analysis of network features of interest, our framework offers potential explanations for which regulatory steps are most important for transcriptional adaptation.

CONCLUSIONS

Our integrative approach identifies several putative hits-genes demonstrating possible transcriptional adaptation-to follow-up on experimentally and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Design principles of multi-map variation in biological systems

Complexity in biology is often described using a multi-map hierarchical architecture, where the genotype, representing the encoded information, is mapped to the functional level, known as the phenotype, which is then connected to a latent phenotype we refer to as fitness. This underlying architecture governs the processes driving evolution. Furthermore, natural selection, along with other neutral forces, can, in turn, modify these maps. At each level, variation is observed. Here, I propose the need to establish principles that can aid in understanding the transformation of variation within this multi-map architecture. Specifically, I will introduce three, related to the presence of modulators, constraints, and the modular channeling of variation. By comprehending these design principles in various biological systems, we can gain better insights into the mechanisms underlying these maps and how they ultimately contribute to evolutionary dynamics.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene can trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, to date, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and contexts. Moreover, how the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear. Here we combine computational analysis of existing datasets with stochastic mathematical modeling and machine learning to uncover the widespread prevalence of transcriptional adaptation in mammalian systems and the diverse single-cell manifestations of minimal compensatory gene networks. Regulon gene expression analysis of a pooled single-cell genetic perturbation dataset recapitulates important model predictions. Our integrative approach uncovers several putative hits-genes demonstrating possible transcriptional adaptation-to follow up on experimentally, and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralogs TDH1 and TDH2. TDH1 and TDH2 are upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

The role of the Aspergillus nidulans high mobility group B protein HmbA, the orthologue of Saccharomyces cerevisiae Nhp6p

The mammalian HMGB1 is a high-mobility-group B protein, which is both an architectural and functional element of chromatin. Nhp6p, the extensively studied fungal homologue of HMGB1 in Saccharomyces cerevisiae has pleiotropic physiological functions. Despite the existence of Nhp6p orthologues in filamentous ascomycetes, little is known about their physiological roles besides their contribution to sexual development. Here we study the function of HmbA, the Aspergillus nidulans orthologue of Nhp6p. We show that HmbA influences the utilization of various carbon- and nitrogen sources, stress tolerance, secondary metabolism, hyphae elongation and maintenance of polarized growth. Additionally, by conducting heterologous expression studies, we demonstrate that HmbA and Nhp6p are partially interchangeable. HmbA restores SNR6 transcription and fitness of nhp6AΔBΔ mutant and reverses its heat sensitivity. Nhp6Ap complements several phenotypes of hmbAΔ, including ascospore formation, utilization of various carbon- and nitrogen-sources, radial growth rate, hypha elongation by polarized growth. However, Nhp6Ap does not complement sterigmatocystin production in a hmbAΔ strain. Finally, we also show that HmbA is necessary for the normal expression of the endochitinase chiA, a cell wall re-modeller that is pivotal for the normal mode of maintenance of polar growth.

Impairment of transcription factor Gcr1p binding motif perturbs OPI3 transcription in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, the transcription factor GCR1 plays a vital role in carbohydrate metabolism and in the current study we tried to elucidate its role in lipid metabolism. In silico analysis revealed the upstream activation sequence (UAS) in the promoter region of OPI3 possessed six conserved recognition sequences for Gcr1p and the ChIP assay confirmed the binding of Gcr1p on the OPI3 promoter region. The real-time quantitative polymerase chain reaction and promoter-reporter activity revealed a substantial reduction in OPI3 expression and was supported with decreased phosphatidylcholine (PC) level that is rescued with exogenous choline supplementation in gcr1∆ cells. Simultaneously, there was an increase in triacylglycerol level, accompanied with increased number and size of lipid droplets in gcr1∆ cells. The expression of pT1, pT2 truncations in opi3∆ cells revealed the -1 to -500 bp in the promoter region is essential for the activation of OPI3 transcription. The mutation specifically at UAS_{CT box} (-265) in the OPI3 promoter region displayed a reduction in the PC level and the additional mutation at UAS_INO (-165) further reduced the PC level. Collectively, our data suggest that the GCR1 transcription factor also regulates the OPI3 expression and has an impact on lipid homeostasis.

Convergent phenotypic evolution towards fosfomycin collateral sensitivity of Pseudomonas aeruginosa antibiotic-resistant mutants

The rise of antibiotic resistance and the reduced amount of novel antibiotics support the need of developing novel strategies to fight infections, based on improving the use of the antibiotics we already have. Collateral sensitivity is an evolutionary trade-off associated with the acquisition of antibiotic resistance that can be exploited to tackle this relevant health problem. However, different works have shown that patterns of collateral sensitivity are not always conserved, thus precluding the exploitation of this evolutionary trade-off to fight infections. In this work, we identify a robust pattern of collateral sensitivity to fosfomycin in Pseudomonas aeruginosa antibiotic-resistant mutants, selected by antibiotics belonging to different structural families. We characterize the underlying mechanism of the collateral sensitivity observed, which is a reduced expression of the genes encoding the peptidoglycan-recycling pathway, which preserves the peptidoglycan synthesis in situations where its de novo synthesis is blocked, and a reduced expression of fosA, encoding a fosfomycin-inactivating enzyme. We propose that the identification of robust collateral sensitivity patterns, as well as the understanding of the molecular mechanisms behind these phenotypes, would provide valuable information to design evolution-based strategies to treat bacterial infections.