Transcriptional Rewiring: The Proof Is in the Eating

Transcriptional activation and coactivator binding by yeast Ino2 and human proto-oncoprotein c-Myc

Basic helix-loop-helix domains in yeast regulatory proteins Ino2 and Ino4 mediate formation of a heterodimer which binds to and activates expression of phospholipid biosynthetic genes. The human proto-oncoprotein c-Myc (Myc) and its binding partner Max activate genes important for cellular proliferation and contain functional domains structure and position of which strongly resembles Ino2 and Ino4. Since Ino2-Myc and Ino4-Max may be considered as orthologs we performed functional comparisons in yeast. We demonstrate that Myc and Max could be stably synthesized in S. cerevisiae and together significantly activated a target gene of Ino2/Ino4 but nevertheless were unable to functionally complement an ino2 ino4 double mutant. We also map two efficient transcriptional activation domains in the N-terminus of Myc (TAD1: aa 1-41 and TAD2: aa 91-140), corresponding to TAD positions in Ino2. We finally show that coactivators such as TFIID subunits Taf1, Taf4, Taf6, Taf10 and Taf12 as well as ATPase subunits of chromatin remodelling complexes Swi2, Sth1 and Ino80 previously shown to interact with TADs of Ino2 were also able to bind TADs of Myc, supporting the view that heterodimers Ino2/Ino4 and Myc/Max are evolutionary related but have undergone transcriptional rewiring of target genes.

Natural trait variation across Saccharomycotina species

Among molecular biologists, the group of fungi called Saccharomycotina is famous for its yeasts. These yeasts in turn are famous for what they have in common-genetic, biochemical, and cell-biological characteristics that serve as models for plants and animals. But behind the apparent homogeneity of Saccharomycotina species lie a wealth of differences. In this review, we discuss traits that vary across the Saccharomycotina subphylum. We describe cases of bright pigmentation; a zoo of cell shapes; metabolic specialties; and species with unique rules of gene regulation. We discuss the genetics of this diversity and why it matters, including insights into basic evolutionary principles with relevance across Eukarya.

The birth, evolution and death of metabolic gene clusters in fungi

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

Regulators of plant biomass degradation in ascomycetous fungi

Fungi play a major role in the global carbon cycle because of their ability to utilize plant biomass (polysaccharides, proteins, and lignin) as carbon source. Due to the complexity and heterogenic composition of plant biomass, fungi need to produce a broad range of degrading enzymes, matching the composition of (part of) the prevalent substrate. This process is dependent on a network of regulators that not only control the extracellular enzymes that degrade the biomass, but also the metabolic pathways needed to metabolize the resulting monomers. This review will summarize the current knowledge on regulation of plant biomass utilization in fungi and compare the differences between fungal species, focusing in particular on the presence or absence of the regulators involved in this process.

Ubiquitin-Dependent Modification of Skeletal Muscle by the Parasitic Nematode, Trichinella spiralis

Trichinella spiralis is a muscle-specific parasitic worm that is uniquely intracellular. T. spiralis reprograms terminally differentiated skeletal muscle cells causing them to de-differentiate and re-enter the cell cycle, a process that cannot occur naturally in mammalian skeletal muscle cells, but one that holds great therapeutic potential. Although the host ubiquitin pathway is a common target for viruses and bacteria during infection, its role in parasite pathogenesis has been largely overlooked. Here we demonstrate that the secreted proteins of T. spiralis contain E2 Ub-conjugating and E3 Ub-ligase activity. The E2 activity is attributed to TsUBE2L3, a novel and conserved T. spiralis enzyme located in the secretory organ of the parasite during the muscle stages of infection. TsUBE2L3 cannot function with any T.spiralis secreted E3, but specifically binds to a panel of human RING E3 ligases, including the RBR E3 ARIH2 with which it interacts with a higher affinity than the mammalian ortholog UbcH7/UBE2L3. Expression of TsUBE2L3 in skeletal muscle cells causes a global downregulation in protein ubiquitination, most predominantly affecting motor, sarcomeric and extracellular matrix proteins, thus mediating their stabilization with regards to proteasomal degradation. This effect is not observed in the presence of the mammalian ortholog, suggesting functional divergence in the evolution of the parasite protein. These findings demonstrate the first example of host-parasite interactions via a parasite-derived Ub conjugating enzyme; an E2 that demonstrates a novel muscle protein stabilization function.

Parasitic worms often establish long-lasting infections in their hosts; tightly regulating their surroundings to strike a delicate balance between host cell modulation and protection that will ensure their replication. This is accomplished via the active secretion of parasite glycolipids and glycoproteins into the host. Trichinella spiralis, a parasitic nematode that infects skeletal muscle of mammals, birds and reptiles, is the only parasitic worm with a true intracellular stage. T. spiralis invade terminally differentiated myotubes, reprogramming them to de-differentiate and re-enter the cell cycle, a process that cannot occur naturally in mammalian skeletal muscle cells, but one that holds great therapeutic potential. We have identified and characterized a novel T. spiralis secreted protein that, despite a high level of sequence identity, appears to have evolved a different function to its host ortholog. This protein is an active Ub conjugating enzyme that binds to a panel of human E3 Ub ligases with higher affinity than the host ortholog. Furthermore, when expressed in skeletal muscle cells in culture, its presence uniquely leads to the stabilization of muscle-specific proteins via the downregulation of their ubiquitination.

The binary response of the GAL/MEL genetic switch of Saccharomyces cerevisiae is critically dependent on Gal80p-Gal4p interaction

Studies on the Saccharomyces cerevisiae GAL/MEL genetic switch have revealed that its bistability is dependent on ultrasensitivity that can be altered or abolished by disabling different combinations of nested feedback loops. In contrast, we have previously demonstrated that weakening of the interaction between Gal80p and Gal4p alone is sufficient to abolish the ultrasensitivity (Das Adhikari et al. 2014). Here, we demonstrate that altering the epistatic interaction between Gal80p and Gal4p also abolishes the bistability, and the switch response to galactose becomes graded instead of binary. However, the GAL/MEL switch of wild-type and epistatically altered strains responded in a graded fashion to melibiose. The properties of the epistatically altered strain resemble Kluyveromyces lactis, which separated from the Saccharomyces lineage 100 mya before whole-genome duplication (WGD). Based on the results reported here, we propose that epistatic interactions played a crucial role in the evolution of the fine regulation of S. cerevisiae GAL/MEL switch following WGD.

High-Throughput Identification of Cis-Regulatory Rewiring Events in Yeast

A coregulated module of genes ("regulon") can have evolutionarily conserved expression patterns and yet have diverged upstream regulators across species. For instance, the ribosomal genes regulon is regulated by the transcription factor (TF) TBF1 in Candida albicans, while in Saccharomyces cerevisiae it is regulated by RAP1. Only a handful of such rewiring events have been established, and the prevalence or conditions conducive to such events are not well known. Here, we develop a novel probabilistic scoring method to comprehensively screen for regulatory rewiring within regulons across 23 yeast species. Investigation of 1,713 regulons and 176 TFs yielded 5,353 significant rewiring events at 5% false discovery rate (FDR). Besides successfully recapitulating known rewiring events, our analyses also suggest TF candidates for certain processes reported to be under distinct regulatory controls in S. cerevisiae and C. albicans, for which the implied regulators are not known: 1) Oxidative stress response (Sc-MSN2 to Ca-FKH2) and 2) nutrient modulation (Sc-RTG1 to Ca-GCN4/Ca-UME6). Furthermore, a stringent screen to detect TF rewiring at individual genes identified 1,446 events at 10% FDR. Overall, these events are supported by strong coexpression between the predicted regulator and its target gene(s) in a species-specific fashion (>50-fold). Independent functional analyses of rewiring TF pairs revealed greater functional interactions and shared biological processes between them (P = 1 × 10(-3)).Our study represents the first comprehensive assessment of regulatory rewiring; with a novel approach that has generated a unique high-confidence resource of several specific events, suggesting that evolutionary rewiring is relatively frequent and may be a significant mechanism of regulatory innovation.

Examining the virulence of Candida albicans transcription factor mutants using Galleria mellonella and mouse infection models

The aim of the present study was to identify Candida albicans transcription factors (TFs) involved in virulence. Although mice are considered the gold-standard model to study fungal virulence, mini-host infection models have been increasingly used. Here, barcoded TF mutants were first screened in mice by pools of strains and fungal burdens (FBs) quantified in kidneys. Mutants of unannotated genes which generated a kidney FB significantly different from that of wild-type were selected and individually examined in Galleria mellonella. In addition, mutants that could not be detected in mice were also tested in G. mellonella. Only 25% of these mutants displayed matching phenotypes in both hosts, highlighting a significant discrepancy between the two models. To address the basis of this difference (pool or host effects), a set of 19 mutants tested in G. mellonella were also injected individually into mice. Matching FB phenotypes were observed in 50% of the cases, highlighting the bias due to host effects. In contrast, 33.4% concordance was observed between pool and single strain infections in mice, thereby highlighting the bias introduced by the "pool effect." After filtering the results obtained from the two infection models, mutants for MBF1 and ZCF6 were selected. Independent marker-free mutants were subsequently tested in both hosts to validate previous results. The MBF1 mutant showed impaired infection in both models, while the ZCF6 mutant was only significant in mice infections. The two mutants showed no obvious in vitro phenotypes compared with the wild-type, indicating that these genes might be specifically involved in in vivo adapt.

Examining the evolution of the regulatory circuit controlling secondary metabolism and development in the fungal genus Aspergillus

Filamentous fungi produce diverse secondary metabolites (SMs) essential to their ecology and adaptation. Although each SM is typically produced by only a handful of species, global SM production is governed by widely conserved transcriptional regulators in conjunction with other cellular processes, such as development. We examined the interplay between the taxonomic narrowness of SM distribution and the broad conservation of global regulation of SM and development in Aspergillus, a diverse fungal genus whose members produce well-known SMs such as penicillin and gliotoxin. Evolutionary analysis of the 2,124 genes comprising the 262 SM pathways in four Aspergillus species showed that most SM pathways were species-specific, that the number of SM gene orthologs was significantly lower than that of orthologs in primary metabolism, and that the few conserved SM orthologs typically belonged to non-homologous SM pathways. RNA sequencing of two master transcriptional regulators of SM and development, veA and mtfA, showed that the effects of deletion of each gene, especially veA, on SM pathway regulation were similar in A. fumigatus and A. nidulans, even though the underlying genes and pathways regulated in each species differed. In contrast, examination of the role of these two regulators in development, where 94% of the underlying genes are conserved in both species showed that whereas the role of veA is conserved, mtfA regulates development in the homothallic A. nidulans but not in the heterothallic A. fumigatus. Thus, the regulation of these highly conserved developmental genes is divergent, whereas-despite minimal conservation of target genes and pathways-the global regulation of SM production is largely conserved. We suggest that the evolution of the transcriptional regulation of secondary metabolism in Aspergillus represents a novel type of regulatory circuit rewiring and hypothesize that it has been largely driven by the dramatic turnover of the target genes involved in the process.

Perturbation of the interaction between Gal4p and Gal80p of the Saccharomyces cerevisiae GAL switch results in altered responses to galactose and glucose

In S. cerevisiae, following the Whole Genome Duplication (WGD), GAL1-encoded galactokinase retained its signal transduction function but lost basal expression. On the other hand, its paralogue GAL3, lost kinase activity but retained its signalling function and basal expression, thus making it indispensable for the rapid induction of the S. cerevisiae GAL switch. However, a gal3Δ strain exhibits delayed growth kinetics due to the redundant signalling function of GAL1. The subfunctionalization between the paralogues GAL1 and GAL3 is due to expression divergence and is proposed to be due to the alteration in the Upstream Activating Sequences (UASG ). We demonstrate that the GAL switch becomes independent of GAL3 by altering the interaction between Gal4p and Gal80p without altering the configuration of UASG . In addition to the above, the altered switch of S. cerevisiae loses ultrasensitivity and stringent glucose repression. These changes caused an increase in fitness in the disaccharide melibiose at the expense of a decrease in fitness in galactose. The above altered features of the ScGAL switch are similar to the features of the GAL switch of K. lactis that diverged from S. cerevisiae before the WGD.

Repeated loss of an anciently horizontally transferred gene cluster in Botrytis

At least five of the six genes of the bikaverin secondary metabolic gene cluster were shown to have undergone horizontal transfer (HGT) from a Fusarium donor to the Botrytis lineage. Of these five, two enzyme-encoding genes are found as pseudogenes in B. cinerea whereas two regulatory genes and the transporter remain intact. To reconstruct the evolutionary events leading to decay of this gene cluster and infer a more precise timing of its transfer, we examined the genomes of nine additional broadly sampled Botrytis species. We found evidence that a Botrytis ancestor acquired the entire gene cluster through an ancient HGT that occurred before the diversification of the genus. During the subsequent evolution and diversification of the genus, four of the 10 genomes appear to have lost the gene cluster, while in the other six the cluster is in various stages of degeneration. Across the Botrytis genomes, the modes of gene decay in the cluster differed between enzyme-encoding genes, which had higher rates of transition to or retention of pseudogenes and were universally inactivated, and regulatory genes (particularly the non-pathway-specific regulator bik4), which more frequently appeared intact. Consistent with these results, the regulatory genes bik4 and bik5 showed stronger evidence of transcriptional expression than other bikaverin genes under multiple conditions in B. cinerea. These results could be explained by pleiotropy in the bikaverin regulatory genes either through rewiring or their interaction with more central pathways or by constraints on the order of gene loss driven by the intrinsic toxicity of the pathway. Our finding that most of the bikaverin pathway genes have been lost or pseudogenized in these Botrytis genomes suggests that the incidence of HGT of gene cluster-encoded metabolic pathways might be higher than what is possible to be inferred from isolated genome analyses.

Upregulation of galactose metabolic pathway by N-acetylglucosamine induced endogenous synthesis of galactose in Candida albicans

N-Acetylglucosamine (GlcNAc) is an important signaling molecule that plays multiple roles in Candida albicans. Induction of galactose metabolic pathway by GlcNAc is an intriguing aspect of C. albicans biology. In order to investigate the role of galactose metabolic genes (GAL genes) in presence of GlcNAc, we created knockouts of galactokinase (GAL1) and UDP galactose epimerase (GAL10) genes. These mutants failed to grow on galactose and also showed lower growth rate in presence of GlcNAc. Interestingly, expression of GAL genes in presence of GlcNAc was higher in gal1Δ strain relative to that of wild type strain. Moreover, no GlcNAc induced upregulation of GAL genes was observed in the gal10Δ strain suggesting that UDP galactose epimerase is essential for GlcNAc induced activation of GAL genes. GlcNAc induced expression of GAL genes was also investigated in GlcNAc metabolic pathway triple mutant N216 (hxk1Δ nag1Δ dac1Δ). Interestingly, in this mutant the GAL genes are neither induced nor repressed and remain derepressed as found on a neutral carbon source such as glycerol, suggesting that catabolism of GlcNAc play an important role in the expression of GAL genes. GC/MS analysis of derivatized metabolites revealed a significant accumulation of galactose in the gal1Δ strain while no galactose was detected in gal10Δ and N216 strain. Solution-state NMR spectroscopy using N-acetyl-¹³C₁-glucosamine confirmed the flow of ¹³C label from GlcNAc to galactose. Thus, internal galactose synthesized via UDP galactose pathway from GlcNAc metabolites acts as the inducer of GAL genes in presence of GlcNAc.

Genome-Wide Scale-Free Network Inference for Candida albicans

Discovery of essential genes in pathogenic organisms is an important step in the development of new medication. Despite a growing number of genome data available, little is known about C. albicans, a major fungal pathogen. Most of the human population carries C. albicans as commensal, but it can cause systemic infection that may lead to the death of the host if the immune system has deteriorated. In many organisms central nodes in the interaction network (hubs) play a crucial role for information and energy transport. Knock-outs of such hubs often lead to lethal phenotypes making them interesting drug targets. To identify these central genes via topological analysis, we inferred gene regulatory networks that are sparse and scale-free. We collected information from various sources to complement the limited expression data available. We utilized a linear regression algorithm to infer genome-wide gene regulatory interaction networks. To evaluate the predictive power of our approach, we used an automated text-mining system that scanned full-text research papers for known interactions. With the help of the compendium of known interactions, we also optimize the influence of the prior knowledge and the sparseness of the model to achieve the best results. We compare the results of our approach with those of other state-of-the-art network inference methods and show that we outperform those methods. Finally we identify a number of hubs in the genome of the fungus and investigate their biological relevance.

Horizontal transfer and death of a fungal secondary metabolic gene cluster

A cluster composed of four structural and two regulatory genes found in several species of the fungal genus Fusarium (class Sordariomycetes) is responsible for the production of the red pigment bikaverin. We discovered that the unrelated fungus Botrytis cinerea (class Leotiomycetes) contains a cluster of five genes that is highly similar in sequence and gene order to the Fusarium bikaverin cluster. Synteny conservation, nucleotide composition, and phylogenetic analyses of the cluster genes indicate that the B. cinerea cluster was acquired via horizontal transfer from a Fusarium donor. Upon or subsequent to the transfer, the B. cinerea gene cluster became inactivated; one of the four structural genes is missing, two others are pseudogenes, and the fourth structural gene shows an accelerated rate of nonsynonymous substitutions along the B. cinerea lineage, consistent with relaxation of selective constraints. Interestingly, the bik4 regulatory gene is still intact and presumably functional, whereas bik5, which is a pathway-specific regulator, also shows a mild but significant acceleration of evolutionary rate along the B. cinerea lineage. This selective preservation of the bik4 regulator suggests that its conservation is due to its likely involvement in other non-bikaverin-related biological processes in B. cinerea. Thus, in addition to novel metabolism, horizontal transfer of wholesale metabolic gene clusters might also be contributing novel regulation.

Modeling the evolution of a classic genetic switch

Background

The regulatory network underlying the yeast galactose-use pathway has emerged as a model system for the study of regulatory network evolution. Evidence has recently been provided for adaptive evolution in this network following a whole genome duplication event. An ancestral gene encoding a bi-functional galactokinase and co-inducer protein molecule has become subfunctionalized as paralogous genes (GAL1 and GAL3) in Saccharomyces cerevisiae, with most fitness gains being attributable to changes in cis-regulatory elements. However, the quantitative functional implications of the evolutionary changes in this regulatory network remain unexplored.

Results

We develop a modeling framework to examine the evolution of the GAL regulatory network. This enables us to translate molecular changes in the regulatory network to changes in quantitative network function. We computationally reconstruct an inferred ancestral version of the network and trace the evolutionary paths in the lineage leading to S. cerevisiae. We explore the evolutionary landscape of possible regulatory networks and find that the operation of intermediate networks leading to S. cerevisiae differs substantially depending on the order in which evolutionary changes accumulate; in particular, we systematically explore evolutionary paths and find that some network features cannot be optimized simultaneously.

Conclusions

We find that a computational modeling approach can be used to analyze the evolution of a well-studied regulatory network. Our results are consistent with several experimental studies of the evolutionary of the GAL regulatory network, including increased fitness in Saccharomyces due to duplication and adaptive regulatory divergence. The conceptual and computational tools that we have developed may be applicable in further studies of regulatory network evolution.

Azole drugs are imported by facilitated diffusion in Candida albicans and other pathogenic fungi

Despite the wealth of knowledge regarding the mechanisms of action and the mechanisms of resistance to azole antifungals, very little is known about how the azoles are imported into pathogenic fungal cells. Here the in-vitro accumulation and import of Fluconazole (FLC) was examined in the pathogenic fungus, Candida albicans. In energized cells, FLC accumulation correlates inversely with expression of ATP-dependent efflux pumps. In de-energized cells, all strains accumulate FLC, suggesting that FLC import is not ATP-dependent. The kinetics of import in de-energized cells displays saturation kinetics with a K_m of 0.64 uM and V_max of 0.0056 pmol/min/10⁸ cells, demonstrating that FLC import proceeds via facilitated diffusion through a transporter rather than passive diffusion. Other azoles inhibit FLC import on a mole/mole basis, suggesting that all azoles utilize the same facilitated diffusion mechanism. An analysis of related compounds indicates that competition for azole import depends on an aromatic ring and an imidazole or triazole ring together in one molecule. Import of FLC by facilitated diffusion is observed in other fungi, including Cryptococcus neoformans, Saccharomyces cerevisiae, and Candida krusei, indicating that the mechanism of transport is conserved among fungal species. FLC import was shown to vary among Candida albicans resistant clinical isolates, suggesting that altered facilitated diffusion may be a previously uncharacterized mechanism of resistance to azole drugs.

Azole antifungals are used to treat a wide variety of fungal infections of humans, animals and plants. A great deal is known about how the azoles interact with their target enzyme within fungal cells and how the azoles are exported from the fungal cell through various efflux pumps. Altered interactions with the target enzyme and altered efflux pump expression are common mechanisms of azole resistance in fungi. However, the mechanism by which azoles enter a fungal cell is not clear—many have assumed that azoles passively diffuse into the cell. This study demonstrates that azoles are not passively diffused, or actively pumped, into the cell. Instead, azoles are imported by facilitated diffusion, mediated by a transporter. Facilitated diffusion of azoles is saturable. All clinically important azoles, and many structurally related compounds, compete for FLC import, while structurally unrelated drugs do not compete. Azole import by facilitated diffusion is shown in four species of fungi, suggesting that it is common for most if not all fungi. Altered facilitated diffusion is observed in a collection of clinical isolates, suggesting that altered import is a previously uncharacterized mechanism of resistance.

From elements to modules: regulatory evolution in Ascomycota fungi

Regulatory divergence is likely a major driving force in evolution. Comparative transcriptomics provides a new glimpse into the evolution of gene regulation. Ascomycota fungi are uniquely suited among eukaryotes for studies of regulatory evolution, because of broad phylogenetic scope, many sequenced genomes, and facility of genomic analysis. Here we review the substantial divergence in gene expression in Ascomycota and how this is reconciled with the modular organization of transcriptional networks. We show that flexibility and redundancy in both cis-regulation and trans-regulation can lead to changes from altered expression of single genes to wholesale rewiring of regulatory modules. Redundancy thus emerges as a major driving force facilitating expression divergence while preserving the coherent functional organization of a transcriptional response.

An integrated approach to identifying cis-regulatory modules in the human genome

In eukaryotic genomes, it is challenging to accurately determine target sites of transcription factors (TFs) by only using sequence information. Previous efforts were made to tackle this task by considering the fact that TF binding sites tend to be more conserved than other functional sites and the binding sites of several TFs are often clustered. Recently, ChIP-chip and ChIP-sequencing experiments have been accumulated to identify TF binding sites as well as survey the chromatin modification patterns at the regulatory elements such as promoters and enhancers. We propose here a hidden Markov model (HMM) to incorporate sequence motif information, TF-DNA interaction data and chromatin modification patterns to precisely identify cis-regulatory modules (CRMs). We conducted ChIP-chip experiments on four TFs, CREB, E2F1, MAX, and YY1 in 1% of the human genome. We then trained a hidden Markov model (HMM) to identify the labels of the CRMs by incorporating the sequence motifs recognized by these TFs and the ChIP-chip ratio. Chromatin modification data was used to predict the functional sites and to further remove false positives. Cross-validation showed that our integrated HMM had a performance superior to other existing methods on predicting CRMs. Incorporating histone signature information successfully penalized false prediction and improved the whole performance. The dataset we used and the software are available at http://nash.ucsd.edu/CIS/.