Conservation and evolution of cis-regulatory systems in ascomycete fungi

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Con7 is a key transcription regulator for conidiogenesis in the plant pathogenic fungus Fusarium graminearum

The mycelium of the plant pathogenic fungus Fusarium graminearum exhibits distinct structures for vegetative growth, asexual sporulation, sexual development, virulence, and chlamydospore formation. These structures are vital for the survival and pathogenicity of the fungus, necessitating precise regulation based on environmental cues. Initially identified in Magnaporthe oryzae, the transcription factor Con7p regulates conidiation and infection-related morphogenesis, but not vegetative growth. We characterized the Con7p ortholog FgCon7, and deletion of FgCON7 resulted in severe defects in conidium production, virulence, sexual development, and vegetative growth. The mycelia of the deletion mutant transformed into chlamydospore-like structures with high chitin level accumulation. Notably, boosting FgABAA expression partially alleviated developmental issues in the FgCON7 deletion mutant. Chromatin immunoprecipitation (ChIP)-quantitative PCR (qPCR) analysis confirmed a direct genetic link between FgABAA and FgCON7. Furthermore, the chitin synthase gene Fg6550 (FGSG_06550) showed significant upregulation in the FgCON7 deletion mutant, and altering FgCON7 expression affected cell wall integrity. Further research will focus on understanding the behavior of the chitin synthase gene and its regulation by FgCon7 in F. graminearum. This study contributes significantly to our understanding of the genetic pathways that regulate hyphal differentiation and conidiation in this plant pathogenic fungus.

IMPORTANCE

The ascomycete fungus Fusarium graminearum is the primary cause of head blight disease in wheat and barley, as well as ear and stalk rot in maize. Given the importance of conidia and ascospores in the disease cycle of F. graminearum, precise spatiotemporal regulation of these biological processes is crucial. In this study, we characterized the Magnaporthe oryzae Con7p ortholog and discovered that FgCon7 significantly influences various crucial aspects of fungal development and pathogenicity. Notably, overexpression of FgABAA partially restored developmental defects in the FgCON7 deletion mutant. ChIP-qPCR analysis confirmed a direct genetic link between FgABAA and FgCON7. Furthermore, our research revealed a clear correlation between FgCon7 and chitin accumulation and the expression of chitin synthase genes. These findings offer valuable insights into the genetic mechanisms regulating conidiation and the significance of mycelial differentiation in this plant pathogenic fungus.

Gene loss and cis-regulatory novelty shaped core histone gene evolution in the apiculate yeast Hanseniaspora uvarum

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage.

An in Silico Approach to Identifying TF Binding Sites: Analysis of the Regulatory Regions of BUSCO Genes from Fungal Species in the Ceratocystidaceae Family

Transcriptional regulation controls gene expression through regulatory promoter regions that contain conserved sequence motifs. These motifs, also known as regulatory elements, are critically important to expression, which is driving research efforts to identify and characterize them. Yeasts have been the focus of such studies in fungi, including in several in silico approaches. This study aimed to determine whether in silico approaches could be used to identify motifs in the Ceratocystidaceae family, and if present, to evaluate whether these correspond to known transcription factors. This study targeted the 1000 base-pair region upstream of the start codon of 20 single-copy genes from the BUSCO dataset for motif discovery. Using the MEME and Tomtom analysis tools, conserved motifs at the family level were identified. The results show that such in silico approaches could identify known regulatory motifs in the Ceratocystidaceae and other unrelated species. This study provides support to ongoing efforts to use in silico analyses for motif discovery.

From a genome assembly to full regulatory network prediction: the case study of Rhodotorula toruloides putative Haa1-regulon

Numerous genomes are sequenced and made available to the community through the NCBI portal. However, and, unlike what happens for gene function annotation, annotation of promoter sequences and the underlying prediction of regulatory associations is mostly unavailable, severely limiting the ability to interpret genome sequences in a functional genomics perspective. Here we present an approach where one can download a genome of interest from NCBI in the GenBank Flat File (.gbff) format and, with a minimum set of commands, have all the information parsed, organized and made available through the platform web interface. Also, the new genomes are compared with a given genome of reference in search of homologous genes, shared regulatory elements and predicted transcription associations. We present this approach within the context of Community YEASTRACT of the YEASTRACT + portal, thus benefiting from immediate access to all the comparative genomics queries offered in the YEASTRACT + portal. Besides the yeast community, other communities can install the platform independently, without any constraints. In this work, we exemplify the usefulness of the presented tool, within Community YEASTRACT, in constructing a dedicated database and analysing the genome of the highly promising oleaginous red yeast species Rhodotorula toruloides currently poorly studied at the genome and transcriptome levels and with limited genome editing tools. Regulatory prediction is based on the conservation of promoter sequences and available regulatory networks. The case-study examined is focused on the Haa1 transcription factor—a key regulator of yeast resistance to acetic acid, an important inhibitor of industrial bioconversion of lignocellulosic hydrolysates. The new tool described here led to the prediction of a RtHaa1 regulon with expected impact in the optimization of R. toruloides robustness for lignocellulosic and pectin-rich residue biorefinery processes.

The master regulator MAT1-1-1 of fungal mating binds to its targets via a conserved motif in the human pathogen Aspergillus fumigatus

Mating-type transcription factors are master regulators of sexually related signal transduction pathways in fungi; however, their recognition of specific DNA sequences from target genes is widely undetermined. Here, we identified and characterized the DNA-binding sequence of the MAT1-1-1 alpha-box domain transcription factor from the human pathogen Aspergillus fumigatus. In order to explore MAT1-1-1 DNA-binding targets, we used the previously reported MAT1-1-1 binding motif from Penicillium chrysogenum, in a bioinformatics approach. We identified 18 A. fumigatus genes carrying the MAT1.1 sequence in their upstream region, among them genes for the α-pheromone precursor (PpgA), G-protein-coupled pheromone receptor (PreA), and for TomA, an unidentified protein. To validate our prediction further, quantification of transcript levels showed a decrease in expression of ppgA, tomA, and others in a MAT1-1 deletion strain. For a functional analysis of the binding sites, truncated variants of the A. fumigatus MAT1-1-1 gene were introduced into Escherichia coli for heterologous expression. The yield of recombinant protein was further optimized for the AfMAT1-1-178-235 variant that harbors an extended alpha-box domain. AfMAT1-1-178-235 bound to a subset of the most strongly upregulated genes: ppgA, preA, and tomA. The DNA-binding specificity was confirmed by testing mutated binding sequences, as well as performing competition experiments with specific and non-specific sequences. Finally, equilibrium dissociation constants of 1.83 ± 0.1 and 1.45 ± 0.26 µM were determined for AfMAT1-1-178-235 and fusion protein GST-AfMAT1-1-178-235. Collectively, these findings provide further insights into AfMAT1-1-1-mediated gene expression and imply that alpha-box domain regulators from other members of Eurotiales control fungal development in a conserved manner.

Function Analysis of MBF1, a Factor Involved in the Response to Amino Acid Starvation and Virulence in Candida albicans

Candida albicans is a commensal of human mucosae, but also one of the most common fungal pathogens of humans. Systemic infections caused by this fungus, mostly affecting immunocompromised patients, are associated to fatality rates as high as 50% despite the available treatments. In order to improve this situation, it is necessary to fully understand how C. albicans is able to cause disease and how it copes with the host defenses. Our previous studies have revealed the importance of the C. albicans gene MBF1 in virulence and ability to colonize internal organs of mammalian and insect hosts. MBF1 encodes a putative transcriptional regulator, and as such it likely has an impact in the regulation of C. albicans gene expression during host infection. Here, recent advances in RNA-seq technologies were used to obtain a detailed analysis of the impact of MBF1 on C. albicans gene expression both in vitro and during infection. MBF1 was involved in the regulation of several genes with a role in glycolysis and response to stress, particularly to nutritional stress. We also investigated whether an interaction existed between MBF1 and GCN4, a master regulator of response to starvation, and found that both genes were needed for resistance to amino acid starvation, suggesting some level of interaction between the two. Reinforcing this idea, we showed that the proteins encoded by both genes could interact. Consistent with the role of MBF1 in virulence, we also established that GCN4 was necessary for virulence in the mouse model of systemic infection as well as in the Galleria mellonella infection model.

Transcriptome-wide differences between Saccharomyces cerevisiae and Saccharomyces cerevisiae var. boulardii: Clues on host survival and probiotic activity based on promoter sequence variability

Although Saccharomyces cerevisiae and S. cerevisiae var. boulardii share more than 95% genome sequence homology, only S. cerevisiae var. boulardii displays probiotic activity. In this study, the transcriptomic differences exhibited by S. cerevisiae and S. cerevisiae var. boulardii in intestinal like medium were evaluated. S. cerevisiae was found to display stress response overexpression, consistent with higher ability of S. cerevisiae var. boulardii to survive within the human host, while S. cerevisiae var. boulardii exhibited transcriptional patterns associated with probiotic activity, suggesting increased acetate biosynthesis. Resorting to the creation of a S. cerevisiae var. boulardii genomic database within Yeastract+, a possible correlation between loss or gain of transcription factor binding sites in S. cerevisiae var. boulardii promoters and the transcriptomic pattern is discussed. This study suggests that S. cerevisiae var. boulardii probiotic activity, when compared to S. cerevisiae, relies, at least partially, on differential expression regulation, based on promoter variability.

Yeast Rpn4 Links the Proteasome and DNA Repair via RAD52 Regulation

Environmental and intracellular factors often damage DNA, but multiple DNA repair pathways maintain genome integrity. In yeast, the 26S proteasome and its transcriptional regulator and substrate Rpn4 are involved in DNA damage resistance. Paradoxically, while proteasome dysfunction may induce hyper-resistance to DNA-damaging agents, Rpn4 malfunction sensitizes yeasts to these agents. Previously, we proposed that proteasome inhibition causes Rpn4 stabilization followed by the upregulation of Rpn4-dependent DNA repair genes and pathways. Here, we aimed to elucidate the key Rpn4 targets responsible for DNA damage hyper-resistance in proteasome mutants. We impaired the Rpn4-mediated regulation of candidate genes using the CRISPR/Cas9 system and tested the sensitivity of mutant strains to 4-NQO, MMS and zeocin. We found that the separate or simultaneous deregulation of 19S or 20S proteasome subcomplexes induced MAG1, DDI1, RAD23 and RAD52 in an Rpn4-dependent manner. Deregulation of RAD23, DDI1 and RAD52 sensitized yeast to DNA damage. Genetic, epigenetic or dihydrocoumarin-mediated RAD52 repression restored the sensitivity of the proteasome mutants to DNA damage. Our results suggest that the Rpn4-mediated overexpression of DNA repair genes, especially RAD52, defines the DNA damage hyper-resistant phenotype of proteasome mutants. The developed yeast model is useful for characterizing drugs that reverse the DNA damage hyper-resistance phenotypes of cancers.

Gene Regulatory Networks of Penicillium echinulatum 2HH and Penicillium oxalicum 114-2 Inferred by a Computational Biology Approach

Penicillium echinulatum 2HH and Penicillium oxalicum 114-2 are well-known cellulase fungal producers. However, few studies addressing global mechanisms for gene regulation of these two important organisms are available so far. A recent finding that the 2HH wild-type is closely related to P. oxalicum leads to a combined study of these two species. Firstly, we provide a global gene regulatory network for P. echinulatum 2HH and P. oxalicum 114-2, based on TF-TG orthology relationships, considering three related species with well-known regulatory interactions combined with TFBSs prediction. The network was then analyzed in terms of topology, identifying TFs as hubs, and modules. Based on this approach, we explore numerous identified modules, such as the expression of cellulolytic and xylanolytic systems, where XlnR plays a key role in positive regulation of the xylanolytic system. It also regulates positively the cellulolytic system by acting indirectly through the cellodextrin induction system. This remarkable finding suggests that the XlnR-dependent cellulolytic and xylanolytic regulatory systems are probably conserved in both P. echinulatum and P. oxalicum. Finally, we explore the functional congruency on the genes clustered in terms of communities, where the genes related to cellular nitrogen, compound metabolic process and macromolecule metabolic process were the most abundant. Therefore, our approach allows us to confer a degree of accuracy regarding the existence of each inferred interaction.

Replicative aging is associated with loss of genetic heterogeneity from extrachromosomal circular DNA in Saccharomyces cerevisiae

Circular DNA can arise from all parts of eukaryotic chromosomes. In yeast, circular ribosomal DNA (rDNA) accumulates dramatically as cells age, however little is known about the accumulation of other chromosome-derived circles or the contribution of such circles to genetic variation in aged cells. We profiled circular DNA in Saccharomyces cerevisiae populations sampled when young and after extensive aging. Young cells possessed highly diverse circular DNA populations but 94% of the circular DNA were lost after ∼15 divisions, whereas rDNA circles underwent massive accumulation to >95% of circular DNA. Circles present in both young and old cells were characterized by replication origins including circles from unique regions of the genome and repetitive regions: rDNA and telomeric Y' regions. We further observed that circles can have flexible inheritance patterns: [HXT6/7circle] normally segregates to mother cells but in low glucose is present in up to 50% of cells, the majority of which must have inherited this circle from their mother. Interestingly, [HXT6/7circle] cells are eventually replaced by cells carrying stable chromosomal HXT6 HXT6/7 HXT7 amplifications, suggesting circular DNAs are intermediates in chromosomal amplifications. In conclusion, the heterogeneity of circular DNA offers flexibility in adaptation, but this heterogeneity is remarkably diminished with age.

Candida glabrata Transcription Factor Rpn4 Mediates Fluconazole Resistance through Regulation of Ergosterol Biosynthesis and Plasma Membrane Permeability

The ability to acquire azole resistance is an emblematic trait of the fungal pathogen Candida glabrata Understanding the molecular basis of azole resistance in this pathogen is crucial for designing more suitable therapeutic strategies. This study shows that the C. glabrata transcription factor (TF) CgRpn4 is a determinant of azole drug resistance. RNA sequencing during fluconazole exposure revealed that CgRpn4 regulates the expression of 212 genes, activating 80 genes and repressing, likely in an indirect fashion, 132 genes. Targets comprise several proteasome and ergosterol biosynthesis genes, including ERG1, ERG2, ERG3, and ERG11 The localization of CgRpn4 to the nucleus increases upon fluconazole stress. Consistent with a role in ergosterol and plasma membrane homeostasis, CgRpn4 is required for the maintenance of ergosterol levels upon fluconazole stress, which is associated with a role in the upkeep of cell permeability and decreased intracellular fluconazole accumulation. We provide evidence that CgRpn4 directly regulates ERG11 expression through the TTGCAAA binding motif, reinforcing the relevance of this regulatory network in azole resistance. In summary, CgRpn4 is a new regulator of the ergosterol biosynthesis pathway in C. glabrata, contributing to plasma membrane homeostasis and, thus, decreasing azole drug accumulation.

Mutational bias and the protein code shape the evolution of splicing enhancers

Exonic splicing enhancers (ESEs) are enriched in exons relative to introns and bind splicing activators. This study considers a fundamental question of co-evolution: How did ESE motifs become enriched in exons prior to the evolution of ESE recognition? We hypothesize that the high exon to intron motif ratios necessary for ESE function were created by mutational bias coupled with purifying selection on the protein code. These two forces retain certain coding motifs in exons while passively depleting them from introns. Through the use of simulations, genomic analyses, and high throughput splicing assays, we confirm the key predictions of this hypothesis, including an overlap between protein and splicing information in ESEs. We discuss the implications of mutational bias as an evolutionary driver in other cis-regulatory systems.

Splicing is regulated by cis-acting elements in pre-mRNAs such as exonic or intronic splicing enhancers and silencers. Here the authors show that exonic splicing enhancers are enriched in exons compared to introns due to mutational bias coupled with purifying selection on the protein code.

Mitochondrial Biogenesis Is Positively Regulated by Casein Kinase I Hrr25 Through Phosphorylation of Puf3 in Saccharomyces cerevisiae

Mitochondrial biogenesis requires coordinated expression of genes encoding mitochondrial proteins, which in Saccharomyces cerevisiae is achieved in part via post-transcriptional control by the Pumilio RNA-binding domain protein Puf3 Puf3 binds to the 3'-UTR of many messenger RNAs (mRNAs) that encode mitochondrial proteins, regulating their turnover, translation, and/or mitochondrial targeting. Puf3 hyperphosphorylation correlates with increased mitochondrial biogenesis; however, the kinase responsible for Puf3 phosphorylation is unclear. Here, we show that the casein kinase I protein Hrr25 negatively regulates Puf3 by mediating its phosphorylation. An hrr25 mutation results in reduced phosphorylation of Puf3 in vivo and a puf3 deletion mutation reverses growth defects of hrr25 mutant cells grown on medium with a nonfermentable carbon source. We show that Hrr25 directly phosphorylates Puf3, and that the interaction between Puf3 and Hrr25 is mediated through the N-terminal domain of Puf3 and the kinase domain of Hrr25 We further found that an hrr25 mutation reduces GFP expression from GFP reporter constructs carrying the 3'-UTR of Puf3 targets. Downregulation of GFP expression due to an hrr25 mutation can be reversed either by puf3Δ or by mutations to the Puf3-binding sites in the 3'-UTR of the GFP reporter constructs. Together, our data indicate that Hrr25 is a positive regulator of mitochondrial biogenesis by phosphorylating Puf3 and inhibiting its function in downregulating target mRNAs encoding mitochondrial proteins.

YEASTRACT+: a portal for cross-species comparative genomics of transcription regulation in yeasts

The YEASTRACT+ information system (http://YEASTRACT-PLUS.org/) is a wide-scope tool for the analysis and prediction of transcription regulatory associations at the gene and genomic levels in yeasts of biotechnological or human health relevance. YEASTRACT+ is a new portal that integrates the previously existing YEASTRACT (http://www.yeastract.com/) and PathoYeastract (http://pathoyeastract.org/) databases and introduces the NCYeastract (Non-Conventional Yeastract) database (http://ncyeastract.org/), focused on the so-called non-conventional yeasts. The information in the YEASTRACT database, focused on Saccharomyces cerevisiae, was updated. PathoYeastract was extended to include two additional pathogenic yeast species: Candida parapsilosis and Candida tropicalis. Furthermore, the NCYeastract database was created, including five biotechnologically relevant yeast species: Zygosaccharomyces baillii, Kluyveromyces lactis, Kluyveromyces marxianus, Yarrowia lipolytica and Komagataella phaffii. The YEASTRACT+ portal gathers 289 706 unique documented regulatory associations between transcription factors (TF) and target genes and 420 DNA binding sites, considering 247 TFs from 10 yeast species. YEASTRACT+ continues to make available tools for the prediction of the TFs involved in the regulation of gene/genomic expression. In this release, these tools were upgraded to enable predictions based on orthologous regulatory associations described for other yeast species, including two new tools for cross-species transcription regulation comparison, based on multi-species promoter and TF regulatory network analyses.

Extreme divergence between one-to-one orthologs: the structure of N15 Cro bound to operator DNA and its relationship to the λ Cro complex

The gene cro promotes lytic growth of phages through binding of Cro protein dimers to regulatory DNA sites. Most Cro proteins are one-to-one orthologs, yet their sequence, structure and binding site sequences are quite divergent across lambdoid phages. We report the cocrystal structure of bacteriophage N15 Cro with a symmetric consensus site. We contrast this complex with an orthologous structure from phage λ, which has a dissimilar binding site sequence and a Cro protein that is highly divergent in sequence, dimerization interface and protein fold. The N15 Cro complex has less DNA bending and smaller DNA-induced changes in protein structure. N15 Cro makes fewer direct contacts and hydrogen bonds to bases, relying mostly on water-mediated and Van der Waals contacts to recognize the sequence. The recognition helices of N15 Cro and λ Cro make mostly nonhomologous and nonanalogous contacts. Interface alignment scores show that half-site binding geometries of N15 Cro and λ Cro are less similar to each other than to distantly related CI repressors. Despite this divergence, the Cro family shows several code-like protein–DNA sequence covariations. In some cases, orthologous genes can achieve a similar biological function using very different specific molecular interactions.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

Evolutionary Transition of GAL Regulatory Circuit from Generalist to Specialist Function in Ascomycetes

The Gal4 transcription factor (TF) controls gene expression by binding the DNA sequence motif CGG(N₁₁)CCG. Well studied versions regulate metabolism of glucose in Candida albicans and galactose in Saccharomyces cerevisiae. Gal4 is also found within Aspergillus species and shows a wide range of potential binding targets. Members of the CTG clade that reassigned CUG codons from leucine to serine lack the Gal80 binding domain of Gal4, and they use the TF to regulate only glycolytic genes. In this clade, the galactose catabolic pathway (also known as the Leloir pathway) genes are regulated by Rtg1/Rtg3. In the WGD species, the complete Gal4/Gal80 module is limited to regulation of the Leloir pathway, while glycolysis is controlled by Gcr1/Gcr2. This shows a switch of Gal4 from a generalist to a specialist within the ascomycetes, and the split of glucose and galactose metabolism into distinct regulatory circuits.

Unique gene Pmhyp controlling melanization of pycnidia in Phoma medicaginis

Phoma medicaginis (syn. Ascochyta medicaginicola Qchen & L. Cai) causes spring black stem and leaf spot of alfalfa and the model legume Medicago truncatula. Phoma medicaginis produces uninucleate conidia in melanized pycnidia and is genetically tractable through Agrobacterium tumefaciens-mediated transformation (ATMT), which can result in insertional mutants. One T-DNA-tagged mutant, P1A17 produced conidia in non-melanized (hyaline) pycnidia. Pycnidial melanization recovered if the mutant was supplemented with melanin precursors or allowed to age. DNA sequences flanking the insertion did not predict any disrupted open reading frames (ORF) unless a Coccidioides prediction algorithm was used. Pmhyp gene was expressed in the wild type, but not the mutant, and has not been annotated in any genomes, to date. Expression of two conserved genes flanking the T-DNA disrupted Pmhyp was unchanged from the wild type. Knockout of Pmhyp strain displayed same cultural phenotype (non-melanized pycnidia). Complementation of Pmhyp strains with wild type PmHYP partially recovered pycnidial melanization. Both knockout and complementation transformants were confirmed using RT-PCR and southern blot analysis. Taken together, PmHYP appears to be a novel regulator of pycnidium specific melanization.

FlbA-Regulated Gene rpnR Is Involved in Stress Resistance and Impacts Protein Secretion when Aspergillus niger Is Grown on Xylose

Proteins are secreted throughout the mycelium of Aspergillus niger except for the sporulating zone. A link between sporulation and repression of protein secretion was underlined by the finding that inactivation of the sporulation gene flbA results in mycelial colonies that secrete proteins throughout the colony. However, ΔflbA strain hyphae also lyse and have thinner cell walls. This pleiotropic phenotype is associated with differential expression of 36 predicted transcription factor genes, one of which, rpnR, was inactivated in this study. Sporulation, biomass, and secretome complexity were not affected in the ΔrpnR deletion strain of the fungus. In contrast, ribosomal subunit expression and protein secretion into the medium were reduced when A. niger was grown on xylose. Moreover, the ΔrpnR strain showed decreased resistance to H₂O₂ and the proteotoxic stress-inducing agent dithiothreitol. Taking the data together, RpnR is involved in proteotoxic stress resistance and impacts protein secretion when A. niger is grown on xylose.IMPORTANCE Aspergillus niger secretes a large amount and diversity of industrially relevant enzymes into the culture medium. This makes the fungus a widely used industrial cell factory. For instance, carbohydrate-active enzymes of A. niger are used in biofuel production from lignocellulosic feedstock. These enzymes represent a major cost factor in this process. Higher production yields could substantially reduce these costs and therefore contribute to a more sustainable economy and less dependence on fossil fuels. Enzyme secretion is inhibited in A. niger by asexual reproduction. The sporulation protein FlbA is involved in this process by impacting the expression of 36 predicted transcription factor genes. Here, we show that one of these predicted transcriptional regulators, RpnR, regulates protein secretion and proteotoxic stress resistance. The gene is thus an interesting target to improve enzyme production in A. niger.

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

Cis-Regulatory Divergence in Gene Expression between Two Thermally Divergent Yeast Species

Gene regulation is a ubiquitous mechanism by which organisms respond to their environment. While organisms are often found to be adapted to the environments they experience, the role of gene regulation in environmental adaptation is not often known. In this study, we examine divergence in cis-regulatory effects between two Saccharomycesspecies, S. cerevisiaeand S. uvarum, that have substantially diverged in their thermal growth profile. We measured allele specific expression (ASE) in the species’ hybrid at three temperatures, the highest of which is lethal to S. uvarumbut not the hybrid or S. cerevisiae. We find that S. uvarumalleles can be expressed at the same level as S. cerevisiaealleles at high temperature and most cis-acting differences in gene expression are not dependent on temperature. While a small set of 136 genes show temperature-dependent ASE, we find no indication that signatures of directional cis-regulatory evolution are associated with temperature. Within promoter regions we find binding sites enriched upstream of temperature responsive genes, but only weak correlations between binding site and expression divergence. Our results indicate that temperature divergence between S. cerevisiaeand S. uvarumhas not caused widespread divergence in cis-regulatory activity, but point to a small subset of genes where the species’ alleles show differences in magnitude or opposite responses to temperature. The difficulty of explaining divergence in cis-regulatory sequences with models of transcription factor binding sites and nucleosome positioning highlights the importance of identifying mutations that underlie cis-regulatory divergence between species.

Comparative Transcriptomics Highlights New Features of the Iron Starvation Response in the Human Pathogen Candida glabrata

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

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

Cancer cell line specific co-factors modulate the FOXM1 cistrome

ChIP-seq has been commonly applied to identify genomic occupation of transcription factors (TFs) in a context-specific manner. It is generally assumed that a TF should have similar binding patterns in cells from the same or closely related tissues. Surprisingly, this assumption has not been carefully examined. To this end, we systematically compared the genomic binding of the cell cycle regulator FOXM1 in eight cell lines from seven different human tissues at binding signal, peaks and target genes levels. We found that FOXM1 binding in ER-positive breast cancer cell line MCF-7 are distinct comparing to those in not only other non-breast cell lines, but also MDA-MB-231, ER-negative breast cancer cell line. However, binding sites in MDA-MB-231 and non-breast cell lines were highly consistent. The recruitment of estrogen receptor alpha (ERα) caused the unique FOXM1 binding patterns in MCF-7. Moreover, the activity of FOXM1 in MCF-7 reflects the regulatory functions of ERα, while in MDA-MB-231 and non-breast cell lines, FOXM1 activities regulate cell proliferation. Our results suggest that tissue similarity, in some specific contexts, does not hold precedence over TF-cofactors interactions in determining transcriptional states and that the genomic binding of a TF can be dramatically affected by a particular co-factor under certain conditions.

Functional regulatory evolution outside of the minimal even-skipped stripe 2 enhancer

Transcriptional enhancers are regions of DNA that drive precise patterns of gene expression. Although many studies have elucidated how individual enhancers can evolve, most of this work has focused on what are called 'minimal' enhancers, the smallest DNA regions that drive expression that approximates an aspect of native gene expression. Here, we explore how the Drosophila erecta even-skipped (eve) locus has evolved by testing its activity in the divergent D. melanogaster genome. We found, as has been reported previously, that the D. erecta eve stripe 2 enhancer (eveS2) fails to drive appreciable expression in D. melanogaster However, we found that a large transgene carrying the entire D. erecta eve locus drives normal eve expression, including in stripe 2. We performed a functional dissection of the region upstream of the D. erecta eveS2 region and found multiple Zelda motifs that are required for normal expression. Our results illustrate how sequences outside of minimal enhancer regions can evolve functionally through mechanisms other than changes in transcription factor-binding sites that drive patterning.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

Recurrent rewiring and emergence of RNA regulatory networks

Alterations in regulatory networks contribute to evolutionary change. Transcriptional networks are reconfigured by changes in the binding specificity of transcription factors and their cognate sites. The evolution of RNA-protein regulatory networks is far less understood. The PUF (Pumilio and FBF) family of RNA regulatory proteins controls the translation, stability, and movements of hundreds of mRNAs in a single species. We probe the evolution of PUF-RNA networks by direct identification of the mRNAs bound to PUF proteins in budding and filamentous fungi and by computational analyses of orthologous RNAs from 62 fungal species. Our findings reveal that PUF proteins gain and lose mRNAs with related and emergent biological functions during evolution. We demonstrate at least two independent rewiring events for PUF3 orthologs, independent but convergent evolution of PUF4/5 binding specificity and the rewiring of the PUF4/5 regulons in different fungal lineages. These findings demonstrate plasticity in RNA regulatory networks and suggest ways in which their rewiring occurs.

Evolution of Brain Active Gene Promoters in Human Lineage Towards the Increased Plasticity of Gene Regulation

Adaptability to a variety of environmental conditions is a prominent feature of Homo sapiens. We hypothesize that this feature can be explained by evolutionary changes in gene promoters active in the brain prefrontal cortex leading to a more flexible gene regulation network. The genotype-dependent range of gene expression can be broader in humans than in other higher primates. Thus, we searched for specific signatures of evolutionary changes in promoter architectures of multiple hominid genes, including the genes active in human cortical neurons that may indicate an increase of variability of gene expression rather than just changes in the level of expression, such as downregulation or upregulation of the genes. We performed a whole-genome search for genetic-based alterations that may impact gene regulation "flexibility" in a process of hominids evolution, such as (i) CpG dinucleotide content, (ii) predicted nucleosome-DNA dissociation constant, and (iii) predicted affinities for TATA-binding protein (TBP) in gene promoters. We tested all putative promoter regions across the human genome and especially gene promoters in active chromatin state in neurons of prefrontal cortex, the brain region critical for abstract thinking and social and behavioral adaptation. Our data imply that the origin of modern man has been associated with an increase of flexibility of promoter-driven gene regulation in brain. In contrast, after splitting from the ancestral lineages of H. sapiens, the evolution of ape species is characterized by reduced flexibility of gene promoter functioning, underlying reduced variability of the gene expression.

Natural gene expression variation studies in yeast

The rise of sequence information across different yeast species and strains is driving an increasing number of studies in the emerging field of genomics to associate polymorphic variants, mRNA abundance and phenotypic differences between individuals. Here, we gathered evidence from recent studies covering several layers that define the genotype-phenotype gap, such as mRNA abundance, allele-specific expression and translation efficiency to demonstrate how genetic variants co-evolve and define an individual's genome. Moreover, we exposed several antecedents where inter- and intra-specific studies led to opposite conclusions, probably owing to genetic divergence. Future studies in this area will benefit from the access to a massive array of well-annotated genomes and new sequencing technologies, which will allow the fine breakdown of the complex layers that delineate the genotype-phenotype map. Copyright © 2016 John Wiley & Sons, Ltd.

Evolution of transcriptional networks in yeast: alternative teams of transcriptional factors for different species

Background

The diversity in eukaryotic life reflects a diversity in regulatory pathways. Nocedal and Johnson argue that the rewiring of gene regulatory networks is a major force for the diversity of life, that changes in regulation can create new species.

Results

We have created a method (based on our new “ping-pong algorithm) for detecting more complicated rewirings, where several transcription factors can substitute for one or more transcription factors in the regulation of a family of co-regulated genes. An example is illustrative. A rewiring has been reported by Hogues et al. that RAP1 in Saccharomyces cerevisiae substitutes for TBF1/CBF1 in Candida albicans for ribosomal RP genes. There one transcription factor substitutes for another on some collection of genes. Such a substitution is referred to as a “rewiring”. We agree with this finding of rewiring as far as it goes but the situation is more complicated. Many transcription factors can regulate a gene and our algorithm finds that in this example a “team” (or collection) of three transcription factors including RAP1 substitutes for TBF1 for 19 genes. The switch occurs for a branch of the phylogenetic tree containing 10 species (including Saccharomyces cerevisiae), while the remaining 13 species (Candida albicans) are regulated by TBF1.

Conclusions

To gain insight into more general evolutionary mechanisms, we have created a mathematical algorithm that finds such general switching events and we prove that it converges. Of course any such computational discovery should be validated in the biological tests. For each branch of the phylogenetic tree and each gene module, our algorithm finds a sub-group of co-regulated genes and a team of transcription factors that substitutes for another team of transcription factors. In most cases the signal will be small but in some cases we find a strong signal of switching. We report our findings for 23 Ascomycota fungi species.

Electronic supplementary material

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

Evolution of New cis-Regulatory Motifs Required for Cell-Specific Gene Expression in Caenorhabditis

Patterning of C. elegans vulval cell fates relies on inductive signaling. In this induction event, a single cell, the gonadal anchor cell, secretes LIN-3/EGF and induces three out of six competent precursor cells to acquire a vulval fate. We previously showed that this developmental system is robust to a four-fold variation in lin-3/EGF genetic dose. Here using single-molecule FISH, we find that the mean level of expression of lin-3 in the anchor cell is remarkably conserved. No change in lin-3 expression level could be detected among C. elegans wild isolates and only a low level of change—less than 30%—in the Caenorhabditis genus and in Oscheius tipulae. In C. elegans, lin-3 expression in the anchor cell is known to require three transcription factor binding sites, specifically two E-boxes and a nuclear-hormone-receptor (NHR) binding site. Mutation of any of these three elements in C. elegans results in a dramatic decrease in lin-3 expression. Yet only a single E-box is found in the Drosophilae supergroup of Caenorhabditis species, including C. angaria, while the NHR-binding site likely only evolved at the base of the Elegans group. We find that a transgene from C. angaria bearing a single E-box is sufficient for normal expression in C. elegans. Even a short 58 bp cis-regulatory fragment from C. angaria with this single E-box is able to replace the three transcription factor binding sites at the endogenous C. elegans lin-3 locus, resulting in the wild-type expression level. Thus, regulatory evolution occurring in cis within a 58 bp lin-3 fragment, results in a strict requirement for the NHR binding site and a second E-box in C. elegans. This single-cell, single-molecule, quantitative and functional evo-devo study demonstrates that conserved expression levels can hide extensive change in cis-regulatory site requirements and highlights the evolution of new cis-regulatory elements required for cell-specific gene expression.

Diversification of mechanisms regulating gene expression of key developmental factors is a major force in the evolution of development. However, in the past, comparisons of gene expression across different species have often been qualitative (i.e. ‘expression is on versus off’ in a certain cell) without precise quantification. New experimental methods now allow us to quantitatively compare the expression of gene homologs across species, with single cell resolution. Moreover, the development of genome editing tools enables the dissection of regulatory DNA sequences that drive gene expression. We use here a well-established “textbook” example of animal organogenesis in the microscopic nematode, Caenorhabditis elegans, focusing on the expression of lin-3, coding for the main inducer of the vulva, in a single cell called the anchor cell. We find that the lin-3 expression level is remarkably conserved, with 20–25 messenger RNAs per anchor cell, in species that are molecularly as distant as fish and mammals. This conservation occurs despite substantial changes and compensation in the regulatory elements required for cell-specific gene expression.

Compartmentalized gene regulatory network of the pathogenic fungus Fusarium graminearum

Head blight caused by Fusarium graminearum threatens world-wide wheat production, resulting in both yield loss and mycotoxin contamination. We reconstructed the global F. graminearum gene regulatory network (GRN) from a large collection of transcriptomic data using Bayesian network inference, a machine-learning algorithm. This GRN reveals connectivity between key regulators and their target genes. Focusing on key regulators, this network contains eight distinct but interwoven modules. Enriched for unique functions, such as cell cycle, DNA replication, transcription, translation and stress responses, each module exhibits distinct expression profiles. Evolutionarily, the F. graminearum genome can be divided into core regions shared with closely related species and variable regions harboring genes that are unique to F. graminearum and perform species-specific functions. Interestingly, the inferred top regulators regulate genes that are significantly enriched from the same genomic regions (P < 0.05), revealing a compartmentalized network structure that may reflect network rewiring related to specific adaptation of this plant pathogen. This first-ever reconstructed filamentous fungal GRN primes our understanding of pathogenicity at the systems biology level and provides enticing prospects for novel disease control strategies involving the targeting of master regulators in pathogens. The program can be used to construct GRNs of other plant pathogens.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi

Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

DOI:http://dx.doi.org/10.7554/eLife.09492.001

Living cells grow and divide with remarkable precision to ensure that their genetic material is faithfully duplicated and distributed equally to the newly formed daughter cells. This precision is achieved through a series of steps known as the cell cycle. The cell cycle is ancient and conserved across all Eukaryotes, including plants, animals and fungi. However, some of the core proteins present in animals and fungi are unrelated. This raises the question as to how a drastic change could have occurred and been tolerated over evolution.

In animals and plants, a protein called E2F controls the expression of genes that are needed to begin the cell cycle. In most fungi, an equivalent protein called SBF performs the same role as E2F, but the two proteins are very different and do not appear to share a common ancestor. This is unexpected given that fungi and animals are more closely related to one another than either is to plants.

Medina et al. searched the genomes of many animals, fungi, plants, algae, and their closest relatives for genes that encoded proteins like E2F and SBF. SBF-like proteins were only found in fungi, yet some fungal groups had cell cycle regulators like those found in animals. Zoosporic fungi, which diverged early from the fungal ancestor, had both SBF- and E2F-like proteins, while many fungi later lost E2F during evolution.

So how did fungi acquire SBF? Medina et al. observed that part of the SBF protein is similar to proteins found in many viruses. The broad distribution of these viral SBF-like proteins suggests that they arose first in viruses, and a fungal ancestor acquired one such protein during a viral infection. As SBF and E2F bind similar DNA sequences, Medina et al. hypothesized that this viral SBF hijacked control of the cell cycle in the fungal ancestor by controlling expression of genes that were originally controlled only by E2F. In support of this idea, experiments showed that many E2F binding sites in modern genes are also SBF binding sites, and that E2F sites can substitute for SBF sites in SBF-controlled genes. Future experiments in zoosporic fungi, which have animal-like and fungal-like features, would provide a glimpse of how a fungal ancestor may have used both SBF and E2F. These experiments may also reveal why most fungi have retained the newer SBF but lost the ancestral and widely conserved E2F protein.

DOI:http://dx.doi.org/10.7554/eLife.09492.002

The Power of Natural Variation for Model Organism Biology

Genetic background effects have long been recognized and, in some cases studied, but they are often viewed as a nuisance by molecular biologists. We suggest that genetic variation currently represents a critical frontier for molecular studies. Human genetics has seen a surge of interest in genetic variation and its contributions to disease, but insights into disease mechanisms are difficult since information about gene function is lacking. By contrast, model organism genetics has excelled at revealing molecular mechanisms of cellular processes, but often de-emphasizes genetic variation and its functional consequences. We argue that model organism biology would benefit from incorporating natural variation, both to capture how well laboratory lines exemplify the species they represent and to inform on molecular processes and their variability. Such a synthesis would also greatly expand the relevance of model systems for studies of complex trait variation, including disease.

Divergent MLS1 Promoters Lie on a Fitness Plateau for Gene Expression

Qualitative patterns of gene activation and repression are often conserved despite an abundance of quantitative variation in expression levels within and between species. A major challenge to interpreting patterns of expression divergence is knowing which changes in gene expression affect fitness. To characterize the fitness effects of gene expression divergence, we placed orthologous promoters from eight yeast species upstream of malate synthase (MLS1) in Saccharomyces cerevisiae. As expected, we found these promoters varied in their expression level under activated and repressed conditions as well as in their dynamic response following loss of glucose repression. Despite these differences, only a single promoter driving near basal levels of expression caused a detectable loss of fitness. We conclude that the MLS1 promoter lies on a fitness plateau whereby even large changes in gene expression can be tolerated without a substantial loss of fitness.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Reprogramming of a gene’s expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100–500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

A map of the evolutionary history of Puf proteins and their RNA targets shows that reprogramming of global gene expression programs via adaptive mutations that affect protein-RNA interactions is an important source of biological diversity.

We set out to trace the evolutionary history of an RNA binding protein and how its interactions with targets change over evolution. Identifying this natural history is a step toward understanding the critical differences between organisms and how gene expression programs are rewired during evolution. Using bioinformatics and experimental approaches, we broadly surveyed the evolution of binding targets of a particular family of RNA binding proteins—the Puf proteins, whose protein sequences and target RNA sequences are relatively well-characterized—across 99 eukaryotic species. We found five groups of species in which targets have been conserved for at least 100 million years and then took advantage of genome sequences from a large number of fungal species to deeply investigate the conservation and changes in Puf proteins and their RNA targets. Our analyses identified multiple and extensive reconfigurations during the natural history of fungi and suggest that RNA binding proteins and their RNA targets are profoundly involved in evolutionary reprogramming of gene expression and help define distinct programs unique to each organism. Continuing to uncover the natural history of RNA binding proteins and their interactions will provide a unique window into the gene expression programs of present day species and point to new ways to engineer gene expression programs.

Comparative analysis of gene regulatory networks: from network reconstruction to evolution

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

Metabolic regulation in model ascomycetes--adjusting similar genomes to different lifestyles

The related yeasts Saccharomyces cerevisiae and Candida albicans have similar genomes but very different lifestyles. These fungi have modified transcriptional and post-translational regulatory processes to adapt their similar genomes to the distinct biological requirements of the two yeasts. We review recent findings comparing the differences between these species, highlighting how they have achieved specialized metabolic capacities tailored to their lifestyles despite sharing similar genomes. Studying this transcriptional and post-transcriptional rewiring may improve our ability to interpret phenotype from genotype.

Investigation of protein secretion and secretion stress in Ashbya gossypii

Background

Ashbya gossypii is a filamentous Saccharomycete used for the industrial production of riboflavin that has been recently explored as a host system for recombinant protein production. To gain insight into the protein secretory pathway of this biotechnologically relevant fungus, we undertook genome-wide analyses to explore its secretome and its transcriptional responses to protein secretion stress.

Results

A computational pipeline was used to predict the inventory of proteins putatively secreted by A. gossypii via the general secretory pathway. The proteins actually secreted by this fungus into the supernatants of submerged cultures in minimal and rich medium were mapped by two-dimensional gel electrophoresis, revealing that most of the A. gossypii secreted proteins have an isoelectric point between 4 and 6, and a molecular mass above 25 kDa. These analyses together indicated that 1-4% of A. gossypii proteins are likely to be secreted, of which less than 33% are putative hydrolases. Furthermore, transcriptomic analyses carried out in A. gossypii cells under recombinant protein secretion conditions and dithiothreitol-induced secretion stress unexpectedly revealed that a conventional unfolded protein response (UPR) was not activated in any of the conditions, as the expression levels of several well-known UPR target genes (e.g. IRE1, KAR2, HAC1 and PDI1 homologs) remained unaffected. However, several other genes involved in protein unfolding, endoplasmatic reticulum-associated degradation, proteolysis, vesicle trafficking, vacuolar protein sorting, secretion and mRNA degradation were up-regulated by dithiothreitol-induced secretion stress. Conversely, the transcription of several genes encoding secretory proteins, such as components of the glycosylation pathway, was severely repressed by dithiothreitol

Conclusions

This study provides the first insights into the secretion stress response of A. gossypii, as well as a basic understanding of its protein secretion potential, which is more similar to that of yeast than to that of other filamentous fungi. Contrary to what has been widely described for yeast and fungi, a conventional UPR was not observed in A. gossypii, but alternative protein quality control mechanisms enabled it to cope with secretion stress. These data will help provide strategies for improving heterologous protein secretion in A. gossypii.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1137) contains supplementary material, which is available to authorized users.

Promoter architectures in the yeast ribosomal expression program

Ribosome biogenesis begins with the orchestrated expression of hundreds of genes, including the three large classes of ribosomal protein, ribosome biogenesis and snoRNA genes. Current knowledge about the corresponding promoters suggests the existence of novel class-specific transcriptional strategies and crosstalk between telomere length and cell growth control.

Inference of expanded Lrp-like feast/famine transcription factor targets in a non-model organism using protein structure-based prediction

Widespread microbial genome sequencing presents an opportunity to understand the gene regulatory networks of non-model organisms. This requires knowledge of the binding sites for transcription factors whose DNA-binding properties are unknown or difficult to infer. We adapted a protein structure-based method to predict the specificities and putative regulons of homologous transcription factors across diverse species. As a proof-of-concept we predicted the specificities and transcriptional target genes of divergent archaeal feast/famine regulatory proteins, several of which are encoded in the genome of Halobacterium salinarum. This was validated by comparison to experimentally determined specificities for transcription factors in distantly related extremophiles, chromatin immunoprecipitation experiments, and cis-regulatory sequence conservation across eighteen related species of halobacteria. Through this analysis we were able to infer that Halobacterium salinarum employs a divergent local trans-regulatory strategy to regulate genes (carA and carB) involved in arginine and pyrimidine metabolism, whereas Escherichia coli employs an operon. The prediction of gene regulatory binding sites using structure-based methods is useful for the inference of gene regulatory relationships in new species that are otherwise difficult to infer.