A Colletotrichum graminicola mutant deficient in the establishment of biotrophy reveals early transcriptional events in the maize anthracnose disease interaction
© Torres et al. 2016
Received: 28 October 2015
Accepted: 26 February 2016
Published: 8 March 2016
Colletotrichum graminicola is a hemibiotrophic fungal pathogen that causes maize anthracnose disease. It progresses through three recognizable phases of pathogenic development in planta: melanized appressoria on the host surface prior to penetration; biotrophy, characterized by intracellular colonization of living host cells; and necrotrophy, characterized by host cell death and symptom development. A “Mixed Effects” Generalized Linear Model (GLM) was developed and applied to an existing Illumina transcriptome dataset, substantially increasing the statistical power of the analysis of C. graminicola gene expression during infection and colonization. Additionally, the in planta transcriptome of the wild-type was compared with that of a mutant strain impaired in the establishment of biotrophy, allowing detailed dissection of events occurring specifically during penetration, and during early versus late biotrophy.
More than 2000 fungal genes were differentially transcribed during appressorial maturation, penetration, and colonization. Secreted proteins, secondary metabolism genes, and membrane receptors were over-represented among the differentially expressed genes, suggesting that the fungus engages in an intimate and dynamic conversation with the host, beginning prior to penetration. This communication process probably involves reception of plant signals triggering subsequent developmental progress in the fungus, as well as production of signals that induce responses in the host. Later phases of biotrophy were more similar to necrotrophy, with increased production of secreted proteases, inducers of plant cell death, hydrolases, and membrane bound transporters for the uptake and egress of potential toxins, signals, and nutrients.
This approach revealed, in unprecedented detail, fungal genes specifically expressed during critical phases of host penetration and biotrophic establishment. Many encoded secreted proteins, secondary metabolism enzymes, and receptors that may play roles in host-pathogen communication necessary to promote susceptibility, and thus may provide targets for chemical or biological controls to manage this important disease. The differentially expressed genes could be used as ‘landmarks’ to more accurately identify developmental progress in compatible versus incompatible interactions involving genetic variants of both host and pathogen.
KeywordsRNA-Seq Plant disease Fungal development Transcriptional profiling Biotrophic development
The fungus Colletotrichum graminicola Ces. Wils. is the causal agent of anthracnose leaf blight and anthracnose stalk rot diseases of maize . ASR is one of the most economically important maize diseases, and is estimated to cause billions of dollars in losses annually in the United States . A better understanding of the molecular mechanisms of C. graminicola pathogenicity to maize leaves and stalks might lead to improved methods for disease management.
Most plant pathogens can be classified into one of two groups based on their feeding strategies. Necrotrophic pathogens induce host cell death in advance of their growth, and then feed on the dead plant tissues . Biotrophic pathogens, in contrast, invade living host cells and reprogram them, using a variety of molecular signals, to divert nutrients for their own use [3, 4]. Some pathogens, including C. graminicola, have a hemibiotrophic lifestyle that appears to be intermediate between these two extremes. After mechanically penetrating the plant epidermis via a melanized appressorium, C. graminicola grows initially as a biotroph, producing thick primary hyphae that invade living host cells, and are separated from the host cytoplasm by a membrane. Later, it switches to necrotrophic growth, producing thinner secondary hyphae that colonize dead cells and are no longer surrounded by a membrane [5–11]. Host tissue collapse and necrotic symptoms occur only during the necrotrophic phase of development [1, 4, 11–14].
The genome of C. graminicola includes a large number of genes that are associated with the production of secondary metabolites (SM) , a feature that is reminiscent of many necrotrophic plant pathogens that secrete phytotoxic compounds to kill plant cells ahead of colonization [16, 17]. However, C. graminicola has not been observed to kill host cells in advance, and each new cell is invaded while it is still alive . This behavior is more like that of biotrophic pathogens, which suppress host cell death and defense pathways by producing a wide variety of secreted proteins known as “effectors” [18–21]. Fungal effectors are typically characterized as small, secreted proteins (SSPs), usually induced in planta, which function to facilitate pathogen colonization [15, 22–25]. Fungal SSP effectors are often enriched in cysteine residues (SSP-CRs): cysteine-rich proteins may be more stable and resistant to plant proteases during infection . The C. graminicola genome encodes hundreds of putative SSP and SSP-CR effectors [15, 27]. However, plant defense responses are reportedly activated during the earliest stages of infection of maize leaves, when the fungus still appears to be growing biotrophically  and cells die rapidly after biotrophic invasion , showing that biotrophy in C. graminicola differs fundamentally from that in true biotrophic pathogens.
Biotrophic invasion by C. graminicola continues at the edges of the expanding lesion, even as the pathogen switches to necrotrophic growth in the center . Thus biotrophy and necrotrophy coexist in C. graminicola lesions. A similar growth pattern has been observed in the closely related sorghum anthracnose pathogen, C. sublineola [29, 30]. Contrasting with these two species, the biotrophic hyphae of C. higginsianum, a pathogen of Arabidopsis thaliana and other Brassicaceae, persist in only one cell before making a complete switch to necrotrophy for invasion of subsequent cells . Analysis of differential gene expression patterns in biotrophic hyphae of C. higginsianum suggested that they functioned primarily to produce SSPs and SMs to support the establishment of biotrophy and the subsequent switch to necrotrophy [15, 31]. A preliminary description of the transcriptome of C. graminicola during the biotrophic phase of development in planta was presented previously . However, in contrast with C. higginsianum, it was not possible to obtain detailed information about differential patterns of gene expression in the biotrophic hyphae of C. graminicola, probably because of the asynchronous nature of the infection.
In the present work, we have continued our analysis of our in planta RNA-Seq transcriptome data from C. graminicola by utilizing improved methods that have allowed us to map more reads to the fungal genome, and by applying a more rigorous statistical analysis of differential gene expression during pre-penetration, biotrophic, and necrotrophic phases of development. We developed a “Mixed Effects” Generalized Linear Model (GLM), using the edgeR  package of Bioconductor . The Bioconductor is an open source project that includes tools for high-throughput genomic data analysis. Bioconductor packages are based on the R language, and support object-oriented frameworks, visualization tools, and vignette documents. R is the programming language of choice for statistical and computing applications , and its data structure and codes are used by Bioconductor packages. We chose the edgeR package due to its well-known capability in RNA-Seq analysis, especially its powerful GLM-based differential expression analysis .
We applied the same methods to compare the in planta transcriptome of the pathogenic wild-type strain of C. graminicola (WT) with that of a non-pathogenic mutant strain (MT) that is impaired specifically during early biotrophic development [11, 36]. This strain has an insertional mutation in a gene (Cpr1) that is predicted to encode one subunit of the signal peptidase complex [36, 37]. The MT germinates and produces appressoria normally on maize leaves and leaf sheaths, but penetration is delayed, and once inside it fails to progress beyond the first biotrophically colonized host cell . Comparison of the MT with the WT allowed us to characterize the transcriptional activity of C. graminicola during penetration, and during early versus later stages of biotrophy, in more detail than previously possible in this species. This analysis enabled us to identify genes that are expressed specifically at these critical points in disease establishment, including genes encoding putative SSP effectors, SM enzymes, membrane receptors, and transcription factors. These genes can be subjected to more detailed functional analyses in the future in order to understand their roles in biotrophic development. The products of some of these genes may provide novel targets for improved control of this important disease.
Results and discussion
Overview of the RNA sequencing results and statistical analysis
A total of ~3.5 × 108 Illumina sequencing reads were obtained from samples consisting of the WT and MT strains at three different stages of in planta development (treatments): pre-penetration appressoria (AP); biotrophy (BT); and necrotrophy (NT: WT only) (Additional file 1: Figure S1; Additional file 2: Table S1). Most of the reads were derived from the plant transcriptome (data not included), but a total of 2.2 × 107 reads (6.2 %) could be mapped to the fungal genome (Genbank accession GSE71919) (Additional file 2: Table S1). The percentage of mapped WT reads in the current study was increased to 9.5 %, compared with only 5.6 % reported for the previous analysis  (Genbank accession PRJNA151285) (Additional file 2: Table S1). It should be noted that the previous study did not include the MT strain.
Identification and validation of differentially expressed genes
Counts Per Million (CPM) is the default data filtering method of the edgeR package. As an alternative to the CPM, we applied our own filtering method, mainly because CPM did not account adequately for our sample replicates. In the first round of our filtering method, we chose genes with a total of at least 20 reads across all replicates of at least one treatment. In the next step, genes that had at least 3 treatments with a total of at least 15 reads across all replicates were selected. We determined that the lists of filtered genes generated by our method and by the default CPM method were more than 96 % identical, and also in each comparison there was more than 93 % agreement. However our method avoided screening out of genes that actually had sufficient read depth to be included, when replicates were taken into account.
Eighty-four percent of the C. graminicola predicted genes (10,028/12,006) had sufficient read depth to include after filtering. Differential expression analysis of these genes via edgeR identified 2412 statistically significant differentially expressed genes across the various comparisons (FDR ≤0.05, Log2 fold change >2) (Additional file 2: Table S2). This is compared with 2619 differentially expressed genes identified in our previous study . However, the former analysis included only P value and Log2 fold change. Application of the more rigorous FDR criterion to the previous dataset reduces the number of statistically significant differentially expressed genes to 1855.
The RNA-Seq data for fourteen differentially expressed genes were validated by using quantitative real-time reverse-transcription polymerase chain reaction (qRT-PCR). Transcript log2 fold-changes (AP vs BT, BT vs NT, AP vs NT) measured by RNA-Seq and qRT-PCR had a positive correlation of R 2 = 0.8604, and a simple linear model slope of y = 1.036 (Additional file 1: Figure S2), supporting the reliability of the RNA-Seq data.
One in five fungal genes is differentially expressed during development of the WT in planta
Summary of the number of genes with significantly different expression (FDR ≤0.05, Log2 fold change >2) between different fungal stage comparisons
WTAP vs. WTBT
WTAP vs. WTNT
WTBT vs. WTNT
MTAP vs. MTBT
WTAP vs. MTAP
WTBT vs. MTBT
Genes that are differentially expressed specifically during the WT biotrophic phase of development
PHI database hit
LIP1 Botrytis cinerea
Hydrophobic surface binding protein A
Bacterial transferase hexapeptide (six repeats) (probable transcription factor)
GzZC087 Gibberella zeae
Fungal trichothecene efflux pump (TRI12)
TRI12 Fusarium sporotrichoides (toxin associated)
AKT1 Alternaria alternate (toxin associated)
CobW/HypB/UreG, nucleotide binding domain
Fungal specific transcription factor domain
GzZC297 Gibberella zeae
Heterokaryon incompatibility protein (HET)
Related to O-methylsterigmatocystin oxidoreductase Gibberella zeae (toxin associated)
Predicted nuclear conserved hypothetical protein
Predicted cytosolic conserved hypothetical protein
Major Facilitator Superfamily
GzCCHC002 Gibberella zeae
Fumarylacetoacetate (FAA) hydrolase family
Alpha/beta hydrolase family
PTH11 Magnaporthe oryzae
ATPase family associated with various cellular activities (AAA)
Major Facilitator Superfamily
MGG_10702 Magnaporthe oryzae
Small secreted protein putative effector, similar to MGG_02647
Pyridine nucleotide disulphide oxidoreductase
DEP4 Alternaria brassicicola (toxin associated)
Secreted WSC domain protein
Sulfate transporter family
Predicted mitochondrial protein, unique to C. graminicola
Short chain dehydrogenase
MFP1 Magnaporthe oryzae
Predicted nuclear protein, unique to C. graminicola
Relatively few genes were differentially expressed in the MT
Genes that are differentially expressed during the transition from AP to BT in the mutant
ChLae1 Cochliobolus heterostrophus
Glycosyl hydrolase family 3 N terminal domain
TOM1 Septoria lycopersici
BcBOT1 (related CND5) Botrytis cinerea
Sedlin, N-terminal conserved regions
Predicted small mitochondrial protein, unique to C. graminicola
Secreted aminopeptidase M28A
Predicted plasma membrane protein
PTH11 Magnaporthe oryzae
Glycosyl hydrolases family 11
XYL2 Cochliobolous carbonum
Secreted Zinc carboxypeptidase M14A
Glycosyl hydrolase family 61
GzOB021 Gibberella zeae
Pregnancy-associated plasma protein-A: secreted metalloprotease M43B
MEP1 Coccidioides posadasii
Extracellular conserved hypothetical protein
Glycosyl hydrolase family 10
Endo-1,4-beta-xylanase M. oryzae PHI_2208
Dienelactone hydrolase family
Alpha/beta hydrolase family
Cazyme, (CE1). SM Cluster 22
Zinc binding dehydrogenase
SM Cluster 22
Fungal Cellulose Binding Domain
FAED1 Gibberella zeae
Conserved SSP, cysteine rich
Aconitase family (aconitase hydratase)
LYSF Aspergillus fumigatus
The Cpr1 gene is not differentially expressed in planta in the MT or WT strains
The MT strain has an insertion of foreign DNA into the 3′ untranslated region of the Cpr1 gene (GLRG_04964), which is predicted to encode a non-catalytic component of the signal peptidase . Analysis of the RNA-Seq data indicated that Cpr1 was expressed at similar levels in the WT and in the MT in planta across all phases of development.
Membrane-protein and secreted-protein genes are over-represented among the differentially expressed genes
Relative representation of differentially-expressed genes encoding potential pathogenicity-associated categories of proteins in different comparisons
Comparison (number of genes)
Total Genome (12,006)
687 (6 %)
3224 (27 %)
110 (1 %)
704 (6 %)
1476 (12 %)
300 (3 %)
WT Early Genes (159)
24 (15 %)
55 (35 %)
2 (1 %)
14 (9 %)
13 (8 %)
14 (9 %)
WT Late Genes (440)
42 (10 %)
144 (33 %)
6 (1 %)
77 (18 %)
48 (11 %)
26 (6 %)
2 (40 %)
2 (40 %)
2 (9 %)
9 (39 %)
1 (4 %)
1 (4 %)
5 (22 %)
43 (15 %)
97 (33 %)
6 (2 %)
25 (9 %)
30 (10 %)
20 (7 %)
62 (13 %)
186 (40 %)
29 (6 %)
84 (18 %)
70 (15 %)
19 (4 %)
WTAP_WTBT Total (760)
105 (14 %)
283 (37 %)
35 (5 %)
109 (14 %)
100 (13 %)
39 (5 %)
62 (20 %)
92 (30 %)
10 (3 %)
36 (12 %)
31 (10 %)
17 (5 %)
73 (11 %)
227 (33 %)
15 (2 %)
112 (17 %)
91 (13 %)
37 (5 %)
WTBT_WTNT Total (992)
135 (14 %)
319 (32 %)
25 (3 %)
148 (15 %)
122 (12 %)
54 (5 %)
17 (23 %)
18 (24 %)
2 (3 %)
9 (12 %)
8 (11 %)
7 (5 %)
41 (29 %)
3 (2 %)
9 (6 %)
15 (10 %)
12 (8 %)
WTAP_MTAP Total (218)
24 (11 %)
59 (27 %)
5 (2 %)
18 (8 %)
23 (11 %)
12 (6 %)
29 (15 %)
77 (40 %)
20 (10 %)
34 (18 %)
28 (15 %)
2 (1 %)
43 (8 %)
194 (37 %)
4 (1 %)
50 (10 %)
63 (12 %)
32 (6 %)
WTBT_MTBT Total (714)
72 (10 %)
271 (38 %)
24 (3 %)
84 (12 %)
91 (13 %)
34 (5 %)
4 (20 %)
10 (50 %)
3 (15 %)
6 (30 %)
2 (10 %)
Prepenetration appressoria rely on stored nutrients, and experience significant oxidative stress
A Gene Ontology (GO) analysis of transcripts that were differentially expressed during the transition from pre-penetration WT appressoria (AP) to biotrophy (BT) was performed by using Blast2GO  (Additional file 2: Table S5). An enhanced expression of genes in AP encoding proteins associated with oxidation and export of potentially toxic compounds and defense against stress, especially oxidative stress, suggests that the fungus is exposed to antimicrobial host defenses even at this early stage. Genes involved in proline metabolism were elevated in AP relative to BT (Additional file 2: Table S5). In C. trifolii, proline was a potent antioxidant under nutrient-limiting conditions . Proline accumulation may protect appressoria from ROS accumulation during pre-penetration stages. ROS production by the host cells can be detected during all stages of development in infected leaf tissues [11, 28], but the specific fungal response to oxidative stress appears to be stage-specific, with different antioxidant genes expressed at different phases of development. ROS are also important regulators of normal fungal development, and induction of some antioxidant genes may be related to developmental regulation of the fungus, rather than a response to host defenses [40, 41].
Before host penetration, plant pathogens face starvation conditions and rely on the metabolism of stored compounds. Carbon and nitrogen starvation have been reported as environmental cues for the expression of many pathogenicity-associated genes [42–44]. Spores of C. graminicola contain large lipid bodies, which provide the nutrients required for germination, appressorial maturation, and host penetration . Expression of isocytrate lyase (ICL) (GLRG_04226), a key enzyme of the glyoxylate cycle, was significantly increased in AP relative to BT and NT. ICL is elevated during early infection stages in several other pathogenic fungi, and it is essential for pathogenicity of the hemibiotrophs Magnaporthe oryzae and Leptosphaeria maculans [46–48]. ICL and other components of the glyoxylate cycle are required for appressorial maturation and function in both C. orbiculare and M. oryzae [49, 50] and play important roles in pathogenicity of many fungal and bacterial pathogens of both plants and animals .
Genes involved in transferase activities, including dimethylallyl tryptophan synthases (DMATs) and other aromatic prenyltransferases, were overrepresented in AP relative to BT (Additional file 2: Table S5). Other SM associated genes, including genes involved in melanin biosynthesis, were also more highly expressed. Some of these SM may play roles in stress mitigation during AP, while others may be important for inducing host susceptibility. Several genes associated with adhesion to hydrophobic surfaces (PF12296) were more highly expressed during AP (Additional file 2: Table S5). These proteins may be involved in appressorial maturation and attachment, which are important for appressorial function and signaling .
The biotrophic phase is characterized by an increase in secretory activity, especially the production of secreted proteases
The transition from AP to BT was characterized by an increase in the expression of genes encoding secreted proteins (Additional file 2: Table S5). These included proteins that have been implicated in interactions with other organisms, and numerous secreted proteases belonging to classes that are known to play important roles in plant pathogenicity (Additional file 2: Table S5) [26, 53–55]. Subtilisins are serine proteases that can degrade cell wall proteins and plant defense proteins . Targeted deletion of the subtilisin gene Spm1 severely compromised pathogenicity of M. oryzae to rice plants  Another up-regulated protease was a metalloprotease belonging to the M. oryzae AVR-Pita avirulence gene family . Carboxypeptidases and aspartic proteases were also represented. Mutant rice plants that transiently expressed a carboxypeptidase inhibitor from potato were highly resistant to infection by the blast fungus M. oryzae and the root pathogen Fusarium verticillioides . Aspartic proteases have been identified in Botrytis cinerea-inoculated carrots, cabbage and grapes . Application of the purified enzymes from B. cinerea induced cell death in carrot cell cultures, and inhibition of the enzymatic activity significantly reduced virulence. It appears that in C. graminicola, production of a wide variety of secreted proteases is an important feature of biotrophic establishment and colonization. Many of these proteases are likely to target host defense mechanisms, as has been observed in other pathosystems .
In both C. higginsianum and C. orbiculare, biotrophic hyphae are proposed to function primarily as secretory organs for the production of SSP effectors and SM for modification of the plant environment, induction of compatibility, and suppression of programmed host cell death (PCD) [31, 60]. Host cell death could be delayed directly, by targeting host PAMP receptors, and also indirectly, by down-regulating the production of potential elicitors of PCD.
Among the 23 genes that were specifically down-regulated during BT in C. graminicola, some were homologs of genes encoding known toxin-associated proteins, eg. TRI12 from F. sporotrichoides , AKT1 from Alternaria alternata , and DEP4 from A. brassicciola  (Table 2). There was also a homolog of a cytochrome P450 monoxygenase gene from F. graminearum (FGSG_00007) that, when knocked out, resulted in greatly increased toxin production and virulence of that pathogen to wheat . It is possible that these genes are down-regulated in biotrophic hyphae to avoid triggering host cell death prematurely. Among the other genes that were specifically reduced in expression were two transcription factors, suggesting the possibility of phase-specific negative regulation of transcription.
The transition to BT was associated with an increase in expression of genes encoding secreted hydrolases, including cutinases, pectate lyases, and chitin deacetylases (Additional file 2: Table S5). Cutinases are involved in the interactions of fungi with host cuticle, and play an important role in signaling to trigger fungal development, and in the activation of host defense . Pectate lyases target pectin, the major polysaccharide comprising the middle lamella, and they are important pathogenicity factors in several pathosystems [66–68]. Chitin deacetylases convert chitin in the cell walls of primary hyphae of Colletotrichum to chitosan, and it has been suggested that this prevents it from being recognized by the host plant and triggering PAMP-mediated resistance during biotrophic colonization [69, 70].
Expression of multiple acid phosphatases was increased in BT relative to AP. Phytases are one class of acid phosphatases that convert organic forms of phosphorous into inorganic phosphate that can be utilized by the fungus . A predicted extracellular phytase (GLRG_06496) and multiple membrane-bound phosphate transporters (GLRG_01384, GLRG_03681, GLRG_006247, GLRG_09610, and GLRG_10529) were significantly increased in expression during both BT and NT relative to AP. This suggests that the in planta environment is relatively poor in phosphorous . A previous study identified phytase as a highly expressed activity in biotrophic hyphae purified from maize stalks by laser capture . The siderophore biosynthesis gene SID1 (GLRG_06540) and the gene encoding NPS6, which is responsible for synthesis of secreted siderophores (GLRG_08065), were both more highly expressed in BT versus AP, rather than being specifically down-regulated as previously reported . Aryl-sulfatase, responsible for sulfur catabolism, was reportedly highly expressed during penetration, and then progressively reduced in expression during biotrophy and necrotrophy in C. gloeosporioides , but our RNA-Seq data did not support differential expression of the homolog of this gene (GLRG_11452) in C. graminicola. The differences in our findings may relate to differences in the strains or species, or in the methods used (RNA-Seq versus semi-quantitative PCR or fluorescent reporters), or to difficulties with precisely and reproducibly identifying the biotrophic phase of development, which is non-synchronous and relatively brief. Our data suggest that biotrophic hyphae of C. graminicola actively assimilate phosphorus, sulfur, and iron from the plant host.
Nearly all of the genes that were increased in expression in BT relative to AP were also increased in NT. Only five genes were increased in expression specifically during BT. One of these was associated with vitamin B6 biosynthesis. Genes involved with biosynthesis of other B vitamins, B1 and B12, were significantly increased during both BT and NT relative to AP. The B vitamins are cofactors of many different enzymes, including some involved in carbon and amino acid metabolism [75, 76], suggesting that there is an increase in these activities during the growth of C. graminicola in planta. Vitamin B1 biosynthetic genes were significantly induced in planta in the biotrophs Puccinia triticina and Uromyces fabae [77, 78]. Vitamin B6, also known as pyridoxine, has also been associated with antioxidant activities and resistance to oxidative stress in plant and fungi [79, 80]. Active detoxification mechanisms could be important to overcome defense mechanisms and establish a successful biotrophic interaction.
Two other biotrophy-specific genes encode members of large classes of membrane proteins that may play roles in signaling. One was a potential G-protein coupled receptor protein in the Pth11 family [81, 82]. This family is very large in C. graminicola, with 50 members. Ten Pth11-family genes were differentially expressed in the WT, with eight transcribed early, and two late during the infection. These CFEM-domain receptors may play roles in the specific recognition of plant signals, and mediation of developmental transitions during the Colletotrichum-maize interaction. The second biotrophy-specific signaling gene encoded a putative heterokaryon incompatibility protein. Two additional HET genes were up-regulated during both BT and NT, relative to AP. HET proteins interact in heterologous pairs to trigger cell death in fungi, including C. graminicola . Expression of HET proteins during growth in planta may serve to regulate intra-mycelial interactions, and protect the colonized tissues from encroachment by other C. graminicola strains.
The transition to necrotrophy was characterized by an increase in production of hydrolytic enzymes, and in the utilization of carbon and nitrogen from the host
The transition to NT in the WT was associated primarily with an increase in activities that have been related to degradation of host cell walls  (Additional file 2: Table S6). These included hydrolase activities (including cellulases, laccases, and peptidases), carbohydrate catabolism, and cellulose binding (Additional file 2: Table S6). Cell wall degrading enzymes (CWDE) are utilized for host penetration and colonization, and also to obtain nutrients from plant polymers . Necrotrophic plant pathogens normally secrete large quantities of CWDE during host colonization, whereas biotrophs produce relatively few [17, 21] Categories related to stress response, particularly hypoxia, were decreased in NT compared with BT, signifying that the necrotrophic hyphae were experiencing less oxygen deprivation than the biotrophic hyphae. This may be related to the relative disorganization of the host tissues in NT versus BT, and lack of a host membrane surrounding the necrotrophic hyphae.
There are a large number of C. graminicola genes that were primarily or only expressed during NT (late genes) (Additional file 2: Table S3). A majority of these encoded CWDE of various classes. This is consistent with cytological evidence that dissolution of host cell walls becomes evident only after the transition to necrotrophy [7, 11]. Among the late genes were two glutamate importers (GLRG_03954 and GLRG_04076), suggesting increased availability of this preferred nitrogen source after induction of host cell death. Availability of preferred carbon sources is also indicated by a significant decrease in the expression of ICL during NT relative to AP and BT. These results indicate that necrotrophic hyphae of C. graminicola are more active than biotrophic hyphae in the uptake and utilization of nutrients from the host tissues. Degradation of the tissue during this phase likely facilitates this uptake by making nutrients more accessible and available to the fungal mycelium.
Number of genes in each comparison that were annotated using Blast2GO
Patterns of expression of predicted effector genes
Annotation of differentially expressed genes encoding secreted proteins
Total (Annotated + Non-annotated)
Categories of fungal effectors and effector families in the C. graminicola genome and transcriptome
GH61 motif SSPa
CFEM motif SSPb
Conserved fungal effectors in the C. graminicola transcriptome
ClH1 LysM (GLRG_02947)
ClH1 LysM (GLRG_07767)
ChNLP2 (NPP1) (GLRG_11600)
Hemibiotrophic pathogens produce SSPs during later phases of infection that induce PCD of host cells, and are thought to mediate the switch to necrotrophy [31, 86, 90]. For example, there are six genes in C. higginsianum encoding proteins that belong to the NPP1 family of PCD-inducing effectors found in Phytophthora species . There are homologs in C. graminicola for all but one of these C. higginsianum proteins (Table 8). C. graminicola homologs of ChNLP2, ChNLP3 and ChNLP5 are differentially expressed. Only CgNLP1 and CgNLP2 in C. graminicola share the amino acids residues crucial for PCD-inducing activity [99, 100]. CgNLP2 is most highly expressed during NT. In contrast, the CgNLP3 and CgNLP5 transcripts are more abundant during AP, like their homologs in C. higginsianum. Another conserved SSP gene, GLRG_01192, is a homolog of CgEC91, a hypersensitive-response inducing protein effector, induced during the switch to necrotrophy in C. higginsianum . Expression of this gene was increased in C. graminicola BT relative to AP.
The patterns of expression of SSPs in C. graminicola suggest that, as in other hemibiotrophic pathogens, “early” SSPs, produced by appressoria and biotrophic hyphae, promote compatibility by subverting host recognition and PCD, while “late” SSPs, produced late during biotrophic development, are involved in the switch to necrotrophy by promoting host PCD.
Patterns of expression of predicted secondary metabolism (SM) clusters
The C. graminicola genome is unusually rich in SM genes . SM genes generally occur as part of co-regulated gene clusters. Forty-two predicted SM clusters are encoded by the C. graminicola genome . The RNA-Seq data provided evidence for co-regulation of only seven of these clusters in planta, defined as having most or all of the genes in the cluster significantly differentially expressed, and with the same pattern of expression (Additional file 2: Table S8).
Expression of SM genes during NT is expected, given that SM are usually associated with necrotrophic plant pathogens. However, induction of SM during earlier, biotrophic stages of development has been described in other intracellular hemibiotrophs including C. higginsianum ; C. orbiculare ; and M. oryzae [101, 102] and it has been suggested that SM expressed during early infection stages act to suppress host defenses, rather than to kill tissues. SM expressed later in the infection process may induce host PCD, or protect necrotic tissues from microbial competitors.
The identities of the products generated by most of the differentially regulated C. graminicola SM clusters can only be guessed, since the clusters do not closely resemble known clusters from other fungi. However, cluster 38 was an exception because it is identical, both in gene content and in gene order, to the RADS cluster of Pochonia chlamydosporia , which is responsible for production of the antifungal SM radicicol (aka. monorden). Monorden, and its biosynthetic intermediates monocillins I, II and III, , are among the few known fungal metabolites that are produced by C. graminicola in maize stalks, and in vitro analyses demonstrated their antifungal activity against other maize stalk-rot and foliar pathogens . Monorden inhibits heat-shock protein (Hsp) 90, by competition with ATP for the binding site required for its activation . It was suggested that C. graminicola may secrete these SM products during early penetration and biotrophic stages, in order to suppress basal host defense responses. However, our transcriptional analysis revealed that cluster 38 is significantly induced during WTNT, suggesting instead that monorden plays a role in defending necrotic tissue from microbial competitors.
Analysis of the MT transcriptome provides additional clues to the nature of biotrophy
Very few biotrophy-specific genes could be identified based on analysis of the WT RNA-Seq data alone. This is probably because BT consists of a mixture of cell types, including numerous pre-penetration appressoria, biotrophic primary hyphae that are just entering living host cells, and intercalary primary hyphae behind the advancing colony front that are occupying cells that are already dead or dying. Nearly all of the genes that were increased in expression in BT relative to AP were also increased in NT. NT is also a mixed culture, with necrotrophic hyphae produced in the center of the colony, but persistence of biotrophy at the colony edges. This lack of synchronicity would be expected to mute potential cell-specific differences in gene expression. Penetration and the establishment of biotrophic hyphae in the living host cell are of great interest as potential determinants of the disease outcome, but these processes are transient and nonsynchronous, and thus the transcriptional activities associated with them are not clearly revealed by the WT RNA-Seq analysis.
The MT strain allowed us to dissect the nature of BT and biotrophic establishment in C. graminicola in more detail. The MT is blocked early during the production of the primary hyphae: it does not advance to adjacent cells and thus, MTBT consists only of the initial invasive primary hyphae, and includes no intercalary hyphae . Only 20 genes were differentially expressed in the transition from MTAP to MTBT, compared with 760 in the WT transition. EdgeR MDS plots indicated that MTBT was more similar to WTAP and MTAP than to WTBT. Patterns of the most highly expressed genes in each condition also supported this conclusion (Additional file 2: Table S7).
The establishment of BT is expected to involve the activities of numerous SSP effectors and SM that facilitate invasion by promoting susceptibility and inhibiting host PCD. Among the 20 genes that were upregulated in MTBT relative to MTAP, two were members of the DMAT SM cluster 22 (Table 3). There were also four SSPs. One is a homolog of the glycosyl hydrolase XYL2 from C. carbonum. In M. oryzae there was a reduction in virulence when the homolog of this gene was knocked out . GLRG_06286 encodes a secreted metalloprotease and is homologous to MEP1 from Coccidioides posadasii. The MEP1 protein prevents host detection by digesting surface antigens from mice cells. When this gene was mutated, virulence was reduced . The product of this gene in C. graminicola may be involved in blocking activity of PAMP receptors of maize. The remaining two differentially expressed SSPs are uncharacterized. One is conserved in multiple Colletotrichum species, while the other is found only in C. graminicola and its close relative C. sublineola. All four SSP genes are also increased in expression in WTBT compared with WTAP, and remained up-regulated in WTNT, suggesting that their function is not specific to penetration or to the establishment of biotrophy. A gene homologous to ChLae1, a master regulator thought to function by altering heterochromatin, and that contributes to host selective toxin production, pathogenicity, and adaptation to oxidative stress in Cochliobolus heterostrophus , was also upregulated in MTBT relative to MTAP. This suggests that the early transition to BT involves chromatin remodeling and a resulting shift toward expression of pathogenicity-specific genes.
Most of these 20 genes were also more highly expressed in WTBT versus WTAP, suggesting that the transition to biotrophy is initiated normally in the MT. There were five exceptions: One of these genes encoded a homolog of LysF from Aspergillus fumigatus, which functions in lysine biosynthesis and was necessary for pathogenicity of that fungus . Another gene encoded a sedlin-domain protein. Sedlin is a component of the secretory pathway that functions in protein transport from the endoplasmic reticulum to the golgi . Increased activity of the secretory pathway in biotrophic hyphae would be consistent with their putative function as secretory cells for protein effectors [31, 60]. It is possible that all five of these genes are only transiently up-regulated during the establishment of biotrophy, and that this was not detected in the more heterogeneous WTBT samples.
Comparison of WT samples to MT samples collected at the same stage of development provided even more detailed information, specifically about events that occur during the crucial penetration and early, versus late, biotrophic phases.
Genes expressed during penetration include numerous hydrolases, phytase, and cysteine-rich SSP effectors
Although germination and appressorial induction and maturation occurred at the same rate, penetration by the MT was delayed compared with the WT . Thus, when the MTAP and WTAP samples were collected, the WTAP had already initiated penetration while the MTAP had not. Comparisons between WTAP versus MTAP indicated that WTAP is enriched in genes encoding secreted hydrolases and a phytase, suggesting an increase in these activities coinciding with the initiation of host wall penetration (Additional file 2: Table S9). Seventeen SSPs, nine of which were cysteine-rich, were more highly expressed in the WTAP versus the MTAP (Table 7). Five of these were also among the 100 most highly expressed genes in WTAP (Table 7). Four were early genes, suggesting that they are transiently up-regulated and may have a specific function in penetration, while nearly all of the remainder were expressed at similar levels in both WTAP and WTBT, and then down-regulated during WTNT, suggesting that they have functions that are important both during penetration and biotrophic colonization. These SSPs include the LysM effectors that are presumed to be important for sequestering chitin, and four LS-SSPs that may be involved in specific interactions with the host.
Only seven SSP effectors were more highly expressed in MTAP, suggesting these were induced earlier, prior to the initiation of penetration: none were LS-SSPs, and only one was cysteine-rich (Table 7). Thus, there appears to be a bias toward increased expression of SSPs, especially SSP-CRs and LS-SSPs, coinciding with the initiation of penetration. It is interesting that six of the seven SSP genes that were higher in MTAP versus WTAP were also significantly increased in expression during WTBT and WTNT compared with WTAP. This suggests that these effectors function during all phases of development, not just in appressoria, and also that the expression of these genes might be down-regulated during penetration and early biotrophy.
MTAP was also more active in expression of eight of the genes contained in SM cluster 18. Cluster 18 as predicted contains 24 genes, and ten of these are homologous to genes in cluster 10 in C. higginsianum . However, the predicted cluster 18 probably represents two separate clusters in C. graminicola, one homologous to the C. higginsianum cluster 10 (which our RNA-Seq and qRT-PCR data both confirm is not expressed in any phase); and a second without a homolog in C. higginsianum. The eight differentially expressed genes in MTAP belong to the second cluster. The backbone of this second cluster is a non-reducing PKS gene, GLRG_08632. This SM cluster is apparently expressed during appressorial formation prior to the initation of penetration, and produces an unknown SM product.
Early versus late biotrophy: late biotrophy is characterized by increased expression of genes encoding secreted proteases and transmembrane receptors and transporters
Relationships among gene expression patterns for genes that are differentially expressed in both the WTAP vs. WTBT comparisons, and in the WTBT vs. MTBT “Up” or “Down” refers in each case to the first term, in the comparison relative to the second. For example, WTAP_WTBT_up means that the genes are more highly expressed in WTAP than in WTBT
In a comparison of GO-Terms associated with WTBT versus MTBT, the WTBT was enriched in extracellular and integral membrane proteins, as well as terms associated with host-pathogen interactions (including several subtilases and proteins involved in peroxide detoxification), transmembrane transport, cellular homeostasis, peptidase activities, and SM (Additional file 2: Table S10). Secreted proteases were over-represented among this “late BT” group, suggesting that these protease functions are more important later in biotrophy versus during the initial invasion of living host cells.
MTBT over-expressed genes involved in the glyoxylate cycle and melanin biosynthesis, an indication that the MTBT sample contained a larger proportion of transcripts from appressoria. This is consistent with the extremely limited development of intracellular hyphae in MTBT compared with WTBT. MTBT over-expressed 43 SSP genes, including 12 SSP-CRs. Twenty-four of the genes are “early” genes that are specifically more highly expressed in WTAP. These are likely to be expressed primarily in the appressoria in the MTBT sample. The other 19 SSP genes were expressed at similar levels in WTAP and WTBT, with five of these decreasing during WTNT, suggesting that they were specific to late AP and early BT.
In comparison with the MTBT, the WTBT was enriched in proteins that might function in communication between host and pathogen, including secreted and integral membrane proteins and proteases potentially involved in evasion of host immune response. Proteins involved in peroxide detoxification that may protect the intracellular hyphae from host ROS, transmembrane transporters to allow uptake and egress of compounds including nutrients or SM, and proteins involved in maintenance of cellular homeostasis, which may preserve hyphal integrity in the dynamic environment of the living host cell, were also over-represented. Twenty-nine SSPs were overexpressed in WTBT including two pectate-lyase domain SSPs that might target the host middle lamella, and the NPP1 homolog ChNLP2, which may be involved in inducing PCD. Cells die rapidly once they have been invaded by C. graminicola primary hyphae, possibly due to activation of these necrosis-inducing proteins and production of endogenous elicitors by disruption of the host cellular structure.
Maize leaf sheaths respond differently to WT and MT infections
It is possible that the delay in penetration by MT appressoria is related to a failure to secrete proteins that are necessary to prepare the epidermal cell for invasion and establishment of biotrophy, and/or to process signals from the host that trigger further development of the pathogen. Host responses to the MT and WT strains at different stages of infection were studied by using qRT-PCR to evaluate expression of maize genes previously associated in the literature with the response to C. graminicola or biotrophic pathogens [6, 28, 112, 113].
In barley, expression of the Bax-family PCD inhibitor Bi1  was induced during infection by the obligate biotrophic pathogen Blumeria graminis f.sp. hordei . Silencing of Bi1 significantly reduced infection, while overexpression increased the susceptibility of barley to biotrophs [112, 115]. The maize homolog was increased in expression during WTAP, suggesting that appressoria may produce factors that induce this gene, thereby inhibiting hypersensitive cell death prior to penetration, and preparing the host cell for biotrophic invasion. The maize pathogenicity-related proteins PR1, PR3, and PR5 are induced as a result of activation of the salicylic acid (SA) defense pathway . The Bowman-birk trypsin inhibitor (Bti) is a marker for induction of jasmonic acid (JA) . SA-dependent pathways are typically deployed against biotrophic plant pathogens, while JA-mediated pathways are usually associated with plant responses against necrotrophic pathogens [118, 119]. The role of SA and JA signaling defense pathways is less clear for hemibiotrophic plant pathogens. In maize seedlings inoculated with C. graminicola, SA accumulated at 36 hpi (during biotrophy), while JA could be detected by 96 hpi (during necrotrophy) . Interestingly, expression of Bti was significantly induced in maize leaves as early as 24 hpi with the biotrophic U. maydis, although silencing of this gene did not affect pathogen colonization . Our data suggest that both defense pathways are activated in leaf sheaths in response to WT appressoria prior to penetration.
Our study indicates that there is a continuum of activities that occur during colonization of maize by C. graminicola, and that the boundaries we have drawn between “AP”, “BT” and “NT” are artificial. About 20 % of the fungal genes are differentially transcribed in “waves” across the processes of appressorial maturation, penetration, and colonization. Expression of genes involved in nutrient uptake and utilization suggest that there is a shift during development from the utilization of stored lipids during appressorial development and early biotrophy toward the uptake and utilization of nutrients including amino acids and sugars from the host during late biotrophy and necrotrophy. Secreted protein, SM protein, and membrane receptor genes are over-represented among the differentially expressed genes, suggesting that the fungus engages in an intimate and dynamic conversation with the host, beginning prior to penetration. This communication process is likely to involve reception of plant signals that trigger subsequent developmental progress in the fungus, as well as the production of signals that induce responses in the host. Later phases of biotrophy are more similar to necrotrophy, with increased production of secreted proteases, inducers of PCD, hydrolases, and membrane bound transporters for the uptake and egress of potential toxins, signals, and nutrients. The initiation of biotrophy on the other hand, when the fungus establishes itself in living host cells, is similar to the penetration phase that occurs late during appressorial development, with increased production of hydrolases targeting the host cell surface and involved in signaling and modification of fungal cell wall to avoid triggering host recognition, as well as the secretion of many SSPs, particularly lineage specific and cysteine-rich SSPs, that are likely to be involved in suppression of host defense responses, including hypersensitive cell death. The mutant provides us with a “snapshot” of the transcriptional events that occur during the critical points of penetration and biotrophic establishment, as it appears to pause the infection process at this phase so that it can be captured by the transcriptome analysis. Identification of a series of highly expressed genes that can serve as markers for different phases of disease development is possible based on the data presented in this study. For example, the following series of effector genes is predicted to represent a temporal sequence that spans the entire infection process: GLRG_08629; GLRG_01735; GLRG_02947; GLRG_00338; GLRG_08975; GLRG_11600; and GLRG_09214. These genes could potentially be used as ‘landmarks’ to identify and compare developmental progress across different host-pathogen strain combinations or different conditions.
These RNA-Seq data do not provide many insights into the root cause of the mutant phenotype, since most genes appear to be transcribed normally in the MT, including the mutant Cpr1 gene itself. The mutant phenotype is probably related to post-transcriptional events, which must await further study. Nonetheless, detailed investigation of the genes that are differentially expressed between the MT and WT can help to identify potential virulence factors (e.g. SM and SSPs), receptors, and transcriptional regulators that play critical roles in the establishment of biotrophy, and can also help to differentiate genes that are expressed during invasion of living host cells from genes that are expressed later in biotrophy, in cells that are dead and dying, in preparation for the switch to necrotrophy. These genes must be subjected to further study to elucidate their precise roles, and to determine if they represent viable targets for disease management.
Plants and fungal strains
The susceptible maize inbred Mo940 was used for this study. Plants were grown to the V3 stage in the greenhouse, with a 14-h day length, in 3.8 × 21 cm plastic Conetainers (Super SC-10 UV stabilized Stuewe & Sons, Inc. Oregon, USA), containing a mixture of three parts Pro-Mix BX (Premiere Horticulture, Ltd, Riviere du Loup, PQ, Canada) to two parts of sterile topsoil. Plants were watered daily to saturation and fertilized two to three times per week with a solution of 150 ppm of Peters 20-10-20 (Scotts-Sierra Horticultural Products Co., Marysville, OH), beginning one week after germination.
C. graminicola strain M1.001 was the wild type (WT) . The mutant (MT) was derived from M1.001 by restriction-enzyme mediated insertional mutagenesis [36, 37]. The MT is nonpathogenic to maize stalks and leaves due to an insertion in the 3′UTR of the Cpr1 gene, which is predicted to encode one component of the microsomal signal peptidase . A complemented strain (Cpr1-C), generated by transformation of the mutant with a 3.6 kb fragment of genomic DNA containing the WT Cpr1 gene, is fully restored in pathogenicity and comparable to the WT strain both in vitro and in planta . The fungal strains were routinely cultured on potato dextrose agar (PDA, Difco) at 23 °C under continuous light.
Sample preparation and RNA extraction
C. graminicola falcate spores were harvested and inoculated on detached maize leaf sheaths as described , except that two 20-μl inoculum drops were applied to each leaf sheath, approximately 1 cm from either end, and at least 3 cm apart. Sheaths containing mature pre-penetration appressoria (WTAP, approximately 20 hpi); intracellular biotrophic hyphae, before symptoms or secondary hyphae were visible (WTBT, approximately 36 hpi); or necrotrophic hyphae, in which initial browning of the tissue and secondary hyphae were visible (WTNT, approximately 60 hpi), were collected. The MT does not progress to the production of necrotrophic hyphae, and so only the AP and BT stages were collected for that strain. Thus, our dataset consisted of five genotype-timepoint combinations (i.e. treatments): WTAP, WTBT, WTNT, MTAP, and MTBT. Each infected leaf sheath was inspected under the microscope to confirm the developmental stage. For the BT and NT samples, each sheath was wiped gently with a moistened sterile cotton swab to remove unattached spores and superficial mycelia. The process of trimming, cleaning, and examination did not take more than two minutes per sheath. Six trimmed leaf sheaths (~100 mg) were pooled into each tube, flash-frozen in liquid nitrogen, and maintained at −80 °C until RNA extraction.
Total RNA was extracted by crushing the frozen tissue, followed by grinding in TRIzol (1 ml per 100 mg sample) (Invitrogen). Samples were incubated for 7 h in isopropanol followed by 2 h in 100 % ethanol, both incubated at −20 °C. Samples were purified and treated with DNAse using the RNeasy Plant Mini Kit (Qiagen), according to the manufacturer’s protocol as modified , and then resuspended in 50 μl of RNAse-free water. Extracts from approximately 30 pieces of leaf sheath (500 mg) were pooled for each biological replicate. Each biological replicate represented an independent inoculation experiment.
Three-hundred micrograms of total RNA from three biological replicates of each treatment (with the exception of MTBT, which had only two usable biological replicates) were submitted for sequencing to AgriLife Genomics and Bioinformatics Services (Texas A&M University). Libraries were prepared for each sample by using the Illumina TruSeqTM RNA Sample Preparation Kit and the manufacturer’s instructions. Data were generated from ten lanes of Illumina GAII sequencing, in two separate runs, with barcoding to multiplex biological replicates. Read lengths were 76 bp (including 7 bp for the barcode adaptor). For the first run, eight lanes of a flowcell were used, and lane five was spiked with 1 % PhiX as a control. For the second run, in which additional data for the WTAP was obtained, two lanes were used. A total of four lanes (i.e. technical replicates) of data were generated for WTAP; one lane each was produced for WTBT and WTNT; and two lanes each were produced for MTAP and MTBT. The technical replicates for each treatment were pooled. Data were processed using the Illumina software CASAVA-1.7.0 for base calling and de-multiplexing, and the final results were stored as individual files for each sample in FASTQ format.
Alignment to the reference genome
The sequencing reads were mapped to the 12,006 annotated genes of C. graminicola strain M1.001 (ACOD01000001) using the CLC Genomics Workbench (GWB) RNA-Seq analysis tool (http://www.clcbio.com). CLC GWB reports gene expression as “Reads Per Kilobase of Transcript per Million Mapped Reads” (RPKM) . The RPKM value has been shown to underperform compared to alternative normalizations that may be applied to “Total Gene Reads” . Pairwise comparison capabilities provided by CLC did not take full advantage of our experimental design. Therefore, we extracted the “Total Gene Reads” field for each replication of each treatment, and subjected this output to further statistical analysis as described below.
We used a “Mixed Effects” Generalized Linear Model (GLM) to account for biological replication, developmental stage, and genotype effects. The mixed-effect GLM allowed us to leverage data from all developmental stages simultaneously in a unified analysis framework. The GLM improves the statistical power and inference capabilities, and can protect against multiple testing as compared to the pairwise testing paradigm through a single prerequisite “no effects” test. We used the GLM implemented in the R package, edgeR (version 3.0.8) . The methods implemented in edgeR assume that the underlying distribution for the discrete count measures of the next generation sequencing (NGS) is a negative binomial (NB). In conjunction with an internal, model-based normalization method that has been shown to be superior to RPKM, edgeR estimates mean and variance of the NB distribution for each gene , and then proceed in a gene by gene manner for experimental design differential expression analysis.
Our GLM specification included main effects and interactions across genotypes (WT and MT) and developmental stage (AP, BT, and NT). Based on these coefficient estimates, six biologically informative contrasts were estimated: (1) WTAP vs. WTBT; (2) WTAP vs. WTNT; (3) WTBT vs. WTNT; (4) WTAP vs. MTAP; (5) WTBT vs. MTBT; (6) MTAP vs. MTBT.
Normalization and statistical testing
The filtered genes and their counts were converted into a DGEList R object. The default normalization method of edgeR – Trimmed Mean of M values (TMM) – was then applied to the data using the “calcNormFactors” EdgeR function. TMM uses one sample as the reference, and assumes that most genes are not differentially expressed. As the result of the normalization, the scaling factors are calculated for each sample library , and these scale factors are then used for model-based normalization by EdgeR.
The edgeR assumption for the underlying distribution of the RNA-Seq data is a negative binomial distribution. The dispersion for each gene is using the quantile-adjusted conditional maximum likelihood (qCML) implemented by the “estimateGLMTagwiseDisp” EdgeR function. Statistical testing was carried out with the “glmFit” and “glmLRT” R functions. The genes were classified as up- or down-regulated using the “decideTestsDGE” EdgeR function. Multiple testing was addressed using the false discovery rate (FDR)  as implemented in the p.adjust R function.
Heatmaps were built by using Genesis (Release 1.7.6) , and represent log2-fold changes of a transcript in each of the three fungal developmental stages, relative to the average expression across all stages.
RNA extraction from inoculated leaf sheaths
Leaf sheaths were inoculated with WT, MT and complemented mutant strains, incubated, trimmed, and flash-frozen for RNA extraction as described above. Additional treatments were also tested in the (q)RT-PCR experiments. Tissue inoculated with the Cpr1-C strain, and sheaths that were mock-inoculated with water, were also collected, trimmed, and frozen for RNA extraction as described previously.
RNA extraction from in vitro appressoria
Appressoria of the WT, MT, and Cpr1-C strains were produced in vitro on polystyrene Petri dishes as described in , with some modifications. C. graminicola spores were collected and washed three times, and 40 ml of a spore suspension at a concentration of 2 × 104 spores/ml was added to each Petri dish. Twenty hours later, each plate was inspected under the microscope to verify the presence of mature melanized appressoria. Appressoria were broken and scraped from the bottoms of ten Petri plates, using a sterile culture spreader, into a total of five ml of TRIzol. Appressoria collected from 40 Petri plates were combined into each replicate. RNA purification was performed as described above.
cDNA synthesis and cycling reactions
The SuperScript II reverse transcriptase kit (Invitrogen) was used to synthesize the first strand of cDNA from 1 μg of DNase and RNase-treated total RNA in a volume of 20 μl. Primers were designed to amplify 100–200 bp fragments, using PrimerQuest (Integrated DNA Technologies) software. The reaction mix for real-time PCR contained 0.4 mM of each primer, 10 μl of SYBR green PCR Master Mix (AppliedBiosystems), 5 μl of a 1:5 dilution of the cDNA product, and DEPC water to a final volume of 20 μl. Cycling conditions were as follows: 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. The reactions were carried out in fast 96-well reaction plates on the ABI 7900HT fast RT-PCR system (Applied Biosystems). Fungal transcript levels were normalized by using the fungal actin gene as an internal standard, and relative expression was calculated using the Pfaffl method . Maize genes were normalized against the maize actin gene, and expression was calculated relative to mock-inoculated plants.
Functional annotation and gene ontology
Nucleotide sequences similar to differentially expressed genes (P < 0.05), were identified by BLASTx searches of the non-redundant database (E value = 1 × 10−3) of the Blast2GO suite . Functional characterization and gene ontology (GO) categories for cellular functions, cellular components, and biological processes, were assigned using the Blast2Go platform. The GOSSIP function was utilized to determine GO term enrichment in different comparisons . Manual annotation of specific genes was performed by using BLAST searches against the NCBI databases (E value = 1 × 10−5) and InterproScan analysis. Comparisons to C. higginsianum were performed by using the tools available on the Broad Institute Colletotrichum database (https://www.broadinstitute.org/scientific-community/science/projects/fungal-genome-initiative/colletotrichumgenome-project). Cellular localization of significantly expressed transcripts was predicted by using WoLF PSORT (http://wolfpsort.org/) . To identify protein families, the Pfam database (http://pfam.sanger.ac.uk/) (E value = 1 × 10−5) was used . Transporters were predicted using the Transporters Classification Database (http://www.tcdb.org) (E value = 1 × 10−5) . The identification of carbohydrate active enzymes (CAZymes) was done using the web resource dbCAN (http://csbl.bmb.uga.edu/dbCAN/annotate.php), an automated CAZyme annotation that is based on the classification scheme of CAZyDB [131, 132]. For the classification of putative secreted proteases, the sequences for the secreted proteins predicted by WoLF PSORT were submitted to MEROPS Batch Blast analysis (http://merops.sanger.ac.uk)  as described in O’Connell et al.. BLAST with an e-value cutoff of 1e−5 was used to identify candidate Colletotrichum pathogenicity proteins with similarity to proteins present in the Pathogen-Host Interaction  .
Illumina sequencing was done at the Texas A&M AgriLife Genomics and Bioinformatics Service. The authors would like to thank Dr. Richard Metz, Ms. Etta Nuckles, and Mr. Doug Brown for excellent technical assistance. This is publication No. 16-12-022 of the Kentucky Agricultural Experiment Station and is published with the approval of the Director. This work was supported by the U.S. Department of Agriculture-Cooperative State Research, Education, and Extension Service (USDA-CSREES) grant 2009-34457-20125 (LJV); by the National Institute of Food and Agriculture, U.S. Department of Agriculture Hatch project 0231781 (LJV); and by a University of Kentucky College of Agriculture, Food, and Environment Research Activity Award (LJV).
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.
- Bergstrom GC, Nicholson RL. The Biology of Corn Anthracnose: Knowledge to Exploit for Improved Management. Plant Dis. 1999;83(7):596–608.View ArticleGoogle Scholar
- Frey TJ et al. Fitness Evaluation of, a Locus that Confers Resistance to (Ces.) GW Wils. Using Near-Isogenic Maize Hybrids. Crop Sci. 2011;51(4):1551–63.View ArticleGoogle Scholar
- Horbach R et al. When and how to kill a plant cell: infection strategies of plant pathogenic fungi. J Plant Physiol. 2011;168(1):51–62.PubMedView ArticleGoogle Scholar
- Bancal M-O, Hansart A, Sache I, Bancal P. Modelling fungal sink competitiveness with grains for assimilates in wheat infected by a biotrophic pathogen. Annals of Botany. 2012;110:113–123.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldwin TK et al. The Pathogen-Host Interactions Database (PHI-base) Provides Insights into Generic and Novel Themes of Pathogenicity. Mol Plant Microbe Interact. 2006;19(12):1451–62.PubMedView ArticleGoogle Scholar
- Balmer D et al. Induced resistance in maize is based on organ-specific defence responses. Plant J. 2013;74(2):213–25.PubMedView ArticleGoogle Scholar
- Mims CW, Vaillancourt LJ. Ultrastructural characterization of infection and colonization of maize leaves by Colletotrichum graminicola, and by a C. graminicola pathogenicity mutant. Phytopathology. 2002;92(7):803–12.PubMedView ArticleGoogle Scholar
- Politis D, Wheeler H. Ultrastructural study of penetration of maize leaves by Colletotrichum graminicola. Physiol Plant Pathol. 1973;3(4):465–71.View ArticleGoogle Scholar
- Venard C, Vaillancourt L. Colonization of fiber cells by Colletotrichum graminicola in wounded maize stalks. Phytopathology. 2007;97(4):438–47.PubMedView ArticleGoogle Scholar
- Venard C, Vaillancourt L. Penetration and colonization of unwounded maize tissues by the maize anthracnose pathogen Colletotrichum graminicola and the related nonpathogen C. sublineolum. Mycologia. 2007;99(3):368–77.PubMedView ArticleGoogle Scholar
- Torres MF, Cuadros DF, Vaillancourt LJ. Evidence for a diffusible factor that induces susceptibility in the Colletotrichum–maize disease interaction. Mol Plant Pathol. 2014;15(1):80–93.PubMedView ArticleGoogle Scholar
- Albarouki E, et al. Biotrophy-specific downregulation of siderophore biosynthesis in Colletotrichum graminicola is required for modulation of immune responses of maize. Mol Microbiol. 2014;92(2):338–55.PubMedPubMed CentralView ArticleGoogle Scholar
- Crouch J, et al. The genomics of Colletotrichum, in Genomics of Plant-Associated Fungi: Monocot Pathogens2014. Heidelberg, New York, Dordrecht, London: Springer. p. 69–102.
- Münch S et al. The hemibiotrophic lifestyle of Colletotrichum species. J Plant Physiol. 2008;165(1):41–51.PubMedView ArticleGoogle Scholar
- O’Connell RJ, et al. Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses. Nat Genet. 2012;44(9):1060–5.PubMedView ArticleGoogle Scholar
- Markham JE, Hille J. Host‐selective toxins as agents of cell death in plant–fungus interactions. Mol Plant Pathol. 2001;2(4):229–39.PubMedView ArticleGoogle Scholar
- Amselem J, et al. Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet. 2011;7(8):e1002230.PubMedPubMed CentralView ArticleGoogle Scholar
- Djamei A, et al. Metabolic priming by a secreted fungal effector. Nature. 2011;478(7369):395–8.PubMedView ArticleGoogle Scholar
- Doehlemann G, et al. Pep1, a secreted effector protein of Ustilago maydis, is required for successful invasion of plant cells. 2009.Google Scholar
- Rafiqi M, et al. Challenges and progress towards understanding the role of effectors in plant–fungal interactions. Curr Opin Plant Biol. 2012;15(4):477–82.PubMedView ArticleGoogle Scholar
- Spanu PD, et al. Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism. Science. 2010;330(6010):1543–6.PubMedView ArticleGoogle Scholar
- de Wit PJGM, et al. Fungal effector proteins: past, present and future. Mol Plant Pathol. 2009;10(6):735–47.PubMedView ArticleGoogle Scholar
- Kämper J, et al. Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis. Nature. 2006;444(7115):97–101.PubMedView ArticleGoogle Scholar
- Lo Presti L, et al. Fungal effectors and plant susceptibility. Annu Rev Plant Biol. 2015;66:513–45.PubMedView ArticleGoogle Scholar
- Rooney HCE, et al. Cladosporium Avr2 inhibits tomato Rcr3 protease required for Cf-2-dependent disease resistance. Science. 2005;308(5729):1783–6.PubMedView ArticleGoogle Scholar
- Kamoun S. A catalogue of the effector secretome of plant pathogenic oomycetes. Annu Rev Phytopathol. 2006;44:41–60.PubMedView ArticleGoogle Scholar
- de Jonge R, et al. Conserved fungal LysM effector Ecp6 prevents chitin-triggered immunity in plants. Science. 2010;329(5994):953–5.PubMedView ArticleGoogle Scholar
- Vargas WA et al. Plant defense mechanisms are activated during biotrophic and necrotrophic development of Colletotricum graminicola in maize. Plant Physiol. 2012;158(3):1342–58.PubMedPubMed CentralView ArticleGoogle Scholar
- Wharton PS, Julian AM. A cytological study of compatible and incompatible interactions between Sorghum bicolor and Colletotrichum sublineolum. New Phytol. 1996;134(1):25–34.View ArticleGoogle Scholar
- Wharton PS, Julian AM, O’Connell RJ. Ultrastructure of the infection of Sorghum bicolor by Colletotrichum sublineolum. Phytopathology. 2001;91(2):149–58.PubMedView ArticleGoogle Scholar
- Kleemann J, et al. Sequential delivery of host-induced virulence effectors by appressoria and intracellular hyphae of the phytopathogen Colletotrichum higginsianum. PLoS Pathog. 2012;8(4):e1002643.PubMedPubMed CentralView ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.PubMedPubMed CentralView ArticleGoogle Scholar
- Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.PubMedPubMed CentralView ArticleGoogle Scholar
- Team, RC. R Language Definition, 2000, Available from CRAN sites.
- Law CW, et al. Voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014;15(2):R29.PubMedPubMed CentralView ArticleGoogle Scholar
- Thon MR, et al. CPR1: a gene encoding a putative signal peptidase that functions in pathogenicity of Colletotrichum graminicola to maize. Mol Plant Microbe Interact. 2002;15(2):120–8.PubMedView ArticleGoogle Scholar
- Thon MR, Nuckles EM, Vaillancourt LJ. Restriction enzyme-mediated integration used to produce pathogenicity mutants of Colletotrichum graminicola. Mol Plant Microbe Interact. 2000;13(12):1356–65.PubMedView ArticleGoogle Scholar
- Conesa, A. and S. Götz, Blast2GO: A comprehensive suite for functional analysis in plant genomics. International journal of plant genomics. 2008;2008:12. doi:10.1155/2008/619832.
- Chen C, Dickman MB. Proline suppresses apoptosis in the fungal pathogen Colletotrichum trifolii. Proc Natl Acad Sci U S A. 2005;102(9):3459–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Egan MJ, et al. Generation of reactive oxygen species by fungal NADPH oxidases is required for rice blast disease. Proc Natl Acad Sci. 2007;104(28):11772–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Fang GC, Hanau RM, Vaillancourt LJ. The SOD2 gene, encoding a manganese-type superoxide dismutase, is up-regulated during conidiogenesis in the plant-pathogenic fungus Colletotrichum graminicola. Fungal Genet Biol. 2002;36(2):155–65.PubMedView ArticleGoogle Scholar
- Marzluf GA. Genetic regulation of nitrogen metabolism in the fungi. Microbiol Mol Biol Rev. 1997;61(1):17–32.PubMedPubMed CentralGoogle Scholar
- Coleman M, et al. Starvation-induced genes of the tomato pathogen Cladosporium fulvum are also induced during growth in planta. Mol Plant Microbe Interact. 1997;10(9):1106–9.PubMedView ArticleGoogle Scholar
- Pellier AL, et al. CLNR1, the AREA/NIT2‐like global nitrogen regulator of the plant fungal pathogen Colletotrichum lindemuthianum is required for the infection cycle. Mol Microbiol. 2003;48(3):639–55.PubMedView ArticleGoogle Scholar
- Schadeck RJG, Leite B, de Freitas Buchi D. Lipid mobilization and acid phosphatase activity in lytic compartments during conidium dormancy and appressorium formation of Colletotrichum graminicola. Cell Struct Funct. 1998;23(6):333–40.PubMedView ArticleGoogle Scholar
- Idnurm A, Howlett BJ. Isocitrate lyase is essential for pathogenicity of the fungus Leptosphaeria maculans to canola (Brassica napus). Eukaryot Cell. 2002;1(5):719–24.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee H-S, et al. Inhibition of the pathogenicity of Magnaporthe grisea by bromophenols, isocitrate lyase inhibitors, from the red alga Odonthalia corymbifera. J Agric Food Chem. 2007;55(17):6923–8.PubMedView ArticleGoogle Scholar
- Wang ZY, et al. The glyoxylate cycle is required for temporal regulation of virulence by the plant pathogenic fungus Magnaporthe grisea. Mol Microbiol. 2003;47(6):1601–12.PubMedView ArticleGoogle Scholar
- Asakura M, Okuno T, Takano Y. Multiple Contributions of Peroxisomal Metabolic Function to Fungal Pathogenicity in Colletotrichum lagenarium. Appl Environ Microbiol. 2006;72(9):6345–54.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Z-Y, et al. Functional analysis of lipid metabolism in Magnaporthe grisea reveals a requirement for peroxisomal fatty acid β-oxidation during appressorium-mediated plant infection. Mol Plant Microbe Interact. 2007;20(5):475–91.PubMedView ArticleGoogle Scholar
- Dunn MF, Ramirez-Trujillo JA, Hernández-Lucas I. Major roles of isocitrate lyase and malate synthase in bacterial and fungal pathogenesis. Microbiology. 2009;155(10):3166–75.PubMedView ArticleGoogle Scholar
- Hamer JE, Talbot NJ. Infection-related development in the rice blast fungus Magnaporthe grisea. Curr Opin Microbiol. 1998;1(6):693–7.PubMedView ArticleGoogle Scholar
- Movahedi S, Heale JB. The roles of aspartic proteinase and endo-pectin lyase enzymes in the primary stages of infection and pathogenesis of various host tissues by different isolates of Botrytis cinerea Pers ex. Pers Physiol Mol Plant Pathology. 1990;36(4):303–24.View ArticleGoogle Scholar
- ten Have A, et al. An aspartic proteinase gene family in the filamentous fungus Botrytis cinerea contains members with novel features. Microbiology. 2004;150(7):2475–89.PubMedView ArticleGoogle Scholar
- Valueva TA, Mosolov VV. Role of inhibitors of proteolytic enzymes in plant defense against phytopathogenic microorganisms. Biochemistry (Moscow). 2004;69(11):1305–9.View ArticleGoogle Scholar
- Olivieri F, et al. Characterization of an extracellular serine protease of Fusarium eumartii and its action on pathogenesis related proteins. Eur J Plant Pathol. 2002;108(1):63–72.View ArticleGoogle Scholar
- Oh Y, et al. Transcriptome analysis reveals new insight into appressorium formation and function in the rice blast fungus Magnaporthe oryzae. Genome Biol. 2008;9(5):R85.PubMedPubMed CentralView ArticleGoogle Scholar
- Khang CH, et al. Genome organization and evolution of the AVR-Pita avirulence gene family in the Magnaporthe grisea species complex. Mol Plant Microbe Interact. 2008;21(5):658–70.PubMedView ArticleGoogle Scholar
- Quilis J, et al. A potato carboxypeptidase inhibitor gene provides pathogen resistance in transgenic rice. Plant Biotechnol J, 2007;5(4):537–53.PubMedView ArticleGoogle Scholar
- Gan P. et al. Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi. New Phytol. 2013;197(4):1236–49.PubMedView ArticleGoogle Scholar
- Alexander NJ, McCormick SP, Hohn TM. TRI12, a trichothecene efflux pump from Fusarium sporotrichioides: gene isolation and expression in yeast. Mol Gen Genet. 1999;261(6):977–84.PubMedView ArticleGoogle Scholar
- Tanaka A, et al. Insertional mutagenesis and cloning of the genes required for biosynthesis of the host-specific AK-toxin in the Japanese pear pathotype of Alternaria alternata. Mol Plant Microbe Interact. 1999;12(8):691–702.PubMedView ArticleGoogle Scholar
- Wight WD, et al. Biosynthesis and role in virulence of the histone deacetylase inhibitor depudecin from Alternaria brassicicola. Mol Plant Microbe Interact. 2009;22(10):1258–67.PubMedView ArticleGoogle Scholar
- Gardiner DM, Kazan K, Manners JM. Novel genes of Fusarium graminearum that negatively regulate deoxynivalenol production and virulence. Mol Plant Microbe Interact. 2009;22(12):1588–600.PubMedView ArticleGoogle Scholar
- Serrano, M, et al. The cuticle and plant defense to pathogens. Frontiers in plant science, 2014;5.
- Ben-Daniel B-H, Bar-Zvi D, Tsror L. Pectate lyase affects pathogenicity in natural isolates of Colletotrichum coccodes and in pelA gene-disrupted and gene-overexpressing mutant lines. Mol Plant Pathol. 2012;13(2):187–97.PubMedView ArticleGoogle Scholar
- Cho Y, et al. A Pectate Lyase-Coding Gene Abundantly Expressed during Early Stages of Infection Is Required for Full Virulence in Alternaria brassicicola. PloS One, 2015. 10(5).
- Cnossen-Fassoni A, et al. The pectate lyase encoded by the pecCl1 gene is an important determinant for the aggressiveness of Colletotrichum lindemuthianum. J Microbiol. 2013;51(4):461–70.PubMedView ArticleGoogle Scholar
- El Gueddari NE, et al. Developmentally regulated conversion of surface‐exposed chitin to chitosan in cell walls of plant pathogenic fungi. New Phytologist. 2002;156(1):103–12.View ArticleGoogle Scholar
- O’Connell RJ, Ride JP. Chemical detection and ultrastructural localization of chitin in cell walls of Colletotrichum lindemuthianum. Physiol Mol Plant Pathology. 1990;37(1):39–53.View ArticleGoogle Scholar
- Abelson PH. A potential phosphate crisis. Science. 1999;283:2015.PubMedView ArticleGoogle Scholar
- Mullaney EJ, Daly CB, Ullah A. Advances in phytase research. Adv Appl Microbiol. 1999;47:157–99.View ArticleGoogle Scholar
- Tang W, et al. The application of laser microdissection to in planta gene expression profiling of the maize anthracnose stalk rot fungus Colletotrichum graminicola. Mol Plant Microbe Interact. 2006;19(11):1240–50.PubMedView ArticleGoogle Scholar
- Goodwin PH, Li J, Jin S. Evidence for sulfate derepression of an arylsulfatase gene of Colletotrichum gloeosporioides f. sp. malvae during infection of round-leaved mallow, Malva pusilla. Physiol Mol Plant Pathology. 2000;57(4):169–76.View ArticleGoogle Scholar
- Sohn J, et al. High level activation of vitamin B1 biosynthesis genes in haustoria of the rust fungus Uromyces fabae. Mol Plant Microbe Interact. 2000;13(6):629–36.PubMedView ArticleGoogle Scholar
- Percudani R, Peracchi A. A genomic overview of pyridoxal‐phosphate‐dependent enzymes. EMBO Rep. 2003;4(9):850–4.PubMedPubMed CentralView ArticleGoogle Scholar
- Thara VK, Fellers JP, Zhou JM. In planta induced genes of Puccinia triticina. Mol Plant Pathol. 2003;4(1):51–6.PubMedView ArticleGoogle Scholar
- Struck C, et al. The Uromyces fabae UfAAT3 gene encodes a general amino acid permease that prefers uptake of in planta scarce amino acids. Mol Plant Pathol. 2004;5(3):183–9.PubMedView ArticleGoogle Scholar
- Bilski P, et al. Vitamin B6 (Pyridoxine) and Its Derivatives Are Efficient Singlet Oxygen Quenchers and Potential Fungal Antioxidants. Photochem Photobiol. 2000;71(2):129–34.PubMedView ArticleGoogle Scholar
- Titiz O, et al. PDX1 is essential for vitamin B6 biosynthesis, development and stress tolerance in Arabidopsis. Plant J. 2006;48(6):933–46.PubMedView ArticleGoogle Scholar
- DeZwaan TM, et al. Magnaporthe grisea pth11p is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues. Plant Cell. 1999;11(10):2013–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Kulkarni RD, et al. Novel G-protein-coupled receptor-like proteins in the plant pathogenic fungus Magnaporthe grisea. Genome Biol. 2005;6(3):R24.PubMedPubMed CentralView ArticleGoogle Scholar
- Vaillancourt LJ, Hanau RM. Nitrate-nonutilizing mutants used to study heterokaryosis and vegetative compatibility in Glomerella graminicola (Colletotrichum graminicola). Exp Mycol. 1994;18(4):311–9.View ArticleGoogle Scholar
- Kubicek CP, Starr TL, Glass NL. Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol. 2014;52:427–51.PubMedView ArticleGoogle Scholar
- Catanzariti A-M, et al. Haustorially Expressed Secreted Proteins from Flax Rust Are Highly Enriched for Avirulence Elicitors. Plant Cell. 2006;18(1):243–56.PubMedPubMed CentralView ArticleGoogle Scholar
- Koeck M, Hardham AR, Dodds PN. The role of effectors of biotrophic and hemibiotrophic fungi in infection. Cell Microbiol. 2011;13(12):1849–57.PubMedPubMed CentralView ArticleGoogle Scholar
- Ökmen B, Doehlemann G. Inside plant: biotrophic strategies to modulate host immunity and metabolism. Curr Opin Plant Biol. 2014;20:19–25.PubMedView ArticleGoogle Scholar
- Kleemann J, et al. Identification of soluble secreted proteins from appressoria of Colletotrichum higginsianum by analysis of expressed sequence tags. Microbiology. 2008;154(4):1204–17.PubMedView ArticleGoogle Scholar
- Guyon K, et al. Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum. BMC Genomics. 2014;15(1):336.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee S-J, Rose JKC. Mediation of the transition from biotrophy to necrotrophy in hemibiotrophic plant pathogens by secreted effector proteins. Plant Signal Behav. 2010;5(6):769–72.PubMedPubMed CentralView ArticleGoogle Scholar
- Schulze-Lefert P, Panstruga R. A molecular evolutionary concept connecting nonhost resistance, pathogen host range, and pathogen speciation. Trends Plant Sci. 2011;16(3):117–25.PubMedView ArticleGoogle Scholar
- Pedersen C, et al. Structure and evolution of barley powdery mildew effector candidates. BMC Genomics. 2012;13(1):694.PubMedPubMed CentralView ArticleGoogle Scholar
- Dong S, et al. Effector specialization in a lineage of the Irish potato famine pathogen. Science. 2014;343(6170):552–5.PubMedView ArticleGoogle Scholar
- Mosquera G, et al. Interaction transcriptome analysis identifies Magnaporthe oryzae BAS1-4 as biotrophy-associated secreted proteins in rice blast disease. Plant Cell. 2009;21(4):1273–90.PubMedPubMed CentralView ArticleGoogle Scholar
- Xue C, et al. Two novel fungal virulence genes specifically expressed in appressoria of the rice blast fungus. Plant Cell. 2002;14(9):2107–19.PubMedPubMed CentralView ArticleGoogle Scholar
- Marshall R, et al. Analysis of two in planta expressed LysM effector homologs from the fungus Mycosphaerella graminicola reveals novel functional properties and varying contributions to virulence on wheat. Plant Physiol. 2011;156(2):756–69.PubMedPubMed CentralView ArticleGoogle Scholar
- Mentlak TA, et al. Effector-mediated suppression of chitin-triggered immunity by Magnaporthe oryzae is necessary for rice blast disease. Plant Cell. 2012;24(1):322–35.PubMedPubMed CentralView ArticleGoogle Scholar
- Kombrink A, Thomma BPHJ. LysM effectors: secreted proteins supporting fungal life. 2013.Google Scholar
- Dong S, et al. The NLP toxin family in Phytophthora sojae includes rapidly evolving groups that lack necrosis-inducing activity. Mol Plant Microbe Interact. 2012;25(7):896–909.PubMedView ArticleGoogle Scholar
- Ottmann C et al. A common toxin fold mediates microbial attack and plant defense. Proc Natl Acad Sci. 2009;106(25):10359–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Fudal I, et al. Expression of Magnaporthe grisea avirulence gene ACE1 is connected to the initiation of appressorium-mediated penetration. Eukaryot Cell. 2007;6(3):546–54.PubMedPubMed CentralView ArticleGoogle Scholar
- Collemare J, et al. Biosynthesis of secondary metabolites in the rice blast fungus Magnaporthe grisea: the role of hybrid PKS-NRPS in pathogenicity. Mycol Res. 2008;112(2):207–15.PubMedView ArticleGoogle Scholar
- Wang S, et al. Functional characterization of the biosynthesis of radicicol, an Hsp90 inhibitor resorcylic acid lactone from Chaetomium chiversii. Chem Biol. 2008;15(12):1328–38.PubMedView ArticleGoogle Scholar
- Zhou H, et al. Insights into radicicol biosynthesis via heterologous synthesis of intermediates and analogs. J Biol Chem. 2010;285(53):41412–21.PubMedPubMed CentralView ArticleGoogle Scholar
- Wicklow DT, Jordan AM, Gloer JB. Antifungal metabolites (monorden, monocillins I, II, III) from Colletotrichum graminicola, a systemic vascular pathogen of maize. Mycol Res. 2009;113(12):1433–42.PubMedView ArticleGoogle Scholar
- Roe SM, et al. Structural basis for inhibition of the Hsp90 molecular chaperone by the antitumor antibiotics radicicol and geldanamycin. J Med Chem. 1999;42(2):260–6.PubMedView ArticleGoogle Scholar
- Nguyen QB, et al. Simultaneous silencing of endo‐β‐1, 4 xylanase genes reveals their roles in the virulence of Magnaporthe oryzae. Mol Microbiol. 2011;81(4):1008–19.PubMedView ArticleGoogle Scholar
- Hung C-Y, et al. A metalloproteinase of Coccidioides posadasii contributes to evasion of host detection. Infect Immun. 2005;73(10):6689–703.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu D, et al. ChLae1 and ChVel1 regulate T-toxin production, virulence, oxidative stress response, and development of the maize pathogen Cochliobolus heterostrophus. PLoS Pathog. 2012;8(2):e1002542.PubMedPubMed CentralView ArticleGoogle Scholar
- Liebmann B, et al. Deletion of the Aspergillus fumigatus lysine biosynthesis gene lysF encoding homoaconitase leads to attenuated virulence in a low-dose mouse infection model of invasive aspergillosis. Arch Microbiol. 2004;181(5):378–83.PubMedView ArticleGoogle Scholar
- Venditti R, et al. Sedlin controls the ER export of procollagen by regulating the Sar1 cycle. Science. 2012;337(6102):1668–72.PubMedPubMed CentralView ArticleGoogle Scholar
- Doehlemann G, et al. Reprogramming a maize plant: transcriptional and metabolic changes induced by the fungal biotroph Ustilago maydis. Plant J. 2008;56(2):181–95.PubMedView ArticleGoogle Scholar
- van der Linde K, et al. Systemic virus‐induced gene silencing allows functional characterization of maize genes during biotrophic interaction with Ustilago maydis. New Phytol. 2011;189(2):471–83.PubMedView ArticleGoogle Scholar
- Hückelhoven R. BAX Inhibitor-1, an ancient cell death suppressor in animals and plants with prokaryotic relatives. Apoptosis. 2004;9(3):299–307.PubMedView ArticleGoogle Scholar
- Eichmann R, et al. The barley apoptosis suppressor homologue BAX inhibitor-1 compromises nonhost penetration resistance of barley to the inappropriate pathogen Blumeria graminis f. sp. tritici. Mol Plant Microbe Interact. 2004;17(5):484–90.PubMedView ArticleGoogle Scholar
- Van Loon LC, Rep M, Pieterse CMJ. Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006;44:135–62.PubMedView ArticleGoogle Scholar
- Rakwal R, Agrawal GK, Jwa N-S. Characterization of a rice (Oryza sativa L.) Bowman–Birk proteinase inhibitor: tightly light regulated induction in response to cut, jasmonic acid, ethylene and protein phosphatase 2A inhibitors. Gene. 2001;263(1):189–98.PubMedView ArticleGoogle Scholar
- Thomma BPHJ, et al. Separate jasmonate-dependent and salicylate-dependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens. Proc Natl Acad Sci. 1998;95(25):15107–11.PubMedPubMed CentralView ArticleGoogle Scholar
- Govrin EM, Levine A. The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea. Curr Biol. 2000;10(13):751–7.PubMedView ArticleGoogle Scholar
- Metz RP, et al. Differential transcriptional regulation by mouse single-minded 2 s. J Biol Chem. 2006;281(16):10839–48.PubMedView ArticleGoogle Scholar
- Mortazavi A, et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.PubMedView ArticleGoogle Scholar
- Dillies M-A, et al. A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief Bioinform. 2013;14(6):671–83.PubMedView ArticleGoogle Scholar
- Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25.PubMedPubMed CentralView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Series B. 1995;57(1):289–300.Google Scholar
- Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics. 2002;18(1):207–8.PubMedView ArticleGoogle Scholar
- Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001;29(9):e45.PubMedPubMed CentralView ArticleGoogle Scholar
- Blüthgen N, et al. Biological Profiling of Gene Groups utilizing Gene Ontology. Genome Inform. 2005;16(1):106–15.PubMedGoogle Scholar
- Horton P, et al. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35 suppl 2:W585–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Punta M, et al. The Pfam protein families database. Nucleic acids research, 2011. doi:10.1093/nar/gkr1065.
- Saier MH, Tran CV, Barabote RD. TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006;42(34 suppl 1):D181–6.View ArticleGoogle Scholar
- Cantarel BL, et al. The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 2009;37 suppl 1:D233–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Yin Y, et al. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40(W1):W445–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A. MEROPS: the peptidase database. Nucleic Acids Res. 2010;38 suppl 1:D227–33.PubMedPubMed CentralView ArticleGoogle Scholar