Transcriptome analysis in response to heat shock and cadmium in the aquatic fungus Blastocladiella emersonii

The global transcriptional response of the chytridiomycete Blastocladiella emersonii to environmental stress conditions was explored by sequencing a large number of expressed sequence tags (ESTs) from three distinct cDNA libraries, constructed with mRNA extracted from cells exposed to heat shock and different concentrations of cadmium chloride. A total of 6,350 high-quality EST sequences were obtained and assembled into 2,326 putative unigenes, 51% of them not previously described in B. emersonii. To approximately 59% of the unigenes it was possible to assign an orthologue in another organism, whereas 41% of them remained without a putative identification, with transcripts related to protein folding and antioxidant activity being highly enriched in the stress libraries. A microarray chip was constructed encompassing 3,773 distinct ESTs from the B. emersonii transcriptome presently available, which correspond to a wide range of biological processes. Global gene expression analysis of B. emersonii cells exposed to stress conditions revealed a large number of differentially expressed genes: 122 up- and 60 downregulated genes during heat shock and 189 up- and 110 downregulated genes during exposure to cadmium. The main functional categories represented among the upregulated genes were protein folding and proteolysis, proteins with antioxidant properties, and cellular transport. Interestingly, in response to cadmium stress, B. emersonii cells induced genes encoding six different glutathione S-transferases and six distinct metacaspases, as well as genes coding for several proteins of sulfur amino acid metabolism, indicating that cadmium causes oxidative stress and apoptosis in this fungus. All sequences described in this study have been submitted to the GenBank EST section with the accession numbers EE 730389 to EE 736848.

Comparative expression analysis of members of the Hsp70 family in the chytridiomycete Blastocladiella emersonii

Sequencing of a large number of expressed sequence tags from Blastocladiella emersonii revealed the presence of ten distinct putative members of the 70 kDa-heat shock protein (Hsp70) family in this fungus. The amino acid sequence deduced from eight of these cDNAs showed significant similarity to members of the Saccharomyces cerevisiae Hsp70 family, and the remaining displayed high sequence homology with hsp70 gene products from other organisms. The hsp70-3 gene was the most highly expressed at normal temperatures and was poorly induced during heat shock. Except for hsp70-4 and hsp70-6, all other hsp70 genes were induced to different degrees upon exposure of B. emersonii cells to heat shock, with hsp70-1 gene presenting the highest transcript levels. Phylogenetic analysis of complete B. emersonii putative Hsp70 protein sequences indicated that Hsp70-1 and Hsp70-3 corresponded to cytosolic proteins, whereas Hsp70-7 and Hsp70-9 are probably localized in the endoplasmic reticulum and mitochondria, respectively.

Comparative EST analysis provides insights into the basal aquatic fungus Blastocladiella emersonii

Background

Blastocladiella emersonii is an aquatic fungus of the Chytridiomycete class, which is at the base of the fungal phylogenetic tree. In this sense, some ancestral characteristics of fungi and animals or fungi and plants could have been retained in this aquatic fungus and lost in members of late-diverging fungal species. To identify in B. emersonii sequences associated with these ancestral characteristics two approaches were followed: (1) a large-scale comparative analysis between putative unigene sequences (uniseqs) from B. emersonii and three databases constructed ad hoc with fungal proteins, animal proteins and plant unigenes deposited in Genbank, and (2) a pairwise comparison between B. emersonii full-length cDNA sequences and their putative orthologues in the ascomycete Neurospora crassa and the basidiomycete Ustilago maydis.

Results

Comparative analyses of B. emersonii uniseqs with fungi, animal and plant databases through the two approaches mentioned above produced 166 B. emersonii sequences, which were identified as putatively absent from other fungi or not previously described. Through these approaches we found: (1) possible orthologues of genes previously identified as specific to animals and/or plants, and (2) genes conserved in fungi, but with a large difference in divergence rate in B. emersonii. Among these sequences, we observed cDNAs encoding enzymes from coenzyme B₁₂-dependent propionyl-CoA pathway, a metabolic route not previously described in fungi, and validated their expression in Northern blots.

Conclusion

Using two different approaches involving comparative sequence analyses, we could identify sequences from the early-diverging fungus B. emersonii previously considered specific to animals or plants, and highly divergent sequences from the same fungus relative to other fungi.

SpotWhatR: a user-friendly microarray data analysis system

SpotWhatR is a user-friendly microarray data analysis tool that runs under a widely and freely available R statistical language (http://www.r-project.org) for Windows and Linux operational systems. The aim of SpotWhatR is to help the researcher to analyze microarray data by providing basic tools for data visualization, normalization, determination of differentially expressed genes, summarization by Gene Ontology terms, and clustering analysis. SpotWhatR allows researchers who are not familiar with computational programming to choose the most suitable analysis for their microarray dataset. Along with well-known procedures used in microarray data analysis, we have introduced a stand-alone implementation of the HTself method, especially designed to find differentially expressed genes in low-replication contexts. This approach is more compatible with our local reality than the usual statistical methods. We provide several examples derived from the Blastocladiella emersonii and Xylella fastidiosa Microarray Projects. SpotWhatR is freely available at http://blasto.iq.usp.br/~tkoide/SpotWhatR, in English and Portuguese versions. In addition, the user can choose between "single experiment" and "batch processing" versions.

Gene discovery and expression profile analysis through sequencing of expressed sequence tags from different developmental stages of the chytridiomycete Blastocladiella emersonii

Blastocladiella emersonii is an aquatic fungus of the chytridiomycete class which diverged early from the fungal lineage and is notable for the morphogenetic processes which occur during its life cycle. Its particular taxonomic position makes this fungus an interesting system to be considered when investigating phylogenetic relationships and studying the biology of lower fungi. To contribute to the understanding of the complexity of the B. emersonii genome, we present here a survey of expressed sequence tags (ESTs) from various stages of the fungal development. Nearly 20,000 cDNA clones from 10 different libraries were partially sequenced from their 5' end, yielding 16,984 high-quality ESTs. These ESTs were assembled into 4,873 putative transcripts, of which 48% presented no matches with existing sequences in public databases. As a result of Gene Ontology (GO) project annotation, 1,680 ESTs (35%) were classified into biological processes of the GO structure, with transcription and RNA processing, protein biosynthesis, and transport as prevalent processes. We also report full-length sequences, useful for construction of molecular phylogenies, and several ESTs that showed high similarity with known proteins, some of which were not previously described in fungi. Furthermore, we analyzed the expression profile (digital Northern analysis) of each transcript throughout the life cycle of the fungus using Bayesian statistics. The in silico approach was validated by Northern blot analysis with good agreement between the two methodologies.

DNA nucleotide sequence homologies between some zoosporic fungi

Representative species from the monoflagellate Blastocladiales and the biflagellate Saprolegniales were studied for their DNA base composition, heterogeneity, nucleotide sequence homology and divergence. Intergeneric, intrageneric and interstrain DNAs of Blastocladiales were heterogenous. The G+C values for their main component (average 64 percent) and two minor ones (average 52 and 44 percent) were found to be significantly higher than the corresponding values from the biflagellate Saprolegnia ferax (55, 46 and 36 percent respectively). In Allomyces species, the two hybrid, male and female strains were found to have closer homology with their parental types than these last between themselves. Among Blastocladiales, interspecific similarities between the epigynous A. macrogynus and the hypogynous A. arbuscula were higher (average 75 percent) than intergeneric similarities between Allomyces and Blastocladiella (average 58 percent). The biflagellate S. ferax was found to be distantly related to the uniflagellate Allomyces (average 48 percent similarity). The nucleotide sequence divergences obtained from thermal elution data correlated the hybridization values.

Blastocladiella emersoniiexpresses a centrin similar toChlamydomonas reinhardtiiisoform not found in late-diverging fungi

Centrins are members of the calcium-binding EF-hand protein superfamily which can be divided into two subfamilies, probably associated with different functions: one related to Chlamydomonas reinhardtii centrin, CrCenp, and the other, represented by Saccharomyces cerevisiae isoform, ScCdc31p. ESTs encoding the two isoforms (BeCen1 and BeCen3) from the chytridiomycete Blastocladiella emersonii were isolated, and expression of the CrCenp-type centrin, BeCen1, was analyzed throughout the fungus life cycle. Becen1 mRNA levels increase transiently during sporulation and protein levels present a similar pattern. Immunolocalization studies seem to localize BeCen1 at the basal body zone and in the cytoplasm surrounding the nuclear cap, a zoospore organelle.

Characterization and expression of two genes encoding isoforms of a putative Na, K-ATPase in the chytridiomycete Blastocladiella emersonii

A P-type ATPase gene (BePAT1) from the aquatic fungus Blastocladiella emersonii, which surprisingly showed high similarity with the alpha-subunit of Na, K-ATPases from animal cells, has been reported recently [Biochim. Biophys. Acta 1383 (1998) 183]. In the present study, we describe the characterization of a second gene, denominated BePAT2, and show that these two genes have a different intron-exon structure but encode putative proteins with greater than 90% amino acid identity. Northern blot and multiplex reverse transcription and polymerase chain reaction (RT-PCR) assays have revealed that BePAT1 and BePAT2 genes have a non-coordinate, developmentally regulated expression during B. emersonii life cycle. Phosphoenzyme formation experiments using the immunopurified enzymes have indicated the presence of a Na, K-ATPase-like activity. Furthermore, immunofluorescence studies using B. emersonii zoospores localized the ATPases on the plasma membrane of these cells.

Structure, expression, and functional analysis of the gene coding for calmodulin in the chytridiomycete Blastocladiella emersonii

The single calmodulin (CaM) gene and the corresponding cDNA of the chytridiomycete Blastocladiella emersonii were isolated and characterized. The CaM gene is interrupted by three introns and transcribed in a single 0.7-kb mRNA species encoding a predicted protein 91% identical to human CaM. B. emersonii CaM has been expressed in Escherichia coli as a fusion protein with gluthatione S-transferase (GST) and purified by affinity chromatography and cleavage from the GST portion using a site-specific protease. In the presence of Ca(2+), B. emersonii CaM exhibited a shift in apparent molecular mass similar to that observed with bovine CaM and was able to activate the autophosphorylation of CaM-dependent protein kinase II (CaMKII) from rat brain. CaM expression is developmentally regulated in B. emersonii, with CaM mRNA and protein concentrations increasing during sporulation to maximum levels observed just prior to the release of the zoospores into the medium. Both CaM protein and mRNA levels decrease drastically at the zoospore stage, increasing again during germination. The CaM antagonists compound 48/80, calmidazolium, and W7 were shown to completely inhibit B. emersonii sporulation when added to the cultures at least 120, 150, and 180 min after induction, respectively. All these drugs also inhibited growth and zoospore production in this fungus. The Ca(2+) channel blocker TMB-8 and the CaMKII inhibitor KN93 completely inhibited sporulation if added up to 60 min after induction of this stage, but only KN93 affected fungal growth. The data presented suggest that the Ca(2+)-CaM complex and CaMKII play an important role during growth and sporulation in B. emersonii.

Characterization and submitochondrial localization of the alpha subunit of the mitochondrial processing peptidase from the aquatic fungus Blastocladiella emersonii

In an effort to investigate the molecular mechanisms responsible for the drastic morphological changes the mitochondria go through during the life cycle of the aquatic fungus Blastocladiella emersonii, the gene encoding the alpha subunit of the mitochondrial processing peptidase (alpha-MPP) was isolated. Nucleotide sequence analysis revealed that the predicted alpha-MPP polypeptide comprises 474 amino acids with a calculated molecular mass of 51,900 Da, presenting a characteristic mitochondrial signal sequence. Northern blot analysis indicated a single 1.4-kb transcript encoding the B. emersonii alpha-MPP, whose levels decrease during sporulation, becoming very low in the zoospore, and increase again during germination. Despite these variations in mRNA concentration, B. emersonii alpha-MPP protein levels do not change significantly during the life cycle of the fungus, as observed in Western blots. Experiments to investigate the submitochondrial localization of B. emersonii alpha-MPP and beta-MPP were also carried out, and the results indicated that both subunits are associated with the mitochondrial inner membrane, possibly as part of the bc1 complex, as described for plants.

Isolation, characterization, and expression of the gene encoding the beta subunit of the mitochondrial processing peptidase from Blastocladiella emersonii

A 2.3-kb BamHI-KpnI fragment was isolated from a partial genomic library and shown by nucleotide sequence analysis to contain the entire coding region of the gene encoding the beta subunit of the Blastocladiella mitochondrial processing peptidase (beta-MPP). The predicted beta-MPP protein has 465 amino acids and a calculated molecular mass of 50.8 kDa. S1 nuclease protection assays revealed an intron, 209 bp in size, interrupting the coding region between the putative signal sequence and the mature protein. Northern blot analysis showed that beta-MPP mRNA levels decrease significantly during B. emersonii sporulation, reaching basal levels in the zoospore stage. The amount of beta-MPP protein, determined in Western blots, unlike its mRNA, does not vary significantly throughout the fungal life cycle.

A P-type ATPase from the aquatic fungus Blastocladiella emersonii similar to animal Na,K-ATPases

We have cloned a P-type ATPase gene from the aquatic fungus Blastocladiella emersonii (BePAT1) using a probe obtained with degenerate oligonucleotides, corresponding to two amino acid sequences highly conserved among all P-type ATPase isoforms, and the polymerase chain reaction technique. Nucleotide sequence analysis revealed a 3.4 kb open reading frame encoding a putative peptide of 1080 amino acid residues with a calculated molecular mass of 119 kDa, which presents all diagnostic features of P-type transporting ATPases. Comparison to other members of the family and phylogenetic analyses have shown that the BePAT1 protein belongs to the subfamily of Na,K- and H,K-ATPases, indicating that the divergence between the alpha-subunit of the Na,K-ATPase and other members of the P-type ATPase family has occurred before the divergence between the animal and fungal lineages in evolution.

Protein factors in Blastocladiella emersonii cell extracts recognize similar sequence elements in the promoters of the genes encoding cAMP-dependent protein kinase subunits

Blastocladiella emersonii contains a single cAMP-dependent protein kinase (PKA), which is similar to the mammalian type II isoforms. Its activity is regulated during development by changes in the levels of the catalytic (C) and regulatory (R) subunits, which occur in parallel with changes in levels of the corresponding mRNAs, suggesting coordinate transcriptional control of the expression of both subunits. Both R and C mRNA levels are low in vegetative cells, rise sharply during sporulation and decrease to basal levels again after germination. To investigate sequence elements common to both Blastocladiella R and C gene promoters, which might be involved in the coordinate regulation of these genes, their 5'-flanking regions were analyzed by gel mobility shift and DNase I footprinting assays. We determined that different DNA-protein complexes are generated when fragments of the R and C gene promoters are incubated with extracts from cells expressing (sporulating cells) or not expressing (vegetative cells) both subunits, and competition experiments suggested that similar protein factors bind to both promoters. DNase I footprinting experiments have indicated that a sequence common to both R and C promoters, and similar to mammalian E-boxes, binds factors present in extracts from vegetative and sporulating cells, whereas sequences flanking the E-boxes in both promoters showed a change in the pattern of DNase I digestion only when the vegetative cell extract was used. This result suggests that the composition of the protein complexes binding to these regions changes during sporulation.

Large-subunit rRNA sequence of the chytridiomycete Blastocladiella emersonii, and implications for the evolution of zoosporic fungi

The 5.8S and 28S ribosomal RNA sequences of the chytridiomycete Blastocladiella emersonii were determined. These data were combined with 18S rRNA sequences in order to carry out a phylogenetic analysis based on distance matrix, parsimony, and maximum likelihood methods. The new data confirmed that chytridiomycetes are true fungi and not protists, as was already suggested on the basis of biochemical, ultrastructural, and 18S rRNA data. Within the fungal clade, B. emersonii formed the first line of divergence. The position of the fungi within the eukaryotic "crown" taxa was also reassessed, and the alveolate-stramenopile cluster appeared as their sister group. The stramenopiles also comprise a number of zoosporic fungi, which resemble chytridiomycetes in so many respects, e.g., production of motile spores, thallus morphology, and absorptive nutrition, that they have been classified together with them in the past. This suggests that the possible common ancestor of the fungi, stramenopiles, and alveolates may have been a zoosporic fungus, which would mean that zoosporic fungi are paraphyletic instead of polyphyletic as previously suggested.

A unique intron-containing hsp70 gene induced by heat shock and during sporulation in the aquatic fungus Blastocladiella emersonii

We have isolated and characterized cDNA and genomic DNA clones encoding the 70-kDa heat-shock protein (Hsp70) from the aquatic fungus Blastocladiella emersonii (Be). Nucleotide (nt) sequence analysis predicts an acidic protein containing 650 amino acids, with a calculated molecular mass of 70.8 kDa. The Be hsp70 gene is induced by heat shock (HS), as well as during sporulation of the fungus, and its coding region is interrupted by a single intron. All the evidence seems to indicate that this is the only hsp70 in Be. S1 nuclease protection assays revealed that splicing of the hsp70 intron is highly thermoresistant; at the lethal temperature of 42 degrees C, only 30% of the hsp70 mRNAs have not been processed. A single transcription start point (tsp), localized about 30 nt downstream from a putative TATA box, was determined both during HS and at normal temperatures. The promoter region presented several NGAAN repeats (where N is any nucleotide) characteristic of HS elements, as well as putative binding sites for ATF, Sp1 and two metal-responsive elements.

Cloning and characterization of the gene for the catalytic subunit of cAMP-dependent protein kinase in the aquatic fungus Blastocladiella emersonii

We have isolated and characterized cDNA and genomic DNA clones encoding the catalytic subunit (C) of cAMP-dependent protein kinase in the aquatic fungus Blastocladiella emersonii. The C-subunit amino acid sequence derived from the nucleotide sequence predicts a basic polypeptide of 424 residues, excluding the initiator methionine, which by amino-terminal sequence analysis has been shown to be absent from the mature protein. The Blastocladiella C presents a 70-amino-acid extension at the amino terminus, when aligned to the mouse C alpha subunit, being one of the largest C subunits already characterized. The B. emersonii C-gene-coding region is interrupted by three introns, ranging in size over 57-69 bp. The positions of the introns are quite different from those found in other species, suggesting a considerable amount of evolutionary drift in the gene structure. The 5'-flanking region lacks recognizable TATA or CCAAT sequences, is remarkably high in GC content (70%), and primer extension experiments indicate that transcription initiates from multiple sites. Several sequence motifs were identified in the promoter region which could be involved in the developmental control of this gene.

Cloning and structural analysis of the gene for the regulatory subunit of cAMP-dependent protein kinase in Blastocladiella emersonii

We have isolated and characterized cDNA and genomic DNA clones encoding the regulatory subunit of cAMP-dependent protein kinase in the aquatic fungus Blastocladiella emersonii. Nucleotide sequence analysis has shown that the predicted protein comprises 403 amino acids with a calculated molecular mass of 44,263 Da and an overall 40% identity to mammalian RII subunits, including a serine in the phosphorylation site, which confirms the Blastocladiella protein as a type II regulatory subunit. The B. emersonii R gene presents two introns, one located in the 5'-noncoding region, whereas the other interrupts the coding region, just after the dimerization domain of the protein. The promoter region does not contain recognizable TATA or CCAAT sequences and is very GC rich, a characteristic shared by mammalian cAMP-dependent protein kinase subunit genes previously analyzed. S1 mapping and primer extension experiments revealed multiple transcription initiation sites. Several sequence motifs were identified in the 5'-flanking region which could be responsible for the regulation of this gene.

Coordinate pretranslational control of cAMP-dependent protein kinase subunit expression during development in the water mold Blastocladiella emersonii

The aquatic fungus Blastocladiella emersonii provides a system for studying the regulation of expression of regulatory (R) and catalytic (C) subunits of cAMP-dependent protein kinase (PKA). Blastocladiella cells contain a single PKA with properties very similar to type II kinases of mammalian tissues. During development cAMP-dependent protein kinase activity and its associated cAMP-binding activity change drastically. We have previously shown that the increase in cAMP-binding activity during sporulation is due to de novo synthesis of R subunit and to an increase in the translatable mRNA coding for R (Marques et al., Eur. J. Biochem. 178, 803, 1989). In the present work we have continued these studies to investigate the mechanism by which the changes in the level of kinase activity take place. The C subunit of Blastocladiella has been purified; antiserum has been raised against it and used to determine amounts of C subunit throughout the fungus' life cycle. A sharp increase in C subunit content occurs during sporulation and peaks at the zoospore stage. Northern blot analyses, using Blastocladiella C and R cDNA probes, have shown that the levels of C and R mRNAs parallel their intracellular protein concentrations. These results indicate a coordinate pretranslational control for C and R subunit expression during differentiation in Blastocladiella.

Developmental regulation of expression of the regulatory subunit of the cAMP-dependent protein kinase of Blastocladiella emersonii

A monospecific polyclonal antiserum to the regulatory subunit (R) of the cAMP-dependent protein kinase of Blastocladiella emersonii has been developed by immunization with purified regulatory subunit. In Western blots, the antiserum displays high affinity and specificity for the intact R monomer of Mr = 58,000, as well as for its proteolytic products of Mr = 43,000 and Mr = 36,000, even though the antiserum has been raised against the Mr = 43,000 fragment. Western blots of cell extracts prepared at different times during the life cycle of the fungus indicate that the increase in cAMP-binding activity occurring during sporulation, as well as its decrease during germination, are associated with the accumulation of the regulatory subunit during sporulation and its disappearance during germination, respectively. Pulse labeling with [35S]methionine and immunoprecipitation indicate that the accumulation of R is due to its increased synthesis during sporulation. Two-dimensional gel electrophoresis of affinity purified cell extracts obtained after [35S]methionine pulse labeling during sporulation confirms de novo synthesis of R during this stage and furthermore shows that the protein is rapidly phosphorylated after its synthesis. In vitro translation studies using RNA isolated from different stages of the life cycle followed by immunoprecipitation have shown that the time course of expression of the mRNA coding for the regulatory subunit parallels the rate of its synthesis in vivo.

Differential expression of heat-shock proteins and spontaneous synthesis of HSP70 during the life cycle of Blastocladiella emersonii

The heat-shock response in Blastocladiella emersonii is dependent on the developmental stage. Cells exposed to elevated temperatures at different stages of the life cycle (sporulation, germination or growth) show a differential synthesis of heat-shock proteins (hsps). Of a total of 22 polypeptides induced, particular subsets of hsps appear in each phase, demonstrating a non-coordinate heat-shock gene expression. In contrast, heat-shock-related proteins (hsp76, hsp70, hsp39a) are spontaneously expressed at a high level during sporulation. By the criteria of two-dimensional gel electrophoresis and partial proteolysis mapping, the 70,000-Da protein, whose synthesis is induced spontaneously during sporulation, is indistinguishable from the heat-inducible hsp70. The techniques of in vitro translation, and Northern analysis using a Drosophila hsp70 probe, demonstrated that enhanced synthesis of hsp70, which occurs during heat-shock treatment and spontaneously during sporulation, is associated with an accumulation of hsp70 mRNA. These observations suggest that hsp70 gene expression is induced during sporulation.

Developmental changes in translatable RNA species and protein synthesis during sporulation in the aquatic fungus Blastocladiella emersonii

Protein synthesis during sporulation in Blastocladiella emersonii is developmentally regulated as revealed using [35S]methionine pulse labeling and two-dimensional gel electrophoresis. A large increase in the synthesis of several proteins is associated with particular stages. A large number of basic proteins are synthesized exclusively during late sporulation. Changes in translatable mRNA species were also detected by two-dimensional gel electrophoresis of the polypeptides produced in a cell-free rabbit reticulocyte lysate primed with RNA prepared at different stages of sporulation. The synthesis of several proteins during sporulation seems to be transcriptionally controlled. Most of the sporulation-specific messages are not present in the mature zoospores.

Nucleotide sequences of 5S ribosomal RNA from four oomycete and chytrid water molds

The nucleotide sequences of the 5S rRNAs of the oomycete water molds Saprolegnia ferax and Pythium hydnosporum and of the chytrid water molds Blastocladiella simplex and Phlyctochytrium irregulare were determined by chemical and enzymatic partial degradation of 3' and 5' end-labelled molecules, followed by gel sequence analysis. The two oomycete sequences differed in 24 positions and the two chytrid sequences differed in 27 positions. These pairs differed in a mean of 44 positions. The chytrid sequences clearly most resemble the sequence from the zygomycete Phycomyces, while the oomycete sequences appear to be allied with those from protozoa and slime molds.