Statistical analysis of sequence features of introns with positive transcriptional regulation in yeast genes

A great deal of experimental studies have shown that many introns of eukaryotic genes function as regulators of transcription. However, comprehensive studies of this problem have not yet been conducted. After checking the transcription frequencies of some Saccharomyces cerevisiae (yeast), genes and their introns, a remarkable phenomenon was discovered that generally the introns of the genes with higher transcription frequencies are longer, and the introns of the genes with lower transcription frequencies are shorter. This suggests that the longer introns of genes with higher transcription frequencies may contain some characteristic sequence structures, which could enhance the transcription of genes. Therefore, two sets of introns of yeast genes were chosen for further study. The transcription frequencies of the first set of genes are higher (>30), and those of the second set of genes are lower (less than or equal to10). Some oligonucleotides are detected by statistically comparative analyses of the occurrence frequencies of oligonucleotides (mainly tetranucleotides and pentanucleotides), whose occurrence frequencies in the first set of introns; are significantly higher than those in the second set of introns, and are also significantly higher than those in the exons flanking the introns of the first set. Some of these extracted oligonucleotides are the same as the regulatory elements of transcription revealed by experimental analyses. Besides, the distributions of these extracted oligonucleotides in the two sets of introns and the exons show that the sequence structures of the first set of introns are favorable for transcription of genes.

Uncovering the rules for protein-protein interactions from yeast genomic data

Identifying protein-protein interactions is crucial for understanding cellular functions. Genomic data provides opportunities and challenges in identifying these interactions. We uncover the rules for predicting protein-protein interactions using a frequent pattern tree (FPT) approach modified to generate a minimum set of rules (mFPT), with rule attributes constructed from the interaction features of the yeast genomic data. The mFPT prediction accuracy is benchmarked against other commonly used methods such as Bayesian networks and logistic regressions under various statistical measures. Our study indicates that mFPT outranks other methods in predicting the protein-protein interactions for the database used. We predict a new protein-protein interaction complex whose biological function is related to premRNA splicing and new protein-protein interactions within existing complexes based on the rules generated.

Cflec-5, a pattern recognition receptor in scallop Chlamys farreri agglutinating yeast Pichia pastoris

C-type lectins are a superfamily of carbohydrate-recognition proteins which play crucial roles as pattern recognition receptors (PRRs) in the innate immunity. In this study, the full-length cDNA of a C-type lectin was cloned from scallop Chlamys farreri (designated as Cflec-5) by expression sequence tag (EST) analysis and rapid amplification of cDNA ends (RACE) approach The full-length cDNA of Cflec-5 was of 1412 bp. The open reading frame encoded a polypeptide of 153 amino acids, including a signal sequence and a conserved carbohydrate-recognition domain with the EPN motif determining the mannose-binding specificity The deduced amino acid sequence of Cflec-5 showed high similarity to members of C-type lectin superfamily. The quantitative real-time PCR was performed to investigate the tissue distribution of Cflec-5 mRNA and its temporal expression profiles in hemocytes post pathogen-associated molecular patterns (PAMPs) stimulation. In healthy scallops, the Cflec-5 mRNA was mainly detected in gill and mantle, and marginally in other tissues The mRNA expression of Cflec-5 could be significantly induced by lipopolysaccharide (LPS) and glucan stimulation and reached the maximum level at 6 h and 12 h, respectively But its expression level did not change significantly during peptidoglycan (PGN) stimulation The function of Cflec-5 was investigated by recombination and expression of the cDNA fragment encoding its mature peptide in Escherichia coli Rosetta Gami (DE3) The recombinant Cflec-5 agglutinated Pichia pastoris in a calcium-independent way The agglutinating activity could be inhibited by D-mannose. LPS and glucan, but not by D-galactose or PGN. These results collectively suggested that Cflec-5 was involved in the innate Immune response of scallops and might contribute to nonself-recognition through its interaction with various PAMPs (C) 2010 Elsevier Ltd All rights reserved

Predicting gene function in Saccharomyces cerevisiae

Clare, A. and King R.D. (2003) Predicting gene function in Saccharomyces cerevisiae. 2nd European Conference on Computational Biology (ECCB '03). (published as a journal supplement in Bioinformatics 19: ii42-ii49)

Integration of relational and hierarchical network information for protein function prediction

BACKGROUND:In the current climate of high-throughput computational biology, the inference of a protein's function from related measurements, such as protein-protein interaction relations, has become a canonical task. Most existing technologies pursue this task as a classification problem, on a term-by-term basis, for each term in a database, such as the Gene Ontology (GO) database, a popular rigorous vocabulary for biological functions. However, ontology structures are essentially hierarchies, with certain top to bottom annotation rules which protein function predictions should in principle follow. Currently, the most common approach to imposing these hierarchical constraints on network-based classifiers is through the use of transitive closure to predictions.RESULTS:We propose a probabilistic framework to integrate information in relational data, in the form of a protein-protein interaction network, and a hierarchically structured database of terms, in the form of the GO database, for the purpose of protein function prediction. At the heart of our framework is a factorization of local neighborhood information in the protein-protein interaction network across successive ancestral terms in the GO hierarchy. We introduce a classifier within this framework, with computationally efficient implementation, that produces GO-term predictions that naturally obey a hierarchical 'true-path' consistency from root to leaves, without the need for further post-processing.CONCLUSION:A cross-validation study, using data from the yeast Saccharomyces cerevisiae, shows our method offers substantial improvements over both standard 'guilt-by-association' (i.e., Nearest-Neighbor) and more refined Markov random field methods, whether in their original form or when post-processed to artificially impose 'true-path' consistency. Further analysis of the results indicates that these improvements are associated with increased predictive capabilities (i.e., increased positive predictive value), and that this increase is consistent uniformly with GO-term depth. Additional in silico validation on a collection of new annotations recently added to GO confirms the advantages suggested by the cross-validation study. Taken as a whole, our results show that a hierarchical approach to network-based protein function prediction, that exploits the ontological structure of protein annotation databases in a principled manner, can offer substantial advantages over the successive application of 'flat' network-based methods.

Helicobacter pylori: comparative genomics and structure-function analysis of the flagellum biogenesis protein HP0958

Helicobacter pylori is a gastric pathogen which infects ~50% of the global population and can lead to the development of gastritis, gastric and duodenal ulcers and carcinoma. Genome sequencing of H. pylori revealed high levels of genetic variability; this pathogen is known for its adaptability due to mechanisms including phase variation, recombination and horizontal gene transfer. Motility is essential for efficient colonisation by H. pylori. The flagellum is a complex nanomachine which has been studied in detail in E. coli and Salmonella. In H. pylori, key differences have been identified in the regulation of flagellum biogenesis, warranting further investigation. In this study, the genomes of two H. pylori strains (CCUG 17874 and P79) were sequenced and published as draft genome sequences. Comparative studies identified the potential role of restriction modification systems and the comB locus in transformation efficiency differences between these strains. Core genome analysis of 43 H. pylori strains including 17874 and P79 defined a more refined core genome for the species than previously published. Comparative analysis of the genome sequences of strains isolated from individuals suffering from H. pylori related diseases resulted in the identification of “disease-specific” genes. Structure-function analysis of the essential motility protein HP0958 was performed to elucidate its role during flagellum assembly in H. pylori. The previously reported HP0958-FliH interaction could not be substantiated in this study and appears to be a false positive. Site-directed mutagenesis confirmed that the coiled-coil domain of HP0958 is involved in the interaction with RpoN (74-284), while the Zn-finger domain is required for direct interaction with the full length flaA mRNA transcript. Complementation of a non-motile hp0958-null derivative strain of P79 with site-directed mutant alleles of hp0958 resulted in cells producing flagellar-type extrusions from non-polar positions. Thus, HP0958 may have a novel function in spatial localisation of flagella in H. pylori

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

The Function and Regulation of Photobodies in Phytochrome Signaling

Light is a critical environmental signal that regulates every phase of the plant life cycle, from germination to floral initiation. Of the many light receptors in the model plant Arabidopsis thaliana, the red- and far-red light-sensing phytochromes (phys) are arguably the best studied, but the earliest events in the phy signaling pathway remain poorly understood. One of the earliest phy signaling events is the translocation of photoactivated phys from the cytoplasm to the nucleus, where they localize to subnuclear foci termed photobodies; in continuous light, photobody localization correlates closely with the light-dependent inhibition of embryonic stem growth. Despite a growing body of evidence supporting the biological significance of photobodies in light signaling, photobodies have also been shown to be dispensable for seedling growth inhibition in continuous light, so their physiological importance remains controversial; additionally, the molecular components that are required for phy localization to photobodies are largely unknown. The overall goal of my dissertation research was to gain insight into the early steps of phy signaling by further defining the role of photobodies in this process and identifying additional intragenic and extragenic requirements for phy localization to photobodies.

Even though the domain structure of phys has been extensively studied, not all of the intramolecular requirements for phy localization to photobodies are known. Previous studies have shown that the entire C-terminus of phys is both necessary and sufficient for their localization to photobodies. However, the importance of the individual subdomains of the C-terminus is still unclear. For example a truncation lacking part of the most C-terminal domain, the histidine kinase-related domain (HKRD), can still localize to small photobodies in the light and behaves like a weak allele. However, a point mutation within the HKRD renders the entire molecule completely inactive. To resolve this discrepancy, I explored the hypothesis that this point mutation might impair the dimerization of the HKRD; dimerization has been shown to occur via the C-terminus of phy and is required for more efficient signaling. I show that this point mutation impairs nuclear localization of phy as well as its subnuclear localization to photobodies. Additionally, yeast-two-hybrid analysis shows that the wild-type HKRD can homodimerize but that the HKRD containing the point mutation fails to dimerize with both itself and with wild-type HKRD. These results demonstrate that dimerization of the HKRD is required for both nuclear and photobody localization of phy.

Studies of seedlings grown in diurnal conditions show that photoactivated phy can persist into darkness to repress seedling growth; a seedling's growth rate is therefore fastest at the end of the night. To test the idea that photobodies could be involved in regulating seedling growth in the dark, I compared the growth of two transgenic Arabidopsis lines, one in which phy can localize to photobodies (PBG), and one in which it cannot (NGB). Despite these differences in photobody morphology, both lines are capable of transducing light signals and inhibiting seedling growth in continuous light. After the transition from red light to darkness, the PBG line was able to repress seedling growth, as well as the accumulation of the growth-promoting, light-labile transcription factor PHYTOCHROME INTERACTING FACTOR 3 (PIF3), for eighteen hours, and this correlated perfectly with the presence of photobodies. Reducing the amount of active phy by either reducing the light intensity or adding a phy-inactivating far-red pulse prior to darkness led to faster accumulation of PIF3 and earlier seedling growth. In contrast, the NGB line accumulated PIF3 even in the light, and seedling growth was only repressed for six hours; this behavior was similar in NGB regardless of the light treatment. These results suggest that photobodies are required for the degradation of PIF3 and for the prolonged stabilization of active phy in darkness. They also support the hypothesis that photobody localization of phys could serve as an instructive cue during the light-to-dark transition, thereby fine-tuning light-dependent responses in darkness.

In addition to determining an intragenic requirement for photobody localization and further exploring the significance of photobodies in phy signaling, I wanted to identify extragenic regulators of photobody localization. A recent study identified one such factor, HEMERA (HMR); hmr mutants do not form large photobodies, and they are tall and albino in the light. To identify other components in the HMR-mediated branch of the phy signaling pathway, I performed a forward genetic screen for suppressors of a weak hmr allele. Surprisingly, the first three mutants isolated from the screen were alleles of the same novel gene, SON OF HEMERA (SOH). The soh mutations rescue all of the phenotypes associated with the weak hmr allele, and they do so in an allele-specific manner, suggesting a direct interaction between SOH and HMR. Null soh alleles, which were isolated in an independent, tall, albino screen, are defective in photobody localization, demonstrating that SOH is an extragenic regulator of phy localization to photobodies that works in the same genetic pathway as HMR.

In this work, I show that dimerization of the HKRD is required for both the nuclear and photobody localization of phy. I also demonstrate a tight correlation between photobody localization and PIF3 degradation, further establishing the significance of photobodies in phy signaling. Finally, I identify a novel gene, SON OF HEMERA, whose product is necessary for phy localization to photobodies in the light, thereby isolating a new extragenic determinant of photobody localization. These results are among the first to focus exclusively on one of the earliest cellular responses to light - photobody localization of phys - and they promise to open up new avenues into the study of a poorly understood facet of the phy signaling pathway.

Arp2p, a potential substrate for NatB acetyltransferase in fission yeast, is important in actin patch nucleation and motility in arp2-1 and arm1 ts mutants

The actin cytoskeleton is a dynamic and complex structure in fission yeast that plays a major function in many cell processes including cellular growth, septa formation, endocytosis and cellular division. Computational studies have shown that Arp2p, which forms part of the Arp2/3 complex, is a potential substrate of NatB acetyltransferase which has specificity for proteins possessing an N-terminal Met-Asp or Met-Glu sequence motif. In arm1- mutants the loss of function of Arm1p, an auxillary subunit required for NatB activity, results in a temperature sensitive phenotype characterized by multiple septa, failure of endocytosis, and the inability to form actin cables. A temperature sensitive mutant of Schizosaccharomyces pombe arp2 gene exhibits a similar phenotype as seen by the formation of improper septa, slow growth, and the delocalization of actin patches. Four expression vectors encoding the open reading frames of arp2 and cdc8 (tropomyosin) were constructed with a modification changing the second residue to a Histidine, believed to mimic the charge distribution of natural acetylation by NatB. Constructs tested in normal yeast strains remained viable and grew normally in the presence of Met-His Arp2p and tropomyosin. Analysis of their ability to suppress the mutant phenotypes of arp2-1 and arm1- mutants is an area of research to be explored in future studies.

Evaluation of predicted network modules in yeast metabolism using NMR-based metabolite profiling

Genome-scale metabolic models promise important insights into cell function. However, the definition of pathways and functional network modules within these models, and in the biochemical literature in general, is often based on intuitive reasoning. Although mathematical methods have been proposed to identify modules, which are defined as groups of reactions with correlated fluxes, there is a need for experimental verification. We show here that multivariate statistical analysis of the NMR-derived intra- and extracellular metabolite profiles of single-gene deletion mutants in specific metabolic pathways in the yeast Saccharomyces cerevisiae identified outliers whose profiles were markedly different from those of the other mutants in their respective pathways. Application of flux coupling analysis to a metabolic model of this yeast showed that the deleted gene in an outlying mutant encoded an enzyme that was not part of the same functional network module as the other enzymes in the pathway. We suggest that metabolomic methods such as this, which do not require any knowledge of how a gene deletion might perturb the metabolic network, provide an empirical method for validating and ultimately refining the predicted network structure.

In vivo and in vitro function of human UDP-galactose 4 '-epimerase variants

Type III galactosemia results from reduced activity of the enzyme UDP-galactose 4'-epimerase. Five disease-associated alleles (G90E, V94M, D103G, N34S and L183P) and three artificial alleles (Y105C, N268D, and M284K) were tested for their ability to alleviate galactose-induced growth arrest in a Saccharomyces cerevisiae strain which lacks endogenous UDP-galactose 4'-epimerase. For all of these alleles, except M284K, the ability to alleviate galactose sensitivity was correlated with the UDP-galactose 4'-epimerase activity detected in cell extracts. The M284K allele, however, was able to substantially alleviate galactose sensitivity, but demonstrated near-zero activity in cell extracts. Recombinant expression of the corresponding protein in Escherichia coil resulted in a protein with reduced enzymatic activity and reduced stability towards denaturants in vitro. This lack of stability may result from the introduction of an unpaired positive charge into a bundle of three alpha-helices near the surface of the protein. The disparities between the in vivo and in vitro data for M284K-hGALE further suggest that there are additional, stabilising factors present in the cell. Taken together, these results reinforce the need for care in the interpretation of in vitro, enzymatic diagnostic tests for type III galactosemia. (C) 2011 Elsevier Masson SAS. All rights reserved.

The GAL genetic switch: visualisation of the interacting proteins by split-EGFP bimolecular fluorescence complementation

A split-EGFP bimolecular fluorescence complementation assay was used to visualise and locate three interacting pairs of proteins from the GAL genetic switch of the budding yeast, Saccharomyces cerevisiae. Both the Gal4p-Gal80p and Gal80p-Gal3p pairs were found to be located in the nucleus under inducing conditions. However, the Gal80p-Gal1p complex was located throughout the cell. These results support recent work establishing an initial interaction between Gal3p and Gal80p occurring in the nucleus. Labelling of all three protein pairs impaired the growth of the yeast strains and resulted in reduced galactokinase activity in cell extracts. The most likely cause of this impairment is decreased dissociation rates of the complexes, caused by the essentially irreversible reassembly of the EGFP fragments. This suggests that a fully functional GAL genetic switch requires dynamic interactions between the protein components. These results also highlight the need for caution in the interpretation of in vivo split-EGFP experiments.

Artemisinins act through at least two targets in a yeast model

Artemisinin and related compounds are potent and widely used antimalarial drugs but their biochemical mode of action is not clear. There is strong evidence that ATP-dependent calcium transporters are a key target in the malarial parasite. However, work using Saccharomyces cerevisiae suggests that disruption of mitochondrial function is critical in the cell killing activity of these compounds. Here it is shown that, in the absence of reducing agents, artemisinin and artesunate targeted the S. cerevisiae calcium channels Pmr1p and Pmc1p. Both compounds affected the growth of yeast on fermentable and nonfermentable media. This growth inhibition was not seen in a yeast strain in which the genes encoding both calcium channels were deleted. In the presence of reducing agents, which break the endoperoxide bridge in the drugs, growth inhibition was only observed in nonfermentable media. This inhibition could be partially relieved by the addition of a free radical scavenger. These results suggest that the drugs have two biochemical modes of action - one acting by specific binding to calcium channels and one involving free radical production in the mitochondria.

Wzx proteins involved in biosynthesis of O antigen function in association with the first sugar of the O-specific lipopolysaccharide subunit

One of the most common pathways for the export of O-specific lipopolysaccharide (LPS) across the plasma membrane requires the participation of the Wzx protein. Wzx belongs to a family of integral membrane proteins that share little conservation in their primary amino acid sequence, making it difficult to delineate functional domains. This paper reports the cloning and expression in Escherichia coli K-12 of various Wzx homologues from different bacteria as FLAG epitope-tagged protein fusions. A reconstitution system for O16 LPS synthesis was used to assess the ability of each Wzx protein to complement an E. coli K-12 Deltawzx mutant. The results demonstrate that Wzx proteins from O-antigen systems that use N-acetylglucosamine or N-acetylgalactosamine for the initiation of the biosynthesis of the O repeat can fully complement the formation of O16 LPS. Partial complementation was seen with Wzx from Pseudomonas aeruginosa, a system that uses N-acetylfucosamine in the initiation reaction. In contrast, there was negligible complementation with the Wzx protein from Salmonella enterica, a system in which galactose is the initiating sugar. These results support a model whereby the first sugar of the O repeat can be recognized by the O-antigen translocation machinery.

Surface expression of O-specific lipopolysaccharide in Escherichia coli requires the function of the TolA protein

We investigated the involvement of Tol proteins in the surface expression of lipopolysaccharide (LPS). tolQ, -R, -A and -B mutants of Escherichia coli K-12, which do not form a complete LPS-containing O antigen, were transformed with the O7+ cosmid pJHCV32. The tolA and tolQ mutants showed reduced O7 LPS expression compared with the respective isogenic parent strains. No changes in O7 LPS expression were found in the other tol mutants. The O7-deficient phenotype in the tolQ and tolA mutants was complemented with a plasmid encoding the tolQRA operon, but not with a similar plasmid containing a frameshift mutation inactivating tolA. Therefore, the reduction in O7 LPS was attributed to the lack of a functional tolA gene, caused either by a direct mutation of this gene or by a polar effect on tolA gene expression exerted by the tolQ mutation. Reduced surface expression of O7 LPS was not caused by changes in lipid A-core structure or downregulation of the O7 LPS promoter. However, an abnormal accumulation of radiolabelled mannose was detected in the plasma membrane. As mannose is a sugar unique to the O7 subunit, this result suggested the presence of accumulated O7 LPS biosynthesis intermediates. Attempts to construct a tolA mutant in the E. coli O7 wild-type strain VW187 were unsuccessful, suggesting that this mutation is lethal. In contrast, a polar tolQ mutation affecting tolA expression in VW187 caused slow growth rate and serum sensitivity in addition to reduced O7 LPS production. VW187 tolQ cells showed an elongated morphology and became permeable to the membrane-impermeable dye propidium iodide. All these phenotypes were corrected upon complementation with cloned tol genes but were not restored by complementation with the tolQRA operon containing the frameshift mutation in tolA. Our results demonstrate that the TolA protein plays a critical role in the surface expression of O antigen subunits by an as yet uncharacterized involvement in the processing of O antigen.