Molecular characterization of ubiquitin genes from Aspergillus nidulans: mRNA expression on different stress and growth conditions

Module-detection approaches for the integration of multilevel omics data highlight the comprehensive response of Aspergillus fumigatus to caspofungin

Background

Omics data provide deep insights into overall biological processes of organisms. However, integration of data from different molecular levels such as transcriptomics and proteomics, still remains challenging. Analyzing lists of differentially abundant molecules from diverse molecular levels often results in a small overlap mainly due to different regulatory mechanisms, temporal scales, and/or inherent properties of measurement methods. Module-detecting algorithms identifying sets of closely related proteins from protein-protein interaction networks (PPINs) are promising approaches for a better data integration.

Results

Here, we made use of transcriptome, proteome and secretome data from the human pathogenic fungus Aspergillus fumigatus challenged with the antifungal drug caspofungin. Caspofungin targets the fungal cell wall which leads to a compensatory stress response. We analyzed the omics data using two different approaches: First, we applied a simple, classical approach by comparing lists of differentially expressed genes (DEGs), differentially synthesized proteins (DSyPs) and differentially secreted proteins (DSePs); second, we used a recently published module-detecting approach, ModuleDiscoverer, to identify regulatory modules from PPINs in conjunction with the experimental data. Our results demonstrate that regulatory modules show a notably higher overlap between the different molecular levels and time points than the classical approach. The additional structural information provided by regulatory modules allows for topological analyses. As a result, we detected a significant association of omics data with distinct biological processes such as regulation of kinase activity, transport mechanisms or amino acid metabolism. We also found a previously unreported increased production of the secondary metabolite fumagillin by A. fumigatus upon exposure to caspofungin. Furthermore, a topology-based analysis of potential key factors contributing to drug-caused side effects identified the highly conserved protein polyubiquitin as a central regulator. Interestingly, polyubiquitin UbiD neither belonged to the groups of DEGs, DSyPs nor DSePs but most likely strongly influenced their levels.

Conclusion

Module-detecting approaches support the effective integration of multilevel omics data and provide a deep insight into complex biological relationships connecting these levels. They facilitate the identification of potential key players in the organism’s stress response which cannot be detected by commonly used approaches comparing lists of differentially abundant molecules.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0620-8) contains supplementary material, which is available to authorized users.

Proteome-wide profiling of protein lysine acetylation in Aspergillus flavus

Protein lysine acetylation is a prevalent post-translational modification that plays pivotal roles in various biological processes in both prokaryotes and eukaryotes. Aspergillus flavus, as an aflatoxin-producing fungus, has attracted tremendous attention due to its health impact on agricultural commodities. Here, we performed the first lysine-acetylome mapping in this filamentous fungus using immune-affinity-based purification integrated with high-resolution mass spectrometry. Overall, we identified 1383 lysine-acetylation sites in 652 acetylated proteins, which account for 5.18% of the total proteins in A. flavus. According to bioinformatics analysis, the acetylated proteins are involved in various cellular processes involving the ribosome, carbon metabolism, antibiotic biosynthesis, secondary metabolites, and the citrate cycle and are distributed in diverse subcellular locations. Additionally, we demonstrated for the first time the acetylation of fatty acid synthase α and β encoded by aflA and aflB involved in the aflatoxin-biosynthesis pathway (cluster 54), as well as backbone enzymes from secondary metabolite clusters 20 and 21 encoded by AFLA_062860 and AFLA_064240, suggesting important roles for acetylation associated with these processes. Our findings illustrating abundant lysine acetylation in A. flavus expand our understanding of the fungal acetylome and provided insight into the regulatory roles of acetylation in secondary metabolism.

Qualitative ubiquitome unveils the potential significances of protein lysine ubiquitination in hyphal growth of Aspergillus nidulans

Protein ubiquitination is an evolutionarily conserved post-translational modification process in eukaryotes, and it plays an important role in many biological processes. Aspergillus nidulans, a model filamentous fungus, contributes to our understanding of cellular physiology, metabolism and genetics, but its ubiquitination is not completely revealed. In this study, the ubiquitination sites in the proteome of A. nidulans were identified using a highly sensitive mass spectrometry combined with immuno-affinity enrichment of the ubiquitinated peptides. The 4816 ubiquitination sites were identified in 1913 ubiquitinated proteins, accounting for 18.1% of total proteins in A. nidulans. Bioinformatic analysis suggested that the ubiquitinated proteins associated with a number of biological functions and displayed various sub-cellular localisations. Meanwhile, seven motifs were revealed from the ubiquitinated peptides, and significantly over-presented in the different pathways. Comparison of the enriched functional catalogues indicated that the ubiquitination functions divergently during growth of A. nidulans and Saccharomyces cerevisiae. Additionally, the proteins in A. nidulans-specific sub-category (cell growth/morphogenesis) were subjected to the protein interaction analysis which demonstrated that ubiquitination is involved in the comprehensive protein interactions. This study presents a first proteomic view of ubiquitination in the filamentous fungus, and provides an initial framework for exploring the physiological roles of ubiquitination in A. nidulans.

A genomic survey of proteases in Aspergilli

BACKGROUND

Proteases can hydrolyze peptides in aqueous environments. This property has made proteases the most important industrial enzymes by taking up about 60% of the total enzyme market. Microorganisms are the main sources for industrial protease production due to their high yield and a wide range of biochemical properties. Several Aspergilli have the ability to produce a variety of proteases, but no comprehensive comparative study has been carried out on protease productivity in this genus so far.

RESULTS

We have performed a combined analysis of comparative genomics, proteomics and enzymology tests on seven Aspergillus species grown on wheat bran and sugar beet pulp. Putative proteases were identified by homology search and Pfam domains. These genes were then clusters based on orthology and extracellular proteases were identified by protein subcellular localization prediction. Proteomics was used to identify the secreted enzymes in the cultures, while protease essays with and without inhibitors were performed to determine the overall protease activity per protease class. All this data was then integrated to compare the protease productivities in Aspergilli.

CONCLUSIONS

Genomes of Aspergillus species contain a similar proportion of protease encoding genes. According to comparative genomics, proteomics and enzymatic experiments serine proteases make up the largest group in the protease spectrum across the species. In general wheat bran gives higher induction of proteases than sugar beet pulp. Interesting differences of protease activity, extracellular enzyme spectrum composition, protein occurrence and abundance were identified for species. By combining in silico and wet-lab experiments, we present the intriguing variety of protease productivity in Aspergilli.

Molecular characterization of the putative transcription factor SebA involved in virulence in Aspergillus fumigatus

Aspergillus fumigatus is a major opportunistic pathogen and allergen of mammals. Nutrient sensing and acquisition mechanisms, as well as the capability to cope with different stressing conditions, are essential for A. fumigatus virulence and survival in the mammalian host. This study characterized the A. fumigatus SebA transcription factor, which is the putative homologue of the factor encoded by Trichoderma atroviride seb1. The ΔsebA mutant demonstrated reduced growth in the presence of paraquat, hydrogen peroxide, CaCl2, and poor nutritional conditions, while viability associated with sebA was also affected by heat shock exposure. Accordingly, SebA::GFP (SebA::green fluorescent protein) was shown to accumulate in the nucleus upon exposure to oxidative stress and heat shock conditions. In addition, genes involved in either the oxidative stress or heat shock response had reduced transcription in the ΔsebA mutant. The A. fumigatus ΔsebA strain was attenuated in virulence in a murine model of invasive pulmonary aspergillosis. Furthermore, killing of the ΔsebA mutant by murine alveolar macrophages was increased compared to killing of the wild-type strain. A. fumigatus SebA plays a complex role, contributing to several stress tolerance pathways and growth under poor nutritional conditions, and seems to be integrated into different stress responses.

Posttranslational modifications of proteins in the pathobiology of medically relevant fungi

Posttranslational modifications of proteins drive a wide variety of cellular processes in eukaryotes, regulating cell growth and division as well as adaptive and developmental processes. With regard to the fungal kingdom, most information about posttranslational modifications has been generated through studies of the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, where, for example, the roles of protein phosphorylation, glycosylation, acetylation, ubiquitination, sumoylation, and neddylation have been dissected. More recently, information has begun to emerge for the medically important fungal pathogens Candida albicans, Aspergillus fumigatus, and Cryptococcus neoformans, highlighting the relevance of posttranslational modifications for virulence. We review the available literature on protein modifications in fungal pathogens, focusing in particular upon the reversible peptide modifications sumoylation, ubiquitination, and neddylation.

Involvement of the Aspergillus nidulans protein kinase C with farnesol tolerance is related to the unfolded protein response

Previously, we demonstrated that the Aspergillus nidulans calC2 mutation in protein kinase C pkcA was able to confer tolerance to farnesol (FOH), an isoprenoid that has been shown to inhibit proliferation and induce apoptosis. Here, we investigate in more detail the role played by A. nidulans pkcA in FOH tolerance. We demonstrate that pkcA overexpression during FOH exposure causes increased cell death. FOH is also able to activate several markers of endoplasmic reticulum (ER) stress and the unfolded protein response (UPR). Our results suggest an intense cross-talk between PkcA and the UPR during FOH-induced cell death. Furthermore, the overexpression of pkcA increases both mRNA accumulation and metacaspases activity, and there is a genetic interaction between PkcA and the caspase-like protein CasA. Mutant analyses imply that MAP kinases are involved in the signal transduction in response to the effects caused by FOH.

The roles played by Aspergillus nidulans apoptosis-inducing factor (AIF)-like mitochondrial oxidoreductase (AifA) and NADH-ubiquinone oxidoreductases (NdeA-B and NdiA) in farnesol resistance

Farnesol (FOH) is a nonsterol isoprenoid produced by dephosphorylation of farnesyl pyrophosphate, a catabolite of the cholesterol biosynthetic pathway. These isoprenoids inhibit proliferation and induce apoptosis. Here, we show that Aspergillus nidulans AifA encoding the apoptosis-inducing factor (AIF)-like mitochondrial oxidoreductase plays a role in the function of the mitochondrial Complex I. Additionally, we demonstrated that ndeA-B and ndiA encode external and internal alternative NADH dehydrogenases, respectively, that have a function in FOH resistance. When exposed to FOH, the ΔaifA and ΔndeA strains have increased ROS production while ΔndeB, ΔndeA ΔndeB, and ΔndiA mutant strains showed the same ROS accumulation than in the absence of FOH. We observed several compensatory mechanisms affecting the differential survival of these mutants to FOH.

Functional characterization of the Aspergillus nidulans methionine sulfoxide reductases (msrA and msrB)

Proteins are subject to modification by reactive oxygen species (ROS), and oxidation of specific amino acid residues can impair their biological function, leading to an alteration in cellular homeostasis. Sulfur-containing amino acids as methionine are the most vulnerable to oxidation by ROS, resulting in the formation of methionine sulfoxide [Met(O)] residues. This modification can be repaired by methionine sulfoxide reductases (Msr). Two distinct classes of these enzymes, MsrA and MsrB, which selectively reduce the two methionine sulfoxide epimers, methionine-S-sulfoxide and methionine-R-sulfoxide, respectively, are found in virtually all organisms. Here, we describe the homologs of methionine sulfoxide reductases, msrA and msrB, in the filamentous fungus Aspergillus nidulans. Both single and double inactivation mutants were viable, but more sensitive to oxidative stress agents as hydrogen peroxide, paraquat, and ultraviolet light. These strains also accumulated more carbonylated proteins when exposed to hydrogen peroxide indicating that MsrA and MsrB are active players in the protection of the cellular proteins from oxidative stress damage.

Transcriptome analysis of the Aspergillus nidulans AtmA (ATM, Ataxia-Telangiectasia mutated) null mutant

ATM is a phosphatidyl-3-kinase-related protein kinase that functions as a central regulator of DNA damage response in eukaryotes. In humans, mutations in ATM cause the devastating neurodegenerative disease Ataxia-Telangiectasia. Previously, we characterized the homologue of ATM (AtmA) in the filamentous fungus Aspergillus nidulans. In addition to its expected role in the DNA damage response, we found that AtmA is also required for polarized hyphal growth. Our results suggested that AtmA probably regulates the function and/or localization of landmark proteins required for the formation of a polarity axis. Here, we extended these studies by investigating which pathways are influenced by AtmA during proliferation and polar growth by comparatively determining the transcriptional profile of A. nidulans wild-type and DeltaatmA mutant strains in different growth conditions. Our results indicate an important role of the pentose phosphate pathway in the fungal proliferation during endogenous DNA damage and polar growth monitored by the AtmA kinase. Furthermore, we identified several genes that have decreased mRNA expression in the DeltaatmA mutant that are involved in the formation of a polarized hyphae and control of polar growth; in the synthesis of phosphatidic acid (e.g. phospholipase D); in the ergosterol biosynthesis (plasma membrane microdomains, lipid rafts); and in intracellular trafficking.

Functional characterization of the putative Aspergillus nidulans DNA damage binding protein homologue DdbA

Nucleotide excision repair (NER) eliminates helix-distorting DNA base lesions. Seven XP-deficient genetic complementation groups (XPA to XPG) have already been identified in mammals, and their corresponding genes have been cloned. Hereditary defects in NER are associated with several diseases, including xeroderma pigmentosum (XP). UV-DDB (XPE) is formed by two associated subunits, DDB1 and DDB2. UV-DDB was identified biochemically as a protein factor that exhibits very strong and specific binding to ultraviolet (UV)-treated DNA. As a preliminary step to characterize the components of the NER in the filamentous fungus Aspergillus nidulans, here we identified a putative DDB1 homologue, DdbA. Deletion and expression analysis indicated that A. nidulans ddbA gene is involved in the DNA damage response, more specifically in the UV light response and 4-nitroquinoline oxide (4-NQO) sensitivity. Furthermore, the DeltaddbA strain cannot self-cross and expression analysis showed that ddbA can be induced by oxidative stress and is developmentally regulated in both asexual and sexual processes. The DeltaddbA mutation can genetically interact with uvsB (ATR), atmA(ATM), nkuA (KU70), H2AX-S129A (a replacement of the conserved serine in the C-terminal of H2AX with alanine), and cshB (a mutation in CSB Cockayne's syndrome protein involved in the transcription-coupled repair subpathway of NER) mutations. Finally, to determine the DdbA cellular localization, we constructed a GFP::DdbA strain. In the presence and absence of DNA damage, DdbA was mostly detected in the nuclei, indicating that DdbA localizes to nuclei and its cellular localization is not affected by the cellular response to DNA damage induced by 4-NQO and UV light.

Purification and characterization of a nitrilase from Aspergillus niger K10

Aspergillus niger K10 cultivated on 2-cyanopyridine produced high levels of an intracellular nitrilase, which was partially purified (18.6-fold) with a 24% yield. The N-terminal amino acid sequence of the enzyme was highly homologous with that of a putative nitrilase from Aspergillus fumigatus Af293. The enzyme was copurified with two proteins, the N-terminal amino acid sequences of which revealed high homology with those of hsp60 and an ubiquitin-conjugating enzyme. The nitrilase exhibited maximum activity (91.6 U mg(-1)) at 45 degrees C and pH 8.0. Its preferred substrates, in the descending order, were 4-cyanopyridine, benzonitrile, 1,4-dicyanobenzene, thiophen-2-acetonitrile, 3-chlorobenzonitrile, 3-cyanopyridine, and 4-chlorobenzonitrile. Formation of amides as by-products was most intensive, in the descending order, for 2-cyanopyridine, 4-chlorobenzonitrile, 4-cyanopyridine, and 1,4-dicyanobenzene. The enzyme stability was markedly improved in the presence of D: -sorbitol or xylitol (20% w/v each). p-Hydroxymercuribenzoate and heavy metal ions were the most powerful inhibitors of the enzyme.

Exposure of Lemna minor to arsenite: expression levels of the components and intermediates of the ubiquitin/proteasome pathway

In animal cells, arsenite has been reported to cause sulfhydryl depletion, generate reactive oxygen species and increase the level of large ubiquitin-protein conjugates. Plant viability tests and DNA laddering experiments have shown that Lemna minor remains viable after exposure to 50 microM NaAsO(2) for periods of at least 6 h. However, protein metabolism is affected in two major ways: the synthesis of an array of stress proteins, which confer thermotolerance; and an increase in the amount of large ubiquitin-protein conjugates, particularly evident after 2-3 h of stress, indicative of a role for the ubiquitin/proteasome pathway. This outcome is primarily attributed to an increased availability of protein substrates during arsenite treatment for three main reasons: an increase in protein carbonyl content after 1-2 h of stress; moderate increments in the transcript levels of the sequences coding for the ubiquitin pathway components chosen as markers (polyubiquitin, E1 and E2, and the beta subunit and the ATPase subunits of the 26S proteasome); the observed increase in ubiquitin conjugates does not depend on de novo protein synthesis. This study is the first report on the involvement of the ubiquitin/proteasome pathway in response to arsenite in plants. In addition, it addresses the simultaneous expression of selected genes encoding the various components of the pathway. The results suggest that in plants, unlike in animals, the response to a relatively low level of arsenite does not induce apoptotic cell death. As a whole, the response to arsenite apparently involves a conjugation of salvage and proteolytic machineries, including heat shock protein synthesis and the ubiquitin/proteasome pathway.

Analysis of expressed sequence tag data and gene expression profiles involved in conidial germination of Fusarium oxysporum

We obtained 3,372 tentative unique transcripts (TUTs) from a cDNA library of Fusarium oxysporum. A cDNA array with 3,158 TUTs was produced to analyze gene expression profiles in conidial germination. It seems that ras and other signaling genes, e.g., ccg, cooperatively initiate conidial germination in Fusarium by increasing protein synthesis.

Ubiquitin depletion as a key mediator of toxicity by translational inhibitors

Cycloheximide acts at the large subunit of the ribosome to inhibit translation. Here we report that ubiquitin levels are critical for the survival of Saccharomyces cerevisiae cells in the presence of cycloheximide: ubiquitin overexpression confers resistance to cycloheximide, while a reduced ubiquitin level confers sensitivity. Consistent with these findings, ubiquitin is unstable in yeast (t(1/2) = 2 h) and is rapidly depleted upon cycloheximide treatment. Cycloheximide does not noticeably enhance ubiquitin turnover, but serves principally to block ubiquitin synthesis. Cycloheximide also induces UBI4, the polyubiquitin gene. The cycloheximide-resistant phenotype of ubiquitin overexpressors is also characteristic of partial-loss-of-function proteasome mutants. Ubiquitin is stabilized in these mutants, which may account for their cycloheximide resistance. Previous studies have reported that ubiquitin is destabilized in the absence of Ubp6, a proteasome-associated deubiquitinating enzyme, and that ubp6 mutants are hypersensitive to cycloheximide. Consistent with the model that cycloheximide-treated cells are ubiquitin deficient, the cycloheximide sensitivity of ubp6 mutants can be rescued either by ubiquitin overexpression or by mutations in proteasome subunit genes. These results also show that ubiquitin wasting in ubp6 mutants is proteasome mediated. Ubiquitin overexpression rescued cells from additional translational inhibitors such as anisomycin and hygromycin B, suggesting that ubiquitin depletion may constitute a widespread mechanism for the toxicity of translational inhibitors.

Transcriptome characterization of the dimorphic and pathogenic fungus Paracoccidioides brasiliensis by EST analysis

Paracoccidioides brasiliensis is a pathogenic fungus that undergoes a temperature-dependent cell morphology change from mycelium (22 degrees C) to yeast (36 degrees C). It is assumed that this morphological transition correlates with the infection of the human host. Our goal was to identify genes expressed in the mycelium (M) and yeast (Y) forms by EST sequencing in order to generate a partial map of the fungus transcriptome. Individual EST sequences were clustered by the CAP3 program and annotated using Blastx similarity analysis and InterPro Scan. Three different databases, GenBank nr, COG (clusters of orthologous groups) and GO (gene ontology) were used for annotation. A total of 3,938 (Y = 1,654 and M = 2,274) ESTs were sequenced and clustered into 597 contigs and 1,563 singlets, making up a total of 2,160 genes, which possibly represent one-quarter of the complete gene repertoire in P. brasiliensis. From this total, 1,040 were successfully annotated and 894 could be classified in 18 functional COG categories as follows: cellular metabolism (44%); information storage and processing (25%); cellular processes-cell division, posttranslational modifications, among others (19%); and genes of unknown functions (12%). Computer analysis enabled us to identify some genes potentially involved in the dimorphic transition and drug resistance. Furthermore, computer subtraction analysis revealed several genes possibly expressed in stage-specific forms of P. brasiliensis. Further analysis of these genes may provide new insights into the pathology and differentiation of P. brasiliensis.

Osmotic stress-coupled maintenance of polar growth in Aspergillus nidulans

Free-living cells monitor extracellular 'osmotic strength' and respond metabolically to offset unfavourable osmotic intracellular solute concentrations. Here, we report the reconstruction of the Aspergillus nidulans salt stress-controlling MAP kinase pathway, based on homology analysis with known yeast genes. In A. nidulans, salt stress HOG genes, such as pbsA, hogA, ptpA and msnA, are upregulated when exposed to high concentrations of salt and, in a hogA deletion mutant (SIK1), the accumulation of pbsA is strongly reduced, suggesting a salt-specific feedback induction mechanism. Growth of SIK1 appears to be unchanged in unstressed cells, but hyphal extension rates are reduced by as much as 60% in the presence of salt. Microscopic observation revealed abnormal hyperbranched hyphal tips, disproportionate accumulation of nuclei and absence of septa. Thus, the inability to maintain turgor pressure depresses cell expansion and results in slower volume increases. In addition, SIK1 fails to partition the apical cell; thus, nuclei are not likely to arrest mitosis in interphase as in normal cells, but continue to divide, accumulating to high levels.

Expression of ubiquitin gene in Microsporum canis and Trichophyton mentagrophytes cultured with fluconazole

The expression of the ubiquitin (Ub) gene in dermatophytes was examined for its relation to resistance against the antifungal drug fluconazole. The nucleotide sequences and the deduced amino acid sequences of the Ub gene in Microsporum canis were proven to be 99% similar to those of the Ub gene in Trichophyton mentagrophytes. Expression of mRNA of Ub in M. canis and T. mentagrophytes was enhanced when the fungi were cultured with fluconazole. The antifungal activity of fluconazole against these dermatophytes was increased in the presence of Ub proteasome inhibitor.