Genome-wide analysis of transcriptional dependence and probable target sites for Abf1 and Rap1 in Saccharomyces cerevisiae

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

A half century of exploring DNA excision repair in chromatin

DNA in eukaryotic cells is packaged into the compact and dynamic structure of chromatin. This packaging is a double-edged sword for DNA repair and genomic stability. Chromatin restricts the access of repair proteins to DNA lesions embedded in nucleosomes and higher order chromatin structures. However, chromatin also serves as a signaling platform in which post-translational modifications of histones and other chromatin-bound proteins promote lesion recognition and repair. Similarly, chromatin modulates the formation of DNA damage, promoting or suppressing lesion formation depending on the chromatin context. Therefore, the modulation of DNA damage and its repair in chromatin is crucial to our understanding of the fate of potentially mutagenic and carcinogenic lesions in DNA. Here, we survey many of the landmark findings on DNA damage and repair in chromatin over the last 50 years (i.e., since the beginning of this field), focusing on excision repair, the first repair mechanism studied in the chromatin landscape. For example, we highlight how the impact of chromatin on these processes explains the distinct patterns of somatic mutations observed in cancer genomes.

SIR telomere silencing depends on nuclear envelope lipids and modulates sensitivity to a lysolipid

The nuclear envelope (NE) is important in maintaining genome organization. The role of lipids in communication between the NE and telomere regulation was investigated, including how changes in lipid composition impact gene expression and overall nuclear architecture. Yeast was treated with the non-metabolizable lysophosphatidylcholine analog edelfosine, known to accumulate at the perinuclear ER. Edelfosine induced NE deformation and disrupted telomere clustering but not anchoring. Additionally, the association of Sir4 at telomeres decreased. RNA-seq analysis showed altered expression of Sir-dependent genes located at sub-telomeric (0-10 kb) regions, consistent with Sir4 dispersion. Transcriptomic analysis revealed that two lipid metabolic circuits were activated in response to edelfosine, one mediated by the membrane sensing transcription factors, Spt23/Mga2, and the other by a transcriptional repressor, Opi1. Activation of these transcriptional programs resulted in higher levels of unsaturated fatty acids and the formation of nuclear lipid droplets. Interestingly, cells lacking Sir proteins displayed resistance to unsaturated-fatty acids and edelfosine, and this phenotype was connected to Rap1.

A novel allelic donkey β-lactogobulin I protein isoform generated by a non-AUG translation initiation codon is associated with a nonsynonymous SNP

β-Lactoglobulin I (β-LG I) is one of the most important whey proteins in donkey milk. However, to our knowledge, there has been no study focusing on the full nucleotide sequences of this gene (BLG I). Current investigation of donkey BLG I gene is very limited with only 2 variants (A and B) characterized so far at the protein level. Recently, a new β-LG I variant, with a significantly higher mass (+1,915 Da) than known variants has been detected. In this study, we report the whole nucleotide sequence of the BLG I gene from 2 donkeys, whose milk samples are characterized by the β-LG I SDS-PAGE band with a normal electrophoretic mobility (18,514.25 Da, β-LG I B1 form) the first, and by the presence of a unique β-LG I band with a higher electrophoretic mobility (20,428.5 Da, β-LG I D form) the latter. A high genetic variability was found all over the 2 sequenced BLG I alleles. In particular, 16 polymorphic sites were found in introns, one in the 5' flanking region, 3 SNPs in the 5' untranslated region and one SNP in the coding region (g.458G > A) located at the 40th nucleotide of exon 2 and responsible for the AA substitutions p.Asp28 > Asn in the mature protein. Two SNPs (g.920-922CAC > TGT and g.1871G/A) were genotyped in 93 donkeys of 2 Italian breeds (60 Ragusana and 33 Amiatina, respectively) and the overall frequencies of g.920-922CAC and g.1871A were 0.3065 and 0.043, respectively. Only the rare allele g.1871A was observed to be associated with the slower migrating β-LG I. Considering this genetic diversity and those found in the database, it was possible to deduce at least 5 different alleles (BLG I A, B, B1, C, D) responsible for 4 potential β-LG I translations. Among these alleles, B1 and D are those characterized in the present research, with the D allele of real novel identification. Haplotype data analysis suggests an evolutionary pathway of donkey BLG I gene and a possible phylogenetic map is proposed. Analyses of mRNA secondary structure showed relevant changes in the structures, as consequence of the g.1871G > A polymorphism, that might be responsible for the recognition of an alternative initiation site providing an additional signal peptide. The extension of 19 AA sequence to the mature protein, corresponding to the canonical signal peptide with an additional alanine residue, is sufficient to provide the observed molecular weight of the slower migrating β-LG I encoded by the BLG I D allele.

Comparative whole transcriptome analysis of Parkinson's disease focusing on the efficacy of zonisamide

OBJECTIVE

Interindividual variations in responsiveness to zonisamide in patients with Parkinson's disease (PD) have been observed in clinical settings. To decipher the molecular mechanisms determining the efficacy of zonisamide, we conducted whole transcriptome sequencing analysis of patients with PD.

METHODS

We selected 23 super-responders (SRs) and 25 non-responders (NRs) to zonisamide from patients with PD who had participated in a previous clinical trial for the approval of zonisamide for the treatment of 'wearing-off'. Whole transcriptome analysis of peripheral blood was conducted on samples taken before and 12 weeks after zonisamide treatment. We performed differential gene expression analysis to compare between the SRs and NRs at each time point.

RESULTS

Differentially expressed genes in the pre-treatment samples were significantly enriched for glutamatergic synapses and insulin-like growth factor binding (P_adj=7.8 × 10^-3 and 0.029, respectively). The gene sets associated with these functions changed more dynamically by treatment in SRs than NRs (p=7.2 × 10^-3 and 8.2 × 10^-3, respectively).

CONCLUSIONS

Our results suggest that the efficacy of zonisamide in PD patients is associated with glutamate-related synaptic modulation and p53-mediated dopaminergic neural loss. Their transcriptomic differences could be captured before treatment, which would lead to the realisation of future personalised treatment.

Abf1 Is an Essential Protein That Participates in Cell Cycle Progression and Subtelomeric Silencing in Candida glabrata

Accurate DNA replication and segregation is key to reproduction and cell viability in all organisms. Autonomously replicating sequence-binding factor 1 (Abf1) is a multifunctional protein that has essential roles in replication, transcription, and regional silencing in the model yeast Saccharomyces cerevisiae. In the opportunistic pathogenic fungus Candida glabrata, which is closely related to S. cerevisiae, these processes are important for survival within the host, for example, the regulation of transcription of virulence-related genes like those involved in adherence. Here, we describe that CgABF1 is an essential gene required for cell viability and silencing near the telomeres, where many adhesin-encoding genes reside. CgAbf1 mediated subtelomeric silencing depends on the 43 C-terminal amino acids. We also found that abnormal expression, depletion, or overexpression of Abf1, results in defects in nuclear morphology, nuclear segregation, and transit through the cell cycle. In the absence of ABF1, cells are arrested in G2 but start cycling again after 9 h, coinciding with the loss of cell viability and the appearance of cells with higher DNA content. Overexpression of CgABF1 causes defects in nuclear segregation and cell cycle progression. We suggest that these effects could be due to the deregulation of DNA replication.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

Histone exchange is associated with activator function at transcribed promoters and with repression at histone loci

Transcription in eukaryotes correlates with major chromatin changes, including the replacement of old nucleosomal histones by new histones at the promoters of genes. The role of these histone exchange events in transcription remains unclear. In particular, the causal relationship between histone exchange and activator binding, preinitiation complex (PIC) assembly, and/or subsequent transcription remains unclear. Here, we provide evidence that histone exchange at gene promoters is not simply a consequence of PIC assembly or transcription but instead is mediated by activators. We further show that not all activators up-regulate gene expression by inducing histone turnover. Thus, histone exchange does not simply correlate with transcriptional activity, but instead reflects the mode of action of the activator. Last, we show that histone turnover is not only associated with activator function but also plays a role in transcriptional repression at the histone loci.

Nucleosome remodeling at origins of global genome-nucleotide excision repair occurs at the boundaries of higher-order chromatin structure

Repair of UV-induced DNA damage requires chromatin remodeling. How repair is initiated in chromatin remains largely unknown. We recently demonstrated that global genome–nucleotide excision repair (GG-NER) in chromatin is organized into domains in relation to open reading frames. Here, we define these domains, identifying the genomic locations from which repair is initiated. By examining DNA damage–induced changes in the linear structure of nucleosomes at these sites, we demonstrate how chromatin remodeling is initiated during GG-NER. In undamaged cells, we show that the GG-NER complex occupies chromatin, establishing the nucleosome structure at these genomic locations, which we refer to as GG-NER complex binding sites (GCBSs). We demonstrate that these sites are frequently located at genomic boundaries that delineate chromosomally interacting domains (CIDs). These boundaries define domains of higher-order nucleosome–nucleosome interaction. We demonstrate that initiation of GG-NER in chromatin is accompanied by the disruption of dynamic nucleosomes that flank GCBSs by the GG-NER complex.

Spt6 Is Required for the Fidelity of Promoter Selection

Spt6 is a conserved factor that controls transcription and chromatin structure across the genome. Although Spt6 is viewed as an elongation factor, spt6 mutations in Saccharomyces cerevisiae allow elevated levels of transcripts from within coding regions, suggesting that Spt6 also controls initiation. To address the requirements for Spt6 in transcription and chromatin structure, we have combined four genome-wide approaches. Our results demonstrate that Spt6 represses transcription initiation at thousands of intragenic promoters. We characterize these intragenic promoters and find sequence features conserved with genic promoters. Finally, we show that Spt6 also regulates transcription initiation at most genic promoters and propose a model of initiation site competition to account for this. Together, our results demonstrate that Spt6 controls the fidelity of transcription initiation throughout the genome.

Nascent RNA signaling to yeast RNA Pol II during transcription elongation

Transcription as the key step in gene expression is a highly regulated process. The speed of transcription elongation depends on the underlying gene sequence and varies on a gene by gene basis. The reason for this sequence dependence is not known in detail. Recently, our group studied the cross talk between the nascent RNA and the transcribing RNA polymerase by screening the Escherichia coli genome for RNA sequences with high affinity to RNA Pol by performing genomic SELEX. This approach led to the identification of RNA polymerase-binding APtamers termed "RAPs". RAPs can have positive and negative effects on gene expression. A subgroup is able to downregulate transcription via the activity of the termination factor Rho. In this study, we used a similar SELEX setup using yeast genomic DNA as source of RNA sequences and highly purified yeast RNA Pol II as bait and obtained almost 1300 yeast-derived RAPs. Yeast RAPs are found throughout the genome within genes and antisense to genes, they are overrepresented in the non-transcribed strand of yeast telomeres and underrepresented in intergenic regions. Genes harbouring a RAP are more likely to show lower mRNA levels. By determining the endogenous expression levels as well as using a reporter system, we show that RAPs located within coding regions can reduce the transcript level downstream of the RAP. Here we demonstrate that RAPs represent a novel type of regulatory RNA signal in Saccharomyces cerevisiae that act in cis and interfere with the elongating transcription machinery to reduce the transcriptional output.

Histone deacetylase-mediated regulation of endolysosomal pH

The pH of the endolysosomal system is tightly regulated by a balance of proton pump and leak mechanisms that are critical for storage, recycling, turnover, and signaling functions in the cell. Dysregulation of endolysosomal pH has been linked to aging, amyloidogenesis, synaptic dysfunction, and various neurodegenerative disorders, including Alzheimer's disease. Therefore, understanding the mechanisms that regulate luminal pH may be key to identifying new targets for managing these disorders. Meta-analysis of yeast microarray databases revealed that nutrient-limiting conditions inhibited the histone deacetylase (HDAC) Rpd3 and thereby up-regulated transcription of the endosomal Na⁺/H⁺ exchanger Nhx1, resulting in vacuolar alkalinization. Consistent with these findings, Rpd3 inhibition by the HDAC inhibitor and antifungal drug trichostatin A induced Nhx1 expression and vacuolar alkalinization. Bioinformatics analysis of Drosophila and mouse databases revealed that caloric control of the Nhx1 orthologs DmNHE3 and NHE6, respectively, is also mediated by HDACs. We show that NHE6 is a target of the transcription factor cAMP-response element-binding protein (CREB), a known regulator of cellular responses to low-nutrient conditions, providing a molecular mechanism for nutrient- and HDAC-dependent regulation of endosomal pH. Of note, pharmacological targeting of the CREB pathway to increase NHE6 expression helped regulate endosomal pH and correct defective clearance of amyloid Aβ in an apoE4 astrocyte model of Alzheimer's disease. These observations from yeast, fly, mouse, and cell culture models point to an evolutionarily conserved mechanism for HDAC-mediated regulation of endosomal NHE expression. Our insights offer new therapeutic strategies for modulation of endolysosomal pH in fungal infection and human disease.

A global function for transcription factors in assisting RNA polymerase II termination

The role of transcription factors (TFs) on nucleosome positioning, RNA polymerase recruitment, and transcription initiation has been extensively characterized. Here, we propose that a subset of TFs such as Reb1, Abf1, Rap1, and TFIIIB also serve a major function in partitioning transcription units by assisting the Nrd1p-Nab3p-Sen1p Pol II termination pathway.

Abf1 and other general regulatory factors control ribosome biogenesis gene expression in budding yeast

Ribosome biogenesis in Saccharomyces cerevisiae involves a regulon of >200 genes (Ribi genes) coordinately regulated in response to nutrient availability and cellular growth rate. Two cis-acting elements called PAC and RRPE are known to mediate Ribi gene repression in response to nutritional downshift. Here, we show that most Ribi gene promoters also contain binding sites for one or more General Regulatory Factors (GRFs), most frequently Abf1 and Reb1, and that these factors are enriched in vivo at Ribi promoters. Abf1/Reb1/Tbf1 promoter association was required for full Ribi gene expression in rich medium and for its modulation in response to glucose starvation, characterized by a rapid drop followed by slow recovery. Such a response did not entail changes in Abf1 occupancy, but it was paralleled by a quick increase, followed by slow decrease, in Rpd3L histone deacetylase occupancy. Remarkably, Abf1 site disruption also abolished Rpd3L complex recruitment in response to starvation. Extensive mutational analysis of the DBP7 promoter revealed a complex interplay of Tbf1 sites, PAC and RRPE in the transcriptional regulation of this Ribi gene. Our observations point to GRFs as new multifaceted players in Ribi gene regulation both during exponential growth and under repressive conditions.

Identification of a transcriptional activation domain in yeast repressor activator protein 1 (Rap1) using an altered DNA-binding specificity variant

Repressor activator protein 1 (Rap1) performs multiple vital cellular functions in the budding yeast Saccharomyces cerevisiae These include regulation of telomere length, transcriptional repression of both telomere-proximal genes and the silent mating type loci, and transcriptional activation of hundreds of mRNA-encoding genes, including the highly transcribed ribosomal protein- and glycolytic enzyme-encoding genes. Studies of the contributions of Rap1 to telomere length regulation and transcriptional repression have yielded significant mechanistic insights. However, the mechanism of Rap1 transcriptional activation remains poorly understood because Rap1 is encoded by a single copy essential gene and is involved in many disparate and essential cellular functions, preventing easy interpretation of attempts to directly dissect Rap1 structure-function relationships. Moreover, conflicting reports on the ability of Rap1-heterologous DNA-binding domain fusion proteins to serve as chimeric transcriptional activators challenge use of this approach to study Rap1. Described here is the development of an altered DNA-binding specificity variant of Rap1 (Rap1AS). We used Rap1AS to map and characterize a 41-amino acid activation domain (AD) within the Rap1 C terminus. We found that this AD is required for transcription of both chimeric reporter genes and authentic chromosomal Rap1 enhancer-containing target genes. Finally, as predicted for a bona fide AD, mutation of this newly identified AD reduced the efficiency of Rap1 binding to a known transcriptional coactivator TFIID-binding target, Taf5. In summary, we show here that Rap1 contains an AD required for Rap1-dependent gene transcription. The Rap1AS variant will likely also be useful for studies of the functions of Rap1 in other biological pathways.

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

Global genome nucleotide excision repair is organized into domains that promote efficient DNA repair in chromatin

The rates at which lesions are removed by DNA repair can vary widely throughout the genome, with important implications for genomic stability. To study this, we measured the distribution of nucleotide excision repair (NER) rates for UV-induced lesions throughout the budding yeast genome. By plotting these repair rates in relation to genes and their associated flanking sequences, we reveal that, in normal cells, genomic repair rates display a distinctive pattern, suggesting that DNA repair is highly organized within the genome. Furthermore, by comparing genome-wide DNA repair rates in wild-type cells and cells defective in the global genome-NER (GG-NER) subpathway, we establish how this alters the distribution of NER rates throughout the genome. We also examined the genomic locations of GG-NER factor binding to chromatin before and after UV irradiation, revealing that GG-NER is organized and initiated from specific genomic locations. At these sites, chromatin occupancy of the histone acetyl-transferase Gcn5 is controlled by the GG-NER complex, which regulates histone H3 acetylation and chromatin structure, thereby promoting efficient DNA repair of UV-induced lesions. Chromatin remodeling during the GG-NER process is therefore organized into these genomic domains. Importantly, loss of Gcn5 significantly alters the genomic distribution of NER rates; this has implications for the effects of chromatin modifiers on the distribution of mutations that arise throughout the genome.

Dynamics of Nucleosome Positioning Maturation following Genomic Replication

Chromatin is thought to carry epigenetic information from one generation to the next, although it is unclear how such information survives the disruptions of nucleosomal architecture occurring during genomic replication. Here, we measure a key aspect of chromatin structure dynamics during replication—how rapidly nucleosome positions are established on the newly replicated daughter genomes. By isolating newly synthesized DNA marked with 5-ethynyl-2′-deoxyuridine (EdU), we characterize nucleosome positions on both daughter genomes of S. cerevisiae during chromatin maturation. We find that nucleosomes rapidly adopt their mid-log positions at highly transcribed genes, which is consistent with a role for transcription in positioning nucleosomes in vivo. Additionally, experiments in hir1Δ mutants reveal a role for HIR in nucleosome spacing. We also characterized nucleosome positions on the leading and lagging strands, uncovering differences in chromatin maturation dynamics at hundreds of genes. Our data define the maturation dynamics of newly replicated chromatin and support a role for transcription in sculpting the chromatin template.

•

Nucleosome positions are determined on newly replicated DNA

•

Transcription reorders nucleosomes in gene bodies after DNA replication

•

The HIR complex tightens nucleosome spacing in gene bodies following replication

•

Nucleosome positions on leading and lagging strands depend on genes’ orientation

A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation

R2R3-MYB transcription factors (TFs) belong to a large and functionally diverse protein superfamily in plants. In this study, we explore the evolution and function of this family in grapevine (Vitis vinifera L.), a high-value fruit crop. We identified and manually curated 134 genes using RNA-Seq data, and named them systematically according to the Super-Nomenclature Committee. We identified novel genes, splicing variants and grapevine/woody-specific duplicated subgroups, suggesting possible neo- and sub-functionalization events. Regulatory network analysis ascribed biological functions to uncharacterized genes and validated those of known genes (e.g. secondary cell wall biogenesis and flavonoid biosynthesis). A comprehensive analysis of different MYB binding motifs in the promoters of co-expressed genes predicted grape R2R3-MYB binding preferences and supported evidence for putative downstream targets. Enrichment of cis-regulatory motifs for diverse TFs reinforced the notion of transcriptional coordination and interaction between MYBs and other regulators. Analysis of the network of Subgroup 2 showed that the resveratrol-related VviMYB14 and VviMYB15 share common co-expressed STILBENE SYNTHASE genes with the uncharacterized VviMYB13. These regulators have distinct expression patterns within organs and in response to biotic and abiotic stresses, suggesting a pivotal role of VviMYB13 in regulating stilbene accumulation in vegetative tissues and under biotic stress conditions.

Multiple roles of the general regulatory factor Abf1 in yeast ribosome biogenesis

In Saccharomyces cerevisiae, the large majority of the genes coding for cytoplasmic ribosomal proteins (RPs) depend on the general regulatory factor Rap1 for their transcription, but a small cohort of them relies on Abf1 regulatory activity. A recent study showed that unlike Rap1, whose association with RP gene promoters is not affected by environmental changes causing RP gene repression/reactivation, Abf1 association with both RP gene and ribosome biogenesis (Ribi) gene promoters dynamically responds to changes in growth conditions. This observation changes the paradigm of general regulatory factors as relatively static DNA-binding proteins constitutively bound to highly active promoters, and point to Abf1, which binds hundreds of non-RPG promoters within the yeast genome, as a possible key regulatory switch in nutrient- and stress-dependent transcriptional modulation. Moreover, the frequent presence of Abf1 binding sites in the promoters of mitochondrial RP genes evokes the possibility that Abf1 might orchestrate still unexplored levels of co-regulation involving growth-related gene networks in yeast cells.

Promoter architecture and transcriptional regulation of Abf1-dependent ribosomal protein genes in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, ribosomal protein gene (RPG) promoters display binding sites for either Rap1 or Abf1 transcription factors. Unlike Rap1-associated promoters, the small cohort of Abf1-dependent RPGs (Abf1-RPGs) has not been extensively investigated. We show that RPL3, RPL4B, RPP1A, RPS22B and RPS28A/B share a common promoter architecture, with an Abf1 site upstream of a conserved element matching the sequence recognized by Fhl1, a transcription factor which together with Ifh1 orchestrates Rap1-associated RPG regulation. Abf1 and Fhl1 promoter association was confirmed by ChIP and/or gel retardation assays. Mutational analysis revealed a more severe requirement of Abf1 than Fhl1 binding sites for RPG transcription. In the case of RPS22B an unusual Tbf1 binding site promoted both RPS22B and intron-hosted SNR44 expression. Abf1-RPG down-regulation upon TOR pathway inhibition was much attenuated at defective mutant promoters unable to bind Abf1. TORC1 inactivation caused the expected reduction of Ifh1 occupancy at RPS22B and RPL3 promoters, but unexpectedly it entailed largely increased Abf1 association with Abf1-RPG promoters. We present evidence that Abf1 recruitment upon nutritional stress, also observed for representative ribosome biogenesis genes, favours RPG transcriptional rescue upon nutrient replenishment, thus pointing to nutrient-regulated Abf1 dynamics at promoters as a novel mechanism in ribosome biogenesis control.

Linking replication stress with heterochromatin formation

The eukaryotic genome can be roughly divided into euchromatin and heterochromatin domains that are structurally and functionally distinct. Heterochromatin is characterized by its high compaction that impedes DNA transactions such as gene transcription, replication, or recombination. Beyond its role in regulating DNA accessibility, heterochromatin plays essential roles in nuclear architecture, chromosome segregation, and genome stability. The formation of heterochromatin involves special histone modifications and the recruitment and spreading of silencing complexes that impact the higher-order structures of chromatin; however, its molecular nature varies between different chromosomal regions and between species. Although heterochromatin has been extensively characterized, its formation and maintenance throughout the cell cycle are not yet fully understood. The biggest challenge for the faithful transmission of chromatin domains is the destabilization of chromatin structures followed by their reassembly on a novel DNA template during genomic replication. This destabilizing event also provides a window of opportunity for the de novo establishment of heterochromatin. In recent years, it has become clear that different types of obstacles such as tight protein-DNA complexes, highly transcribed genes, and secondary DNA structures could impede the normal progression of the replisome and thus have the potential to endanger the integrity of the genome. Multiple studies carried out in different model organisms have demonstrated the capacity of such replisome impediments to favor the formation of heterochromatin. Our review summarizes these reports and discusses the potential role of replication stress in the formation and maintenance of heterochromatin and the role that silencing proteins could play at sites where the integrity of the genome is compromised.

Telomere DNA recognition in Saccharomycotina yeast: potential lessons for the co-evolution of ssDNA and dsDNA-binding proteins and their target sites

In principle, alterations in the telomere repeat sequence would be expected to disrupt the protective nucleoprotein complexes that confer stability to chromosome ends, and hence relatively rare events in evolution. Indeed, numerous organisms in diverse phyla share a canonical 6 bp telomere repeat unit (5'-TTAGGG-3'/5'-CCCTAA-3'), suggesting common descent from an ancestor that carries this particular repeat. All the more remarkable, then, are the extraordinarily divergent telomere sequences that populate the Saccharomycotina subphylum of budding yeast. These sequences are distinguished from the canonical telomere repeat in being long, occasionally degenerate, and frequently non-G/C-rich. Despite the divergent telomere repeat sequences, studies to date indicate that the same families of single-strand and double-strand telomere binding proteins (i.e., the Cdc13 and Rap1 families) are responsible for telomere protection in Saccharomycotina yeast. The recognition mechanisms of the protein family members therefore offer an informative paradigm for understanding the co-evolution of DNA-binding proteins and the cognate target sequences. Existing data suggest three potential, inter-related solutions to the DNA recognition problem: (i) duplication of the recognition protein and functional modification; (ii) combinatorial recognition of target site; and (iii) flexibility of the recognition surfaces of the DNA-binding proteins to adopt alternative conformations. Evidence in support of these solutions and the relevance of these solutions to other DNA-protein regulatory systems are discussed.

Chromatin mediation of a transcriptional memory effect in yeast

Previous studies have described a transcriptional "memory effect," whereby transcript levels of many Abf1-regulated genes in the budding yeast Saccharomyces cerevisiae are undiminished even after Abf1 has dissociated from its regulatory sites. Here we provide additional support for this effect and investigate its molecular basis. We show that the effect is observed in a distinct abf1 ts mutant from that used in earlier studies, demonstrating that it is robust, and use chromatin immunoprecipitation to show that Abf1 association is decreased similarly from memory effect and transcriptionally responsive genes at the restrictive temperature. We also demonstrate that the association of TATA-binding protein and Pol II decreases after the loss of Abf1 binding for transcriptionally responsive genes but not for memory effect genes. Examination of genome-wide nucleosome occupancy data reveals that although transcriptionally responsive genes exhibit increased nucleosome occupancy in abf1 ts yeast, the promoter regions of memory effect targets show no change in abf1 ts mutants, maintaining an open chromatin conformation even after Abf1 eviction. This contrasting behavior reflects different inherent propensity for nucleosome formation between the two classes, driven by the presence of A/T-rich sequences upstream of the Abf1 site in memory effect gene promoters. These sequence-based differences show conservation in closely related fungi and also correlate with different gene expression noise, suggesting a physiological basis for greater access to "memory effect" promoter regions. Thus, our results establish a conserved mechanism underlying a transcriptional memory effect whereby sequences surrounding Abf1 binding sequences affect local nucleosome occupancy following loss of Abf1 binding. Furthermore, these findings demonstrate that sequence-based differences in the propensity for nucleosome occupancy can influence the transcriptional response of genes to an altered regulatory signal.

GENETIC MUTATIONS ARE CHARACTERIZED BY INCREASE IN ENTROPY AT THE TRANSCRIPTIONAL LEVEL

Predicting the genomic and phenotypic re-programming in organisms undergoing genetic perturbations is a challenging task in modern biology. It is hypothesized that genomic alterations perturb the dynamics of biological information flow. In the present study, a statistical data analysis framework was designed and the network entropy concept was employed to quantify the level of disorder at the transcriptional level as a result of the genomic re-programming of S. cerevisiae cells under genetic perturbations. The customized re-programming in transcription levels to different genetic modifications was observed and genetic mutations were characterized by enhanced network entropies, which revealed higher degree of randomness in mRNA expression levels. To our knowledge, this study constitutes the first numerical demonstration on the conservative energetic state of the microorganisms against genetic perturbations.

Histone hypoacetylation-activated genes are repressed by acetyl-CoA- and chromatin-mediated mechanism

Transcriptional activation is typically associated with increased acetylation of promoter histones. However, this paradigm does not apply to transcriptional activation of all genes. In this study we have characterized a group of genes that are repressed by histone acetylation. These histone hypoacetylation-activated genes (HHAAG) are normally repressed during exponential growth, when the cellular level of acetyl-CoA is high and global histone acetylation is also high. The HHAAG are induced during diauxic shift, when the levels of acetyl-CoA and global histone acetylation decrease. The histone hypoacetylation-induced activation of HHAAG is independent of Msn2/Msn4. The repression of HSP12, one of the HHAAG, is associated with well-defined nucleosomal structure in the promoter region, while histone hypoacetylation-induced activation correlates with delocalization of positioned nucleosomes or with reduced nucleosome occupancy. Correspondingly, unlike the majority of yeast genes, HHAAG are transcriptionally upregulated when expression of histone genes is reduced. Taken together, these results suggest a model in which histone acetylation is required for proper positioning of promoter nucleosomes and repression of HHAAG.

Silencing subtelomeric VSGs by Trypanosoma brucei RAP1 at the insect stage involves chromatin structure changes

Trypanosoma brucei causes human African trypanosomiasis and regularly switches its major surface antigen variant surface glycoprotein (VSG) to evade mammalian host immune responses at the bloodstream form (BF) stage. Monoallelic expression of BF Expression Site (BES)-linked VSGs and silencing of metacyclic VSGs (mVSGs) in BF cells are essential for antigenic variation, whereas silencing of both BES-linked and mVSGs in the procyclic form (PF) cells is important for cell survival in the midgut of its insect vector. We have previously shown that silencing BES-linked VSGs in BF cells depends on TbRAP1. We now show that TbRAP1 silences both BES-linked and mVSGs at both BF and PF stages. The strength of TbRAP1-mediated BES-linked VSG silencing is stronger in the PF cells than that in BF cells. In addition, Formaldehyde-Assisted Isolation of Regulatory Elements analysis and MNase digestion demonstrated that depletion of TbRAP1 in PF cells led to a chromatin structure change, which is significantly stronger at the subtelomeric VSG loci than at chromosome internal loci. On the contrary, no significant chromatin structure changes were detected on depletion of TbRAP1 in BF cells. Our observations indicate that TbRAP1 helps to determine the chromatin structure at the insect stage, which likely contributes to its strong silencing effect on VSGs.

Nucleotide excision repair in cellular chromatin: studies with yeast from nucleotide to gene to genome

Here we review our development of, and results with, high resolution studies on global genome nucleotide excision repair (GGNER) in Saccharomyces cerevisiae. We have focused on how GGNER relates to histone acetylation for its functioning and we have identified the histone acetyl tranferase Gcn5 and acetylation at lysines 9/14 of histone H3 as a major factor in enabling efficient repair. We consider results employing primarily MFA2 as a model gene, but also those with URA3 located at subtelomeric sequences. In the latter case we also see a role for acetylation at histone H4. We then go on to outline the development of a high resolution genome-wide approach that enables one to examine correlations between histone modifications and the nucleotide excision repair (NER) of UV-induced cyclobutane pyrimidine dimers throughout entire genomes. This is an approach that will enable rapid advances in understanding the complexities of how compacted chromatin in chromosomes is processed to access DNA damage and then returned to its pre-damaged status to maintain epigenetic codes.

Chromatin and transcription in yeast

Understanding the mechanisms by which chromatin structure controls eukaryotic transcription has been an intense area of investigation for the past 25 years. Many of the key discoveries that created the foundation for this field came from studies of Saccharomyces cerevisiae, including the discovery of the role of chromatin in transcriptional silencing, as well as the discovery of chromatin-remodeling factors and histone modification activities. Since that time, studies in yeast have continued to contribute in leading ways. This review article summarizes the large body of yeast studies in this field.

Functions of protosilencers in the formation and maintenance of heterochromatin in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, transcriptionally silent heterochromatin at HML and HMR loci is established by silencers that recruit SIR complex and promote its propagation along chromatin. Silencers consist of various combinations of two or three binding sites for origin recognition complex (ORC), Abf1 and Rap1. A single ORC, Abf1 or Rap1 site cannot promote silencing, but can enhance silencing by a distant silencer, and is called a protosilencer. The mechanism of protosilencer function is not known. We examine the functions of ORC, Abf1 and Rap1 sites as components of the HMR-E silencer, and as protosilencers. We find that the Rap1 site makes a larger and unique contribution to HMR-E function compared to ORC and Abf1 sites. On the other hand, Rap1 site does not act as a protosilencer to assist HML-E silencer in forming heterochromatin, whereas ORC and Abf1 sites do. Therefore, different mechanisms may be involved in the roles of Rap1 site as a component of HMR-E and as a protosilencer. Heterochromatin formed by ORC or Abf1 site in collaboration with HML-E is not as stable as that formed by HMR-E and HML-E, but increasing the copy number of Abf1 site enhances heterochromatin stability. ORC and Abf1 sites acting as protosilencers do not modulate chromatin structure in the absence of SIR complex, which argues against the hypothesis that protosilencers serve to create a chromatin structure favorable for SIR complex propagation. We also investigate the function of ARS1 containing an ORC site and an Abf1 site as a protosilencer. We find that ARS1 inserted at HML enhances heterochromatin stability, and promotes de novo formation of a chromatin structure that partially resembles heterochromatin in an S phase dependent manner. Taken together, our results indicate that protosilencers aid in the formation and maintenance of heterochromatin structure.

The orientation of the C-terminal domain of the Saccharomyces cerevisiae Rap1 protein is determined by its binding to DNA

Rap1 is an essential DNA-binding factor from the yeast Saccharomyces cerevisiae involved in transcription and telomere maintenance. Its binding to DNA targets Rap1 at particular loci, and may optimize its ability to form functional macromolecular assemblies. It is a modular protein, rich in large potentially unfolded regions, and comprising BRCT, Myb and RCT well-structured domains. Here, we present the architectures of Rap1 and a Rap1/DNA complex, built through a step-by-step integration of small angle X-ray scattering, X-ray crystallography and nuclear magnetic resonance data. Our results reveal Rap1 structural adjustment upon DNA binding that involves a specific orientation of the C-terminal (RCT) domain with regard to the DNA binding domain (DBD). Crystal structure of DBD in complex with a long DNA identifies an essential wrapping loop, which constrains the orientation of the RCT and affects Rap1 affinity to DNA. Based on our structural information, we propose a model for Rap1 assembly at telomere.

A MULTI-STRATEGY APPROACH TO INFORMATIVE GENE IDENTIFICATION FROM GENE EXPRESSION DATA

An unsupervised multi-strategy approach has been developed to identify informative genes from high throughput genomic data. Several statistical methods have been used in the field to identify differentially expressed genes. Since different methods generate different lists of genes, it is very challenging to determine the most reliable gene list and the appropriate method. This paper presents a multi-strategy method, in which a combination of several data analysis techniques are applied to a given dataset and a confidence measure is established to select genes from the gene lists generated by these techniques to form the core of our final selection. The remainder of the genes that form the peripheral region are subject to exclusion or inclusion into the final selection. This paper demonstrates this methodology through its application to an in-house cancer genomics dataset and a public dataset. The results indicate that our method provides more reliable list of genes, which are validated using biological knowledge, biological experiments, and literature search. We further evaluated our multi-strategy method by consolidating two pairs of independent datasets, each pair is for the same disease, but generated by different labs using different platforms. The results showed that our method has produced far better results.

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.

Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization

The packaging of eukaryotic genomes into nuclesomes plays critical roles in chromatin organization and gene regulation. Studies in Saccharomyces cerevisiae indicate that nucleosome occupancy is partially encoded by intrinsic antinucleosomal DNA sequences, such as poly(A) sequences, as well as by binding sites for trans-acting factors that can evict nucleosomes, such as Reb1 and the Rsc3/30 complex. Here, we use genome-wide nucleosome occupancy maps in 13 Ascomycota fungi to discover large-scale evolutionary reprogramming of both intrinsic and trans determinants of chromatin structure. We find that poly(G)s act as intrinsic antinucleosomal sequences, comparable to the known function of poly(A)s, but that the abundance of poly(G)s has diverged greatly between species, obscuring their antinucleosomal effect in low-poly(G) species such as S. cerevisiae. We also develop a computational method that uses nucleosome occupancy maps for discovering trans-acting general regulatory factor (GRF) binding sites. Our approach reveals that the specific sequences bound by GRFs have diverged substantially across evolution, corresponding to a number of major evolutionary transitions in the repertoire of GRFs. We experimentally validate a proposed evolutionary transition from Cbf1 as a major GRF in pre-whole-genome duplication (WGD) yeasts to Reb1 in post-WGD yeasts. We further show that the mating type switch-activating protein Sap1 is a GRF in S. pombe, demonstrating the general applicability of our approach. Our results reveal that the underlying mechanisms that determine in vivo chromatin organization have diverged and that comparative genomics can help discover new determinants of chromatin organization.

GAS3, a developmentally regulated gene, encodes a highly mannosylated and inactive protein of the Gas family of Saccharomyces cerevisiae

The multigene GAS family of Saccharomyces cerevisiae is constituted by five genes encoding GPI-anchored proteins required for cell wall or spore wall assembly. GAS1 and GAS5 are expressed in vegetative growth and repressed during sporulation, whereas GAS2 and GAS4 exhibit the opposite expression pattern. This study focuses on GAS3, a still poorly characterized member of the family. To date, attempts to reveal the glucan elongase activity typical of Gas proteins have been unsuccessful, suggesting that Gas3p is the only inactive member of the family. Here, we compared the mRNA levels of GAS1, GAS3 and GAS5 and demonstrate that GAS3 is the weakest-expressed paralogue in vegetative growth. Moreover, GAS3 mRNA increased during sporulation, showing a bimodal profile typical of the early-middle meiotic genes. GAS3 product was identified as a low-abundance, polydisperse mannoprotein. Loss of Gas3p did not affect growth and sporulation. The overexpression of GAS3, driven by the GAS1 promoter, slightly reduced growth rate in a wild-type strain and led to hyperaccumulation of Gas3p in the membranes and in the cell wall. To determine whether GAS3 could replace GAS1 function in vivo, GAS3 was also overexpressed in a gas1Delta mutant. Increased amounts of Gas3p were not only unable to complement the defects of the gas1Delta cells but exacerbated them. A mutated Gas3p-E283Q, where one of the catalytic glutamate residues essential for GH72 enzyme activity was replaced by glutamine, was also noxious to gas1Delta cells, indicating that the increased expression of Gas3p, rather than a potential activity, is deleterious for gas1Delta cells.

Extensive role of the general regulatory factors, Abf1 and Rap1, in determining genome-wide chromatin structure in budding yeast

The packaging of eukaryotic DNA into chromatin has profound consequences for gene regulation, as well as for other DNA transactions such as recombination, replication and repair. Understanding how this packaging is determined is consequently a pressing problem in molecular genetics. DNA sequence, chromatin remodelers and transcription factors affect chromatin structure, but the scope of these influences on genome-wide nucleosome occupancy patterns remains uncertain. Here, we use high resolution tiling arrays to examine the contributions of two general regulatory factors, Abf1 and Rap1, to nucleosome occupancy in Saccharomyces cerevisiae. These factors have each been shown to bind to a few hundred promoters, but we find here that thousands of loci show localized regions of altered nucleosome occupancy within 1 h of loss of Abf1 or Rap1 binding, and that altered chromatin structure can occur via binding sites having a wide range of affinities. These results indicate that DNA-binding transcription factors affect chromatin structure, and probably dynamics, throughout the genome to a much greater extent than previously appreciated.

The role of nucleosome positioning in the evolution of gene regulation

Chromatin organization plays a major role in gene regulation and can affect the function and evolution of new transcriptional programs. However, it can be difficult to decipher the basis of changes in chromatin organization and their functional effect on gene expression. Here, we present a large-scale comparative genomic analysis of the relationship between chromatin organization and gene expression, by measuring mRNA abundance and nucleosome positions genome-wide in 12 Hemiascomycota yeast species. We found substantial conservation of global and functional chromatin organization in all species, including prominent nucleosome-free regions (NFRs) at gene promoters, and distinct chromatin architecture in growth and stress genes. Chromatin organization has also substantially diverged in both global quantitative features, such as spacing between adjacent nucleosomes, and in functional groups of genes. Expression levels, intrinsic anti-nucleosomal sequences, and trans-acting chromatin modifiers all play important, complementary, and evolvable roles in determining NFRs. We identify five mechanisms that couple chromatin organization to evolution of gene regulation and have contributed to the evolution of respiro-fermentation and other key systems, including (1) compensatory evolution of alternative modifiers associated with conserved chromatin organization, (2) a gradual transition from constitutive to trans-regulated NFRs, (3) a loss of intrinsic anti-nucleosomal sequences accompanying changes in chromatin organization and gene expression, (4) re-positioning of motifs from NFRs to nucleosome-occluded regions, and (5) the expanded use of NFRs by paralogous activator-repressor pairs. Our study sheds light on the molecular basis of chromatin organization, and on the role of chromatin organization in the evolution of gene regulation.

TOPTMH: topology predictor for transmembrane alpha-helices

Alpha-helical transmembrane proteins mediate many key biological processes and represent 20%-30% of all genes in many organisms. Due to the difficulties in experimentally determining their high-resolution 3D structure, computational methods to predict the location and orientation of transmembrane helix segments using sequence information are essential. We present TOPTMH, a new transmembrane helix topology prediction method that combines support vector machines, hidden Markov models, and a widely used rule-based scheme. The contribution of this work is the development of a prediction approach that first uses a binary SVM classifier to predict the helix residues and then it employs a pair of HMM models that incorporate the SVM predictions and hydropathy-based features to identify the entire transmembrane helix segments by capturing the structural characteristics of these proteins. TOPTMH outperforms state-of-the-art prediction methods and achieves the best performance on an independent static benchmark.

Uncovering transcriptional interactions via an adaptive fuzzy logic approach

BACKGROUND

To date, only a limited number of transcriptional regulatory interactions have been uncovered. In a pilot study integrating sequence data with microarray data, a position weight matrix (PWM) performed poorly in inferring transcriptional interactions (TIs), which represent physical interactions between transcription factors (TF) and upstream sequences of target genes. Inferring a TI means that the promoter sequence of a target is inferred to match the consensus sequence motifs of a potential TF, and their interaction type such as AT or RT is also predicted. Thus, a robust PWM (rPWM) was developed to search for consensus sequence motifs. In addition to rPWM, one feature extracted from ChIP-chip data was incorporated to identify potential TIs under specific conditions. An interaction type classifier was assembled to predict activation/repression of potential TIs using microarray data. This approach, combining an adaptive (learning) fuzzy inference system and an interaction type classifier to predict transcriptional regulatory networks, was named AdaFuzzy.

RESULTS

AdaFuzzy was applied to predict TIs using real genomics data from Saccharomyces cerevisiae. Following one of the latest advances in predicting TIs, constrained probabilistic sparse matrix factorization (cPSMF), and using 19 transcription factors (TFs), we compared AdaFuzzy to four well-known approaches using over-representation analysis and gene set enrichment analysis. AdaFuzzy outperformed these four algorithms. Furthermore, AdaFuzzy was shown to perform comparably to 'ChIP-experimental method' in inferring TIs identified by two sets of large scale ChIP-chip data, respectively. AdaFuzzy was also able to classify all predicted TIs into one or more of the four promoter architectures. The results coincided with known promoter architectures in yeast and provided insights into transcriptional regulatory mechanisms.

CONCLUSION

AdaFuzzy successfully integrates multiple types of data (sequence, ChIP, and microarray) to predict transcriptional regulatory networks. The validated success in the prediction results implies that AdaFuzzy can be applied to uncover TIs in yeast.

Cis-regulatory elements in plant cell signaling

Plant cell signaling pathways are in part dependent on transcriptional regulatory networks comprising circuits of transcription factors (TFs) and regulatory DNA elements that control the expression of target genes. Here, we describe experimental and bioinformatic approaches for identifying potential cis-regulatory elements. We also discuss recent integrative genomics studies aimed at elucidating the functions of cis-regulatory elements in aspects of plant biology, including the circadian clock, interactions with the environment, stress responses, and regulation of growth and development by phytohormones. Finally, we discuss emerging technologies and approaches that offer great potential for accelerating the discovery and functional characterization of cis-elements and interacting TFs--which will help realize the promise of systems biology.

Mediator complex association with constitutively transcribed genes in yeast

Mediator is a large, multisubunit complex that is essential for transcription of mRNA by RNA Pol II in eukaryotes and is believed to bridge transcriptional activators and the general transcription machinery. However, several recent studies suggest that the requirement for Mediator during transcriptional activation is not universal, but rather activator dependent, and may be indirect for some genes. Here we have investigated Mediator association with several constitutively transcribed genes in yeast by comparing a yeast strain that harbors a temperature-sensitive mutation in an essential Mediator subunit, Srb4, with its wild-type (WT) counterpart. We find modest association of Mediator with constitutively active genes and show that this association is strongly decreased in srb4 ts yeast, whereas association with a nontranscribed region or repressed gene promoters is lower and unaffected in the mutant yeast. The tail module of Mediator remains associated with ribosomal protein (RP) gene promoters in srb4 ts yeast, while subunits from the head and middle modules are lost. Tail module association at Rap1-dependent gene promoters is lost in rap1 ts yeast, indicating that Rap1 is required for Mediator recruitment at these gene promoters and that its recruitment occurs via the tail module. Pol II association is also rapidly and severely affected in srb4 ts yeast, indicating that Mediator is directly required for pol II association at constitutively transcribed genes. Our results are consistent with Mediator functioning as a general transcription factor in yeast.

Distinct differences in chromatin structure at subtelomeric X and Y' elements in budding yeast

In Saccharomyces cerevisiae, all ends of telomeric DNA contain telomeric repeats of (TG(1-3)), but the number and position of subtelomeric X and Y' repeat elements vary. Using chromatin immunoprecipitation and genome-wide analyses, we here demonstrate that the subtelomeric X and Y' elements have distinct structural and functional properties. Y' elements are transcriptionally active and highly enriched in nucleosomes, whereas X elements are repressed and devoid of nucleosomes. In contrast to X elements, the Y' elements also lack the classical hallmarks of heterochromatin, such as high Sir3 and Rap1 occupancy as well as low levels of histone H4 lysine 16 acetylation. Our analyses suggest that the presence of X and Y' elements govern chromatin structure and transcription activity at individual chromosome ends.

VisHiC--hierarchical functional enrichment analysis of microarray data

Measuring gene expression levels with microarrays is one of the key technologies of modern genomics. Clustering of microarray data is an important application, as genes with similar expression profiles may be regulated by common pathways and involved in related functions. Gene Ontology (GO) analysis and visualization allows researchers to study the biological context of discovered clusters and characterize genes with previously unknown functions. We present VisHiC (Visualization of Hierarchical Clustering), a web server for clustering and compact visualization of gene expression data combined with automated function enrichment analysis. The main output of the analysis is a dendrogram and visual heatmap of the expression matrix that highlights biologically relevant clusters based on enriched GO terms, pathways and regulatory motifs. Clusters with most significant enrichments are contracted in the final visualization, while less relevant parts are hidden altogether. Such a dense representation of microarray data gives a quick global overview of thousands of transcripts in many conditions and provides a good starting point for further analysis. VisHiC is freely available at http://biit.cs.ut.ee/vishic.

Progress and challenges in predicting protein-protein interaction sites

The identification of protein-protein interaction sites is an essential intermediate step for mutant design and the prediction of protein networks. In recent years a significant number of methods have been developed to predict these interface residues and here we review the current status of the field. Progress in this area requires a clear view of the methodology applied, the data sets used for training and testing the systems, and the evaluation procedures. We have analysed the impact of a representative set of features and algorithms and highlighted the problems inherent in generating reliable protein data sets and in the posterior analysis of the results. Although it is clear that there have been some improvements in methods for predicting interacting sites, several major bottlenecks remain. Proteins in complexes are still under-represented in the structural databases and in particular many proteins involved in transient complexes are still to be crystallized. We provide suggestions for effective feature selection, and make it clear that community standards for testing, training and performance measures are necessary for progress in the field.

Using triplet periodicity of nucleotide sequences for finding potential reading frame shifts in genes

We introduce a novel approach for the detection of possible mutations leading to a reading frame (RF) shift in a gene. Deletions and insertions of DNA coding regions are considerable events for genes because an RF shift results in modifications of the extensive region of amino acid sequence coded by a gene. The suggested method is based on the phenomenon of triplet periodicity (TP) in coding regions of genes and its relative resistance to substitutions in DNA sequence. We attempted to extend 326 933 regions of continuous TP found in genes from the KEGG databank by considering possible insertions and deletions. We revealed totally 824 genes where such extension was possible and statistically significant. Then we generated amino acid sequences according to active (KEGG's) and hypothetically ancient RFs in order to find confirmation of a shift at a protein level. Consequently, 64 sequences have protein similarities only for ancient RF, 176 only for active RF, 3 for both and 581 have no protein similarity at all. We aimed to have revealed lower bound for the number of genes in which a shift between RF and TP is possible. Further ways to increase the number of revealed RF shifts are discussed.

Asymmetric relationships between proteins shape genome evolution

An investigation of metabolic networks in E. coli and S. cerevisiae reveals that asymmetric protein interactions affect gene expression, the relative effect of gene-knockouts and genome evolution.

Background

The relationships between proteins are often asymmetric: one protein (A) depends for its function on another protein (B), but the second protein does not depend on the first. In metabolic networks there are multiple pathways that converge into one central pathway. The enzymes in the converging pathways depend on the enzymes in the central pathway, but the enzymes in the latter do not depend on any specific enzyme in the converging pathways. Asymmetric relations are analogous to the “if->then” logical relation where A implies B, but B does not imply A (A->B).

Results

We show that the majority of relationships between enzymes in metabolic flux models of metabolism in Escherichia coli and Saccharomyces cerevisiae are asymmetric. We show furthermore that these asymmetric relationships are reflected in the expression of the genes encoding those enzymes, the effect of gene knockouts and the evolution of genomes. From the asymmetric relative dependency, one would expect that the gene that is relatively independent (B) can occur without the other dependent gene (A), but not the reverse. Indeed, when only one gene of an A->B pair is expressed, is essential, is present in a genome after an evolutionary gain or loss, it tends to be the independent gene (B). This bias is strongest for genes encoding proteins whose asymmetric relationship is evolutionarily conserved.

Conclusions

The asymmetric relations between proteins that arise from the system properties of metabolic networks affect gene expression, the relative effect of gene knockouts and genome evolution in a predictable manner.