A systems approach to delineate functions of paralogous transcription factors: role of the Yap family in the DNA damage response

Yap5 Competes With Hap4 for the Regulation of Iron Homeostasis Genes in the Human Pathogen Candida glabrata

The CCAAT-binding complex (CBC) is a conserved heterotrimeric transcription factor which, in fungi, requires additional regulatory subunits to act on transcription. In the pathogenic yeast Candida glabrata, CBC has a dual role. Together with the Hap4 regulatory subunit, it activates the expression of genes involved in respiration upon growth with non-fermentable carbon sources, while its association with the Yap5 regulatory subunit is required for the activation of iron tolerance genes in response to iron excess. In the present work, we investigated further the interplay between CBC, Hap4 and Yap5. We showed that Yap5 regulation requires a specific Yap Response Element in the promoter of its target gene GRX4 and that the presence of Yap5 considerably strengthens the binding of CBC to the promoters of iron tolerance genes. Chromatin immunoprecipitation (ChIP) and transcriptome experiments showed that Hap4 can also bind these promoters but has no impact on the expression of those genes when Yap5 is present. In the absence of Yap5 however, GRX4 is constitutively regulated by Hap4, similarly to the genes involved in respiration. Our results suggest that the distinction between the two types of CBC targets in C. glabrata is mainly due to the dependency of Yap5 for very specific DNA sequences and to the competition between Hap4 and Yap5 at the promoter of the iron tolerance genes.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.

DNA damage induces Yap5-dependent transcription of ECO1/CTF7 in Saccharomyces cerevisiae

Yeast Eco1 (ESCO2 in humans) acetyltransferase converts chromatin-bound cohesins to a DNA tethering state, thereby establishing sister chromatid cohesion. Eco1 establishes cohesion during DNA replication, after which Eco1 is targeted for degradation by SCF E3 ubiquitin ligase. SCF E3 ligase, and sequential phosphorylations that promote Eco1 ubiquitination and degradation, remain active throughout the M phase. In this way, Eco1 protein levels are high during S phase, but remain low throughout the remaining cell cycle. In response to DNA damage during M phase, however, Eco1 activity increases-providing for a new wave of cohesion establishment (termed Damage-Induced Cohesion, or DIC) which is critical for efficient DNA repair. To date, little evidence exists as to the mechanism through which Eco1 activity increases during M phase in response to DNA damage. Possibilities include that either the kinases or E3 ligase, that target Eco1 for degradation, are inhibited in response to DNA damage. Our results reveal instead that the degradation machinery remains fully active during M phase, despite the presence of DNA damage. In testing alternate models through which Eco1 activity increases in response to DNA damage, the results reveal that DNA damage induces new transcription of ECO1 and at a rate that exceeds the rate of Eco1 turnover, providing for rapid accumulation of Eco1 protein. We further show that DNA damage induction of ECO1 transcription is in part regulated by Yap5-a stress-induced transcription factor. Given the role for mutated ESCO2 (homolog of ECO1) in human birth defects, this study highlights the complex nature through which mutation of ESCO2, and defects in ESCO2 regulation, may promote developmental abnormalities and contribute to various diseases including cancer.

Genome-Wide Identification and Expression Analysis of the bZIP Transcription Factors in the Mycoparasite Coniothyrium minitans

The basic leucine zipper (bZIP) proteins family is one of the largest and most diverse transcription factors, widely distributed in eukaryotes. However, no information is available regarding the bZIP gene family in Coniothyrium minitans, an important biocontrol agent of the plant pathogen Sclerotinia sclerotiorum. In this study, we identified 34 bZIP genes from the C. minitans genome, which were classified into 8 groups based on their phylogenetic relationships. Intron analysis showed that 28 CmbZIP genes harbored a variable number of introns, and 15 of them shared a feature that intron inserted into the bZIP domain. The intron position in bZIP domain was highly conserved, which was related to recognize the arginine (R) and could be treated as a genomic imprinting. Expression analysis of the CmbZIP genes in response to abiotic stresses indicated that they might play distinct roles in abiotic stress responses. Results showed that 22 CmbZIP genes were upregulated during the later stage of conidial development. Furthermore, transcriptome analysis indicated that CmbZIP genes are involved in different stages of mycoparasitism. Among deletion mutants of four CmbZIPs (CmbZIP07, -09, -13, and -16), only ΔCmbZIP16 mutants significantly reduced its tolerance to the oxidative stress. The other mutants exhibited no significant effects on colony morphology, mycelial growth, conidiation, and mycoparasitism. Taken together, our results suggested that CmbZIP genes play important roles in the abiotic stress responses, conidial development, and mycoparasitism. These results provide comprehensive information of the CmbZIP gene family and lay the foundation for further research on the bZIP gene family regarding their biological functions and evolutionary history.

The fungal CCAAT-binding complex and HapX display highly variable but evolutionary conserved synergetic promoter-specific DNA recognition

To sustain iron homeostasis, microorganisms have evolved fine-tuned mechanisms for uptake, storage and detoxification of the essential metal iron. In the human pathogen Aspergillus fumigatus, the fungal-specific bZIP-type transcription factor HapX coordinates adaption to both iron starvation and iron excess and is thereby crucial for virulence. Previous studies indicated that a HapX homodimer interacts with the CCAAT-binding complex (CBC) to cooperatively bind bipartite DNA motifs; however, the mode of HapX-DNA recognition had not been resolved. Here, combination of in vivo (genetics and ChIP-seq), in vitro (surface plasmon resonance) and phylogenetic analyses identified an astonishing plasticity of CBC:HapX:DNA interaction. DNA motifs recognized by the CBC:HapX protein complex comprise a bipartite DNA binding site 5′-CSAATN₁₂RWT-3′ and an additional 5′-TKAN-3′ motif positioned 11–23 bp downstream of the CCAAT motif, i.e. occasionally overlapping the 3′-end of the bipartite binding site. Phylogenetic comparison taking advantage of 20 resolved Aspergillus species genomes revealed that DNA recognition by the CBC:HapX complex shows promoter-specific cross-species conservation rather than regulon-specific conservation. Moreover, we show that CBC:HapX interaction is absolutely required for all known functions of HapX. The plasticity of the CBC:HapX:DNA interaction permits fine tuning of CBC:HapX binding specificities that could support adaptation of pathogens to their host niches.

The bZIP transcription factor FpAda1 is essential for fungal growth and conidiation in Fusarium pseudograminearum

Fusarium pseudograminearum is an important pathogen of Fusarium crown rot and Fusarium head blight, which is able to infect wheat and barley worldwide, causing great economic losses. Transcription factors (TFs) of the basic leucine zipper (bZIP) protein family control important processes in all eukaryotes. In this study, we identified a gene, designated FpAda1, encoding a bZIP TF in F. pseudograminearum. The homolog of FpAda1 is also known to affect hyphal growth in Neurospora crassa. Deletion of FpAda1 in F. pseudograminearum resulted in defects in hyphal growth, mycelial branching and conidia formation. Pathogenicity assays showed that virulence of the Δfpada1 mutant was dramatically decreased on wheat coleoptiles and barley leaves. However, wheat coleoptile inoculation assay showed that Δfpada1 could penetrate and proliferate in wheat cells. Moreover, the FpAda1 was required for abnormal nuclear morphology in conidia and transcription of FpCdc2 and FpCdc42. Taken together, these results indicate that FpAda1 is an important transcription factor involved in growth and development in F. pseudograminearum.

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.

The regulation of iron homeostasis in the fungal human pathogen Candida glabrata

Iron is an essential element to most microorganisms, yet an excess of iron is toxic. Hence, living cells have to maintain a tight balance between iron uptake and iron consumption and storage. The control of intracellular iron concentrations is particularly challenging for pathogens because mammalian organisms have evolved sophisticated high-affinity systems to sequester iron from microbes and because iron availability fluctuates among the different host niches. In this review, we present the current understanding of iron homeostasis and its regulation in the fungal pathogen Candida glabrata. This yeast is an emerging pathogen which has become the second leading cause of candidemia, a life-threatening invasive mycosis. C. glabrata is relatively poorly studied compared to the closely related model yeast Saccharomyces cerevisiae or to the pathogenic yeast Candida albicans. Still, several research groups have started to identify the actors of C. glabrata iron homeostasis and its transcriptional and post-transcriptional regulation. These studies have revealed interesting particularities of C. glabrata and have shed new light on the evolution of fungal iron homeostasis.

Doxorubicin induces an extensive transcriptional and metabolic rewiring in yeast cells

Doxorubicin is one of the most effective chemotherapy drugs used against solid tumors in the treatment of several cancer types. Two different mechanisms, (i) intercalation of doxorubicin into DNA and inhibition of topoisomerase II leading to changes in chromatin structure, (ii) generation of free radicals and oxidative damage to biomolecules, have been proposed to explain the mode of action of this drug in cancer cells. A genome-wide integrative systems biology approach used in the present study to investigate the long-term effect of doxorubicin in Saccharomyces cerevisiae cells indicated the up-regulation of genes involved in response to oxidative stress as well as in Rad53 checkpoint sensing and signaling pathway. Modular analysis of the active sub-network has also revealed the induction of the genes significantly associated with nucleosome assembly/disassembly and DNA repair in response to doxorubicin. Furthermore, an extensive re-wiring of the metabolism was observed. In addition to glycolysis, and sulfate assimilation, several pathways related to ribosome biogenesis/translation, amino acid biosynthesis, nucleotide biosynthesis, de novo IMP biosynthesis and one-carbon metabolism were significantly repressed. Pentose phosphate pathway, MAPK signaling pathway biological processes associated with meiosis and sporulation were found to be induced in response to long-term exposure to doxorubicin in yeast cells.

The CCAAT-Binding Complex Controls Respiratory Gene Expression and Iron Homeostasis in Candida Glabrata

The CCAAT-binding complex (CBC) is a heterotrimeric transcription factor which is widely conserved in eukaryotes. In the model yeast S. cerevisiae, CBC positively controls the expression of respiratory pathway genes. This role involves interactions with the regulatory subunit Hap4. In many pathogenic fungi, CBC interacts with the HapX regulatory subunit to control iron homeostasis. HapX is a bZIP protein which only shares with Hap4 the Hap4Like domain (Hap4L) required for its interaction with CBC. Here, we show that CBC has a dual role in the pathogenic yeast C. glabrata. It is required, along with Hap4, for the constitutive expression of respiratory genes and it is also essential for the iron stress response, which is mediated by the Yap5 bZIP transcription factor. Interestingly, Yap5 contains a vestigial Hap4L domain. The mutagenesis of this domain severely reduced Yap5 binding to its targets and compromised its interaction with Hap5. Hence, Yap5, like HapX in other species, acts as a CBC regulatory subunit in the regulation of iron stress response. This work reveals new aspects of iron homeostasis in C. glabrata and of the evolution of the role of CBC and Hap4L-bZIP proteins in this process.

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.

Autoregulation of ZEB2 expression for zearalenone production in Fusarium graminearum

Several Fusarium species produce the polyketide mycotoxin zearalenone (ZEA), a causative agent of hyperestrogenic syndrome in animals that is often found in F. graminearum-infected cereals in temperate regions. The ZEA biosynthetic cluster genes PKS4, PKS13, ZEB1 and ZEB2 encode a reducing polyketide synthase, a non-reducing polyketide synthase, an isoamyl alcohol oxidase and a transcription factor respectively. In this study, the production of two isoforms (ZEB2L and ZEB2S) from the ZEB2 gene in F. graminearum via an alternative promoter was characterized. ZEB2L contains a basic leucine zipper (bZIP) DNA-binding domain at the N-terminus, whereas ZEB2S is an N-terminally truncated form of ZEB2L that lacks the bZIP domain. Interestingly, ZEA triggers the induction of both ZEB2L and ZEB2S transcription. ZEB2L and ZEB2S interact with each other to form a heterodimer that regulates ZEA production by reducing the binding affinity of ZEB2L for the ZEB2L gene promoter. Our study provides insight into the autoregulation of ZEB2 expression by alternative promoter usage and a feedback loop during ZEA production; this regulatory mechanism is similar to that observed in higher eukaryotes.

Pervasive variation of transcription factor orthologs contributes to regulatory network evolution

Differences in transcriptional regulatory networks underlie much of the phenotypic variation observed across organisms. Changes to cis-regulatory elements are widely believed to be the predominant means by which regulatory networks evolve, yet examples of regulatory network divergence due to transcription factor (TF) variation have also been observed. To systematically ascertain the extent to which TFs contribute to regulatory divergence, we analyzed the evolution of the largest class of metazoan TFs, Cys2-His2 zinc finger (C2H2-ZF) TFs, across 12 Drosophila species spanning ~45 million years of evolution. Remarkably, we uncovered that a significant fraction of all C2H2-ZF 1-to-1 orthologs in flies exhibit variations that can affect their DNA-binding specificities. In addition to loss and recruitment of C2H2-ZF domains, we found diverging DNA-contacting residues in ~44% of domains shared between D. melanogaster and the other fly species. These diverging DNA-contacting residues, found in ~70% of the D. melanogaster C2H2-ZF genes in our analysis and corresponding to ~26% of all annotated D. melanogaster TFs, show evidence of functional constraint: they tend to be conserved across phylogenetic clades and evolve slower than other diverging residues. These same variations were rarely found as polymorphisms within a population of D. melanogaster flies, indicating their rapid fixation. The predicted specificities of these dynamic domains gradually change across phylogenetic distances, suggesting stepwise evolutionary trajectories for TF divergence. Further, whereas proteins with conserved C2H2-ZF domains are enriched in developmental functions, those with varying domains exhibit no functional enrichments. Our work suggests that a subset of highly dynamic and largely unstudied TFs are a likely source of regulatory variation in Drosophila and other metazoans.

Repositioning monocyte TFRN into fibroblasts

Repositioning of cell-specific Transcription Factor Regulatory Networks (TFRNs) into another cell type is a pragmatic approach to elicit desired functions in impaired cells and also to investigate transcription regulations underlying cellular functions. However, a systematic method to identify a set of transcription factors (TFs) which can efficiently promote TFRN in a foreign cell type is lacking. Here, we demonstrate a perturbation-matrix method for the identification and repositioning of monocyte TFRNs in dermal fibroblasts.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Transcriptional regulation in Saccharomyces cerevisiae: transcription factor regulation and function, mechanisms of initiation, and roles of activators and coactivators

Here we review recent advances in understanding the regulation of mRNA synthesis in Saccharomyces cerevisiae. Many fundamental gene regulatory mechanisms have been conserved in all eukaryotes, and budding yeast has been at the forefront in the discovery and dissection of these conserved mechanisms. Topics covered include upstream activation sequence and promoter structure, transcription factor classification, and examples of regulated transcription factor activity. We also examine advances in understanding the RNA polymerase II transcription machinery, conserved coactivator complexes, transcription activation domains, and the cooperation of these factors in gene regulatory mechanisms.

Determining PTEN functional status by network component deduced transcription factor activities

PTEN-controlled PI3K-AKT-mTOR pathway represents one of the most deregulated signaling pathways in human cancers. With many small molecule inhibitors that target PI3K-AKT-mTOR pathway being exploited clinically, sensitive and reliable ways of stratifying patients according to their PTEN functional status and determining treatment outcomes are urgently needed. Heterogeneous loss of PTEN is commonly associated with human cancers and yet PTEN can also be regulated on epigenetic, transcriptional or post-translational levels, which makes the use of simple protein or gene expression-based analyses in determining PTEN status less accurate. In this study, we used network component analysis to identify 20 transcription factors (TFs) whose activities deduced from their target gene expressions were immediately altered upon the re-expression of PTEN in a PTEN-inducible system. Interestingly, PTEN controls the activities (TFA) rather than the expression levels of majority of these TFs and these PTEN-controlled TFAs are substantially altered in prostate cancer mouse models. Importantly, the activities of these TFs can be used to predict PTEN status in human prostate, breast and brain tumor samples with enhanced reliability when compared to straightforward IHC-based or expression-based analysis. Furthermore, our analysis indicates that unique sets of PTEN-controlled TFAs significantly contribute to specific tumor types. Together, our findings reveal that TFAs may be used as “signatures” for predicting PTEN functional status and elucidate the transcriptional architectures underlying human cancers caused by PTEN loss.

Evolution of biological interaction networks: from models to real data

We are beginning to uncover common mechanisms leading to the evolution of biological networks. The driving force behind these advances is the increasing availability of comparative data in several species.

Curated collection of yeast transcription factor DNA binding specificity data reveals novel structural and gene regulatory insights

Background

Transcription factors (TFs) play a central role in regulating gene expression by interacting with cis-regulatory DNA elements associated with their target genes. Recent surveys have examined the DNA binding specificities of most Saccharomyces cerevisiae TFs, but a comprehensive evaluation of their data has been lacking.

Results

We analyzed in vitro and in vivo TF-DNA binding data reported in previous large-scale studies to generate a comprehensive, curated resource of DNA binding specificity data for all characterized S. cerevisiae TFs. Our collection comprises DNA binding site motifs and comprehensive in vitro DNA binding specificity data for all possible 8-bp sequences. Investigation of the DNA binding specificities within the basic leucine zipper (bZIP) and VHT1 regulator (VHR) TF families revealed unexpected plasticity in TF-DNA recognition: intriguingly, the VHR TFs, newly characterized by protein binding microarrays in this study, recognize bZIP-like DNA motifs, while the bZIP TF Hac1 recognizes a motif highly similar to the canonical E-box motif of basic helix-loop-helix (bHLH) TFs. We identified several TFs with distinct primary and secondary motifs, which might be associated with different regulatory functions. Finally, integrated analysis of in vivo TF binding data with protein binding microarray data lends further support for indirect DNA binding in vivo by sequence-specific TFs.

Conclusions

The comprehensive data in this curated collection allow for more accurate analyses of regulatory TF-DNA interactions, in-depth structural studies of TF-DNA specificity determinants, and future experimental investigations of the TFs' predicted target genes and regulatory roles.

YeTFaSCo: a database of evaluated yeast transcription factor sequence specificities

The yeast Saccharomyces cerevisiae is a prevalent system for the analysis of transcriptional networks. As a result, multiple DNA-binding sequence specificities (motifs) have been derived for most yeast transcription factors (TFs). However, motifs from different studies are often inconsistent with each other, making subsequent analyses complicated and confusing. Here, we have created YeTFaSCo (The Yeast Transcription Factor Specificity Compendium, http://yetfasco.ccbr.utoronto.ca/), an extensive collection of S. cerevisiae TF specificities. YeTFaSCo differs from related databases by being more comprehensive (including 1709 motifs for 256 proteins or protein complexes), and by evaluating the motifs using multiple objective quality metrics. The metrics include correlation between motif matches and ChIP-chip data, gene expression patterns, and GO terms, as well as motif agreement between different studies. YeTFaSCo also features an index of ‘expert-curated’ motifs, each associated with a confidence assessment. In addition, the database website features tools for motif analysis, including a sequence scanning function and precomputed genome-browser tracks of motif occurrences across the entire yeast genome. Users can also search the database for motifs that are similar to a query motif.

Transcriptional regulation and the diversification of metabolism in wine yeast strains

Transcription factors and their binding sites have been proposed as primary targets of evolutionary adaptation because changes to single transcription factors can lead to far-reaching changes in gene expression patterns. Nevertheless, there is very little concrete evidence for such evolutionary changes. Industrial wine yeast strains, of the species Saccharomyces cerevisiae, are a geno- and phenotypically diverse group of organisms that have adapted to the ecological niches of industrial winemaking environments and have been selected to produce specific styles of wine. Variation in transcriptional regulation among wine yeast strains may be responsible for many of the observed differences and specific adaptations to different fermentative conditions in the context of commercial winemaking. We analyzed gene expression profiles of wine yeast strains to assess the impact of transcription factor expression on metabolic networks. The data provide new insights into the molecular basis of variations in gene expression in industrial strains and their consequent effects on metabolic networks important to wine fermentation. We show that the metabolic phenotype of a strain can be shifted in a relatively predictable manner by changing expression levels of individual transcription factors, opening opportunities to modify transcription networks to achieve desirable outcomes.

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.

Environment-responsive transcription factors bind subtelomeric elements and regulate gene silencing

Chromosome position analysis of ChIP-chip data revealed that several carbon source and stress-responsive yeast transcription factors conditionally bind subtelomeric X elements.

Integration of several microarray gene expression data sets showed that, in this context, the factors conditionally control the boundaries and strength of subtelomeric silencing.

Regulation of silencing by a fatty acid-responsive factor was found to be dependent on Sir2p and independent of Hda1p.

These findings provide a critical link for establishing the mechanisms by which telomere biology is coordinated with other cellular processes including responses to environmental stimuli, aging and adaptation.

It is well established that environmental conditions modulate gene expression through local binding of a variety of conditionally active transcription factors, each responsive to specific environmental cues. However, another prevalent mechanism of gene regulation in eukaryotic cells is the long-range control of groups of genes by chromatin modifications or other position-dependent mechanisms. One such phenomenon, gene silencing, is an important and evolutionarily conserved mode of regulation that controls expression of subtelomeric genes. These genes are enriched for stress response and metabolic genes and their regulation is controlled by the spreading of silencing molecules from chromosome ends (telomeres) into subtelomeric regions. Levels of subtelomeric silencing have been linked to cellular lifespan, and study of the regulation of silencing is fundamental to our understanding of human aging. The spread of silencing in subtelomeric regions is discontinuous, and is controlled by various genomic elements that can either relay and enhance silencing from telomeres (proto-silencing) or create boundaries that protect some genomic regions from silencing. In yeast, every subtelomeric region contains an X element that proto-silences centromere-proximal genes, and also insulates telomere-proximal genes from silencing.

In this paper, we identify a regulatory mechanism to control X element-mediated proto-silencing and insulating activities in response to environmental cues. The mechanism was identified using chromosome position analysis of microarray-based chromatin immunoprecipitation (ChIP-chip) data for environment-responsive TFs and genome-wide gene expression data under the same conditions. The mechanism involves the conditional association of environment-responsive transcription factors to X elements. The binding at X elements results in regulation of proto-silencing of centromere-proximal genes, or insulation of telomere-proximal genes (depending on the factor) in response to environmental stimuli related to stress response and metabolism. One example is shown below (Figure 4B). Transcription factor, Oaf1p, conditionally binds X elements in the presence of fatty acids and enhances proto-silencing specifically under this condition. Oaf1p and several other factors implicated here are known to control adjacent genes at intrachromosomal positions, suggesting their dual functionality in both gene-specific transcriptional regulation, and long-range position-dependent mechanism. Investigation of this mechanism during the response to fatty acid exposure showed that conditional proto-silencing activity is dependent on Sir2p, a molecule known to be involved in subtelomeric silencing related to aging. This study reveals a path cells can use to coordinate subtelomeric silencing related to aging with cellular environment, and with the activities of other cellular processes.

Subtelomeric chromatin is subject to evolutionarily conserved complex epigenetic regulation and is implicated in numerous aspects of cellular function including formation of heterochromatin, regulation of stress response pathways and control of lifespan. Subtelomeric DNA is characterized by the presence of specific repeated segments that serve to propagate silencing or to protect chromosomal regions from spreading epigenetic control. In this study, analysis of genome-wide chromatin immunoprecipitation and expression data, suggests that several yeast transcription factors regulate subtelomeric silencing in response to various environmental stimuli through conditional association with proto-silencing regions called X elements. In this context, Oaf1p, Rox1p, Gzf1p and Phd1p control the propagation of silencing toward centromeres in response to stimuli affecting stress responses and metabolism, whereas others, including Adr1p, Yap5p and Msn4p, appear to influence boundaries of silencing, regulating telomere-proximal genes in Y′ elements. The factors implicated here are known to control adjacent genes at intrachromosomal positions, suggesting their dual functionality. This study reveals a path for the coordination of subtelomeric silencing with cellular environment, and with activities of other cellular processes.

Exploring the bZIP transcription factor regulatory network in Neurospora crassa

Transcription factors (TFs) are key nodes of regulatory networks in eukaryotic organisms, including filamentous fungi such as Neurospora crassa. The 178 predicted DNA-binding TFs in N. crassa are distributed primarily among six gene families, which represent an ancient expansion in filamentous ascomycete genomes; 98 TF genes show detectable expression levels during vegetative growth of N. crassa, including 35 that show a significant difference in expression level between hyphae at the periphery versus hyphae in the interior of a colony. Regulatory networks within a species genome include paralogous TFs and their respective target genes (TF regulon). To investigate TF network evolution in N. crassa, we focused on the basic leucine zipper (bZIP) TF family, which contains nine members. We performed baseline transcriptional profiling during vegetative growth of the wild-type and seven isogenic, viable bZIP deletion mutants. We further characterized the regulatory network of one member of the bZIP family, NCU03905. NCU03905 encodes an Ap1-like protein (NcAp-1), which is involved in resistance to multiple stress responses, including oxidative and heavy metal stress. Relocalization of NcAp-1 from the cytoplasm to the nucleus was associated with exposure to stress. A comparison of the NcAp-1 regulon with Ap1-like regulons in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Candida albicans and Aspergillus fumigatus showed both conservation and divergence. These data indicate how N. crassa responds to stress and provide information on pathway evolution.

How Saccharomyces cerevisiae copes with toxic metals and metalloids

Toxic metals and metalloids are widespread in nature and can locally reach fairly high concentrations. To ensure cellular protection and survival in such environments, all organisms possess systems to evade toxicity and acquire tolerance. This review provides an overview of the molecular mechanisms that contribute to metal toxicity, detoxification and tolerance acquisition in budding yeast Saccharomyces cerevisiae. We mainly focus on the metals/metalloids arsenic, cadmium, antimony, mercury, chromium and selenium, and emphasize recent findings on sensing and signalling mechanisms and on the regulation of tolerance and detoxification systems that safeguard cellular and genetic integrity.

Comparative transcriptome profiling analyses during the lag phase uncover YAP1, PDR1, PDR3, RPN4, and HSF1 as key regulatory genes in genomic adaptation to the lignocellulose derived inhibitor HMF for Saccharomyces cerevisiae

BACKGROUND

The yeast Saccharomyces cerevisiae is able to adapt and in situ detoxify lignocellulose derived inhibitors such as furfural and HMF. The length of lag phase for cell growth in response to the inhibitor challenge has been used to measure tolerance of strain performance. Mechanisms of yeast tolerance at the genome level remain unknown. Using systems biology approach, this study investigated comparative transcriptome profiling, metabolic profiling, cell growth response, and gene regulatory interactions of yeast strains and selective gene deletion mutations in response to HMF challenges during the lag phase of growth.

RESULTS

We identified 365 candidate genes and found at least 3 significant components involving some of these genes that enable yeast adaptation and tolerance to HMF in yeast. First, functional enzyme coding genes such as ARI1, ADH6, ADH7, and OYE3, as well as gene interactions involved in the biotransformation and inhibitor detoxification were the direct driving force to reduce HMF damages in cells. Expressions of these genes were regulated by YAP1 and its closely related regulons. Second, a large number of PDR genes, mainly regulated by PDR1 and PDR3, were induced during the lag phase and the PDR gene family-centered functions, including specific and multiple functions involving cellular transport such as TPO1, TPO4, RSB1, PDR5, PDR15, YOR1, and SNQ2, promoted cellular adaptation and survival in order to cope with the inhibitor stress. Third, expressed genes involving degradation of damaged proteins and protein modifications such as SHP1 and SSA4, regulated by RPN4, HSF1, and other co-regulators, were necessary for yeast cells to survive and adapt the HMF stress. A deletion mutation strain Δrpn4 was unable to recover the growth in the presence of HMF.

CONCLUSIONS

Complex gene interactions and regulatory networks as well as co-regulations exist in yeast adaptation and tolerance to the lignocellulose derived inhibitor HMF. Both induced and repressed genes involving diversified functional categories are accountable for adaptation and energy rebalancing in yeast to survive and adapt the HMF stress during the lag phase of growth. Transcription factor genes YAP1, PDR1, PDR3, RPN4, and HSF1 appeared to play key regulatory rules for global adaptation in the yeast S. cerevisiae.

Coevolution within a transcriptional network by compensatory trans and cis mutations

Transcriptional networks have been shown to evolve very rapidly, prompting questions as to how such changes arise and are tolerated. Recent comparisons of transcriptional networks across species have implicated variations in the cis-acting DNA sequences near genes as the main cause of divergence. What is less clear is how these changes interact with trans-acting changes occurring elsewhere in the genetic circuit. Here, we report the discovery of a system of compensatory trans and cis mutations in the yeast AP-1 transcriptional network that allows for conserved transcriptional regulation despite continued genetic change. We pinpoint a single species, the fungal pathogen Candida glabrata, in which a trans mutation has occurred very recently in a single AP-1 family member, distinguishing it from its Saccharomyces ortholog. Comparison of chromatin immunoprecipitation profiles between Candida and Saccharomyces shows that, despite their different DNA-binding domains, the AP-1 orthologs regulate a conserved block of genes. This conservation is enabled by concomitant changes in the cis-regulatory motifs upstream of each gene. Thus, both trans and cis mutations have perturbed the yeast AP-1 regulatory system in such a way as to compensate for one another. This demonstrates an example of "coevolution" between a DNA-binding transcription factor and its cis-regulatory site, reminiscent of the coevolution of protein binding partners.

Yeast biological networks unfold the interplay of antioxidants, genome and phenotype, and reveal a novel regulator of the oxidative stress response

BACKGROUND

Identifying causative biological networks associated with relevant phenotypes is essential in the field of systems biology. We used ferulic acid (FA) as a model antioxidant to characterize the global expression programs triggered by this small molecule and decipher the transcriptional network controlling the phenotypic adaptation of the yeast Saccharomyces cerevisiae.

METHODOLOGY/PRINCIPAL FINDINGS

By employing a strict cut off value during gene expression data analysis, 106 genes were found to be involved in the cell response to FA, independent of aerobic or anaerobic conditions. Network analysis of the system guided us to a key target node, the FMP43 protein, that when deleted resulted in marked acceleration of cellular growth (∼15% in both minimal and rich media). To extend our findings to human cells and identify proteins that could serve as drug targets, we replaced the yeast FMP43 protein with its human ortholog BRP44 in the genetic background of the yeast strain Δfmp43. The conservation of the two proteins was phenotypically evident, with BRP44 restoring the normal specific growth rate of the wild type. We also applied homology modeling to predict the 3D structure of the FMP43 and BRP44 proteins. The binding sites in the homology models of FMP43 and BRP44 were computationally predicted, and further docking studies were performed using FA as the ligand. The docking studies demonstrated the affinity of FA towards both FMP43 and BRP44.

CONCLUSIONS

This study proposes a hypothesis on the mechanisms yeast employs to respond to antioxidant molecules, while demonstrating how phenome and metabolome yeast data can serve as biomarkers for nutraceutical discovery and development. Additionally, we provide evidence for a putative therapeutic target, revealed by replacing the FMP43 protein with its human ortholog BRP44, a brain protein, and functionally characterizing the relevant mutant strain.

Identification of response-modulated genetic interactions by sensitivity-based epistatic analysis

Background

High-throughput genomics has enabled the global mapping of genetic interactions based on the phenotypic impact of combinatorial genetic perturbations. An important next step is to understand how these networks are dynamically remodelled in response to environmental stimuli. Here, we report on the development and testing of a method to identify such interactions. The method was developed from first principles by treating the impact on cellular growth of environmental perturbations equivalently to that of gene deletions. This allowed us to establish a novel neutrality function marking the absence of epistasis in terms of sensitivity phenotypes rather than fitness. We tested the method by identifying fitness- and sensitivity-based interactions involved in the response to drug-induced DNA-damage of budding yeast Saccharomyces cerevisiae using two mutant libraries - one containing transcription factor deletions, and the other containing deletions of DNA repair genes.

Results

Within the library of transcription factor deletion mutants, we observe significant differences in the sets of genetic interactions identified by the fitness- and sensitivity-based approaches. Notably, among the most likely interactions, only ~50% were identified by both methods. While interactions identified solely by the sensitivity-based approach are modulated in response to drug-induced DNA damage, those identified solely by the fitness-based method remained invariant to the treatment. Comparison of the identified interactions to transcriptional profiles and protein-DNA interaction data indicate that the sensitivity-based method improves the identification of interactions involved in the DNA damage response. Additionally, for the library containing DNA repair mutants, we observe that the sensitivity-based method improves the grouping of functionally related genes, as well as the identification of protein complexes, involved in DNA repair.

Conclusion

Our results show that the identification of response-modulated genetic interactions can be improved by incorporating the effect of a changing environment directly into the neutrality function marking the absence of epistasis. We expect that this extension of conventional epistatic analysis will facilitate the development of dynamic models of gene networks from quantitative measurements of genetic interactions. While the method was developed for growth phenotype, it should apply equally well for other phenotypes, including the expression of fluorescent reporters.

An integrative multi-network and multi-classifier approach to predict genetic interactions

Genetic interactions occur when a combination of mutations results in a surprising phenotype. These interactions capture functional redundancy, and thus are important for predicting function, dissecting protein complexes into functional pathways, and exploring the mechanistic underpinnings of common human diseases. Synthetic sickness and lethality are the most studied types of genetic interactions in yeast. However, even in yeast, only a small proportion of gene pairs have been tested for genetic interactions due to the large number of possible combinations of gene pairs. To expand the set of known synthetic lethal (SL) interactions, we have devised an integrative, multi-network approach for predicting these interactions that significantly improves upon the existing approaches. First, we defined a large number of features for characterizing the relationships between pairs of genes from various data sources. In particular, these features are independent of the known SL interactions, in contrast to some previous approaches. Using these features, we developed a non-parametric multi-classifier system for predicting SL interactions that enabled the simultaneous use of multiple classification procedures. Several comprehensive experiments demonstrated that the SL-independent features in conjunction with the advanced classification scheme led to an improved performance when compared to the current state of the art method. Using this approach, we derived the first yeast transcription factor genetic interaction network, part of which was well supported by literature. We also used this approach to predict SL interactions between all non-essential gene pairs in yeast (http://sage.fhcrc.org/downloads/downloads/predicted_yeast_genetic_interactions.zip). This integrative approach is expected to be more effective and robust in uncovering new genetic interactions from the tens of millions of unknown gene pairs in yeast and from the hundreds of millions of gene pairs in higher organisms like mouse and human, in which very few genetic interactions have been identified to date.

Finding regulatory DNA motifs using alignment-free evolutionary conservation information

As an increasing number of eukaryotic genomes are being sequenced, comparative studies aimed at detecting regulatory elements in intergenic sequences are becoming more prevalent. Most comparative methods for transcription factor (TF) binding site discovery make use of global or local alignments of orthologous regulatory regions to assess whether a particular DNA site is conserved across related organisms, and thus more likely to be functional. Since binding sites are usually short, sometimes degenerate, and often independent of orientation, alignment algorithms may not align them correctly. Here, we present a novel, alignment-free approach for using conservation information for TF binding site discovery. We relax the definition of conserved sites: we consider a DNA site within a regulatory region to be conserved in an orthologous sequence if it occurs anywhere in that sequence, irrespective of orientation. We use this definition to derive informative priors over DNA sequence positions, and incorporate these priors into a Gibbs sampling algorithm for motif discovery. Our approach is simple and fast. It requires neither sequence alignments nor the phylogenetic relationships between the orthologous sequences, yet it is more effective on real biological data than methods that do.

An atlas of combinatorial transcriptional regulation in mouse and man

Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

Mechanisms of transcription factor selectivity

The initiation of transcription is regulated by transcription factors (TFs) binding to DNA response elements (REs). How do TFs recognize specific binding sites among the many similar ones available in the genome? Recent research has illustrated that even a single nucleotide substitution can alter the selective binding of TFs to coregulators, that prior binding events can lead to selective DNA binding, and that selectivity is influenced by the availability of binding sites in the genome. Here, we combine structural insights with recent genomics screens to address the problem of TF-DNA interaction specificity. The emerging picture of selective binding site sequence recognition and TF activation involves three major factors: the cellular network, protein and DNA as dynamic conformational ensembles and the tight packing of multiple TFs and coregulators on stretches of regulatory DNA. The classification of TF recognition mechanisms based on these factors impacts our understanding of how transcription initiation is regulated.

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.

A multiparameter network reveals extensive divergence between C. elegans bHLH transcription factors

Differences in expression, protein interactions, and DNA binding of paralogous transcription factors ("TF parameters") are thought to be important determinants of regulatory and biological specificity. However, both the extent of TF divergence and the relative contribution of individual TF parameters remain undetermined. We comprehensively identify dimerization partners, spatiotemporal expression patterns, and DNA-binding specificities for the C. elegans bHLH family of TFs, and model these data into an integrated network. This network displays both specificity and promiscuity, as some bHLH proteins, DNA sequences, and tissues are highly connected, whereas others are not. By comparing all bHLH TFs, we find extensive divergence and that all three parameters contribute equally to bHLH divergence. Our approach provides a framework for examining divergence for other protein families in C. elegans and in other complex multicellular organisms, including humans. Cross-species comparisons of integrated networks may provide further insights into molecular features underlying protein family evolution. For a video summary of this article, see the PaperFlick file available with the online Supplemental Data.

Dynamic and complex transcription factor binding during an inducible response in yeast

Complex biological processes are often regulated, at least in part, by the binding of transcription factors to their targets. Recently, considerable effort has been made to analyze the binding of relevant factors to the suite of targets they regulate, thereby generating a regulatory circuit map. However, for most studies the dynamics of binding have not been analyzed, and thus the temporal order of events and mechanisms by which this occurs are poorly understood. We globally analyzed in detail the temporal order of binding of several key factors involved in the salt response of yeast to their target genes. Analysis of Yap4 and Sko1 binding to their target genes revealed multiple temporal classes of binding patterns: (1) constant binding, (2) rapid induction, (3) slow induction, and (4) transient induction. These results demonstrate that individual transcription factors can have multiple binding patterns and help define the different types of temporal binding patterns used in eukaryotic gene regulation. To investigate these binding patterns further, we also analyzed the binding of seven other key transcription factors implicated in osmotic regulation, including Hot1, Msn1, Msn2, Msn4, Skn7, and Yap6, and found significant coassociation among the different factors at their gene targets. Moreover, the binding of several key factors was correlated with distinct classes of Yap4- and Sko1-binding patterns and with distinct types of genes. Gene expression studies revealed association of Yap4, Sko1, and other transcription factor-binding patterns with different gene expression patterns. The integration and analysis of binding and expression information reveals a complex dynamic and hierarchical circuit in which specific combinations of transcription factors target distinct sets of genes at discrete times to coordinate a rapid and important biological response.

Adaptable gene-specific dye bias correction for two-channel DNA microarrays

DNA microarray technology is a powerful tool for monitoring gene expression or for finding the location of DNA-bound proteins. DNA microarrays can suffer from gene-specific dye bias (GSDB), causing some probes to be affected more by the dye than by the sample. This results in large measurement errors, which vary considerably for different probes and also across different hybridizations. GSDB is not corrected by conventional normalization and has been difficult to address systematically because of its variance. We show that GSDB is influenced by label incorporation efficiency, explaining the variation of GSDB across different hybridizations. A correction method (Gene- And Slide-Specific Correction, GASSCO) is presented, whereby sequence-specific corrections are modulated by the overall bias of individual hybridizations. GASSCO outperforms earlier methods and works well on a variety of publically available datasets covering a range of platforms, organisms and applications, including ChIP on chip. A sequence-based model is also presented, which predicts which probes will suffer most from GSDB, useful for microarray probe design and correction of individual hybridizations. Software implementing the method is publicly available.

A fitness-based interferential genetics approach using hypertoxic/inactive gene alleles as references

Genetics, genomics, and biochemistry have all been of immense help in characterizing macromolecular cell entities and their interactions. Still, obtaining an overall picture of the functioning of even a simple unicellular species has remained a challenging task. One possible way to obtain a comprehensive picture has been described: by capitalizing on the observation that the overexpression on a multicopy plasmid of apparently any wild-type gene in yeast can lead to some negative effect on cell fitness (referring to the concept of "gene toxicity"), the FIG (fitness-based interferential genetics) approach was devised for selecting normal genes that are in antagonistic (and potentially also agonistic) relationship with a particular gene used as a reference. Herein, we take a complementary approach to FIG, by first selecting a "hypertoxic" allele of the reference gene--which easily provides the general possibility of obtaining gene products with the remarkable property of being inactive without altering their macromolecular interactivity--and then looking for the genes that interact functionally with this reference. Thus, FIG and the present approach (Trap-FIG), both taking advantage of the negative effects on cell fitness induced by various quantitative modulations in cellular networks, could potentially pave the way for the emergence of efficient in situ biochemistry.

Characterization of the DNA-binding motif of the arsenic-responsive transcription factor Yap8p

Saccharomyces cerevisiae uses several mechanisms for arsenic detoxification including the arsenate reductase Acr2p and the arsenite efflux protein Acr3p. ACR2 and ACR3 are transcribed in opposite directions from the same promoter and expression of these genes is regulated by the AP-1 (activator protein 1)-like transcription factor Yap8p. Yap8p has been shown to permanently associate with this promoter and to stimulate ACR2/ACR3 expression in response to arsenic. In the present study we characterized the DNA sequence that is targeted by Yap8p. We show that Yap8p binds to a pseudo-palindromic TGATTAATAATCA sequence that is related to, but distinct from, the sequence recognized by other fungal AP-1 proteins. Probing the promoter by mutational analysis, we confirm the importance of the TTAATAA core element and pin-point nucleotides that flank this element as crucial for Yap8p binding and in vivo activation of ACR3 expression. A genome-wide search for this element combined with global gene expression analysis indicates that the principal function of Yap8p is to control expression of ACR2 and ACR3. We conclude that Yap8p and other yeast AP-1 proteins require distinct DNA-binding motifs to induce gene expression and propose that this fact contributed towards a separation of function between AP-1 proteins during evolution.