Comparative transcriptome analysis of microsclerotia development in Nomuraea rileyi
© Song et al.; licensee BioMed Central Ltd. 2013
Received: 15 January 2013
Accepted: 7 June 2013
Published: 19 June 2013
Nomuraea rileyi is used as an environmental-friendly biopesticide. However, mass production and commercialization of this organism are limited due to its fastidious growth and sporulation requirements. When cultured in amended medium, we found that N. rileyi could produce microsclerotia bodies, replacing conidiophores as the infectious agent. However, little is known about the genes involved in microsclerotia development. In the present study, the transcriptomes were analyzed using next-generation sequencing technology to find the genes involved in microsclerotia development.
A total of 4.69 Gb of clean nucleotides comprising 32,061 sequences was obtained, and 20,919 sequences were annotated (about 65%). Among the annotated sequences, only 5928 were annotated with 34 gene ontology (GO) functional categories, and 12,778 sequences were mapped to 165 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) database. Furthermore, we assessed the transcriptomic differences between cultures grown in minimal and amended medium. In total, 4808 sequences were found to be differentially expressed; 719 differentially expressed unigenes were assigned to 25 GO classes and 1888 differentially expressed unigenes were assigned to 161 KEGG pathways, including 25 enrichment pathways. Subsequently, we examined the up-regulation or uniquely expressed genes following amended medium treatment, which were also expressed on the enrichment pathway, and found that most of them participated in mediating oxidative stress homeostasis. To elucidate the role of oxidative stress in microsclerotia development, we analyzed the diversification of unigenes using quantitative reverse transcription-PCR (RT-qPCR).
Our findings suggest that oxidative stress occurs during microsclerotia development, along with a broad metabolic activity change. Our data provide the most comprehensive sequence resource available for the study of N. rileyi. We believe that the transcriptome datasets will serve as an important public information platform to accelerate studies on N. rileyi microsclerotia.
KeywordsTranscriptome Oxidative Stress Microsclerotia Nomuraea Rileyi
Nomuraea rileyi is a well-known fungal pathogen thriving in a range of obligate hosts, particularly the noctuids. It affects polyphagous species of Heliothis, Spodoptera, Pseudoplusia, Trichoplusia, Plutella, and Rachiplusia[1, 2], and hence, can be used as a mycoinsecticide. However, the lack of reliable cost-effective protocols for mass production of this entomopathogenic fungus limits its commercialization. The present mass productions methods are not cost effective and require special growth conditions (stimulatory light and maltose) for sporulation . As a result, there is a need to find an active agent for commercial production of N. rileyi. Microsclerotium (MS), with a diameter of 200–600 μm, is a pseudoparenchymatous aggregation of hyphae, comprising only a few cells. It is produced by many phytopathogenic fungi for persistence in the soil and decaying plant material. The structures of fungal phytopathogen, including the infective and spreading structures, help them to survive for long periods in the environment [4, 5]. MS production was successfully induced in phytopathogenic fungus [6, 7] and entomopathogenic fungus Metarhizium anisopliae during submerged liquid culture fermentation. Our laboratory successfully realized MS production by N. rileyi on liquid amended medium (AM), which exhibited insecticidal efficacy on Spodoptera litura, and thus can be used as an active agent (data not shown).
The morphological events that occur during microsclerotia production have been studied extensively [9, 10]. A mature MS has been observed to exhibit three distinct layers: a pigmented rind, a thin-walled cortex, and a large central medulla [11, 12]. During the initial stages of development, the hyphae swell and aggregate. The aggregations subsequently enlarge and the pigment emerges. In the final phase, the pigment gets deposited within the cell wall. The mature microsclerotia can disperse in the environment as persistent infective propagules, affecting a broad range of herbaceous and woody hosts .
The potentially important genes in microsclerotia development have already been investigated. Two expressed sequence tag (EST) libraries developed from different cultures of a tomato isolate strain Verticillium dahliae, were found to exhibit the most unique sequences, of which 55% had protein sequences similar to those presented in the database . Based on the EST databases, VDH1, a hydrophobin gene, was confirmed to be involved in microsclerotia development [15, 16]. Subsequently, using the Agrobacterium tumefaciens-mediated transformation (ATMT) method, the role of the gene in microsclerotia development was identified. In addition to targeted gene disruption by homologous recombination, the ATMT has been successfully exploited for large-scale forward genetic screen to create insertional mutants. A gene that encodes glutamic acid-rich protein and affects microsclerotia formation and pathogenicity was identified by the ATMT . It has been reported that the physiological and morphological differences during microsclerotia development should be, at least in part, due to numerous genes expressed at various development stages [18–21]. Transcript profiling is an important strategy for studying the expression of large genes. Although EST sequencing has been used for the detection of reference transcripts, it has some inherent limitations, such as low throughput, high cost, and long experimental cycle. Recently, next-generation sequencing (NGS) technology has emerged for high-throughput sequence determination, which is quite economical . The RNA-Seq technology has enabled us to investigate the transcriptome for various gene expression studies without reference genome sequences .
In the present study, the NGS technology was used to examine the transcriptome of MS of N. rileyi. We constructed two libraries sequenced using Illumina HiSeq™ 2000 based on N. rileyi cultures grown in liquid AM and minimal medium (MM), which produced abundant MS and mycelium, respectively. Subsequently, we analyzed the expression of 32,061 unigenes, found the differentially expressed genes related to MS development, compared the up-regulation of the expression of unigenes in the AM libraries, and validated some unigenes based on different development stages by using RT-qPCR. The two assembled and annotated transcriptome sequences provide an invaluable resource for the identification of N. rileyi involved in MS development.
Illumina sequencing and reads assembly
Summary of the sequence assembly after Illumina sequencing
Average read length
Mean length of contigs
Mean length of unigenes
Mean length of sequences
A total of 32,061 sequences were assembled, with an average of 933 bp and an N50 of 1522 bp (Table 1). The size distribution of those sequences is shown in Additional file 1: Figure S1. Among the sequences, 19,788 were longer than 500 bp, 11,626 were longer than 1000 bp, and 3896 were longer than 2000 bp (Additional file 1: Figure S1B).
All-unigenes annotation provides functional information, including protein sequence similarities, clusters of orthologous groups (COG), gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) pathway information. For annotation, the sequences were compared with the protein databases (NCBI Nr, SwissProt, KEGG, and COG) using BLASTX (E-value <10−5), and the protein functions were predicted from the annotations of the most similar proteins. The sequences were first searched using BLASTX against the NCBI Nr database, and 20,919 sequences (65.25%) had hits that exceeded the E-value cutoff (Additional file 2: Figure S2). A total of 13,385 sequences (41.75%) were identified from the SwissProt database.
The 20,919 nr hits were mapped to the reference canonical pathways in KEGG. In total, 12,778 sequences were assigned to 165 KEGG pathways. To identify the biological pathways that are active in different ways among the two treatments, a total of 1888 differentially expressed sequences was assigned to 161 KEGG pathways. With significant enrichment pathways, we could discern the main biochemical and signal transduction pathways. After multiple testing corrections, we chose pathways with Qvalue ≤ 0.05 as significantly enriched among the differentially expressed genes. In total, 26 significant enrichment pathways were detected. The enrichment pathways that were most represented were the biosynthesis of secondary metabolites (305 members), metabolic pathways (721 members), and microbial metabolism in diverse environments (300 members). These annotations are valuable resource for the comparison of processes, functions, and pathways in MS research.
Detection of upregulated genes related to MS development
Following AM and MM treatments, a total of 4808 differentially expressed unigenes were identified. Among the 4808 unigenes, 1538 were up-regulated and 322 were induced to express during AM treatment. Sequence analyses revealed that out of the 1538 unigenes, 691 (44.93%) showed no homology with known sequences and were classified as “no hits”. Among the 322 unigenes, 165 (51.24%) were found to have known functions. We identified a number of unigenes that were induced during AM treatment as well as expressed in the enrichment pathways, such as bifunctional P-450:NADPH-P450 reductase (unigene24005), FAD binding domain protein (unigene23494), glyceraldehyde 3-phosphate dehydrogenase (unigene23268), heat-shock protein (unigene23729, unigene23337, unigene22930, and unigene22875), SUR2-type hydroxylase/desaturase (unigene23801), delta-9 desaturase isoenzyme B (unigene23839), acyl-CoA desaturase (unigene23046), and NADP-dependent mannitol dehydrogenase (unigene23414). Although these unigenes were noted to play different roles in various metabolisms, they were found to be commonly involved in oxidative stress detoxification. This indicates that oxidative stress may occur during MS development [27, 28].
To confirm the results of the Solexa/Illumina sequencing, the following 17 unigenes, up- regulated or uniquely expressed following AM treatment, were selected for RT- qPCR assays: superoxide dismutase (SOD; unigene5568), catalase (CAT; unigene2410), heat-shock protein (SSC1; unigene23729), glutathione synthase (GS; unigene5050), glutathione S transferase (GSTs; unigene13979), glutathione reductase (GR; unigene22615), polyketone synthase (Pks; unigene22495), delta-9 desaturase isoenzyme B (DB; unigene23839), 2OG-Fe (II) oxygenase family oxidoreductase (FO; unigene1892), pyruvate carboxylase (PYC; unigene3898), pyruvate decarboxylase (PDC; unigene9437), aspartate aminotransferase (AST; unigene12484), alanine aminotransferase (MAC; unigene9129), aspartate-tRNA ligase (ARS; unigene18089), prolyl-tRNA synthase (PRS; unigene16762), acetoacetyl-CoA synthase (ACS; unigene22462), and vacuolar ATP synthase catalytic subunit A (ATP-synt A; unigene18015).
High-throughput mRNA sequencing technology is particularly suitable for gene expression profiling in non-model organisms without prior genome annotation . In the present study, we used NGS technology, based on the Illumina HiSeq™ 2000 platform, obtained 9.2 Gb of coverage with 52,200,794 clean sequencing reads, and identified 32,061 all-unigenes from two libraries with a mean length of 933. Among them, 20,919 (65.25%) unigenes were successfully annotated. Among the two transcriptomes, approximately 4800 unigenes were found to be differentially expressed. The results of the present study expand on those reported in microsclerotia studies based on EST sequencing, and provide more comprehensive genetic and genomic information [14, 29].
During the undifferentiated (UD) stage of MS, similar to other phytopathogenic fungi, there is abundant hyphal growth, accompanied by material and energy metabolism. In the present study, numerous related genes were found to be up-regulated during this stage (Figure 3). Furthermore, the natural byproducts of mitochondrial respiration and oxygen consumption were noted to be converted to ROS, such as superoxide anion radical (O2−). This is the initial step in the formation of ROS by donating a single electron to the oxygen molecule (O2), as well as a product of important enzymes such as xanthine oxidase (XO) and NADPH oxidase [30, 31]. However, O2− radical is not highly oxidizing, and subsequently produces H2O2 to exert its toxic role. H2O2 can capture an electron from another O2− or free ferrous (Fe2+) and cuprous ions, forming hydroxyl radical (OH), which is highly reactive and causes lipid peroxidation; one of the major end products of lipid peroxidation had been shown to be involved in signal transduction . Furthermore, H2O2 is stable and can cross cell membrane. Thus, although H2O2 is a harmful metabolic byproduct, it also functions as a ubiquitous intracellular messenger. Different levels of H2O2 can induce distinct responses within a cell, especially for the proper development and proliferation of cells [25, 33, 34].
Enzymatic defenses against O2− include SOD that dismute O2− to H2O2 in fungi , and the organisms defend against H2O2 mainly by CAT to catalytically decompose it to H2O and O2. In the present study, unigene5568, annotated as sod, was highly expressed in the AM library. In particular, the expression level increased between the development and maturation stages (Figure 3). An earlier study also reported a similar finding during MS development in Sclerotinia sclerotiorum. In the present study, maximum consumption of O2 was observed at the MS initiation (SI) period (data not shown), demonstrating that O2− was continuously produced and the alteration occurred between the UD and SI period. The unigenes (cat, unigene11956, unigene13044, unigene13045, and unigene19140) were highly expressed in the AM library, and their expression was the highest during the SI period (Figure 3). Among the basal salts in the AM culture, Fe2+ was noted to promote MS biogenesis (data not shown), which could catalyze the formation of OH. radicals. In the present study, the expression of unigene1892 (fo) was not only high in the AM library, but was also relatively highest during the SI period (Figure 3). Thus, the generation of ROS may play an important role in the induction of MS formation, as observed in other filamentous fungi [38, 39]. Furthermore, enzymes involved in the detoxification of free radicals also appear to have important role during MS development.
Generation of ROS has forced many organisms to devise complex antioxidant defense, such as thioredoxin (TRX), glyceraldehyde 3-phosphate dehydrogenase (G3PDH), and reduced glutathione (GSH). The GSH, in particular, has been shown to be the main non-enzymatic antioxidant defense, helping to maintain the reduced environment of cytosol, and can be used as a substrate to detoxify the cell from oxidant xenobiotic molecules using GSTs to form oxidized glutathione (GSSG) . GR is a ubiquitous flavor enzyme of disulfide reductase family, which catalyzes the NADPH-dependent reduction of GSSG to GSH. These enzymes maintain the balance of the redox couple (GSSG/GSH). The expression of unigenes (gs, gr, and gsts) of functional redox systems involved in the balance of thiol and oxydic glutathione was high in the AM library, and highest during the SI period (Figure 3).This shows the increase in the oxidative stress associated reduced thiol when the undifferentiated mycelium enters the differentiated state , which can combat the incremental oxidative stress during MS development. Other stress regulators, such as heat-shock proteins, play important role during oxidative stress, and have well-established roles during different stress conditions in the organism and are up-regulated during stress [41, 42]. The unigene (ssc1) showed prominent multiple peaks in RT-qPCR and illustrated regulatory activity under oxidative stress (Figure 3).
ROS are not only inevitable byproducts of oxygen metabolism, but also play a role in cellular signaling, which mediates or augments the effects of growth factors, cell differentiation, and programmed cell death . Based on the AM library, the dynamic activity of MS was found to be cell differentiation during the SI period (Figure 1), and the doses of H2O2 were presumed to reach a certain level (Figure 3) or “sub-poison” concentration [33, 34]. The levels of ROS accumulated during stress condition and their subcellular source were noted to determine the expression pattern of specific genes and induction of stress responsive pathways . Much of this response appeared to be designed to decrease energy production and biochemical processes. In this regard, the TCA cycle was reported to be a significant source of reducing equivalents. The regulation of the two genes, pyc and pdc, related to pyruvate carboxylation, was low (5 and 13 folds, respectively; Figure 3), and the induction of genes of amino acid and fatty acid metabolism was 10 and 30 folds (mac, ast, acs), respectively. These findings indicate that low-level metabolism may be important to balance the generation of proton. In addition, several other genes (prs, ars, ATP-synt A) were also significantly down-regulated during this period (Figure 3).
The pigment granules were found to be deposited on the rind during Nomuraea rileyi MS formation, but not during melanin production [6–8]. An earlier study examined the biosynthetic pathway of the pigment in other fungi . In our transcriptome, the unigene22495, identified as pks, was observed to be involved in the biosynthesis of the pigment. The abundance of this gene increased significantly between the SI and MS development (SM) period (Figure 3), and the change in genotype was in accordance with the phenotype. It has been demonstrated that pigment production can be uncoupled from other morphological changes associated with MS development [14, 46]. Two unigenes (acs and db) were noted to be involved in the metabolism of lipid acid. The db gene increased significantly during the SI period (Figure 3), indicating that the unsaturated lipid acid participates in resisting oxidative stress. Similarly, the acs gene was maintained at high levels during the UD period and remained at certain expression level during other periods (Figure 3), implying that the gene not only takes part in lipid metabolism, but is also related to the biosynthetic pathway of the pigment .
In the present study, a large number of differentially expressed genes involved in MS development were found. We analyzed the expression pattern of these differentially expressed genes, particularly, their up-regulation and unique unigenes in the AM libraries. These unigenes were found to primarily code for stress-responsive genes or a part of the cell machinery underlying the biogenesis of various bodies; however, numerous genes with unknown functions were also noted .These results indicate that oxidative stress takes place during MS differentiation, and the signaling genes should be further identified. Overall, this study suggests that fungi employ various pathways and regulatory networks of genes in response to different culture conditions, and demonstrated at the genetic level that MS differentiation in N. rileyi is related to high oxidative stress.
Using the Illumina sequencing, we surveyed the two transcriptomes of N. rileyi under different culture conditions, assembled 32,061 unigenes, annotated 20,919 of those unigenes, and found 4808 differentially expressed unigenes among the two libraries. Our findings substantially contribute to the existing sequence resources for N. rileyi, and will certainly accelerate the research on MS development. In addition, the de novo transcriptome analysis of non-model organism without prior genome annotation has again been successfully demonstrated in this study. The genes data will provide a portfolio of candidate genes for further research on the gene expression, genomics, and functional genomics of N. rileyi.
Strain and growth conditions
N. rileyi strain was obtained from the Genetic Engineering Research Center, Chongqing, China. The strain was grown on Sabouraud maltose agar fortified with 1% yeast extract (SMAY) for 14 days, under continuous light at 25°C. The conidia suspension, suspended in sterile deionized water with 0.5% Tween 80 (Sigma-Aldrich, America) at a concentration of 1 × 108 conidia/mL, was inoculated in a 100-mL liquid medium (in flask) and incubated in a shaker at 25°C and 200 rpm. The MM (consisting of 40 g/L of glucose, 2.5 g/L of peptone, and 5 g/L of yeast extract) was designed to reflect the nutritional conditions of hypha development. The AM (comprising 40 g/L of glucose, 2.5 g/L of peptone, 5 g/L of yeast extract, 4.0 g/L of KH2PO4, 0.8 g/L of CaCl2.2H2O, 0.6 g/L of MgSO4.7H2O, 0.1 g/L of FeSO4.7H2O, 37 mg/L of CoCl2.6H2O, 16 mg/L of MnSO4.H2O, and 14 mg/L of ZnSO4.7H2O) contained substances that could induce the development of MS. The development stages, including UD, SI, MS development (SD), and SM, of N. rileyi were confirmed in 2, 3.5, 4, and 5 days, respectively. Two samples, separately cultured in the liquid MM and AM for 3.5 days, were collected using Whitman filter paper. The mycelia tissues were separated from MS using forceps and washed twice with distilled autoclaved water. Subsequently, the two samples were frozen at −80°C until further use.
RNA isolation and library preparation for transcriptome analysis
Total RNA was isolated using the TRIzol® reagent (Invitrogen, USA). The RNA samples were then digested using DNase I for 30 min at 37°C to remove potential genomic DNA. RNA integrity was confirmed using the 2100 Bioanalyzer. The mRNA was purified from total RNA using oligo (dT) magnetic beads. Subsequently, the mRNA was fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were used for first-strand cDNA synthesis by employing random primers. The second-strand cDNA synthesis was carried out using DNA polymerase I and RNase H. These cDNA fragments were subjected to an end repair process and ligation of adapters. Subsequently, these products were enriched by PCR to create the final libraries (AM and MM).
Analysis of Illumina sequencing results
The sequencing library was sequenced using the HiSeqTM2000 platform (Illumina). Over 26 million sequence raw reads were generated for the analysis. The adaptor fragments were removed from the raw reads to yield the clean reads required for the analysis. De novo transcriptome assembly of these short reads was carried out with short reads assembling program –Trinity . Trinity first combined the reads with certain length of overlap to form longer fragments without N, using “N” to represent unknown sequences. These contigs were utilized for further process of cluster using sequence clustering software to form longer sequences without N . These sequences were defined as unigenes. Unigenes from each sample’s assembly were subjected to further sequence splicing and redundancy removal using sequence clustering software to acquire non-redundant unigenes as long as possible. These were designated as all-unigenes.
For further analysis, we used BLASTX (E-value < 10−5) to compare between unigenes and protein databases such as NCBI nr, Swiss-Prot, KEGG, and COG, and the best aligning results were used to decide the sequence direction of the unigenes. Unigene annotation provided information regarding the expression and functional annotation of the identified genes. With nr annotation, we used Blast2GO program to obtain GO annotation of the unigenes. After acquiring the GO annotation for every unigene, we used WEGO software to carry out GO functional classification for all-unigenes and understand the distribution of gene functions of the species at the macro level. The COG and KEGG pathways annotation was carried out using Blastall software against COG (http://www.ncbi.nlm.nih.gov/COG) and KEGG (http://www.genome.jp/kegg/) database .
Gene validation and expression analysis
One differentially expressed gene with potential roles in MS development was chosen for validation using RT-qPCR (the primers designed for RT-qPCR analyses are shown in Additional file 5: Figure S2). The samples, cultured in the liquid AM for 2, 3.5, and 5days, were respectively collected, and the mycelia tissues were separated from MS using forceps. Subsequently, the mycelia were washed with sterile distilled water for at least two times. The total RNA was extracted using the TRIzol® reagent (Invitrogen, USA) and purified with RNA purification kit (Qiagen, USA). The first-strand cDNA fragment was synthesized from purified RNA. A total of 17 genes were chosen for RT-qPCR. All the genes were first subjected to RT-PCR analysis, and the gel was run to obtain suitable conditions for RT-qPCR analysis. The RT-qPCR reaction was performed according to the protocol of Green Mastermix (Qiagen, USA). Three biological duplicates of each sample and triplicates of each reaction were acquired. The genes tef and tub were used for calibration in all the experiments. The expression ratios were calculated from cycle threshold values using the 2−△△CT method.
This study was financially supported by Special Fund for Agro-scientific in the Public Interest (Project 201103002) and State 863 Project of China (Project 2011AA10A201).
- Boucias DG, Bradford DL, Barfield CS: Susceptibility of the velvetbean caterpillar and soybean looper (Lepidoptera: Noctuidae) toNomuraea rileyi: effects of pathotype, dosage, temperature and host age.J. Econ. Entomol. 1984, 77: 247-253.View ArticleGoogle Scholar
- Chen H, Yin Y, Li Y, Mahmud MS, Wang ZK: Identification and analysis of gene differentially expressed in the Spodoptera litura fat body in response to the biocontrol fungus, Nomuraea rileyi. Comp Biochem Physiol B Biochem Mol Biol. 2012, 163 (2): 203-210. 10.1016/j.cbpb.2012.05.018.View ArticlePubMedGoogle Scholar
- De Faria MR, Wraight SP: Mycoinsecticides and Mycoacaricides: a comprehensive list with worldwide coverage and international classification of formulation types. Biol Control. 2007, 43: 237-256. 10.1016/j.biocontrol.2007.08.001.View ArticleGoogle Scholar
- Boland GJ, Hall R: Index of plant hosts of Sclerotinia sclerotiorum. Canadian J Plant Pathol. 1995, 16: 93-108.View ArticleGoogle Scholar
- Huang HC, Miao HJ, Erickson RS, Jiang DH, Xiao YN, LiGQ: Biological control of sclerotinia diseases of rapeseed by aerial applications of the mycoparasite Coniothyrium minitans. Eur J Plant Pathol. 2006, 114: 345-355. 10.1007/s10658-005-2232-6.View ArticleGoogle Scholar
- Jackson MA, Schisler DA: Liquid culture production of microsclerotia of Colletotrichum truncatum for use as bioherbicidal propagules. Mycol Res. 1995, 99 (7): 879-884. 10.1016/S0953-7562(09)80745-4.View ArticleGoogle Scholar
- Shearer JF, Jackson MA: Liquid culturing of microsclerotia of Mycoleptodiscus terrestris, a potential biological control agent for the management of Hydrilla. Biol Control. 2006, 38: 298-306. 10.1016/j.biocontrol.2006.04.012.View ArticleGoogle Scholar
- Jackson MA, Jaronski ST: Production of microsclerotia of the fungal entomopathogen Metarhizium anisopliae and their potential for use as a biocontrol agent for soil-inhabiting insects. Mycol Res. 2009, 113: 842-850. 10.1016/j.mycres.2009.03.004.View ArticlePubMedGoogle Scholar
- Griffiths DA: The fine structure of developing microsclerotia of Verticillium dahliae Kleb. Arch Microbiol. 1970, 74: 207-212.Google Scholar
- Perry JW, Evert RF: Structure of microsclerotia of Verticillium dahliaein roots of ‘Russett Burbank’ potatoes. Can J Bot. 1982, 62: 396-401.View ArticleGoogle Scholar
- Jones D: Ultrastructure and composition of the cell walls of Sclerotinia sclerotiorum. Trans Br Mycol Soc. 1970, 54: 351-360. 10.1016/S0007-1536(70)80148-6.View ArticleGoogle Scholar
- Colotelo N: A scanning electron microscope study of developing sclerotinia of Sclerotinia sclerotiorum. Can J Bot. 1974, 52: 1127-1130. 10.1139/b74-143.View ArticleGoogle Scholar
- López-Escudero FJ, Mwanza CMA, Blanco-López MA: Production of homogeneous and viable Verticillium dahliae microsclerotia effective for Verticillium wilt studies. Biotechnology. 2006, 5 (4): 421-428.View ArticleGoogle Scholar
- Neumann MJ, Dobinson KF: Sequence tag analysis of gene expression in the vascular wilt pathogen Verticillium dahliae during microsclerotia development and pathogenic growth. Fungal GenetBiol. 2003, 38: 54-62. 10.1016/S1087-1845(02)00507-8.View ArticleGoogle Scholar
- Klimes A, Dobinson KF: A hydrophobin gene, VDH1, is involved in microsclerotial development and spore viability in the plant pathogen Verticillium dahliae. Fungal Genet Biol. 2006, 43: 283-294. 10.1016/j.fgb.2005.12.006.View ArticlePubMedGoogle Scholar
- Klimes A, Amyotte SG, Grant S, Kang S, Dobinson KF: Microsclerotia development inVerticillium dahliae: Regulation and differential expression of the hydrophobin geneVDH1.Fungal Genet Biol. 2008, 45: 1525-1532. 10.1016/j.fgb.2008.09.014.View ArticlePubMedGoogle Scholar
- Gao F, Zhou BJ, Li GY, Jia PS, Li H, Zhao YL, Zhao P, Xia GX, Guo HS: A glutamic acid-rich protein identified in Verticillium dahliae from an insertional mutagenesis affects microsclerotial formation and pathogenicity. PLoS One. 2010, 5 (12): e15319-10.1371/journal.pone.0015319.PubMed CentralView ArticlePubMedGoogle Scholar
- Rauyaree P, Ospina-Giraldo MD, Kang S, Bhat RG, Subbarao KV, Grant SJ, Dobinson KF: Mutations in VMK1, a mitogen-activated protein kinase gene, affect microsclerotia formation and pathogenicity in Verticillium dahliae. Curr Genet. 2005, 48: 109-116. 10.1007/s00294-005-0586-0.View ArticlePubMedGoogle Scholar
- Duran RM, Cary JW, Calvo AM: Production of cyclopiazonic acid, aflatrem, and aflatoxin by Aspergillus flavusis regulated by veA, a gene necessary for sclerotial formation. Appl Microbiol Biotechnol. 2007, 73: 1158-1168.View ArticlePubMedGoogle Scholar
- Liu YY, Guo SX: Involvement of Ca2+ channel signaling in sclerotial formation of Polyporus umbellatus. Mycopathologia. 2010, 169: 139-150. 10.1007/s11046-009-9238-0.View ArticlePubMedGoogle Scholar
- Jin FJ, Takahashi T, Matsushima KI, Hara S, Shinohara Y, Maruyama JI, Kitamoto K, Koyama Y: SclR, a basic helix-loop-helix transcription factor, regulates hyphal morphology and promotes sclerotial formation in Aspergillus oryzae. Eukaryotic Cell. 2011, 10 (7): 945-955. 10.1128/EC.00013-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster SC: Next-generation sequencing transforms today’s biology. Nat Methods. 2008, 5: 16-18.View ArticlePubMedGoogle Scholar
- Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS: De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010, 11: 400-10.1186/1471-2164-11-400.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Gōtz S, García-Gómez JM, Terol J, Talón M, Robles M, Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-297. 10.1093/nar/gkl031.PubMed CentralView ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.PubMedGoogle Scholar
- Gerogiou CD, Patsoukis N, Papapostolou I, Zervoudakis G: Sclerotial metamorphosis in filamentous fungi is induced by oxidative stress. Integr Comp Biol. 2006, 46 (6): 691-712. 10.1093/icb/icj034.View ArticleGoogle Scholar
- Patsoukis N, Georgiou CD: Thiol redox state and oxidative stress affect sclerotial differentiation of the phytopathogenic fungi Sclerotium rolfsii and Sclerotinia sclerotiorum. J Appl Microbiol. 2008, 104: 42-50.PubMedGoogle Scholar
- Oh Y, Donofrio N, Pan H, Coughlan S, Brown DE, Meng SW, Mitchell T, Dean RA: Transcriptome analysis reveals new insight into appressorium formation and function in the rice blast fungus Magnaporthe oryzae. Genome Biol. 2008, 9: R85-10.1186/gb-2008-9-5-r85.PubMed CentralView ArticlePubMedGoogle Scholar
- Egan MJ, Wang ZY, Jones MA, Smirnoff N, Talbot NJ: Generation of reactive oxygen species by fungal NADPH oxidases is required for rice blast disease. PNAS. 2007, 104 (28): 11772-11777. 10.1073/pnas.0700574104.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahashi M: Oxidative stress and redox regulation on In Vitro development of mammalian embryos. J Reprod Dev. 2012, 58 (1): 1-9. 10.1262/jrd.11-138N.View ArticlePubMedGoogle Scholar
- Yang Y, Sharma R, Sharma A, Awasthi S, Awasthi YC: Lipid peroxidation and cell cycle signaling: 4-hydroxynonenal, a key molecule in stress mediated signaling. Acta Biochemica Polonica. 2003, 50 (2): 319-336.Google Scholar
- Srinivasa K, Kim J, Yee S, Kim W, Choi W: A MAPK kinase pathway is implicated in the pseudohyphal induction by hydrogen peroxide in Candidaal bicans. Molecules and Cells. 2012, 33 (29): 183-193.PubMed CentralView ArticlePubMedGoogle Scholar
- Aguirre J, Momberg MR, Hewitt D, Hansberg W: Reactive oxygen species and development in microbial eukaryotes. Trends Microbiol. 2005, 13 (3): 111-118. 10.1016/j.tim.2005.01.007.View ArticlePubMedGoogle Scholar
- Xie XQ, Li F, Ying SH, Feng MG: Additive contributions of two manganese-cored superoxide dismutases (MnSODs) to antioxidation, UV tolerance and virulence of Beauveria bassiana. PLoS One. 2012, 7 (1): e30298-10.1371/journal.pone.0030298.PubMed CentralView ArticlePubMedGoogle Scholar
- Reth M: Hydrogen peroxide as second messenger in lymphocyte activation. Nat Immunol. 2002, 3 (12): 1129-1134. 10.1038/ni1202-1129.View ArticlePubMedGoogle Scholar
- Liang Y, Rahman MH, Strelkov SE, Kav NNV: Developmentally induced changes in the sclerotial proteome of Sclerotinia sclerotiorum. Fungal Biology. 2010, 114: 619-627. 10.1016/j.funbio.2010.05.003.View ArticlePubMedGoogle Scholar
- Papapostolou I, Georgiou CD: Superoxide radical is involved in the sclerotial differentiation of filamentous phytopathogenic fungi: identification of a fungal xanthine oxidase. Fungal Biology. 2010, 114: 387-395. 10.1016/j.funbio.2010.01.010.View ArticlePubMedGoogle Scholar
- Papapostolou I, Georgiou CD: Hydrogen peroxide is involved in the sclerotial differentiation of filamentous phytopathogenic fungi. Appl Microbiol. 2010, 109: 1926-1936.View ArticleGoogle Scholar
- Patsoukis N, Georgiou CD: Thiol redox state and related enzymes in sclerotium-forming filamentous phytopathogenic fungi. Mycol Res. 2008, 112: 602-610. 10.1016/j.mycres.2007.10.006.View ArticlePubMedGoogle Scholar
- Singh K, Nizam S, Sinha M, Verma PK: Comparative transcriptome analysis of the necrotrophic fungus Ascochyta rabiei during oxidative stress: insight for fungal survival in the host plant. PLoS One. 2012, 7 (3): e33128-10.1371/journal.pone.0033128.PubMed CentralView ArticlePubMedGoogle Scholar
- Mathioni SM, Beló A, Rizzo CJ, Dean RA, Donofrio NM: Transcriptome profiling of the rice blast fungus during invasive plant infection and in vitro stresses. BMC Genomics. 2011, 12: 49-10.1186/1471-2164-12-49.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartosz G: Reactive oxygen species: destroyers or messengers?. Biochem Pharmacol. 2009, 77: 1303-1315. 10.1016/j.bcp.2008.11.009.View ArticlePubMedGoogle Scholar
- Kim D, Kim MS, Cho KH: The core regulation module of stress-responsive regulatory networks in yeast. Nucleic Acids Res. 2012, 40 (18): 8793-8802. 10.1093/nar/gks649.PubMed CentralView ArticlePubMedGoogle Scholar
- Lan LF, Cheng A, Dunman PM, Missiakas D, He C: Golden pigment production and virulence gene expression are affected by metabolisms in Staphylococcus aureus. J Bacteriol. 2010, 192 (12): 3068-3077. 10.1128/JB.00928-09.PubMed CentralView ArticlePubMedGoogle Scholar
- ABO Ellil AHA: Sclerotial development, melanin production and lipid peroxidation by Sclerotium rolfsii. Folia Microbiol. 1999, 44 (2): 181-186. 10.1007/BF02816239.View ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hachen N, Geirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.PubMed CentralView ArticlePubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19 (5): 651-652. 10.1093/bioinformatics/btg034.View ArticlePubMedGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-484.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.