Secretomes of medically important fungi reflect morphological and phylogenetic diversity

Secretome represents a main target for understanding the mechanisms of fungal adaptation. In the present study, we focus on the secretomes of fungi associated with infections in humans and other mammals in order to explore relationships between the diverse morphological and phylogenetic groups. Almost all the mammalian pathogenic fungi analyzed have secretome sizes smaller than 1000 proteins and, secreted proteins comprise between 5% and 10% of the total proteome. As expected, the correlation pattern between the secretome size and the total proteome was similar to that described in previous secretome studies of fungi. With regard to the morphological groups, minimum secretome sizes of less than 250 secreted proteins and low values for the fraction of secreted proteins are shown in mammalian pathogenic fungi with reduced proteomes such as microsporidia, atypical fungi and some species of yeasts and yeast-like fungi (Malassezia). On the other hand, filamentous fungi have significantly more secreted proteins and the highest numbers are present in species of filamentous fungi that also are plant or insect pathogens (Fusarium verticilloides, Fusarium oxysporum and Basidiobolus meristosporus). With respect to phylogeny, there are also variations in secretome size across fungal subphyla: Microsporidia, Taphrinomycotina, Ustilagomycotina and Saccharomycotina contain small secretomes; whereas larger secretomes are found in Agaricomycotina, Pezizomycotina, Mucoromycotina and Entomophthoromycotina. Finally, principal component analysis (PCA) was conducted on the complete secretomes. The PCA results revealed that, in general, secretomes of fungi belonging to the same morphological group or subphyla cluster together. In conclusion, our results point out that in medically important fungi there is a relationship between the secretome and the morphological group or phylogenetic classification.

VerSeDa: vertebrate secretome database

Based on the current tools, de novo secretome (full set of proteins secreted by an organism) prediction is a time consuming bioinformatic task that requires a multifactorial analysis in order to obtain reliable in silico predictions. Hence, to accelerate this process and offer researchers a reliable repository where secretome information can be obtained for vertebrates and model organisms, we have developed VerSeDa (Vertebrate Secretome Database). This freely available database stores information about proteins that are predicted to be secreted through the classical and non-classical mechanisms, for the wide range of vertebrate species deposited at the NCBI, UCSC and ENSEMBL sites. To our knowledge, VerSeDa is the only state-of-the-art database designed to store secretome data from multiple vertebrate genomes, thus, saving an important amount of time spent in the prediction of protein features that can be retrieved from this repository directly.

Database URL: VerSeDa is freely available at http://genomics.cicbiogune.es/VerSeDa/index.php

Expansion of Signal Transduction Pathways in Fungi by Extensive Genome Duplication

Plants and fungi use light and other signals to regulate development, growth, and metabolism. The fruiting bodies of the fungus Phycomyces blakesleeanus are single cells that react to environmental cues, including light, but the mechanisms are largely unknown [1]. The related fungus Mucor circinelloides is an opportunistic human pathogen that changes its mode of growth upon receipt of signals from the environment to facilitate pathogenesis [2]. Understanding how these organisms respond to environmental cues should provide insights into the mechanisms of sensory perception and signal transduction by a single eukaryotic cell, and their role in pathogenesis. We sequenced the genomes of P. blakesleeanus and M. circinelloides and show that they have been shaped by an extensive genome duplication or, most likely, a whole-genome duplication (WGD), which is rarely observed in fungi [3-6]. We show that the genome duplication has expanded gene families, including those involved in signal transduction, and that duplicated genes have specialized, as evidenced by differences in their regulation by light. The transcriptional response to light varies with the developmental stage and is still observed in a photoreceptor mutant of P. blakesleeanus. A phototropic mutant of P. blakesleeanus with a heterozygous mutation in the photoreceptor gene madA demonstrates that photosensor dosage is important for the magnitude of signal transduction. We conclude that the genome duplication provided the means to improve signal transduction for enhanced perception of environmental signals. Our results will help to understand the role of genome dynamics in the evolution of sensory perception in eukaryotes.

Genomewide analysis of phytochrome proteins in the phylum Basidiomycota

Phytochromes are photoreceptor proteins involved in the detection of the red and far-red regions of the visible light spectrum. Fungal phytochromes are hybrid histidine kinases with a conserved domain architecture composed of an N-terminal photosensory module and a C-terminal regulatory output module that includes the histidine kinase and response regulator receiver domains. In this study, we have analyzed the distribution, domain architecture, and phylogenetic analysis of phytochrome proteins in 47 published genome sequences among the phylum Basidiomycota. Genome analysis revealed that almost every genome of basidiomycetes contained at least one gene encoding a phytochrome protein. Domain architecture of fungal phytochromes was completely conserved in the identified phytochromes of basidiomycetes, and phylogenetic analysis clustered these proteins into clades related with the phylogenetic classification of this fungal phylum.

Highly expressed captured genes and cross-kingdom domains present in Helitrons create novel diversity in Pleurotus ostreatus and other fungi

BACKGROUND

Helitrons are class-II eukaryotic transposons that transpose via a rolling circle mechanism. Due to their ability to capture and mobilize gene fragments, they play an important role in the evolution of their host genomes. We have used a bioinformatics approach for the identification of helitrons in two Pleurotus ostreatus genomes using de novo detection and homology-based searching. We have analyzed the presence of helitron-captured genes as well as the expansion of helitron-specific helicases in fungi and performed a phylogenetic analysis of their conserved domains with other representative eukaryotic species.

RESULTS

Our results show the presence of two helitron families in P. ostreatus that disrupt gene colinearity and cause a lack of synteny between their genomes. Both putative autonomous and non-autonomous helitrons were transcriptionally active, and some of them carried highly expressed captured genes of unknown origin and function. In addition, both families contained eukaryotic, bacterial and viral domains within the helitron's boundaries. A phylogenetic reconstruction of RepHel helicases using the Helitron-like and PIF1-like helicase conserved domains revealed a polyphyletic origin for eukaryotic helitrons.

CONCLUSION

P. ostreatus helitrons display features similar to other eukaryotic helitrons and do not tend to capture host genes or gene fragments. The occurrence of genes probably captured from other hosts inside the helitrons boundaries pose the hypothesis that an ancient horizontal transfer mechanism could have taken place. The viral domains found in some of these genes and the polyphyletic origin of RepHel helicases in the eukaryotic kingdom suggests that virus could have played a role in a putative lateral transfer of helitrons within the eukaryotic kingdom. The high similarity of some helitrons, along with the transcriptional activity of its RepHel helicases indicates that these elements are still active in the genome of P. ostreatus.

Two-component signal transduction in Agaricus bisporus: a comparative genomic analysis with other basidiomycetes through the web-based tool BASID2CS

Two-component systems (TCSs) are signal transduction mechanisms present in many eukaryotes, including fungi that play essential roles in the regulation of several cellular functions and responses. In this study, we carry out a genomic analysis of the TCS proteins in two varieties of the white button mushroom Agaricus bisporus. The genomes of both A. bisporus varieties contain eight genes coding for TCS proteins, which include four hybrid Histidine Kinases (HKs), a single histidine-containing phosphotransfer (HPt) protein and three Response Regulators (RRs). Comparison of the TCS proteins among A. bisporus and the sequenced basidiomycetes showed a conserved core complement of five TCS proteins including the Tco1/Nik1 hybrid HK, HPt protein and Ssk1, Skn7 and Rim15-like RRs. In addition, Dual-HKs, unusual hybrid HKs with 2 HK and 2 RR domains, are absent in A. bisporus and are limited to various species of basidiomycetes. Differential expression analysis showed no significant up- or down-regulation of the Agaricus TCS genes in the conditions/tissue analyzed with the exception of the Skn7-like RR gene (Agabi_varbisH97_2|198669) that is significantly up-regulated on compost compared to cultured mycelia. Furthermore, the pipeline web server BASID2CS (http://bioinformatics.unavarra.es:1000/B2CS/BASID2CS.htm) has been specifically designed for the identification, classification and functional annotation of putative TCS proteins from any predicted proteome of basidiomycetes using a combination of several bioinformatic approaches.

The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi

Brown rot decay removes cellulose and hemicellulose from wood--residual lignin contributing up to 30% of forest soil carbon--and is derived from an ancestral white rot saprotrophy in which both lignin and cellulose are decomposed. Comparative and functional genomics of the "dry rot" fungus Serpula lacrymans, derived from forest ancestors, demonstrated that the evolution of both ectomycorrhizal biotrophy and brown rot saprotrophy were accompanied by reductions and losses in specific protein families, suggesting adaptation to an intercellular interaction with plant tissue. Transcriptome and proteome analysis also identified differences in wood decomposition in S. lacrymans relative to the brown rot Postia placenta. Furthermore, fungal nutritional mode diversification suggests that the boreal forest biome originated via genetic coevolution of above- and below-ground biota.

Comparative genomic analysis of two-component regulatory proteins in Pseudomonas syringae

Background

Pseudomonas syringae is a widespread bacterial plant pathogen, and strains of P. syringae may be assigned to different pathovars based on host specificity among different plant species. The genomes of P. syringae pv. syringae (Psy) B728a, pv. tomato (Pto) DC3000 and pv. phaseolicola (Pph) 1448A have been recently sequenced providing a major resource for comparative genomic analysis. A mechanism commonly found in bacteria for signal transduction is the two-component system (TCS), which typically consists of a sensor histidine kinase (HK) and a response regulator (RR). P. syringae requires a complex array of TCS proteins to cope with diverse plant hosts, host responses, and environmental conditions.

Results

Based on the genomic data, pattern searches with Hidden Markov Model (HMM) profiles have been used to identify putative HKs and RRs. The genomes of Psy B728a, Pto DC3000 and Pph 1448A were found to contain a large number of genes encoding TCS proteins, and a core of complete TCS proteins were shared between these genomes: 30 putative TCS clusters, 11 orphan HKs, 33 orphan RRs, and 16 hybrid HKs. A close analysis of the distribution of genes encoding TCS proteins revealed important differences in TCS proteins among the three P. syringae pathovars.

Conclusion

In this article we present a thorough analysis of the identification and distribution of TCS proteins among the sequenced genomes of P. syringae. We have identified differences in TCS proteins among the three P. syringae pathovars that may contribute to their diverse host ranges and association with plant hosts. The identification and analysis of the repertoire of TCS proteins in the genomes of P. syringae pathovars constitute a basis for future functional genomic studies of the signal transduction pathways in this important bacterial phytopathogen.

Molecular and insecticidal characterization of a Cry1I protein toxic to insects of the families Noctuidae, Tortricidae, Plutellidae, and Chrysomelidae

The most notable characteristic of Bacillus thuringiensis is its ability to produce insecticidal proteins. More than 300 different proteins have been described with specific activity against insect species. We report the molecular and insecticidal characterization of a novel cry gene encoding a protein of the Cry1I group with toxic activity towards insects of the families Noctuidae, Tortricidae, Plutellidae, and Chrysomelidae. PCR analysis detected a DNA sequence with an open reading frame of 2.2 kb which encodes a protein with a molecular mass of 80.9 kDa. Trypsin digestion of this protein resulted in a fragment of ca. 60 kDa, typical of activated Cry1 proteins. The deduced sequence of the protein has homologies of 96.1% with Cry1Ia1, 92.8% with Cry1Ib1, and 89.6% with Cry1Ic1. According to the Cry protein classification criteria, this protein was named Cry1Ia7. The expression of the gene in Escherichia coli resulted in a protein that was water soluble and toxic to several insect species. The 50% lethal concentrations for larvae of Earias insulana, Lobesia botrana, Plutella xylostella, and Leptinotarsa decemlineata were 21.1, 8.6, 12.3, and 10.0 microg/ml, respectively. Binding assays with biotinylated toxins to E. insulana and L. botrana midgut membrane vesicles revealed that Cry1Ia7 does not share binding sites with Cry1Ab or Cry1Ac proteins, which are commonly present in B. thuringiensis-treated crops and commercial B. thuringiensis-based bioinsecticides. We discuss the potential of Cry1Ia7 as an active ingredient which can be used in combination with Cry1Ab or Cry1Ac in pest control and the management of resistance to B. thuringiensis toxins.

Extracytoplasmic function sigma factors in Pseudomonas syringae

Genome analyses of the plant pathogens Pseudomonas syringae pv. tomato DC3000, pv. syringae B728a and pv. phaseolicola 1448A reveal fewer extracytoplasmic function (ECF) sigma factors than in related Pseudomonads with different lifestyles. We highlight the presence of a P. syringae-specific ECF sigma factor that is an interesting target for future studies because of its potential role in the adaptation of P. syringae to its specialized phytopathogenic lifestyle.

The VirPphA/AvrPtoB family of type III effectors in Pseudomonas syringae

The VirPphA/AvrPtoB family of type III effector proteins from the phytopathogenic bacterium Pseudomonas syringae is one of the models providing insights into the molecular mechanisms conferring plant disease resistance and pathogenesis. In this review we summarize recent advances concerning the VirPphA/AvrPtoB family of effectors involved in the elicitation and suppression of plant defense responses.

Pseudomonas syringae pv. phaseolicola can be separated into two genetic lineages distinguished by the possession of the phaseolotoxin biosynthetic cluster

The bean (Phaseolus spp.) plant pathogen Pseudomonas syringae pv. phaseolicola is characterized by the ability to produce phaseolotoxin (Tox(+)). We recently reported that the majority of the Spanish P. syringae pv. phaseolicola population is unable to synthesize this toxin (Tox(-)). These Tox(-) isolates appear to lack the entire DNA region for the biosynthesis of phaseolotoxin (argK-tox gene cluster), as shown by PCR amplification and DNA hybridization using DNA sequences specific for separated genes of this cluster. Tox(+) and Tox(-) isolates also showed genomic divergence that included differences in ERIC-PCR and arbitrarily primed-PCR profiles. Tox(+) isolates showed distinct patterns of IS801 genomic insertions and contained a chromosomal IS801 insertion that was absent from Tox(-) isolates. Using a heteroduplex mobility assay, sequence differences were observed only among the intergenic transcribed spacer of the five rDNA operons of the Tox(-) isolates. The techniques used allowed the unequivocal differentiation of isolates of P. syringae pv. phaseolicola from the closely related soybean (Glycine max) pathogen, P. syringae pv. glycinea. Finally, a pathogenicity island that is essential for the pathogenicity of P. syringae pv. phaseolicola on beans appears to be conserved among Tox(+), but not among Tox(-) isolates, which also lacked the characteristic large plasmid that carries this pathogenicity island. It is proposed that the results presented here justify the separation of the Tox(+) and Tox(-) P. syringae pv. phaseolicola isolates into two distinct genetic lineages, designated Pph1 and Pph2, respectively, that show relevant genomic differences that include the pathogenicity gene complement.

Use of a PCR method on fecal samples for diagnosis of sheep paratuberculosis

The high sensitivity of PCR compared to the difficulties of fecal culture in sheep prompted the development of PCR protocols for detection of Mycobacterium avium subsp. paratuberculosis DNA in sheep feces. Although the PCR itself is well developed, and does not pose large technical problems, concentrating the bacteria from samples that may contain low numbers of bacilli using practical methods is still the main difficulty for the use of this technique. In this study, we describe an extraction protocol for the concentration and purification of M. avium subsp. paratuberculosis DNA from fecal samples and we compare it with other methods. The diagnostic performance of the freeze-boiling method was evaluated using a reference method [Vet. Rec. 134 (1994) 95] on fecal samples from a group of selected sheep from different flocks of known individual serological, pathological, and cultural paratuberculosis status. Using, as a reference, a combination of results in those conventional methods, the freeze-boiling PCR protocol showed a sensitivity of 94.1%, and a specificity of 92.3%.

Involvement of the sigmaN DNA-binding domain in open complex formation

sigmaN (sigma54) RNA polymerase holoenzyme closed complexes isomerize to open complexes in a reaction requiring nucleoside triphosphate hydrolysis by enhancer binding activator proteins. Here, we characterize Klebsiella pneumoniae sigmaN mutants, altered in the carboxy DNA-binding domain (F354A/F355A, F402A, F403A and F402A/F403A), that fail in activator-dependent transcription. The mutant holoenzymes have altered activator-dependent interactions with promoter sequences that normally become melted. Activator-dependent stable complexes accumulated slowly in vitro (F402A) and to a reduced final level (F403A, F402A/F403A, F354A/F355A). Similar results were obtained in an assay of activator-independent stable complex formation. Premelted templates did not rescue the mutants for stable preinitiation complex formation but did for deleted region I sigmaN, suggesting different defects. The DNA-binding domain substitutions are within sigmaN sequences previously shown to be buried upon formation of the wild-type holoenzyme or closed complex, suggesting that, in the mutants, alteration of the sigmaN-core and sigmaN-DNA interfaces has occurred to change holoenzyme activity. Core-binding assays with the mutant sigmas support this view. Interestingly, an internal deletion form of sigmaN lacking the major core binding determinant was able to assemble into holoenzyme and, although unable to support activator-dependent transcription, formed a stable activator-independent holoenzyme promoter complex on premelted DNA templates.

DNA-binding domain mutants of sigma-N (sigmaN, sigma54) defective between closed and stable open promoter complex formation

The sigmaN RNA polymerase binds promoters in a transcriptionally inactive form. Activation by enhancer binding positive control proteins results in the formation of an open promoter complex. In the closed complex, DNA sequences melted upon activation are close contacted by the sigmaN C-terminal DNA-binding domain. Conserved phenylalanine residues within the DNA-binding domain were mutated to examine their contribution to sigmaN function. Mutants defective in supporting sigmaN-dependent growth and in vivo promoter activation were obtained. The mutant proteins were able to bind promoter DNA and to form an RNA polymerase holoenzyme closed complex in vitro. However, they were defective in response to activator in vitro. They failed in the formation of heparin-stable promoter complexes characteristic of open promoter complexes. The sigmaN mutant forms, displaying good promoter occupancy but poor open complex formation, appear defective for some function of the holoenzyme required after initial promoter recognition. The possibilities that the defect could be located in a DNA contact important for DNA melting or is associated with activator interaction and conformational change in sigmaN are discussed.

Transcriptional analysis of the sigA and sigB genes of Brevibacterium lactofermentum

Transcription of the sigA and sigB genes of Brevibacterium lactofermentum, encoding the principal sigma factors of RNA polymerase, has been investigated by Northern blot and primer extension analysis. Northern hybridizations revealed that sigA is transcribed as a monocistronic mRNA of 1.7 kb and sigB gives two transcripts of 1.1 and 1.5 kb. Similar transcription patterns of sigA and sigB genes in nutrient-rich medium were observed; transcripts of both genes occurred simultaneously throughout the exponential growth phase and decayed clearly when the stationary phase was reached. Primer extension analysis of B. lactofermentum RNA showed that the sigA and sigB transcription initiation sites are located 17 bp and 24 bp upstream from the first nucleotide of the respective translation initiation codons. Alignment of the sigA and sigB promoter regions provided evidence for highly conserved sequences in both of them.

The galE gene encoding the UDP-galactose 4-epimerase of Brevibacterium lactofermentum is coupled transcriptionally to the dmdR gene

The galE gene of Brevibacterium lactofermentum, encoding UDP-galactose 4-epimerase (EC 5.1.3.2), has been identified by DNA sequencing downstream from the orf1-sigB-dmdR region. The arrangement of the sigB-dtxR-galE cluster is also conserved in Corynebacterium diphtheriae. The deduced galE product was a protein of 329 aa residues (35.4 kDa) that shared a high degree of identity to known UDP-galactose 4-epimerase proteins from Gram-positive microorganisms (Streptomyces lividans and Streptococcus thermophilus). Transcriptional analysis of the dmdR and galE genes in nutrient-rich medium showed that these genes are part of an operon, that is actively transcribed as a bicistronic mRNA during the exponential growth phase, but transcription of the operon is decreased during the stationary growth phase. In addition, the dmdR gene was also expressed as a monocistronic 0.7-kb transcript during the active growth phase.

Multiple sigma factor genes in Brevibacterium lactofermentum: characterization of sigA and sigB

Four rpoD hybridizing signals have been identified in the chromosome of Brevibacterium lactofermentum. Two rpoD-like genes, sigA and sigB, have been cloned and sequenced, and they encode principal sigma factors of the RNA polymerase. The deduced amino acid sequences of SigA and SigB showed very high similarities to those of Mycobacterium smegmatis MysA and MysB proteins, respectively, and also to those of HrdB proteins from different Streptomyces species. SigA and SigB maintain the conserved motifs of sigma 70-like principal sigma factors. sigB is closely linked to the dtxR gene (encoding a repressor of iron-regulated promoters homologous to the diphtheria toxin repressor from Corynebacterium diphtheriae.

Three genes hrdB, hrdD and hrdT of Streptomyces griseus IMRU 3570, encoding sigma factor-like proteins, are differentially expressed under specific nutritional conditions

Three genes (hrd) homologous to the rpoD gene of Escherichia coli, that encode sigma factor-like proteins, have been cloned from DNA of the candicidin-producing strain Streptomyces griseus IMRU 3570. They are located in different regions of the chromosome. Sequence analysis showed that the first one is analogous to the hrdB gene of S. coelicolor. The second showed high similarity to the hrdD gene of S. coelicolor and S. aureofaciens and is linked, as in S. coelicolor, to a N-acetyltransferase-encoding gene (nat) distantly related to the pat and bar genes that encode resistance to bialafos. The third showed no close homology with other known hrd genes from actinomycetes and has been named hrdT. Functional domains in the three S. griseus Hrd proteins are highly conserved in relation to those of the sigma 70 protein family. Northern analysis showed that hrdB is expressed as a 1.9-kb transcript during active growth in phosphate-rich medium, but it is less efficiently transcribed under sporulation conditions (phosphate-starved) or after a heat-shock treatment. Two other shorter transcripts of 1.2 and 0.7 kb were also detected with the same probe. The hrdD gene is transcribed as a single 1.1-kb transcript under sporulation conditions following nutritional shiftdown and, to a lower extent, during growth conditions in phosphate-rich medium. The hrdT gene is weakly transcribed (1.5-kb RNA) under all conditions tested. The hrd-encoded sigma factors probably recognize actinomycetes promoters (SEP type) with E. coli-like consensus sequences.

Molecular cloning, DNA sequence analysis, and characterization of the Corynebacterium diphtheriae dtxR homolog from Brevibacterium lactofermentum

A homolog of the Corynebacterium diphtheriae dtxR gene was isolated from Brevibacterium lactofermentum. The product of the B. lactofermentum dtxR gene was immunoreactive with polyclonal anti-DtxR antibodies and functioned as an iron-activated repressor capable of regulating the expression of beta-galactosidase from a diphtheria tox promoter/operator transcriptional fusion in recombinant Escherichia coli. The extents of induction by increasing concentrations of the chelator 2,2'-dipyridyl were identical in cells expressing DtxR from either C. diphtheriae or B. lactofermentum.

A gene encoding arginyl-tRNA synthetase is located in the upstream region of the lysA gene in Brevibacterium lactofermentum: regulation of argS-lysA cluster expression by arginine

The Brevibacterium lactofermentum argS gene, which encodes an arginyl-tRNA synthetase, was identified in the upstream region of the lysA gene. The cloned gene was sequenced; it encodes a 550-amino-acid protein with an M(r) of 59,797. The deduced amino acid sequence showed 28% identical and 49% similar residues when compared with the sequence of the Escherichia coli arginyl-tRNA synthetase. The B. lactofermentum enzyme showed the highly conserved motifs of class I aminoacyl-tRNA synthetases. Expression of the argS gene in B. lactofermentum and E. coli resulted in an increase in aminoacyl-tRNA synthetase activity, correlated with the presence in sodium dodecyl sulfate-polyacrylamide gels of a clear protein band that corresponds to this enzyme. One single transcript of about 3,000 nucleotides and corresponding to the B. lactofermentum argS-lysA operon was identified. The transcription of these genes is repressed by lysine and induced by arginine, showing an interesting pattern of biosynthetic interlock between the pathways of both amino acids in corynebacteria.

A cluster of three genes (dapA, orf2, and dapB) of Brevibacterium lactofermentum encodes dihydrodipicolinate synthase, dihydrodipicolinate reductase, and a third polypeptide of unknown function

The dapA and dapB genes, encoding, respectively, dihydrodipicolinate synthase and dihydrodipicolinate reductase, the two first enzymes of the lysine branch of the aspartic amino acid family, were cloned from the DNA of the amino acid-producing bacterium Brevibacterium lactofermentum. The two genes were clustered in a 3.5-kb Sau3AI-BamHI fragment but were separated by an open reading frame of 750 nucleotides. The protein encoded by this open reading frame had little similarity to any protein in the data banks, and its function remains unknown. The three genes were translated in Escherichia coli, giving the corresponding polypeptides.