A drastic improvement in the analysis of gene expression has lead to new discoveries in bioinformatics research. In order to analyse the gene expression data, fuzzy clustering algorithms are widely used. However, the resulting analyses from these specific types of algorithms may lead to confusion in hypotheses with regard to the suggestion of dominant function for genes of interest. Besides that, the current fuzzy clustering algorithms do not conduct a thorough analysis of genes with low membership values. Therefore, we present a novel computational framework called the "multi-stage filtering-Clustering Functional Annotation" (msf-CluFA) for clustering gene expression data. The framework consists of four components: fuzzy c-means clustering (msf-CluFA-0), achieving dominant cluster (msf-CluFA-1), improving confidence level (msf-CluFA-2) and combination of msf-CluFA-0, msf-CluFA-1 and msf-CluFA-2 (msf-CluFA-3). By employing double filtering in msf-CluFA-1 and apriori algorithms in msf-CluFA-2, our new framework is capable of determining the dominant clusters and improving the confidence level of genes with lower membership values by means of which the unknown genes can be predicted.
The protease activity of different isolates of dermatophytes representing different ecological groups namely geophilic, zoopahilic and anthropophilic, in their vegetative and sporulation growth phases were compared. Unlike their geophilic and zoophilic counterparts, all the isolates of anthropophilic dermatophytes viz. Trichophyton rubrum, T. mentagrophytes, T. tonsurans, T. violaceum and Epidermophyton floccosum recorded reduced protease activity during artificially induced sporulation phase in comparison to their vegetative growth phase. Even among the anthropophilic group, a classical moderation of protease activity was recorded in Trichyphyton rubrum which also correlates to its clinical manifestation. This enzyme moderation could also be an evolutionary adaptation of the anthropization of these species.
The psychrophilic yeast Glaciozyma antarctica demonstrated high antifreeze activity in its culture filtrate. The culture filtrate exhibited both thermal hysteresis (TH) and ice recrystallization inhibition (RI) properties. The TH of 0.1 °C was comparable to that previously reported for bacteria and fungi. A genome sequence survey of the G. antarctica genome identified a novel antifreeze protein gene. The cDNA encoded a 177 amino acid protein with 30 % similarity to a fungal antifreeze protein from Typhula ishikariensis. The expression levels of AFP1 were quantified via real time-quantitative polymerase chain reaction (RT-qPCR), and the highest expression levels were detected within 6 h of growth at -12 °C. The cDNA of the antifreeze protein was cloned into an Escherichia coli expression system. Expression of recombinant Afp1 in E. coli resulted in the formation of inclusion bodies that were subsequently denatured by treatment with urea and allowed to refold in vitro. Activity assays of the recombinant Afp1 confirmed the antifreeze protein properties with a high TH value of 0.08 °C.
A lectin gene from the Tiger Milk Mushroom Lignosus rhinocerus TM02® was successfully cloned and expressed via vector pET28a in Escherichia coli BL21(DE3). The recombinant lectin, Rhinocelectin, with a predicted molecular mass of 22.8 kDa, was overexpressed in water-soluble form without signal peptide and purified via native affinity chromatography Ni-NTA agarose. Blast protein analysis indicated the lectin to be homologous to jacalin-related plant lectin. In its native form, Rhinocelectin exists as a homo-tetramer predicted with four chains of identical proteins consisting of 11 beta-sheet structures with only one alpha-helix structure. The antiproliferative activity of the Rhinocelectin against human cancer cell lines was concentration dependent and selective. The IC50 values against triple negative breast cancer cell lines MDA-MB-231 and breast cancer MCF-7 are 36.52 ± 13.55 μg mL-1 and 53.11 ± 22.30 μg mL-1 , respectively. Rhinocelectin is only mildly cytotoxic against the corresponding human nontumorigenic breast cell line 184B5 with IC50 value at 142.19 ± 36.34 μg mL-1 . The IC50 against human lung cancer cell line A549 cells is 46.14 ± 7.42 μg mL-1 while against nontumorigenic lung cell line NL20 is 41.33 ± 7.43 μg mL-1 . The standard anticancer drug, Doxorubicin exhibited IC50 values mostly below 1 μg mL-1 for the cell lines tested. Flow cytometry analysis showed the treated breast cancer cells were arrested at G0/G1 phase and apoptosis induced. Rhinocelectin agglutinated rat and rabbit erythrocytes at a minimal concentration of 3.125 μg mL-1 and 6.250 μg mL-1 , respectively.
The aims of the present study were to determine whether Allium sativum (garlic) extract has any effect on the morphology transformation of Candida albicans, and to investigate whether it could alter the gene expression level of SIR2, a morphogenetic control gene and SAP4, a gene encoding secreted aspartyl proteinase.
Lipid biosynthesis produces glycerol, which is important in fueling turgor pressure necessary for germination and penetration of plant host by fungi. As the relationship between pathogenicity and the lipid biosynthetic pathway is not fully understood, we have elucidated the role of the fatty acid synthase beta subunit dehydratase (FAS1) gene in lipid biosynthesis. The FAS1 gene was silenced through homologous double crossover in Magnaporthe oryzae strain S6 to study the effect on lipid biosynthesis. The vegetative growth of Δfas1 mutants show the highest drop on oleic acid (between 10 and 50%), while the mycelial dry weight of mutants dropped significantly on all media. Conidiation of FAS1 mutants show a ~10- and ~5-fold reduction on oatmeal and Potato Dextrose Agar (PDA), respectively. Mutants formed mycelium that were mildly pigmented, indicating that the deletion of FAS1 may have affected melanin biosynthesis. Biochemical and gene expression studies concluded that the fatty acid degradation pathway might have been interrupted by FAS1 deletion. FAS1 mutants showed no enzyme activity on glucose or olive oil, suggesting that the mutants may lack functional peroxisomes and be defective in β-oxidation of fatty acids, hence explaining the reduced lipid deposits in the spores.
Eukaryotic cells typically respond to stress conditions by inhibiting global protein synthesis. The initiation phase is the main target of regulation and represents a key control point for eukaryotic gene expression. In Saccharomyces cerevisiae and mammalian cells this is achieved by phosphorylation of eukaryotic initiation factor 2 (eIF2α). We have examined how the fungal pathogen Candida albicans responds to oxidative stress conditions and show that oxidants including hydrogen peroxide, the heavy metal cadmium and the thiol oxidant diamide inhibit translation initiation. The inhibition in response to hydrogen peroxide and cadmium largely depends on phosphorylation of eIF2α since minimal inhibition is observed in a gcn2 mutant. In contrast, translation initiation is inhibited in a Gcn2-independent manner in response to diamide. Our data indicate that all three oxidants inhibit growth of C. albicans in a dose-dependent manner, however, loss of GCN2 does not improve growth in the presence of hydrogen peroxide or cadmium. Examination of translational activity indicates that these oxidants inhibit translation at a post-initiation phase which may account for the growth inhibition in a gcn2 mutant. As well as inhibiting global translation initiation, phosphorylation of eIF2α also enhances expression of the GCN4 mRNA in yeast via a well-known translational control mechanism. We show that C. albicans GCN4 is similarly induced in response to oxidative stress conditions and Gcn4 is specifically required for hydrogen peroxide tolerance. Thus, the response of C. albicans to oxidative stress is mediated by oxidant-specific regulation of translation initiation and we discuss our findings in comparison to other eukaryotes including the yeast S. cerevisiae.
Proper management of agricultural disease is important to ensure sustainable food security. Staple food crops like rice, wheat, cereals, and other cash crops hold great export value for countries. Ensuring proper supply is critical; hence any biotic or abiotic factors contributing to the shortfall in yield of these crops should be alleviated. Rhizoctonia solani is a major biotic factor that results in yield losses in many agriculturally important crops. This paper focuses on genome informatics of our Malaysian Draft R. solani AG1-IA, and the comparative genomics (inter- and intra- AG) with four AGs including China AG1-IA (AG1-IA_KB317705.1), AG1-IB, AG3, and AG8. The genomic content of repeat elements, transposable elements (TEs), syntenic genomic blocks, functions of protein-coding genes as well as core orthologous genic information that underlies R. solani's pathogenicity strategy were investigated. Our analyses show that all studied AGs have low content and varying profiles of TEs. All AGs were dominant for Class I TE, much like other basidiomycete pathogens. All AGs demonstrate dominance in Glycoside Hydrolase protein-coding gene assignments suggesting its importance in infiltration and infection of host. Our profiling also provides a basis for further investigation on lack of correlation observed between number of pathogenicity and enzyme-related genes with host range. Despite being grouped within the same AG with China AG1-IA, our Draft AG1-IA exhibits differences in terms of protein-coding gene proportions and classifications. This implies that strains from similar AG do not necessarily have to retain similar proportions and classification of TE but must have the necessary arsenal to enable successful infiltration and colonization of host. In a larger perspective, all the studied AGs essentially share core genes that are generally involved in adhesion, penetration, and host colonization. However, the different infiltration strategies will depend on the level of host resilience where this is clearly exhibited by the gene sets encoded for the process of infiltration, infection, and protection from host.
This study was carried out to examine the efficacy of two biocontrol agents (Clonostachys rosea 016, BCA1; Gram-negative bacterium, BCA5) for control of FUM1 gene expression and fumonisin B1 (FB1) production by F. verticillioides FV1 on maize cobs of different ripening stages: R3, Milk (0.985 aw); R4, Dough (0.976 aw); R5, Dent (0.958 aw). Initially, temporal studies on FUM1 gene expression and FB1 production were performed on maize kernels for up to 14days. This revealed that day 10 was optimum for both parameters, and was used in the biocontrol studies. Maize cobs were inoculated with 50:50 mixtures of the pathogen:antagonist inoculum and incubated in environmental chambers to maintain the natural aw conditions for ten days at 25 and 30°C. The growth rates of F. verticillioides FV1, the relative expression of the FUM1 gene and FB1 production were quantified. It was found that, aw×temp had significant impacts on growth, FUM1 gene expression and FB1 production by F. verticillioides FV1 on maize cobs of different maturities. The fungal antagonist (BCA1) significantly reduced FB1 contamination on maize cobs by >70% at 25°C, and almost 60% at 30°C regardless of maize ripening stage. For the bacterial antagonist (BCA5) however, FB1 levels on maize cobs were significantly decreased only in some treatments. These results suggest that efficacy of antagonists to control mycotoxin production in ripening maize cobs needs to take account of the ecophysiology of the pathogen and the antagonists, as well as the physiological status of the maize during silking to ensure effective control.
Squalene synthase (SS) is the key precursor and first committed enzyme of the sterol biosynthesis pathway. In a previous work, SS has been identified as one of the immunogenic proteins that could be a potential diagnostic candidate for the pathogenic fungus Candida tropicalis. In this study, SS from C. tropicalis was cloned and expressed as recombinant protein in Pichia pastoris to investigate its reactivity with serum antibodies. ERG9 gene that encodes for SS was amplified by PCR and cloned in-frame into pPICZB expression vector. The recombinant construct was then transformed into P. pastoris GS115 host strain. Expression of the recombinant protein was confirmed by SDS-PAGE and Western blot analysis using anti-His tag probe. Optimal protein production was achieved by cultivating the culture with 1.0% methanol for 72h. The recombinant protein was purified to approximately 97% pure in a single step immobilized metal affinity chromatography with a yield of 70.3%. Besides, the purified protein exhibited specific reactivity with immune sera on Western blot. This is the first report on heterologous expression of antigenic SS from C. tropicalis in P. pastoris which can be exploited for large-scale production and further research. The results also suggested that the protein might be of great value as antigen candidate for serodiagnosis of Candida infection.
Lignocellulosic biomass dedicated to bioethanol production usually contains pentoses and inhibitory compounds such as furfural that are not well tolerated by Saccharomyces cerevisiae. Thus, S. cerevisiae strains with the capability of utilizing both glucose and xylose in the presence of inhibitors such as furfural are very important in industrial ethanol production. Under the synergistic conditions of transaldolase (TAL) and alcohol dehydrogenase (ADH) overexpression, S. cerevisiae MT8-1X/TAL-ADH was able to produce 1.3-fold and 2.3-fold more ethanol in the presence of 70 mM furfural than a TAL-expressing strain and a control strain, respectively. We also tested the strains' ability by mimicking industrial ethanol production from hemicellulosic hydrolysate containing fermentation inhibitors, and ethanol production was further improved by 16% when using MT8-1X/TAL-ADH compared to the control strain. Transcript analysis further revealed that besides the pentose phosphate pathway genes TKL1 and TAL1, ADH7 was also upregulated in response to furfural stress, which resulted in higher ethanol production compared to the TAL-expressing strain. The improved capability of our modified strain was based on its capacity to more quickly reduce furfural in situ resulting in higher ethanol production. The co-expression of TAL/ADH genes is one crucial strategy to fully utilize undetoxified lignocellulosic hydrolysate, leading to cost-competitive ethanol production.
The aims of our research were to investigate the gene expression of the multidrug efflux transporter, CDR1 and the major drug facilitator superfamily transporter, MDR1 gene in azole drug-resistant Candida albicans and Candida glabrata clinical isolates recovered from vaginitis patients; and to identify hotspot mutations that may be present in the C. albicans CaCDR1 gene that could be associated with drug-resistance. The relative expression of the CDR1 and MDR1 transcripts in ketoconazole and clotrimazole-resistant isolates and drug-susceptible ATCC strains were determined by semi-quantitative reverse transcription-polymerase chain reaction. Expression of CaCDR1 transcript was upregulated to varying extents in all three azole-resistant C. albicans isolates studied (1.6-, 3.7- and 3.9-fold) and all three C. glabrata isolates tested (at 1.9-, 2.3- and 2.7-fold). The overexpression level of CaCDR1 in the isolates correlated with the degree of resistance as reflected by the minimum inhibitory concentration (MIC) of the drugs. The messenger RNA for another efflux pump, MDR1, was also overexpressed in one of the azole-resistant C. albicans isolates that overexpressed CDR1. This finding suggests that drug-resistance may involve synergy between energy-dependent drug efflux pumps CDR1p and MDR1p in some but not all isolates. Interestingly, DNA sequence analysis of the promoter region of the CaCDR1 gene revealed several point mutations in the resistant clinical isolates compared to the susceptible isolates at 39, 49 and 151 nucleotides upstream from the ATG start codon. This finding provides new information on point mutations in the promoter region which may be responsible for the overexpression of CDR1 in drug-resistant isolates.
Candida albicans is capable of undergoing yeast-hypha transition to attain pathogenicity in humans. In this study, we investigated the differential expression of CaSIR2 via quantitative real-time PCR (qPCR), during yeast-hypha transition with and without the presence of 2-dodecanol. SIR2 transcript levels were found to be significantly enhanced after hyphal induction as compared to the yeast form. This study found that 2-dodecanol is able to inhibit hyphal development and block SIR2 up-regulation, even in hyphal-inducing growth conditions. We suggest that SIR2 may be involved in Candida albicans quorum-sensing and serum-induced yeast-hyphae transition via the Ras1-cAMP-Efg1 signalling cascade.
Trichoderma is a genus of soil-borne fungus with an abundance of reports of its economic importance in the agriculture industry. Thus, the correct identification of Trichoderma species is necessary for its commercial purposes. Globally, Trichoderma species are routinely identified from micro-morphological descriptions which can be tedious and prone to errors. Thus, we emphasize that the accurate identification of Trichoderma strains requires a three-pronged approach i.e. based on its morphological characteristics, multilocus gene sequences of the rDNA [internal transcribed spacer (ITS) 1 and 2 regions], translation elongation factor 1-α (TEF-1α), Calmodulin (CAL) and its lignocellulolytic activities. We used this approach to identify a total of 53 Trichoderma strains which were isolated from a wet paddy field located at Tuaran, Sabah, Malaysia. The 53 strains were positively identified as belonging to three Trichoderma species, namely T. asperellum (43 strains), T. harzianum (9 strains), and T. reesei (one strain) on the basis of its morphological characteristics and multilocus gene sequences. Phylogenetic trees constructed based on the UPGMA method of the ITS 1 and 2 regions of the rDNA, TEF-1α and CAL revealed three distinct groups with the T. asperellum, T. harzianum and T. reesei strains placed under the section of Trichoderma, Pachybasium and Longibrachiatum, respectively. In addition, the lignocellulolytic activities of the isolates were measured based on the diameters of the halo zones produced when degrading cellulose, lignin, and starch, respectively. This diagnostic assay can be used to identify Trichoderma as it produces polyphenol oxidase when Tannic Acid Media is used for the lignin test, endoglucanases when Jensen media is used for cellulose, and it hydrolyzes starch to glucose when the modified Melin-Nokrans media is used for the starch test. Accurate identification of Trichoderma species is needed as these strains can potentially be used as a biocontrol agent to prevent diseases and to increase yield in agriculture crops.
Clustering is used widely in 'omics' studies and is often tackled with standard methods, e.g. hierarchical clustering. However, the increasing need for integration of multiple data sets leads to a requirement for clustering methods applicable to mixed data types, where the straightforward application of standard methods is not necessarily the best approach. A particularly common problem involves clustering entities characterized by a mixture of binary data (e.g. presence/absence of mutations, binding, motifs and epigenetic marks) and continuous data (e.g. gene expression, protein abundance, metabolite levels). Here, we present a generic method based on a probabilistic model for clustering this type of data, and illustrate its application to genetic regulation and the clustering of cancer samples. We show that the resulting clusters lead to useful hypotheses: in the case of genetic regulation these concern regulation of groups of genes by specific sets of transcription factors and in the case of cancer samples combinations of gene mutations are related to patterns of gene expression. The clusters have potential mechanistic significance and in the latter case are significantly linked to survival. The method is available as a stand-alone software package (GNU General Public Licence) from http://github.com/BioToolsLeeds/FlexiCoClusteringPackage.git.
The ability of eukaryotes to adapt to an extreme range of temperatures is critically important for survival. Although adaptation to extreme high temperatures is well understood, reflecting the action of molecular chaperones, it is unclear whether these molecules play a role in survival at extremely low temperatures. The recent genome sequencing of the yeast Glaciozyma antarctica, isolated from Antarctic sea ice near Casey Station, provides an opportunity to investigate the role of molecular chaperones in adaptation to cold temperatures. We isolated a G. antarctica homologue of small heat shock protein 20 (HSP20), GaSGT1, and observed that the GaSGT1 mRNA expression in G. antarctica was markedly increased following culture exposure at low temperatures. Additionally, we demonstrated that GaSGT1 overexpression in Escherichia coli protected these bacteria from exposure to both high and low temperatures, which are lethal for growth. The recombinant GaSGT1 retained up to 60 % of its native luciferase activity after exposure to luciferase-denaturing temperatures. These results suggest that GaSGT1 promotes cell thermotolerance and employs molecular chaperone-like activity toward temperature assaults.
The rice blast fungus Magnaporthe oryzae causes plant disease via specialised infection structures called appressoria. These dome-shaped cells are able to generate enormous internal pressure, which enables penetration of rice tissue by invasive hyphae. Previous studies have shown that mobilisation of lipid bodies and subsequent lipid metabolism are essential pre-requisites for successful appressorium-mediated plant infection, which requires autophagic recycling of the contents of germinated spores and germ tubes to the developing appressorium. Here, we set out to identify putative regulators of lipid metabolism in the rice blast fungus. We report the identification of FAR1 and FAR2, which encode highly conserved members of the Zn2-Cys6 family of transcriptional regulators. We generated Δfar1, Δfar2 and Δfar1Δfar2 double mutants in M. oryzae and show that these deletion mutants are deficient in growth on long chain fatty acids. In addition, Δfar2 mutants are also unable to grow on acetate and short chain fatty acids. FAR1 and FAR2 are necessary for differential expression of genes involved in fatty acid β-oxidation, acetyl-CoA translocation, peroxisomal biogenesis, and the glyoxylate cycle in response to the presence of lipids. Furthermore, FAR2 is necessary for expression of genes associated with acetyl-CoA synthesis. Interestingly, Δfar1, Δfar2 and Δfar1Δfar2 mutants show no observable delay or reduction in lipid body mobilisation during plant infection, suggesting that these transcriptional regulators control lipid substrate utilization by the fungus but not the mobilisation of intracellular lipid reserves during infection-related morphogenesis.
The small ubiquitin-like modifier (SUMO) is implicated in various cellular activities, including transcriptional regulation. We previously showed that the yeast activator Gcn4 becomes sumoylated during activation, facilitating its eventual promoter eviction and transcriptional shut off. Here we show that the corepressor Tup1 is sumoylated, at two specific lysines, under various stress conditions. Mutation of these sites has no effect on Tup1 recruitment or RNAP II promoter occupancy immediately following induction. However, Tup1 levels subsequently decrease, while RNAP II and transcription increase in Tup1 mutant cells. Consistent with this, a Tup1 mutant displaying increased sumoylation led to reduced transcription. We also show that coordinated sumoylation of Gcn4 and Tup1 enhances Gcn4 promoter eviction and that multiple Tup1-interacting proteins become sumoylated after stress. Together, our studies provide evidence that coordinated sumoylation of Gcn4, Tup1 and likely other factors dampens activated transcription by stabilizing Tup1 binding and stimulating Gcn4 and RNAP II removal.
Comparative genomic studies have reported widespread variation in levels of gene expression within and between species. Using these data to infer organism-level trait divergence has proven to be a key challenge in the field. We have used a wild Malaysian population of S. cerevisiae as a test bed in the search to predict and validate trait differences based on observations of regulatory variation. Malaysian yeast, when cultured in standard medium, activated regulatory programs that protect cells from the toxic effects of high iron. Malaysian yeast also showed a hyperactive regulatory response during culture in the presence of excess iron and had a unique growth defect in conditions of high iron. Molecular validation experiments pinpointed the iron metabolism factors AFT1, CCC1, and YAP5 as contributors to these molecular and cellular phenotypes; in genome-scale sequence analyses, a suite of iron toxicity response genes showed evidence for rapid protein evolution in Malaysian yeast. Our findings support a model in which iron metabolism has diverged in Malaysian yeast as a consequence of a change in selective pressure, with Malaysian alleles shifting the dynamic range of iron response to low-iron concentrations and weakening resistance to extreme iron toxicity. By dissecting the iron scarcity specialist behavior of Malaysian yeast, our work highlights the power of expression divergence as a signpost for biologically and evolutionarily relevant variation at the organismal level. Interpreting the phenotypic relevance of gene expression variation is one of the primary challenges of modern genomics.
The infection of Cryptococcus neoformans is acquired through the inhalation of desiccated yeast cells and basidiospores originated from the environment, particularly from bird's droppings and decaying wood. Three environmental strains of C. neoformans originated from bird droppings (H4, S48B and S68B) and C. neoformans reference clinical strain (H99) were used for intranasal infection in C57BL/6 mice. We showed that the H99 strain demonstrated higher virulence compared to H4, S48B and S68B strains. To examine if gene expression contributed to the different degree of virulence among these strains, a genome-wide microarray study was performed to inspect the transcriptomic profiles of all four strains. Our results revealed that out of 7,419 genes (22,257 probes) examined, 65 genes were significantly up-or down-regulated in H99 versus H4, S48B and S68B strains. The up-regulated genes in H99 strain include Hydroxymethylglutaryl-CoA synthase (MVA1), Mitochondrial matrix factor 1 (MMF1), Bud-site-selection protein 8 (BUD8), High affinity glucose transporter 3 (SNF3) and Rho GTPase-activating protein 2 (RGA2). Pathway annotation using DAVID bioinformatics resource showed that metal ion binding and sugar transmembrane transporter activity pathways were highly expressed in the H99 strain. We suggest that the genes and pathways identified may possibly play crucial roles in the fungal pathogenesis.