- Research article
- Open Access
The wheat pathogen Zymoseptoria tritici senses and responds to different wavelengths of light
BMC Genomics volume 21, Article number: 513 (2020)
The ascomycete fungus Zymoseptoria tritici (synonyms: Mycosphaerella graminicola, Septoria tritici) is a major pathogen of wheat that causes the economically important foliar disease Septoria tritici blotch. Despite its importance as a pathogen, little is known about the reaction of this fungus to light. To test for light responses, cultures of Z. tritici were grown in vitro for 16-h days under white, blue or red light, and their transcriptomes were compared with each other and to those obtained from control cultures grown in darkness.
There were major differences in gene expression with over 3400 genes upregulated in one or more of the light conditions compared to dark, and from 1909 to 2573 genes specifically upregulated in the dark compared to the individual light treatments. Differences between light treatments were lower, ranging from only 79 differentially expressed genes in the red versus blue comparison to 585 between white light and red. Many of the differentially expressed genes had no functional annotations. For those that did, analysis of the Gene Ontology (GO) terms showed that those related to metabolism were enriched in all three light treatments, while those related to growth and communication were more prevalent in the dark. Interestingly, genes for effectors that have been shown previously to be involved in pathogenicity also were upregulated in one or more of the light treatments, suggesting a possible role of light for infection.
This analysis shows that Z. tritici can sense and respond to light with a huge effect on transcript abundance. High proportions of differentially expressed genes with no functional annotations illuminates the huge gap in our understanding of light responses in this fungus. Differential expression of genes for effectors indicates that light could be important for pathogenicity; unknown effectors may show a similar pattern of transcription. A better understanding of the effects of light on pathogenicity and other biological processes of Z. tritici could help to manage Septoria tritici blotch in the future.
Light is essential for many biological processes and is an important environmental cue. Fungi have multiple responses to light, which can vary from species to species. Much of the research done on light responses in ascomycetes has been with the model filamentous fungus Neurospora crassa [1,2,3,4,5,6,7,8,9]. This species has clearly delineated morphology between growth that occurs in the light versus the dark [1,2,3,4,5,6].
Multiple species in a variety of families rely on light to regulate the machinery that helps them cope with stresses. UV radiation damages DNA, and the cells respond by producing photolyase proteins to repair the damage. The expression of photolyases is often induced by light, as it is only required after exposure to UV [10,11,12,13,14]. Oxidative stress, which can be caused by light, relies on light-sensing genes to induce production of proteins required for cellular repair [15, 16]. Many fungi are protected from photobiological damage by pigments, such as melanin and various carotenoids, production of which can be induced by light [2, 16,17,18,19,20,21].
Light can have major effects on fungal growth and morphology. In some species of fungi, light induces asexual reproduction via the formation and germination of conidia, and represses sexual reproduction [22,23,24,25,26,27]. Yet the opposite pattern is found in other species [17, 28,29,30]. Growth in day:night cycles can lead to fungal colonies showing concentric circles of differing morphologies on agar plates [17, 31, 32]. Presumably this also would occur within plant hosts and may explain the spread of lesions in concentric circles seen in fungi such as many species in the genus Alternaria, some of which show maximum growth in culture under alternating cycles of light and dark .
In some fungal species, light regulates when toxins and other secondary metabolites are produced [30, 33,34,35,36,37,38]. In Cercospora species that produce the light-activated phytotoxin cercosporin, light induces production of the toxin, as that is when it is most effective [39, 40], and no cercosporin is produced in the dark. A similar phenomenon occurs in Alternaria alternata where production of the mycotoxins altertoxin and alternariol are induced exclusively or have expression increased greatly under blue light in contrast to dark . Similarly, the production of aflatoxin in Aspergillus species can be affected by both the color and intensity of light, although the conditions under which mycotoxins are produced at the highest rates are not uniform within the genus [35, 42]. Alternatively, in Fusarium graminearum, light represses the production of the trichothecene mycotoxins deoxynivalenol and 15-acetyl-deoxynivalenol .
Curiously, the ability to sense light also is required for pathogenicity and virulence in some pathogenic fungal species. For example, in the plant pathogens Botryis cinerea and Cercospora zea-maydis, knocking out wc1, the light-sensing component of the proteins comprising the White Collar Complex (WCC), leads to lowered virulence or a complete lack of pathogenicity [15, 44]. In Aspergillus flavus, the deletion of wc1 or velvet homologs drastically reduce the ability of the fungus to infiltrate corn kernels, peanuts, and cotton bolls [30, 37, 38]. However, in the rice pathogen Magnaporthe oryzae, light represses infection, and leads to much lower disease severity . A wc1 knockout in M. oryzae showed greater disease severity, as it could infect in the light as well as in the dark, compared to wild type which only infects in the dark .
Zymoseptoria tritici (synonyms: Mycosphaerella graminicola, Septoria tritici) is the causal agent of Septoria tritici blotch, an economically important disease of wheat. Losses due to this disease can reach up to 50% in epidemic years, and often vary between 5 and 20% depending on the environment and the cultivar of wheat; it has been estimated that up to 70% of fungicide use in Europe is to control this disease [46,47,48]. Spores of Z. tritici are splash dispersed during rainstorms, and need humid conditions for successful infection . After landing on a leaf, the spores germinate, and invade the wheat plants via the stomata . Initial growth appears to be biotrophic, but the fungus rapidly switches to necrotrophic growth beginning 12–14 days after penetration . Controlling Z. tritici is becoming more difficult, as resistance to the strobilurin (quinone-outside inhibitor) fungicides has become widespread in Europe, and also has appeared in North America [48, 52]. Other effective fungicides are the demethylation inhibitors (azoles), the SDHIs (succinate dehydrogenase inhibitors), and multi-site fungicides . Resistance to the azoles has already begun to spread, and the SDHI fungicides are at a medium to high risk of fungal populations developing resistance [48, 53,54,55,56]. Understanding more about how Z. tritici infects wheat and what conditions are necessary for it to reproduce are critical for developing better methods of disease control.
Currently, little is known about how Z. tritici senses and responds to light. Light does have minor effects on the growth and development of Z. tritici, stimulating formation of aerial hyphae, and prolonging the time growing in yeast-like form before transitioning to hyphal form [57, 58]. Light is a very important environmental cue for some fungi in the Dothideomycetes, the class that contains Z. tritici and many other important plant pathogens [59, 60]. For example, in the genus Cercospora, another genus in the same taxonomic order (Capnodiales), light represses asexual sporulation, and melanization is controlled by circadian rhythms [17, 44]. The disruption of CRP1, a homolog of wc1 in N. crassa, eliminates the stomatal tropism that C. zea-maydis needs to infect maize . The conidiation and growth of two species in the order Pleosporales, Alternaria alternata and Exserohilum turcicum, are also regulated by light [41, 61]. Additionally, as mentioned previously, light regulates toxin production in multiple Dothideomycetes species [33, 40, 41, 62].
The Z. tritici genome  contains genes coding for a number of homologs to photoreceptors in other fungi. These include homologs for the two White Collar Complex genes (wc1 and wc2), VIVID, a blue-light-sensing cryptochrome, a photolyase plus two putative photolyase genes and a red-light-sensing phytochrome [10, 25, 64, 65]. It also has homologs of some genes that may not be involved in sensing light, but do respond to light in other species, such as frq, a circadian rhythm gene that codes for the frequency clock protein, and velvet, which responds to light and has been shown to react with a phytochrome and the white collar complex in Aspergillus nidulans [66,67,68,69]. Knocking out function of the velvet gene, MVE1, in Z. tritici had large effects on growth and development, including increased sensitivity to stresses, reduced melanin production, and blindness to light-induced aerial hyphae formation .
Between the evidence that closely related fungi can sense light and the presence of putative photoreceptors and photoreceptor homologs, it seems highly likely that Z. tritici can sense and respond to light, but this has not been tested specifically. The goal of this research was to test the hypothesis that Z. tritici can sense and respond to light. Secondary goals were to identify genes that are highly regulated by light for future study and to augment the annotation of the reference genome by analyzing the expression of genes induced by different wavelengths of light.
RNA samples extracted from mycelia grown in three biological replicates of four light growth conditions, 16 h per day of white, blue (400–530 nm), or red (600–700 nm) light and continuous darkness, were sequenced and the poly A, unstranded reads were filtered for quality, mapped, and analyzed. On average, 1.35% of all sequence reads were too short after trimming and were removed from each of the samples, and an average of 98.7% of the reads mapped to the Z. tritici genome annotations (Table 1). The remaining 1.3% of the reads mostly appeared to be contaminants. When the read counts for each gene were plotted between the three biological replicates of each condition, the reproducibility was very high (Additional File 1A). The lowest variability appeared to occur between the three white light replicates, while the greatest was between the first and third dark replicates (Additional File 1A); however, even this was low compared to the variability among the different treatments.
Out of 13,522 total genes called in the reference genome, 12,079 (89%) had greater than 10 reads in each replicate and more than 100 reads across all twelve replicates, and were used in analysis with DESeq2. The cleaned data were distributed normally, both within individual replicates and as a whole (Additional File 1B).
Differential gene expression in response to light
There were major differences in the numbers of differentially expressed genes (adjusted p value < 0.05 and an absolute log2 fold change (LFC) > 2) between the dark condition and any of the three light treatments, with all three of the light conditions versus dark having over 2000 differentially expressed genes each, while white versus red was the only other comparison to break 100 differentially expressed genes (Table 2). The white light versus red light comparison was distinct from the white light versus blue light and red light versus blue light comparisons, as well as from the light vs dark comparisons (Table 2). All six comparisons between the treatments showed the expected concavity in the volcano plots, but the magnitudes of the adjusted p values were much greater in the three light versus dark comparisons, with the lowest magnitude of adjusted p values in white light versus blue light (Additional File 1C).
In total, 4187 unique genes were differentially expressed between all six comparisons. This is just under a third (31%) of all genes present in the Z. tritici genome. In the three light versus dark comparisons, 4019 unique genes were differentially expressed, which is very near the total unique genes, and still represents 30% of all genes in Z. tritici. In the three light versus light comparisons, 804 unique genes were differentially expressed, which is only about 6% of the genes in the Z. tritici genome (Additional File 2).
The dark condition had the highest number of unique differentially expressed genes in all three comparisons with the light treatments (Fig. 1). Blue light consistently had the fewest unique differentially expressed genes compared to all of the other treatments (Fig. 1). Many of the genes that are differentially expressed between any light condition and full darkness are differentially expressed between all light conditions and full darkness, with white light having the highest number of unique differentially expressed genes (Fig. 1).
Receptors play an important role in sensing and responding to light. Many of the potential light-sensing proteins in Z. tritici were not differentially expressed in any comparison (Fig. 2a). Multiple genes with photolyase-like regions were differentially expressed under various conditions. The cryptochrome/photolyase gene CRY was upregulated in white light and blue light, but not in red light or darkness. The rhodopsin-like gene NOP-1 was upregulated in all types of light versus dark. The known light-signaling gene MVE1 was upregulated in the dark relative to white light and blue light (Fig. 2b). Interestingly, four genes for effector proteins were differentially expressed under different light conditions, two LysM effectors, Avr3D1, and AvrStb6 (Fig. 2a).
There are three MAPK (mitogen-activated protein kinase)-encoding genes in Z. tritici [71,72,73]. MgHog1 (ZtritIPO323_04g02798) was strongly downregulated in light, especially in red light (Fig. 2b). Another MAPK-encoding gene, MgSlt2 (ZtritIPO323_04g00461), was upregulated in light (Fig. 2b). A third MAPK-encoding gene, MgFus3 (ZtritIPO323_04g10805), was not differentially expressed in any comparison (Fig. 2b). All three of these MAPK-encoding genes are essential for full virulence of Z. tritici as well as the production of melanin [71,72,73], a critical component of photoprotection in this fungus. Another melanization-related gene, PKS1, is downregulated in white light versus blue and red light, which indicated that it also may be regulated by a light color such as green, or by an interplay of the blue-light and red-light sensing genes (Fig. 2b) .
Three small, secreted proteins (SSP) are highly differentially expressed in various comparisons (Fig. 2a) . ZtritIPO323_04g09921 is strongly upregulated in blue and white light versus red light and dark and is downregulated in red light versus dark. This indicated that expression of this gene might be increased by white and blue light but suppressed by red light relative to dark. Another SSP, ZtritIPO323_04g11421, was highly upregulated in blue and white light versus red light and dark, indicating that its expression may be stimulated by the blue light spectrum. The third SSP, ZtritIPO323_04g03806, is strongly downregulated in all light conditions versus darkness, so may be repressed by all types of light. All three of these SSPs are functionally unannotated, and had no fungal matches in the NCBI NR database. One of the few SSPs with a functional annotation, ZtritIPO323_04g10536, has a non-orthologous group (NOG) of “chitin binding peritrophin-A” domains (10PPR@NOG according to the eggNOG annotation) , and is slightly downregulated in light (log2 fold change in white versus dark: − 2.16, blue versus dark: − 2.01, red versus dark: − 1.88).
One protease had a highly significant expression pattern (Fig. 2a) . MgAsp2 (ZtritIPO323_04g06056) was downregulated in all but one comparison, the exception being blue versus red light. The NOG associated with this gene is 03JP4@ascNOG, a secreted aspartic protease orthologous group, and blastp links it to other aspartic proteases . Many other proteases had significant differences in transcript abundance between treatments, with approximately a third of the proteases upregulated, a third downregulated, and a third not differentially expressed in the light versus dark comparisons (Additional File 3D).
No obvious patterns were apparent in differential expression of genes for secreted proteins from other functional classes. Very few secreted lipase genes were differentially expressed in the various light versus dark comparisons (Additional File 3B). Most secreted peroxidases were more differentially expressed in white light rather than blue or red (Additional File 3E). Almost half of the plant cell wall degrading enzymes (PCWDEs) were downregulated by light, but with no consistent pattern (Additional File 3C).
The genes related to secondary metabolism  tend to be upregulated by light if they are differentially expressed. Among all the genes, nearly half of them are upregulated, and only a quarter of them are downregulated in light versus dark comparisons (Fig. 2c). In each of the categories of secondary metabolism genes, most had more genes upregulated than downregulated or not differentially expressed (Fig. 2c). The exception is the polyketide synthase-like genes (PKS-L), where there was only one, and it was downregulated only in white light versus dark (Fig. 2c).
Gene ZtritIPO323_04g04132 showed very high differential expression (Fig. 2a). The adjusted p value for this gene in the white versus dark and blue versus dark comparisons was 0, with a log2 fold change of 8.7 and 8.9, respectively, and it had an adjusted p value of 1.1e− 214 and a log2 fold change of 6.7 in the red versus dark comparison. It was also less dramatically differentially expressed in comparisons between white light versus red and blue versus red. The only functional annotations for this gene are an alpha-beta hydrolase fold (03NR3@ascNOG) and a methyl ester carboxylesterase conserved domain according to the NCBI Conserved Domain Database; it was not annotated as being in a class related to pathogenicity by Palma-Guerrero et al. [75, 76, 78,79,80,81].
Among all differentially expressed genes, approximately 3% of those in the light treatments relative to dark were located on dispensable chromosomes that are not present in all isolates of the pathogen ; all of the rest occurred on the core chromosome set. Very few genes on dispensable chromosomes were differentially expressed when the three light treatments were compared to each other, ranging from 1 in the blue versus red light comparison to 57 between white light and red (Table 3).
From 48 to 63% of the genes with significant differential expression in each comparison had no corresponding functional annotation (Table 2, Additional File 2). This includes KOG classes and GO terms in the frozen gene catalog created as part of the original annotation by the Joint Genome Institute , as well as EggNOG classifications  generated with the analysis of additional RNA sequences by King . On average, 45.5% of genes could be assigned functional protein annotations, while 27.7% could be assigned GO terms. In the full genome, 55.3% of genes have annotations, and 32.8% have GO terms assigned. There is a distinct division between the three light/dark comparisons and the light/light comparisons; the light/light comparisons have a higher percentage of genes annotated to the full genome (average of 59.3%), while the light/dark comparisons have fewer genes annotated than the full genome (average of 49.7%).
The enriched GO terms were quite diverse between the comparisons. The numbers of GO terms enriched between the three light treatments and the dark treatment were higher than those between the light conditions (Table 4). This was similar to the raw numbers of differentially expressed genes (Table 2).
In pairwise comparisons between the light conditions, white light had fewer enriched GO terms than blue or red (Table 4). Compared to blue, white light had only one enriched GO term, protein phosphorylation (GO:0006468), that was also enriched in the dark in all three of the light versus dark comparisons. The gene this GO term was associated with was a serine/threonine kinase. More enriched GO terms were identified when cultures were grown under blue or red light than in the wider-range, white-light condition.
The dark condition had a wider range of enriched GO terms relative to the light conditions (Fig. 3). The variety of enriched GO terms was especially notable in the dark versus blue light comparison (Fig. 3). The dark treatment had more GO terms that were related to growth and development than did cultures grown in the light, where those GO terms occurred rarely. This was very apparent in the blue light and red light versus dark comparisons, where multiple GO terms for growth were enriched in the dark in each comparison. These GO terms included the high-level growth GO term (GO:0040007), as well as filamentous growth (GO:0030447), chromosome segregation (GO:0007059), cellular component biosynthesis (GO:0044085), and numerous cellular components related to mitotic division (GO:0000776, GO:0000794, GO:0000942, GO:0000793), among others.
The light conditions versus dark had many more enriched GO terms related to metabolism (Fig. 3 in red) as well as some that were related to transportation of substances in the cell. The dark versus white and dark versus blue comparisons were enriched in GO terms related to growth while the white versus dark comparison contained several GO terms related to communication.
The comparisons between the three light treatments had fewer enriched GO terms compared to the light versus dark comparisons. The exception to this was the genes that were downregulated in white versus red light, which looked more like a comparison with a dark treatment than did the white versus red and blue versus red comparisons (Fig. 3).
Many of the light comparisons have enriched GO terms for responses to oxidative stress, including those that are enriched on both sides of a comparison, such as the general response to oxidative stress (GO:0006979), which was enriched in both white and red light. White and blue light had more of these GO terms enriched than red light, and dark had none. These include response to oxidative stress (GO:0006979), a response to hydrogen peroxide (GO:0042542), peroxidase activity (GO:0004601), and multiple oxidoreductase-related GO terms (GO:0016684, GO:0016705, GO:0016634, GO:0006733, GO:0016634, GO:0016722).
The KOG functional annotations had similar results to the GO annotations, showing that more genes related to the transport and metabolism of major substances are upregulated in the light rather than the dark. The class “defense mechanisms” also had more genes upregulated in the light than in the dark. The opposite trend occurred with some KOG classes related to growth, which had more genes upregulated in the dark. These KOG classes include “cell cycle control, cell division, chromosome partitioning”, “cell wall/membrane/envelope biogenesis”, “signal transduction mechanisms”, “replication, recombination and repair”, and “intracellular trafficking, secretion, and vesicular transport”. These differences did not appear in the light versus light comparisons. Despite these changes in transcription of genes involved in growth and development and previous studies [57, 58], there were no obvious differences in the morphologies of cultures grown in the dark versus any of the light conditions (Fig. 4).
The RNAseq analyses show clearly that Z. tritici can sense and respond to light. The high numbers of differentially expressed genes between the comparisons indicate that Z. tritici alters its transcriptome in the dark compared to when it is growing in the light and between different wavelengths of light. While this is not surprising, as many other fungi have been shown to sense and respond to light, this is the first confirmation that the hemibiotrophic wheat pathogen Z. tritici also uses light to trigger large differences in gene transcription, which presumably has a correspondingly large effect on metabolic processes [15, 20, 34, 83, 84].
Approximately half of the differentially expressed genes (48–63% depending on the comparison) do not have any functional annotations so this is the first time they have been implicated in possible responses to light or dark. The proportions of unannotated genes among those that were differentially expressed in the light and dark treatments were moderately higher than their representations in the genome in general. This could be in part because no other analyses of differential expression in response to light have been performed on Z. tritici, but also could indicate that this fungus has different light responses than those that have been seen in other species. Many prior analyses of light responses in fungi have been with the model fungus N. crassa in the class Sordariomycetes. Light responses of fungi in the class Dothideomycetes to which Z. tritici belongs are less well understood, and it is possible that this species may use different genes or have different responses compared to other fungi. The differentially expressed genes with no functional annotations could indicate that novel pathways are involved in the responses of Dothideomycetes to light. This hypothesis could be tested in the future by knocking out those genes to identify any altered phenotypes in response to different light conditions.
Light affected expression of genes on all 21 chromosomes in the Z. tritici genome, including all eight dispensable (or accessory) chromosomes, indicating that light most likely has a significant effect on growth and development of this fungus. The percentages of genes that were differentially expressed in any of the comparisons agrees with what was found previously in N. crassa, indicating that while the response was large, this is likely correct . One interesting result was that genes with annotations relating to oxidative stress, metabolism and transportation were upregulated primarily in the light, while the fungus would be growing on its host during the day, while genes with annotations relating to growth and general cell maintenance were more common primarily in the dark, and therefore would be expressed at night. Further experiments are needed to test whether these in vitro experiments reflect gene expression on the host. Unlike the pattern seen in Z. tritici, in Trichoderma atroviride, transport-related genes were repressed in the wild type when grown under constant light and were expressed in a wc1 knockout mutant . However, our findings on cellular metabolism are similar to what was found in N. crassa previously [16, 86].
The upregulation of genes related to oxidative stress under light indicates another common way that fungi can respond to radiation-induced stresses. Oxidoreductases are involved in reducing free radicals generated by UV radiation so their increased expression in response to light would be adaptive for the fungus. Light has been shown to increase a response to oxidative stresses in other organisms [15, 16, 20, 85] and seems to induce similar responses in Z. tritici.
Expression patterns for many genes were as expected based on their predicted biology. For example, cryptochromes are involved in sensing and responding to blue light in other organisms [10, 87]. The primary cryptochrome gene in Z. tritici, Cry, is upregulated very specifically in white and blue light as would be expected from its biological function. Photolyases repair DNA damage caused by exposure to light, particularly UV . One of the two putative Z. tritici photolyase genes analyzed was strongly upregulated in response to blue and white light, which would be the closest to UV, while the other was upregulated at a low level. The main photolyase gene, Phr, was not differentially expressed in response to light so may have a different function or could be expressed constitutively to protect against DNA damage in general.
In contrast, many genes with light-related biological functions in other organisms were not differentially expressed in the Z. tritici experiments. This included light-sensing genes such as wc1 and the phytochrome genes (Phy-1/2). While they may not be differentially expressed under long exposure to light, such as that used here, it is possible that other results would be obtained during transitory periods of light, such as during dawn or dusk, which affects wc-1 in N. crassa . Many of these genes may be regulated in a similar fashion, with higher differential expression when light conditions are changing. Another possibility is that the major light-sensing proteins are expressed constitutively so that they are available to detect changes when they occur. Additional experiments covering diurnal variations are required to answer these questions.
Proteins likely to be involved in fungal virulence showed some interesting expression patterns. The most obvious of these were the increased transcription of LysM domain-containing genes in response to light. These code for effector proteins that can minimize wheat host defenses in response to Z. tritici during the initial invasion phase of infection by binding chitin, and protecting the cell wall from hydrolysis . Mg1LysM was upregulated by any light, while Mg3LysM was upregulated primarily by red light. Mg3LysM is the effector primarily responsible for blocking the wheat host from activating defenses , so it is interesting that it would be upregulated in red light, rather than be expressed constitutively. Similarly, while Mg1LysM is not actively responsible for blocking a defense response, it is not clear why it would be expressed in the light rather than constitutively unless it has another role related to growth or pathogenicity.
Some other genes that could be linked to pathogenicity, including proteases and SSPs, also showed differential expression in response to light. Very few of the differentially expressed SSPs had any functional annotations other than being flagged as such by Palma-Guerrero et al. . One of those with an annotation that was downregulated in response to light, ZtritIPO323_04g10536, had a chitin-binding peritrophin-A domain. These domains are found primarily in insects , and most of the non-orthologous group annotations were from Drosophila species. Preventing the plant from sensing chitin is critical for avoiding PAMP-triggered immunity and other responses against fungal invaders, and is the major effect of the LysM effector proteins secreted by Z. tritici during infection. Curiously, the expression patterns of ZtritIPO323_04g10536 were opposite those of the LysM genes, which were induced by light. This could be an indication that this is another chitin-binding effector that works opposite to the LysM effectors.
The proteases were fairly evenly divided between being upregulated by light, downregulated by light, or not differentially expressed. Expression of proteases is consistent with the conclusion of Goodwin et al.  that the genome of Z. tritici was more similar to those of endophytes rather than other pathogens and that pathogenicity might involve catalysis of proteins rather than carbohydrates during the early stages of infection. However, why most of those genes would have higher expression in response to light is not known. In contrast to most proteases, the one with the highest differential expression, aspartic protease MgAsp2, is part of a family of genes that is potentially involved in pathogenicity in B. cinerea , yet had much higher expression in dark than in light.
Pigmentation is a critical component of photoprotection, or is controlled by circadian rhythms to produce pigments during the day [15, 17, 20, 44, 91, 92]. Despite this tendency, the pigment-related genes in Z. tritici varied in whether they were up or down regulated in the light. The three MAPK-encoding genes provide a good example. All three are required for melanization [71,72,73], yet each has a different expression pattern. MgFus3 was not differentially expressed, MgHog1 was downregulated by light, and MgSlt2 was upregulated by white and blue light (Fig. 2b). These genes are also required for full virulence, via different parts of the infection cycle [71,72,73]. Another gene that is known to be directly involved in the synthesis of melanin, the polyketide synthase PKS1 , was downregulated in white light only. It was not differentially expressed in any other comparison, indicating that it may not be photoregulated or is regulated by a color of light outside of blue and red. We found that melanization of Z. tritici hyphae occurs independent of light conditions, and this further supports the hypothesis that melanization is not regulated by light in this fungus.
Other genes in the secondary metabolism pathways also are differentially expressed and tend to be upregulated in the light. Nearly half of them were upregulated in the light, and the remaining half were divided between being downregulated and not differentially expressed (Fig. 2c). This correlates with the findings from the enriched GO terms, where metabolism and secondary metabolism GO classes were expressed more in the light than in the dark (Fig. 3).
There are many other classes of genes that have been linked to pathogenicity in other species of fungal pathogens, including the PCWDEs, peroxidases, and lipases. The effect of light on these genes varied within classes. There were some general trends, such as nearly half of the PCWDEs were downregulated in light, but due to a lack of detailed annotations, more granular regulation trends are difficult to ascertain. A better understanding of how these genes are expressed in the host under various light conditions is essential for a complete picture of gene expression during infection.
This experiment provided a single snapshot of gene expression under different light conditions, rather than analyzing changes over time or during the transitions from dark to light and vice versa. As such, some genes that may be differentially expressed only during a short period after light exposure would not be found. Previous RNAseq experiments in N. crassa show that there are genes that are differentially expressed, both up and down regulated, during the 15 to 240 min after initial exposure to light that then return to basal dark expression . It is highly likely that additional Z. tritici genes would have been differentially expressed if samples had been collected during the first few minutes or hours after the transitions to light or dark. For instance, white collar 1 in N. crassa was highly differentially expressed during the first 15 min following exposure to white light, but returned to dark levels of transcription after an hour .
It was not possible to conclude which color of light stimulates the greatest response from Z. tritici. The expression patterns under red and blue light were very similar, with only a small percentage of genes differing between these two light-versus-dark comparisons. The white light versus dark comparison had the highest number of differentially expressed genes among the three light treatments, but since blue and red are both components of white light this does not help to identify which color can be sensed most efficiently. The genome sequence of Z. tritici  appears to have more genes for sensing and responding to blue light, including the WCC, VIVID, cryptochrome, and photolyase genes compared to only a single phytochrome gene for red. Biologically it may make sense for the fungus to have a heightened ability to sense and respond to blue light, as this would most likely co-occur with damaging UV, while red light alone would be present in nature mostly when the light is changing near sunrise or sunset and the amount of UV would be lower.
These limited initial experiments on gene expression show that Z. tritici can sense and respond to light, with profound effects on growth, development and metabolism. The large number of functionally unannotated genes involved in light responses indicates a huge gap in our knowledge that must be filled before we can fully understand the effects of light on this fungus. We still know nothing about the minimum intensity of light that can trigger a response, how long light must be present before a response is initiated, or whether the fungus can specifically sense UV. Differential expression of various genes potentially or proven to be involved with infection indicates that light could be important for pathogenicity, and other effectors may show a similar pattern of transcription. A better understanding of the effects of light on pathogenicity and other biological processes of this fungus could provide the basis for development of improved disease management strategies in the future.
Fungal growth and light treatments
The Zymoseptoria tritici isolate IPO323 was used for these analyses, as it has been sequenced completely, grows readily in culture and, as an excellent reference genome, has been the subject of much prior research . This isolate was grown on sealed yeast-sucrose agar plates (YSA; 10.0 g/L yeast extract, 10.0 g/L sucrose, 15.0 g/L BactoAgar) under filters to control the wavelengths of light received by each culture. The treatments were white light (no filter), blue (400–530 nm) and red (600–700 nm) light using colored acrylic glass filters (American Acrylics, Skokie, IL) . Full dark was achieved by wrapping the plates in aluminum foil. Light for the cultures was provided by 32-W fluorescent bulbs. Each condition had three biological replicates consisting of one 9-cm Petri dish each. The three light conditions had 16:8 day:night cycles. The cultures were initiated with plugs taken from the margins of an actively growing mycelium and grown for a week at 23 °C.
RNA extraction and sequencing
Fungal tissue was quickly scraped from the agar, frozen immediately in liquid nitrogen for grinding, and total RNA from the cultures was extracted using a QIAGEN RNeasy Plant Mini Kit (Catalog Number 74903), following the manufacturer’s recommendations. The RNA was then sent to the Purdue Genomics Core Facility (West Lafayette, IN) to be processed and sequenced on an Illumina HiSeq 2500. The data were downloaded to the Halstead Computing Cluster of the Purdue Rosen Center for Advanced Computing for analysis.
RNAseq and statistical analyses
Trimmomatic (version 0.36) was used to remove leading and trailing bases with a phred score lower than 30, and reads that were shorter than 40 base pairs long after base removal . HISAT2 (version 2.0.5) was used to map the remaining reads to the Z. tritici reference genome [94,95,96]. Samtools (version 1.4) was used to sort the mapped reads, and HTSeq-count (version 0.7.0) was used with the King Rothamsted reference annotation to obtain a counts table for genes [82, 97, 98].
The creation of the reads library was done by using the paired-end reads from each of three biological replicates for each of the four light conditions. The table of gene read counts was exported to R (version 3.4.0). The gene read counts were cleaned by removing any genes that had a rowsum of fewer than 100 reads totaled over all samples, and any genes where any one replicate had fewer than 10 reads, in that order. This conservative approach was used to limit the analysis to genes with reliable data for all replicates. DESeq2 (version 1.15.51) was used to calculate differential gene expression and to obtain adjusted p values and log2 fold changes. UpSetR (version 1.4.0) was used to generate the UpSet plot .
A custom Perl script was used to pull gene sequences from the reference annotation and the gene sequences were analyzed by eggnog-mapper (version 4.5.1) and InterProScan (version 5.24–63.0) [76, 99]. Another custom Perl script was used to collect each GO term assigned to each gene in a gene-to-GO table, which was fed in to TopGO (version 2.28.0) . These Perl scripts are available on request or from GitHub. Lists of differentially expressed genes from DESeq2 were used to determine overexpressed GO terms found in each comparison.
GOOSE (GO Online SQL Environment) was used to query the Gene Ontology (GO) database to generate levels [101, 102]. The maximum tree depth was used to organize the enriched GO terms for visual analysis. CateGOrizer was used to map the biological processes of each enriched GO term to a parent term in the standard GO_slim . ggplot2 was used to generate stacked bar charts from the CateGOrizer data .
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the SRA repository, under accession number SRP151591.
Custom Perl scripts are available from the authors upon request, or from GitHub under project name Z.triticiLightRNAseq at home page https://github.com/CBMcCorison/Z.triticiLightRNAseq.
White Collar Complex
Succinate dehydrogenase inhibitors
Zymoseptoria tritici Velvet gene
Mitogen-activated protein kinase
Small secreted protein
Plant cell wall degrading enzymes
Nonribosomal peptide synthetase
Sargent ML, Briggs WR. The effects of light on a circadian rhythm of conidiation in Neurospora. Plant Physiol. 1967;42:1504–10. https://doi.org/10.1104/pp.42.11.1504.
Harding RW, Turner RV. Photoregulation of the carotenoid biosynthetic pathway in albino and White collar mutants of Neurospora crassa. Plant Physiol. 1981;68:745–9. https://doi.org/10.1104/pp.68.3.745.
Degli-Innocenti F, Russo VEA. Isolation of new White collar mutants of Neurospora crassa and studies on their behavior in the blue light-induced formation of protoperithecia. J Bacteriol. 1984;159:757–61 http://jb.asm.org/content/159/2/757.full.pdf. Accessed 7 Jun 2018.
Nelson MA, Morelli G, Carattoli A, Romano N, Macino G. Molecular cloning of a Neurospora crassa carotenoid biosynthetic gene (albino-3) regulated by blue light and the products of the white collar genes. Mol Cell Biol. 1989;9:1271–6. https://doi.org/10.1128/mcb.9.3.1271.
Schmidhauser TJ, Lauter FR, Russo VE, Yanofsky C. Cloning, sequence, and photoregulation of al-1, a carotenoid biosynthetic gene of Neurospora crassa. Mol Cell Biol. 1990;10:5064–70. https://doi.org/10.1128/mcb.10.10.5064.
Lauter FR, Russo VEA. Blue light induction of conidiation-specific genes in Neurospora crassa. Nucleic Acids Res. 1991;19:6883–6. https://doi.org/10.1093/nar/19.24.6883.
Arpaia G, Loros JJ, Dunlap JC, Morelli G, Macino G. Light induction of the clock-controlled gene ccg-1 is not transduced through the circadian clock in Neurospora crassa. Mol Gen Genet. 1995;247:157–63. https://doi.org/10.1007/BF00705645.
Dasgupta A, Chen CH, Lee CH, Gladfelter AS, Dunlap JC, Loros JJ. Biological significance of photoreceptor Photocycle length: VIVID Photocycle governs the dynamic VIVID-White collar complex Pool mediating photo-adaptation and response to changes in light intensity. PLoS Genet. 2015;11:e1005215. https://doi.org/10.1371/journal.pgen.1005215.
Wang B, Zhou X, Loros JJ, Dunlap JC. Alternative use of DNA binding domains by the Neurospora White collar complex dictates circadian regulation and light responses. Mol Cell Biol. 2016;36:781–93. https://doi.org/10.1128/mcb.00841-15.
Bluhm BH, Dunkle LD. PHL1 of Cercospora zeae-maydis encodes a member of the photolyase/cryptochrome family involved in UV protection and fungal development. Fungal Genet Biol. 2008;45:1364–72. https://doi.org/10.1016/j.fgb.2008.07.005.
Berrocal-Tito G, Sametz-Baron L, Eichenberg K, Horwitz BA, Herrera-Estrella A. Rapid blue light regulation of a Trichoderma harzianum photolyase gene. J Biol Chem. 1999;274:14288–94. https://doi.org/10.1074/jbc.274.20.14288.
Veluchamy S, Rollins JA. A CRY-DASH-type photolyase/cryptochrome from Sclerotinia sclerotiorum mediates minor UV-A-specific effects on development. Fungal Genet Biol. 2008;45:1265–76.
Brancini GTP, Bachmann L, Eliana M, Ferreira S, Rangel DEN, Braga GÚL. Exposing Metarhizium acridum mycelium to visible light up-regulates a photolyase gene and increases photoreactivating ability. J Invertebr Pathol. 2018;152:35–7. https://doi.org/10.1016/j.jip.2018.01.007.
Alejandre-Durán E, Roldán-Arjona T, Ariza RR, Ruiz-Rubio M. The photolyase gene from the plant pathogen Fusarium oxysporum f. sp. lycopersici is induced by visible light and α-tomatine from tomato plant. Fungal Genet Biol. 2003;40:159–65.
Canessa P, Schumacher J, Hevia MA, Tudzynski P, Larrondo LF. Assessing the effects of light on differentiation and virulence of the plant pathogen Botrytis cinerea: characterization of the white collar complex. PLoS One. 2013;8.
Wu C, Yang F, Smith KM, Peterson M, Dekhang R, Zhang Y, et al. Genome-wide characterization of light-regulated genes in Neurospora crassa. G3-Genes Genom Genet. 2014;4:1731–45.
Bluhm BH, Burnham AM, Dunkle LD. A circadian rhythm regulating hyphal melanization in Cercospora kikuchii. Mycologia. 2010;102:1221–8. https://doi.org/10.3852/09-041.
Ruiz-Roldán MC, Garre V, Guarro J, Mariné M, Roncero MIG. Role of the White collar 1 photoreceptor in carotenogenesis, UV resistance, hydrophobicity, and virulence of Fusarium oxysporum. Eukaryot Cell. 2008;7:1227–30.
Krobanan K, Liang S-W, Chiu H-C, Shen W-C. The blue-light photoreceptor Sfwc-1 gene regulates the phototropic response and fruiting-body development in the homothallic ascomycete Sordaria fimicola. Appl Environ Microbiol. 2019;85. https://doi.org/10.1128/AEM.02206-18.
Fuller KK, Ringelberg CS, Loros JJ, Dunlap JC. The fungal pathogen Aspergillus fumigatus regulates growth, metabolism, and stress resistance in response to light. MBio. 2013;4:4–6.
Tang Y, Zhu P, Lu Z, Qu Y, Huang L, Zheng N, et al. The photoreceptor components FaWC1 and FaWC2 of Fusarium asiaticum cooperatively regulate light responses but play independent roles in virulence expression. Microorganisms. 2020;8:365.
Mooney JL, Yager LN. Light is required for conidiation in Aspergillus nidulans. Genes Dev. 1990;4:1473–82. https://doi.org/10.1101/gad.4.9.1473.
Brandhoff B, Simon A, Dornieden A, Schumacher J. Regulation of conidiation in Botrytis cinerea involves the light-responsive transcriptional regulators BcLTF3 and BcREG1. Curr Genet. 2017;63:931–49. https://doi.org/10.1007/s00294-017-0692-9.
Suzuki T, Nishimura S, Yagi K, Nakamura R, Takikawa Y, Matsuda Y, et al. Effects of light quality on conidiophore formation of the melon powdery mildew pathogen Podosphaera xanthii. Phytoparasitica. 2017. https://doi.org/10.1007/s12600-017-0631-9.
Blumenstein A, Vienken K, Tasler R, Purschwitz J, Veith D, Frankenberg-Dinkel N, et al. The Aspergillus nidulans phytochrome FphA represses sexual development in red light. Curr Biol. 2005;15:1833–8.
Wang Z, Li N, Li J, Dunlap JC, Trail F, Townsend JP. The fast-evolving phy-2 gene modulates sexual development in response to light in the model fungus Neurospora crassa. MBio. 2016;7. https://doi.org/10.1128/mBio.02148-15.
Wang Z, Wang J, Li N, Li J, Trail F, Dunlap JC, et al. Light sensing by opsins and fungal ecology: NOP-1 modulates entry into sexual reproduction in response to environmental cues. Mol Ecol. 2017;27:216–32. https://doi.org/10.1111/mec.14425.
Arraiano LS, Brown JKM. Sources of resistance and susceptibility to Septoria tritici blotch of wheat. Mol Plant Pathol. 2016;18:276–92. https://doi.org/10.1111/mpp.12482.
Arazoe T, Miyoshi K, Yamato T, Ogawa T, Ohsato S, Arie T, et al. Tailor-made CRISPR/Cas system for highly efficient targeted gene replacement in the rice blast fungus. Biotechnol Bioeng. 2015;112:2543–9.
Amaike S, Keller NP. Distinct roles for VeA and LaeA in development and pathogenesis of Aspergillus flavus. Eukaryot Cell. 2009;8:1051–60. https://doi.org/10.1128/EC.00088-09.
Das J, Busse H-G. Light-driven diurnal zonation in the filamentous fungus Fusarium solani. Int J Dev Biol. 1990;34:319–22 http://www.ijdb.ehu.es/web/paper.php?doi=2386731.
Hubballi M, Nakkeeran S, Raguchander T, Anand T, Samiyappan R. Effect of environmental conditions on growth of Alternaria alternata causing leaf blight of noni. World J Agric Sci. 2010;6:171–7.
Fanelli F, Reveglia P, Masi M, Mulè G, Zonno MC, Cimmino A, et al. Influence of light on the biosynthesis of ophiobolin a by Bipolaris maydis. Nat Prod Res. 2017;31:909–17. https://doi.org/10.1080/14786419.2016.1253084.
Fanelli F, Schmidt-Heydt M, Haidukowski M, Geisen R, Logrieco A, Mulè G. Influence of light on growth, fumonisin biosynthesis and FUM1 gene expression by Fusarium proliferatum. Int J Food Microbiol. 2012;153:148–53.
Suzuki T. Light-irradiation wavelength and intensity changes influence aflatoxin synthesis in fungi. Toxins (Basel). 2018;10:1–12.
Kim H, Ridenour JB, Dunkle LD, Bluhm BH. Regulation of pathogenesis by light in Cercospora zeae-maydis: an updated perspective. Plant Pathol J. 2011;27:103–9. https://doi.org/10.5423/PPJ.2011.27.2.103.
Duran RM, Cary JW, Calvo AM. The role of veA in Aspergillus flavus infection of peanut, corn and cotton. Open Mycol J. 2009;3:27–36. https://doi.org/10.1016/j.nedt.2013.03.020.
Kale SP, Milde L, Trapp MK, Frisvad JC, Keller NP, Bok JW. Requirement of LaeA for secondary metabolism and sclerotial production in Aspergillus flavus. Fungal Genet Biol. 2008;45:1422–9. https://doi.org/10.1016/j.fgb.2008.06.009.
Daub ME, Ehrenshaft M. The photoactivated Cercospora toxin cercosporin: contributions to plant disease and fundamental biology. Annu Rev Phytopathol. 2000;38:461–90. https://doi.org/10.1146/annurev.phyto.38.1.461.
You B-J, Lee M-H, Chung K-R. Production of cercosporin toxin by the phytopathogenic Cercospora fungi is affected by diverse environmental signals. Can J Microbiol. 2008;54:259–69. https://doi.org/10.1139/W08-002.
Pruß S, Fetzner R, Seither K, Herr A, Pfeiffer E, Metzler M, et al. Role of the Alternaria alternata blue-light receptor LreA (White-collar 1) in spore formation and secondary metabolism. Appl Environ Microbiol. 2014;80:2582–91. https://doi.org/10.1128/AEM.00327-14.
García-Cela E, Marin S, Sanchis V, Crespo-Sempere A, Ramos AJ. Effect of ultraviolet radiation a and B on growth and mycotoxin production by Aspergillus carbonarius and Aspergillus parasiticus in grape and pistachio media. Fungal Biol. 2015;119:67–78.
Kim H, Son H, Lee Y-W. Effects of light on secondary metabolism and fungal development of Fusarium graminearum. J Appl Microbiol. 2014;116:380–9. https://doi.org/10.1111/jam.12381.
Kim H, Ridenour JB, Dunkle LD, Bluhm BH. Regulation of stomatal tropism and infection by light in Cercospora zeae-maydis: evidence for coordinated host/pathogen responses to photoperiod? PLoS Pathog. 2011;7.
Kim S, Singh P, Park J, Park S, Friedman A, Zheng T, et al. Genetic and molecular characterization of a blue light photoreceptor MGWC-1 in Magnaporthe oryzae. Fungal Genet Biol. 2011;48:400–7. https://doi.org/10.1016/j.fgb.2011.01.004.
Fones H, Gurr S. The impact of Septoria tritici blotch disease on wheat: an EU perspective. Fungal Genet Biol. 2015;79:3–7. https://doi.org/10.1016/j.fgb.2015.04.004.
Eyal Z, Scharen AL, Prescott JM, van Ginkel M. The Septoria diseases of wheat: concepts and methods of disease management. Mexico: CIMMYT; 1987.
Torriani SFF, Melichar JPE, Mills C, Pain N, Sierotzki H, Courbot M. Zymoseptoria tritici: a major threat to wheat production, integrated approaches to control. Fungal Genet Biol. 2015;79:8–12.
Shaw MW, Royle DJ. Factors determining the severity of epidemics of Mycosphaerella graminicola (Septoria tritici) on winter wheat in the UK. Plant Pathol. 1993;42:882–99. https://doi.org/10.1111/j.1365-3059.1993.tb02674.x.
Kema GHJ, Yu D, Rijkenberg FHJ, Shaw MW, Baayen RP. Histology of pathogenesis of Mycosphaerella graminicola in wheat. Phytopathology. 1996;7:777–86.
Ponomarenko A, Goodwin SB, Kema GHJ. Septoria tritici blotch (STB) of wheat. Plant Health Instr. 2011;2:1–7.
Cheval P, Siah A, Bomble M, Popper AD, Reignault P, Halama P. Evolution of QoI resistance of the wheat pathogen Zymoseptoria tritici in northern France. Crop Prot. 2017;92:131–3. https://doi.org/10.1016/j.cropro.2016.10.017.
Wieczorek TM, Berg G, Semaškienė R, Mehl A, Sierotzki H, Stammler G, et al. Impact of DMI and SDHI fungicides on disease control and CYP51 mutations in populations of Zymoseptoria tritici from northern Europe. Eur J Plant Pathol. 2015;143:861–71.
Estep LK, Torriani SFF, Zala M, Anderson NP, Flowers MD, McDonald BA, et al. Emergence and early evolution of fungicide resistance in North American populations of Zymoseptoria tritici. Plant Pathol. 2015;64:961–71.
Fraaije BA, Bayon C, Atkins S, Cools HJ, Lucas JA, Fraaije MW. Risk assessment studies on succinate dehydrogenase inhibitors, the new weapons in the battle to control Septoria leaf blotch in wheat. Mol Plant Pathol. 2012;13:263–75.
Rehfus A, Strobel D, Bryson R, Stammler G. Mutations in sdh genes in field isolates of Zymoseptoria tritici and impact on the sensitivity to various succinate dehydrogenase inhibitors. Plant Pathol. 2018;67:175–80.
Choi Y-E, Goodwin SB. MVE1, encoding the velvet gene product homolog in Mycosphaerella graminicola, is associated with aerial mycelium formation, melanin biosynthesis, hyphal swelling, and light signaling. Appl Environ Microbiol. 2011;77:942–53.
Tiley AMM, White HJ, Foster GD, Bailey AM. The ZtvelB gene is required for vegetative growth and sporulation in the wheat pathogen Zymoseptoria tritici. Front Microbiol. 2019;10:2210. https://doi.org/10.3389/fmicb.2019.02210.
Quaedvlieg W, Kema GHJ, Groenewald JZ, Verkley GJM, Seifbarghi S, Razavi M, et al. Zymoseptoria gen. Nov.: a new genus to accommodate Septoria-like species occurring on graminicolous hosts. Persoonia. 2011;26:57–69.
Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, et al. Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 2012;8.
Flaherty JE, Dunkle LD. Identification and expression analysis of regulatory genes induced during conidiation in Exserohilum turcicum. Fungal Genet Biol. 2005;42:471–81. https://doi.org/10.1016/j.fgb.2005.02.001.
Daub ME, Herrero S, Chung KR. Photoactivated perylenequinone toxins in fungal pathogenesis of plants. FEMS Microbiol Lett. 2005;252:197–206.
Goodwin SB, Ben M’BS, Dhillon B, AHJ W, Crane CF, Hane JK, et al. Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011;7.
Smith KM, Sancar G, Dekhang R, Sullivan CM, Li S, Tag AG, et al. Transcription factors in light and circadian clock signaling networks revealed by genomewide mapping of direct targets for Neurospora white collar complex. Eukaryot Cell. 2010;9:1549–56. https://doi.org/10.1128/EC.00154-10.
Gin E, Diernfellner A, Brunner M, Höfer T. The Neurospora photoreceptor VIVID exerts negative and positive control on light sensing to achieve adaptation. Mol Syst Biol. 2013;9:1–11. https://doi.org/10.1038/msb.2013.24.
Bayram Ö, Braus GH. Coordination of secondary metabolism and development in fungi: the velvet family of regulatory proteins. FEMS Microbiol Rev. 2012;36:1–24. https://doi.org/10.1111/j.1574-6976.2011.00285.x.
Froehlich AC, Liu Y, Loros JJ, Dunlap JC. White Collar-1, a circadian blue light photoreceptor, binding to the frequency promoter. Science. 2002;297:815–9.
Lombardi LM, Brody S. Circadian rhythms in Neurospora crassa: clock gene homologues in fungi. Fungal Genet Biol. 2005;42:887–92.
Hunt SM, Thompson S, Elvin M, Heintzen C. VIVID interacts with the WHITE COLLAR complex and FREQUENCY-interacting RNA helicase to alter light and clock responses in Neurospora. Proc Natl Acad Sci U S A. 2010;107:16709–14. https://doi.org/10.1073/pnas.1009474107.
Conway JR, Lex A, Gehlenborg N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics. 2017;33:2938–40. https://doi.org/10.1093/bioinformatics/btx364.
Cousin A, Mehrabi R, Guilleroux M, Dufresne M, Van Der Lee T, Waalwijk C, et al. The MAP kinase-encoding gene MgFus3 of the non-appressorium phytopathogen Mycosphaerella graminicola is required for penetration and in vitro pycnidia formation. Mol Plant Pathol. 2006;7:269–78. https://doi.org/10.1111/j.1364-3703.2006.00337.x.
Mehrabi R, Zwiers L-H, de Waard MA, Kema GHJ. MgHog1 regulates dimorphism and pathogenicity in the fungal wheat pathogen Mycosphaerella graminicola. Mol Plant-Microbe Interact. 2006;19:1262–9. https://doi.org/10.1094/mpmi-19-1262.
Mehrabi R, van der Lee T, Waalwijk C, Kema GHJ. MgSlt2, a cellular integrity MAP kinase gene of the fungal wheat pathogen Mycosphaerella graminicola, is dispensable for penetration but essential for invasive growth. Mol Plant-Microbe Interact. 2006;19:389–98. https://doi.org/10.1094/mpmi-19-0389.
Lendenmann MH, Croll D, Stewart EL, McDonald BA. Quantitative trait locus mapping of melanization in the plant pathogenic fungus Zymoseptoria tritici. G3 genes, genomes. Genet. 2014;4:2519–33. https://doi.org/10.1534/g3.114.015289.
Palma-Guerrero J, Ma X, Torriani SFF, Zala M, Francisco CS, Hartmann FE, et al. Comparative transcriptome analyses in Zymoseptoria tritici reveal significant differences in gene expression among strains during plant infection. Mol Plant-Microbe Interact. 2017;30:231–44. https://doi.org/10.1094/mpmi-07-16-0146-r.
Huerta-Cepas J, Szklarczyk D, Forslund K, Cook H, Heller D, Walter MC, et al. EGGNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 2016;44:D286–93. https://doi.org/10.1093/nar/gkv1248.
Cairns T, Meyer V. In silico prediction and characterization of secondary metabolite biosynthetic gene clusters in the wheat pathogen Zymoseptoria tritici. BMC Genomics. 2017;18. https://doi.org/10.1186/s12864-017-3969-y.
Marchler-Bauer A, Bryant SH. CD-Search: Protein domain annotations on the fly. Nucleic Acids Res. 2004;32 WEB SERVER ISS:W327–31. https://doi.org/10.1093/nar/gkh454.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011;39(SUPPL. 1):D225–9. https://doi.org/10.1093/nar/gkq1189.
Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015;43:D222–6. https://doi.org/10.1093/nar/gku1221.
Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45:D200–3. https://doi.org/10.1093/nar/gkw1129.
King RC. Annotation of Zymoseptoria tritici merging JGI and custom annotations at Rothamsted. https://doi.org/10.6084/m9.figshare.4753708.v1. Accessed 15 Apr 2016.
Lee K, Singh P, Chung W, Ash J, Kim TS, Hang L, et al. Light regulation of asexual development in the rice blast fungus, Magnaporthe oryzae. Fungal Genet Biol. 2006;43:694–706. https://doi.org/10.1016/j.fgb.2006.04.005.
Ballario P, Vittorioso P, Magrelli A, Talora C, Cabibbo A, Macino G. White collar-1, a central regulator of blue light responses in Neurospora, is a zinc finger protein. EMBO J. 1996;15:1650–7.
Cetz-Chel JE, Balcázar-López E, Esquivel-Naranjo EU, Herrera-Estrella A. The Trichoderma atroviride putative transcription factor Blu7 controls light responsiveness and tolerance. BMC Genomics. 2016;17. https://doi.org/10.1186/s12864-016-2639-9.
Chen C-H, Ringelberg CS, Gross RH, Dunlap JC, Loros JJ. Genome-wide analysis of light-inducible responses reveals hierarchical light signalling in Neurospora. EMBO J. 2009;28:1029–42. https://doi.org/10.1038/emboj.2009.54.
Castrillo M, García-Martínez J, Avalos J. Light-dependent functions of the Fusarium fujikuroi CryD DASH cryptochrome in development and secondary metabolism. Appl Environ Microbiol. 2013;79:2777–88.
Sancar A. Structure and function of DNA photolyase. Biochemistry. 1994;33:2–9. https://doi.org/10.1021/bi00167a001.
Marshall R, Kombrink A, Motteram J, Loza-Reyes E, Lucas J, Hammond-Kosack KE, 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:756–69. https://doi.org/10.1104/pp.111.176347.
ten Have A, Dekkers E, Kay J, Phylip LH, van Kan JAL. An aspartic proteinase gene family in the filamentous fungus Botrytis cinerea contains members with novel features. Microbiology. 2004;150:2475–89. https://doi.org/10.1099/mic.0.27058-0.
Avalos J, Estrada AF. Regulation by light in Fusarium. Fungal Genet Biol. 2010;47:930–8. https://doi.org/10.1016/j.fgb.2010.05.001.
De Fabo EC, Harding RW, Shropshire W. Action spectrum between 260 and 800 nanometers for the photoinduction of carotenoid biosynthesis in Neurospora crassa. Plant Physiol. 1976;57:440–5. https://doi.org/10.1104/pp.57.3.440.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Sirén J, Välimäki N, Mäkinen V. Indexing graphs for path queries with applications in genome research. IEEE/ACM Trans Comput Biol Bioinforma. 2014;11:375–88. https://doi.org/10.1109/TCBB.2013.2297101.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60. https://doi.org/10.1038/nmeth.3317.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9. https://doi.org/10.1038/nmeth.1923.
Anders S, Pyl PT, Huber W. HTSeq-A Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9. https://doi.org/10.1093/bioinformatics/btu638.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9. https://doi.org/10.1093/bioinformatics/btp352.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40. https://doi.org/10.1093/bioinformatics/btu031.
Alexa A, Rahnenfuhrer J. topGO: Enrichment analysis for gene ontology; 2016.
Alterovitz G, Xiang M, Mohan M, Ramoni MF. GO PaD: The Gene Ontology partition database. Nucleic Acids Res. 2007;35(SUPPL. 1). https://doi.org/10.1093/nar/gkl799.
Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S, et al. AmiGO: online access to ontology and annotation data. Bioinformatics. 2009;25:288–9. https://doi.org/10.1093/bioinformatics/btn615.
Hu Z-L, Bao J, Reecy J. CateGOrizer: A web-based program to batch analyze Gene Ontology classification categories. Online J Bioinforma. 2008;9 January:108–12. http://users.comcen.com.au/~journals/geneontologyabs2008.htm.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer-Verlag; 2009.
This research was supported in part through computational resources provided by Information Technology at Purdue University, West Lafayette, Indiana. We thank Ian Thompson for growing the cultures, Jessica Cavaletto for extracting the RNA and sending materials to the Purdue Genomics Facility for sequencing, and Michael Gribskov for suggesting improved methods of analysis.
This work was funded by USDA-ARS research project 3602–22000-017-00D. The funding body had no role in designing the experiments, collecting the data or writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
png Quality-control checks of the read libraries. A, To the right of the treatment diagonal: Scatterplots comparing two replicates on a log10 scale; To the left of the treatment diagonal: Individual histograms showing the gene read count distribution over each replicate. Treatments of white, blue and red light, and dark are indicated in the diagonal by W, B, R and D, respectively. Replication number is indicated by an integer from 1 to 3, e.g., W2 is the second replication of the white light treatment. B. Histogram of the gene read sums across all replicates. The x axis is the log10 read count and the y axis is the number of genes. C. Volcano plots of the log2 fold changes versus the adjusted p values.
Additional file 2.
All genes that were differentially expressed in the six possible comparisons between the dark and three light treatments. Each tab contains the differentially expressed genes in a single comparison, where the first treatment listed is the one showing up regulation, e.g., White versus dark shows all genes that were significantly upregulated under white light compared to the dark at a significance level of p = 0.05 or less.
Additional file 3.
Heatmaps of modified adjusted p values for selected Zymoseptoria tritici genes possibly involved in pathogenicity and light sensing and response. The values are the inverse of the adjusted p value multiplied by the direction of the fold change. Comparisons are indicated at the top of each column. The classes A-E are from Palma-Guerrero et al. , where: A is small secreted proteins; B are secreted lipases; C shows plant cell wall degrading enzymes; D summarizes protease genes; and E indicates secreted peroxidases.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
McCorison, C.B., Goodwin, S.B. The wheat pathogen Zymoseptoria tritici senses and responds to different wavelengths of light. BMC Genomics 21, 513 (2020). https://doi.org/10.1186/s12864-020-06899-y
- Gene expression
- Mycosphaerella graminicola
- Zymoseptoria tritici