Genome adaptation to chemical stress: clues from comparative transcriptomics in Saccharomyces cerevisiae and Candida glabrata

ifDEEPre: large protein language-based deep learning enables interpretable and fast predictions of enzyme commission numbers

Accurate understanding of the biological functions of enzymes is vital for various tasks in both pathologies and industrial biotechnology. However, the existing methods are usually not fast enough and lack explanations on the prediction results, which severely limits their real-world applications. Following our previous work, DEEPre, we propose a new interpretable and fast version (ifDEEPre) by designing novel self-guided attention and incorporating biological knowledge learned via large protein language models to accurately predict the commission numbers of enzymes and confirm their functions. Novel self-guided attention is designed to optimize the unique contributions of representations, automatically detecting key protein motifs to provide meaningful interpretations. Representations learned from raw protein sequences are strictly screened to improve the running speed of the framework, 50 times faster than DEEPre while requiring 12.89 times smaller storage space. Large language modules are incorporated to learn physical properties from hundreds of millions of proteins, extending biological knowledge of the whole network. Extensive experiments indicate that ifDEEPre outperforms all the current methods, achieving more than 14.22% larger F1-score on the NEW dataset. Furthermore, the trained ifDEEPre models accurately capture multi-level protein biological patterns and infer evolutionary trends of enzymes by taking only raw sequences without label information. Meanwhile, ifDEEPre predicts the evolutionary relationships between different yeast sub-species, which are highly consistent with the ground truth. Case studies indicate that ifDEEPre can detect key amino acid motifs, which have important implications for designing novel enzymes. A web server running ifDEEPre is available at https://proj.cse.cuhk.edu.hk/aihlab/ifdeepre/ to provide convenient services to the public. Meanwhile, ifDEEPre is freely available on GitHub at https://github.com/ml4bio/ifDEEPre/.

Is There a Relationship Between Mating and Pathogenesis in Two Human Fungal Pathogens, Candida albicans and Candida glabrata?

Purpose of Review

Human fungal pathogens are rapidly increasing in incidence and readily able to evade the host immune responses. Our ability to study the genetic behind this has been limited due to the apparent lack of a sexual cycle and forward genetic tools. In this review, we discuss the evolution of mating, meiosis, and pathogenesis and if these processes are advantageous to pathogens.

Recent Findings

This review summarises what is currently known about the sexual cycles of two important human fungal pathogens, Candida albicans and Candida glabrata. This includes the identification of parasexual cycle in C. albicans and the observed low levels of recombination in C. glabrata populations.

Summary

In this review, we present what is currently known about the mating types and mating/sexual cycles of two clinically important human fungal pathogens, Candida albicans and Candida glabrata. We discuss the evolution of meiosis using the knowledge that has been amassed from the decades of studying Saccharomyces cerevisiae and how this can be applied to fungal pathogens. We further discuss how the evolution of pathogenesis has played a role in influencing mating processes in human fungal pathogens and compare sexual cycles between C. albicans and C. glabrata, highlighting knowledge gaps and suggesting how these two fungi have evolved distinct mating niches to allow the development of disease in a human host.

Transcriptional Response of Candida auris to the Mrr1 Inducers Methylglyoxal and Benomyl

Candida auris is an urgent threat to human health due to its rapid spread in health care settings and its repeated development of multidrug resistance. Diseases that increase risk for C. auris infection, such as diabetes, kidney failure, or immunocompromising conditions, are associated with elevated levels of methylglyoxal (MG), a reactive dicarbonyl compound derived from several metabolic processes. In other Candida species, expression of MG reductase enzymes that catabolize and detoxify MG are controlled by Mrr1, a multidrug resistance-associated transcription factor, and MG induces Mrr1 activity. Here, we used transcriptomics and genetic assays to determine that C. aurisMRR1a contributes to MG resistance, and that the main Mrr1a targets are an MG reductase and MDR1, which encodes a drug efflux protein. The C. auris Mrr1a regulon is smaller than Mrr1 regulons described in other species. In addition to MG, benomyl (BEN), a known Mrr1 stimulus, induces C. auris Mrr1 activity, and characterization of the MRR1a-dependent and -independent transcriptional responses revealed substantial overlap in genes that were differentially expressed in response to each compound. Additionally, we found that an MRR1 allele specific to one C. auris phylogenetic clade, clade III, encodes a hyperactive Mrr1 variant, and this activity correlated with higher MG resistance. C. aurisMRR1a alleles were functional in Candida lusitaniae and were inducible by BEN, but not by MG, suggesting that the two Mrr1 inducers act via different mechanisms. Together, the data presented in this work contribute to the understanding of Mrr1 activity and MG resistance in C. auris.

IMPORTANCECandida auris is a fungal pathogen that has spread since its identification in 2009 and is of concern due to its high incidence of resistance against multiple classes of antifungal drugs. In other Candida species, the transcription factor Mrr1 plays a major role in resistance against azole antifungals and other toxins. More recently, Mrr1 has been recognized to contribute to resistance to methylglyoxal (MG), a toxic metabolic product that is often elevated in different disease states. MG can activate Mrr1 and its induction of Mdr1 which can protect against diverse challenges. The significance of this work lies in showing that MG is also an inducer of Mrr1 in C. auris, and that one of the major pathogenic C. auris lineages has an activating Mrr1 mutation that confers protection against MG.

Erg25 Controls Host-Cholesterol Uptake Mediated by Aus1p-Associated Sterol-Rich Membrane Domains in Candida glabrata

The uptake of cholesterol from the host is closely linked to the proliferation of pathogenic fungi and protozoa during infection. For some pathogenic fungi, cholesterol uptake is an important strategy for decreasing susceptibility to antifungals that inhibit ergosterol biosynthesis. In this study, we show that Candida glabrata ERG25, which encodes an enzyme that demethylates 4,4-dimethylzymosterol, is required for cholesterol uptake from host serum. Based on the screening of C. glabrata conditional knockdown mutants for each gene involved in ergosterol biosynthesis, ERG25 knockdown was found to decrease lethality of infected mice. ERG25 knockdown impairs the plasma membrane localization of the sterol importer Aus1p, suggesting that the accumulated 4,4-dimethylzymosterol destabilizes the lipid domain with which Aus1p functionally associates. ERG25 knockdown further influences the structure of the membrane compartment of Can1p (MCC)/eisosomes (ergosterol-rich lipid domains), but not the localization of the membrane proteins Pma1p and Hxt1p, which localize to sterol-poor domains. In the sterol-rich lipid domain, Aus1p-contining domain was mostly independent of MCC/eisosomes, and the nature of these domains was also different: Ausp1-contining domain was a dynamic network-like domain, whereas the MCC/eisosomes was a static dot-like domain. However, deletion of MCC/eisosomes was observed to influence the localization of Aus1p after Aus1p was transported from the endoplasmic reticulum (ER) through the Golgi apparatus to the plasma membrane. These findings suggest that ERG25 plays a key role in stabilizing sterol-rich lipid domains, constituting a promising candidate target for antifungal therapy.

Yeast AP-1 like transcription factors (Yap) and stress response: a current overview

Yeast adaptation to stress has been extensively studied. It involves large reprogramming of genome expression operated by many, more or less specific, transcription factors. Here, we review our current knowledge on the function of the eight Yap transcription factors (Yap1 to Yap8) in Saccharomyces cerevisiae, which were shown to be involved in various stress responses. More precisely, Yap1 is activated under oxidative stress, Yap2/Cad1 under cadmium, Yap4/Cin5 and Yap6 under osmotic shock, Yap5 under iron overload and Yap8/Arr1 by arsenic compounds. Yap3 and Yap7 seem to be involved in hydroquinone and nitrosative stresses, respectively. The data presented in this article illustrate how much knowledge on the function of these Yap transcription factors is advanced. The evolution of the Yap family and its roles in various pathogenic and non-pathogenic fungal species is discussed in the last section.

Pixel: a content management platform for quantitative omics data

BACKGROUND

In biology, high-throughput experimental technologies, also referred as "omics" technologies, are increasingly used in research laboratories. Several thousands of gene expression measurements can be obtained in a single experiment. Researchers are routinely facing the challenge to annotate, store, explore and mine all the biological information they have at their disposal. We present here the Pixel web application (Pixel Web App), an original content management platform to help people involved in a multi-omics biological project.

METHODS

The Pixel Web App is built with open source technologies and hosted on the collaborative development platform GitHub (https://github.com/Candihub/pixel). It is written in Python using the Django framework and stores all the data in a PostgreSQL database. It is developed in the open and licensed under the BSD 3-clause license. The Pixel Web App is also heavily tested with both unit and functional tests, a strong code coverage and continuous integration provided by CircleCI. To ease the development and the deployment of the Pixel Web App, Docker and Docker Compose are used to bundle the application as well as its dependencies.

RESULTS

The Pixel Web App offers researchers an intuitive way to annotate, store, explore and mine their multi-omics results. It can be installed on a personal computer or on a server to fit the needs of many users. In addition, anyone can enhance the application to better suit their needs, either by contributing directly on GitHub (encouraged) or by extending Pixel on their own. The Pixel Web App does not provide any computational programs to analyze the data. Still, it helps to rapidly explore and mine existing results and holds a strategic position in the management of research data.

Microevolution of the pathogenic yeasts Candida albicans and Candida glabrata during antifungal therapy and host infection

Infections by the pathogenic yeasts Candida albicans and Candida glabrata are among the most common fungal diseases. The success of these species as human pathogens is contingent on their ability to resist antifungal therapy and thrive within the human host. C. glabrata is especially resilient to azole antifungal treatment, while C. albicans is best known for its wide array of virulence features. The core mechanisms that underlie antifungal resistance and virulence in these pathogens has been continuously addressed, but the investigation on how such mechanisms evolve according to each environment is scarcer. This review aims to explore current knowledge on micro-evolution experiments to several treatment and host-associated conditions in C. albicans and C. glabrata. The analysis of adaptation strategies that evolve over time will allow to better understand the mechanisms by which Candida species are able to achieve stable phenotypes in real-life scenarios, which are the ones that should constitute the most interesting drug targets.

Identifying novel genetic determinants for oxidative stress tolerance in Candida glabrata via adaptive laboratory evolution

Candida glabrata (C glabrata) is an important yeast of industrial and medical significance. Resistance to oxidative stress is an important trait affecting its robustness as a production host or virulence as a pathogenic agent, but current understanding of resistance mechanisms is still limited in this fungus. In this study, we rapidly evolved C glabrata population to adapt to oxidative challenge (from 80mM to 350mM of H₂ O₂ ) through short-term adaptive laboratory evolution. Adaptive mutants were isolated from evolved populations and subjected to phenotypic and omics analyses to identify potential mechanisms of tolerance to H₂ O₂ . Phenotypic characterizations revealed faster detoxification of H₂ O₂ and ability to initiate growth at a higher concentration of the oxidant in the isolated adaptive mutants compared with the wild type. Genome resequencing and genome-wide transcriptome analysis revealed multiple genetic determinants (eg, CAGL0E01243g, CAGL0F06831g, and CAGL0C00385g) that potentially contribute to enhanced H₂ O₂ resistance. Subsequent experimental verification confirmed that CgCth2 (CAGL0E01243g) and CgMga2 (CAGL0F06831g) are important in C glabrata tolerance to oxidative stress. Transcriptome profiling of adaptive mutants and bioinformatic analysis suggest that NADPH regeneration, modulation of membrane composition, cell wall remodeling, and/or global regulatory changes are involved in C glabrata tolerance to H₂ O₂ .

The Concept of Resilience in Ecological Risk Assessment: Scientific and Regulatory Issues

Resilience represents one of the key components of the vulnerability of ecological systems and may refer to different levels of biological organization, from populations to the biosphere. A short description is given on the concept of resilience applied to the levels that are directly involved in ecological risk assessment (ERA): populations, communities, and ecosystems. The opportunities and challenges for measuring and quantifying resilience are discussed. Finally, some suggestions for introducing the resilience concept in regulatory ERA are proposed. Integr Environ Assess Manag 2018;14:581-585. © 2018 SETAC.

The CgHaa1-Regulon Mediates Response and Tolerance to Acetic Acid Stress in the Human Pathogen Candida glabrata

To thrive in the acidic vaginal tract, Candida glabrata has to cope with high concentrations of acetic acid. The mechanisms underlying C. glabrata tolerance to acetic acid at low pH remain largely uncharacterized. In this work, the essential role of the CgHaa1 transcription factor (encoded by ORF CAGL0L09339g) in the response and tolerance of C. glabrata to acetic acid is demonstrated. Transcriptomic analysis showed that CgHaa1 regulates, directly or indirectly, the expression of about 75% of the genes activated under acetic acid stress. CgHaa1-activated targets are involved in multiple physiological functions including membrane transport, metabolism of carbohydrates and amino acids, regulation of the activity of the plasma membrane H⁺-ATPase, and adhesion. Under acetic acid stress, CgHaa1 increased the activity and the expression of the CgPma1 proton pump and contributed to increased colonization of vaginal epithelial cells by C. glabrata. CgHAA1, and two identified CgHaa1-activated targets, CgTPO3 and CgHSP30, are herein demonstrated to be determinants of C. glabrata tolerance to acetic acid. The protective effect of CgTpo3 and of CgHaa1 was linked to a role of these proteins in reducing the accumulation of acetic acid inside C. glabrata cells. In response to acetic acid stress, marked differences were found in the regulons controlled by CgHaa1 and by its S. cerevisiae ScHaa1 ortholog, demonstrating a clear divergent evolution of the two regulatory networks. The results gathered in this study significantly advance the understanding of the molecular mechanisms underlying the success of C. glabrata as a vaginal colonizer.

Membrane Proteomics Analysis of the Candida glabrata Response to 5-Flucytosine: Unveiling the Role and Regulation of the Drug Efflux Transporters CgFlr1 and CgFlr2

Resistance to 5-flucytosine (5-FC), used as an antifungal drug in combination therapy, compromises its therapeutic action. In this work, the response of the human pathogen Candida glabrata to 5-FC was evaluated at the membrane proteome level, using an iTRAQ-based approach. A total of 32 proteins were found to display significant expression changes in the membrane fraction of cells upon exposure to 5-FC, 50% of which under the control of CgPdr1, the major regulator of azole drug resistance. These proteins cluster into functional groups associated to cell wall assembly, lipid metabolism, amino acid/nucleotide metabolism, ribosome components and translation machinery, mitochondrial function, glucose metabolism, and multidrug resistance transport. Given the obtained indications, the function of the drug:H+ antiporters CgFlr1 (ORF CAGL0H06017g) and CgFlr2 (ORF CAGL0H06039g) was evaluated. The expression of both proteins, localized to the plasma membrane, was found to confer flucytosine resistance. CgFlr2 further confers azole drug resistance. The deletion of CgFLR1 or CgFLR2 was seen to increase the intracellular accumulation of 5-FC, or 5-FC and clotrimazole, suggesting that these transporters play direct roles in drug extrusion. The expression of CgFLR1 and CgFLR2 was found to be controlled by the transcription factors CgPdr1 and CgYap1, major regulator of oxidative stress resistance.

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.

A Network of Paralogous Stress Response Transcription Factors in the Human Pathogen Candida glabrata

The yeast Candida glabrata has become the second cause of systemic candidemia in humans. However, relatively few genome-wide studies have been conducted in this organism and our knowledge of its transcriptional regulatory network is quite limited. In the present work, we combined genome-wide chromatin immunoprecipitation (ChIP-seq), transcriptome analyses, and DNA binding motif predictions to describe the regulatory interactions of the seven Yap (Yeast AP1) transcription factors of C. glabrata. We described a transcriptional network containing 255 regulatory interactions and 309 potential target genes. We predicted with high confidence the preferred DNA binding sites for 5 of the 7 CgYaps and showed a strong conservation of the Yap DNA binding properties between S. cerevisiae and C. glabrata. We provided reliable functional annotation for 3 of the 7 Yaps and identified for Yap1 and Yap5 a core regulon which is conserved in S. cerevisiae, C. glabrata, and C. albicans. We uncovered new roles for CgYap7 in the regulation of iron-sulfur cluster biogenesis, for CgYap1 in the regulation of heme biosynthesis and for CgYap5 in the repression of GRX4 in response to iron starvation. These transcription factors define an interconnected transcriptional network at the cross-roads between redox homeostasis, oxygen consumption, and iron metabolism.

Evolution of regulatory networks in Candida glabrata: learning to live with the human host

The opportunistic human fungal pathogen Candida glabrata is second only to C. albicans as the cause of Candida infections and yet is more closely related to Saccharomyces cerevisiae. Recent advances in functional genomics technologies and computational approaches to decipher regulatory networks, and the comparison of these networks among these and other Ascomycete species, have revealed both unique and shared strategies in adaptation to a human commensal/opportunistic pathogen lifestyle and antifungal drug resistance in C. glabrata. Recently, several C. glabrata sister species in the Nakeseomyces clade representing both human associated (commensal) and environmental isolates have had their genomes sequenced and analyzed. This has paved the way for comparative functional genomics studies to characterize the regulatory networks in these species to identify informative patterns of conservation and divergence linked to phenotypic evolution in the Nakaseomyces lineage.

Saccharomyces cerevisiae as a model in ecotoxicological studies: A post-genomics perspective

The budding yeast Saccharomyces cerevisiae represents a well-consolidated and widely used eukaryotic model, with a number of features that make it an ideal organism to carry out functional toxicological studies. Several advantages are permitted by the use of yeast cells, as the possibility to identify molecular biomarkers, unknown mechanisms of action and novel potential targets. Thanks to the evolutionary conservation, yeast can provide also useful clues allowing the prioritization of more complex analyses and toxicity predictions in higher eukaryotes. The last two decades were incredibly fruitful for yeast "omics", but referring to the analysis of the effects of pesticides on yeast much still remains to be done. Furthermore, a deeper knowledge of the effects of environmental pollutants on biotechnological processes associated with the use of yeasts is to be hoped.

Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance

In light of ongoing climate changes in wine-growing regions, the selection of drought-tolerant rootstocks is becoming a crucial factor for developing a sustainable viticulture. In this study, M4, a new rootstock genotype that shows tolerance to drought, was compared from a genomic and transcriptomic point of view with the less drought-tolerant genotype 101.14. The root and leaf transcriptome of both 101.14 and the M4 rootstock genotype was analysed, following exposure to progressive drought conditions. Multifactorial analyses indicated that stress treatment represents the main factor driving differential gene expression in roots, whereas in leaves the genotype is the prominent factor. Upon stress, M4 roots and leaves showed a higher induction of resveratrol and flavonoid biosynthetic genes, respectively. The higher expression of VvSTS genes in M4, confirmed by the accumulation of higher levels of resveratrol in M4 roots compared with 101.14, was coupled to an up-regulation of several VvWRKY transcription factors. Interestingly, VvSTS promoter analyses performed on both the resequenced genomes highlighted a significantly higher number of W-BOX elements in the tolerant genotype. It is proposed that the elevated synthesis of resveratrol in M4 roots upon water stress could enhance the plant's ability to cope with the oxidative stress usually associated with water deficit.

Stress response and expression of fluconazole resistance associated genes in the pathogenic yeast Candida glabrata deleted in the CgPDR16 gene

In yeasts, the PDR16 gene encodes a phosphatidylinositol transfer protein which belongs to the Sec14 homologue (SFH) family and localizes to lipid droplets, microsomes and at the cell periphery. The loss of its function alters the lipid droplet metabolism and plasma membrane properties, and renders yeast cells more sensitive to azole antimycotics. In this study, the entire chromosomal CgPDR16 ORF was replaced by the ScURA3 gene both in azole sensitive and azole resistant strains of Candida glabrata bearing a gain-of-function mutation in the CgPDR1 gene, and their responses to different stresses were assessed. The CgPDR16 deletion was found to sensitize the mutant strains to azole antifungals without changes in their osmo- and halotolerance. Fluconazole treated pdr16Δ mutant strains displayed a reduced expression of several genes involved in azole tolerance. The gain-of-function CgPDR1 allele as well as the cycloheximide and hydrogen peroxide treatments of cells enhanced the expression of the CgPDR16 gene. The results indicate that CgPDR16 belongs to genes whose expression is induced by chemical and oxidative stresses. The loss of its function can attenuate the expression of drug efflux pump encoding genes that might also contribute to the decreased azole tolerance in pdr16Δ mutant cells.

bPeaks: a bioinformatics tool to detect transcription factor binding sites from ChIPseq data in yeasts and other organisms with small genomes

Peak calling is a critical step in ChIPseq data analysis. Choosing the correct algorithm as well as optimized parameters for a specific biological system is an essential task. In this article, we present an original peak-calling method (bPeaks) specifically designed to detect transcription factor (TF) binding sites in small eukaryotic genomes, such as in yeasts. As TF interactions with DNA are strong and generate high binding signals, bPeaks uses simple parameters to compare the sequences (reads) obtained from the immunoprecipitation (IP) with those from the control DNA (input). Because yeasts have small genomes (<20 Mb), our program has the advantage of using ChIPseq information at the single nucleotide level and can explore, in a reasonable computational time, results obtained with different sets of parameter values. Graphical outputs and text files are provided to rapidly assess the relevance of the detected peaks. Taking advantage of the simple promoter structure in yeasts, additional functions were implemented in bPeaks to automatically assign the peaks to promoter regions and retrieve peak coordinates on the DNA sequence for further predictions of regulatory motifs, enriched in the list of peaks. Applications of the bPeaks program to three different ChIPseq datasets from Saccharomyces cerevisiae, Candida albicans and Candida glabrata are presented. Each time, bPeaks allowed us to correctly predict the DNA binding sequence of the studied TF and provided relevant lists of peaks. The bioinformatics tool bPeaks is freely distributed to academic users. Supplementary data, together with detailed tutorials, are available online: http://bpeaks.gene-networks.net.

Improvement of oxidative stress tolerance in Saccharomyces cerevisiae through global transcription machinery engineering

Excessive oxidative stress poses significant damage to yeast cells during fermentation process, and finally affects fermentation efficiency and the quality of products. In this paper, global transcription machinery engineering was employed to elicit Saccharomyces cerevisiae phenotypes of higher tolerance against oxidative stress caused by H2O2. Two strains from two plasmid-based mutagenesis libraries (Spt15 and Taf25), which exhibited significant increases in oxidative stress tolerance, were successfully isolated. At moderate H2O2 shock (≤3.5 mM), a positive correlation was found between the outperformance in cell growth of the oxidation-tolerate strains and H2O2 concentration. Several mutations were observed in the native transcription factors, which resulted in a different transcriptional profile compared with the control. Catalase and superoxide dismutase activities of the two mutants increased under H2O2 stress conditions. Fermentation experiments revealed that the mutant strain taf25-3 has a shorter lag phase compared to the control one, indicating that taf25-3 had improved adaptation ability to H2O2-induced oxidative stress and higher fermentation efficiency. Our study demonstrated that several amino acid substitutions in general transcription factors (Spt15 and Taf25) could modify the cellular oxidation defense systems and improve the anti-oxidation ability of S. cerevisiae. It could make the industrial ethanol fermentation more efficient and cost-effective by using the strain of higher stress tolerance.

Mutation of the CgPDR16 gene attenuates azole tolerance and biofilm production in pathogenic Candida glabrata

The PDR16 gene encodes the homologue of Sec14p, participating in protein secretion, regulation of lipid synthesis and turnover in vivo and acting as a phosphatidylinositol transfer protein in vitro. This gene is also involved in the regulation of multidrug resistance in Saccharomyces cerevisiae and pathogenic yeasts. Here we report the results of functional analysis of the CgPDR16 gene, whose mutation has been previously shown to enhance fluconazole sensitivity in Candida glabrata mutant cells. We have cloned the CgPDR16 gene, which was able to complement the pdr16Δ mutation in both C. glabrata and S. cerevisiae. Along with fluconazole, the pdr16Δ mutation resulted in increased susceptibility of mutant cells to several azole antifungals without changes in sensitivity to polyene antibiotics, cycloheximide, NQO, 5-fluorocytosine and oxidants inducing the intracellular formation of reactive oxygen species. The susceptibility of the pdr16Δ mutant strain to itraconazole and 5-fluorocytosine was enhanced by CTBT [7-chlorotetrazolo(5,1-c)benzo(1,2,4)triazine] inducing oxidative stress. The pdr16Δ mutation increased the accumulation of rhodamine 6G in mutant cells, decreased the level of itraconazole resistance caused by gain-of-function mutations in the CgPDR1 gene, and reduced cell surface hydrophobicity and biofilm production. These results point to the pleiotropic phenotype of the pdr16Δ mutant and support the role of the CgPDR16 gene in the control of drug susceptibility and virulence in the pathogenic C. glabrata.

Evolutionary principles of modular gene regulation in yeasts

Divergence in gene regulation can play a major role in evolution. Here, we used a phylogenetic framework to measure mRNA profiles in 15 yeast species from the phylum Ascomycota and reconstruct the evolution of their modular regulatory programs along a time course of growth on glucose over 300 million years. We found that modules have diverged proportionally to phylogenetic distance, with prominent changes in gene regulation accompanying changes in lifestyle and ploidy, especially in carbon metabolism. Paralogs have significantly contributed to regulatory divergence, typically within a very short window from their duplication. Paralogs from a whole genome duplication (WGD) event have a uniquely substantial contribution that extends over a longer span. Similar patterns occur when considering the evolution of the heat shock regulatory program measured in eight of the species, suggesting that these are general evolutionary principles.

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

The incredible diversity of living creatures belies the fact that their genes are quite similar. In the 1970s Mary-Claire King and Allan Wilson proposed that a process called gene regulation—which determines when, where and how genes are expressed as proteins—is responsible for this diversity. Four decades later, the central role of gene regulation in evolution has been confirmed in a wide range of species including bacteria, fungi, flies and mammals, although the details remain poorly understood. In recent years it has been suggested that the duplication of genes—and sometimes the duplication of whole genomes—has had a crucial influence on the part played by gene regulation in the evolution of many different species.

Ascomycota fungi are uniquely suited to the study of genetics and evolution because of their diversity—they include C. albicans, a fungus that is found in the human mouth and gut, and various species of yeast—and because many of their genomes have already been sequenced. Moreover, their genomes are relatively small, which simplifies the task of working out how it has changed over the course of evolution. It is also known that species in this branch of the tree of life diverged before and after an event in which a whole genome was duplicated.

Ascomycota fungi use glucose as a source of carbon in different ways during aerobic growth. Most, including C. albicans, are respiratory and rely on oxidative phosphorylation processes to produce energy. However, a small number—including S. cerevisiae and S. pombe, two types of yeast that are widely used as model organisms—prefer to ferment glucose, even when oxygen is available. Species that favor the latter respiro-fermentative lifestyle have evolved independently at least twice: once after the whole genome duplication event that lead to S. cerevisiae, and once when S. pombe and the other fission yeasts evolved.

Thompson et al. have measured mRNA profiles in 15 different species of yeast and reconstructed how the regulation of groups of genes (modules) have evolved over a period of more than 300 million years. They found that modules have diverged proportionally to evolutionary time, with prominent changes in gene regulation being associated with changes in lifestyle (especially changes in carbon metabolism) and a whole genome duplication event.

Gene duplication events result in gene paralogs—identical genes at different places in the genome—and these have made significant contributions to the evolution of different forms of gene regulation, especially just after the duplication event. Moreover, the paralogs produced in whole genome duplication events have resulted in bigger changes over longer periods of time. Similar patterns were observed in the regulation of the genes involved in the response to heat shock in eight of the species, which suggests that these are general evolutionary principles.

The changes in gene expression associated with the respiro-fermentative lifestyle may also have implications for our understanding of cancer: healthy cells rely on oxidative phosphorylation to produce energy whereas, similar to yeast cells, most cancerous cells rely on respiro-fermentation. Furthermore, yeast cells and cancer cells both support their rapid growth and proliferation by using glucose for biosynthesis to support cell division, although this process is not fully understood. Normal cells, on the other hand, use glucose primarily for energy and tend not to divide rapidly.

Thompson et al. found that the genes encoding enzymes in two biosynthetic pathways—one that produces the nucleotides necessary for DNA replication, and one that synthesizes glycine—are induced in respiro-fermentative yeasts but repressed in respiratory yeast cells. The fact that similar changes are observed in the same two pathways when normal cells become cancer cells suggests that these pathways have an important role in the development of cancer. The framework developed by Thompson et al. could also be used to explore the evolution of gene regulation in other species and biological processes.

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

Antifungal activity of fused Mannich ketones triggers an oxidative stress response and is Cap1-dependent in Candida albicans

We investigated the antifungal activity of fused Mannich ketone (FMK) congeners and two of their aminoalcohol derivatives. In particular, FMKs with five-membered saturated rings were shown to have minimum inhibitory concentration (MIC90s) ranging from 0.8 to 6 µg/mL toward C. albicans and the closely related C. parapsilosis and C. krusei while having reduced efficacy toward C. glabrata and almost no efficacy against Aspergillus sp. Transcript profiling of C. albicans cells exposed for 30 or 60 min to 2-(morpholinomethyl)-1-indanone, a representative FMK with a five-membered saturated ring, revealed a transcriptional response typical of oxidative stress and similar to that of a C. albicans Cap1 transcriptional activator. Consistently, C. albicans lacking the CAP1 gene was hypersensitive to this FMK, while C. albicans strains overexpressing CAP1 had decreased sensitivity to 2-(morpholinomethyl)-1-indanone. Quantitative structure-activity relationship studies revealed a correlation of antifungal potency and the energy of the lowest unoccupied molecular orbital of FMKs and unsaturated Mannich ketones thereby implicating redox cycling-mediated oxidative stress as a mechanism of action. This conclusion was further supported by the loss of antifungal activity upon conversion of representative FMKs to aminoalcohols that were unable to participate in redox cycles.

Comparative gene expression between two yeast species

Background

Comparative genomics brings insight into sequence evolution, but even more may be learned by coupling sequence analyses with experimental tests of gene function and regulation. However, the reliability of such comparisons is often limited by biased sampling of expression conditions and incomplete knowledge of gene functions across species. To address these challenges, we previously systematically generated expression profiles in Saccharomyces bayanus to maximize functional coverage as compared to an existing Saccharomyces cerevisiae data repository.

Results

In this paper, we take advantage of these two data repositories to compare patterns of ortholog expression in a wide variety of conditions. First, we developed a scalable metric for expression divergence that enabled us to detect a significant correlation between sequence and expression conservation on the global level, which previous smaller-scale expression studies failed to detect. Despite this global conservation trend, between-species gene expression neighborhoods were less well-conserved than within-species comparisons across different environmental perturbations, and approximately 4% of orthologs exhibited a significant change in co-expression partners. Furthermore, our analysis of matched perturbations collected in both species (such as diauxic shift and cell cycle synchrony) demonstrated that approximately a quarter of orthologs exhibit condition-specific expression pattern differences.

Conclusions

Taken together, these analyses provide a global view of gene expression patterns between two species, both in terms of the conditions and timing of a gene's expression as well as co-expression partners. Our results provide testable hypotheses that will direct future experiments to determine how these changes may be specified in the genome.

Milbemycins: more than efflux inhibitors for fungal pathogens

Existing antifungal agents are still confronted to activities limited to specific fungal species and to the development of resistance. Several improvements are possible either by tackling and overcoming resistance or exacerbating the activity of existing antifungal agents. In Candida glabrata, azole resistance is almost exclusively mediated by ABC transporters (including C. glabrata CDR1 [CgCDR1] and CgCDR2) via gain-of-function mutations in the transcriptional activator CgPDR1 or by mitochondrial dysfunctions. We also observed that azole resistance was correlating with increasing virulence and fitness of C. glabrata in animal models of infection. This observation motivated the re-exploitation of ABC transporter inhibitors as a possible therapeutic intervention to decrease not only the development of azole resistance but also to interfere with the virulence of C. glabrata. Milbemycins are known ABC transporter inhibitors, and here we used commercially available milbemycin A3/A4 oxim derivatives to verify this effect. As expected, the derivatives were inhibiting C. glabrata efflux with the highest activity for A3 oxim below 1 μg/ml. More surprising was that oxim derivatives had intrinsic fungicidal activity above 3.2 μg/ml, thus highlighting effects additional to the efflux inhibition. Similar values were obtained with C. albicans. Our data show that the fungicidal activity could be related to reactive oxygen species formation in these species. Transcriptional analysis performed both in C. glabrata and C. albicans exposed to A3 oxim highlighted a core of commonly regulated genes involved in stress responses, including genes involved in oxidoreductive processes, protein ubiquitination, and vesicle trafficking, as well as mitogen-activated protein kinases. However, the transcript profiles contained also species-specific signatures. Following these observations, experimental treatments of invasive infections were performed in mice treated with the commercial A3/A4 oxim preparation alone or in combination with fluconazole. Tissue burden analysis revealed that oxims on their own were able to decrease fungal burdens in both Candida species. In azole-resistant isolates, oxims acted synergistically in vivo with fluconazole to reduce fungal burden to levels of azole-susceptible isolates. In conclusion, we show here the potential of milbemycins not only as drug efflux inhibitors but also as effective fungal growth inhibitors in C. glabrata and C. albicans.

Divergence of the yeast transcription factor FZF1 affects sulfite resistance

Changes in gene expression are commonly observed during evolution. However, the phenotypic consequences of expression divergence are frequently unknown and difficult to measure. Transcriptional regulators provide a mechanism by which phenotypic divergence can occur through multiple, coordinated changes in gene expression during development or in response to environmental changes. Yet, some changes in transcriptional regulators may be constrained by their pleiotropic effects on gene expression. Here, we use a genome-wide screen for promoters that are likely to have diverged in function and identify a yeast transcription factor, FZF1, that has evolved substantial differences in its ability to confer resistance to sulfites. Chimeric alleles from four Saccharomyces species show that divergence in FZF1 activity is due to changes in both its coding and upstream noncoding sequence. Between the two closest species, noncoding changes affect the expression of FZF1, whereas coding changes affect the expression of SSU1, a sulfite efflux pump activated by FZF1. Both coding and noncoding changes also affect the expression of many other genes. Our results show how divergence in the coding and promoter region of a transcription factor alters the response to an environmental stress.

Changes in gene regulation are thought to play an important role in evolution. While variation in gene expression between species is common, it is hard to identify the phenotypic consequences of this variation since many changes in gene expression may have subtle or no phenotypic effects. In this study, we investigate changes in sulfite resistance and gene expression caused by the transcription factor, FZF1, that has evolved rapidly during the divergence of related yeast species. We find that divergence in the ability of FZF1 to confer sulfite resistance is mediated by changes in its expression as well as changes in its protein structure, both of which cause changes in the expression of other genes. Our results show how the combination of multiple changes within a transcription factor can produce substantial changes in phenotype and the expression of many genes.

Extensive divergence of yeast stress responses through transitions between induced and constitutive activation

Closely related species show a high degree of differences in gene expression, but the functional significance of these differences remains unclear. Similarly, stress responses in yeast typically involve differential expression of numerous genes, and it is unclear how many of these are functionally significant. To address these issues, we compared the expression programs of four yeast species under different growth conditions, and found that the response of these species to stress has diverged extensively. On an individual gene basis, most transcriptional responses are not conserved in any pair of species, and there are very limited common responses among all four species. We present evidence that many evolutionary changes in stress responses are compensated either (i) by the response of related genes or (ii) by changes in the basal expression levels of the genes whose responses have diverged. Thus, stress-related genes are often induced upon stress in some species but maintain high levels even in the absence of stress at other species, indicating a transition between induced and constitutive activation. In addition, ~15% of the stress responses are specific to only one of the four species, with no evidence for compensating effects or stress-related annotations, and these may reflect fortuitous regulation that is unimportant for the stress response (i.e., biological noise). Frequent compensatory changes and biological noise may explain how diverged expression responses support similar physiological responses.

The evolution of gene expression regulatory networks in yeasts

Gene regulation is a major source of phenotypic diversity between and within species. This aspect of evolution has long been addressed from the sole point of view of the genome sequence. The incredible development of transcriptomics approaches now allows one to actually study the topology and the properties of regulatory networks on an evolutionary perspective. This new discipline is called comparative functional genomics or comparative transcriptomics. This article reviews some of the main advances made in this field, using yeast species, and especially the species sequenced in the frame of the Genolevures program, as a model.

The reconstruction of condition-specific transcriptional modules provides new insights in the evolution of yeast AP-1 proteins

AP-1 proteins are transcription factors (TFs) that belong to the basic leucine zipper family, one of the largest families of TFs in eukaryotic cells. Despite high homology between their DNA binding domains, these proteins are able to recognize diverse DNA motifs. In yeasts, these motifs are referred as YRE (Yap Response Element) and are either seven (YRE-Overlap) or eight (YRE-Adjacent) base pair long. It has been proposed that the AP-1 DNA binding motif preference relies on a single change in the amino acid sequence of the yeast AP-1 TFs (an arginine in the YRE-O binding factors being replaced by a lysine in the YRE-A binding Yaps). We developed a computational approach to infer condition-specific transcriptional modules associated to the orthologous AP-1 protein Yap1p, Cgap1p and Cap1p, in three yeast species: the model yeast Saccharomyces cerevisiae and two pathogenic species Candida glabrata and Candida albicans. Exploitation of these modules in terms of predictions of the protein/DNA regulatory interactions changed our vision of AP-1 protein evolution. Cis-regulatory motif analyses revealed the presence of a conserved adenine in 5' position of the canonical YRE sites. While Yap1p, Cgap1p and Cap1p shared a remarkably low number of target genes, an impressive conservation was observed in the YRE sequences identified by Yap1p and Cap1p. In Candida glabrata, we found that Cgap1p, unlike Yap1p and Cap1p, recognizes YRE-O and YRE-A motifs. These findings were supported by structural data available for the transcription factor Pap1p (Schizosaccharomyces pombe). Thus, whereas arginine and lysine substitutions in Cgap1p and Yap1p proteins were reported as responsible for a specific YRE-O or YRE-A preference, our analyses rather suggest that the ancestral yeast AP-1 protein could recognize both YRE-O and YRE-A motifs and that the arginine/lysine exchange is not the only determinant of the specialization of modern Yaps for one motif or another.

Development of a novel multiplex DNA microarray for Fusarium graminearum and analysis of azole fungicide responses

Background

The toxigenic fungal plant pathogen Fusarium graminearum compromises wheat production worldwide. Azole fungicides play a prominent role in controlling this pathogen. Sequencing of its genome stimulated the development of high-throughput technologies to study mechanisms of coping with fungicide stress and adaptation to fungicides at a previously unprecedented precision. DNA-microarrays have been used to analyze genome-wide gene expression patterns and uncovered complex transcriptional responses. A recently developed one-color multiplex array format allowed flexible, effective, and parallel examinations of eight RNA samples.

Results

We took advantage of the 8 × 15 k Agilent format to design, evaluate, and apply a novel microarray covering the whole F. graminearum genome to analyze transcriptional responses to azole fungicide treatment. Comparative statistical analysis of expression profiles uncovered 1058 genes that were significantly differentially expressed after azole-treatment. Quantitative RT-PCR analysis for 31 selected genes indicated high conformity to results from the microarray hybridization. Among the 596 genes with significantly increased transcript levels, analyses using GeneOntology and FunCat annotations detected the ergosterol-biosynthesis pathway genes as the category most significantly responding, confirming the mode-of-action of azole fungicides. Cyp51A, which is one of the three F. graminearum paralogs of Cyp51 encoding the target of azoles, was the most consistently differentially expressed gene of the entire study. A molecular phylogeny analyzing the relationships of the three CYP51 proteins in the context of 38 fungal genomes belonging to the Pezizomycotina indicated that CYP51C (FGSG_11024) groups with a new clade of CYP51 proteins. The transcriptional profiles for genes encoding ABC transporters and transcription factors suggested several involved in mechanisms alleviating the impact of the fungicide. Comparative analyses with published microarray experiments obtained from two different nutritional stress conditions identified subsets of genes responding to different types of stress. Some of the genes that responded only to tebuconazole treatment appeared to be unique to the F. graminearum genome.

Conclusions

The novel F. graminearum 8 × 15 k microarray is a reliable and efficient high-throughput tool for genome-wide expression profiling experiments in fungicide research, and beyond, as shown by our data obtained for azole responses. The array data contribute to understanding mechanisms of fungicide resistance and allow identifying fungicide targets.

COMODO: an adaptive coclustering strategy to identify conserved coexpression modules between organisms

Increasingly large-scale expression compendia for different species are becoming available. By exploiting the modularity of the coexpression network, these compendia can be used to identify biological processes for which the expression behavior is conserved over different species. However, comparing module networks across species is not trivial. The definition of a biologically meaningful module is not a fixed one and changing the distance threshold that defines the degree of coexpression gives rise to different modules. As a result when comparing modules across species, many different partially overlapping conserved module pairs across species exist and deciding which pair is most relevant is hard. Therefore, we developed a method referred to as conserved modules across organisms (COMODO) that uses an objective selection criterium to identify conserved expression modules between two species. The method uses as input microarray data and a gene homology map and provides as output pairs of conserved modules and searches for the pair of modules for which the number of sharing homologs is statistically most significant relative to the size of the linked modules. To demonstrate its principle, we applied COMODO to study coexpression conservation between the two well-studied bacteria Escherichia coli and Bacillus subtilis. COMODO is available at: http://homes.esat.kuleuven.be/∼kmarchal/Supplementary_Information_Zarrineh_2010/comodo/index.html.

Comparative transcriptomics of drought responses in Populus: a meta-analysis of genome-wide expression profiling in mature leaves and root apices across two genotypes

Background

Comparative genomics has emerged as a promising means of unravelling the molecular networks underlying complex traits such as drought tolerance. Here we assess the genotype-dependent component of the drought-induced transcriptome response in two poplar genotypes differing in drought tolerance. Drought-induced responses were analysed in leaves and root apices and were compared with available transcriptome data from other Populus species.

Results

Using a multi-species designed microarray, a genomic DNA-based selection of probesets provided an unambiguous between-genotype comparison. Analyses of functional group enrichment enabled the extraction of processes physiologically relevant to drought response. The drought-driven changes in gene expression occurring in root apices were consistent across treatments and genotypes. For mature leaves, the transcriptome response varied weakly but in accordance with the duration of water deficit. A differential clustering algorithm revealed similar and divergent gene co-expression patterns among the two genotypes. Since moderate stress levels induced similar physiological responses in both genotypes, the genotype-dependent transcriptional responses could be considered as intrinsic divergences in genome functioning. Our meta-analysis detected several candidate genes and processes that are differentially regulated in root and leaf, potentially under developmental control, and preferentially involved in early and long-term responses to drought.

Conclusions

In poplar, the well-known drought-induced activation of sensing and signalling cascades was specific to the early response in leaves but was found to be general in root apices. Comparing our results to what is known in arabidopsis, we found that transcriptional remodelling included signalling and a response to energy deficit in roots in parallel with transcriptional indices of hampered assimilation in leaves, particularly in the drought-sensitive poplar genotype.

Cross-species queries of large gene expression databases

MOTIVATION

Expression databases, including the Gene Expression Omnibus and ArrayExpress, have experienced significant growth over the past decade and now hold hundreds of thousands of arrays from multiple species. Since most drugs are initially tested on model organisms, the ability to compare expression experiments across species may help identify pathways that are activated in a similar way in humans and other organisms. However, while several methods exist for finding co-expressed genes in the same species as a query gene, looking at co-expression of homologs or arbitrary genes in other species is challenging. Unlike sequence, which is static, expression is dynamic and changes between tissues, conditions and time. Thus, to carry out cross-species analysis using these databases, we need methods that can match experiments in one species with experiments in another species.

RESULTS

To facilitate queries in large databases, we developed a new method for comparing expression experiments from different species. We define a distance metric between the ranking of orthologous genes in the two species. We show how to solve an optimization problem for learning the parameters of this function using a training dataset of known similar expression experiments pairs. The function we learn outperforms previous methods and simpler rank comparison methods that have been used in the past for single species analysis. We used our method to compare millions of array pairs from mouse and human expression experiments. The resulting matches can be used to find functionally related genes, to hypothesize about biological response mechanisms and to highlight conditions and diseases that are activating similar pathways in both species.

AVAILABILITY

Supporting methods, results and a Matlab implementation are available from http://sb.cs.cmu.edu/ExpQ/.

Evolutionary divergence in the fungal response to fluconazole revealed by soft clustering

Background

Fungal infections are an emerging health risk, especially those involving yeast that are resistant to antifungal agents. To understand the range of mechanisms by which yeasts can respond to anti-fungals, we compared gene expression patterns across three evolutionarily distant species - Saccharomyces cerevisiae, Candida glabrata and Kluyveromyces lactis - over time following fluconazole exposure.

Results

Conserved and diverged expression patterns were identified using a novel soft clustering algorithm that concurrently clusters data from all species while incorporating sequence orthology. The analysis suggests complementary strategies for coping with ergosterol depletion by azoles - Saccharomyces imports exogenous ergosterol, Candida exports fluconazole, while Kluyveromyces does neither, leading to extreme sensitivity. In support of this hypothesis we find that only Saccharomyces becomes more azole resistant in ergosterol-supplemented media; that this depends on sterol importers Aus1 and Pdr11; and that transgenic expression of sterol importers in Kluyveromyces alleviates its drug sensitivity.

Conclusions

We have compared the dynamic transcriptional responses of three diverse yeast species to fluconazole treatment using a novel clustering algorithm. This approach revealed significant divergence among regulatory programs associated with fluconazole sensitivity. In future, such approaches might be used to survey a wider range of species, drug concentrations and stimuli to reveal conserved and divergent molecular response pathways.

Genome-wide characterisation of the Gcn5 histone acetyltransferase in budding yeast during stress adaptation reveals evolutionarily conserved and diverged roles

Background

Gcn5 is a transcriptional coactivator with histone acetyltransferase activity that is conserved with regard to structure as well as its histone substrates throughout the eukaryotes. Gene regulatory networks within cells are thought to be evolutionarily diverged. The use of evolutionarily divergent yeast species, such as S. cerevisiae and S. pombe, which can be studied under similar environmental conditions, provides an opportunity to examine the interface between conserved regulatory components and their cellular applications in different organisms.

Results

We show that Gcn5 is important for a common set of stress responses in evolutionarily diverged yeast species and that the activity of the conserved histone acetyltransferase domain is required. We define a group of KCl stress response genes in S. cerevisiae that are specifically dependent on Gcn5. Gcn5 is localised to many Gcn5-dependent genes including Gcn5 repressed targets such as FLO8. Gcn5 regulates divergent sets of KCl responsive genes in S. cerevisiae and S. pombe. Genome-wide localization studies showed a tendency for redistribution of Gcn5 during KCl stress adaptation in S. cerevisiae from short genes to the transcribed regions of long genes. An analogous redistribution was not observed in S. pombe.

Conclusions

Gcn5 is required for the regulation of divergent sets of KCl stress-response genes in S. cerevisiae and S. pombe even though it is required a common group of stress responses, including the response to KCl. Genes that are physically associated with Gcn5 require its activity for their repression or activation during stress adaptation, providing support for a role of Gcn5 as a corepressor as well as a coactivator. The tendency of Gcn5 to re-localise to the transcribed regions of long genes during KCl stress adaptation suggests that Gcn5 plays a specific role in the expression of long genes under adaptive conditions, perhaps by regulating transcriptional elongation as has been seen for Gcn5 in S. pombe. Interestingly an analogous redistribution of Gcn5 is not seen in S. pombe. The study thus provides important new insights in relation to why coregulators like Gcn5 are required for the correct expression of some genes but not others.

Heat-shock responsive genes identified and validated in Atlantic cod (Gadus morhua) liver, head kidney and skeletal muscle using genomic techniques

BACKGROUND

Daily and seasonal changes in temperature are challenges that fish within aquaculture settings cannot completely avoid, and are known to elicit complex organismal and cellular stress responses. We conducted a large-scale gene discovery and transcript expression study in order to better understand the genes that are potentially involved in the physiological and cellular aspects of stress caused by heat-shock. We used suppression subtractive hybridization (SSH) cDNA library construction and characterization to identify transcripts that were dysregulated by heat-shock in liver, skeletal muscle and head kidney of Atlantic cod. These tissues were selected due to their roles in metabolic regulation, locomotion and growth, and immune function, respectively. Fish were exposed for 3 hours to an 8 degrees C elevation in temperature, and then allowed to recover for 24 hours at the original temperature (i.e. 10 degrees C). Tissue samples obtained before heat-shock (BHS), at the cessation of heat-shock (CS), and 3, 12, and 24 hours after the cessation of heat-shock (ACS), were used for reciprocal SSH library construction and quantitative reverse transcription - polymerase chain reaction (QPCR) analysis of gene expression using samples from a group that was transferred but not heat-shocked (CT) as controls.

RESULTS

We sequenced and characterized 4394 ESTs (1524 from liver, 1451 from head kidney and 1419 from skeletal muscle) from three "forward subtracted" libraries (enriched for genes up-regulated by heat-shock) and 1586 from the liver "reverse subtracted" library (enriched for genes down-regulated by heat-shock), for a total of 5980 ESTs. Several cDNAs encoding putative chaperones belonging to the heat-shock protein (HSP) family were found in these libraries, and "protein folding" was among the gene ontology (GO) terms with the highest proportion in the libraries. QPCR analysis of HSP90alpha and HSP70-1 (synonym: HSPA1A) mRNA expression showed significant up-regulation in all three tissues studied. These transcripts were more than 100-fold up-regulated in liver following heat-shock. We also identified HSP47, GRP78 and GRP94-like transcripts, which were significantly up-regulated in all 3 tissues studied. Toll-like receptor 22 (TLR22) transcript, found in the liver reverse SSH library, was shown by QPCR to be significantly down-regulated in the head kidney after heat-shock.

CONCLUSION

Chaperones are an important part of the cellular response to stress, and genes identified in this work may play important roles in resistance to thermal-stress. Moreover, the transcript for one key immune response gene (TLR22) was down-regulated by heat-shock, and this down-regulation may be a component of heat-induced immunosuppression.

Fungal regulatory evolution: cis and trans in the balance

Regulatory divergence is likely a major driving force in evolution. Comparative genomics is being increasingly used to infer the evolution of gene regulation. Ascomycota fungi are uniquely suited among eukaryotes for regulatory evolution studies, due to broad phylogenetic scope, many sequenced genomes, and tractability of genomic analysis. Here we review recent advances in the identification of the contribution of cis- and trans-factors to expression divergence. Whereas current strategies have led to the discovery of surprising signatures and mechanisms, we still understand very little about the adaptive role of regulatory evolution. Empirical studies including experimental evolution, comparative functional genomics and hybrid and engineered strains are showing early promise toward deciphering the contribution of regulatory divergence to adaptation.

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.

Cross species analysis of microarray expression data

MOTIVATION

Many biological systems operate in a similar manner across a large number of species or conditions. Cross-species analysis of sequence and interaction data is often applied to determine the function of new genes. In contrast to these static measurements, microarrays measure the dynamic, condition-specific response of complex biological systems. The recent exponential growth in microarray expression datasets allows researchers to combine expression experiments from multiple species to identify genes that are not only conserved in sequence but also operated in a similar way in the different species studied.

RESULTS

In this review we discuss the computational and technical challenges associated with these studies, the approaches that have been developed to address these challenges and the advantages of cross-species analysis of microarray data. We show how successful application of these methods lead to insights that cannot be obtained when analyzing data from a single species. We also highlight current open problems and discuss possible ways to address them.