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

Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Candida glabrata is an emerging human fungal pathogen that is frequently drug tolerant, resulting in difficulties in treatment and a higher mortality in immunocompromised patients. The calcium-activated protein phosphatase calcineurin plays critical roles in controlling drug tolerance, hyphal growth, and virulence in diverse fungal pathogens via distinct mechanisms involving survival in serum or growth at host temperature (37° and higher). Here, we comprehensively studied the calcineurin signaling cascade in C. glabrata and found novel and uncharacterized functions of calcineurin and its downstream target Crz1 in governing thermotolerance, intracellular architecture, and pathogenesis in murine ocular, urinary tract, and systemic infections. This represents a second independent origin of a role for calcineurin in thermotolerant growth of a major human fungal pathogen, distinct from that which arose independently in Cryptococcus neoformans. Calcineurin also promotes survival of C. glabrata in serum via mechanisms distinct from C. albicans and thereby enables establishment of tissue colonization in a murine systemic infection model. To understand calcineurin signaling in detail, we performed global transcript profiling analysis and identified calcineurin- and Crz1-dependent genes in C. glabrata involved in cell wall biosynthesis, heat shock responses, and calcineurin function. Regulators of calcineurin (RCN) are a novel family of calcineurin modifiers, and two members of this family were identified in C. glabrata: Rcn1 and Rcn2. Our studies demonstrate that Rcn2 expression is controlled by calcineurin and Crz1 to function as a feedback inhibitor of calcineurin in a circuit required for calcium tolerance in C. glabrata. In contrast, the calcineurin regulator Rcn1 activates calcineurin signaling. Interestingly, neither Rcn1 nor Rcn2 is required for virulence in a murine systemic infection model. Taken together, our findings show that calcineurin signaling plays critical roles in thermotolerance and virulence, and that Rcn1 and Rcn2 have opposing functions in controlling calcineurin signaling in C. glabrata.

As a tool in directed genome manipulations, site-specific recombination is a double-edged sword. Exquisite specificity, while highly desirable, makes it imperative that the target site be first inserted at the desired genomic locale before it can be manipulated. We describe a combination of computational and experimental strategies, based on the tyrosine recombinase Flp and its target site FRT, to overcome this impediment. We document the systematic evolution of Flp variants that can utilize, in a bacterial assay, two sites from the human interleukin 10 gene, IL10, as recombination substrates. Recombination competence on an end target site is acquired via chimeric sites containing mixed sequences from FRT and the genomic locus. This is the first time that a tyrosine site-specific recombinase has been coaxed successfully to perform DNA exchange within naturally occurring sequences derived from a foreign genomic context. We demonstrate the ability of an Flp variant to mediate integration of a reporter cassette in Escherichia coli via recombination at one of the IL10-derived sites.