- Research article
- Open Access
Dual species transcript profiling during the interaction between banana (Musa acuminata) and the fungal pathogen Fusarium oxysporum f. sp. cubense
© The Author(s). 2019
- Received: 21 January 2019
- Accepted: 12 June 2019
- Published: 24 June 2019
Banana wilt disease, caused by Fusarium oxysporum f. sp. cubense Tropical Race 4 (Foc TR4), is one of the most devastating diseases in banana (Musa spp.). Foc is a soil borne pathogen that causes rot of the roots or wilt of leaves by colonizing the xylem vessels. The dual RNA sequencing is used to simultaneously assess the transcriptomes of pathogen and host. This method greatly helps to understand the responses of pathogen and host to each other and discover the potential pathogenic mechanism.
Plantlets of two economically important banana cultivars, Foc TR4 less susceptible cultivar NK and susceptible cultivar BX, were used to research the Foc-banana interaction mechanism. Notably, the infected NK had more significantly up-regulated genes on the respiration machinery including TCA cycle, glyoxylate, glycerol, and glycolysis compared to BX at 27 h post inoculation (hpi). In addition, genes involved in plant-pathogen interaction, starch, sucrose, linolenic acid and sphingolipid metabolisms were uniquely more greatly induced in BX than those in NK during the whole infection. Genes related to the biosynthesis and metabolism of SA and JA were greatly induced in the infected NK; while auxin and abscisic acid metabolisms related genes were strongly stimulated in the infected BX at 27 hpi. Furthermore, most of fungal genes were more highly expressed in the roots of BX than in those of NK. The fungal genes related to pathogenicity, pectin and chitin metabolism, reactive oxygen scavenging played the important roles during the infection of Foc. CCP1 (cytochrome c peroxidase 1) was verified to involve in cellulose utilization, oxidative stress response and pathogenicity of fungus.
The transcriptome indicated that NK had much faster defense response against Foc TR4 than BX and the expression levels of fungal genes were higher in BX than those in NK. The metabolisms of carbon, nitrogen, and signal transduction molecular were differentially involved in pathogen infection in BX and NK. Additionally, the putative virulence associated fungal genes involved in colonization, nutrition acquirement and transport provided more insights into the infection process of Foc TR4 in banana roots.
- Musa spp.
- Fusarium oxysporum
- Fusarium wilt
- Cytochrome c peroxidase
Banana (Musa spp.) ranks fourth globally in terms of gross value production, after rice, wheat and maize . Fusarium wilt in banana, also called Panama disease, is caused by the fungus Fusarium oxysporum f. sp. cubense (Foc). The most virulent race of Foc for ‘Cavendish’ banana is Foc tropical race 4 (Foc TR4). This fungus severely limits the productivity of the leading Cavendish banana cultivar, resulting in drastic economic losses throughout the banana-producing regions . The vascular wilt fungus Fusarium oxysporum is an asexual, soil inhabiting facultative parasite . Fusarium oxysporum has different asexual spores: microconidia directly penetrate the living roots, macroconidia usually colonize the surface of dead plants, and chlamydospores lie dormant in soil for decades [4, 5]. Foc causes severe diseases in many economically important crops, including tomato (Solanum lycopersicum L.) , cucumber (Cucumis sativus L.) , and cotton (Gossypium hirsutum L.) . A seedling plant infected by Foc shows drastic physiological and morphological alterations, including slowed growth, leaf chlorosis, and finally, whole plant wilt. The fungal colonization disrupts the transportation of nutrition and water of plants through vascular tissue. In response to this biotic stress, plants induce a range of immune responses, including physical barriers (e.g., cell wall-associated compounds)  and molecular responses (e.g., hypersensitive response, production of reactive oxygen species, activation of signal molecules, and expression of pathogen-related genes) .
High-throughput sequencing technology, especially the RNA-seq, facilitates the more precise molecular changes of plants that accompany biotic and abiotic stress. In previous researches, many banana plant genes related to defense-related pathways mainly including phenylpropanoid biosynthesis, phenylalanine metabolism, carbon metabolisms, amino acid recycling, hormone signal transduction and plant-pathogen interaction were greatly induced in the face of Foc TR4 [11–14]. Meanwhile, many fungal genes related to virulence, transporters and transcription factors for toxin and nutrient, signaling pathways, and adaptation to host were strongly activated in the co-culture with banana roots . Recently, the dual RNA-seq has become a powerful method to comprehensively understand the interaction between host and pathogen in vivo . This dual sequencing uncovers biosynthetic and metabolic pathways of cross-talk from participants simultaneously and specifically links the dynamic expression profiles of genes with the interaction [17, 18]. However, the mutual response of attack and counterattack in vivo between banana and fungus has been poorly known till now.
In the present work, the susceptible cultivar ‘Baxi’ (BX) and its less susceptible mutant cultivar ‘Nongke No. 1’ (NK) to Foc TR4 were employed for Foc infection [12, 19]. NK is bred from BX clones and has the similar growth period and morphology of fruits with BX; however, it has the higher soluble sugar and less susceptible to Foc TR4 than BX . To research the infection response at the early stage, the banana sampling times of 27 and 51 h post inoculation (hpi) were determined according to our previous research .
This study begins with a comprehensive transcriptomic analysis of the infected banana roots and Foc TR4 in vivo using the dual RNA-seq method. The differentially expressed banana genes from the infected BX and NK indicated the potential response mechanisms of different banana cultivars to Foc. Moreover, the expression levels of fungal genes in different hosts at two time points were investigated at an unprecedented depth. The faster response of banana defense genes and the lower fold change of fungal genes in NK compared to BX might determine the less susceptibility of NK. Furthermore, fungal gene encoding cytochrome c peroxidase (CCP1) was investigated to play a role in cellulose utilization. The low disease incidence in BX plantlets infected by the CCP1 mutant strains indicated that CCP1 might play an important role in fungal pathogenicity.
Overall, this transcriptome dataset allows us to gain insights into the complex banana-Foc TR4 interaction and to investigate the potential pathogenesis of Foc. Moreover, it would benefit to attenuate/control fungal pathogenicity or breed the resistant cultivars.
Symptoms of the infected banana plants
Banana root transcriptomes during the early stages of Foc TR4 infection
Number of differentially expressed genes (DEGs) in susceptible and less susceptible banana cultivars during Foc infection
Banana roots showed up-regulation of transcripts related to energy metabolism during infection
The GO term analysis also uncovered the differences between the infected BX and NK under q value < 0.05. For instance, the DEGs with organic substance metabolic process and molecular functions were uniquely induced in NK, as eleven GO terms were not distributed by NK at either time point. In addition, the down-regulated DEGs from two cultivars distributed in different terms except for ‘heme and tetrapyrole binding’ (Additional file 4: Figure S2).
We also investigated the function enrichment of the DEGs between mock-inoculated BX and NK cultivars (DEMG) through KEGG and GO. The up-regulated DEMGs were mainly related to citrate cycle and phenylalamine metabolism, cutin, suberine, wax, steroid biosynthesis through KEGG analysis (Additional file 5: Figure S3). The down-regulated DEMGs were mainly clustered in cysteine, methionmine, glycerolipid metabolism and fatty acid elongation. Under GO analysis, the up-regulated DEMGs were only clustered at 27 hpi, and they greatly distributed in biological, metabolic processes, and catalytic activity terms. The enriched up-regulated DEMGs at 51 hpi uniquely distributed into the terms of molecular function except for the cofactor binding (Additional file 6: Figure S4). DEGs from the infected banana roots, related to ‘plant-pathogen interaction’, ‘hormone signal transduction’ or ‘phenylpropanoid biosynthesis’, were not found in these DEMGs.
Insight into signal molecules metabolism in infected banana plants
Salicylic acid and Jasmonic acid
More DEGs involved in SA biosynthesis and metabolism were greatly up-regulated in NK compared to those in BX at 27 hpi, including phenylalanine ammonia-lyase (PAL) and pathogenesis-related protein 1 (PR-1) genes (Fig. 4). By contrast, more DEGs were greatly up-regulated in BX than those in NK at 51 hpi. As for JA, all jasmonate ZIM domain-containing protein (JAZ) genes, a negative regulator in JA transduction, had the higher expression changes in NK than those in BX at 27 hpi. The 12-oxophytodienoate reductase (OTR) and allene oxide synthase (AOS) genes, responsible for the biosynthesis of JA, had higher expression fold changes at 51 hpi compared to those at 27 hpi in both cultivars. These genes expression profiles were consistent with our previous data: the infected NK had a higher concentration of SA and a lower concentration of JA than the infected BX at 27 hpi .
Primary auxin-responsive DEGs fall into three major classes: auxin responsive GH3 gene family (GH3), auxin/indole-3-acetic acid protein (AUX/IAA), and small auxin-up RNA (SAUR). Most of three families genes were greatly induced in the infected BX at 27 hpi, and half of them were depressed at 51 hpi. In our previous research, the level of 3-indoleacetic acid was lower in BX roots at 51 hpi than that at 27 hpi . The depressed AUX/IAA genes might indicate the lower level of 3-indoleacetic acid in NK at 27 hpi than that in BX.
The DEGs encoding abscisic acid receptor (PYR/PYL) and protein phosphatase 2C (PP2C) were more induced in BX compared to in NK at 27 hpi. At 51 hpi, more than half of these genes were depressed in NK, suggesting that the metabolism of ABA might be inhibited in NK (Fig. 4). The higher fold changes of these DEGs at 27 hpi than those at 51 hpi also verified the peak level of ABA at 27 hpi in BX .
Some DEGs (mean FPKM > 10) related to cytokinin (A-ARR family), gibberellin (DELLA protein), brassinosteroid (shaggy-related protein kinase eta-like isoform X1, BIN2), and ethylene (ethylene-responsive transcription factor 1, ERF1/2) biosynthesis and metabolism were also found in both banana cultivars (Additional file 8: Table S4).
Identification and functional analysis of fungal transcripts present in banana roots
The RNA-seq data was also mapped to the Foc TR4 genome with 18,065 potential protein coding genes. About 8% of the reads from the infected samples (14,459 transcripts) were mapped to the Foc TR4 genome (Additional file 9: Table S5). When the data were tallied under a less stringent threshold (> 0.5 FPKM), 10,777 fungal transcripts accounting for about 60% of the predicted fungal genes were obtained in this research. It was notable that there were more fungal genome mapped-reads present in BX roots than in NK (6.29% vs. 3.23% at 27 hpi and 2.23% vs. 1.59% at 51 hpi). Additionally, more up-regulated differentially expressed fungal genes (DEFGs) were from BX roots compared to those from NK roots at both time points (Additional file 10: Table S6).
Characterizing the fungal genes expressed within plants
During its colonization in banana, fungus expressed a number of genes involved in detoxification and stress tolerance . For instance, catalase peroxidase (CAT), superoxide dismutase (SOD) and cytochrome c peroxidase (CCP) genes were distinctively expressed in both hosts. Furthermore, two CATs were significantly stimulated from BX, indicating that the fungus mounted more defenses in BX than in NK during the early infection stage. In addition, seventeen glutathione S-transferase (GST) family members, which combine glutathione with an electrophilic substance, and seven cytochrome P450 (CTP) members, which are indispensable for pathogens adapting to a hostile environment , were highly expressed in both hosts (Fig. 6). The fungus also induced many transporter genes to facilitate the infection of the plant. For instance, the ATP binding cassette superfamily (ABC) transporters, which mainly function in pathogen defense and transport of virulence factors , were also induced in both cultivars. The transcriptome data of all fungal genes were shown in Additional file 11: Table S7.
Evidence of the nutritional-acquisition strategy of Foc TR4 fungus
Expression levels of fungal genes related to nutrition utilization from banana roots
Nitrogen assimilation transcription factor
Potential function of fungal genes CCP 1
One of the earliest defense reactions activated in plant tissues in response to pathogen attack is the accumulation of reactive oxygen species (ROS) [23, 28]. Necrotrophic pathogens even stimulate H2O2 production [29, 30]. Catalase (CAT) from septoria tritici in wheat effectively eliminates H2O2 to maintain the intercellular redox balance . We also verified that CAT was greatly related to the pathogenicity of Foc TR4 during its infection into banana roots . Cytochrome c peroxidase (CCP) is an antioxidant enzyme in the mitochondria; it is a possible sensor in the oxidative stress response in fungi ; and it is also one of the key enzymes for camalexin biosynthesis in plants . To seek the role of CCP named as CCP1 in Foc TR4, wild type (WT), mutant (ccp1) and complemented fungal strains (Comp) were generated to study their colonization, development, and pathogenicity in banana plants.
The ccp1 and the Comp strains were constructed through PCR overlap and homologous recombination. The ccp1 strains carried a cassette of green fluorescent protein and hygromycin, and the Comp strains carried a neomycin cassette for their selection and verification (Additional file 12: Figure S5). The ccp1 strains showed strong green fluorescence and the multiple PCRs indicated that CCP1 was successfully knocked out from the strains (Additional file 13: Figure S6A, B, C and D). The Comp strains were verified by green fluorescence deletion and the PCR products of about 820 bp. The PCR products were amplified with the CCP1-specific primer pair F/R and the PCR product has the same sequence with that of CCP1 (Additional file 13: Fig. S6).
Plants and pathogens constantly compete for survival, and the physiological and molecular change results in dynamic but regulated alterations in gene expression occurring in both organisms. In this respect, we present a dual high throughput RNA-seq on banana roots and Foc TR4 to unveil the dynamic changes in both species. It will provide a broader understanding of the banana defense response triggered by Foc TR4, or identify a potential pathogenic mechanism of Foc TR4 in both hosts with similar genetic background.
Since the transcriptomic profile reflects the culture condition and the treatment method, the test tissue and the fungi are crucial to this dual RNA-seq approach. In our research, the wounded banana roots wrapped with fungal agar block had more hyphea or conidia than those dipped into fungal suspension . Our infection method resulted in some different expression profiles of genes from the previous researches [11, 12, 14]. For example, the number of up-regulated DEGs was more than that of the down-regulated DEGs in BX at both time points, while it showed the contrary in BX with roots soaking at 3 days after infection . In addition, more up-regulated DEGs were found in NK at 27 hpi than in BX, although they were contrary at 51 hpi. Based on our dual RNA transcriptomic analysis, we try to present a model that shows the major molecular metabolism occurring in the test host and fungi (Additional file 15: Fig. S8 and data mainly from Additional files 3 and 16: Table S2 and S8). Briefly, the DEGs associated with glycolysis and glyoxylate cycle were widely involved in the defense response in the infected banana. Furthermore, the higher fold changes of these genes in the two pathways were found in NK than in BX at 27 hpi. Most DEGs required for the tricarboxylic acid cycle (TCA) were up-regulated only in the infected NK at 27 hpi, indicating that energy consumption was accelerated in the less susceptible cultivar. Together, many genes encoding trehalose phosphate phosphatase were more highly up-regulated in NK than in BX at 27 hpi, and a shift in metabolism could provide extra energy for banana plants to maintain growth and metabolism. These increasing in metabolism in NK at 27 hpi indicated that more energy would be available to resist the biotic stress at the infection early stage. However, it was likely that the glycerol metabolism were inhibited in NK due to more down-regulated genes at 27 hpi than in BX. Curiously, most of DEGs committed into the fatty acid degradation were only found in BX at 27 hpi. Based on our transcriptome, DEGs with the degradation of carbohydrates, including sucrose synthase, trehalose phosphate phosphatase, and α-amylase, were induced in both banana cultivars.
In addition, SA and JA are important signal molecules during the plant response to biotic and abiotic stresses [39, 40]. In our study, SA and JA were more greatly induced in the infected NK with the higher gene fold change than in the infected BX at 27 hpi. It suggested that SA and JA might be more early induced in the infected NK than in the infected BX. On the contrary, the fold change of DEGs involved in AUX and ABA metabolism showed higher in BX than in NK at 27 hpi. Phytohormone had cross-talk during modulating disease respond in plant . The emerging evidence verifies that auxin is widely involved in plant disease susceptibility . GH3 increases the combination of IAA with amino acids to inactivate IAA, which stimulates the expression of PR-1 and CBF in response to stress . The levels of ABA are usually negatively associated with plant resistance due to its antagonistic interaction with the SA signaling pathway [44–46]. The induced DEGs related to ABA metabolism might contribute to the lower concentration of SA in BX than in NK.
The defense-related genes encoding PR1, thaumatin-like proteins (PR-5) peroxidase (PR-9), and chitinase play important roles in pathogen-infected plants or challenge with virulence factors from pathogens [9, 47]. These genes were more highly induced in NK than in BX at 27 hpi, which might be intimately related to the resistance of NK at the early infection stage (Additional file 3).
It was obvious that the chitin and pectin metabolism in Foc TR4 from the infected BX were higher compared to those from the infected NK. Chitin is an essential structural component of fungal cell walls, and fungus has multiple classes of CHS (chitin synthase) that deacetylate chitin to chitosan. The roles of some fungal CHS genes have been investigated in Candida albicans, Aspergillus fumigatus, Aspergillus nidulans and Wangiella dermatiditis [48, 49]. In our manuscript, more CHS genes were found compared to the previous research , and they were distributed in six classes according to the phylogenetic analysis [50, 51] (Additional file 17: Figure S9). It was notable that FOC4_g10006856 was clustered with FoCHSV that was required for pathogenesis of Fusarium oxysporum [52, 53], and it had the highest expression level (mean FPKM > 600 in BX and > 450 in NK) among all CHS genes. PL (pectin lyase) and PEL (pectate lyase) are the first two enzymes secreted by fungal pathogens when they directly degrade plant cell walls [54–56]. The high expression levels of PL and PEL genes might mean higher pathogenicity in the infected BX than in the infected NK.
Meanwhile, many fungal genes benefitting Foc TR4 infection and colonization were greatly induced in the infected banana roots. For instance, fungal effectors, defined as pathogen proteins and small molecules, facilitate infection and trigger defense response in fungi . The Secreted In Xylem (SIX) gene has been identified as one family of effectors in Fusarium oxysporum f sp. cubense [58, 59]. The increasing expression level of SIX from 27 to 51 hpi (FPKM = 31 to 262.3) in the infected BX suggested that fungus might be more virulent to the susceptible BX, since it was hardly detected in the infected NK at 51 hpi. In agreement with this finding, 60 and 55% of the identified fungal genes were > 1-fold higher in BX than those in NK at 27 hpi and 51 hpi, respectively (data not shown). The higher expression levels of fungal genes in the BX might more suppress or postpone defense response in the susceptible BX than in the less susceptible NK. Many potential effectors including endoglucanase, elongation factor, and endopolygalacturonase found in our research needed to be investigated in the future (Additional file 18: Table S9). The fungal catalase, superoxide dismutase, cytochrome c peroxidase, glutathione S-transferase, and cytochrome P450 were widely induced in both hosts, which greatly contribute to remove active oxygen species from host, reduce the harm of reactive oxygen, and contribute to adaptation of fungi in host [60, 61].
The cytochrome c peroxidase (CCP) is a water-soluble heme-containing enzyme of the peroxidase family. It catalyzes the oxidation of organic substrates by H2O2 to resist peroxide-induce cellular stress in eukaryotes and prokaryotes [62–64]. In addition to its peroxidase activity, it acts as a sensor and signaling molecule to exogenous H2O2, which activates mitochondrial catalase activity . Campylobacter jejuni carrying a mutation in the Cj0020c (CCP) gene are unable to effectively colonise chicks . In our report, the mutant (ccp1) showed the greatly reduced virulence in banana and this is a unique result in our study, as Cryotococcus neoformans or the plant pathogen Septoria tritici lacking the CCP gene did not show diminished virulence [31, 67].
However, the fungal transcriptome also indicated that the lower mapped reads at 51 hpi than those at 27 hpi (Additional file 9: Table S5), which was inconsistent with the amount of hyphae or spores in the roots. It was reported that approximately 86.06 and 81.08% of genes in Foc TR4 were transcribed at vegetative stage and 48 h post inoculation , suggesting that the transcription of genes in Foc TR4 should be inhibited when facing to host. We speculated that the lower mapped reads at 51 hpi than at 27 hpi from Foc might be due to the inhibition from banana increasing immune response, which was consistent with the increased mapped reads from banana. However, the phenomenon might be temporary at the early infection stage, and it need to be investigated further.
Taken together, our research reported the intriguing interaction between Fusarium oxysporum f. sp. cubense and banana plant, which make wilt disease for further in-depth studies in future.
1. Through the dual RNA-seq method, the dynamic transcript profiles at the early stages occurred in banana roots responding to Foc TR4 were explored in unprecedented perspectives. About 60 and 8% of the reads were mapped to the banana and fungal genome; however, 97% (35,276) of banana genes and 60% (6690) of the fungal genes were obtained in our research.
2. Infection with fungus triggers massive metabolic reprogramming in the banana roots using KEGG and GO analysis. More DEGs related to glycolysis, tricarboxylic acid cycle, and degradation of carbohydrates in NK at 27 hpi than in BX might be related to the less susceptibility of NK to Foc TR4.
3. The DEGs involved in the biosynthesis and metabolism of salicylic acid, jasmonic acid, auxin, and abscisic acid showed different expression changes between both cultivars, which contribute to the immune response of banana to Foc TR4.
4. The genes encoding the potential antivirulence, virulence, and effectors proteins were identified from Foc fungus. The high expression levels of these fungal genes in BX discovered that the fungus was more virulent to BX than NK, which might result in the susceptibility of BX to Foc TR4.
Plant growth conditions and fungus inoculation
The susceptible BX and the less susceptible NK of banana (Musa acuminata) cultivars were from the Tropical Crop Germplasm Research Institute, Chinese Academy of Tropical Agricultural Sciences (Danzhou city, Hainan, China). Tissue culture-derived plantlets were grown in pots under 16-h light/8-h dark at 28 °C with 100 μmol m− 2 s− 1 light intensity for 90 days. The Foc TR4 strains (VCG01213/16) isolated from Hainan in China by Dr. Junsheng Huang (Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou, China) were used in our research. The fungus was transformed with GFP according to the reported method  and then used for the infection processing detection in banana roots. The GFP-expressing Foc TR4 was newly transferred to the potato dextrose agar (PDA) media with 100 μg/ml hygromycin B 4 to 5 days before plant inoculation.
The root epidermis was artificially damaged (about 0.5 cm2) and then covered with a freshly prepared Foc agar block of about 0.5 cm2 for fungus inoculation . The treated banana roots above 0.5 cm of the infection point were observed with a laser confocal microscope (OLYMPUS, FV10-ASW) using GFP filters .
The number of hyphea or conidia in the vascular tissue was counted through the green fluorescent light points in an intact root transverse section. Six banana plants with 4 infected roots for every treatment samples and a total of about 60–80 roots sections were obtained from three replications. The statistical significance between different samples was carried out using Student’s test.
Plant sampling and sample sequencing
About 10 roots of three banana plants inoculated by wild Foc TR4 for every treatment were collected at 27 h and 51 h post inoculation (hpi), respectively. RNA (3 μg) was extracted from the treated parts of banana roots using an RNA extraction kit (OMEGA, USA). Poly (A) -RNA enrichment and a strand-specific RNA-seq library were prepared using the NEBNext Ultra™ Directional RNA Library Prep Kit for Illumina (NEB, USA) following manufacturer’s recommendations. Library quality was assessed on the Agilent Bioanalyzer 2100 system (Agilent Technologies, USA). Libraries were sequenced on an Illumina Hiseq 2500 platform with 125-bp paired-end reads.
Read mapping and gene expression quantification
After removing low-quality sequences containing uncalled bases (Ns), Bowtie v2.0.6  and TopHat v2.0.9  were used to align the RNA-seq reads against 36,542 gene models from M. acuminata (version 2, https://banana-genome-hub.southgreen.fr/) and 18,037 genes models from Foc (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nuccore/AMGQ00000000). Cuffdiff (v2.2.1) was used to calculate FPKMs (fragments per kilobase of transcript per million mapped reads) and to determine differential expression in each sample . Hierarchical clustering for banana and fungal genes was performed with MultiExperiment Viewer (MeV) (version 4.8.1) from the TM4 suite (http://www.tm4.org) . Genes were clustered using the Euclidean Distance measure and the complete linkage clustering method. PCA was conducted using the prcomp command with default parameters in the R software package. The transcriptome datasets are available at the NCBI Sequence Read Archive (SRA) with the accession number from SRR8661622 to 45. The expression levels of all banana plant transcripts are shown in Additional file 19: Table S10.
Gene functional classification
The functions of genes were annotated using BLASTx with an e-value threshold of 10− 5 against the NCBI non-redundant (NR) protein database, the KEGG (Kyoto Encyclopedia of Genes and Genomes) protein database, COGs (NCBI phylogenetic classification of proteins encoded in complete genomes), and the Swiss-Prot database.
GO and KEGG enrichment analyses on differentially expressed genes
Gene Ontology (GO) enrichment analysis of differentially expressed genes was conducted using the GOseq R package. GO terms with corrected p-values less than 0.05 were considered significantly enriched. KOBAS software was used to test the statistical enrichment of DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.
Comparing RNA-seq and qPCR results
The qPCR assays were performed to confirm RNA-seq results by an independent technique (Additional file 2: Table S1). The expression levels of 30 banana and 10 Foc genes were analyzed by qPCR from the RNA sampled above. The qPCR experiments were conducted on a Step One Real-Time PCR system (Applied Biosystems) using SYBR Green I (Takara, Japan). Each reaction was performed in a final volume of 20 μl, containing 10 μl of 2 × SYBR Green PCR Master Mix (Takara, Japan), 200 nM each gene-specific primer, and 50 ng cDNA template. No-template reactions were included as negative controls for each set of primers used. The thermal cycling conditions were 95 °C for 30 s, followed by 40 cycles of 5 s at 95 °C, 20 s at 58–62 °C depending on primer melting temperature, and 20 s at 72 °C, with fluorescence detection at the end of each cycle. The amplification of a single product per reaction was confirmed by melting curve analysis. All reactions were performed in technical triplicates. Banana data were normalized using two reference genes that showed little variation in the RNA-seq analysis (Actin and glyceraldehydes-3-phosphate dehydrogenase 2, GAPDH). Expression levels of fungal genes were given in relation to the fungal β-actin and IF3b (transcription initiation factor) genes. Primers used in these experiments were designed by Primer Premier 6.0 software.
Construction of mutant and complementation cassettes targeting cytochrome c peroxidase (CCP1)
The mutant of CCP1 deletion (ccp1) was performed by homologous recombination. In brief, the 5′- and 3′-flanking sequences of the CCP1 gene were amplified using the genomic DNA as a template with primer pairs F1/F2 and F3/F4, respectively (Additional file 12: Figure S5, upper). The hygromycin (HYG) and green fluorescent protein (GFP) cassette was amplified with primer pairs M1/M2 from the binary vector pCT74. Three respective products were ligated through three rounds of PCR for the ccp1 strain construction.
For the mutant complementation (Comp) construction, the ORF of the CCP1 gene was amplified with the 5′-flanking sequences using primer pair F1/CR from the Foc genomic DNA. The product was overlapped with a neomycin cassette (amplified using primers C1/C2), and the 3′-flanking sequences (amplified using primers CF3/F4) through PCR using primer pairs F1/F4 (Additional file 12: Figure S5, lower). Briefly, in the first round, three fragments were amplified and extracted. In the second round, all three amplified products with the same quantity of total DNA (> 800 ng) were amplified using LA Taq polymerase (Takara, Japan) without primers to produce the full-length fusion PCR products. In the third round, primers F1 and F4 were used to yield the mutant and the complementation construction. All primers are shown in Additional file 20: Table S11.
The first- and third-round amplification reactions consisted of 35 cycles of 10 s at 98 °C, 5 s at 55 °C, and 72 °C, 25 s. The second-round reactions consisted of 31 cycles of 5 min at 94 °C, 30 s at 94 °C, 55 °C, 30 s and a final extension of 7 min at 72 °C. The final reaction was the same with the first round except for the 30 s cycle length at 72 °C due to the longer products.
PEG-mediated transformation of Foc protoplast
About 600 mg of fresh-cultured hyphae from Foc TR4 was catalyzed with 5 ml of mixed lysate (40 mg driselase, 40 mg lyticase and 20 mg snailase in 2 ml NaCl (0.8 M)) at 28 °C, for 2–3 h. After cleavage, the precipitate was collected at 2800 rpm for 10 min, washed with 40 ml of STC (1.2 mM sorbitol, 50 mM CaCl2, 10 mM Tris, pH 7.5), centrifuged at 3200 rpm for 10 min, and then redissolved in 200 μl of STC for transformation.
Each fusion PCR product (4 μg), i.e. the mutant or complementation cassette, was added into 200 μl STC containing catalyzed Foc TR4, and let to stand for 20 min. After that, the solution was added with 1.25 ml PTC (1000 mg PEG4000 in 2 ml STC) under gently shaking, and then again set for 20 min. The mutant and complemented strains grew on PDA with hygromycin (100 μg/ml) and neomycin (100 μg/ml), respectively, at 28 °C for 4–6 d.
Analysis of CCP1 gene function in the ccp1 and comp strains of Foc
The penetration abilities of all strains were checked by growth on a cellophane membrane  that covered a semicircle of the pathogen plaque, while the other half remained uncovered, at 28 °C for 5–7 d.
Cellulose-Congo Red medium (Sangon, Shanghai, China) was used to track the utilization of cellulose by all strains . Cellulose Congo Red Medium (25.3 g) and agar were dissolved into 1 L dH2O and autoclaved at 121 °C for 15 min.
Oxidative stress was applied by adding H2O2 to a final concentration of 0.1, 0.5%, or 1% to the PDA medium.
Osmotic stress was applied at two concentrations, 1 M or 2 M NaCl, in the PDA medium.
Cell wall selection pressure was tested with two concentrations of sorbitol (1 M and 2 M) in the PDA medium.
The phenotype of banana infected by the ccp1 strains
Healthy roots of the BX were inoculated by all strains according to our previous method . The banana plants were planted as description before for 45 d. We collected 50 banana plants per treatment for phenotypic analysis. The incidence of wilt was the ratio of plants with wilt leaves and brown spots in pseudostem to all treated plants. Statistical significance between treatments was carried out using Student’s test.
Illumina sequencing was completed by Beijing Novogene Bioinformatics Technology Co. Ltd.
This study was supported by the Program of Hainan Association for Science and Technology Plans to Youth R & D Innovation (HAST201630) and Central Public-interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (No. 1630052016005 and 1630052018003).
Conceived and designed the experiments: MP. Performed the experiments and analyzed the data: WBL. Constructed mutants and complemented strains of Foc: XLW and CQL. Conducted GFP-expression Foc inoculation: JBS. Investigated expression levels of genes: SXL.
Ethics approval and consent to participate
Our research is not involved in human subjects, human material, or human data that must be approved by any ethics committee.
Consent for publication
The work described has not been submitted elsewhere for publication, in whole or in part, and all authors listed have approved the manuscript for submission.
The authors declare no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Moffat AS. Crop engineering Goes south. Science. 1999;285(5426):370.PubMedGoogle Scholar
- Ploetz RC. Fusarium wilt of Banana. Phytopathology. 2015;105(12):1512–21.PubMedGoogle Scholar
- Alaouvette C, Couteaudier Y, Louvet J. Research on the resistance of soils to diseases. IX. Population dynamics of Fusarium spp. and Fusarium oxysporum f. sp. melonis in a soil resistant to and a soil conducive to vascular Fusarium disease. Agronomie. 1984;4(8):729–33.Google Scholar
- Ploetz RC. Panama disease: a classic and destructive disease of banana. Plant Health Progress. 2000;65(3):18–22.Google Scholar
- Wang ZY, Li DB. Genetic diversity in Fusarium oxysporum. Acta Phytopathologica Sinica. 2000;30(3):193–9.Google Scholar
- Chakraborty N, Chandra S, Acharya K. Biochemical basis of improvement of defense in tomato plant against Fusarium wilt by CaCl2. Physiol Mol Biol Plants. 2017;23(3):581–96.PubMedPubMed CentralGoogle Scholar
- Singh PP, Shin YC, Park CS, Chung YR. Biological control of fusarium wilt of cucumber by chitinolytic bacteria. Phytopathology. 1999;89(1):92–9.PubMedGoogle Scholar
- Wang P, Su L, Qin L, Hu B, Guo W, Zhang T. Identification and molecular mapping of a Fusarium wilt resistant gene in upland cotton. Theor Appl Genet. 2009;119(4):733–9.PubMedGoogle Scholar
- Berg NVD, Berger DK, Hein I, Birch PRJ, Wingfield MJ, Viljoen A. Tolerance in banana to fusarium wilt is associated with early up-regulation of cell wall-strengthening genes in the roots. Mol Plant Pathol. 2007;8(3):333–41.Google Scholar
- Berrocal-Lobo M, Molina A. Arabidopsis defense response against Fusarium oxysporum. Trends Plant Sci. 2008;13(3):145–50.PubMedGoogle Scholar
- Bai TT, Xie WB, Zhou PP, Wu ZL, Xiao WC, Zhou L, Sun J, Ruan XL, Li HP. Transcriptome and expression profile analysis of highly resistant and susceptible banana roots challenged with Fusarium oxysporum f. sp. cubense tropical race 4. PLoS One. 2013;8(9):e73945.PubMedPubMed CentralGoogle Scholar
- Li CY, Deng GM, Yang J, Viljoen A, Jin Y, Kuang RB, Zuo CW, Lv ZC, Yang QS, Sheng O, et al. Transcriptome profiling of resistant and susceptible Cavendish banana roots following inoculation with Fusarium oxysporum f. sp. cubense tropical race 4. BMC Genomics. 2012;13(1):374.PubMedPubMed CentralGoogle Scholar
- Li CQ, Shao JF, Wang YJ, Li WB, Guo DJ, Yan B, Xia YJ, Peng M. Analysis of banana transcriptome and global gene expression profiles in banana roots in response to infection by race 1 and tropical race 4 of Fusarium oxysporum f. sp. cubense. BMC Genomics. 2013;14(1):851.PubMedPubMed CentralGoogle Scholar
- Wang Z, Zhang J, Jia C, Liu J, Li Y, Yin X, Xu B, Jin Z. De novo characterization of the banana root transcriptome and analysis of gene expression under Fusarium oxysporum f. sp. cubense tropical race 4 infection. BMC Genomics. 2012;13(1):650.PubMedPubMed CentralGoogle Scholar
- Guo L, Han L, Yang L, Zeng H, Fan D, Zhu Y, Feng Y, Wang G, Peng C, Jiang X, et al. Genome and transcriptome analysis of the fungal pathogen Fusarium oxysporum f. sp. cubense causing banana vascular wilt disease. PLoS One. 2014;9(4):e95543.PubMedPubMed CentralGoogle Scholar
- Westermann AJ, Gorski SA, Vogel J. Dual RNA-seq of pathogen and host. Nat Rev Microbiol. 2012;10(9):618–30.PubMedGoogle Scholar
- Teixeira PJPL, Thomazella DPT, Reis O, do Prado PFV, do Rio MCS, Fiorin GL, Jose J, Costa GGL, Negri VA, Mondego JMC, et al. High resolution transcript profiling of the atypical biotrophic interaction between Theobroma cacao and the fungal pathogen Moniliophthora perniciosa. Plant Cell. 2014;26(11):4245–69.PubMedPubMed CentralGoogle Scholar
- Burgess DJ. Gene expression: host-pathogen duels revealed by dual RNA-seq in vivo. Nat Rev Genet. 2017;18(3):143.PubMedGoogle Scholar
- Hwang SC, Ko WH. Cavendish banana cultivars resistant to Fusarium wilt accquired through somaclonal variation in Taiwan. Plant Dis. 2004;88(6):580–8.PubMedGoogle Scholar
- Liu SQ, Liang ZH, Huang ZH, Huang YX: Breeding of a new banana variety Nongke no. 1 with resistant to Fusarium wilt (Foc TR4). Guangdong agricultural sciences 2007, 1:30–32.Google Scholar
- Xiao R, Zhu Y, Li Y, Huang BL. Green fluorescent protein gene transformation on Fusarium oxysporum f. sp. niveum strain, FOV-135. Fujian Journal of Agricultural Sciences. 2009;24(6):521–4.Google Scholar
- Li WB, Li CQ, Li SX, Peng M. Long noncoding RNAs that respond to Fusarium oxysporum infection in ‘Cavendish’ banana (Musa acuminata). Sci Rep. 2017;7(1):16939.PubMedPubMed CentralGoogle Scholar
- Li W, Li C, Sun J, Peng M. Metabolomic, biochemical and gene expression analyses reveal the underlying responses of resistant and susceptible banana species during early infection with Fusarium oxysporum f. sp. cubense. Plant Dis. 2017;101(4):534–43.PubMedGoogle Scholar
- Ling J, Zeng F, Cao Y, Zhang J, Chen G, Mao Z. Identification of a class of CFEM proteins containing a new conserved motif in Fusarium oxysporum. Physiol mol plant. 2015;89(1):41–8.Google Scholar
- Apel K, Hirt H. Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004;55:373–99.PubMedGoogle Scholar
- Moktali V, Park J, Fedorova-Abrams ND, Park B, Choi J, Lee YH, Kang S. Systematic and searchable classification of cytochrome P450 proteins encoded by fungal and oomycete genomes. BMC Genomics. 2012;13(1):525.PubMedPubMed CentralGoogle Scholar
- Coleman JJ, Mylonakis E. Efflux in fungi: la piece de resistance. PLoS Pathog. 2009;5(6):e100486.Google Scholar
- Bolwell GP. Role of active oxygen species and NO in plant defence responses. Curr Opin Plant Bio. 1999;2(4):287–94.Google Scholar
- Able AJ. Role of reactive oxygen species in the response of barley to necrotrophic pathogens. Protoplasma. 2003;221(1–2):137–43.PubMedGoogle Scholar
- Kumar J, Huckelhoven R, Beckhove U, Nagarajan S, Kogel K-H. A compromised Mlo pathway affects the response of barley to the necrotrophic fungus Bipolaris sorokiniana (teleomorph: Cochliobolus sativus) and its toxins. Phytopathology. 2001;91(2):127–33.PubMedGoogle Scholar
- Shetty NP, Mehrabi R, Lutken H, Haldrup A, Kema GHJ, Collinge DB, Jorgensen HH. Role of hydrogen peroxide during the interaction between the hemibiotrophic fungal pathogen Septoria tritici and wheat. New Phytol. 2007;174(3):637–47.PubMedGoogle Scholar
- Wang XL, Li CQ, Yang JH, Li WB, Sun JB, Peng M. Gene knockout and phenotype analysis of the catalase-1 in Fusarium oxysporum f. sp. cubense race 4. Chinese Journal of Tropical Crops. 2017;38(2):335–42.Google Scholar
- Charizanis C, Juhnke H, Krems B, Entian KD. The mitochondrial cytochrome c peroxidase Ccp1 of Saccharomyces cerevisiae is involved in conveying an oxidative stress siganl to the transcription factor Pos9 (Skn7). Mol Gen Genet. 1999;262(3):437–47.PubMedGoogle Scholar
- Zhou N, Tootle TL, Glazebrook J. Arabidopsis PAD3, a gene required for camalexin biosynthesis, encodes a putative cytochrome P450 monoxygenase. Plant Cell. 1999;11(12):2419–28.PubMedPubMed CentralGoogle Scholar
- Rodriguez MA, Cabrera G, Gozzo FC, Eberlin MN, Godeas A. Clonostachys rosea BAFC3874 as a Sclerotinia sclerotiorum antagonist: mechanisms involved and potential as a biocontrol agent. J Appl Microbiol. 2011;110(5):1177–86.PubMedGoogle Scholar
- Kullnig C, Mach RL, Lorito M, Kubicek CP. Enzyme diffusion from Trichoderma atroviride (=T. harzianum P1) to Rhizoctonia solani is a prerequisite for triggering of Trihoderma ech42 gene expression before mycoparasitic contact. Appl Environ Microb. 2000;66(5):2232–4.Google Scholar
- Manca C, Paul S, Iii CEB, Freedman VH, Kaplan G. Mycobacterium tuberculosis catalase and peroxidase activities and resistance to oxidative killing in human monocytes in vitro. Infect Immun. 1999;67(1):74–9.PubMedPubMed CentralGoogle Scholar
- Seib KL, Wu HJ, Kidd SP, Apicella MA, Jennings MP, Mcewan AG. Defense against oxidative stress in Neisseria gonorrhoeae: a system tailored for a challenging environment. Microbiol Mol Biol Rev. 2006;70(2):344.PubMedPubMed CentralGoogle Scholar
- Vlot AC, Dempsey DA, Klessiq DF. Salicylic acid, a multifaceted hormone to combat disease. Annu Rev Phytopathol. 2009;47(1):177–206.PubMedGoogle Scholar
- Thaler JS, Humphrey PT, Whiteman NK. Evolution of jasmonate and salicylate signal crosstalk. Trends Plant Sci. 2012;17(5):260–70.PubMedGoogle Scholar
- Robert-Seilaniantz A, Grant M, Jones JD. Hormone crosstalk in plant diease and defense: more than just jasmonate-salicylate antagonism. Annu Rev Phytopathol. 2011;49(1):317–43.PubMedGoogle Scholar
- Chen Z, Kloek AP, Cuzick A, Moeder W, Tang D, Innes RW, Klessig DF, McDowell JM, Kunkel BN. The Pseudomonas syringae type III effector AvrRpt2 functions downstream or independently of SA to promote virulence on Arabidopsis thaliana. Plant J. 2004;37(4):494–504.PubMedGoogle Scholar
- Park JE, Park JY, Kim YS, Staswick PE, Jeon J, Y J, Kim SY, Kim J, Lee YH, Park CM. GH3-mediated auxin homeostasis links growth regulation with stress adaptation response in Arabidopsis. J Biol Chem. 2007;13(282):10036–46.Google Scholar
- Xu J, Audenaert K, Hofte M, De Vleesschauwer D. Abscisic acid promotes susceptibility to the rice leaf blight pathogen Xanthomonas oryzae pv oryzae by suppressing salicylic acid-mediated defenses. PLoS One. 2013;8(6):e67413.PubMedPubMed CentralGoogle Scholar
- Jiang CJ, Shimono M, Sugano S, Kojima M, Yazawa K, Yoshida R, Inoue H, Hayashi N, Sakakibara H, Takatsuji H. Abscisic acid interacts antagonistically with salicylic acid signaling pathway in rice-Magnaporthe grisea interaction. Mol Plant-Microbe Interact. 2010;23(6):791–8.PubMedGoogle Scholar
- Meguro A, Sato Y. Salicylic acid antagonizes abscisic acid inhibition of shoot growth and cell cycle progression in rice. Sci Rep. 2014;4:4555.PubMedPubMed CentralGoogle Scholar
- Lin KC, Bushnell WR, Szabo LJ, Smith AG. Isolation and expression of a host response gene family encoding thaumatin-like proteins in incompatible oat-stem rust fungus interactions. Mol Plant-Microbe Interact. 1996;9(6):511–22.PubMedGoogle Scholar
- Specht CA, Liu YL, Robbins PW, Bulawa CE, Iartchouk N, Winter KR, Riggle PJ, Rhodes JC, Dodge CL, Culp DW, et al. The chsD and chsE genes of Aspergillus nidulas and their roles in chitin synthesis. Fungal Genet Biol. 1996;20(2):153–67.PubMedGoogle Scholar
- Mellado E, Aufauvrebrown A, Gow NAR, Holden DW. The Aspergillus fumigatus chsC and chsG genes encode class III chitin synthases with different functions. Mol Microbiol. 2010;20(3):667–79.Google Scholar
- Bowen AR, Chen-wu JL, Momany M, Young R, Szaniszlo PJ, Robbins PW. Classification of fungal chitin synthases. P Natl Acad Sci USA. 1992;89(2):519–23.Google Scholar
- Munro CA, Gow NAR. Chitin synthesis in human pathogenic fungi. Med Mycol. 2001;39(Suppl 1):41–53.PubMedGoogle Scholar
- Madrid MP, Di PA, Roncero MI. Class V chitin synthase determines pathogenesis in the vascular wilt fungus Fusarium oxysporum and mediates resistance to plant defence compounds. Mol Microbiol. 2003;47(1):257–66.PubMedGoogle Scholar
- Takeshita N, Yamashita S, Ohta A, Horiuchi H. Aspergillus nidulans class V and VI chitin synthases CsmA and CsmB, each with a myosin motor-like domain, perform compensatory functions that are essential for hyphal tip growth. Mol Microbiol. 2006;59(5):1380–94.PubMedGoogle Scholar
- Idnurm A, Howlett BJ. Pathogenicity genes of phyto-pathogenic fungi. Mol Plant Pathol. 2001;2(4):241–55.PubMedGoogle Scholar
- Niture SK, Kumar AR, Pant A. Role of glucose in production and repression of polygalacturonase and pectate lyase from phytopathogenic fungus Fusarium moniliforme, NCIM 1276. World J Microb Biot. 2006;22(9):893–9.Google Scholar
- Guo W, Gonzalez-Candelas L, Kolattukudy PE. Cloning of a novel constitutively expressed pectate lyase gene pelB from Fusarium solani f. sp. pisi (Nectria haematococca, mating type VI) and characterization of the gene product expressed in Pichia pastoris. J Bacteriol. 1995;177(24):7070–7.PubMedPubMed CentralGoogle Scholar
- Hogenhout SA, Van der Hoorn RA, Terauchi R, Kamoun S: Emerging concepts in effector biology of plant-associated organisms. Mol Plant Microbe In 2009, 22(2):115–122.PubMedGoogle Scholar
- Czislowski E, Fraser-Smith S, Zander M, O'Neil WT, Meldrum RA, Tran-Nguyen LTT, Batley J, Aiken EAB. Investigation of the diversity of effector genes in the banana pathogen, Fusarium oxysporum f. sp. cubense, reveals evidence of horizontal gene transfer. Mol Plant Pathol. 2017;19(5):1155–71.PubMedGoogle Scholar
- Meldrum RA, Fraser-Smith S, Tran-Nguyen LTT, Daly AM, Aitken EAB. Presence of putative pathogenicity genes in isolates of Fusarium oxysporum fsp cubense from Australia. Australs, Plant Pathology. 2012;41(5):551–7.Google Scholar
- Gao S, Gold SE, Glenn AE. Characterization of two catalase-peroxidase-encoding genes in Fusarium verticillioides reveals differential responses to in vitro versus in planta oxidative challenges. Mol Plant Pathol. 2017;19(5):1127–39.PubMedGoogle Scholar
- Levine A, Tenhaken R, Dixon R, Lamb C. H2O2 from the oxidative burst orchestrates the plant hypersensitive desease resistance response. Cell. 1994;79(4):583–93.PubMedGoogle Scholar
- Yin V, Shaw GS, Konermann L. Cytochrome c as a peroxidase: activation of the precatalytic aative state by H2O2-induced covalent modifications. J Am Chem Soc. 2017;139(44):15701–9.PubMedGoogle Scholar
- Kaya A, Mariotti M, Gladyshev VN. Cytochrome c peroxidase facilitates the beneficial use of H2O2 in prokaryotes. Proc Natl Acad Sci U S A. 2017;114(33):8678–80.PubMedPubMed CentralGoogle Scholar
- Kathiresan M, Martins D, English AM. Respiration triggers heme transfer from cytochrome c peroxidase to catalase in yeast mitochondria. Proc Natl Acad Sci U S A. 2014;111(49):17468–73.PubMedPubMed CentralGoogle Scholar
- Martins D, Kathiresan M, English AM. Cytochrome c peroxidase is a mitochondrial heme-based H2O2 sensor that modulates antioxidant defense. Free Radic Biol Med. 2013;65(6):541–51.PubMedGoogle Scholar
- Hendrixson DR, Dirita VJ. Identification of Campylobacter jejuni genes involved in commensal colonization of the chick gastrointestinal tract. Mol Microb. 2004;52(2):471–84.Google Scholar
- Giles SS, Perfect JR, Cox GM. Cytochrome c peroxidase contributes to the antioxidant defense of Cryptococcus neoformans. Fungal Genet Biol. 2005;42(1):20–9.PubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.PubMedPubMed CentralGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.PubMedPubMed CentralGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.PubMedPubMed CentralGoogle Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34(2):374–8.PubMedGoogle Scholar
- Araki F, Miyagi Y. Effects of fungicides on penetration by Pyricularia oryzae as evaluated by an improved cellophane method. J Pestic Sci. 1977;2(4):457–61.Google Scholar
- Ye JY. A new differential medium for cellulose decomposing microorganisms. Microbiology. 1997;24(4):251–2.Google Scholar