Peter J. Schaap

Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88

The filamentous fungus Aspergillus niger is widely exploited by the fermentation industry for the production of enzymes and organic acids, particularly citric acid. We sequenced the 33.9-megabase genome of A. niger CBS 513.88, the ancestor of currently used enzyme production strains. A high level of synteny was observed with other aspergilli sequenced. Strong function predictions were made for 6,506 of the 14,165 open reading frames identified. A detailed description of the components of the protein secretion pathway was made and striking differences in the hydrolytic enzyme spectra of aspergilli were observed. A reconstructed metabolic network comprising 1,069 unique reactions illustrates the versatile metabolism of A. niger. Noteworthy is the large number of major facilitator superfamily transporters and fungal zinc binuclear cluster transcription factors, and the presence of putative gene clusters for fumonisin and ochratoxin A synthesis.

Comparative genomics of citric-acid-producing Aspergillus niger ATCC 1015 versus enzyme-producing CBS 513.88

The filamentous fungus Aspergillus niger exhibits great diversity in its phenotype. It is found globally, both as marine and terrestrial strains, produces both organic acids and hydrolytic enzymes in high amounts, and some isolates exhibit pathogenicity. Although the genome of an industrial enzyme-producing A. niger strain (CBS 513.88) has already been sequenced, the versatility and diversity of this species compel additional exploration. We therefore undertook whole-genome sequencing of the acidogenic A. niger wild-type strain (ATCC 1015) and produced a genome sequence of very high quality. Only 15 gaps are present in the sequence, and half the telomeric regions have been elucidated. Moreover, sequence information from ATCC 1015 was used to improve the genome sequence of CBS 513.88. Chromosome-level comparisons uncovered several genome rearrangements, deletions, a clear case of strain-specific horizontal gene transfer, and identification of 0.8 Mb of novel sequence. Single nucleotide polymorphisms per kilobase (SNPs/kb) between the two strains were found to be exceptionally high (average: 7.8, maximum: 160 SNPs/kb). High variation within the species was confirmed with exo-metabolite profiling and phylogenetics. Detailed lists of alleles were generated, and genotypic differences were observed to accumulate in metabolic pathways essential to acid production and protein synthesis. A transcriptome analysis supported up-regulation of genes associated with biosynthesis of amino acids that are abundant in glucoamylase A, tRNA-synthases, and protein transporters in the protein producing CBS 513.88 strain. Our results and data sets from this integrative systems biology analysis resulted in a snapshot of fungal evolution and will support further optimization of cell factories based on filamentous fungi.

RNA Targeting by the Type III-A CRISPR-Cas Csm Complex of Thermus thermophilus

CRISPR-Cas is a prokaryotic adaptive immune system that provides sequence-specific defense against foreign nucleic acids. Here we report the structure and function of the effector complex of the Type III-A CRISPR-Cas system of Thermus thermophilus: the Csm complex (TtCsm). TtCsm is composed of five different protein subunits (Csm1-Csm5) with an uneven stoichiometry and a single crRNA of variable size (35-53 nt). The TtCsm crRNA content is similar to the Type III-B Cmr complex, indicating that crRNAs are shared among different subtypes. A negative stain EM structure of the TtCsm complex exhibits the characteristic architecture of Type I and Type III CRISPR-associated ribonucleoprotein complexes. crRNA-protein crosslinking studies show extensive contacts between the Csm3 backbone and the bound crRNA. We show that, like TtCmr, TtCsm cleaves complementary target RNAs at multiple sites. Unlike Type I complexes, interference by TtCsm does not proceed via initial base pairing by a seed sequence.

3DM: Systematic analysis of heterogeneous superfamily data to discover protein functionalities

Ten years of experience with molecular class–specific information systems (MCSIS) such as with the hand‐curated G protein–coupled receptor database (GPCRDB) or the semiautomatically generated nuclear receptor database has made clear that a wide variety of questions can be answered when protein‐related data from many different origins can be flexibly combined. MCSISes revolve around a multiple sequence alignment (MSA) that includes “all” available sequences from the entire superfamily, and it has been shown at many occasions that the quality of these alignments is the most crucial aspect of the MCSIS approach. We describe here a system called 3DM that can automatically build an entire MCSIS. 3DM bases the MSA on a multiple structure alignment, which implies that the availability of a large number of superfamily members with a known three‐dimensional structure is a requirement for 3DM to succeed well. Thirteen MCSISes were constructed and placed on the Internet for examination. These systems have been instrumental in a large series of research projects related to enzyme activity or the understanding and engineering of specificity, protein stability engineering, DNA‐diagnostics, drug design, and so forth. Proteins 2010.

An iterative workflow for mining the human intestinal metaproteome.

BackgroundPeptide spectrum matching (PSM) is the standard method in shotgun proteomics data analysis. It relies on the availability of an accurate and complete sample proteome that is used to make interpretation of the spectra feasible. Although this procedure has proven to be effective in many proteomics studies, the approach has limitations when applied on complex samples of microbial communities, such as those found in the human intestinal tract. Metagenome studies have indicated that the human intestinal microbiome contains over 100 times more genes than the human genome and it has been estimated that this ecosystem contains over 5000 bacterial species. The genomes of the vast majority of these species have not yet been sequenced and hence their proteomes remain unknown. To enable data analysis of shotgun proteomics data using PSM, and circumvent the lack of a defined matched metaproteome, an iterative workflow was developed that is based on a synthetic metaproteome and the developing metagenomic databases that are both representative for but not necessarily originating from the sample of interest.ResultsTwo human fecal samples for which metagenomic data had been collected, were analyzed for their metaproteome using liquid chromatography-mass spectrometry and used to benchmark the developed iterative workflow to other methods. The results show that the developed method is able to detect over 3,000 peptides per fecal sample from the spectral data by circumventing the lack of a defined proteome without naive translation of matched metagenomes and cross-species peptide identification.ConclusionsThe developed iterative workflow achieved an approximate two-fold increase in the amount of identified spectra at a false discovery rate of 1% and can be applied in metaproteomic studies of the human intestinal tract or other complex ecosystems.

Oxaloacetate Hydrolase, the C–C Bond Lyase of Oxalate Secreting Fungi

Oxalate secretion by fungi is known to be associated with fungal pathogenesis. In addition, oxalate toxicity is a concern for the commercial application of fungi in the food and drug industries. Although oxalate is generated through several different biochemical pathways, oxaloacetate acetylhydrolase (OAH)-catalyzed hydrolytic cleavage of oxaloacetate appears to be an especially important route. Below, we report the cloning of the Botrytis cinerea oahA gene and the demonstration that the disruption of this gene results in the loss of oxalate formation. In addition, through complementation we have shown that the intact B. cinerea oahA gene restores oxalate production in an Aspergillus niger mutant strain, lacking a functional oahA gene. These observations clearly indicate that oxalate production in A. niger and B. cinerea is solely dependent on the hydrolytic cleavage of oxaloacetate catalyzed by OAH. In addition, the B. cinera oahA gene was overexpressed in Escherichia coli and the purified OAH was used to define catalytic efficiency, substrate specificity, and metal ion activation. These results are reported along with the discovery of the mechanism-based, tight binding OAH inhibitor 3,3-difluorooxaloacetate (Ki = 68 nm). Finally, we propose that cellular uptake of this inhibitor could reduce oxalate production.

A comprehensive metatranscriptome analysis pipeline and its validation using human small intestine microbiota datasets

BackgroundNext generation sequencing (NGS) technologies can be applied in complex microbial ecosystems for metatranscriptome analysis by employing direct cDNA sequencing, which is known as RNA sequencing (RNA-seq). RNA-seq generates large datasets of great complexity, the comprehensive interpretation of which requires a reliable bioinformatic pipeline. In this study, we focus on the development of such a metatranscriptome pipeline, which we validate using Illumina RNA-seq datasets derived from the small intestine microbiota of two individuals with an ileostomy.ResultsThe metatranscriptome pipeline developed here enabled effective removal of rRNA derived sequences, followed by confident assignment of the predicted function and taxonomic origin of the mRNA reads. Phylogenetic analysis of the small intestine metatranscriptome datasets revealed a strong similarity with the community composition profiles obtained from 16S rDNA and rRNA pyrosequencing, indicating considerable congruency between community composition (rDNA), and the taxonomic distribution of overall (rRNA) and specific (mRNA) activity among its microbial members. Reproducibility of the metatranscriptome sequencing approach was established by independent duplicate experiments. In addition, comparison of metatranscriptome analysis employing single- or paired-end sequencing methods indicated that the latter approach does not provide improved functional or phylogenetic insights. Metatranscriptome functional-mapping allowed the analysis of global, and genus specific activity of the microbiota, and illustrated the potential of these approaches to unravel syntrophic interactions in microbial ecosystems.ConclusionsA reliable pipeline for metatransciptome data analysis was developed and evaluated using RNA-seq datasets obtained for the human small intestine microbiota. The set-up of the pipeline is very generic and can be applied for (bacterial) metatranscriptome analysis in any chosen niche.

Identification of a Gatekeeper Residue That Prevents Dehydrogenases from Acting as Oxidases

The oxygen reactivity of flavoproteins is poorly understood. Here we show that a single Ala to Gly substitution in l-galactono-γ-lactone dehydrogenase (GALDH) turns the enzyme into a catalytically competent oxidase. GALDH is an aldonolactone oxidoreductase with a vanillyl-alcohol oxidase (VAO) fold. We found that nearly all oxidases in the VAO family contain either a Gly or a Pro at a structurally conserved position near the C4a locus of the isoalloxazine moiety of the flavin, whereas dehydrogenases prefer another residue at this position. Mutation of the corresponding residue in GALDH (Ala-113 → Gly) resulted in a striking 400-fold increase in oxygen reactivity, whereas the cytochrome c reductase activity is retained. The activity of the A113G variant shows a linear dependence on oxygen concentration (kox = 3.5 × 105 m-1 s-1), similar to most other flavoprotein oxidases. The Ala-113 → Gly replacement does not change the reduction potential of the flavin but creates space for molecular oxygen to react with the reduced flavin. In the wild-type enzyme, Ala-113 acts as a gatekeeper, preventing oxygen from accessing the isoalloxazine nucleus. The presence of such an oxygen access gate seems to be a key factor for the prevention of oxidase activity within the VAO family and is absent in members that act as oxidases.

An inventory of the Aspergillus niger secretome by combining in silico predictions with shotgun proteomics data

BackgroundThe ecological niche occupied by a fungal species, its pathogenicity and its usefulness as a microbial cell factory to a large degree depends on its secretome. Protein secretion usually requires the presence of a N-terminal signal peptide (SP) and by scanning for this feature using available highly accurate SP-prediction tools, the fraction of potentially secreted proteins can be directly predicted. However, prediction of a SP does not guarantee that the protein is actually secreted and current in silico prediction methods suffer from gene-model errors introduced during genome annotation.ResultsA majority rule based classifier that also evaluates signal peptide predictions from the best homologs of three neighbouring Aspergillus species was developed to create an improved list of potential signal peptide containing proteins encoded by the Aspergillus niger genome. As a complement to these in silico predictions, the secretome associated with growth and upon carbon source depletion was determined using a shotgun proteomics approach. Overall, some 200 proteins with a predicted signal peptide were identified to be secreted proteins. Concordant changes in the secretome state were observed as a response to changes in growth/culture conditions. Additionally, two proteins secreted via a non-classical route operating in A. niger were identified.ConclusionsWe were able to improve the in silico inventory of A. niger secretory proteins by combining different gene-model predictions from neighbouring Aspergilli and thereby avoiding prediction conflicts associated with inaccurate gene-models. The expected accuracy of signal peptide prediction for proteins that lack homologous sequences in the proteomes of related species is 85%. An experimental validation of the predicted proteome confirmed in silico predictions.

Correlated mutation analyses on super‐family alignments reveal functionally important residues

Correlated mutation analyses (CMA) on multiple sequence alignments are widely used for the prediction of the function of amino acids. The accuracy of CMA‐based predictions is mainly determined by the number of sequences, by their evolutionary distances, and by the quality of the alignments. These criteria are best met in structure‐based sequence alignments of large super‐families. So far, CMA‐techniques have mainly been employed to study the receptor interactions. The present work shows how a novel CMA tool, called Comulator, can be used to determine networks of functionally related residues in enzymes. These analyses provide leads for protein engineering studies that are directed towards modification of enzyme specificity or activity. As proof of concept, Comulator has been applied to four enzyme super‐families: the isocitrate lyase/phoshoenol‐pyruvate mutase super‐family, the hexokinase super‐family, the RmlC‐like cupin super‐family, and the FAD‐linked oxidases super‐family. In each of those cases networks of functionally related residue positions were discovered that upon mutation influenced enzyme specificity and/or activity as predicted. We conclude that CMA is a powerful tool for redesigning enzyme activity and selectivity. Proteins 2009.