The function and properties of the Azf1 transcriptional regulator change with growth conditions in Saccharomyces cerevisiae

The pleiotropic phenotype of FlbA of Aspergillus niger is explained in part by the activity of seven of its downstream-regulated transcription factors

Inactivation of flbA in Aspergillus niger results in thinner cell walls, increased cell lysis, abolished sporulation, and an increased secretome complexity. A total of 36 transcription factor (TF) genes are differentially expressed in ΔflbA. Here, seven of these genes (abaA, aslA, aslB, azf1, htfA, nosA, and srbA) were inactivated. Inactivation of each of these genes affected sporulation and, with the exception of abaA, cell wall integrity and protein secretion. The impact on secretion was strongest in the case of ΔaslA and ΔaslB that showed increased pepsin, cellulase, and amylase activity. Biomass was reduced of agar cultures of ΔabaA, ΔaslA, ΔnosA, and ΔsrbA, while biomass was higher in liquid shaken cultures of ΔaslA and ΔaslB. The ΔaslA and ΔhtfA strains showed increased resistance to H2O2, while ΔaslB was more sensitive to this reactive oxygen species. Together, inactivation of the seven TF genes impacted biomass formation, sporulation, protein secretion, and stress resistance, and thereby these genes explain at least part of the pleiotropic phenotype of ΔflbA of A. niger.

Systems level analysis of time and stimuli specific signaling through PKA

It is well known that eukaryotic cells create gradients of cAMP across space and time to regulate the cAMP dependent protein kinase (PKA) and, in turn, growth and metabolism. However, it is unclear how PKA responds to different concentrations of cAMP. Here, to address this question, we examine PKA signaling in Saccharomyces cerevisiae in different conditions, timepoints, and concentrations of the chemical inhibitor 1-NM-PP1, using phosphoproteomics. These experiments show that there are numerous proteins that are only phosphorylated when cAMP and PKA activity are at/near their maximum level, while other proteins are phosphorylated even when cAMP levels and PKA activity are low. The data also show that PKA drives cells into distinct growth states by acting on proteins with different thresholds for phosphorylation in different conditions. Analysis of the sequences surrounding the 118 PKA-dependent phosphosites suggests that the phosphorylation thresholds are set, at least in part, by the affinity of PKA for each site.

Overexpression of the Transcription Factor Azf1 Reveals Novel Regulatory Functions and Impacts β-Glucosidase Production in Trichoderma reesei

The fungus Trichoderma reesei is an essential producer of enzymes that degrade lignocellulosic biomass to produce value-added bioproducts. The cellulolytic system of T. reesei is controlled by several transcription factors (TFs) that efficiently regulate the production of these enzymes. Recently, a new TF named Azf1 was identified as a positive regulator of cellulase expression. Here, we investigated novel regulatory functions of Azf1 by its overexpression. In the mutant strain OEazf1, overexpression of azf1 was achieved under both repression and induction conditions. Although azf1 was more abundant in transcript and protein, overexpression of this TF did not activate transcription of the cellulase gene in the presence of the repressor glucose, suggesting that Azf1 may be subject to posttranslational regulation. In cellulose, the expression of swo, encoding the accessory protein swollenin, and the β-glucosidases cel1a, cel1b, cel3b, and cel3g increases in the early stages of cultivation. The increased production of these β-glucosidases increases the hydrolysis rate of cellobiose and sophorose, which activates carbon catabolite repression (CCR) and causes repression of cellulase genes and the regulator Xyr1 in the later stages of cultivation. Moreover, overexpression of azf1 led to increased cellulase activity in T. reesei during long-term cultivation in cellulose and sugarcane bagasse. Our results provide new insights into the mechanisms regulating Azf1 and novel genes that are important targets of this TF. This work contributes to a better understanding of the complex mechanisms regulating cellulase expression in T. reesei. It will contribute to the development of strains with higher production of these essential enzymes.

The intricate role of Sir2 in oxidative stress response during the post-diauxic phase in Saccharomyces cerevisiae

Silent information regulator 2 (Sir2) is a conserved NAD+-dependent histone deacetylase crucial for regulating cellular stress response and the aging process in Saccharomyces cerevisiae. In this study, we investigated the molecular mechanism underlying how the absence of Sir2 can lead to altered stress susceptibilities in S. cerevisiae under different environmental and physiological conditions. In a glucose-complex medium, the sir2Δ strain showed increased sensitivity to H2O2 compared to the wild-type strain during the post-diauxic phase. In contrast, it displayed increased resistance during the exponential growth phase. Transcriptome analysis of yeast cells in the post-diauxic phase indicated that the sir2Δ mutant expressed several oxidative defense genes at lower levels than the wild-type, potentially accounting for its increased susceptibility to H2O2. Interestingly, however, the sir2Δras2Δ double mutant exhibited greater resistance to H2O2 than the ras2Δ single mutant counterpart. We found that the expression regulation of the cytoplasmic catalase encoded by CTT1 was critical for the increased resistance to H2O2 in the sir2Δras2Δ strain. The expression of the CTT1 gene was influenced by the combined effect of RAS2 deletion and the transcription factor Azf1, whose level was modulated by Sir2. These findings provide insights into the importance of understanding the intricate interactions among various factors contributing to cellular stress response.

FUN-PROSE: A deep learning approach to predict condition-specific gene expression in fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE-a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

AnAzf1 acts as a positive regulator of ochratoxin A biosynthesis in Aspergillus niger

Aspergillus niger produces genotoxic and carcinogenic ochratoxin A (OTA) that severely threatens human and animal health. Transcription factor Azf1 is essential in regulating fungal cell development and primary metabolism. However, its effect and mechanism on secondary metabolism are unclear. Here, we characterized and deleted a Azf1 homolog gene, An15g00120 (AnAzf1), in A. niger, which completely blocked OTA production, and repressed the OTA cluster genes, p450, nrps, hal, and bzip at the transcriptional level. The results indicated that AnAzf1 was a positive regulator of OTA biosynthesis. Transcriptome sequencing results showed that the AnAzf1 deletion significantly upregulated antioxidant genes and downregulated oxidative phosphorylation genes. Enzymes involved in reactive oxygen species (ROS) scavenging, including catalase (CAT) and peroxidase (POD) were increased, and the corresponding ROS levels were decreased. Upregulation of genes (cat, catA, hog1, and gfd) in the MAPK pathway and downregulation of genes in iron homeostasis were associated with decreased ROS levels, linking the altered MAPK pathway and iron homeostasis to lower ROS levels caused by AnAzf1 deletion. Additionally, enzymes including complex I (NADH-ubiquinone oxidoreductase), and complex V (ATP synthase), as well as ATP levels, were significantly decreased, indicating impaired oxidative phosphorylation caused by the AnAzf1-deletion. During lower ROS levels and impaired oxidative phosphorylation, OTA was not produced in ∆AnAzf1. Together, these results strongly suggested that AnAzf1 deletion blocked OTA production in A. niger by a synergistic interference of ROS accumulation and oxidative phosphorylation. KEY POINTS: • AnAzf1 positively regulated OTA biosynthesis in A. niger. • Deletion of AnAzf1 decreased ROS levels and impaired oxidative phosphorylation. • An altered MAPK pathway and iron homeostasis were associated with lower ROS levels.

The Putative C2H2 Transcription Factor VadH Governs Development, Osmotic Stress Response, and Sterigmatocystin Production in Aspergillus nidulans

The VosA-VelB hetero-dimeric complex plays a pivotal role in regulating development and secondary metabolism in Aspergillus nidulans. In this work, we characterize a new VosA/VelB-activated gene called vadH, which is predicted to encode a 457-amino acid length protein containing four adjacent C₂H₂ zinc-finger domains. Mutational inactivation of vosA or velB led to reduced mRNA levels of vadH throughout the lifecycle, suggesting that VosA and VelB have a positive regulatory effect on the expression of vadH. The deletion of vadH resulted in decreased asexual development (conidiation) but elevated production of sexual fruiting bodies (cleistothecia), indicating that VadH balances asexual and sexual development in A. nidulans. Moreover, the vadH deletion mutant exhibited elevated susceptibility to hyperosmotic stress compared to wild type and showed elevated production of the mycotoxin sterigmatocystin (ST). Genome-wide expression analyses employing RNA-Seq have revealed that VadH is likely involved in regulating more genes and biological pathways in the developmental stages than those in the vegetative growth stage. The brlA, abaA, and wetA genes of the central regulatory pathway for conidiation are downregulated significantly in the vadH null mutant during asexual development. VadH also participates in regulating the genes, mat2, ppgA and lsdA, etc., related to sexual development, and some of the genes in the ST biosynthetic gene cluster. In summary, VadH is a putative transcription factor with four C₂H₂ finger domains and is involved in regulating asexual/sexual development, osmotic stress response, and ST production in A. nidulans.

Three-Dimensional Genome Map of the Filamentous Fungus Penicillium oxalicum

Higher-order spatial organization of the chromatin in the nucleus plays crucial roles in the maintenance of cell functions and the regulation of gene expression. Three-dimensional (3D) genome sequencing has been used to great effect in mammal and plants, but the availability of 3D genomes of filamentous fungi is severely limited. Here, we performed a chromosome-level genome assembly of Penicillium oxalicum through single-molecule real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C), with a scaffold N50 of 4.07 Mb and a contig N50 of 3.81 Mb, and further elucidated the 3D genome architecture of P. oxalicum. High-frequency interchromosomal contacts occurred within the centromeres and telomeres, as well as within individual chromosomes. There were 12,203 cis-interactions and 7,884 trans-interactions detected at a resolution of 1 kb. Moreover, a total of 1,099 topologically associated domains (or globules) were found, ranging in size from 2.0 to 76.0 kb. Interestingly, transcription factor-bound motifs were enriched in the globule boundaries. All the cellulase and xylanase genes were discretely distributed in the 3D model of the P. oxalicum genome as a result of few cis- and trans-interactions. Our results from this study provide a global view of chromatin interactions in the P. oxalicum genome and will act as a resource for studying spatial regulation of gene expression in filamentous fungi. IMPORTANCE The spatial structure of chromatin plays important roles in normal cell functions and the regulation of gene expression. The three-dimensional (3D) architectures of the genomes of many mammals and plants have been elucidated, but corresponding studies on filamentous fungi, which play vital roles as decomposers of organic matter in the soil, are very limited. Penicillium oxalicum is one of the predominant cellulolytic aerobic fungi in subtropical and tropical forest soils and can secrete integrative cellulase and xylanase under integrated regulatory control, degrading plant biomass highly efficiently. In the present study, we employed Hi-C technology to construct the 3D genome model of P. oxalicum strain HP7-1 and to further investigate cellulase and xylanase as well as transcription factor genes in 3D genome. These results provide a resource to achieve a deeper understanding of cell function and the regulation of gene expression in filamentous fungi.

Defining the role of the polyasparagine repeat domain of the S. cerevisiae transcription factor Azf1p

Across eukaryotes, homopolymeric repeats of amino acids are enriched in regulatory proteins such as transcription factors and chromatin remodelers. These domains play important roles in signaling, binding, prion formation, and functional phase separation. Azf1p is a prion-forming yeast transcription factor that contains two homorepeat domains, a polyglutamine and a polyasparagine domain. In this work, we report a new phenotype for Azf1p and identify a large set of genes that are regulated by Azf1p during growth in glucose. We show that the polyasparagine (polyN) domain plays a subtle role in transcription but is dispensable for Azf1p localization and prion formation. Genes upregulated upon deletion of the polyN domain are enriched in functions related to carbon metabolism and storage. This domain may therefore be a useful target for engineering yeast strains for fermentation applications and small molecule production. We also report that both the polyasparagine and polyglutamine domains vary in length across strains of S. cerevisiae and propose a model for how this variation may impact protein function.

DIA-based systems biology approach unveils E3 ubiquitin ligase-dependent responses to a metabolic shift

The yeast Saccharomyces cerevisiae is a powerful model system for systems-wide biology screens and large-scale proteomics methods. Nearly complete proteomics coverage has been achieved owing to advances in mass spectrometry. However, it remains challenging to scale this technology for rapid and high-throughput analysis of the yeast proteome to investigate biological pathways on a global scale. Here we describe a systems biology workflow employing plate-based sample preparation and rapid, single-run, data-independent mass spectrometry analysis (DIA). Our approach is straightforward, easy to implement, and enables quantitative profiling and comparisons of hundreds of nearly complete yeast proteomes in only a few days. We evaluate its capability by characterizing changes in the yeast proteome in response to environmental perturbations, identifying distinct responses to each of them and providing a comprehensive resource of these responses. Apart from rapidly recapitulating previously observed responses, we characterized carbon source-dependent regulation of the GID E3 ligase, an important regulator of cellular metabolism during the switch between gluconeogenic and glycolytic growth conditions. This unveiled regulatory targets of the GID ligase during a metabolic switch. Our comprehensive yeast system readout pinpointed effects of a single deletion or point mutation in the GID complex on the global proteome, allowing the identification and validation of targets of the GID E3 ligase. Moreover, this approach allowed the identification of targets from multiple cellular pathways that display distinct patterns of regulation. Although developed in yeast, rapid whole-proteome-based readouts can serve as comprehensive systems-level assays in all cellular systems.

The secondary metabolite regulator, BbSmr1, is a central regulator of conidiation via the BrlA-AbaA-WetA pathway in Beauveria bassiana

The filamentous fungus Beauveria bassiana, an insect fungal pathogen, is widely used for pest biocontrol. Aerial conidia are infectious propagules, and their yield and viability greatly affect the field application of this fungus; however, little is known about the molecular regulatory mechanism of the triggered conidiation. In the present study, we find that the secondary metabolite regulator BbSmr1 is involved in the regulation of asexual conidiation development and stress response in B. bassiana. A deficiency in Bbsmr1 results in a prominent fluffy-like phenotype on solid medium, decreased conidial yield, accelerated conidial germination, as well as increased tolerance to H₂ O₂ stress and cell wall inhibitors. The deletion of Bbsmr1 also leads to thickened conidial cell walls and changed cell epitopes. Overexpressing either BbbrlA or BbabaA in the ∆Bbsmr1 strain can rescue the phenotypes of conidial development and stress response. BbSmr1 activates BbbrlA transcription by directly binding to the A4GA3 sequence of the BbbrlA promoter. BbBrlA in turn binds to the promoter of Bbsmr1 and negatively regulates the expression of Bbsmr1. These results indicate that BbSmr1 positively regulates conidial development in B. bassiana by activating the central development pathway BrlA-AbaA-WetA and provides insights into the developmental regulatory mechanism of entomopathogenic fungi.

Investigating the Influence of Glycerol on the Utilization of Glucose in Yarrowia lipolytica Using RNA-Seq-Based Transcriptomics

Glycerol is considered as a promising substrate for biotechnological applications and the non-conventional yeast Yarrowia lipolytica has been used extensively for the valorization of this compound. Contrary to S. cerevisiae, Y. lipolytica seems to prefer glycerol over glucose and it has been reported previously that the presence of glycerol can suppress the consumption of glucose in co-substrate fermentations. Based on these observations, we hypothesized glycerol repression-like effects in Y. lipolytica, which are converse to well described carbon repression mechanisms ensuring the prioritized use of glucose (e.g., in S. cerevisiae). We therefore aimed to investigate this effect on the level of transcription. Strains varying in the degree of glucose suppression were chosen and characterized in high-resolution growth screenings, resulting in the detection of different growth phenotypes under glycerol-glucose mixed conditions. Two strains, IBT and W29, were selected and cultivated in chemostats using glucose, glycerol and glucose/glycerol as carbon sources, followed by an RNA-Seq-based transcriptome analysis. We could show that several transporters were significantly higher expressed in W29, which is potentially related to the observed physiological differences. However, most of the expression variation between the strains were regardless of the carbon source applied, and cross-comparisons revealed that the strain-specific carbon source responses underwent in the opposite direction. A deeper analysis of the substrate specific carbon source response led to the identification of several differentially expressed genes with orthologous functions related to signal transduction and transcriptional regulation. This study provides an initial investigation on potentially novel carbon source regulation mechanisms in yeasts.

A Novel Cys2His2 Zinc Finger Homolog of AZF1 Modulates Holocellulase Expression in Trichoderma reesei

In this work, we used a systems biology approach to map new regulatory interactions in Trichoderma reesei controlling the expression of genes encoding cellulase and hemicellulase. By integrating transcriptomics related to complex biomass degradation, we were able to identify a novel transcriptional regulator which is able to activate the expression of these genes in response to two different cellulose sources. In vivo experimental validation confirmed the role of this new regulator in several other processes related to carbon source utilization and nutrient transport. Therefore, this work revealed novel forms of regulatory interaction in this model system for plant biomass deconstruction and also represented a new approach that could be easy applied to other organisms.

Filamentous fungi are remarkable producers of enzymes dedicated to the degradation of sugar polymers found in the plant cell wall. Here, we integrated transcriptomic data to identify novel transcription factors (TFs) related to the control of gene expression of lignocellulosic hydrolases in Trichoderma reesei and Aspergillus nidulans. Using various sets of differentially expressed genes, we identified some putative cis-regulatory elements that were related to known binding sites for Saccharomyces cerevisiae TFs. Comparative genomics allowed the identification of six transcriptional factors in filamentous fungi that have corresponding S. cerevisiae homologs. Additionally, a knockout strain of T. reesei lacking one of these TFs (S. cerevisiaeAZF1 homolog) displayed strong reductions in the levels of expression of several cellulase-encoding genes in response to both Avicel and sugarcane bagasse, revealing a new player in the complex regulatory network operating in filamentous fungi during plant biomass degradation. Finally, RNA sequencing (RNA-seq) analysis showed the scope of the AZF1 homologue in regulating a number of processes in T. reesei, and chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) provided evidence for the direct interaction of this TF in the promoter regions of cel7a, cel45a, and swo. Therefore, we identified here a novel TF which plays a positive effect in the expression of cellulase-encoding genes in T. reesei.

IMPORTANCE In this work, we used a systems biology approach to map new regulatory interactions in Trichoderma reesei controlling the expression of genes encoding cellulase and hemicellulase. By integrating transcriptomics related to complex biomass degradation, we were able to identify a novel transcriptional regulator which is able to activate the expression of these genes in response to two different cellulose sources. In vivo experimental validation confirmed the role of this new regulator in several other processes related to carbon source utilization and nutrient transport. Therefore, this work revealed novel forms of regulatory interaction in this model system for plant biomass deconstruction and also represented a new approach that could be easy applied to other organisms.

Rewired cellular signaling coordinates sugar and hypoxic responses for anaerobic xylose fermentation in yeast

Microbes can be metabolically engineered to produce biofuels and biochemicals, but rerouting metabolic flux toward products is a major hurdle without a systems-level understanding of how cellular flux is controlled. To understand flux rerouting, we investigated a panel of Saccharomyces cerevisiae strains with progressive improvements in anaerobic fermentation of xylose, a sugar abundant in sustainable plant biomass used for biofuel production. We combined comparative transcriptomics, proteomics, and phosphoproteomics with network analysis to understand the physiology of improved anaerobic xylose fermentation. Our results show that upstream regulatory changes produce a suite of physiological effects that collectively impact the phenotype. Evolved strains show an unusual co-activation of Protein Kinase A (PKA) and Snf1, thus combining responses seen during feast on glucose and famine on non-preferred sugars. Surprisingly, these regulatory changes were required to mount the hypoxic response when cells were grown on xylose, revealing a previously unknown connection between sugar source and anaerobic response. Network analysis identified several downstream transcription factors that play a significant, but on their own minor, role in anaerobic xylose fermentation, consistent with the combinatorial effects of small-impact changes. We also discovered that different routes of PKA activation produce distinct phenotypes: deletion of the RAS/PKA inhibitor IRA2 promotes xylose growth and metabolism, whereas deletion of PKA inhibitor BCY1 decouples growth from metabolism to enable robust fermentation without division. Comparing phosphoproteomic changes across ira2Δ and bcy1Δ strains implicated regulatory changes linked to xylose-dependent growth versus metabolism. Together, our results present a picture of the metabolic logic behind anaerobic xylose flux and suggest that widespread cellular remodeling, rather than individual metabolic changes, is an important goal for metabolic engineering.

Genomic and transcriptomic comparison of Aspergillus oryzae strains: a case study in soy sauce koji fermentation

The filamentous fungus Aspergillus oryzae is used in soy sauce koji making due to its high productivity of hydrolytic enzymes. In this study, we compared the genomes and transcriptomes of an industrial strain RD2 and a strain with decreased fermentation performance TS2, aiming to explain their phenotypic differences at the molecular level. Under the regulation of conidiation and fermentation conditions, the enhanced hydrolytic enzyme production and flavor precursor formation in RD2 described a complete expression profile necessary to maintain desirable fermentation performance. By contrast, central carbon metabolism was up-regulated in TS2 for fast growth, suggesting a conflicting relationship between mycelium growth and fermentation performance. Accumulation of mutations also lowered the fermentation performance of TS2. Our study has deepened the understanding of the metabolism and related regulatory mechanisms in desirable koji fermentation. A list of potential molecular markers identified here could facilitate targeted strain maintenance and improvement for better koji fermentation.

Understanding the tolerance of the industrial yeast Saccharomyces cerevisiae against a major class of toxic aldehyde compounds

Development of the next-generation biocatalyst is vital for fermentation-based industrial applications and a sustainable bio-based economy. Overcoming the major class of toxic compounds associated with lignocellulose-to-biofuels conversion is one of the significant challenges for new strain development. A significant number of investigations have been made to understand mechanisms of the tolerance for industrial yeast. It is humbling to learn how complicated the cell's response to the toxic chemicals is and how little we have known about yeast tolerance in the universe of the living cell. This study updates our current knowledge on the tolerance of industrial yeast against aldehyde inhibitory compounds at cellular, molecular and the genomic levels. It is comprehensive yet specific based on reproducible evidence and cross confirmed findings from different investigations using varied experimental approaches. This research approaches a rational foundation toward a more comprehensive understanding on the yeast tolerance. Discussions and perspectives are also proposed for continued exploring the puzzle of the yeast tolerance to aid the next-generation biocatalyst development.

Genome sequencing and transcriptome analysis of Trichoderma reesei QM9978 strain reveals a distal chromosome translocation to be responsible for loss of vib1 expression and loss of cellulase induction

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.

A Multistate Toggle Switch Defines Fungal Cell Fates and Is Regulated by Synergistic Genetic Cues

Heritable epigenetic changes underlie the ability of cells to differentiate into distinct cell types. Here, we demonstrate that the fungal pathogen Candida tropicalis exhibits multipotency, undergoing stochastic and reversible switching between three cellular states. The three cell states exhibit unique cellular morphologies, growth rates, and global gene expression profiles. Genetic analysis identified six transcription factors that play key roles in regulating cell differentiation. In particular, we show that forced expression of Wor1 or Efg1 transcription factors can be used to manipulate transitions between all three cell states. A model for tristability is proposed in which Wor1 and Efg1 are self-activating but mutually antagonistic transcription factors, thereby forming a symmetrical self-activating toggle switch. We explicitly test this model and show that ectopic expression of WOR1 can induce white-to-hybrid-to-opaque switching, whereas ectopic expression of EFG1 drives switching in the opposite direction, from opaque-to-hybrid-to-white cell states. We also address the stability of induced cell states and demonstrate that stable differentiation events require ectopic gene expression in combination with chromatin-based cues. These studies therefore experimentally test a model of multistate stability and demonstrate that transcriptional circuits act synergistically with chromatin-based changes to drive cell state transitions. We also establish close mechanistic parallels between phenotypic switching in unicellular fungi and cell fate decisions during stem cell reprogramming.

Light governs asexual differentiation in the grey mould fungus Botrytis cinerea via the putative transcription factor BcLTF2

Botrytis cinerea is a plant pathogenic fungus known for its utilization of light as environmental cue to regulate asexual differentiation: conidia are formed in the light, while sclerotia are formed in the dark. As no orthologues of known regulators of conidiation (e.g., Aspergillus nidulans BrlA, Neurospora crassa FL) exist in the Leotiomycetes, we initiated a de novo approach to identify the functional counterpart in B. cinerea. The search revealed the light-responsive C2H2 transcription factor BcLTF2 whose expression - usually restricted to light conditions - is necessary and sufficient to induce conidiation and simultaneously to suppress sclerotial development. Light-induced expression of bcltf2 is mediated via a so far unknown pathway, and is attenuated in a (blue) light-dependent fashion by the White Collar complex, BcLTF1 and the VELVET complex. Mutation of either component leads to increased bcltf2 expression and causes light-independent conidiation (always conidia phenotype). Hence, the tight regulation of bcltf2 governs the balance between vegetative growth that allows for the colonization of the substrate and subsequent reproduction via conidia in the light. The orthologue ssltf2 in the closely related species Sclerotinia sclerotiorum is not significantly expressed suggesting that its deregulation may cause the lack of the conidiation program in this fungus.

Disruption of a cystine transporter downregulates expression of genes involved in sulfur regulation and cellular respiration

Cystine and cysteine are important molecules for pathways such as redox signaling and regulation, and thus identifying cellular deficits upon deletion of the Saccharomyces cerevisiae cystine transporter Ers1p allows for a further understanding of cystine homeostasis. Previous complementation studies using the human ortholog suggest yeast Ers1p is a cystine transporter. Human CTNS encodes the protein Cystinosin, a cystine transporter that is embedded in the lysosomal membrane and facilitates the export of cystine from the lysosome. When CTNS is mutated, cystine transport is disrupted, leading to cystine accumulation, the diagnostic hallmark of the lysosomal storage disorder cystinosis. Here, we provide biochemical evidence for Ers1p-dependent cystine transport. However, the accumulation of intracellular cystine is not observed when the ERS1 gene is deleted from ers1-Δ yeast, supporting the existence of modifier genes that provide a mechanism in ers1-Δ yeast that prevents or corrects cystine accumulation. Upon comparison of the transcriptomes of isogenic ERS1+ and ers1-Δ strains of S. cerevisiae by DNA microarray followed by targeted qPCR, sixteen genes were identified as being differentially expressed between the two genotypes. Genes that encode proteins functioning in sulfur regulation, cellular respiration, and general transport were enriched in our screen, demonstrating pleiotropic effects of ers1-Δ. These results give insight into yeast cystine regulation and the multiple, seemingly distal, pathways that involve proper cystine recycling.

Summary: We identify genes that are differentially expressed in yeast lacking vacuolar cystine transporter Ers1p in order to find pathways, such as respiration and sulfur regulation, that are associated with cystine homeostasis.

The Trichoderma atroviride putative transcription factor Blu7 controls light responsiveness and tolerance

Background

Most living organisms use sunlight as a source of energy and/or information about their environment. Consequently, they have developed mechanisms to sense light quality and quantity. In the fungus Trichoderma atroviride blue-light is perceived through the Blue Light Regulator Complex, which in turn up-regulates a set of genes (blu) and down-regulates another set (bld), triggering asexual reproduction. To gain insight into this process, we characterized the blu7 gene, which encodes a protein containing a C2H2 zinc finger domain.

Results

Δblu7 mutants show reduced conidiation at low light fluences, which is still clear even when exposed to saturating light. For the first time we show a genome wide survey of light regulated gene expression in T. atroviride, including RNA-seq analyses of the wild type and the Δblu7 strains after brief exposure to blue-light. Our data show a reduction in the number of induced genes and an increase in down-regulated genes in the mutant. Light activates stress responses and several metabolic processes in the wild type strain that are no longer activated in the mutant. In agreement with the misregulation of metabolic processes, continuous exposure to white light strongly inhibited growth of the ∆blu7 mutant, in a carbon source dependent fashion. RNA-seq analyses under constant white light using glucose as sole carbon source revealed that localization and transport process present the opposite regulation pattern in the ∆blu7 and wild type strains. Genes related to amino acid, sugar and general transporters were enriched in the induced genes in the mutant and the repressed genes of the wild type. Peptone supplemented in the media restored growth of the ∆blu7 mutant in constant light, suggesting a role of Blu7 in the regulation of nitrogen metabolism in the presence of light.

Conclusions

Blu7 appears to regulate light sensitivity in terms of induction of conidiation, and to play a major role in supporting growth under continuous exposure to light. The diminished conidiation observed in ∆blu7 mutants is likely due to misregulation of the cAMP signaling pathway and ROS production, whereas their low tolerance to continuous exposure to light indicates that Blu7 is required for adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2639-9) contains supplementary material, which is available to authorized users.

Genomic and transcriptome analyses reveal that MAPK- and phosphatidylinositol-signaling pathways mediate tolerance to 5-hydroxymethyl-2-furaldehyde for industrial yeast Saccharomyces cerevisiae

The industrial yeast Saccharomyces cerevisiae is a traditional ethanologenic agent and a promising biocatalyst for advanced biofuels production using lignocellulose materials. Here we present the genomic background of type strain NRRL Y-12632 and its transcriptomic response to 5-hydroxymethyl-2-furaldehyde (HMF), a commonly encountered toxic compound liberated from lignocellulosic-biomass pretreatment, in dissecting the genomic mechanisms of yeast tolerance. Compared with the genome of laboratory model strain S288C, we identified more than 32,000 SNPs in Y-12632 with 23,000 missense and nonsense SNPs. Enriched sequence mutations occurred for genes involved in MAPK- and phosphatidylinositol (PI)- signaling pathways in strain Y-12632, with 41 and 13 genes containing non-synonymous SNPs, respectively. Many of these mutated genes displayed consistent up-regulated signature expressions in response to challenges of 30 mM HMF. Analogous single-gene deletion mutations of these genes showed significantly sensitive growth response on a synthetic medium containing 20 mM HMF. Our results suggest at least three MAPK-signaling pathways, especially for the cell-wall integrity pathway, and PI-signaling pathways to be involved in mediation of yeast tolerance against HMF in industrial yeast Saccharomyces cerevisiae. Higher levels of sequence variations were also observed for genes involved in purine and pyrimidine metabolism pathways.

Comprehensive analysis of forty yeast microarray datasets reveals a novel subset of genes (APha-RiB) consistently negatively associated with ribosome biogenesis

BACKGROUND

The scale and complexity of genomic data lend themselves to analysis using sophisticated mathematical techniques to yield information that can generate new hypotheses and so guide further experimental investigations. An ensemble clustering method has the ability to perform consensus clustering over the same set of genes from different microarray datasets by combining results from different clustering methods into a single consensus result.

RESULTS

In this paper we have performed comprehensive analysis of forty yeast microarray datasets. One recently described Bi-CoPaM method can analyse expressions of the same set of genes from various microarray datasets while using different clustering methods, and then combine these results into a single consensus result whose clusters' tightness is tunable from tight, specific clusters to wide, overlapping clusters. This has been adopted in a novel way over genome-wide data from forty yeast microarray datasets to discover two clusters of genes that are consistently co-expressed over all of these datasets from different biological contexts and various experimental conditions. Most strikingly, average expression profiles of those clusters are consistently negatively correlated in all of the forty datasets while neither profile leads or lags the other.

CONCLUSIONS

The first cluster is enriched with ribosomal biogenesis genes. The biological processes of most of the genes in the second cluster are either unknown or apparently unrelated although they show high connectivity in protein-protein and genetic interaction networks. Therefore, it is possible that this mostly uncharacterised cluster and the ribosomal biogenesis cluster are transcriptionally oppositely regulated by some common machinery. Moreover, we anticipate that the genes included in this previously unknown cluster participate in generic, in contrast to specific, stress response processes. These novel findings illuminate coordinated gene expression in yeast and suggest several hypotheses for future experimental functional work. Additionally, we have demonstrated the usefulness of the Bi-CoPaM-based approach, which may be helpful for the analysis of other groups of (microarray) datasets from other species and systems for the exploration of global genetic co-expression.

Characterization of global gene expression during assurance of lifespan extension by caloric restriction in budding yeast

Caloric restriction (CR) is the best-studied intervention known to delay aging and extend lifespan in evolutionarily distant organisms ranging from yeast to mammals in the laboratory. Although the effect of CR on lifespan extension has been investigated for nearly 80years, the molecular mechanisms of CR are still elusive. Consequently, it is important to understand the fundamental mechanisms of when and how lifespan is affected by CR. In this study, we first identified the time-windows during which CR assured cellular longevity by switching cells from culture media containing 2% or 0.5% glucose to water, which allows us to observe CR and non-calorically-restricted cells under the same conditions. We also constructed time-dependent gene expression profiles and selected 646 genes that showed significant changes and correlations with the lifespan-extending effect of CR. The positively correlated genes participated in transcriptional regulation, ribosomal RNA processing and nuclear genome stability, while the negatively correlated genes were involved in the regulation of several metabolic pathways, endoplasmic reticulum function, stress response and cell cycle progression. Furthermore, we discovered major upstream regulators of those significantly changed genes, including AZF1 (YOR113W), HSF1 (YGL073W) and XBP1 (YIL101C). Deletions of two genes, AZF1 and XBP1 (HSF1 is essential and was thus not tested), were confirmed to lessen the lifespan extension mediated by CR. The absence of these genes in the tor1Δ and ras2Δ backgrounds did show non-overlapping effects with regard to CLS, suggesting differences between the CR mechanism for Tor and Ras signaling.

Systematic measurement of transcription factor-DNA interactions by targeted mass spectrometry identifies candidate gene regulatory proteins

Regulation of gene expression involves the orchestrated interaction of a large number of proteins with transcriptional regulatory elements in the context of chromatin. Our understanding of gene regulation is limited by the lack of a protein measurement technology that can systematically detect and quantify the ensemble of proteins associated with the transcriptional regulatory elements of specific genes. Here, we introduce a set of selected reaction monitoring (SRM) assays for the systematic measurement of 464 proteins with known or suspected roles in transcriptional regulation at RNA polymerase II transcribed promoters in Saccharomyces cerevisiae. Measurement of these proteins in nuclear extracts by SRM permitted the reproducible quantification of 42% of the proteins over a wide range of abundances. By deploying the assay to systematically identify DNA binding transcriptional regulators that interact with the environmentally regulated FLO11 promoter in cell extracts, we identified 15 regulators that bound specifically to distinct regions along ∼600 bp of the regulatory sequence. Importantly, the dataset includes a number of regulators that have been shown to either control FLO11 expression or localize to these regulatory regions in vivo. We further validated the utility of the approach by demonstrating that two of the SRM-identified factors, Mot3 and Azf1, are required for proper FLO11 expression. These results demonstrate the utility of SRM-based targeted proteomics to guide the identification of gene-specific transcriptional regulators.

Functional Implications of Local DNA Structures in Regulatory Motifs

The three-dimensional structure of DNA has been proposed to be a major determinant for functional transcription factors (TFs) and DNA interaction. Here, we use hydroxyl radical cleavage pattern as a measure of local DNA structure. We compared the conservation between DNA sequence and structure in terms of information content and attempted to assess the functional implications of DNA structures in regulatory motifs. We used statistical methods to evaluate the structural divergence of substituting a single position within a binding site and applied them to a collection of putative regulatory motifs. The following are our major observations: (i) we observed more information in structural alignment than in the corresponding sequence alignment for most of the transcriptional factors; (ii) for each TF, majority of positions have more information in the structural alignment as compared to the sequence alignment; (iii) we further defined a DNA structural divergence score (SD score) for each wild-type and mutant pair that is distinguished by single-base mutation. The SD score for benign mutations is significantly lower than that of switch mutations. This indicates structural conservation is also important for TFBS to be functional and DNA structures will provide previously unappreciated information for TF to realize the binding specificity.

Reconstructing differentially co-expressed gene modules and regulatory networks of soybean cells

BACKGROUND

Current experimental evidence indicates that functionally related genes show coordinated expression in order to perform their cellular functions. In this way, the cell transcriptional machinery can respond optimally to internal or external stimuli. This provides a research opportunity to identify and study co-expressed gene modules whose transcription is controlled by shared gene regulatory networks.

RESULTS

We developed and integrated a set of computational methods of differential gene expression analysis, gene clustering, gene network inference, gene function prediction, and DNA motif identification to automatically identify differentially co-expressed gene modules, reconstruct their regulatory networks, and validate their correctness. We tested the methods using microarray data derived from soybean cells grown under various stress conditions. Our methods were able to identify 42 coherent gene modules within which average gene expression correlation coefficients are greater than 0.8 and reconstruct their putative regulatory networks. A total of 32 modules and their regulatory networks were further validated by the coherence of predicted gene functions and the consistency of putative transcription factor binding motifs. Approximately half of the 32 modules were partially supported by the literature, which demonstrates that the bioinformatic methods used can help elucidate the molecular responses of soybean cells upon various environmental stresses.

CONCLUSIONS

The bioinformatics methods and genome-wide data sources for gene expression, clustering, regulation, and function analysis were integrated seamlessly into one modular protocol to systematically analyze and infer modules and networks from only differential expression genes in soybean cells grown under stress conditions. Our approach appears to effectively reduce the complexity of the problem, and is sufficiently robust and accurate to generate a rather complete and detailed view of putative soybean gene transcription logic potentially underlying the responses to the various environmental challenges. The same automated method can also be applied to reconstruct differentially co-expressed gene modules and their regulatory networks from gene expression data of any other transcriptome.

Conserved and specific responses to hypoxia in Aspergillus oryzae and Aspergillus nidulans determined by comparative transcriptomics

Hypoxia imposes stress on filamentous fungi that require oxygen to proliferate. Global transcription analysis of Aspergillus oryzae grown under hypoxic conditions found that the expression of about 50% of 4,244 affected genes was either induced or repressed more than 2-fold. A comparison of these genes with the hypoxically regulated genes of Aspergillus nidulans based on their predicted amino acid sequences classified them as bi-directional best hit (BBH), one-way best hit (extra homolog, EH), and no-hit (non-syntenic genes, NSG) genes. Clustering analysis of the BBH genes indicated that A. oryzae and A. nidulans down-regulated global translation and transcription under hypoxic conditions, respectively. Under hypoxic conditions, both fungi up-regulated genes for alcohol fermentation and the γ-aminobutyrate shunt of the tricarboxylate cycle, whereas A. oryzae up-regulated the glyoxylate pathway, indicating that both fungi eliminate NADH accumulation under hypoxic conditions. The A. oryzae NS genes included specific genes for secondary and nitric oxide metabolism under hypoxic conditions. This comparative transcriptomic analysis discovered common and strain-specific responses to hypoxia in hypoxic Aspergillus species.

Epigenetic transgenerational actions of vinclozolin on promoter regions of the sperm epigenome

Previous observations have demonstrated that embryonic exposure to the endocrine disruptor vinclozolin during gonadal sex determination promotes transgenerational adult onset disease such as male infertility, kidney disease, prostate disease, immune abnormalities and tumor development. The current study investigates genome-wide promoter DNA methylation alterations in the sperm of F3 generation rats whose F0 generation mother was exposed to vinclozolin. A methylated DNA immunoprecipitation with methyl-cytosine antibody followed by a promoter tilling microarray (MeDIP-Chip) procedure was used to identify 52 different regions with statistically significant altered methylation in the sperm promoter epigenome. Mass spectrometry bisulfite analysis was used to map the CpG DNA methylation and 16 differential DNA methylation regions were confirmed, while the remainder could not be analyzed due to bisulfite technical limitations. Analysis of these validated regions identified a consensus DNA sequence (motif) that associated with 75% of the promoters. Interestingly, only 16.8% of a random set of 125 promoters contained this motif. One candidate promoter (Fam111a) was found to be due to a copy number variation (CNV) and not a methylation change, suggesting initial alterations in the germline epigenome may promote genetic abnormalities such as induced CNV in later generations. This study identifies differential DNA methylation sites in promoter regions three generations after the initial exposure and identifies common genome features present in these regions. In addition to primary epimutations, a potential indirect genetic abnormality was identified, and both are postulated to be involved in the epigenetic transgenerational inheritance observed. This study confirms that an environmental agent has the ability to induce epigenetic transgenerational changes in the sperm epigenome.

Comparative analyses of time-course gene expression profiles of the long-lived sch9Delta mutant

In an attempt to elucidate the underlying longevity-promoting mechanisms of mutants lacking SCH9, which live three times as long as wild type chronologically, we measured their time-course gene expression profiles. We interpreted their expression time differences by statistical inferences based on prior biological knowledge, and identified the following significant changes: (i) between 12 and 24 h, stress response genes were up-regulated by larger fold changes and ribosomal RNA (rRNA) processing genes were down-regulated more dramatically; (ii) mitochondrial ribosomal protein genes were not up-regulated between 12 and 60 h as wild type were; (iii) electron transport, oxidative phosphorylation and TCA genes were down-regulated early; (iv) the up-regulation of TCA and electron transport was accompanied by deep down-regulation of rRNA processing over time; and (v) rRNA processing genes were more volatile over time, and three associated cis-regulatory elements [rRNA processing element (rRPE), polymerase A and C (PAC) and glucose response element (GRE)] were identified. Deletion of AZF1, which encodes the transcriptional factor that binds to the GRE element, reversed the lifespan extension of sch9Δ. The significant alterations in these time-dependent expression profiles imply that the lack of SCH9 turns on the longevity programme that extends the lifespan through changes in metabolic pathways and protection mechanisms, particularly, the regulation of aerobic respiration and rRNA processing.

High-resolution analysis of condition-specific regulatory modules in Saccharomyces cerevisiae

A novel approach for identifying condition-specific regulatory modules in yeast reveals functionally distinct coregulated submodules.

We present an approach for identifying condition-specific regulatory modules by using separate units of gene expression profiles along with ChIP-chip and motif data from Saccharomyces cerevisiae. By investigating the unique and common features of the obtained condition-specific modules, we detected several important properties of transcriptional network reorganization. Our approach reveals the functionally distinct coregulated submodules embedded in a coexpressed gene module and provides an effective method for identifying various condition-specific regulatory events at high resolution.

Protein kinase A, TOR, and glucose transport control the response to nutrient repletion in Saccharomyces cerevisiae

Nutrient repletion leads to substantial restructuring of the transcriptome in Saccharomyces cerevisiae. The expression levels of approximately one-third of all S. cerevisiae genes are altered at least twofold when a nutrient-depleted culture is transferred to fresh medium. Several nutrient-sensing pathways are known to play a role in this process, but the relative contribution that each pathway makes to the total response has not been determined. To better understand this, we used a chemical-genetic approach to block the protein kinase A (PKA), TOR (target of rapamycin), and glucose transport pathways, alone and in combination. Of the three pathways, we found that loss of PKA produced the largest effect on the transcriptional response; however, many genes required both PKA and TOR for proper nutrient regulation. Those genes that did not require PKA or TOR for nutrient regulation were dependent on glucose transport for either nutrient induction or repression. Therefore, loss of these three pathways is sufficient to prevent virtually the entire transcriptional response to fresh medium. In the absence of fresh medium, activation of the cyclic AMP/PKA pathway does not induce cellular growth; nevertheless, PKA activation induced a substantial fraction of the PKA-dependent genes. In contrast, the absence of fresh medium strongly limited gene repression by PKA. These results account for the signals needed to generate the transcriptional responses to glucose, including induction of growth genes required for protein synthesis and repression of stress genes, as well as the classical glucose repression and hexose transporter responses.

Defining Candida albicans stationary phase by cellular and DNA replication, gene expression and regulation

Stationary phase Candida albicans yeast cells harbour properties of better adherence, virulence and elevated drug resistance. C. albicans stationary phase is not well characterized in vitro either physiologically or molecularly. C. albicans yeast cells were grown in rich medium with 2% glucose. Based on growth and DNA profiles of cells, and by measurement of glucose and ethanol in the medium, we defined the timing of C. albicans entry into different growth transitions. We found that, compared with 24 h cells, mRNA content was less abundant in post-diauxic shift phase and even less in stationary phase C. albicans cells. Further analysis of the C. albicans transcriptome with oligonucleotide-based microarrays revealed that although the overall mRNA content had decreased, transcripts of many genes increased in post-diauxic shift phase as well as stationary phase. Genes involved in processes such as gluconeogenesis, stress resistance, adherence, DNA repair and ageing were expressed at higher levels at and beyond post-diauxic shift phase. Many C. albicans genes associated with virulence, drug resistance and cell-wall biosynthesis were expressed only at stationary phase. By screening 108 C. albicans transcription factor and cell-wall mutants we identified 17 genes essential for either entry or survival in stationary phase at 30 degrees C.