Molecular evolution of the HSP70 multigene family

Eukaryotic genomes encode multiple 70-kDa heat-shock proteins (HSP70s). The Saccharomyces cerevisiae HSP70 family is comprised of eight members. Here we present the nucleotide sequence of the SSA3 and SSB2 genes, completing the nucleotide sequence data for the yeast HSP70 family. We have analyzed these yeast sequences as well as 29 HSP70s from 24 additional eukaryotic and prokaryotic species. Comparison of the sequences demonstrates the extreme conservation of HSP70s; proteins from the most distantly related species share at least 45% identity and more than one-sixth of the amino acids are identical in the aligned region (567 amino acids) among all proteins analyzed. Phylogenetic trees constructed by two independent methods indicate that ancient molecular and cellular events have given rise to at least four monophyletic groups of eukaryotic HSP70 proteins. Each group of evolutionarily similar HSP70s shares a common intracellular localization and is presumed to be comprised of functional homologues; these include heat-shock proteins of the cytoplasm, endoplasmic reticulum, mitochondria, and chloroplasts. HSP70s localized in mitochondria and plastids are most similar to the DnaK HSP70 homologues in purple bacteria and cyanobacteria, respectively, which is consistent with the proposed prokaryotic origin of these organelles. The analyses indicate that the major eukaryotic HSP70 groups arose prior to the divergence of the earliest eukaryotes, roughly 2 billion years ago. In some cases, as exemplified by the SSA genes encoding the cytoplasmic HSP70s of S. cerevisiae, more recent duplication events have given rise to subfamilies within the major groups. The S. cerevisiae SSB proteins comprise a unique subfamily not identified in other species to date. This subfamily appears to have resulted from an ancient gene duplication that occurred at approximately the same time as the origin of the major eukaryotic HSP70 groups. Correspondence to: E.A. Craig     

Molecular population genetics and evolution of a prion-like protein in Saccharomyces cerevisiae

The prion-like behavior of Sup35p, the eRF3 homolog in the yeast Saccharomyces cerevisiae, mediates the activity of the cytoplasmic nonsense suppressor known as [PSI(+)]. Sup35p is divided into three regions of distinct function. The N-terminal and middle (M) regions are required for the induction and propagation of [PSI(+)] but are not necessary for translation termination or cell viability. The C-terminal region encompasses the termination function. The existence of the N-terminal region in SUP35 homologs of other fungi has led some to suggest that this region has an adaptive function separate from translation termination. To examine this hypothesis, we sequenced portions of SUP35 in 21 strains of S. cerevisiae, including 13 clinical isolates. We analyzed nucleotide polymorphism within this species and compared it to sequence divergence from a sister species, S. paradoxus. The N domain of Sup35p is highly conserved in amino acid sequence and is highly biased in codon usage toward preferred codons. Amino acid changes are under weak purifying selection based on a quantitative analysis of polymorphism and divergence. We also conclude that the clinical strains of S. cerevisiae are not recently derived and that outcrossing between strains in S. cerevisiae may be relatively rare in nature.     

Transfer RNA methytransferases and their corresponding modifications in budding yeast and humans: activities, predications, and potential roles in human health

Throughout the kingdoms of life, transfer RNA (tRNA) undergoes over 100 enzyme-catalyzed, methyl-based modifications. Although a majority of the methylations are conserved from bacteria to mammals, the functions of a number of these modifications are unknown. Many of the proteins responsible for tRNA methylation, named tRNA methyltransferases (Trms), have been characterized in Saccharomyces cerevisiae. In contrast, only a few human Trms have been characterized. A BLAST search for human homologs of each S. cerevisiae Trm revealed a total of 34 human proteins matching our search criteria for an S. cerevisiae Trm homolog candidate. We have compiled a database cataloging basic information about each human and yeast Trm. Every S. cerevisiae Trm has at least one human homolog, while several Trms have multiple candidates. A search of cancer cell versus normal cell mRNA expression studies submitted to Oncomine found that 30 of the homolog genes display a significant change in mRNA expression levels in at least one data set. While 6 of the 34 human homolog candidates have confirmed tRNA methylation activity, the other candidates remain uncharacterized. We believe that our database will serve as a resource for investigating the role of human Trms in cellular stress signaling.     

Analysis of the pdx-1 (snz-1/sno-1) region of the Neurospora crassa genome: correlation of pyridoxine-requiring phenotypes with mutations in two structural genes

We report the analysis of a 36-kbp region of the Neurospora crassa genome, which contains homologs of two closely linked stationary phase genes, SNZ1 and SNO1, from Saccharomyces cerevisiae. Homologs of SNZ1 encode extremely highly conserved proteins that have been implicated in pyridoxine (vitamin B6) metabolism in the filamentous fungi Cercospora nicotianae and in Aspergillus nidulans. In N. crassa, SNZ and SNO homologs map to the region occupied by pdx-1 (pyridoxine requiring), a gene that has been known for several decades, but which was not sequenced previously. In this study, pyridoxine-requiring mutants of N. crassa were found to possess mutations that disrupt conserved regions in either the SNZ or SNO homolog. Previously, nearly all of these mutants were classified as pdx-1. However, one mutant with a disrupted SNO homolog was at one time designated pdx-2. It now appears appropriate to reserve the pdx-1 designation for the N. crassa SNZ homolog and pdx-2 for the SNO homolog. We further report annotation of the entire 36,030-bp region, which contains at least 12 protein coding genes, supporting a previous conclusion of high gene densities (12,000-13,000 total genes) for N. crassa. Among genes in this region other than SNZ and SNO homologs, there was no evidence of shared function. Four of the genes in this region appear to have been lost from the S. cerevisiae lineage.     

Characterization of α-factor pheromone and pheromone receptor genes of Ashbya gossypii

The genome of Ashbya gossypii contains homologs of most of the genes that are part of the Saccharomyces cerevisiae pheromone-signal transduction cascade. However, we currently lack understanding of a potential sexual cycle for this pre-whole genome duplication hemiascomycete. The sequenced strain bears three identical copies encoding MATa. We show that the syntenic A. gossypii homolog of MFα1 (AFL062w) does not encode a mature α-factor peptide, but identified another gene, AAR163c, which encodes a candidate α-specific mating pheromone and is thus reannotated as AgMFα2. The expression of the AgSTE2α-factor receptor in an Scste2 S. cerevisiae MATa strain resulted in dosage-dependent growth arrest upon exposure to A. gossypiiα-factor, which indicated that the pheromone response was effectively coupled to the S. cerevisiae signal transduction cascade. Comparison of α-pheromones and α-pheromone receptors showed greater conservation between Eremothecium cymbalariae and S. cerevisiae than between A. gossypii and E. cymbalariae. We constructed A. gossypii strains deleted for the STE2 and STE3 pheromone receptors. These strains showed no phenotypic abnormalities and an ste2, ste3 double mutant is still able to sporulate. The deletion of STE12 as the downstream target of pheromone signalling, however, led to a hypersporulation phenotype.     

Mapping by sequencing the Pneumocystis genome using the ordering DNA sequences V3 tool

A bioinformatics tool called ODS3 has been created for mapping by sequencing. The tool allows the creation of integrated genomic maps from genetic, physical mapping, and sequencing data and permits an integrated genome map to be stored, retrieved, viewed, and queried in a stand-alone capacity, in a client/server relationship with the Fungal Genome Database (FGDB), and as a web-browsing tool for the FGDB. In that ODS3 is programmed in Java, the tool promotes platform independence and supports export of integrated genome-mapping data in the extensible markup language (XML) for data interchange with other genome information systems. The tool ODS3 is used to create an initial integrated genome map of the AIDS-related fungal pathogen, Pneumocystis carinii. Contig dynamics would indicate that this physical map is approximately 50% complete with approximately 200 contigs. A total of 10 putative multigene families were found. Two of these putative families were previously characterized in P. carinii, namely the major surface glycoproteins (MSGs) and HSP70 proteins; three of these putative families (not previously characterized in P. carinii) were found to be similar to families encoding the HSP60 in Schizosaccharomyces pombe, the heat-shock psi protein in S. pombe, and the RNA synthetase family (i.e., MES1) in Saccharomyces cerevisiae. Physical mapping data are consistent with the 16S, 5.8S, and 26S rDNA genes being single copy in P. carinii. No other fungus outside this genus is known to have the rDNA genes in single copy.     

Nonlinear relationships among evolutionary rates identify regions of functional divergence in heat-shock protein 70 genes

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP- dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily C (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.      

Pneumocystis carinii sterol 14α-demethylase activity in Saccharomyces cerevisiae erg11 knockout mutant: sterol biochemistry

Pneumocystis carinii is an unusual fungus that can cause pneumonitis in immunosuppressed laboratory rats. Reactions in sterol biosynthesis are attractive targets for development of antimycotic drugs. A key enzyme in sterol biosynthesis is sterol 14α-demethylase (14DM), which is coded by the erg11 gene. Here we describe detailed sterol analysis of wild-type Saccharomyces cerevisiae and in an erg11 knockout mutant expressing either P. carinii or S. cerevisiae 14DM from a plasmid-borne cDNA. Sterols of the three strains were qualitatively and quantitatively analyzed using thin-layer chromatography, high-performance liquid chromatography, and gas-liquid chromatography and mass spectrometry and nuclear magnetic resonance spectroscopy. Biochemical evidence for functional complementation was provided by detecting the same major sterols in all three strains with ergosterol being by far the most abundant. A total of 25 sterols was identified, 16 of which were identified in all three strains. The ratios of lanosterol:14-desmethyllanosterol in the three strains indicate that the mutant transformed with erg11 showed more 14DM activity than wild-type yeast. The sterol analyses also indicated that the P. carinii 14DM can utilize the sterol substrates used by the S. cerevisiae 14DM and suggested that the yeast 14DM in the yeast cell utilizes 4α-methyl sterols better than the P. carinii enzyme.     

Evolution of HSP70 gene and its implications regarding relationships between archaebacteria,eubacteria, and eukaryotes

The 70-kDa heat-shock protein (HSP70) constitutes the most conserved protein present in all organisms that is known to date. Based on global alignment of HSP70 sequences from organisms representing all three domains, numerous sequence signatures that are specific for prokaryotic and eukaryotic homologs have been identified. HSP70s from the two archaebacterial species examined (viz., Halobacterium marismortui and Methanosarcina mazei) have been found to contain all eubacterial but no eukaryotic signature sequences. Based on several novel features of the HSP70 family of proteins (viz., presence of tandem repeats of a 9-amino-acid [a.a.] polypeptide sequence and structural similarity between the first and second quadrants of HSP70, homology of the N-terminal half of HSP70 to the bacterial MreB protein, presence of a conserved insert of 23–27 a.a. in all HSP70s except those from archaebacteria and gram-positive eubacteria) a model for the evolution of HSP70 gene from an early stage is proposed. The HSP70 homologs from archaebacteria and gram-positive bacteria lacking the insert in the N-terminal quadrants are indicated to be the ancestral form of the protein. Detailed phylogenetic analyses of HSP70 sequence data (viz., by bootstrap analyses, maximum parsimony, and maximum likelihood methods) provide evidence that archaebacteria are not monophyletic and show a close evolutionary linkage with the gram-positive eubacteria. These results do not support the traditional archaebacterial tree, where a close relationship between archaebacterial and eukaryotic homologs is observed. To explain the phylogenies based on HSP70 and other gene sequences, a model for the origin of eukaryotic cells involving fusion between archaebacteria and gram-negative eubacteria is proposed. Correspondence to: R. S. Gupta     

Regulation of RNA processing and transport by a nuclear guanine nucleotide release protein and members of the Ras superfamily.

The RCC1 gene of mammals encodes a guanine nucleotide release protein (GNRP). RCC1 and a homolog in Saccharomyces cerevisiae (MTR1/PRP20/SRM1) have previously been implicated in control of mRNA metabolism and export from the nucleus. We here demonstrate that a temperature-sensitive fission yeast mutant which has a mutation in a homologous gene, and two of three additional (mtr1/prp20/srm1) mutants accumulate nuclear poly(A)+ RNA at 37 degrees C. In S.cerevisiae, maturation of rRNA and tRNA is also inhibited at 37 degrees C. Nevertheless, studies with the corresponding BHK-21 cell mutant indicate that protein import into the nucleus continues. MTR1 homologs regulate RNA processing at a point which is distinct from their regulation of chromosome condensation since: (i) poly(A)+ RNA accumulation in the fission yeast mutant precedes chromosome condensation, and (ii) unlike chromosome condensation, accumulation of nuclear poly(A)+ RNA does not require p34cdc28 kinase activation or protein synthesis. Moreover, experiments involving inhibition of DNA synthesis indicate that the S.cerevisiae homolog does not govern cell cycle checkpoint control. Since RCC1p acts as GNRP for Ran, a small nuclear GTPase of the ras superfamily, we have identified two homologs of Ran in S.cerevisiae (CNR1 and CNR2). Only CNR1 is essential, but both code for proteins extremely similar to Ran and can suppress mtr1 mutations in allele-specific fashion. Thus, MTR1 and its homologs appear to act as GNRPs for a family of conserved GTPases in controlling RNA metabolism and transport. Their role in governing checkpoint control appears to be restricted to higher eukaryotes.     

BiP, HSP70, NDK and PDI in wheat endosperm. II. Effects of high temperature on protein and mRNA accumulation

The effects of high temperature on accumulation of the 70‐kDa heat shock protein (HSP70) and nucleoside diphosphate kinase (NDK) as well as two other proteins that have roles in the biosynthesis of storage proteins were examined during grain development. An HSP70 homolog and a 17‐kDa NDK were co‐purified from wheat endosperm, their identity verified, and a cDNA for an HSP70 expressed in endosperm was isolated. Wheat plants ( Triticum aestivum , cvs Butte and Vulcan) were heat shocked at 40°C or exposed to maximum daily temperatures of 37 or 40°C during early or mid‐grain fill. Antibodies and cDNA probes for BiP, HSP70, NDK and PDI were used to examine the effect of high temperatures on the accumulation of protein and mRNA in the endosperm. HSP70 mRNA levels increased substantially when plants were exposed to heat shock or to a 1‐day gradual increase to 40°C. The effects of a 5‐day heat treatment on mRNA levels were more complicated and depended on the developmental stage of the grain. A treatment that began at 7 days post‐anthesis (DPA) decreased the level of mRNA for HSP70, BiP, PDI and NDK, whereas a treatment that began at 14 DPA slightly increased mRNA levels. The same treatments increased the accumulation of HSP70 but did not affect BiP, PDI, or NDK protein levels. This is the first detailed report on the effects of heat on mRNA and protein levels for HSP70 in a developing seed storage tissue.     

Identification of an Arabidopsis thaliana cDNA encoding a HSP70-related protein belonging to the HSP110/SSE1 subfamily

Heat-shock protein 70 (HSP70)-related proteins are classified in two main subfamilies: the DnaK subfamily and the HSP110/SSE1 subfamily. We have characterized the first plant member of the HSP110/SSE1 subfamily, HSP91. At least two, tightly linked genes encoding HSP91 are present per haploid Arabidopsis genome. HSP91 is constitutively expressed in non-stressed Arabidopsis plants and is transiently induced by heat shock.     

Cloning and Identification of Arp1, an Actin-Related Protein from Pneumocystis carinii

ABSTRACT. The complete Pneumocystis carinii Arp1 gene has been sequenced from two cDNA clones. The gene encodes a protein 385 bp in length with an estimated size of 45,000 kD. The A + T% for the Arp1 gene and a 900-bp sequence upstream of the gene were 63.7% and 70.3%, respectively. These values are consistant with A + T codon preference displayed by P. carinii and are similar to values reported for other P. carinii genes. The predicted amino acid sequence of the P. carinii Arp1 protein had a similarity of 87.6% with Neurospora crassa Arp1, 82.1% similarity with vertebrate centractin, and 71.2% similarity with the Saccharomyces cerevisiae Act5p. Expression of Arp1 mRNA in P. carinii was detectable via synthesis of cDNA and subsequent PCR amplification. Affinity purified antibodies against S. cerevisiae Act5p, and canine centractin recognized both the recombinantly expressed protein and a 45,000 kD protein in P. carinii nuclear extracts. The Arp1 gene is the second member of the actin multigene family that has been identified in P. carinii .