ABSTRACT
G protein-coupled receptors (GPCRs) are transmembrane receptors that relay signals from the external environment inside the cell, allowing an organism to adapt to its surroundings. They are known to detect a vast array of ligands, including sugars, amino acids, pheromone peptides, nitrogen sources, oxylipins, and light. Despite their prevalence in fungal genomes, very little is known about the functions of filamentous fungal GPCRs. Here we present the first full-genome assessment of fungal GPCRs through characterization of null mutants of all 15 GPCRs encoded by the aflatoxin-producing fungus Aspergillus flavus. All strains were assessed for growth, development, ability to produce aflatoxin, and response to carbon sources, nitrogen sources, stress agents, and lipids. Most GPCR mutants were aberrant in one or more response processes, possibly indicative of cross talk in downstream signaling pathways. Interestingly, the biological defects of the mutants did not correspond with assignment to established GPCR classes; this is likely due to the paucity of data for characterized fungal GPCRs. Many of the GPCR transcripts were differentially regulated under various conditions as well. The data presented here provide an extensive overview of the full set of GPCRs encoded by A. flavus and provide a framework for analysis in other fungal species.
IMPORTANCE Aspergillus flavus is an opportunistic pathogen of crops and animals, including humans, and it produces a carcinogenic toxin called aflatoxin. Because of this, A. flavus accounts for food shortages and economic losses in addition to sickness and death. Effective means of combating this pathogen are needed to mitigate its deleterious effects. G protein-coupled receptors (GPCRs) are often used as therapeutic targets due to their signal specificity, and it is estimated that half of all drugs target GPCRs. In fungi such as A. flavus, GPCRs are likely necessary for sensing the changes in the environment, including food sources, developmental signals, stress agents, and signals from other organisms. Therefore, elucidating their functions in A. flavus could identify ideal receptors against which to develop antagonists.
INTRODUCTION
Aspergillus flavus is a soilborne fungus that acts as a saprophyte, decomposing dead organic matter. A. flavus is also an opportunistic pathogen of both plants and animals. It is capable of infecting a range of crops, both pre- and postharvest, including corn, peanuts, and cottonseed, where it produces the carcinogenic mycotoxin aflatoxin (AF). Ingestion of low doses of AF can lead to hepatocellular carcinoma, while high doses are toxic and lethal. The loss of crops infected with this fungus leads to food shortages and economic losses in the billion-dollar range in the United States alone (1). Humans and animals can also be infected by A. flavus via the inhalation of spores, and it is the second leading cause of aspergillosis after Aspergillus fumigatus. It is commonly associated with manifestations of aspergillosis in the skin, oral mucosa, and subcutaneous tissues (2). In order to thrive in such vastly different environments, the fungus must be able to compete with other organisms for nutrition and defend itself against attacks from other microbes and the plant and animal hosts it invades. Therefore, the ability to sense and respond to the various environments in which it lives is critical for A. flavus’ survival.
G protein-coupled receptors (GPCRs) are transmembrane (TM) proteins that detect external signals and transmit that information inside the cell to elicit some type of response. They contain 7 TM helices that are joined by internal and external loops. The N terminus resides outside the cell, while the C terminus is inside. GPCRs bind heterotrimeric G proteins, which are composed of an α, β, and γ subunit, and the α subunit is bound to GTP or GDP and is capable of GTP hydrolysis. Upon activation by their cognate ligands, GPCRs undergo a conformational change. This triggers the exchange of GTP for GDP on the Gα subunit, resulting in dissociation of the G protein. Both the Gα subunit and the Gβγ dimer can initiate downstream signaling pathways. Additional proteins that impact this reaction are regulators of G protein signaling (RGS), which quench the signal by accelerating GTPase activity of the Gα subunit. When GTP is hydrolyzed to GDP, the G protein reassociates in its inactive heterotrimeric state. Signaling is also terminated by internalization of the GPCR, which is regulated by phosphorylation and ubiquitination of the GPCR C terminus, triggering endocytosis and degradation of the receptor (reviewed in reference 3).
Both the Gα and Gβγ subunits can act as effectors of downstream signaling cascades. These have been well studied in the model yeast Saccharomyces cerevisiae. Sugars stimulate a cAMP-protein kinase A (PKA) pathway in which Gα activates adenylate cyclase, which converts ATP to cAMP. The increased concentration of cAMP results in the release of inhibitory subunits from PKA, thereby activating the enzyme. PKA then phosphorylates downstream proteins to impact the cellular response. Following yeast pheromone sensing, the Gβγ dimer acts as the effector for the mitogen-activated protein kinase (MAPK) pathway, which consists of a kinase cascade that regulates mating (4). Components of both the cAMP-PKA and MAPK pathways have been identified in the aspergilli and found to regulate germination, sporulation, mycotoxin production, and stress tolerance (reviewed in references 5 and 6).
Based on their potential importance to fungal growth and survival, GPCRs have been the subject of numerous bioinformatics studies. As a result, the entire set of GPCRs encoded by various fungi has been predicted for Aspergillus nidulans, Aspergillus fumigatus, Aspergillus oryzae, Magnaporthe grisea, Cryptococcus neoformans, Neurospora crassa, Verticillium spp., and Trichoderma spp. (7–12), with numbers ranging from 7 (C. neoformans) to 22 (Verticillium dahliae), in contrast to the three GPCRs in the S. cerevisiae genome. However, the functions of almost all of these GPCRs remain unknown. In this study, we identified the putative GPCRs encoded by A. flavus and carried out functional studies to investigate the roles of these GPCRs in growth, development, response to various carbon and nitrogen sources, stress agents, and fatty acid-derived signals. From this, we can begin to classify A. flavus’ GPCRs by their biological functions. We found that many receptors appear to be involved in multiple processes, highlighting the cross talk that occurs between different signal transduction pathways.
RESULTS
A. flavus encodes 15 putative GPCRs.We began our search for A. flavus GPCRs by examining the set of 16 GPCRs reported for A. nidulans (7). In an effort to focus on the canonical GPCRs, our study did not include the PTH11-like membrane-spanning proteins, of which A. nidulans is predicted to encode at least 25 (7, 8). A BLAST search of the amino acid sequences of each of the 16 A. nidulans GPCRs against the A. flavus genome yielded 77 total hits, 57 of which were unique. GPCRs pass the membrane seven times, so two protein topology prediction tools, TMHMM 2.0 (13, 14) and TopPred 1.10 (15, 16), were used to determine which of these 57 proteins were likely GPCRs. Twenty-nine of the fifty-seven sequences were not predicted to encode any transmembrane (TM) domains, so they were eliminated. Of the remaining 28 proteins, 13 were predicted to have 7 TM domains by at least one of the prediction tools. Two additional proteins that were homologous to A. nidulans GprF and GprO were predicted to have fewer than 7 TM domains, but since their A. nidulans orthologs were predicted to have 7 TM domains, they were also included with the rest of the A. flavus GPCRs, bringing the total to 15 (Table 1).
Overview of A. flavus GPCRs
We also compared the A. flavus genome to the reported A. oryzae GPCRs (7) since these two species are believed to be the same species, with A. oryzae constituting a domesticated clade (17). Similar to what was reported for A. oryzae, no orthologs for A. nidulans GprE, GprI, or GprN were found in the A. flavus genome. A novel A. oryzae GPCR not found in A. nidulans was identified and named GprQ. Its A. flavus ortholog (AFLA_132040; 99% identity, with an E value of 7e−115) had only four (TMHMM) or three (TopPred) predicted TM domains. Since A. oryzae GprQ was reported to encode only five predicted TM domains, the A. flavus homolog was not included in the set of GPCRs.
The BLAST search also identified two GPCRs that were not previously identified in the aspergilli, GprR and GprS. Like GprK, GprR harbors an RGS (regulator of G protein signaling) domain at the C-terminal end. This type of GPCR was first discovered in Arabidopsis thaliana and called A. thaliana RGS1 (AtRGS1). AtRGS1 regulates cell proliferation as well as sensitivity to sugars via its functional RGS domain. Unlike canonical GPCRs, it does not trigger the exchange of GTP for GDP on a G protein. Upon sensing glucose, AtRGS1 interacts with the constitutively active Gα subunit A. thaliana GPA1 (AtGPA1), resulting in hydrolysis of GTP and subsequent deactivation of the G protein (18–20). GPCR-RGS hybrids have since been found in several species of filamentous ascomycete fungi in addition to the aspergilli, including N. crassa, M. grisea, Fusarium graminearum, and several Trichoderma species (7, 10, 12). GprS contains a PQ loop repeat, grouping it with the class IV receptors GprF, GprG, and GprJ. The function of the PQ loop repeat domain is unknown, although this class contains the Schizosaccharomyces pombe nitrogen starvation sensor Stm1 (7).
Based on their domain structures, A. flavus’ 15 putative GPCRs were assigned to nine different classes. These classes have been previously described (7, 21), and a few notes on each class are included in Table 1. The amino acid sequences for all 15 GPCRs, as well as the sequences for GPCRs from A. nidulans, A. fumigatus, A. oryzae, F. graminearum, M. grisea, N. crassa, Trichoderma reesei, and V. dahliae, were aligned using the Clustal Omega software tool (22), and a phylogenetic tree was generated using the Phylogeny.fr software tool (23) (see Fig. S1 in the supplemental material). From this analysis, it was apparent that the A. flavus GPCRs GprF, GprH, GprS, and GprR are not orthologous to any of A. oryzae’s predicted GPCRs, despite the high degree of genetic similarity between the two fungi. Interestingly, of the two GPCRs that were not previously found in the aspergilli, GprS was most closely related to two GPCRs from V. dahliae, while GprR was unique among the fungi analyzed here.
Figure S1
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
In order to verify the coding sequences of the A. flavus GPCRs, 13 of the 15 GPCRs were amplified, cloned, and sequenced from A. flavus mRNA. cDNA for gprJ and gprM could not be amplified. Two GPCRs, gprF and gprO, appeared to be misannotated. In both cases, the actual sequence improved the proteins’ alignments with their A. nidulans orthologs, and for GprO, it changed the number of predicted TM domains from six to seven (TopPred), with the updated prediction reported in Table 1. GprF was only predicted to encode four (TMHMM) or five (TopPred) TM domains, but since A. nidulans GprF was predicted to encode seven TM domains, it was included with the set of A. flavus GPCRs.
Individual GPCRs deleted in A. flavus.All 15 of the GPCRs were individually disrupted in the CA14 ΔpyrG Δku70 strain of A. flavus by replacing the GPCR-encoding gene with pyrG. Deletion mutants were confirmed by PCR (data not shown) and Southern blotting (see Fig. S2 in the supplemental material), and confirmed isolates were checked for marker gene effects by growing on media with and without uridine and uracil (UU) (see >Fig. S3 to S5; also data not shown) (24). The ΔgprG, ΔgprK, ΔgprM, and ΔgprR mutants all exhibited growth defects that were at least partially remediated by supplementation of UU.
Figure S2
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
Figure S3
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
To investigate the impact of UU supplementation on a phenotypic assay, a subset of mutants was examined in greater detail. First, gprR was disrupted in strain NRRL3357 of A. flavus. Like the CA14 ΔgprR mutant, the NRRL3357 ΔgprR mutant was severely restricted in growth in the absence of UU, demonstrating that the need for pyrimidine supplementation was conserved in different genetic backgrounds (see Fig. S3 in the supplemental material). Strains were then exposed to two different stressors: alkaline pH (pH 8) and cell wall stress (via Congo red), with and without UU. The wild types for both the CA14 and NRRL3357 backgrounds were inhibited to similar degrees regardless of UU supplementation, whereas the addition of UU was required to observe stress-induced inhibition for both ΔgprR mutants (see Fig. S3). Multiple independent isolates of a mutant in the same genetic background were similarly examined. The ΔgprM mutant had shown a marker gene effect for all three confirmed transformants, and without UU, the effects of the stresses were more dramatic (see Fig. S4). As a final test of the effects of UU supplementation, multiple isolates of two mutants that had not shown marker gene effects, the ΔgprD and ΔgprH mutants, were analyzed. Independent isolates of both mutations showed similar patterns of inhibition regardless of UU supplementation (see Fig. S5). The impacts of alkaline pH and cell wall stress were quantified for all of the isolates of the ΔgprR, ΔgprM, ΔgprD, and ΔgprH mutants, revealing that multiple isolates of the same mutant behaved similarly (see Table S1). Therefore, a single isolate for each GPCR disruption was chosen, and UU was included in the media for all subsequent experiments.
Figure S4
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
Figure S5
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
Table S1
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
GPCRs important in basic A. flavus development and aflatoxin synthesis.Various fungal developmental processes have previously been linked to G protein signaling. To elucidate which GPCRs may be involved in mediating these signals, we measured the germination rate, AF biosynthesis, sporulation under both light and dark conditions, and production of spores and sclerotia at high and low cell densities, described below.
(i) Germination.Germination studies in A. nidulans have revealed several components of G protein signaling to be important in that process. The heterotrimeric G protein GanB-SfaD-GpgA activates a cAMP pathway in response to glucose to drive germination, and this is abrogated by RGS RgsA (25–27). There is also evidence of a small GTPase, RasA, controlling germination independently of cAMP (25, 28). The only GPCRs implicated in germination thus far are A. nidulans GprD (21) and A. fumigatus GprC and GprD (29). Therefore, germination of the entire set of A. flavus Δgpr mutants was monitored from 4 to 9 h postinoculation. During this time period, the wild type went from 1% to 77% germinated spores (data not shown). The germination rates of the Δgpr mutants were measured, and a heat map was generated with the percentage of mutant germination versus wild-type germination for each time point (Table 2). Aside from a slightly enhanced germination rate at 6 h by the ΔgprO mutant, the rest of the mutants were either equivalent or impaired in germination compared to the wild type. The most severely deficient mutants were the ΔgprA and ΔgprK strains, with <50% of the wild-type germination rate at nearly all time points measured. Many of the other mutants, including the ΔgprB, ΔgprC, ΔgprH, ΔgprJ, ΔgprP, ΔgprR, and ΔgprS strains, were impaired in germination as well.
Developmental phenotypes of A. flavus Δgpr mutantsa
(ii) Light-induced sporulation.Rhodopsin is well characterized as a receptor of light, and A. flavus encodes one GPCR, NopA, that belongs to the class IX bacteriorhodopsin-like GPCR family. Since light is an inducer of sporulation, the fungi were exposed to continuous light or continuous dark, and their spores were counted. The wild type and most of the mutants produced significantly more spores in the light versus production in the dark. Two mutants, the ΔgprF and ΔgprR mutants, showed no significant difference between the two conditions. Notably, the ΔnopA mutant exhibited the wild-type response to light (Table 2).
(iii) Quorum sensing.A. flavus undergoes a quorum-sensing-mediated developmental shift in which at low inoculum density, it produces many sclerotia and very few spores. The opposite is observed at high inoculum density (30). Because an exogenous signal mediates this phenomenon, GPCRs are plausible targets for relaying the signal. In fact, two GPCRs in the A. flavus strain NRRL3357, encoded by gprC and gprD, have already been implicated as quorum-sensing receptors based on their inability to respond to a high-density extract (31). Surprisingly, all 15 Δgpr mutants produced more sclerotia and fewer spores at high density than the wild type (data not shown). However, the trend of decreasing sclerotia and increasing spore production as inoculum density increased was still observed for most of the mutants, with the exceptions of the ΔgprB and ΔgprF mutants There was no significant difference in the amount of spores produced at high and low densities for the ΔgprB mutant, nor was there a significant difference in sclerotia production at the two densities for the ΔgprF mutant (Table 2).
(iv) Aflatoxin.G protein signaling has been shown to regulate synthesis of AF and its precursor, sterigmatocystin (ST). The Gα subunit FadA and its RGS FlbA are part of an adenylate cyclase/cAMP/protein kinase A (PKA) (PkaA) pathway that controls ST production in A. nidulans via transcriptional and posttranscriptional regulation of aflR (32–34). The orthologous pathway regulates AF biosynthesis in Aspergillus parasiticus (35) and A. flavus (J. W. Bok and N. P. Keller, unpublished). Nothing is known about which GPCR(s) might be initiating this pathway, so the panel of Δgpr mutants was tested for their ability to make AF. AF was extracted and measured following growth on AF-inducing and AF-repressing media (YES [2% yeast extract and 6% sucrose, pH 5.8] and YEP [2% yeast extract and 6% peptone, pH 5.8], respectively). No mutant gained the ability to produce detectable AF on YEP. However, on YES medium, the ΔgprA mutant made ~1.5× more AF than the wild type, and the ΔgprP mutant made ~2.3× more AF than the wild type (Table 2). None of the other mutants were affected in their ability to produce AF.
GPCRs impact growth on various carbon and nitrogen sources.In addition to developmental cues, GPCRs are able to detect the presence or absence of important nutrients. GprC and GprD belong to class III, which is the group that also contains the sugar receptors S. cerevisiae Gpr1 and S. pombe Git3. Gpr1 senses glucose and sucrose (with a stronger affinity for sucrose) and transmits the signal through a G protein/cAMP/PKA pathway (36–38). Other Gpr1 homologs are also carbon source sensors, such as S. pombe Git3 (reviewed in reference 39) and N. crassa GPR-4 (40), although the role of C. albicans Gpr1 in carbon sensing is ambiguous (41, 42). While a carbon-sensing GPCR has not been identified in the aspergilli, it has been shown that carbon sensing is mediated by the heterotrimeric G protein GanB-SfaD-GpgA in A. nidulans (27). To identify A. flavus GPCRs involved in carbon sensing, strains were grown on a variety of carbon and nitrogen sources. These included glucose, galactose, xylose, sucrose, and corn oil for carbon sources. Corn oil is mainly composed of C18 and C16 fatty acids (43). The ΔgprA, ΔgprC, ΔgprJ, ΔgprK, and ΔgprR mutants were impaired in growth on several of the carbon sources tested (Table 3). The ΔgprR mutant was restricted in growth on every carbon source tested except galactose and corn oil.
Growth of A. flavus Δgpr mutants on various carbon and nitrogen sourcesa
Class IV receptors are grouped together based on their relatedness to the S. pombe nitrogen starvation sensor Stm1, and A. flavus has four receptors in this group. Nitrogen starvation results in cell cycle arrest in S. pombe, which allows the yeast to remain viable either until nutrition becomes available or mating partners are present. stm1 is expressed at low, moderate, and high levels on peptone, ammonium, and proline, respectively. Deletion of stm1 causes premature cell cycle arrest under nitrogen starvation, while overexpression of stm1 represses vegetative growth under nitrogen-replete conditions (44). Therefore, the mutants were tested for their ability to grow on peptone, ammonium chloride, and proline. While no mutant was restricted on peptone, mutants carrying both ΔgprC and ΔgprD deletions—the class III receptors—were inhibited on proline. The ΔgprR mutant was impaired on ammonium chloride and proline (Table 3).
GPCRs contribute to fungal stress responses.Fungi transmit the myriad stress signals they encounter through several conserved, interconnected signal transduction pathways. Multiple stimuli trigger each pathway, and the pathways themselves engage in cross talk with one another. For example, the cell wall integrity (CWI) MAPK pathway is triggered in response not only to cell wall stress but to osmoregulatory and oxidative stresses as well, among others (reviewed in reference 45). In addition to the CWI pathway, osmotic stress is also managed by another MAPK pathway, the high-osmolarity glycerol (HOG) pathway (reviewed in reference 46). pH tolerance is maintained by the Pal/Rim alkaline response pathway. It transmits its signal via PalH, a 7-TM-domain protein that interacts with an arrestin protein to initiate endocytosis, ultimately activating the PacC transcription factor (reviewed in reference 47). Besides their overlapping roles in stress tolerance, what these pathways have in common is the prevalence of G protein signaling components. To examine whether any of A. flavus’ GPCRs mediate stress responses, the mutants were exposed to a variety of stressors that inhibit the radial growth of the wild type. The percent inhibition of growth in the presence versus absence of the stressor was compared between the wild type and each mutant, and the changes in inhibition between wild type and the mutants are presented as a heat map (Table 4).
Stress responses of A. flavus Δgpr mutantsa
Hydrogen peroxide was used at a concentration of 3 mM to assay oxidative stress. Loss of gprH rendered the fungus more resistant to hydrogen peroxide, while the ΔgprG mutant was more sensitive. Cell wall stress was created by growing the fungi on Congo red, which weakens the cell wall by binding to growing chitin chains, preventing their attachment to other cell wall components (48, 49). The ΔgprK, and ΔgprS mutants were both more sensitive to Congo red. The impact of osmotic pressure was measured with 1 M sodium chloride. Three mutants, carrying the ΔgprK, ΔgprM, and ΔgprR deletions, were more sensitive than the wild type to the hyperosmotic condition. Finally, the strains were grown on media buffered to pH 4 and pH 8. The ΔgprM mutant was more sensitive to both conditions, while the ΔgprK and ΔgprR mutants were more sensitive to acidic and alkaline pHs, respectively. Alternatively, the ΔgprD, ΔgprF, and ΔgprG mutants were more resistant to acidic pH.
GPCRs are involved in response to lipids and oxylipins.It has long been speculated that oxylipins are detected by fungal GPCRs, since this is their mode of perception in mammalian cells (reviewed in references 3, 5, and 50). To address this, strains were exposed to the oxylipins methyl jasmonate (MeJA) and 13-hydroperoxy-octadecadienoic acid [13(S)-HpODE] and to the fatty acid linoleic acid. Any mutant that deviated from the wild-type response to these chemicals is indicated in Table 5 with a shaded box.
A. flavus Δgpr mutant responses to fatty acids and oxylipinsa
The effect of MeJA on AF production varies depending on the experimental setup and concentration of MeJA used (51). Under certain conditions, MeJA has been shown to repress AF production by A. flavus (52), but under others, it has enhanced AF synthesis by A. parasiticus (53). For this study, an AF-repressive protocol was followed (51) in which cultures were exposed to 10−4 M MeJA, and AF was measured after 3 days. AF repression was observed for the wild type and five of the Δgpr mutants: the ΔgprC, ΔgprD, ΔgprH, ΔgprJ, and ΔnopA mutants. The rest of the mutants showed no change in AF production except for the ΔgprK mutant, which actually made more AF when exposed to MeJA (Table 5). A lack of response to MeJA cannot be attributed to a basic biosynthetic defect because, as previously shown, none of the mutants were impaired in AF biosynthesis on YES media (Table 2).
13(S)-HpODE, as well as its precursor fatty acid, linoleic acid, has been shown to induce sporulation in A. flavus (54). Stimulation of conidiation by 13(S)-HpODE was lost for the ΔgprA, ΔgprC, ΔgprD, ΔgprF, ΔgprG, ΔgprJ, ΔgprO, ΔgprP, and ΔgprR mutants. To test whether this represented a general sporulation defect or a specific lack of response to 13(S)-HpODE, sporulation was also monitored in response to linoleic acid. All of the mutants except the ΔgprR mutant responded as did the wild type to linoleic acid, indicating that the lack of response to 13(S)-HpODE was specific for all of the above mutants except the ΔgprR mutant (Table 5).
GPCRs are differentially expressed under disparate conditions.To further understand potential roles of A. flavus GPCRs, their expression levels under different conditions from various published studies were investigated. The results are compiled in Table 6, which revealed several GPCRs to be differentially expressed on maize versus liquid culture mycelia, and these GPCRs do not overlap with those that are differentially expressed on maize versus wheat (55). The corn kernel germ is mainly composed of lipids, while the endosperm is predominantly starch, and A. flavus preferentially grows in the lipid-rich portions of the kernel (56). Interestingly, four GPCRs, encoded by gprC, gprH, gprM, and gprS, were differentially expressed on the germ versus the endosperm (55). Two separate studies compared genomic expression from mycelial growth versus sclerotial production, though the methods used were quite different between the two studies (55, 57). The strains of A. flavus, media, and growth conditions (i.e., liquid versus solid culture) that were used varied, as did the methods for measuring gene expression (i.e., microarray versus transcriptome sequencing [RNA-seq]). Both studies showed that many of A. flavus’ GPCRs were differentially expressed depending on whether the fungus was growing as mycelia or producing sclerotia, although the expression patterns differed greatly between the two studies. 5-Azacytidine (5-AC) is an inhibitor of AF production and sporulation (58), but application of 5-AC to cultures of A. flavus did not impact expression of any of its GPCRs (59).
Expression of GPCR-encoding genes under various conditionsa
DISCUSSION
Several works have bioinformatically identified the GPCRs encoded by various fungi. These include A. nidulans, A. fumigatus, A. oryzae, M. grisea, C. neoformans, N. crassa, Verticillium spp., and Trichoderma spp. (7–12). Functions have been assigned for a few of these GPCRs. However, this is the first study to identify functions for a fungus’ entire set of GPCRs, providing leads on most of the GPCRs that can be pursued further in future studies.
An attempt was made to identify the activating ligands of A. flavus’ GPCRs through heterologous expression in S. cerevisiae. Despite testing more than 60 potential ligands, no ligand-receptor pair was elucidated (data not shown). Individually deleting each GPCR, however, led to an abundance of findings that highlight the complexity of signal transduction networks. Of the 15 deletion mutants, all but two (the ΔgprO and ΔnopA mutants) exhibited aberrancies in multiple processes. This was not unexpected, since there is overlap in many of the pathways that drive various aspects of fungal development. For example, germination in A. nidulans is triggered by a carbon source (28), so without the receptor(s) that detects that signal, the fungus would likely be impaired in both germination and growth on various carbons sources, as was the case with the ΔgprA, ΔgprC, ΔgprJ, ΔgprK, ΔgprP, and ΔgprR mutants.
Stress pathways are also interconnected among themselves and with other developmental pathways, making it difficult to decipher which GPCRs are tied to which responses. Nitrogen starvation, for example, creates a stressful environment for cells. Chung et al. (44) found that the nitrogen starvation sensor Stm1 was important for growth on different nitrogen sources and that it also appeared to stimulate a stress response pathway. These are just a couple examples of cross talk between different signaling and development pathways, though many more such instances exist, creating a domino effect when a component of one pathway is disrupted.
Table 7 lists each GPCR’s predicted function based on its sequence alignment as well as the various processes it may be involved in based on the functional studies reported here. Interestingly, none of the predicted functions (Table 1) were validated by the experimental results as presented here, which aligns with the paucity of data on the functions of filamentous fungus GPCRs. In the case of the yeast pheromone receptor Ste2 and Ste3 homologs GprA and GprB, a role in mating was not evaluated in this study. However, in A. nidulans, gprA and gprB are required for fruiting-body formation during self-fertilization (60). A. flavus is a heterothallic fungus, and it undergoes sexual reproduction when the appropriate mating type is present (61). Therefore, GprA and GprB are likely to be important for detecting the presence of mating partners and triggering sexual reproduction.
Potential roles of A. flavus GPCRs
GprC and GprD were predicted to be carbon source sensors based on their homology with the S. pombe glucose receptor Git3. However, only GprC displayed any growth defects on different carbon sources, suggesting that despite their high degree of sequence similarity, GprC and GprD have distinct roles. In addition to GprC, several other putative GPCRs may be involved in carbon source sensing. Deletion of gprA, gprJ, gprK, and gprR also resulted in impaired growth on multiple carbon sources. It was shown in S. cerevisiae that its carbon sensor Gpr1 was partially redundant with Ras2, a small GTPase, such that a Δgpr1 mutant grew normally but a Δgpr1 Δras2 strain had a severe growth defect (36). Perhaps several of A. flavus’ receptors have overlapping functions in carbon sensing with each other and/or with other proteins, such as small GTPases, which would explain why all of the single GPCR mutants were still viable. It is also possible that several of the GPCRs work as heterodimers, as is found in other systems (62), and thus their full role in fungal development will be uncovered only in multiple deletion strains.
Class IV GPCRs are S. pombe Stm1 homologs and are therefore predicted to have a role in nitrogen sensing. A. flavus’ four class IV GPCRs are GprF, GprG, GprJ, and GprS. However, none of these mutants showed any deviation from the wild-type phenotype on the three nitrogen sources tested. Instead, the ΔgprC, ΔgprD, and ΔgprR mutants were restricted on proline, with the ΔgprR mutant also showing restriction on ammonium chloride. This is reminiscent of the S. pombe Δstm1 mutant showing premature cell cycle arrest on proline, a poor nitrogen source (44). Further experimentation is required to definitively classify these receptors as nitrogen sensors, though it is interesting that GprC and GprD may have overlapping roles in nitrogen sensing and divergent roles in carbon sensing.
Not much is known about class V receptors. They are predicted to be cAMP sensors, though there is no evidence that this occurs in Aspergillus. A. flavus encodes only one class V receptor, GprH. It is homologous to C. neoformans Gpr4, which is a methionine sensor. The Gpr4 protein is internalized in response to methionine, and a Δgpr4 mutant does not produce a spike in cAMP in the presence of methionine as is seen in the wild type. The Δgpr4 mutant is defective in capsule production and mating (9). Heterologous expression of A. flavus GprH in yeast revealed it to be active in the absence of any added ligand (data not shown), although efforts to determine whether this was due to methionine in the media were inconclusive.
NopA belongs to the bacteriorhodopsin-like class of GPCRs (class IX). Rhodopsins consist of an opsin protein and a retinal chromophore and act as light-responsive GPCRs or proton pumps (reviewed in reference 63). In spite of the wealth of knowledge about rhodopsins, the function of many opsins, including NopA, remains enigmatic. The first opsin identified in fungi was N. crassa Nop-1 (64), which binds retinal (65). It has no obvious function, though its expression is induced during conidiation, which is a light-responsive process (66). The Fusarium fujikuroi Nop-1 ortholog also lacks an obvious function, although it appears to play a role in the regulation of retinal biosynthesis (67).
Based on its sequence, A. flavus NopA is an opsin-related protein instead of a rhodopsin, meaning that it lacks a critical lysine residue required for binding retinal. Opsin-related proteins are found throughout the fungal kingdom, though nothing is known about what cofactors, if any, they bind to and what their functions may be (68). Other than a slight germination defect, the ΔnopA mutant did not appear different from the wild type in any of the development, growth, stress, or lipid response assays described here. Either NopA is functionally redundant with another receptor, or else it has a role that was not addressed by these experiments.
Five GPCRs from three different classes (classes VI, VII, and VIII) had no predicted functions. GprK and GprR both belong to class VI, characterized by hybrid GPCR-RGS proteins. Although the ΔgprK and ΔgprR mutants shared some phenotypes, the ΔgprR mutant was distinct for its defective growth/response in nearly every assay performed. This suggests that GprR might regulate a process that is required for growth, and without it, the mutant is impaired in its response to many other stimuli. GprO and GprP are both class VII receptors. They may share overlapping functions in oxylipin sensing, but GprP clearly has additional roles in germination and AF production.
In this study, the Δgpr mutants were evaluated for their basic growth and development, growth on a variety of carbon and nitrogen sources, stress responses, and responses to lipids and oxylipins. An area that was not explored but is likely to be extremely important for A. flavus is the roles of its GPCRs during interaction with other organisms. A. flavus is naturally found in the soil, where it encounters many different microbes, and it is also capable of infecting plants and humans. These interactions are mediated by signaling between organisms, and it is quite possible that some of these signals are transduced via A. flavus’ GPCRs.
GPCRs are an attractive antifungal target for their lock-and-key specificity. In fact, approximately half of all drugs target GPCRs (69). Finding antagonists for GPCRs that are required for fungal invasion of a host, for example, could be a very effective way of combating infection. No single GPCR was identified here as being absolutely required for growth in culture, but perhaps certain receptors are required to grow or produce toxins in a host and/or function as heterodimers. By unraveling the functions of A. flavus’ GPCRs, new therapeutic targets may emerge for reducing the impact of this pathogen.
MATERIALS AND METHODS
Culture conditions.Strains (see Table S2 in the supplemental material) were grown on glucose minimal medium (GMM) (33) with ammonium salts instead of nitrate salts unless otherwise mentioned. In some cases, 5 mM uridine and 5 mM uracil were added, and this is denoted as “plus UU.” For genomic DNA (gDNA) extraction, strains were grown on liquid GMM (NH4+) with 0.5% yeast extract added.
Table S2
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
Strain construction.Strain genotypes and sources are summarized in Table S2 in the supplemental material. All primers used for strain construction are listed in Table S3. GPCRs were disrupted using homologous recombination to replace each gene with pyrG in the parental strain CA14Δku70ΔpyrG (70). Table S3 indicates which primers were used to amplify the GPCR 5′ and 3′ flanks. The pyrG marker gene was amplified from A. fumigatus genomic DNA using primers 61 and 62. In all cases, the 5′ flank forward primer and 3′ flank reverse primer were used to amplify the entire double-joint cassette. CA14Δku70ΔpyrG was also complemented to prototrophy to generate an isogenic “wild type” for comparison with the other mutants. To achieve this, A. fumigatus pyrG was targeted to the nkuA locus. Primers 63 and 64 and primers 65 and 66 were used to amplify the nkuA 5′ and 3′ flanks, respectively. The entire construct was amplified with primers 63 and 66.
Table S3
Copyright © 2014 Affeldt et al.This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.
All of the above constructs were transformed individually into parental strain CA14Δku70ΔpyrG (70), with the exception of ΔgprR, which was also transformed into strain 3357.5 (71). Transformation of the fungus was carried out as described previously (31). At least 12 independent isolates were screened by PCR, followed by Southern analysis to confirm PCR-positive mutants. The same primers that were used to amplify the full double-joint cassette were used to amplify Southern probes from wild-type gDNA. Confirmed mutants were grown with and without UU to check for marker gene effects (24), and a single isolate of each mutant was selected for subsequent analyses.
Bioinformatic characterization of GPCRs.The 15 putative A. flavus GPCRs were identified by a BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) search of the A. nidulans GPCRs (7) against the A. flavus genome. The topology prediction software programs TMHMM 2.0 (13, 14) and TopPred 1.10 (15, 16) were used to predict the number of TM domains for all unique hits. GPCRs were aligned using the software program Clustal Omega (22), and a phylogenetic tree was generated using Phylogeny.fr (23).
Fungal development assays. (i) Germination assay.Germination was assayed using a 48-well plate, with 400 µl of GMM(NH4+) plus UU plus 0.1% yeast extract containing 106 spores/ml dispensed per well (three replicates per strain). The plate was incubated at 29°C, and germination was measured at 4, 5, 6, 7, 8, and 9 h postinoculation. The plate was imaged on an Eclipse TS100 inverted microscope (Nikon) with an EOS Rebel T3 camera (Canon) fitted into the eyepiece. The Cell Counter tool in the ImageJ program (72) was used to count the number of germinated spores per 100 spores. A spore was counted as germinated when the germ tube extended for at least the length of the spore.
(ii) Aflatoxin analysis.Strains were grown on 10-cm-diameter plates containing 25 ml YES (2% yeast extract and 6% sucrose, pH 5.8) or YEP (2% yeast extract and 6% peptone, pH 5.8) with uridine and uracil (plus UU) and 1.5% agar. Spores (103) in 1 µl of 0.1% Tween 20 were point inoculated onto the center of each plate. Three replicates per strain were placed in the dark at 30°C for 5 days. A 15-mm-diameter core was punched from each plate and homogenized in 3 ml 0.01% Tween 20.
To extract AF, 3 ml of ethyl acetate was added to each tube, and the tubes were shaken vigorously and spun at 3,000 rpm for 10 min. The organic layer was removed, dried down, and resuspended in 1 ml 50:40:10 water-methanol-acetonitrile. Samples were filtered through an Acrodisc syringe filter with a nylon membrane (0.45 µm; Pall Corporation) and run on a PerkinElmer Flexar instrument equipped with a Zorbax Eclipse XDB C18 column (150 by 4.6 mm, 5 μm, 100 Å; Agilent). The column was equilibrated in the running solvent (50:40:10 water-methanol-acetonitrile), and 20 µl was injected and run isocratically for 11 min with 100% running solvent at a flow rate of 0.8 ml/min. AF was detected by a fluorescent detector with an emission wavelength of 455 nm and excitation wavelength of 365 nm.
(iii) Impact of light and dark on sporulation.Strains were grown on 10-cm plates containing 10 ml GMM(NH4+) plus UU (15 g/liter agar). Five hundred spores in 1 µl 0.1% Tween 20 were point inoculated into the center of each plate, and plates were incubated for 3 days at 30°C in either constant light or constant dark. There were three replicates per strain per condition. Following incubation, a 15-mm core was punched from the center of each plate and homogenized in 3 ml 0.1% Tween 20. A 100-µl aliquot was removed and diluted, and spores were enumerated using a hemocytometer.
(iv) Density dependence assay.The impact of GPCRs on density dependence was assayed with the method of Horowitz Brown et al. (30) with some modifications. A 35.3 mM concentration of glutamine (Gln) was used as a nitrogen source since NH4+ does not support the normal density-dependent development pattern. Ten-centimeter plates containing 10 ml GMM(Gln) plus UU plus 2% sorbitol (15 g/liter agar) were overlaid with 103 (low density) or 107 (high density) spores in 3 ml of the base medium containing 7.5 g/liter agar. Plates were incubated in the dark at 30°C for 6 days. Following incubation, three replicates were processed for sclerotia production, and the other three replicates were processed for spore production.
To measure sclerotia, plates were sprayed with 70% ethanol to wash off conidia. The sclerotia were scraped from the plates, frozen in liquid nitrogen, and lyophilized. The dry weight was recorded. To measure spore production, three 15-mm cores were punched from each plate and homogenized in 3 ml 0.1% Tween 20. A 100-µl aliquot was removed, diluted, and counted on a hemocytometer.
Growth on alternative carbon and nitrogen sources.To assess the effects of a variety of carbon sources, strains were grown on a variety of media, the base being MM(NH4+) plus UU (15 g/liter agar) as described above. The various carbon sources were added at the following concentrations: glucose, galactose, and xylose, 1%; sucrose, 0.95%. To test alternative nitrogen sources, GMM plus UU (15 g/liter agar) with nitrogen-free salts added was the base medium. Three grams per liter peptone or 8.13 g/liter proline was added, and these media were compared to GMM(NH4+) plus UU (15 g/liter agar). Forty milliliters of medium was dispensed into square 10-cm by 10-cm plates and inoculated with 500 spores of each strain in 1 µl 0.1% Tween 20. Plates were incubated at 30°C for 3 days with three replicate plates per condition.
Fungal stress assays.Several different media were used to assess the roles of GPCRs in stress responses. The base for these media was GMM(NH4+) plus UU (15 g/liter agar). The following stressors were added to the media after autoclaving: hydrogen peroxide (3 mM), Congo red (200 µg/ml), and sodium chloride (1 M). Media were also buffered to two different pHs before autoclaving: pH 4 (100 mM NaH2PO4 and 100 mM NaCl) and pH 8 (100 mM Na2HPO4). The pH stress media contained 30 g/liter agar, so they were compared to GMM(NH4+) plus UU with 30 g/liter agar. Forty milliliters of medium was dispensed into square 10-cm by 10-cm plates and inoculated with 500 spores of each strain in 1 µl 0.1% Tween 20. Plates were incubated at 30°C for 3 days with three replicate plates per condition.
Fatty acid assays. (i) Impact of methyl jasmonate on aflatoxin production.To test the impact of methyl jasmonate (MeJA) on the ability of the GPCR mutants to synthesize aflatoxin, the method of Meimaroglou et al. (51) was followed with some modifications. One hundred spores were inoculated into 10-cm plates containing 10 ml YES (2% yeast extract and 6% sucrose, pH 5.8) plus UU liquid medium with or without 0.1 mM MeJA (Sigma-Aldrich). Three replicates per condition per strain were incubated in the dark at 30°C for 3 days. The contents of each plate were transferred to a 50-ml conical tube and homogenized.
To extract aflatoxin, 5 ml ethyl acetate was added to each tube for 15 min with periodic shaking. The mixtures were spun at 3,000 rpm for 5 min, and 3 ml of the organic layer was transferred to a glass vial, dried down, and resuspended in 1 ml 50:40:10 water-methanol-acetonitrile. They were analyzed by high-performance liquid chromatography (HPLC) as described above.
(ii) Impact of fatty acids on sporulation.Disk assays were carried out to measure the impact of linoleic acid and 13(S)-HpODE (Cayman Chemical) on sporulation, following the method of Calvo et al. (54) with some modifications. Ten-centimeter plates containing 10 ml YGT (0.5% yeast extract, 2% glucose, and 1 ml/liter trace elements) plus UU with 1.5% agar were overlaid with 104 spores per plate in 3 ml YGT plus UU with 0.75% agar. One-centimeter filter disks (Whatman, grade 1) were soaked with either 1 mg linoleic acid in 8 µl ethanol, 1 µg 13(S)-HpODE in 8 µl ethanol, or 8 µl ethanol by itself. The disks were placed at the center of the plates (one disk per plate), and the plates were incubated at 30°C for 3 days (linoleic acid plates) or 4 days [13(S)-HpODE] along with ethanol control plates at each time point.
Following incubation, an 11-mm core containing the disk was removed and discarded. A 23-mm ring was punched out around the removed portion and homogenized in 3 ml Tween H2O. A 100-µl aliquot was removed and diluted, and spores were counted on a hemocytometer.
Expression analysis.Publicly available RNA-seq and microarray expression data for A. flavus under various conditions was accessed from the publication itself or from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/), depending on where the data were deposited. For the GEO data, the GEO2R Web tool was used to compare different samples within an experiment. The log2-fold change (FC) values and P values for A. flavus’ 15 GPCRs were collected from the various data sets.
Statistical analysis.Statistically significant differences were determined by using an unpaired Student’s t test with a two-tailed distribution and a P value of <0.05.
ACKNOWLEDGMENTS
This work was supported by NSF grant IOS-0965649 (to N.P.K. and K.J.A.).
FOOTNOTES
- Received 17 June 2014
- Accepted 15 September 2014
- Published 14 October 2014
- Copyright © 2014 Affeldt et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.