- 1Institut de Recherche en Horticulture et Semences – UMR 1345, INRA, Université d’Angers, Agrocampus-Ouest, SFR 4207 QUASAV, Angers, France
- 2Université de Lille, INRA, ISA, Université d’Artois, Université du Littoral Côte d’Opale, EA 7394 - ICV - Institut Charles Viollette, Lille, France
- 3Departamento de Agronomía, ETSIAM, Universidad de Córdoba, Córdoba, Spain
- 4Department of Plant Sciences, University of Agronomic Sciences and Veterinary Medicine of Bucharest, Bucharest, Romania
- 5Equipe “Sucres & Echanges Végétaux-Environnement,” UMR CNRS 7267 EBI Ecologie et Biologie des Interactions, Université de Poitiers, Poitiers, France
Alternaria brassicicola is a necrotrophic fungus causing black spot disease and is an economically important seed-borne pathogen of cultivated brassicas. Seed transmission is a crucial component of its parasitic cycle as it promotes long-term survival and dispersal. Recent studies, conducted with the Arabidopsis thaliana/A. brassicicola pathosystem, showed that the level of susceptibility of the fungus to water stress strongly influenced its seed transmission ability. In this study, we gained further insights into the mechanisms involved in the seed infection process by analyzing the transcriptomic and metabolomic responses of germinated spores of A. brassicicola exposed to water stress. Then, the repertoire of putative hydrophilins, a group of proteins that are assumed to be involved in cellular dehydration tolerance, was established in A. brassicicola based on the expression data and additional structural and biochemical criteria. Phenotyping of single deletion mutants deficient for fungal hydrophilin-like proteins showed that they were affected in their transmission to A. thaliana seeds, although their aggressiveness on host vegetative tissues remained intact.
Introduction
The fungus Alternaria brassicicola causes black spot disease and is an economically important seed-borne pathogen of Brassicaceae species. This necrotrophic fungus strongly depends on seed transmission process for its long-term survival and dispersal (van den Bosch et al., 2010). However, fungal and plant factors that impact seed transmission are still poorly described. Such knowledge is crucial to propose strategies for improving the seed health, which remains a major issue for seed companies. Recent studies conducted with the Arabidopsis thaliana/A. brassicicola pathosystem showed that the level of susceptibility of the fungus to water stress strongly influenced its seed transmission ability. For instance, two osmosensitive fungal mutants, defective for the class III Histidine kinase (HK) AbNik1 (Pochon et al., 2013) and the MAP kinase AbHog1 (unpublished result), respectively, were highly jeopardized in their ability to colonize seeds. Consistently, Iacomi-Vasilescu et al. (2008) had previously reported that A. brassicicola spontaneous phenylpyrrole resistant mutants with non-functional class III HK were found to be strongly impaired in their ability to infect radish seeds in field conditions, indicating that a functional high osmolarity pathway is required for efficient infection of seeds. Pochon et al. (2013) also showed that two dehydrin-like proteins were required for effective seed colonization by A. brassicicola. During colonization of maturing seeds, the fungi are exposed to severe water and osmotic constraints, and their ability to cope with such stresses is probably a key factor that determines their transmission to seeds.
The strategies to withstand low water availability have been mainly investigated in some organisms called anhydrobiotes belonging to kingdoms of bacteria, fungi and plants, which have the ability to resist desiccation, i.e., to survive to nearly total dehydration (Calahan et al., 2011; Dupont et al., 2014). Various mechanical, structural and oxidative constraints are induced during dehydration, which may lead to damage to membranes, proteins and DNA. Studies in the yeast Saccharomyces cerevisiae have shown that the protective system is based on a set of constitutive and inducible mechanisms. For instance, the synthesis of polyols occurs to maintain the osmotic balance. In yeast, the major polyol concerned with osmotic adjustment is glycerol and its synthesis, and hence its intracellular accumulation, is under control of the high osmolarity glycerol (HOG) pathway mediated by phosphorylation of the HOG1 MAP kinase. Some other molecules, such as heat shock proteins, hydrophilins or the non-reducing disaccharides trehalose and sucrose have been identified to be synthesized during desiccation and may be involved in preserving native protein structures or stabilizing membranes (Gadd et al., 1987; Garay-Arroyo et al., 2000; Sales et al., 2000; Tapia and Koshland, 2014).
Hydrophilins are a group of proteins that are present in organisms such as plants, fungi and bacteria, and that are biochemically defined as glycine-rich and highly hydrophilic disordered proteins. They are assumed to be involved in cellular dehydration tolerance, even if their specific role is still unclear. In Sa. cerevisiae, 12 genes, whose expression were induced in response to water deficit, encode proteins with specific hydrophilin characteristics (Garay-Arroyo et al., 2000). The deletion (or disruption) and overexpression of these genes showed that only some of them contributed to the desiccation tolerance of yeast (Dang and Hincha, 2011; López-Martínez et al., 2012; Rodríguez-Porrata et al., 2012). A first beneficial effect of hydrophilins on yeast viability during water deficit may be due to their antioxidant capacity which minimizes the accumulation of cellular ROS (López-Martínez et al., 2012). Another role of hydrophilins would be the protection of enzymes and the stabilization of membranes during drying (Sales et al., 2000; Goyal et al., 2005; Reyes et al., 2005). In virtue of their structural plasticity, which is due in part to the presence of disordered regions, hydrophilins may act as molecular shields or directly interact with various molecular partners (Tunnacliffe and Wise, 2007).
Except for the yeast hydrophilins, very little information is available on the presence and role of this class of proteins in filamentous fungi. To date, most studies have focused on a particular class of fungal hydrophilins called dehydrins. Dehydrin-encoding genes were first identified in Tuber borchii during searches for genes controlling fruiting body maturation or conidial dormancy (Abba’ et al., 2006). Then, Wong Sak Hoi et al. (2011, 2012) identified three dehydrin-encoding genes in Aspergillus fumigatus. Two of them were found to be involved in the protection against oxidative, osmotic and pH stress and the third dehydrin regulated freezing tolerance. A farnesol-induced dehydrin-like protein, that contributed to the high tolerance of resting conidia against oxidative and heat stress, was also identified in Aspergillus nidulans (Wartenberg et al., 2012). In A. brassicicola, three proteins harboring the typical asparagine-proline-arginine (DPR) signature pattern and sharing the characteristic features of fungal dehydrin-like proteins have been identified (Pochon et al., 2013). The expression of these genes was found to be regulated by the AbHog1 mitogen-activated protein kinase (MAPK) pathway. Phenotyping of single deletion mutants showed that dehydrin-like proteins have an impact mainly on oxidative stress tolerance and on conidial survival upon exposure to high and freezing temperatures. Moreover, a double deletion mutant was strongly affected with respect to conidiation and showed a highly compromised pathogenicity with a lower aggressiveness on Brassica oleracea leaves and a reduced capacity to be transmitted to Arabidopsis seeds via siliques. These results support the importance of dehydrin proteins and, more globally, of hydrophilins, with respect to the ability of A. brassicicola to efficiently accomplish key steps of its pathogen life cycle.
To identify new fungal molecular determinants that may be involved in the silique and seed infection process, the present study aimed at analyzing the transcriptomic and metabolomic responses of A. brassicicola exposed to water stress. Then, the repertoire of putative hydrophilins in this fungus was established based on expression data and structural and biochemical criteria. Finally, functional approaches were initiated to investigate the role of some of these putative hydrophilins.
Results
Application of Water Stresses
Several preliminary tests were conducted to evaluate the level of stress for each treatment applied to fungal cells. Considering that the reference strain is naturally tolerant to water stress (Iacomi-Vasilescu et al., 2008), we additionally analyzed the effect of treatments on three mutant strains (Δabnik1, Δabhog1 and Δabsch9) derived from the wild-type Abra43 strain, for which we hypothesized a higher sensitivity to water stress. As expected, the Abra43 reference strain exhibited high sorbitol tolerance with only 5% inhibition at a concentration of 1.2 M. At this concentration, a higher susceptibility was observed for the mutant strains (Table 1), confirming that the application of such a concentration indeed constitutes a stress for the fungal cells. In the same way, when a drastic decrease in relative humidity (RH) was applied to the germlings, no significant difference in mortality between the treated and the silica gel-free conditions was observed for the reference strain. On the contrary, the results showed a very high susceptibility of Δabnik1 and Δabhog1 to the desiccation stress, while the AbSch9 mutant exhibited an intermediate susceptibility (Table 1).
Table 1. Susceptibility of Alternaria brassicicola wild-type, ΔAbnik1,ΔAbhog1 and ΔAsch9 strains to sorbitol and desiccation stresses.
Based on these results, we chose to analyze the responses of the fungus when exposed to 1.2 M sorbitol for 0.5 and 2 h and when exposed to desiccation stress for 1 and 4 h. Regarding the latter stress, the chosen times implied that germlings exposed to silica gel beads were subjected to a RH value of 1% while control germlings were harvested at RH values of 70% and 90%, respectively.
In addition, we analyzed other phenotypic characteristics of these mutants: the susceptibility to oxidative stress (exposure to 10 mM H2O2 and 20 mM menadione), the ability to transmit to the A. thaliana seeds and the aggressiveness on cabbage leaves. Δabnik1 and Δabhog1 mutants were more susceptible to oxidative stress, and, consistently, had reduced ability to colonize A. thaliana seeds compared to the wild-type strain (Table 2). The Δabsch9 mutant was not susceptible to oxidative stress and its seed transmission capacity was not altered. Compared to seed transmission capacity, the ability of leaf infection was not or less dependent on the sensitivity of fungal genotypes to oxidative stress, since Δabnik1 and Δabsch9 strains did not have reduced aggressiveness on cabbage leaves.
Modulation of Fungal Transcriptomes in Response to Water Stress
In order to identify potential effectors of the response of A. brassicicola to dehydration stress, we first focused on the analysis of the transcriptional response in germinating conidia exposed to previously defined water stress (i.e., exposure to sorbitol 1.2 M for 0.5 and 2 h and exposure to silica gel beads for 1 h and 4 h). Each treated sample was compared to control sample and, thus, we produced four transcriptome data sets (each set was obtained from three biological replicates). We used A. brassicicola Nimblegen microarrays bearing one probe for each of the 10 633 ORFs predicted in the A. brassicicola automatically annotated genome database (JGI Genome Portal1). Only probes with a P-values ≤0.05 and a log2 ratio ≥0.7 or ≤−0.7 were considered as differentially expressed. The numbers of genes that were induced or repressed under each condition, together with a summary of those regulated in more than one condition, are shown in Figure 1. To validate our microarray results, the regulation of several genes by sorbitol or silica gel exposure was confirmed by quantitative PCR (Supplementary Figure S1).
Figure 1. Global modulation of gene expression in Alternaria brassicicola in response to applied treatments. Sorbitol and desiccation treatments were compared to non-treated cultures at two time points 0.5 h, 2 h and 1 h, 4 h, respectively. Genes with a P-values ≤0.05 and a log2 ratio ≥0.7 or ≤–0.7 were considered as differentially expressed. Intersection with less than 10 genes were omitted. Colors displayed regulation compared to the control; red = up-regulated in all conditions, green = down-regulated in all conditions, orange = opposite regulation according to the considered condition.
Our data showed that, for all the treatments, with the exception of the 4 h-silica gel (4 h-desiccation) treatment, the number of up-regulated and down-regulated genes was less than 500 and 300, respectively. Most of these genes were found to be regulated in more than one condition, 37 were found to be induced in all 4 stress conditions while none were found to be repressed. The 4h-desiccation stress led to more profound transcriptomic modifications since 578 and 917 were induced or repressed compared to control condition, respectively. Moreover, this transcriptomic response seemed rather specific since more than 800 genes (268 up-regulated and 601 down-regulated genes) were specifically regulated after the fungus was exposed to this constraint. It should be noted that 157 genes, that were induced by the sorbitol exposure, were repressed following the silica gel treatment.
The GO enrichment analyses showed that the enriched GO terms were not generally shared within the gene sets representative of each type of treatment (sorbitol or desiccation) (Figure 2). One exception was the GO categories corresponding to the response to oxidative and osmotic stresses, which were enriched in all the lists of induced genes. Another configuration was that enriched categories emerging from the 4h-desiccation up-regulated genes were also enriched categories emerging from the down-regulated gene lists after sorbitol exposure. These results seem to point out that the two types of treatment did not cause the same effects at the cellular level and induced quite different transcriptomic responses.
Figure 2. Relevant Gene Ontology Biological Process (GO-BP) enrichments obtained from modulated gene sets after exposure to water stresses.
Enriched categories emerging from genes regulated following silica gel exposure were mainly related to the translation process, mitochondrial functions, and carboxylic acid and amino acids metabolic processes. Within the up-regulated genes, the main enriched GO terms identified in the sorbitol data sets were predominantly associated with a response to oxidative stress and to sugar metabolic process. Within the sorbitol repressed gene list, enriched GO terms referred to a wider range of cellular processes such as gene expression (Figure 2).
Metabolomic Responses to Water Stress
As several enriched GO categories emerging from transcriptomic data were linked to amino acid metabolism and catabolism, we investigated whether the response to water stress in A. brassicicola resulted in an alteration of the cellular content in amino acids. Regardless of the hydric stress applied, the total amino acid concentrations were significantly lower for the treatments conditions compared to the control conditions (Welch test P-value = 5.955E–08, Supplementary Figure S2). Indeed, considering the 24 amino acids and derivatives that were quantified in this study, the mean ratio of amino acid concentration in the treated conditions compared to the respective controls decreased significantly for 17 and 12 amino acids (9 in common), regardless of time exposure for the sorbitol and the desiccation treatment, respectively (Figure 3A). Amino acid profiles that were differentially less accumulated compared to controls were quite similar after 0.5 and 2 h of sorbitol exposure with 9 shared amino acids (over 13). Glutamine, glutamate, alpha alanine, lysine and arginine concentration explained 70% [187.8/270.0 μmol.g–1 DW (dry weight)] and 78% (156.3/199.8 μmol.g–1 DW) of the amino acid variations for the sorbitol treatment at 0.5 and 2 h, respectively. For the 4 h desiccation condition, glutamine, GABA, alpha alanine and lysine explained 70% (134.7/194.6 μmol.g–1 DW) of the amino acid variations. No significant concentration differences were observed for the 1 h desiccation treatment.
Figure 3. Impact of sorbitol and desiccation treatments on amino acid (A), saccharides and sugar alcohols (B) contents of fungal cultures. The heat map background displays the mean control contents. Internal colored boxes display the mean difference between treatments and controls. Quantifications were performed from three separate samples. Asterisks indicate a significant difference between the mutant and the parental isolate (Welch test, ∗p ≤ 0.05).
Then, we focused on sugar content, since their intracellular accumulation is a well-known protective mechanism to cope with the water constraints (Dupont et al., 2014). Eight saccharides and sugar alcohols were quantified in treated and control conditions (Figure 3B). The total concentration of those compounds was found to be significantly different only in the case of comparison between sorbitol-treated and untreated conditions (Welch test P-value = 2.05E–5, Supplementary Figure S2) regardless of the time exposure. Unexpectedly, we observed a strong accumulation of some sugars (arabitol, trehalose, mannitol) in the control of the desiccation condition. This result could explain that no significant difference on the sugar cellular contents was observed between treated and untreated samples, except for an arabitol accumulation after a 4 h-treatment.
As expected, glycerol was accumulated in fungal cells exposed to sorbitol. Two other polyols, arabitol and dulcitol, were also strongly accumulated under these conditions. Conversely, the amounts of other polyols and the disaccharide trehalose were lower in cells treated with sorbitol than in untreated conditions.
Hydrophilins Identification Through in silico Analysis
Hydrophilins could be defined as glycine-rich, highly hydrophilic disordered proteins. This criterion was used to search in A. brassicicola for putative hydrophilins among the whole protein set deduced from the recently published genome of A. brassicicola strain Abra43 (12,456 entries, Belmas et al., 2018). A total of 732 predicted proteins with an average hydrophilicity above one (GRAVY index lower than −1) were first selected and analyzed for their glycine content. The 210 predicted proteins with at least 8% glycine were retained and subjected to secondary structure prediction analysis. Almost three quarters of them (165) were predicted to adopt a coil configuration over more than 50% of their length. After removing mis-annotated and incomplete proteins and those that had an obvious match upon Blast analysis with proteins of known functions (except “stress response”), a final set of 107 candidate hydrophilins (c.a. 0.9% of the whole predicted proteome) was selected (Table 3). Interestingly, among proteins that had physicochemical features of hydrophilins but were discarded after Blast analysis, several displayed similarities with proteins involved in nucleic acids metabolic processes (Supplementary Table S1). Each of the 107 candidate hydrophilins was characterized by an index of significance, calculated based on the following formula: [- (gravy index −1) X (% G – 8) X (% rc – 50)] X 100, whose value varied from 3 to 36.512. Six putative homologs of fungal hydrophilins were found within the whole set of A. brassicicola hydrophilin-like proteins: HSP12 (AB01868/Abra01064 and AB10131/Abra12207), STF2 (AB03803/Abra04061) from Sa. cerevisiae, CON10 (AB02056/Abra00844 and AB01782/Abra01163) and CON6 (Abra02927) from Neurospora crassa, and the dehydrin-like proteins from A. brassicicola AbDhn1 (AB02513/Abra00274) and AbDhn3 (AB05365, AB05364/Abra08036). None of the other yeast hydrophilins had putative homolog in the A. brassicicola genome except NOP6 (AB08803/Abra04570) and WWM1 (AB05865/Abra11690) that were not selected due to a glycine content below 8% or a gravy index above −1, respectively. Similarly, due to an average hydrophilicity below the threshold value, the other A. brassicicola dehydrin-like protein AbDhn2 (AB08993/Abra04787) was not selected within the candidate hydrophilins set. This protein set also contains the putative homolog DNA-damage responsive DDR48 (AB03330/Abra05531) protein from yeast.
Table 3. Set of candidate hydrophilins identified from the Alternaria brassicicola genome sequence through in silico analysis.
Expression of Hydrophilins During Osmotic and Hydric Stress
Then, we explored the expression pattern of genes encoding putative hydrophilins from transcriptome data sets generated after exposing germinating conidia of A. brassicicola to 1.2 M sorbitol and to low relative humidity (silica gel beads). These data sets were obtained using microarrays designed based on the genome of the A. brassicicola ATCC 96836 strain available on the JGI Genome Portal (see footnote 1). This genomic sequence contained 10,514 predicted genes and the proportion of gaps between scaffolding boards in this sequence was relatively high (estimated to 5.2% by Dang et al., 2015). As described above, the in silico analysis was performed from the recently published genome of A. brassicicola strain Abra43 (Belmas et al., 20182 ,3), which is a higher quality sequence containing 12,456 predicted protein-coding genes. Considering these proteome differences, expression patterns were available only for 78 putative hydrophilins that matched with a predicted ORF from the ATCC 96836 strain proteome. Under these conditions, 17 of the 78 (21.8%) candidate hydrophilin genes were up-regulated by the sorbitol treatment (Figure 4A). For all of them, except AB03330 and AB00441, an induced expression was observed at 0.5 h (Figure 4B). Moreover, 13 candidate hydrophilin genes were up-regulated by the desiccation treatment, 10 of them being also up-regulated upon the sorbitol treatment. All of these proteins were named with the abbreviation “SIH” for “stress-induced hydrophilin-like” (Table 3).
Figure 4. Gene expression of Alternaria brassicicola hydrophilins in response to sorbitol and desiccation treatments. Sorbitol and desiccation treatments were compared to non-treated cultures at two time points 0.5 h, 2 h and 1 h, 4 h, respectively. (A) For each gene, numerical data were transformed into color-grid representations in which the fold gene expression induction (log2 values) is represented by a colored scale. In addition, hydrophilin gene expression were evaluated in mutant strains Δabnik1, Δabhog1 and Δabsch9 exposed to 1.2 M sorbitol for 0.5 h. The left part of the panel (A) highlights the possible regulation of hydrophilin gene expression by one or several of these protein kinases. (B) Volcano plot showing fold inductions of genes in the wild-type strain exposed to 1.2 M sorbitol for 0.5 h. The hydrophilin encoding genes are labeled and identified with red dots. The threshold values (P-values ≤0.05 and a log2 ratio ≥0.7 or ≤–0.7) are indicated by dash lines.
To gain insight on the possible regulation of hydrophilin genes expression by different kinases known to regulate the osmostress response, the expression profiles of hydrophilin genes in the Δabhog1, Δabnik1 and Δabsch9 mutants were compared with those in the wild-type parental strain under both control condition and after exposure to sorbitol for 0.5 h (Figure 4). This analysis revealed that, among the sorbitol-inducible hydrophilin genes, 6 were dependent of at least two of the three kinases for their expression. Two were regulated by AbHog1 and AbNik1, two others by AbHog1- and AbSch9 and three others by AbHog1, AbNik1 and AbSch9. The gene encoding AbDhn2 was also induced by the treatment and was found Hog1-Nilk1-Sch9-dependent.
Role of Selected Hydrophilins in A. brassicicola
Two putative hydrophilins (AbSih3/AB01868 and Absih15/AB05148) were selected for further studies aiming to investigate their function in A. brassicicola. Knockout mutants deficient for AbSih3 and Absih15 (called Δabsih3 and Δabsih15), respectively, were constructed by replacing the respective ORFs with a hygromycin B resistance cassette (Supplementary Figure S3). For each targeted gene, two replacement mutants were selected using a PCR screen and further purified by two rounds of single-spore isolation. In all further experiments, the phenotypic characters for transformants of the same genotype were not found to be significantly different and the phenotyping values that we displayed were calculated by taking into account both transformants.
Although the expression of these genes was induced in germinating conidia exposed to previously defined water stress (i.e., exposure to 1.2 M sorbitol and to silica gel beads), the knock-out approach showed that only the AbSih3 deletion had an impact on phenotype, mainly on oxidative stress tolerance. The corresponding deficient single mutants were indeed characterized by a slightly increased susceptibility toward oxidative stress generated by exposure to menadione and H2O2 compared to the wild-type (Figure 5). No significant difference in susceptibility to other tested stresses (exposure to PEG −0.7 MPa, 350 mM NaCl, 1.2 M sorbitol) was found compared to the wild-type (data not shown).
Figure 5. Growth inhibition rates of the wild-type strain (WT) and hydrophilin deficient mutants in the presence of H2O2 (10 mM) and menadione (20 mM). The results are expressed as the percentage of inhibition in treated samples compared to the control without additive. Conidia were used to inoculate microplate wells containing standard PDB medium that was supplemented with the appropriate test substance. Growth was automatically recorded for 30 h at 25°C using a nephelometric reader (see section “Materials and Methods”). Each genotype was analyzed in triplicate and the experiments were repeated at least three times per growth condition. Error bars indicate standard deviations. Asterisks indicate a significant difference between the mutant and the parental isolate (Dunn control test, ∗p ≤ 0.01).
Consistently, a similar higher susceptibility to H2O2 of the Sa. cerevisiae hsp12 null mutant was observed (Figure 6). To confirm that AbSih3 correspond to true functional homolog of the yeast hsp12 gene, complementation of hsp12 mutant null mutant was carried out with relevant recombinant pYES plasmids. As shown in Figure 6, AbSih3 successfully complemented the Sa. cerevisiae hsp12 phenotypes, even if the growth of Δhsp12 cells expressing AbSih3 was not restored to a level comparable to that of the WT strain containing the empty vector pYES2 on media supplemented with 1 mM H2O2.
Figure 6. Complementation of yeast Δhsp12 mutants with AbSih3 coding sequences. (A) Growth curves of Saccharomyces cerevisiae (parental WT strain and Δhsp12 mutants) cells transformed with empty vector (pYES2CT) or recombinant vectors (pYAbsih3) exposed to 1 mM H2O2 or to water. (B) Comparisons of growth inhibition rates obtained for each yeast strain (Wilcoxon test, ∗p < 0.001).
Then, the virulence of the A. brassicicola wild-type, Δabsih3 and Δabsih15 strains were compared on B. oleracea and A. thaliana host plants during host vegetative tissue infection or during the seed transmission process. Brassica oleracea leaves were inoculated with drops of conidia suspension (105, 104, or 103 conidia/mL). As shown in Figure 7, the wild-type and KO strains were all able to produce typical symptoms and their aggressiveness was not impaired. Regardless of the inoculated strain, first necrotic symptoms appeared on leaves at 3 days post-inoculation (dpi). Necrotic areas enlarged and surrounded themselves with chlorotic halos at 6 dpi. During late stages of infection, the fungus produced conidia on the surface of necrosis.
Figure 7. Representative symptoms at 5 dpi obtained by inoculation of wildtype, Δabsch9, Δabsih3, and Δabsih15 on cabbage leaves. Leaves were inoculated with 5 μL drops of conidia suspension (105, 104 or 103 conidia/mL in water). As shown on the right side, mutants were inoculated on the right part of the central vein and compared on the same leaf with the parental strain (inoculated on the left part of the central vein). Bars = 2 cm.
Using the model pathosystem previously described by Pochon et al. (2012) for investigating seed transmission in Arabidopsis plants, the ability of the hydrophilin mutants to transmit to seeds was compared with that of the wild-type. As a positive control, we also analyzed the behavior of the osmosensitive fungal mutants Δabnik1 and Δabhog1, which were previously described as being strongly altered in their seed transmission capacity (Pochon et al., 2012). Silique inoculation with the wild-type and hydrophilin mutant strains resulted in the development of mycelium and typical lesions covered with conidia on siliques within a few days after inoculation. Seeds were then individually harvested and plated on PDA medium for a microbiological examination. As expected, a strong decrease in seed transmission rate was observed for Δabnik1 and Δabhog1 compared to the wild-type parental strain (Figure 8). Interestingly, the seed transmission capacity of the Δabsch9 mutant was not altered although this strain was found to be sensitive to sorbitol and desiccation. While the deletion of AbSih3 had a slight significant impact on seed transmission, a stronger decrease was observed for the Δabsih15 mutant, for which the seed transmission rate was reduced by a factor two compared to the wild-type strain.
Figure 8. Transmission capacity of Alternaria brassicicola wild-type (WT) and null strains to Arabidopsis thaliana seeds (Ler ecotype). The seed transmission capacity were measured as described by Pochon et al. (2012). The five youngest siliques of at least five plants were inoculated with each fungal genotype and the experiment was repeated twice. Contaminated siliques were harvested 10 dpi. After dissection, seeds were incubated separately on PDA medium for 2 days. A seed was considered contaminated when incubation resulted in typical Alternaria brassicicola colony development. For each inoculated fungal genotype, the seed infection probability was evaluated from at least 1000 seeds. Values represent infection probabilities with 95% confidence interval.
Discussion
The ability of A. brassicicola to cope with water-stress constitutes a key factor for successful seed colonization and thus accomplishment of its entire infectious cycle (Iacomi-Vasilescu et al., 2008; Pochon et al., 2012, 2013). Highly hydrophilic and disordered glycin-rich proteins containing repeated DPR-motifs called dehydrins were identified in A. brassicicola and might correspond to effectors of the response to water stress (Pochon et al., 2013). Based on their physicochemical features, dehydrins are considered as members of a larger protein family called hydrophilins initially described in Sa. cerevisiae (Garay-Arroyo et al., 2000). In an attempt to identify additional effectors of the response of A. brassicicola to water stress, the proteome deduced from the recently published genome of A. brassicicola strain Abra43 (Belmas et al., 2018) was first screened for proteins with characteristic features of hydrophilins. Candidate hydrophilins accounted for 0.9% of the whole predicted proteome. In comparison, hydrophilins represented only 0.2% of the yeast proteome (Garay-Arroyo et al., 2000) although the selection criteria were less stringent, i.e., proteins with lower glycine content were included. This difference might reflect a family expansion associated with environmental adaptation and/or fungal pathogenesis as it was proposed for other sets of proteins, e.g., secreted proteins, enzymes for secondary metabolites production (Soanes et al., 2008; Ohm et al., 2012). However, this assumption remains to be demonstrated and would require a more extensive search for hydrophilin repertoires in fungi with various lifestyles.
The fact that two of the three previously described dehydrins and additional homologs of fungal hydrophilins were included in the repertoire of candidate hydrophilins validated our screening process. Garay-Arroyo et al. (2000) proposed that hydrophilins that were not induced by water-stress conditions should not be considered as true hydrophilins. Thus, we analyzed changes in transcript levels in A. brassicicola germlings exposed to two kinds of water-stress: hyperosmosis (exposure to 1.2 M sorbitol) and desiccation (exposure to silica gel beads). Using gene expression array, which was designed from a previous genome sequence with a lower quality, we defined the expression patterns of 78 genes (over 107) encoding these putative hydrophilins. Twenty candidate hydrophilin genes were up-regulated either by the sorbitol treatment or by the desiccation stress, and, among them, ten were induced by both treatments. So, according to the hydrophilin definition given by Garay-Arroyo et al. (2000), we identified at least 20 hydrophilins in A. brassicicola. However, we can reasonably assume that their number is higher since we did not have access to the gene expression of 29 other genes encoding putative hydrophilins.
In this study, strains deficient for AbSih3 exhibited a slightly lower susceptibility to oxidative stress compared to the wild-type strain. No significant difference in susceptibility to other tested stresses (exposure to PEG, NaCl and sorbitol) was found for Δabsih3 and Δabsih15. Given the large number of hydrophilins that were co-expressed in this organism, functional redundancy may occur between some of these proteins. This process would explain the scarcity of significant phenotypes that we obtained from single deficient mutants since the presence of other hydrophilin genes would compensates for the loss. For instance, AbSih3 displays high sequence similarity to AbSih18 (AB10131), and both may be considered as homologs of the Sa. cerevisiae HSP12 hydrophilin, even if AbSih3 only partially compensated the HSP12 loss in yeast. The existence of such a functional redundancy was also previously observed between two dehydrins AbDhn1 and AbDhn2 in A. brassicicola. The double mutant ΔΔabdhn1-abdhn2 was severely compromised in its pathogenicity while the respective single mutants were not affected (Pochon et al., 2013).
Despite this probable functional compensation, we showed here that the hydrophilins Absih3 and Absih15 could be considered as pathogenic factors since the corresponding null strains were affected in its seed transmission ability (more strongly for the Δabsih15 strain). This result confirmed our initial assumption that an efficient fungal adaptive response to severe hydric stress conditions or to oxidative stress strongly influenced the seed transmission ability. This is consistent with the fact that, during seed colonization, fungi are exposed to a gradual decrease in the water potential occurring in maturing reproductive organs and that they have to cope with this particular stress to be efficiently transmitted to seeds and complete their infection cycle. In the same vein, we showed in this study or in previous works that the class III HK AbNik1 and the MAP kinase AbHog1, two components of the high osmolarity pathway, were required for an effective seed colonization by A. brassicicola (Pochon et al., 2012). The Δabsch9 mutant was also found to be sensitive to sorbitol and desiccation but, unexpectedly, its seed transmission capacity was not altered. The fact that this mutant strain did not exhibit a high susceptibility to oxidants suggested that the fungal ability to overcome oxidative stress would be an additional or the main key mechanism for an efficient seed infection process. This hypothesis is consistent with the fact that oxidative damages produced by free radical species are considered to be one of the major causes of desiccation injuries (Dupont et al., 2014).
A final point raised from these results is that several A. brassicicola mutants (Δabnik1, Δabsih3 and Δabsih15) were affected in their transmission to A. thaliana seeds, although their aggressiveness on host vegetative tissues remained intact. This fact confirms the specificity of some molecular mechanisms controlling transgenerational transmission of fungal pathogens (Pochon et al., 2012) and validates the strategy to gain further insight into these mechanisms for the development of specific seed disease control strategies.
In addition to the induction of genes encoding hydrophilins, we identified other metabolic clues related to the activation of an adaptive response to water stress in A. brassicicola. First, the cellular content in amino acids was strongly disturbed following the application of both stresses. The amino acid levels were indeed significantly reduced in treated cells compared to controls. In line with this result, several enriched GO categories emerging from transcriptomic data were linked to amino acid metabolism and catabolism. Moreover, amino acid synthesis was also found to be strongly down-regulated in response to the osmotic stress treatment in another Dothideomycete fungus, the wheat pathogen Stagonospora nodorum (Lowe et al., 2008). As reported in plants (Hildebrandt et al., 2015), the degradation of amino acids may lead to the production of an additional energetic resource, that could be useful to cope with adverse conditions, by providing intermediates of the tricarboxylic acid cycle and finally contributing to ATP production. Moreover, the oxidation of some amino acids feeds electrons into the mitochondrial electron transport chain. Most reaction steps of the catabolic pathways occur either in the mitochondria or cytosol. Consistent with this hypothesis, several enriched GO terms referred to mitochondria organization, carboxylic acid metabolism and to oxidation-reduction process. Regarding the response to desiccation, as a GO enrichment of functions related to the genesis and transport of ribosomes and, more generally, to translation was observed, we could assume that a strong induction of the translation process would trigger a depletion of the cellular amino acid pool. So, the mechanisms leading to the decrease of the amino acid could depend on the type of applied stress.
An intriguing result from our data was that no significant difference in the amino acid concentrations was observed for the 1 h-desiccation treatment compared to the respective control. A more precise inspection showed that, for this condition, both control and treated samples showed a similar low amino acid content (Supplementary Figure S2). In this assay, germlings grown on cellophane were taken off from the PDA plate and directly transferred to sealed boxes (without additional agar medium) with or without silica gel beads. In the control box (without silica gel), we suspect that germlings were also exposed to a strong constraint in the early time of exposure, as a period of 4 h was required to gradually reach a RH value of about 90%. Therefore, although this experimental device triggered a specific regulation at the transcriptome level (Figures 1, 2), this does not appear to be sufficient to reveal significant differences in amino acid accumulation between control and treated samples at this analysis time.
In the sorbitol-exposed samples, a strong increase in global sugar contents was observed compared to sorbitol control. Interestingly, the accumulation of non-reducing disaccharides, such as trehalose, did not seem to be favored in A. brassicicola germlings under the sorbitol stress conditions. On the contrary, the contents of three polyols, glycerol, arabitol and dulcitol, were mainly impacted by this stress, suggesting that they function as intracellular osmolytes to counteract the effects of decreases in the cell volume and the loss of turgor pressure (Sánchez-Fresneda et al., 2013; Dupont et al., 2014). This result was consistent with the increased sensitivity of the Δabhog1 mutant to osmotic stress since the MAP kinase Hog1 is known to be required for endogenous glycerol production in several fungal species (Brewster and Gustin, 2014; Hohmann, 2015). However, its involvement in the synthesis of other protective compounds, such as arabitol or dulcitol, is still unclear (Kayingo and Wong, 2005; Sánchez-Fresneda et al., 2013). Arabitol has been described as part of a set of polyols which is developed to counteract environmental challenges in several yeasts and filamentous fungi (Pascual et al., 2003; Sánchez-Fresneda et al., 2013). For instance, arabitol, and to a lesser extent glycerol, were found to accumulate in response to the osmotic stress treatment in the wheat pathogen St. nodorum (Lowe et al., 2008) or in the rice pathogen Magnaporthe grisea (Dixon et al., 1999). Conversely, the role of dulcitol in stress responses remains highly confidential and has been reported on rare occasions, for instance in the marine fungus Cirrenalia pygmea (Ravishankar et al., 2006).
In addition to its involvement in the synthesis of compatible solutes, we showed here that the HOG pathway, mainly regulated by the MAP kinase AbHog1 and the upstream HK AbNik1, also acted as an expression regulator of some hydrophilin encoding genes. The expression profiles of hydrophilin genes in the Δabhog1, Δabnik1, and Δabsch9 mutants were indeed compared with those in the wild-type parental strain after exposure to sorbitol for 0.5 h. Of the 17 proteins that were induced under these conditions, 11 and 5 were found to be under the control of AbHog1 and AbNik1, respectively. This result, associated with the weakness of Δabnik1 and Δabhog1 to survive on reproductive organs, highlighted that the HOG pathway is a major target for controlling both the synthesis of osmolytes and hydrophilins, and thus limiting the transmission of fungal pathogens to seeds. Unexpectedly, AbDhn1 and AbDhn3 were not part of this set of HOG-regulated proteins although Pochon et al. (2013) previously reported that these two proteins (in addition to AbDhn2) were targets for the AbHog1 signaling cascade when the fungus was exposed to a saline (NaCl) stress. This suggests a possible differential regulation of hydrophilins depending on the applied constraint.
Following a 0.5 h-sorbitol exposure, six other hydrophilins were also found to be AbSch9-dependent. As a major downstream effector of the target of rapamycin (TOR) complex 1, the Sch9 protein kinase was proposed to play multiple roles in stress resistance and nutrient sensing in Sa. cerevisiae (Pascual-Ahuir and Proft, 2007; Huber et al., 2009). However, the functions of Sch9 orthologs in filamentous fungi have been poorly characterized to date. In Fusarium graminearum, Gu et al. (2015) showed that FgSch9 served as a mediator of the TOR and HOG pathways and played important roles in regulating vegetative differentiation, secondary metabolism and stress responses. We showed here that both hydrophilins AbSih3 and AbSih15 were necessary factors for an optimal A. brassicicola transmission to seeds but their expression was not found to be regulated by AbSch9 (regulation by AbHog1 and AbHog1-AbNik1, respectively). This result, associated with the fact that the seed transmission capacity of the Δabsch9 mutant was not affected, suggested that this kinase (unlike AbHog1 and AbNik1), although involved in resistance to sorbitol and desiccation, is not a major player in transgenerational transmission of fungal pathogens.
Materials and Methods
Strains and Culture Conditions
The A. brassicicola wild-type strain Abra43 used in this study was previously described (Sellam et al., 2007). For routine culture, A. brassicicola was grown and maintained on potato dextrose agar (PDA). The genome of Abra43 strain was recently published and is available following this link https://bbric-pipelines.toulouse.inra.fr/myGenomeBrowser? (Carrere and Gouzy, 2017).
For transcriptomic analyses, exposure to sorbitol (osmotic stress) or to low relative humidity (desiccation stress) were achieved as follows. Conidia were first incubated for 24 h on cellophane disks overlaid on solid PDA medium. Germlings and the underlying cellophanes were then taken off from the agar plate and transferred either to a second PDA medium supplemented or not with 1.2 M sorbitol for 0.5 or 2 h.
The desiccation stress was achieved by applying a drastic decrease in relative humidity (RH) to the germlings (obtained as described below) by exposure to silica gel beads in sealed boxes for 0.5 or 4 h. Immediately after deposition of germlings (and the underlying cellophanes) in boxes, a strong RH increase was observed. In the control box (without silica gel), a gradual increase up to a RH value of about 90% in 4 h was observed. With silica gel beads, the opening of the boxes resulted in a rapid RH increase up to 20% followed by a decrease to a RH value of 1% in less than 0.5 h.
To evaluate the level of stress for each treatment applied to fungal cells and considering that the reference strain is naturally tolerant to water stress (Iacomi-Vasilescu et al., 2008), we additionally analyzed the effect of treatments on three mutant strains (derived from the wild-type Abra43 strain) for which we hypothesized a higher sensitivity to water stress. Two mutant strains (Δabnik1 and Δabhog1) were described in previous studies (Dongo et al., 2009; Joubert et al., 2011) and were deficient for genes encoding key regulators of the HOG pathway, the HK AbNik1 and the MAP kinase AbHog1, respectively. The third mutant was obtained in this study and was deficient for the kinase AbSch9, which is also known to be involved in the fungal response to osmotic constraints (Pascual-Ahuir and Proft, 2007). Under non-stress conditions, no significant differences in mycelium growth were detected in any of the tested mutants as compared to the wild-type (data not shown). To assess the impact of low RH exposure on A. brassicicola, viability tests were conducted on young hyphae of wild-type and mutant strains. About 200 conidia were first incubated for 24 h on cellophane disks overlaid on solid PDA medium. Then, germinated conidia were transferred for 4 days to sealed boxes containing or not silica gel beads. The cellophane membranes were then removed from the boxes and deposited again on PDA medium. The survival rate was calculated by counting the number of mycelial colonies (derived from treated or untreated germlings) after 3 days of incubation at 24°C. As previouslys reported (de Lima Alves et al., 2015; Wang et al., 2017), exposure to these two treatments induces water stress in fungal cells. Exposure of wild-type and mutant strains to 1.2 M sorbitol were performed following the procedure described below. The radial growth was calculated after 6 days of incubation.
To study hyphal growth in liquid media, conidial suspensions (105 conidia/mL, final concentration) were inoculated onto microplate wells containing the considered substance (10 mM H2O2 or 20 mM menadione) at the desired concentrations in potato dextrose broth (PDB, Difco) in a total volume of 300 μl. Fungal growth was monitored automatically over a 30 h period with a laser-based microplate nephelometer (NEPHELOstar, BMG Labtech) as described by Joubert et al. (2010). For each condition, the area under the growth curve representative of the lag phase was used to calculate the percentages of growth inhibition under stress conditions compared with standard growth conditions (PDB medium without additive), as described by Gaucher et al. (2013). At least three replicates were conducted per treatment.
RNA Extraction and Microarray Analysis
Total RNA was extracted using the nucleospin RNAplant kit (Macherey Nagel, Düren, Germany) and 400 ng were amplified using the Ambion messageAmp II (Ambion, Austin, TX, United States) according to the manufacturer’s instructions. Five micrograms of amplified RNA were retro-transcribed with 400 U of Superscript II reverse-transcriptase (Invitrogen Corp., Carlsbad, CA, United States) and labeled with 1.5 mmol of Cyanine-3 (Cy3) or Cyanine-5 (Cy5) (Interchim, France) and then purified using NucleoSpin Gel and PCR Clean-up column kits (Macherey-Nagel, GmbH and Co., KG, Germany). Purified and labeled cDNA were quantified using a NanoDrop ND-1000. Corresponding Cy3- and Cy5 labeled samples (30 pmol) were combined and cohybridized to the Abra _v1.0 12 × 135K arrays on a NimbleGen Hybridization System 4 (mix mode B) at 42°C overnight. The Abra_v1.0 chip was in situ synthesized by Nimblegen (Madison, WI) and contained 10,633 60-mer oligoprobes that were designed from the A. brassicicola automatically annotated genome database4, containing one probe per sequence. Slides were washed, dried, and scanned at 2 μm resolution and high sensitivity with a Roche-NimbleGen MS200. DEVA Software version 1.2.1 was used to extract pair data files from the scanned images. Three biological replicates were analyzed per comparison using the dye-switch method. Statistical analyses of the gene expression data were performed using the R language, version 2.14, and the Linear Models for Microarray Analysis package (Smyth et al., 2005) from the Bioconductor project. For the preprocessing step, data were normalized by the LOWESS method. Log ratio and log intensity were calculated before differential expression analyses and performed using the lmFit function and the Bayes-moderated Student’s t-test in a linear model for microarray analysis. Probes with P-value <0.05 were considered as differentially expressed. Gene expression datasets were deposited in the Gene Expression Omnibus (GEO)5.
Quantitative PCR
Total RNA was extracted as described above. Amplification experiments were conducted as previously described (Joubert et al., 2011). The relative quantification analysis was performed using the ΔΔCt method (Winer et al., 1999). All the primer sequences used in real-time quantitative PCR are presented in the Supplementary Table S2.
Generation of Targeted Gene Knockout Mutants and Fusion Strains
We followed the procedure described by Pigné et al. (2017). Gene replacement cassettes were generated with the Hph gene cassette from pCB1636 (Sweigard et al., 1995) conferring resistance to hygromycin B. The sets of primers used to amplify the 5′ and 3′ flanking regions of each targeted gene are presented in the Supplementary Table S2. The replacement cassettes were used to transform A. brassicicola protoplasts as described by Cho et al. (2006). Primer combinations used to confirm integration of the replacement cassette at the targeted locus were presented in Supplementary Figure S3.
Yeast Complementation Assays
Complementation assays were performed in Sa. cerevisiae cells (strains: parental BY4743 Acc#Y20000, ΔHSP12 Acc#Y27070, Euroscarf, Germany) as described by Joubert et al. (2011). Sequences encoding AbSih3 were amplified by PCR with primer pairs AbSih3-Bam (5′- GGATCCAAAAATGACCGACGCATTCCGC-3′) – AbSih3-Eco (5′- GAATTCCCATTCCAAGGGCGTTTTTGG-3′). Unique BamH1-EcoRI restriction sites (indicated in boldface) were thus introduced at the upstream start and downstream of the coding sequences, respectively and were used to clone the PCR product into pGEM-T (Promega) and then into the yeast expression vector pYES2CT (Invitrogen) to generate the recombinant plasmid pYAbhsp12. Sensitivity to 1 mM H2O2 was analyzed based on microcultivation in 300 μl liquid medium. Growth in microplates was automatically recorded using a spectrophotometer (Spectrostar Nano, BMG Labtech) at 30°C.
Infection Assays
Leaf infection assays were performed on B. oleracea plants (var. Bartolo) as described by Joubert et al. (2011). The experiment was repeated twice and involved at least 10 leaves inoculated by each fungal genotype.
Seed infection assays and seed contamination assessments were performed as described by Pochon et al. (2012). At least five plants per fungal genotype were inoculated and the experiment was repeated twice.
Metabolic Extraction and Amino Acids and Sugar Analysis
Metabolic extractions were made on 15 mg of dry powdered samples. Procedures for methanol–chloroform–water-based extractions as well as determination of amino acid and carbohydrate contents were described by Thouvenot et al. (2015).
Author Contributions
GN’G carried out the sample preparation and microarray analyses. RR, BI, and MM involved in the sample preparation for metabolic analysis and performed the quantitative PCR. NB-S carried out the yeast complementation assays. CM, MM, and CA-B involved in the construction and phenotyping of Alternaria brassicicola mutant strains. PS and TG conceived the study and participated in the design of the experiments as well as in the analysis of the results. RL and BP involved in sugar analyses from fungal samples. SP and J-PR involved in the gene expression analyses. CC, MM, FB, BH, AK, and BI involved in the pathological tests.
Funding
This research was supported in part by grants awarded by the Region des Pays de la Loire (FUNHY project). This research was also supported by post-doctoral fellowships (CA-B) provided by the University of Angers.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank the Metabolic Profiling and Metabolomic (P2M2) Platform from the French public institute IGEPP for assistance during metabolic phenotyping.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01969/full#supplementary-material
FIGURE S1 | Comparison of transcript levels obtained by microarray or RT-qPCR methods. The expression of 9 genes was monitored in response to sorbitol and desiccation stresses. The fold inductions (log2 values) are presented in (A) and the correlation plot in (B).
FIGURE S2 | Total organic contents in a 1 mg sample of dried fungal cultures sample exposed or not to sorbitol and desiccation treatments. (A) Total contents of amino acids. (B) Total content of saccharides and sugar alcohols. Colors display the treatments. Box plots show the median, 1st and 3rd quantiles (n = 3). ∗denotes a significant difference (P-value <0.001, Welsh test).
FIGURE S3 | Generation of null mutants by homologous recombination. (A) Schematic representation of the deleted locus in the wild-type and the replacement construct with the Hyg B resistance cassette (Hph gene). Arrows indicate the position of primers used for the cassette generation and for PCR screening of mutants. (B) Gel electrophoresis of PCR products obtained from template DNA of the wild-type, Δabsch9, Δabsih3 and Δabsih15 strains with the indicated primer pairs. Molecular sizes (kb) were estimated based on a 1 kb ladder (lane M, NEW ENGLAND BioLabs® Inc.). A supplementary PCR was done with the primer pair ITS1/ITS4 as a positive control to assess the quality of the wild-type template DNA.
TABLE S1 | List of predicted sequences from the Alternaria brassicicola genome exhibiting physicochemical features of hydrophilins and an obvious match upon Blast analysis with proteins of known functions.
TABLE S2 | List of primers used in this study.
Footnotes
- ^ https://genome.jgi.doe.gov/Altbr1/Altbr1.home.html
- ^ https://www.ncbi.nlm.nih.gov/nuccore/PHFN00000000
- ^ https://bbric-pipelines.toulouse.inra.fr/myGenomeBrowser?public=1
- ^ https://genome.jgi.doe.gov/Altbr1/Altbr1.home.html
- ^ https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE133507
References
Abba’, S., Ghignone, S., and Bonfante, P. (2006). A dehydration-inducible gene in the truffle Tuber borchii identifies a novel group of dehydrins. BMC Genomics 7:39. doi: 10.1186/1471-2164-7-39
Belmas, E., Briand, M., Kwasiborski, A., Colou, J., N’Guyen, G., Iacomi, B., et al. (2018). Genome sequence of the necrotrophic plant pathogen Alternaria brassicicola Abra43. Genome Announc. 6:e01559-17. doi: 10.1128/genomeA.01559-17
Brewster, J. L., and Gustin, M. C. (2014). Hog1: 20 years of discovery and impact. Sci. Signal. 7:re7. doi: 10.1126/scisignal.2005458
Calahan, D., Dunham, M., DeSevo, C., and Koshland, D. E. (2011). Genetic analysis of desiccation tolerance in Sachharomyces cerevisiae. Genetics 189, 507–519. doi: 10.1534/genetics.111.130369
Carrere, S., and Gouzy, J. (2017). myGenomeBrowser: building and sharing your own genome browser. Bioinformatics 33, 1255–1257. doi: 10.1093/bioinformatics/btw800
Cho, Y., Davis, J. W., Kim, K.-H., Wang, J., Sun, Q.-H., Cramer, R. A. Jr., et al. (2006). A high throughput targeted gene disruption method for Alternaria brassicicola functional genomics using linear minimal element (LME) constructs. Mol. Plant Microbe Interact. 19, 7–15. doi: 10.1094/MPMI-19-0007
Dang, H. X., Pryor, B., Peever, T., and Lawrence, C. B. (2015). The Alternaria genomes database: a comprehensive resource for a fungal genus comprised of saprophytes, plant pathogens, and allergenic species. BMC Genomics 16:239. doi: 10.1186/s12864-015-1430-7
Dang, N. X., and Hincha, D. K. (2011). Identification of two hydrophilins that contribute to the desiccation and freezing tolerance of yeast (Saccharomyces cerevisiae) cells. Cryobiology 62, 188–193. doi: 10.1016/j.cryobiol.2011.03.002
de Lima Alves, F., Stevenson, A., Baxter, E., Gillion, J. L., Hejazi, F., Hayes, S., et al. (2015). Concomitant osmotic and chaotropicity-induced stresses in Aspergillus wentii: compatible solutes determine the biotic window. Curr. Genet. 61, 457–477. doi: 10.1007/s00294-015-0496-8
Dixon, K. P., Xu, J. R., Smirnoff, N., and Talbot, N. J. (1999). Independent signaling pathways regulate cellular turgor during hyperosmotic stress and appressorium-mediated plant infection by Magnaporthe grisea. Plant Cell 11, 2045–2058. doi: 10.1105/tpc.11.10.2045
Dongo, A., Bataillé-Simoneau, N., Campion, C., Guillemette, T., Hamon, B., Iacomi-Vasilescu, B., et al. (2009). The group III two-component histidine kinase of filamentous fungi is involved in the fungicidal activity of the bacterial polyketide ambruticin. Appl. Environ. Microbiol. 75, 127–134. doi: 10.1128/AEM.00993-08
Dupont, S., Rapoport, A., Gervais, P., and Beney, L. (2014). Survival kit of Saccharomyces cerevisiae for anhydrobiosis. Appl. Microbiol. Biotechnol. 98, 8821–8834. doi: 10.1007/s00253-014-6028-5
Gadd, G. M., Chalmers, K., and Reed, R. H. (1987). The role of trehalose in dehydration resistance of Saccharomyces cerevisiae. FEMS Microbiol. Lett. 48, 249–254. doi: 10.1111/j.1574-6968.1987.tb02551.x
Garay-Arroyo, A., Colmenero-Flores, J. M., Garciarrubio, A., and Covarrubias, A. A. (2000). Highly hydrophilic proteins in prokaryotes and eukaryotes are common during conditions of water deficit. J. Biol. Chem. 275, 5668–5674. doi: 10.1074/jbc.275.8.5668
Gaucher, M., Dugé de Bernonville, T., Lohou, D., Guyot, S., Guillemette, T., Brisset, M.-N., et al. (2013). Histolocalization and physico-chemical characterization of dihydrochalcones: insight into the role of apple major flavonoids. Phytochemistry 90, 78–89. doi: 10.1016/j.phytochem.2013.02.009
Goyal, K., Walton, L. J., and Tunnacliffe, A. (2005). LEA proteins prevent protein aggregation due to water stress. Biochem. J. 388, 151–157. doi: 10.1042/BJ20041931
Gu, Q., Zhang, C., Yu, F., Yin, Y., Shim, W.-B., and Ma, Z. (2015). Protein kinase FgSch9 serves as a mediator of the target of rapamycin and high osmolarity glycerol pathways and regulates multiple stress responses and secondary metabolism in Fusarium graminearum. Environ. Microbiol. 17, 2661–2676. doi: 10.1111/1462-2920.12522
Hildebrandt, T. M., Nunes Nesi, A., Araújo, W. L., and Braun, H.-P. (2015). Amino acid catabolism in plants. Mol. Plant 8, 1563–1579. doi: 10.1016/j.molp.2015.09.005
Hohmann, S. (2015). An integrated view on a eukaryotic osmoregulation system. Curr. Genet. 61, 373–382. doi: 10.1007/s00294-015-0475-0
Huber, A., Bodenmiller, B., Uotila, A., Stahl, M., Wanka, S., Gerrits, B., et al. (2009). Characterization of the rapamycin-sensitive phosphoproteome reveals that Sch9 is a central coordinator of protein synthesis. Genes Dev. 23, 1929–1943. doi: 10.1101/gad.532109
Iacomi-Vasilescu, B., Bataille-Simoneau, N., Campion, C., Dongo, A., Laurent, E., Serandat, I., et al. (2008). Effect of null mutations in the AbNIK1 gene on saprophytic and parasitic fitness of Alternaria brassicicola isolates highly resistant to dicarboximide fungicides. Plant Pathol. 57, 937–947. doi: 10.1111/j.1365-3059.2008.01864.x
Joubert, A., Bataille-Simoneau, N., Campion, C., Guillemette, T., Hudhomme, P., Iacomi-Vasilescu, B., et al. (2011). Cell wall integrity and high osmolarity glycerol pathways are required for adaptation of Alternaria brassicicola to cell wall stress caused by brassicaceous indolic phytoalexins. Cell. Microbiol. 13, 62–80. doi: 10.1111/j.1462-5822.2010.01520.x
Joubert, A., Calmes, B., Berruyer, R., Pihet, M., Bouchara, J.-P., Simoneau, P., et al. (2010). Laser nephelometry applied in an automated microplate system to study filamentous fungus growth. Biotechniques 48, 399–404. doi: 10.2144/000113399
Kayingo, G., and Wong, B. (2005). The MAP kinase Hog1p differentially regulates stress-induced production and accumulation of glycerol and D-arabitol in Candida albicans. Microbiology 151, 2987–2999. doi: 10.1099/mic.0.28040-0
López-Martínez, G., Rodríguez-Porrata, B., Margalef-Català, M., and Cordero-Otero, R. (2012). The STF2p hydrophilin from Saccharomyces cerevisiae is required for dehydration stress tolerance. PLoS One 7:e33324. doi: 10.1371/journal.pone.0033324
Lowe, R. G. T., Lord, M., Rybak, K., Trengove, R. D., Oliver, R. P., and Solomon, P. S. (2008). A metabolomic approach to dissecting osmotic stress in the wheat pathogen Stagonospora nodorum. Fungal Genet. Biol. 45, 1479–1486. doi: 10.1016/j.fgb.2008.08.006
Ohm, R. A., Feau, N., Henrissat, B., Schoch, C. L., Horwitz, B. A., Barry, K. W., et al. (2012). Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 8:e1003037. doi: 10.1371/journal.ppat.1003037
Pascual, S., Melgarejo, P., and Magan, N. (2003). Water availability affects the growth, accumulation of compatible solutes and the viability of the biocontrol agent Epicoccum nigrum. Mycopathologia 156, 93–100.
Pascual-Ahuir, A., and Proft, M. (2007). The Sch9 kinase is a chromatin-associated transcriptional activator of osmostress-responsive genes. EMBO J. 26, 3098–3108. doi: 10.1038/sj.emboj.7601756
Pigné, S., Zykwinska, A., Janod, E., Cuenot, S., Kerkoud, M., Raulo, R., et al. (2017). A flavoprotein supports cell wall properties in the necrotrophic fungus Alternaria brassicicola. Fungal Biol. Biotechnol. 4:1. doi: 10.1186/s40694-016-0029-3
Pochon, S., Simoneau, P., Pigné, S., Balidas, S., Bataillé-Simoneau, N., Campion, C., et al. (2013). Dehydrin-like proteins in the necrotrophic fungus Alternaria brassicicola have a role in plant pathogenesis and stress response. PLoS One 8:e75143. doi: 10.1371/journal.pone.0075143
Pochon, S., Terrasson, E., Guillemette, T., Iacomi-Vasilescu, B., Georgeault, S., Juchaux, M., et al. (2012). The Arabidopsis thaliana-Alternaria brassicicola pathosystem: a model interaction for investigating seed transmission of necrotrophic fungi. Plant Methods 8:16. doi: 10.1186/1746-4811-8-16
Ravishankar, J. P., Suryanarayanan, T. S., and Muruganandam, V. (2006). Strategies for osmoregulation in the marine fungus Cirrenalia pygmea Kohl.(Hyphomycetes). Indian J. Geo Mar. Sci. 35, 351–358.
Reyes, J. L., Rodrigo, M.-J., Colmenero-Flores, J. M., Gil, J.-V., Garay-Arroyo, A., Campos, F., et al. (2005). Hydrophilins from distant organisms can protect enzymatic activities from water limitation effects in vitro. Plant Cell Environ. 28, 709–718. doi: 10.1111/j.1365-3040.2005.01317.x
Rodríguez-Porrata, B., Carmona-Gutierrez, D., Reisenbichler, A., Bauer, M., Lopez, G., Escoté, X., et al. (2012). Sip18 hydrophilin prevents yeast cell death during desiccation stress. J. Appl. Microbiol. 112, 512–525. doi: 10.1111/j.1365-2672.2011.05219.x
Sales, K., Brandt, W., Rumbak, E., and Lindsey, G. (2000). The LEA-like protein HSP 12 in Saccharomyces cerevisiae has a plasma membrane location and protects membranes against desiccation and ethanol-induced stress. Biochim. Biophys. Acta 1463, 267–278. doi: 10.1016/s0005-2736(99)00215-1
Sánchez-Fresneda, R., Guirao-Abad, J. P., Argüelles, A., González-Párraga, P., Valentín, E., and Argüelles, J.-C. (2013). Specific stress-induced storage of trehalose, glycerol and d-arabitol in response to oxidative and osmotic stress in Candida albicans. Biochem. Biophys. Res. Commun. 430, 1334–1339. doi: 10.1016/j.bbrc.2012.10.118
Sellam, A., Iacomi-Vasilescu, B., Hudhomme, P., and Simoneau, P. (2007). In vitro antifungal activity of brassinin, camalexin and two isothiocyanates against the crucifer pathogens Alternaria brassicicola and Alternaria brassicae. Plant Pathol. 56, 296–301. doi: 10.1111/j.1365-3059.2006.01497.x
Smyth, G. K., Michaud, J., and Scott, H. S. (2005). Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics 21, 2067–2075. doi: 10.1093/bioinformatics/bti270
Soanes, D. M., Alam, I., Cornell, M., Wong, H. M., Hedeler, C., Paton, N. W., et al. (2008). Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis. PLoS One 3:e2300. doi: 10.1371/journal.pone.0002300
Sweigard, J. A., Carroll, A. M., Kang, S., Farrall, L., Chumley, F. G., and Valent, B. (1995). Identification, cloning, and characterization of PWL2, a gene for host species specificity in the rice blast fungus. Plant Cell 7, 1221–1233. doi: 10.1105/tpc.7.8.1221
Tapia, H., and Koshland, D. E. (2014). Trehalose is a versatile and long-lived chaperone for desiccation tolerance. Curr. Biol. 24, 2758–2766. doi: 10.1016/j.cub.2014.10.005
Thouvenot, L., Deleu, C., Berardocco, S., Haury, J., and Thiébaut, G. (2015). Characterization of the salt stress vulnerability of three invasive freshwater plant species using a metabolic profiling approach. J. Plant Physiol. 175, 113–121. doi: 10.1016/j.jplph.2014.11.007
Tunnacliffe, A., and Wise, M. J. (2007). The continuing conundrum of the LEA proteins. Naturwissenschaften 94, 791–812. doi: 10.1007/s00114-007-0254-y
van den Bosch, F., Fraaije, B. A., van den Berg, F., and Shaw, M. W. (2010). Evolutionary bi-stability in pathogen transmission mode. Proc. Biol. Sci. 277, 1735–1742. doi: 10.1098/rspb.2009.2211
Wang, S., Liu, Y., and Xu, J. (2017). Comparison of different drying methods for recovery of mushroom DNA. Sci. Rep. 7:3008. doi: 10.1038/s41598-017-03570-7
Wartenberg, D., Vödisch, M., Kniemeyer, O., Albrecht-Eckardt, D., Scherlach, K., Winkler, R., et al. (2012). Proteome analysis of the farnesol-induced stress response in Aspergillus nidulans. The role of a putative dehydrin. J. Proteomics 75, 4038–4049. doi: 10.1016/j.jprot.2012.05.023
Winer, J., Jung, C. K., Shackel, I., and Williams, P. M. (1999). Development and validation of real-time quantitative reverse transcriptase-polymerase chain reaction for monitoring gene expression in cardiac myocytes in vitro. Anal. Biochem. 270, 41–49. doi: 10.1006/abio.1999.4085
Wong Sak Hoi, J., Beau, R., and Latgé, J.-P. (2012). A novel dehydrin-like protein from Aspergillus fumigatus regulates freezing tolerance. Fungal Genet. Biol. 49, 210–216. doi: 10.1016/j.fgb.2012.01.005
Keywords: hydrophilins, plant pathogenic fungus, seed transmission, dehydration, Alternaria brassicicola
Citation: N’Guyen GQ, Raulo R, Marchi M, Agustí-Brisach C, Iacomi B, Pelletier S, Renou J-P, Bataillé-Simoneau N, Campion C, Bastide F, Hamon B, Mouchès C, Porcheron B, Lemoine R, Kwasiborski A, Simoneau P and Guillemette T (2019) Responses to Hydric Stress in the Seed-Borne Necrotrophic Fungus Alternaria brassicicola. Front. Microbiol. 10:1969. doi: 10.3389/fmicb.2019.01969
Received: 18 September 2018; Accepted: 09 August 2019;
Published: 30 August 2019.
Edited by:
Christophe Le May, Agrocampus Ouest, FranceReviewed by:
Syama Chatterton, Agriculture and Agri-Food Canada (AAFC), CanadaAntoine Gravot, University of Rennes 1, France
Copyright © 2019 N’Guyen, Raulo, Marchi, Agustí-Brisach, Iacomi, Pelletier, Renou, Bataillé-Simoneau, Campion, Bastide, Hamon, Mouchès, Porcheron, Lemoine, Kwasiborski, Simoneau and Guillemette. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thomas Guillemette, thomas.guillemette@univ-angers.fr
†Present address: Guillaume Quang N’Guyen, Département de Biologie, PROTEO, Institut de Biologie Intégrative et des Systèmes, Université Laval, Pavillon Charles-Eugène-Marchand, Québec, QC, Canada