#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Examining the Evolution of the Regulatory Circuit Controlling Secondary Metabolism and Development in the Fungal Genus


Filamentous fungi produce a highly diverse cadre of secondary metabolites, small molecules whose potent toxic activities are integral to the fungal lifestyle. Most secondary metabolites are narrowly taxonomically distributed, whereas the transcriptional regulators that control their production, alongside with controlling other key processes such as development, are broadly conserved. To gain insight into the evolution of the regulatory circuit governing secondary metabolism and development, we examined the evolution of the genes and pathways underlying these processes as well as the evolution of their transcriptional regulation in the filamentous fungal genus Aspergillus, a prolific and important producer of secondary metabolites. We discovered that although secondary metabolic genes and pathways are poorly conserved across Aspergillus, their regulation is largely conserved. In contrast, the regulation of highly conserved developmental genes is divergent. These results point to a new type of rewiring that has occurred during the evolution of the regulatory circuit governing secondary metabolism and development in Aspergillus in which conserved regulators control a conserved biological process (secondary metabolism), even though the underlying genes and pathways that make up the biological process are not themselves conserved.


Published in the journal: . PLoS Genet 11(3): e32767. doi:10.1371/journal.pgen.1005096
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1005096

Summary

Filamentous fungi produce a highly diverse cadre of secondary metabolites, small molecules whose potent toxic activities are integral to the fungal lifestyle. Most secondary metabolites are narrowly taxonomically distributed, whereas the transcriptional regulators that control their production, alongside with controlling other key processes such as development, are broadly conserved. To gain insight into the evolution of the regulatory circuit governing secondary metabolism and development, we examined the evolution of the genes and pathways underlying these processes as well as the evolution of their transcriptional regulation in the filamentous fungal genus Aspergillus, a prolific and important producer of secondary metabolites. We discovered that although secondary metabolic genes and pathways are poorly conserved across Aspergillus, their regulation is largely conserved. In contrast, the regulation of highly conserved developmental genes is divergent. These results point to a new type of rewiring that has occurred during the evolution of the regulatory circuit governing secondary metabolism and development in Aspergillus in which conserved regulators control a conserved biological process (secondary metabolism), even though the underlying genes and pathways that make up the biological process are not themselves conserved.

Introduction

Filamentous fungi produce diverse repertoires of small molecules known as secondary metabolites (SMs) [1]. SMs include widely used pharmaceuticals such the antibiotic penicillin [2], the cholesterol-reducing drug lovastatin [3], and the immunosuppressant cyclosporin [4], as well as potent mycotoxins, such as aflatoxin [5] and fumonisin [6,7]. SMs play key ecological roles in territory establishment and defense, communication, and virulence [812].

The genes involved in fungal SM pathways are often physically linked in the genome, forming contiguous SM gene clusters [13]. These gene clusters are typically characterized by a backbone gene, such as those encoding nonribosomal peptide synthetases (NRPSs), polyketide synthases (PKSs), hybrid NRPS-PKS enzymes, and prenyltransferases, whose protein products are responsible for synthesizing the proto-SM. Additional genetic components of SM gene clusters include genes for one or more tailoring enzymes that chemically modify SM precursors, transporter genes responsible for exporting the final product, and transcription factors that drive expression of the remaining genes in the gene cluster. For example, the gene cluster responsible for the synthesis of the mycotoxin gliotoxin in the opportunistic human pathogen Aspergillus fumigatus contains 13 genes including a non-ribosomal peptide synthase (gliP), multiple tailoring enzymes (gliI, gliJ, gliC, gliM, gliG, gliN, gliF), a transporter gene (gliA), a transcription factor (gliZ), and a gliotoxin oxidase gene that protects the fungus from the harmful effects of gliotoxin (gliT) [14,15].

Filamentous fungi exhibit a huge amount of SM biochemical diversity. Individual SMs are often known to be produced by only one or a handful of species, and the SM chemotypic profiles of closely related fungi are typically non-overlapping [1,16]. For example, the meroterpenoid fumagillin, originally isolated from A. fumigatus, has only been detected in A. fumigatus and some isolates of Penicillium raistrickii [17,18]. The gene cluster required for its production appears to be conserved in the A. fumigatus close relative, Aspergillus fischerianus, though only intermediate compounds have been detected from cultures of this and other closely related species [1921]. In some genera, including Aspergillus [22], the extent of fungal SM distribution is so taxonomically narrow that SM chemotypic profiles have been used as unequivocal species-level identifiers.

As might be expected given their key roles in fungal ecology, SM production–and as a consequence SM gene cluster transcriptional activity–is tightly controlled by a complex network of master SM regulators triggered by a wide variety of environmental cues such as temperature, light, pH, and nutrient availability [23]. Among the master SM regulators identified to date are members of the fungal-specific Velvet protein family, which regulate SM production in a light-dependent manner in the model filamentous fungus Aspergillus nidulans [2427]. The founding member of the Velvet family, VeA, stimulates production of diverse types of SMs in various fungal genomes under dark conditions, and has been shown to regulate gliotoxin, fumagillin, fumitremorgin G, and fumigaclavine C gene cluster expression and metabolite production in A. fumigatus [28]. Recently, a VeA-dependent regulator of secondary metabolism, MtfA, was identified in A. nidulans, which–unlike VeA–is localized in the nucleus regardless of light conditions [29]. MtfA regulates terrequinone, sterigmatocystin, and penicillin in A. nidulans; in A. fumigatus, MtfA is necessary for normal protease activity, and virulence assays using the moth Galleria mellonella suggest it plays a role in pathogenicity [30].

In addition to regulating SM, both of these regulators have been linked to the regulation of asexual and sexual development. Timing of SM production with developmental changes is well established in filamentous fungi, and the presence/absence of certain SMs has been linked with developmental changes [3133]. It has been suggested that regulators that coordinate SM and development allow filamentous fungi to support more “complex” lifestyles through the production of a much greater diversity of natural products than their unicellular yeast relatives, which lack veA as well as backbone synthesis genes necessary for SM production [3335].

Remarkably, both veA and mtfA appear to be broadly conserved in filamentous fungi with non-overlapping SM profiles [29,36]. We used four well-studied organisms from the fungal genus Aspergillus, a highly diverse genus and producer of some of the most iconic SMs, including gliotoxin and penicillin, to investigate the evolutionary variability in the distribution of SM gene clusters and its interaction with these two broadly conserved global transcriptional regulators that differ in their response to light, veA and mtfA. Our evolutionary analyses show that although both the SM gene clusters as well as their gene content are poorly conserved between Aspergillus species, explaining the narrow taxonomic distribution and distinctiveness of their SM profiles, the effects of the global transcriptional regulators on SM production in response to environmental cues are largely conserved across these same species. In contrast, examination of the role of veA and mtfA in development, a process that involves genes that are highly conserved between the two species and whose regulation is intimately linked to SM regulation, yields a different pattern; whereas the role of veA is conserved, mtfA regulates development in the homothallic A. nidulans but not in the heterothallic A. fumigatus.

Although rewiring has been well documented in yeast and animal regulatory networks [3744] these studies largely concern rewirings of regulatory signals between otherwise conserved regulators and target genes. Our finding that VeA, and to a certain extent MtfA, regulate non-homologous SM pathways in both A. fumigatus and A. nidulans suggests a new type of rewiring in which conserved regulators control a conserved biological process, even though the underlying genes and pathways that make up the biological process are not themselves conserved.

Results

The majority of SM gene clusters in Aspergillus are species-specific

The genomes of A. fumigatus, A. nidulans, Aspergillus oryzae, and Aspergillus niger contain 317, 498, 725, and 584 secondary metabolic genes, respectively, which are organized in 37, 70, 75, and 78 corresponding secondary metabolic gene clusters (S1 Table) [45]. We considered SM gene clusters to be conserved between species if greater than half of the genes in the larger gene cluster were orthologous to greater than half of the genes in the smaller gene cluster. Even with this very liberal definition of gene cluster conservation, we found that no SM gene clusters were conserved across all four species. Moreover, 91.9–96.1% of SM gene clusters were specific to each species, with only 7 SM gene clusters conserved between any species (Fig. 1A, S2 Table). Only one two-gene cluster is completely conserved between two species (A. fumigatus 27 and A. oryzae 9). Although none of these conserved SM gene clusters have chemically characterized products, the fact that only one of the 7 gene clusters appear to be 100% conserved in only two of the four Aspergillus species, suggests that the SM chemotypic profiles of the four species are non-overlapping.

Fig. 1. SM gene clusters in Aspergillus are not evolutionarily conserved.
SM gene clusters in <i>Aspergillus</i> are not evolutionarily conserved.
A. Venn diagram showing homologous SM gene clusters between A. fumigatus, A. nidulans, A. niger, and A. oryzae. Two SM gene clusters were considered homologous if greater than 50% of their genes were orthologs. Numbers in parenthesis indicate the total number of SM gene clusters present in each species. The asterisk (*) denotes that two SM gene clusters in A. oryzae are homologous to one gene cluster in A. niger (S2 Table). B. Circos plot showing the evolutionary conservation of the A. nidulans sterigmatocystin gene cluster, one of the most highly conserved SM gene clusters in our study, across Aspergillus. in A. fumigatus, A. nidulans, A. niger, and A. oryzae. The outer black track shows the relative SM gene counts in A. fumigatus, A. nidulans, A. niger, and A. oryzae. SM gene clusters in each species are indicated by the alternating light and dark grey wedges of the inner track; wedge thickness is proportional to number of clustered genes. The sterigmatocystin gene cluster in A. nidulans is colored red. Links indicate SM clusters containing one or more genes assigned to the same orthogroup as gene(s) in the sterigmatocystin gene cluster (S3 Table); link color indicates the number of shared genes.

While only one (at the 100% level) or very few (at the 50% level) conserved SM gene clusters can be identified in comparisons between any of these four species, SM gene clusters do contain genes whose orthologs are parts of other, non-homologous, SM gene clusters. For example, the 25 genes in the sterigmatocystin gene cluster in A. nidulans, one of the largest SM gene clusters present in the genomes analyzed, have orthologs in 25 SM gene clusters in the other three species as well as inparalogs in 8 other A. nidulans SM gene clusters (Fig. 1B, S3 Table). However, in all but one case, less than 20.0% (5 genes) of the sterigmatocystin gene cluster is present in the other gene cluster. The only exception is the truncated aflatoxin gene cluster of A. oryzae, which shares 11 orthologs with the ST gene cluster. Although the A. oryzae aflatoxin gene cluster is non-functional [4648], the evolutionary conservation between the aflatoxin and sterigmatocystin gene clusters is reflected in the fact that sterigmatocystin is the penultimate precursor product of the aflatoxin biosynthetic pathway [49].

Aspergillus SM genes are significantly less conserved than genes for primary metabolism

Given the remarkable lack of conservation of SM gene clusters between Aspergillus species, we next tested whether the genes belonging to these clusters were also less conserved by comparing their degree of evolutionary conservation to that of genes involved in primary metabolism. To determine the percentage of species-specific orthogroups (see Methods) involved in SM and primary metabolism, we first determined the number of orthogroups with genes annotated to these functional categories that contained no genes in any other species (i.e., all orthogroups that contain only a single gene or only putative inparalogs). To determine the conservation of gene function across species, we determined the number of orthogroups in each genome that are either entirely species-specific or whose members in other genomes are annotated to different functional categories.

We found that SM orthogroups were significantly far less conserved than primary metabolic orthogroups in all four genomes examined (adjusted P < 1e-10 for all combinations; S4 Table). Specifically, the percentage of species-specific primary metabolic orthogroups ranged between 7.5 and 15.4% (Fig. 2). When we considered conservation of function (i.e. whether orthogroups contained genes annotated to primary metabolism), we saw a slight and non-significant increase in the percentage of species-specific orthogroups (8.9–18%; Fig. 2). In contrast, 25.0% to 40.5% of orthogroups containing SM backbone genes were species-specific; considering functional conservation had a negligible impact on these percentages (Fig. 2). Orthogroups containing genes in SM clusters exhibited similar percentage ranges for species-specificity (26.0–38.1%; Fig. 2). Strikingly, the orthologs of many SM genes in a given species were not in SM gene clusters in any of the other species; specifically, 53.7% (in A. fumigatus, with 37 SM clusters) to 74.7% (in A. niger, with 77 SM gene clusters) of orthogroups containing an SM gene from a particular genome had no orthologs in SM clusters in any of the other genomes.

Fig. 2. Many more Aspergillus SM genes are species-specific than primary metabolic genes.
Many more <i>Aspergillus</i> SM genes are species-specific than primary metabolic genes.
For each species, dark blue bars indicate the percentage of primary metabolism orthogroups that is species-specific when compared to all genes annotated as participating in primary metabolism in the other three species; light blue bars indicate the percentage of each species’ primary metabolism orthogroups that is species-specific compared to all genes, irrespective of their annotation, in the other three species. Similarly, light orange bars indicate the percentage of each species’ SM backbone synthesis orthogroups that is species-specific when compared to all genes annotated as SM backbone synthesis genes in the other three species; dark orange bars indicate the percentage of each species’ SM backbone synthesis orthogroups that is species-specific compared to all other genes. Light red bars indicate the percentage of each species’ clustered SM orthogroups that is species-specific when compared to all genes annotated as clustered SM genes in the other three species; dark red bars indicate the percentage of each species’ clustered SM orthogroups that is species-specific compared to all other genes. Asterisks indicate statistically significant differences based on a P-value ≤ 0.01 (*), ≤ 0.001 (**), or ≤ 0.0001 (***) in a two-tailed Fisher’s exact test (S4 Table).

VeA regulates the same biological processes as well as the same fraction of the genome in both A. nidulans and A. fumigatus

We next examined the function of the conserved secondary metabolic regulator VeA by performing RNA sequencing [19,48,50] of ΔveA and wild-type (WT) A. fumigatus strains TSD1.15 [51] and CEA10 and A. nidulans strains TXFp2.1 and TRV50.2 [52] to identify genes and biological processes that are differentially regulated in ΔveA vs WT in the two species. Of the 9,783 transcribed genes in the A. fumigatus genome, 1,546 (15.8%) were over-expressed and 1,555 (15.9%) were under-expressed in the ΔveA vs WT analysis in A. fumigatus (Tables 1, S5). We observed very similar numbers of genes differentially regulated in the A. nidulans ΔveA vs WT analysis; out of 10,709 genes in the A. nidulans genome, 1,165 (10.9%) were over-expressed and 1,671 genes (15.6%) were under-expressed. In total, approximately 32% and 26% of protein coding genes were differentially regulated in ΔveA compared to WT in A. fumigatus and A. nidulans, respectively.

Tab. 1. Differentially expressed genes in ΔveA vs WT and ΔmtfA vs WT comparisons for A. fumigatus and A. nidulans.
Differentially expressed genes in Δ<i>veA</i> vs WT and Δ<i>mtfA</i> vs WT comparisons for <i>A</i>. <i>fumigatus</i> and <i>A</i>. <i>nidulans</i>.
aNumber of differentially expressed genes relative to wild type

To characterize the broad biological processes that these differentially regulated genes are involved with, we performed GO term enrichment analysis using the Aspergillus GOSlim term hierarchy [53,54]. Remarkably, the same four GO terms, namely secondary metabolic process, carbohydrate metabolic process, oxidoreductase activity, and extracellular region, were significantly enriched in under-expressed genes in both A. nidulans and A. fumigatus, showing that VeA is a positive regulator of the same processes in both species (Fig. 3). Over-expressed genes in A. fumigatus were significantly enriched for twelve GO terms potentially related to cell growth, namely ribosome biogenesis, cellular amino acid metabolic process, translation, rna metabolic process, structural molecule activity, helicase activity, rna binding, transferase activity, nucleus, nucleolus, ribosome, and cytosol. Five of these twelve terms were also significantly enriched in A. nidulans (ribosome biogenesis, cellular amino acid metabolic process, rna metabolic process, nucleus, and nucleolus). Over-expressed genes were present in the remaining seven terms in A. nidulans but did not show statistically significant enrichment (S6 Table).

Fig. 3. GO term enrichment analysis of genes differentially expressed in ΔveA and ΔmtfA relative to wild type in A. fumigatus and A. nidulans.
GO term enrichment analysis of genes differentially expressed in Δ<i>veA</i> and Δ<i>mtfA</i> relative to wild type in <i>A</i>. <i>fumigatus</i> and <i>A</i>. <i>nidulans</i>.
Statistically overrepresented gene ontology (GO) categories in under-expressed (red) and over-expressed (blue) gene sets in ΔveA and ΔmtfA relative to wild type. Arrows point to GO term ancestors. Horizontal bars show the percentage of each gene set assigned to a particular GO term with black bars indicating significant enrichment (Benjamini & Hochberg adjusted P-value ≤ 0.05 in a hypergeometric test; S6 Table); white bars indicate no significant enrichment.

In addition to the enrichment of the GO term secondary metabolic process in the differentially expressed genes from the ΔveA vs WT comparison in both A. nidulans and A. fumigatus, a large portion of the 317 A. fumigatus SM cluster genes and 498 A. nidulans SM cluster genes was also differentially expressed in the ΔveA strains. In A. fumigatus, 98 genes (30.9%) were under-expressed and 38 (12.0%) were over-expressed; in A. nidulans, 184 genes (36.9%) were under-expressed and 67 (13.5%) were over-expressed (S5 Table). Interestingly, all constituent genes of several SM gene clusters were differentially expressed. For example, all genes in the A. fumigatus pseurotin A gene cluster and all genes in the A. nidulans asperthicin cluster were under-expressed (S10 Table).

MtfA’s regulatory role in secondary metabolism is smaller in scope in A. fumigatus compared to A. nidulans

We next examined the role of the recently identified SM regulator MtfA [29,30] in A. fumigatus and A. nidulans by performing RNA sequencing and differential gene expression analysis of ΔmtfA vs WT strains of both species (A. fumigatus tTDS4.1 ΔmtfA [30] and CEA10, A. nidulans TRVp ΔmtfA and TRV50.2 [29]). In contrast to our findings with veA, we found an approximately 10-fold difference in the percentage of genes regulated in both species (Tables 1, S5). Thirty-six genes were over-expressed (0.4%) and 63 (0.6%) were under-expressed in the A. fumigatus ΔmtfA vs WT analysis, whereas in the A. nidulans ΔmtfA vs WT analysis 400 genes were over-expressed (3.7%) and 568 were under-expressed (5.3%).

To determine the functional categories impacted by mtfA deletion both species, we performed GO term enrichment analysis on the genes differentially expressed between ΔmtfA and WT strains. Under-expressed as well as over-expressed genes in A. nidulans were significantly enriched for secondary metabolic process, toxin metabolic process and oxidoreductase activity, suggesting that MtfA is involved in positive and negative regulation of different secondary metabolites (Fig. 3). Over-expressed genes in A. nidulans were also significantly enriched for asexual developmental processes, namely developmental process and asexual sporulation.

Under-expressed genes in A. fumigatus were significantly enriched for two of the three processes as in A. nidulans, namely secondary metabolic process and oxidoreductase activity. However, over-expressed genes in A. fumigatus were not significantly enriched for any GO terms; some over-expressed genes were present in the secondary metabolic process term, though this was not statistically significant (S6 Table).

Examination of the 317 A. fumigatus SM cluster genes and 498 A. nidulans SM cluster genes showed that they too were also differentially expressed in both A. fumigatus and A. nidulans ΔmtfA mutants compared to the WT strains. In A. nidulans, 107 of the 498 SM cluster genes (21.5%) were under-expressed and 32 (6.4%) were over-expressed (S5 Table). In some SM gene clusters, such as in the dba cluster, all genes were under-expressed in ΔmtfA A. nidulans relative to WT (S10 Table). In contrast, many fewer SM cluster genes were differentially expressed in A. fumigatus ΔmtfA; 32 out of the 317 genes (10.1%) were under-expressed and 4 (1.3%) were over-expressed (S5 Table). Despite this much smaller number of differentially expressed SM genes, at least one entire SM gene cluster (the pseurotin A cluster) was under-expressed in A. fumigatus ΔmtfA relative to ST (S10 Table).

Transcriptional regulation of similar processes regardless of gene conservation

To examine whether gene conservation and involvement in a biological process (secondary metabolism and development) correlated with conservation of regulation by VeA and MtfA, we tested whether orthologous and species-specific genes in A. nidulans and A. fumigatus showed the same or different responses in ΔveA vs WT and ΔmtfA vs WT analyses (Fig. 4). Even though both species contain large numbers of species-specific genes in their SM gene clusters, the proportions of differentially expressed SM cluster genes in both ΔveA A. nidulans and ΔveA A. fumigatus were remarkably similar (Figs. 1A, 4A; S5 Table).

Fig. 4. Orthology of SM and development genes differentially expressed in ΔveA and ΔmtfA relative to wild type in A. fumigatus and A. nidulans.
Orthology of SM and development genes differentially expressed in Δ<i>veA</i> and Δ<i>mtfA</i> relative to wild type in <i>A</i>. <i>fumigatus</i> and <i>A</i>. <i>nidulans</i>.
Circos plots of genes in SM clusters (A,B) and developmental genes (C,D) showing change in gene expression patterns in ΔveA vs WT (A,C) and in ΔmtfA vs WT (B,D) experiments. In all plots, A. fumigatus genes are shown on the left and A. nidulans genes on the right. Inner tracks show the relative number of under-expressed genes (red), over-expressed genes (blue) and not differentially expressed genes (grey) in each species. Links between the two species indicate orthologous genes that are both under-expressed (red links), both over-expressed (blue links) and both not differentially expressed (light grey links); purple links indicate that the orthologous genes have conflicting expression patterns. In both A and B, the 317 SM cluster genes present in A. fumigatus are shown on the left and the 498 A. nidulans SM cluster genes are shown on the right. Many SM cluster genes are differentially expressed in the ΔveA vs. WT comparison in both species (A). Many SM cluster genes are also differentially expressed in the ΔmtfA vs. WT comparison in A. nidulans; the proportion of differentially expressed genes in the ΔmtfA vs. WT comparison in A. fumigatus is much smaller (B). The absence of links between most of the SM cluster genes in the two species indicates that most SM cluster genes are species-specific. For those genes that share orthology in the two species, they often differ in their expression in response to a master transcriptional regulator. In C and D, the 478 genes annotated to the GO term developmental process in A. fumigatus are shown on the left side of the plot and the 490 genes annotated to the same GO term in A. nidulans are shown on the right side. The presence of links between most genes in these two species indicates a high degree of gene conservation. The expression patterns of all individual genes, including those present in SM gene clusters and those involved in developmental processes, are shown in S5 Table. Expression patterns of the SM cluster genes organized by individual SM pathways are shown in S10 Table.

Examination of SM cluster genes that are differentially expressed suggests that even when genes in SM clusters have orthologs, these orthologs were often differentially regulated by ΔveA. Entire SM pathways that are over- or under-expressed in one species may have orthologs with completely different expression patterns in the other species, as in the case of the dba SM gene cluster in A. nidulans (S2 Fig). Specifically, of the 184 under-expressed SM cluster genes in ΔveA A. nidulans, 64 genes (34.8%) had an ortholog in A. fumigatus (S7 Table). Of these 64 genes, 45 (70.3%) had at least one differentially expressed ortholog in A. fumigatus, and 37 (57.8%) had at least one similarly under-expressed ortholog in A. fumigatus (S7 Table). Fewer SM genes were over-expressed in either ΔveA A. nidulans or A. fumigatus; of the 67 over-expressed genes in A. nidulans, only 14 (20.9%) had orthologs in A. fumigatus. Of these 14 genes, 5 had at least one differentially expressed ortholog in A. fumigatus, and 3 had at least one similarly over-expressed ortholog in A. fumigatus (S7 Table).

When mtfA was deleted, fewer SM cluster genes were differentially expressed in A. fumigatus than in A. nidulans. Of the 107 under-expressed genes in ΔmtfA A. nidulans, 36 (33.6%) had an ortholog in A. fumigatus (Fig. 4B, S7 Table). Unlike veA, however, only 6 of these conserved genes had differentially expressed orthologs in A. fumigatus. Finally, of the 32 over-expressed genes in ΔmtfA A. nidulans, 2 of the 12 genes with orthologs in A. fumigatus had orthologs that were differentially expressed.

Apart from their involvement in the global regulation of SM, both VeA and MtfA are also involved in the regulation of asexual and sexual development. In contrast to genes involved in SM, genes involved in asexual and sexual development in Aspergillus have been shown to be highly conserved across the genus [55]. Of the 490 genes annotated to the GO term developmental process in A. nidulans, 462 have at least one ortholog among the 478 genes annotated to this term in A. fumigatus. In ΔveA A. nidulans, 72 developmental genes are under-expressed and 32 are over-expressed. Of the 72 under-expressed genes, 66 (91.7%) have an ortholog in A. fumigatus; of these 66 orthologs, 30 were differentially expressed in both species (Fig. 4C, S7 Table). There were fewer over-expressed developmental genes in ΔveA A. nidulans, but they showed similar trends; 31 of the 32 over-expressed genes have an ortholog, 15 of which had differentially expressed orthologs in A. fumigatus and 11 of which had over-expressed orthologs. In contrast with veA, many more developmental genes were differentially expressed in A. nidulans ΔmtfA (35) than in A. fumigatus ΔmtfA (1). While 4 of the 6 under-expressed genes and 28 of the 29 over-expressed genes in A. nidulans had orthologs in A. fumigatus, none of these orthologs were differentially expressed (Fig. 4D, S7 Table).

Discussion

To gain insight into the evolution of the regulatory circuit governing SM and development in filamentous fungi, we examined the evolution of their underlying genes and pathways as well as the evolution of their regulation by two global transcriptional regulators in Aspergillus, a genus whose members prolifically produce diverse secondary metabolites. We discovered that although secondary metabolic genes and pathways are poorly conserved, their regulation by VeA, and to a certain extent by MtfA, is largely conserved. In contrast, the MtfA regulation of highly conserved developmental genes is divergent, whereas that of VeA is conserved. Below, we discuss the significance of our results with respect to the evolution of SM gene clusters as well as to the evolution of the regulatory circuit controlling secondary metabolism and development in filamentous fungi.

Aspergillus secondary metabolic genes and gene clusters are largely species-specific

All but one of the SM gene clusters present in the four Aspergillus species we examined were species-specific. Even when we used a very low threshold of 50% evolutionary conservation, we found that there were no clusters conserved in all four species, one cluster conserved in three species, and a very small number conserved between pairs of species (Fig. 1A). Consistent with our results, an examination of the close relatives A. fumigatus, A. fischerianus, and Aspergillus clavatus, using an 80% threshold of evolutionary conservation, also found relatively small numbers of conserved SM gene clusters [56]. Remarkably, variation in SM gene cluster content can sometimes be even strain-specific; for example, a single SM gene cluster has been shown to vary in its presence in A. fumigatus isolates [57], and several SM gene clusters also appear to vary between isolates of A. niger [58].

Surveys of fungal genomes show that for any given fungus there are many more gene clusters than known SMs, suggesting that the currently characterized SMs might be only a small fraction of the SMs that a fungus can produce [56,59]. For example, 33 of the 37 putative SM gene clusters in A. fumigatus have no characterized products, despite evidence from metabolomics surveys suggesting that the fungus produces many SMs [19,45,60]. Interestingly, most known SMs are produced by only a handful of species [16] and SM profiles between closely related species are quite distinct [17,61], both indications that the typical taxonomic distribution of SMs is narrow. Thus, the observed extremely high degree of taxonomic narrowness of SM gene clusters is not only consistent with the distribution of known SMs, but also suggests that this distribution is likely to be typical of all the SMs produced by a fungus.

SM gene clusters are largely species-specific, so one hypothesis is that their constituent genes are also species-specific. Another alternative is that their constituent genes are conserved but have been reshuffled and become members of different SM gene clusters in different species. Our analysis shows support for both hypotheses. For example, we found that a much larger fraction of the genes comprising SM gene clusters is species-specific than of genes participating in primary metabolism (Fig. 2), which is likely explained by extensive gene duplication and loss [35,62,63] de novo gene emergence [64], horizontal gene transfer [6569], as well as by very high sequence divergence. However, a considerable fraction of SM genes was not species-specific. In several cases, the orthologs of SM genes in one Aspergillus species were found residing in non-homologous SM pathways in another species; for example, the 25 genes in the sterigmatocystin gene cluster in A. nidulans have orthologs in 25 distinct SM gene clusters in the other three species (Fig. 1B). Surprisingly, the orthologs of many genes that resided within an SM gene cluster in one species were not in an SM gene cluster in another species (Fig. 2). The majority of these orthologs either lack annotation or are annotated as being involved in primary metabolism, consistent with a model in which the gene content of SM gene clusters is formed or altered through the recruitment and incorporation of native genes involved in primary metabolism and other essential cellular processes.

The evolution of the circuit regulating secondary metabolism and development in Aspergillus

By comparing genome-wide gene expression of deletion mutants of veA and mtfA with wild-type strains in both A. fumigatus and A. nidulans we assessed the degree to which their regulatory roles in controlling secondary metabolism and development are conserved (Fig. 4A and 4C). The involvement of VeA in regulating both processes is highly conserved in both A. fumigatus and A. nidulans. Remarkably, however, the downstream SM genes regulated by VeA are different between the two species; this difference is not fully accounted by the fact that VeA regulates many species-specific SM genes. Many genes that are regulated by VeA in one species are present but are not regulated by VeA in the other (Fig. 4A).

This conservation in regulatory logic (i.e., conservation in the regulation of secondary metabolism) despite the dramatic change in the genes involved may be explained by considering the number and complexity of VeA’s interacting partners as well as VeA’s putative transcription factor function, which in combination offer many degrees of freedom for changes in the regulation of specific SM genes and pathways in either species. Specifically, VeA is known to have many interacting partners [33] including the Velvet family protein VelB, which it transports from the cytoplasm to the nucleus, where both proteins form a trimeric complex with LaeA that regulates secondary metabolism production and development [24]. Furthermore, VeA also functions outside of this complex [24,28]; it interacts with red light-sensing proteins in the nucleus [27] and with other methyltransferases [70], and it has been hypothesized that it may also act as a scaffold protein for the recruitment of additional transcriptional regulators [33,36]. Finally, recent analysis has shown that the Velvet domain is a DNA-binding domain, and that Velvet family proteins may act as direct transcriptional regulators [71].

The involvement of MtfA in regulating SM in both A. fumigatus and A. nidulans is also conserved, although the numbers of SM genes that are under its control differ considerably between the two species. As with veA, comparing expression patterns between ΔmtfA and WT identified differentially expressed species-specific and conserved SM genes. Again like veA, many conserved genes that are regulated by MtfA in one species are not regulated in the other (Fig. 4B). This rewiring could either be due to a divergence in the signal that targets these genes for regulation by MtfA, a C2H2 zinc finger transcription factor [29], or its interacting partners. Although not much is known about MtfA’s interacting partners, indirect support for their presence comes from our finding that genes involved in SM are significantly over-represented in both the under-expressed and the over-expressed gene sets in A. nidulans (Fig. 3), suggesting that the regulatory effect of MtfA in certain regulatory partnerships is positive and in others negative. Interestingly, our results also suggest that MtfA genetically interacts with (and acts downstream of) VeA in A. nidulans, but not in A. fumigatus [2830], as its expression is decreased in ΔveA vs WT in A. nidulans but not in ΔveA vs WT in A. fumigatus (S5 Table). Finally, MtfA appears to regulate development in the homothallic A. nidulans, but not in the heterothallic A. fumigatus (Figs. 3,4). Unlike veA, deleting mtfA influenced the expression of developmental genes only in A. nidulans but not in A. fumigatus. This is consistent with the loss in A. fumigatus or gain in A. nidulans of the signal that directs MtfA or its downstream targets to regulate developmental processes.

VeA’s central role in coordinating SM and development under dark conditions as well as its large number of interacting partners likely explain in part why many more genes are differentially expressed in its absence in both A. nidulans and A. fumigatus than in the absence of MtfA. Gaining a more complete understanding of the molecular mechanisms through which VeA and MtfA regulate downstream targets is an interesting line of future inquiry that harbors significant promise in elucidating how VeA and MtfA have been rewired in different fungal species. However, it is likely that VeA and MtfA globally regulate SM and developmental processes not only directly but also indirectly, through interactions with other regulatory proteins; thus, obtaining a full mechanistic understanding of the rewiring of this regulatory circuit will likely also require characterization of VeA and MtfA’s interacting partners and their regulatory functions.

A novel type of regulatory circuit rewiring?

It is abundantly clear that changes in transcriptional regulation, also known as rewiring, are a major driver of phenotypic divergence [40,43]. It is also becoming increasingly clear that rewiring of regulatory circuits can also take place in the absence of phenotypic change [42,72,73]. Irrespective of whether it leads to phenotypic change or not, rewiring is typically thought to occur either through changes to the transcriptional regulators or through changes to the regulatory signals between the regulators and the target genes.

An example of rewiring due to changes to the transcriptional regulator is offered by the galactose pathway, which is controlled by Gal4p in the baker’s yeast Saccharomyces cerevisiae, but by the non-related Cph1p in the human commensal Candida albicans [37,74]. The acquisition or loss of specific motifs in otherwise conserved transcriptional regulators can also lead to rewiring [75,76], as can acquisition or loss of interacting regulatory partners [41]. An example of changes in regulatory signals is offered by the small percentage of Mcm1 regulated genes shared between S. cerevisiae and C. albicans [77], which can be largely attributed to high rates of cis-regulatory sequence gain and loss that Mcm1 binds to. This is likely to be the major explanation in cases that involve conserved regulators, conserved target genes, and conserved phenotypes [78,79].

Although the evidence for the regulatory impact of VeA and MtfA on secondary metabolism and development is genetic in our case (i.e., it is not known whether VeA or MtfA directly bind to target gene cis-regulatory elements or whether they control their expression indirectly), the evolution of VeA and MtfA regulation on development is consistent with rewiring that involves changes either to the transcriptional regulator or the regulatory signals between the regulator and the target genes (Fig. 5A,B). In the case of VeA, the rewiring is not associated with phenotypic change–genes involved in development are regulated by VeA in both A. fumigatus and A. nidulans, whereas in the case of MtfA the rewiring might be associated with the differences in sexual and asexual development between the two species [55,80].

Fig. 5. Model of the evolution of the regulatory circuit controlling secondary metabolism and development in Aspergillus.
Model of the evolution of the regulatory circuit controlling secondary metabolism and development in <i>Aspergillus</i>.
Generalized gene regulatory networks for development (A,B) and secondary metabolism (C,D) for the transcriptional regulators VeA (A,C) and MtfA (B,D) in A. fumigatus and A. nidulans. As master transcriptional regulators, VeA and MtfA can promote the expression of a gene (indicated by an arrow →) or suppress it (indicated by a bar ⟞). Target genes are either conserved (i.e., present in both species; green squares) or not (i.e., species-specific; blue squares). Similarly, the signal emanating from each transcriptional regulator that coordinates target gene regulation is either conserved (i.e., acting in the same fashion on genes conserved in both species; green arrows/bars) or not (i.e., acting in different fashion on genes conserved in both species or acting on species-specific genes; blue arrows/bars). Both transcriptional regulators contribute to the regulation of secondary metabolism in both A. fumigatus and A. nidulans even though the underlying genes in the two species are largely non-homologous (C,D). In contrast, the role of these two regulators in development, where most underlying genes are conserved in both species, is different; the role of veA is conserved, but mtfA regulates development in the homothallic A. nidulans but not in the heterothallic A. fumigatus (A,B).

In contrast to the rewiring examples mentioned so far, the major driver of change in the evolution of VeA and MtfA regulation on SM appears to be the dramatic turnover in the target genes involved in the process, which has resulted in their high degree of species specificity and their even higher species-specific presence in SM gene clusters (Figs. 2 and 5C,D). Although it seems reasonable to exclude the transcriptional regulators as potential co-drivers of the rewiring because the key domains of both VeA and MtfA are conserved across Aspergillus [29,36], it is possible that the regulatory signals or the differential presence or absence of interacting regulatory partners may have also contributed to rewiring. However, given the vital ecological importance of secondary metabolites in fungal ecology, it seems reasonable to hypothesize that the evolution of the regulatory circuit governing secondary metabolism in filamentous fungi is primarily driven by the likely extreme evolutionary pressure imposed on each fungus to produce its own unique blend of secondary metabolites.

Materials and Methods

Genome sequences and orthogroup definitions for A. nidulans, A. fumigatus, A. oryzae, and A. niger

All genome sequences and annotations for A. nidulans FGSC A4 s10-m02-r03, A. fumigatus AF293 s03-m04-r11, A. oryzae RIB40 s01-m08-r21 and A. niger CBS 513.88 s01-m06-r10 were taken from the Aspergillus Genomes Database (AspGD) [54]. Groups of orthologous genes (orthogroups) for these four genomes were taken from AspGD’s orthology assignments for 16 Aspergillus species, which were generated using a Jaccard clustering approach [81]. AspGD orthogroups contain groups of genes that are thought to have descended from the Aspergillus common ancestor; genes from the same species that are part of a given orthogroup are defined as in-paralogs that have duplicated at some later point after the species diverged from the Aspergillus common ancestor. Species-specific genes, which were absent from AspGD orthogroups, were organized into species-specific orthogroups using the MCL algorithm in combination with all-versus-all protein BLAST search [82]. Proteins with BLAST hits with 60% query and subject coverage, an e-value of less than 1e-5, and a percent identity of greater than 60% were subsequently clustered in MCL with an inflation parameter of 2 and were considered species-specific orthogroups. Proteins that did not pass the BLAST cutoffs were considered single-gene, species-specific orthogroups.

Gene category definitions

Genes involved in secondary metabolism were taken from a previous study that expertly annotated secondary metabolic gene clusters in the four species under study [45]. Manually curated gene cluster boundaries were used when available. Primary metabolism genes were annotated using a previously described enzyme classification pipeline which utilizes KEGG Enzyme Commission annotations [62]. Genes involved in development were determined from all genes in A. fumigatus and A. nidulans annotated to the GO term developmental process (GO:0032502) in AmiGO [83]. This data was accessed on 2014–07–19.

Strains and culture conditions

The strains used in this study include A. fumigatus CEA10, TSD1.15(ΔveA) and TTDS4.1(ΔmtfA) [30,51] and A. nidulans TRV50.2 [52], TXFp2.1(ΔveA) generated in this study, and TRVpΔmtfA [29]. Deletion and wild-type strains presented isogenic genetic backgrounds, and all strains used in this study are prototrophs. Many A. nidulans studies have used a veA partial deletion [84]. For the present study we generated a strain with a complete deletion of the veA coding region, TXFp2.1(ΔveA). This strain was constructed as follows. First, The veA deletion cassette was obtained by fusion PCR as previously described [85]. A 1.4 kb 5’ UTR and a 1 kb 3’ UTR veA flanking regions were PCR amplified from wild type FGSC4 genomic DNA with primers veA_comF and AnidveA_p2, and ANVeASTagP3 and ANVeASTagP4 primers sets, respectively (S8 Table). The A. fumigatus pyrG (pyrGA.fum) selectable marker was amplified with AnidveA_p5 and ANVeASTagP6 primers from plasmid p1439. The 5’ and 3’ UTR fragments were then PCR fused to pyrGA.fum to generate the veA replacement construct using primers AnidveA_P7 and AnidveA_P8. The deletion cassette was transformed into A. nidulans RJMP1.49 strain [86]. The resulting colonies were then transformed with the pSM3 plasmid containing the A. nidulans pyroA to generate a prototroph with a ΔveA background. This strain was confirmed by DNA analysis and designated as TXFp2.1.

All strains were grown in liquid stationary cultures in Czapek-Dox medium (Difco) in the dark. The experiments were carried out with two replicates. After 72 hours of incubation at 37°C mycelial samples were harvested, immediately frozen in liquid nitrogen and lyophilized.

RNA extraction

Total RNA was isolated from lyophilized mycelia using the directzol RNA MiniPrep Kit (Zymo) according to the manufacturer’s instructions. RNA then was quantified using a nanodrop instrument. Expression patterns of veA and mtfA were verified in the A. fumigatus and A. nidulans wild types as well as in the deletion mutants by qRT-PCR prior to RNA sequencing (S1 Fig), confirming the absence of transcripts in the deletion mutants.

RNA sequencing

RNA-Seq libraries were constructed and sequenced at the Vanderbilt Technologies for Advanced Genomics Core Facility at Vanderbilt University using the Illumina Tru-seq RNA sample prep kit as previously described [28,48,50]. In brief, total RNA quality was assessed via Bioanalyzer (Agilent). Upon passing quality control, poly-A RNA was purified from total RNA and the second strand cDNA was synthesized from mRNA. cDNA ends were then blunt repaired and given an adenylated 3’ end. Next, barcoded adapters were ligated to the adenylated ends and the libraries were PCR enriched, quantified, pooled and sequenced an on Illumina HiSeq 2500 sequencer. Two biological replicates were generated for each strain sequenced.

RNA-seq read alignment and differential gene expression

Raw RNA-seq reads were trimmed of low-quality reads and adapter sequences using Trimmomatic using the suggested parameters for single-end read trimming [87]. After read trimming, all samples contained between 20–30 million reads. The smallest sample contained 19.9 million and the largest contained 30.6 million reads; the average sample contained 24.0 million reads (S9 Table). Trimmed reads were aligned to A. nidulans and A. fumigatus genomes using Tophat2 using the reference gene annotation to guide alignment and without attempting to detect novel transcripts (parameter—no-novel-juncs) [88]. Reads aligning to each gene were counted using HTSeq-count with the intersection-strict mode [89]. Differential expression between ΔveA vs WT and ΔmtfA vs WT strains of A. fumigatus and A. nidulans were determined using the DESeq2 software, which normalizes read counts by library depth [90]. Genes were considered differentially expressed if their adjusted P-value was less than 0.1 and their log2 fold change was greater than 1 or less than-1.

Statistical analyses

GO term enrichment was determined for over- and under-expressed genes in all four conditions tested (A. nidulans and A. fumigatus ΔveA vs WT and ΔmtfA vs WT) using the Cytoscape plugin Bingo [91,92]. To allow for a high-level view of the types of differentially expressed gene sets, the Aspergillus GOSlim term subset developed by AspGD was used. The Benjamani-Hochberg multiple testing correction was applied, and terms were considered significantly enriched if the adjusted P-value was less than 0.05.

Fisher’s exact tests were performed using the R function fisher.test with a two-sided alternative hypothesis [93]. P-values were adjusted for multiple comparisons using the R function p.adjust with the Benjamini-Hochberg multiple testing correction [94]. Figures were created using the R plotting system ggplot2 [95] and circos [96].

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11

Attachment 12


Zdroje

1. Keller NP, Turner G, Bennett JW (2005) Fungal secondary metabolism—from biochemistry to genomics. Nat Rev Microbiol 3: 937–947. doi: 10.1038/nrmicro1286 16322742

2. Brakhage AA (1998) Molecular Regulation of β-Lactam Biosynthesis in Filamentous Fungi. Microbiol Mol Biol Rev 62: 547–585. 9729600

3. Kennedy J, Auclair K, Kendrew SG, Park C, Vederas JC, et al. (1999) Modulation of polyketide synthase activity by accessory proteins during lovastatin biosynthesis. Science 284: 1368–1372. 10334994

4. Weber G, Schörgendorfer K, Schneider-Scherzer E, Leitner E (1994) The peptide synthetase catalyzing cyclosporine production in Tolypocladium niveum is encoded by a giant 45.8-kilobase open reading frame. Curr Genet 26: 120–125. doi: 10.1007/BF00313798 8001164

5. Chang PK, Ehrlich KC, Yu J, Bhatnagar D, Cleveland TE (1995) Increased expression of Aspergillus parasiticus aflR, encoding a sequence-specific DNA-binding protein, relieves nitrate inhibition of aflatoxin biosynthesis. Appl Environ Microbiol 61: 2372–2377. 7793958

6. Brown DW, Butchko RAE, Busman M, Proctor RH (2007) The Fusarium verticillioides FUM gene cluster encodes a Zn(II)2Cys6 protein that affects FUM gene expression and fumonisin production. Eukaryot Cell 6: 1210–1218. doi: 10.1128/EC.00400–06 17483290

7. Proctor RH, Brown DW, Plattner RD, Desjardins AE (2003) Co-expression of 15 contiguous genes delineates a fumonisin biosynthetic gene cluster in Gibberella moniliformis. Fungal Genet Biol 38: 237–249. 12620260

8. Yim G, Wang HH, Davies J (2007) Antibiotics as signalling molecules. Philos Trans R Soc Lond B Biol Sci 362: 1195–1200. doi: 10.1098/rstb.2007.2044 17360275

9. Rohlfs M, Albert M, Keller NP, Kempken F (2007) Secondary chemicals protect mould from fungivory. Biol Lett 3: 523–525. doi: 10.1098/rsbl.2007.0338 17686752

10. Vining LC (1990) Functions of secondary metabolites. Annu Rev Microbiol 44: 395–427. doi: 10.1146/annurev.mi.44.100190.002143 2252388

11. Demain AL, Fang A (2000) The natural functions of secondary metabolites. Adv Biochem Eng Biotechnol 69: 1–39. 11036689

12. Rohlfs M, Churchill ACL (2011) Fungal secondary metabolites as modulators of interactions with insects and other arthropods. Fungal Genet Biol 48: 23–34. doi: 10.1016/j.fgb.2010.08.008 20807586

13. Keller NP, Hohn TM (1997) Metabolic Pathway Gene Clusters in Filamentous Fungi. Fungal Genet Biol 21: 17–29. 9073477

14. Gardiner DM, Howlett BJ (2005) Bioinformatic and expression analysis of the putative gliotoxin biosynthetic gene cluster of Aspergillus fumigatus. FEMS Microbiol Lett 248: 241–248. doi: 10.1016/j.femsle.2005.05.046 15979823

15. Scharf DH, Heinekamp T, Remme N, Hortschansky P, Brakhage A a, et al. (2012) Biosynthesis and function of gliotoxin in Aspergillus fumigatus. Appl Microbiol Biotechnol 93: 467–472. doi: 10.1007/s00253–011–3689–1 22094977

16. Bennett J, Bentley R (1989) What’s in a name?—Microbial secondary metabolism. Adv Appl Microbiol 34. 2672724

17. Christensen M, Frisvad JC, Tuthill D (1999) Taxonomy of the Penicillium miczynskii group based on morphology and secondary metabolites. Mycol Res 103: 527–541. doi: 10.1017/S0953756298007515

18. McCowen MC, Callender ME, Lawlis JF (1951) Fumagillin (H-3), a new antibiotic with amebicidal properties. Science 113: 202–203. 14809278

19. Lin H-C, Chooi Y-H, Dhingra S, Xu W, Calvo AM, et al. (2013) The Fumagillin Biosynthetic Gene Cluster in Aspergillus fumigatus Encodes a Cryptic Terpene Cyclase Involved in the Formation of β-trans-Bergamotene. J Am Chem Soc 135: 4616–4619. doi: 10.1021/ja312503y 23488861

20. Asami Y, Kakeya H, Onose R, Chang Y-H, Toi M, et al. (2004) RK-805, an endothelial-cell-growth inhibitor produced by Neosartorya sp., and a docking model with methionine aminopeptidase-2. Tetrahedron 60: 7085–7091. doi: 10.1016/j.tet.2003.09.104

21. Wiemann P, Guo C-J, Palmer JM, Sekonyela R, Wang CCC, et al. (2013) Prototype of an intertwined secondary-metabolite supercluster. Proc Natl Acad Sci U S A 110: 17065–17070. doi: 10.1073/pnas.1313258110 24082142

22. Frisvad JC, Andersen B, Thrane U (2008) The use of secondary metabolite profiling in chemotaxonomy of filamentous fungi. Mycol Res 112: 231–240. doi: 10.1016/j.mycres.2007.08.018 18319145

23. Brakhage A a (2013) Regulation of fungal secondary metabolism. Nat Rev Microbiol 11: 21–32. doi: 10.1038/nrmicro2916 23178386

24. Bayram Ö, Krappmann S, Ni M, Bok JW, Helmstaedt K, et al. (2008) VelB/VeA/LaeA complex coordinates light signal with fungal development and secondary metabolism. Science (80-) 320: 1504–1506. doi: 10.1126/science.1155888

25. Kato N, Brooks W, Calvo AM (2003) The expression of sterigmatocystin and penicillin genes in Aspergillus nidulans is controlled by veA, a gene required for sexual development. Eukaryot Cell 2: 1178–1186. 14665453

26. Stinnett SM, Espeso EA, Cobeño L, Araújo-Bazán L, Calvo AM (2007) Aspergillus nidulans VeA subcellular localization is dependent on the importin alpha carrier and on light. Mol Microbiol 63: 242–255. doi: 10.1111/j.1365–2958.2006.05506.x 17163983

27. Purschwitz J, Müller S, Kastner C, Schöser M, Haas H, et al. (2008) Functional and physical interaction of blue- and red-light sensors in Aspergillus nidulans. Curr Biol 18: 255–259. doi: 10.1016/j.cub.2008.01.061 18291652

28. Dhingra S, Lind AL, Lin H-C, Tang Y, Rokas A, et al. (2013) The fumagillin gene cluster, an example of hundreds of genes under veA control in Aspergillus fumigatus. PLoS One 8: e77147. doi: 10.1371/journal.pone.0077147 24116213

29. Ramamoorthy V, Dhingra S, Kincaid A, Shantappa S, Feng X, et al. (2013) The Putative C2H2 Transcription Factor MtfA Is a Novel Regulator of Secondary Metabolism and Morphogenesis in Aspergillus nidulans. PLoS One 8: e74122. doi: 10.1371/journal.pone.0074122 24066102

30. Smith TD, Calvo AM (2014) The mtfA Transcription Factor Gene Controls Morphogenesis, Gliotoxin Production, and Virulence in the Opportunistic Human Pathogen Aspergillus fumigatus. Eukaryot Cell 13: 766–775. doi: 10.1128/EC.00075–14 24728192

31. Tsai HF, Chang YC, Washburn RG, Wheeler MH, Kwon-Chung KJ (1998) The developmentally regulated alb1 gene of Aspergillus fumigatus: its role in modulation of conidial morphology and virulence. J Bacteriol 180: 3031–3038. 9620950

32. Calvo AM, Wilson RA, Bok JW, Keller NP (2002) Relationship between secondary metabolism and fungal development. Microbiol Mol Biol Rev 66: 447–459, table of contents. 12208999

33. Calvo AM (2008) The VeA regulatory system and its role in morphological and chemical development in fungi. Fungal Genet Biol 45: 1053–1061. doi: 10.1016/j.fgb.2008.03.014 18457967

34. Roze L V, Chanda A, Linz JE (2011) Compartmentalization and molecular traffic in secondary metabolism: a new understanding of established cellular processes. Fungal Genet Biol 48: 35–48. doi: 10.1016/j.fgb.2010.05.006 20519149

35. Kroken S, Glass NL, Taylor JW, Yoder OC, Turgeon BG (2003) Phylogenomic analysis of type I polyketide synthase genes in pathogenic and saprobic ascomycetes. Proc Natl Acad Sci U S A 100: 15670–15675. doi: 10.1073/pnas.2532165100 14676319

36. Bayram O, Braus GH (2012) Coordination of secondary metabolism and development in fungi: the velvet family of regulatory proteins. FEMS Microbiol Rev 36: 1–24. doi: 10.1111/j.1574–6976.2011.00285.x 21658084

37. Martchenko M, Levitin A, Hogues H, Nantel A, Whiteway M (2007) Transcriptional rewiring of fungal galactose-metabolism circuitry. Curr Biol 17: 1007–1013. doi: 10.1016/j.cub.2007.05.017 17540568

38. Cain CW, Lohse MB, Homann OR, Sil A, Johnson AD (2012) A conserved transcriptional regulator governs fungal morphology in widely diverged species. Genetics 190: 511–521. doi: 10.1534/genetics.111.134080 22095082

39. Carroll SB (2005) Evolution at two levels: on genes and form. PLoS Biol 3: e245. doi: 10.1371/journal.pbio.0030245 16000021

40. Prud’homme B, Gompel N, Carroll SB (2007) Emerging principles of regulatory evolution. Proc Natl Acad Sci U S A 104 Suppl: 8605–8612. doi: 10.1073/pnas.0700488104 17494759

41. Tuch BB, Li H, Johnson AD (2008) Evolution of eukaryotic transcription circuits. Science (80-) 319: 1797–1799. doi: 10.1126/science.1152398

42. Tanay A, Regev A, Shamir R (2005) Conservation and evolvability in regulatory networks: the evolution of ribosomal regulation in yeast. Proc Natl Acad Sci U S A 102: 7203–7208. doi: 10.1073/pnas.0502521102 15883364

43. Carroll SB (2008) Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution. Cell 134: 25–36. doi: 10.1016/j.cell.2008.06.030 18614008

44. Thompson DA, Roy S, Chan M, Styczynsky MP, Pfiffner J, et al. (2013) Evolutionary principles of modular gene regulation in yeasts. Elife 2: e00603. doi: 10.7554/eLife.00603 23795289

45. Inglis DO, Binkley J, Skrzypek MS, Arnaud MB, Cerqueira GC, et al. (2013) Comprehensive annotation of secondary metabolite biosynthetic genes and gene clusters of Aspergillus nidulans, A. fumigatus, A. niger and A. oryzae. BMC Microbiol 13: 91. doi: 10.1186/1471–2180–13–91 23617571

46. Machida M, Asai K, Sano M, Tanaka T, Kumagai T, et al. (2005) Genome sequencing and analysis of Aspergillus oryzae. Nature 438: 1157–1161. doi: 10.1038/nature04300 16372010

47. Kiyota T, Hamada R, Sakamoto K, Iwashita K, Yamada O, et al. (2011) Aflatoxin non-productivity of Aspergillus oryzae caused by loss of function in the aflJ gene product. J Biosci Bioeng 111: 512–517. doi: 10.1016/j.jbiosc.2010.12.022 21342785

48. Gibbons JG, Salichos L, Slot JC, Rinker DC, McGary KL, et al. (2012) The evolutionary imprint of domestication on genome variation and function of the filamentous fungus Aspergillus oryzae. Curr Biol 22: 1403–1409. doi: 10.1016/j.cub.2012.05.033 22795693

49. Brown DW, Yu JH, Kelkar HS, Fernandes M, Nesbitt TC, et al. (1996) Twenty-five coregulated transcripts define a sterigmatocystin gene cluster in Aspergillus nidulans. Proc Natl Acad Sci U S A 93: 1418–1422. 8643646

50. Gibbons JG, Beauvais A, Beau R, McGary KL, Latgé J-P, et al. (2012) Global transcriptome changes underlying colony growth in the opportunistic human pathogen Aspergillus fumigatus. Eukaryot Cell 11: 68–78. doi: 10.1128/EC.05102–11 21724936

51. Dhingra S, Andes D, Calvo AM (2012) VeA Regulates Conidiation, Gliotoxin Production and Protease Activity in the Opportunistic Human Pathogen Aspergillus fumigatus. Eukaryot Cell 11: 1531–1543. doi: 10.1128/EC.00222–12 23087369

52. Ramamoorthy V, Shantappa S, Dhingra S, Calvo AM (2012) veA-dependent RNA-pol II transcription elongation factor-like protein, RtfA, is associated with secondary metabolism and morphological development in Aspergillus nidulans. Mol Microbiol 85: 795–814. doi: 10.1111/j.1365–2958.2012.08142.x 22783880

53. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25: 25–29. doi: 10.1038/75556 10802651

54. Arnaud MB, Chibucos MC, Costanzo MC, Crabtree J, Inglis DO, et al. (2010) The Aspergillus Genome Database, a curated comparative genomics resource for gene, protein and sequence information for the Aspergillus research community. Nucleic Acids Res 38: D420–D427. doi: 10.1093/nar/gkp751 19773420

55. Galagan JE, Calvo SE, Cuomo C, Ma L-J, Wortman JR, et al. (2005) Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature 438: 1105–1115. doi: 10.1038/nature04341 16372000

56. Khaldi N, Seifuddin FT, Turner G, Haft D, Nierman WC, et al. (2010) SMURF: Genomic mapping of fungal secondary metabolite clusters. Fungal Genet Biol 47: 736–741. doi: 10.1016/j.fgb.2010.06.003 20554054

57. Fedorova ND, Khaldi N, Joardar VS, Maiti R, Amedeo P, et al. (2008) Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genet 4: e1000046. doi: 10.1371/journal.pgen.1000046 18404212

58. Andersen MR, Salazar MP, Schaap PJ, van de Vondervoort PJI, Culley D, et al. (2011) Comparative genomics of citric-acid-producing Aspergillus niger ATCC 1015 versus enzyme-producing CBS 513.88. Genome Res 21: 885–897. doi: 10.1101/gr.112169.110 21543515

59. Hoffmeister D, Keller NP (2007) Natural products of filamentous fungi: enzymes, genes, and their regulation. Nat Prod Rep 24: 393–416. doi: 10.1039/b603084j 17390002

60. Frisvad JC, Rank C, Nielsen KF, Larsen TO (2009) Metabolomics of Aspergillus fumigatus. Med Mycol 47 Suppl 1: S53–S71. doi: 10.1080/13693780802307720 18763205

61. Geiser DM, Klich MA, Frisvad JC, Peterson SW, Varga J, et al. (2007) The current status of species recognition and identification in Aspergillus. Stud Mycol 59: 1–10. doi: 10.3114/sim.2007.59.01 18490947

62. Wisecaver JH, Slot JC, Rokas A (2014) The Evolution of Fungal Metabolic Pathways. PLoS Genet. 10: e1004816 doi: 10.1371/journal.pgen.1004816 25474404

63. Bushley KE, Turgeon BG (2010) Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships. BMC Evol Biol 10: 26. doi: 10.1186/1471–2148–10–26 20100353

64. Carvunis A-R, Rolland T, Wapinski I, Calderwood MA, Yildirim MA, et al. (2012) Proto-genes and de novo gene birth. Nature 487: 370–374. doi: 10.1038/nature11184 22722833

65. Khaldi N, Collemare J, Lebrun M-H, Wolfe KH (2008) Evidence for horizontal transfer of a secondary metabolite gene cluster between fungi. Genome Biol 9: R18. doi: 10.1186/gb-2008–9–1-r18 18218086

66. Slot JC, Rokas A (2011) Horizontal transfer of a large and highly toxic secondary metabolic gene cluster between fungi. Curr Biol 21: 134–139. doi: 10.1016/j.cub.2010.12.020 21194949

67. Campbell M a, Rokas A, Slot JC (2012) Horizontal transfer and death of a fungal secondary metabolic gene cluster. Genome Biol Evol 4: 289–293. doi: 10.1093/gbe/evs011 22294497

68. Proctor RH, Van Hove F, Susca A, Stea G, Busman M, et al. (2013) Birth, death and horizontal transfer of the fumonisin biosynthetic gene cluster during the evolutionary diversification of Fusarium. Mol Microbiol 90: 290–306. doi: 10.1111/mmi.12362 23937442

69. Khaldi N, Wolfe KH (2011) Evolutionary Origins of the Fumonisin Secondary Metabolite Gene Cluster in Fusarium verticillioides and Aspergillus niger. Int J Evol Biol 2011: 423821. doi: 10.4061/2011/423821 21716743

70. Sarikaya-Bayram O, Bayram O, Feussner K, Kim J-H, Kim H-S, et al. (2014) Membrane-bound methyltransferase complex VapA-VipC-VapB guides epigenetic control of fungal development. Dev Cell 29: 406–420. doi: 10.1016/j.devcel.2014.03.020 24871947

71. Ahmed YL, Gerke J, Park H-S, Bayram Ö, Neumann P, et al. (2013) The Velvet Family of Fungal Regulators Contains a DNA-Binding Domain Structurally Similar to NF-κB. PLoS Biol 11: e1001750. doi: 10.1371/journal.pbio.1001750 24391470

72. Tsong AE, Tuch BB, Li H, Johnson AD (2006) Evolution of alternative transcriptional circuits with identical logic. Nature 443: 415–420. doi: 10.1038/nature05099 17006507

73. Rokas A (2006) Evolution: different paths to the same end. Nature 443: 401–402. doi: 10.1038/443401a 17006502

74. Rokas A, Hittinger CT (2007) Transcriptional rewiring: the proof is in the eating. Curr Biol 17: R626–R628. doi: 10.1016/j.cub.2007.06.025 17714646

75. Galant R, Carroll SB (2002) Evolution of a transcriptional repression domain in an insect Hox protein. Nature 415: 910–913. doi: 10.1038/nature717 11859369

76. Ronshaugen M, McGinnis N, McGinnis W (2002) Hox protein mutation and macroevolution of the insect body plan. Nature 415: 914–917. doi: 10.1038/nature716 11859370

77. Tuch BB, Galgoczy DJ, Hernday AD, Li H, Johnson AD (2008) The evolution of combinatorial gene regulation in fungi. PLoS Biol 6: e38. doi: 10.1371/journal.pbio.0060038 18303948

78. Odom DT, Dowell RD, Jacobsen ES, Gordon W, Danford TW, et al. (2007) Tissue-specific transcriptional regulation has diverged significantly between human and mouse. Nat Genet 39: 730–732. doi: 10.1038/ng2047 17529977

79. Borneman AR, Gianoulis TA, Zhang ZD, Yu H, Rozowsky J, et al. (2007) Divergence of transcription factor binding sites across related yeast species. Science 317: 815–819. doi: 10.1126/science.1140748 17690298

80. O’Gorman CM, Fuller HT, Dyer PS (2009) Discovery of a sexual cycle in the opportunistic fungal pathogen Aspergillus fumigatus. Nature 457: 471–474. doi: 10.1038/nature07528 19043401

81. Crabtree J, Angiuoli S V, Wortman JR, White OR (2007) Sybil: methods and software for multiple genome comparison and visualization. Methods Mol Biol 408: 93–108. doi: 10.1007/978–1–59745–547–3_6 18314579

82. Enright AJ, Van Dongen S, Ouzounis CA (2002) An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res 30: 1575–1584. 11917018

83. Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, et al. (2009) AmiGO: online access to ontology and annotation data. Bioinformatics 25: 288–289. doi: 10.1093/bioinformatics/btn615 19033274

84. Kim H-S, Han K-Y, Kim K-J, Han D-M, Jahng K-Y, et al. (2002) The veA gene activates sexual development in Aspergillus nidulans. Fungal Genet Biol 37: 72–80. doi: 10.1016/S1087–1845(02)00029–4 12223191

85. Szewczyk E, Nayak T, Oakley CE, Edgerton H, Xiong Y, et al. (2006) Fusion PCR and gene targeting in Aspergillus nidulans. Nat Protoc 1: 3111–3120. doi: 10.1038/nprot.2006.405 17406574

86. Shaaban M, Palmer JM, El-Naggar WA, El-Sokkary MA, Habib E-SE, et al. (2010) Involvement of transposon-like elements in penicillin gene cluster regulation. Fungal Genet Biol 47: 423–432. doi: 10.1016/j.fgb.2010.02.006 20219692

87. Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics: btu170 –. doi: 10.1093/bioinformatics/btu170 24695404

88. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, et al. (2013) TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14: R36. doi: 10.1186/gb-2013–14–4-r36 23618408

89. Anders S, Pyl PT, Huber W (2014) HTSeq—A Python framework to work with high-throughput sequencing data. bioRxiv. doi: 10.1101/002824

90. Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. bioRxiv. doi: 10.1101/002832

91. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, et al. (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13: 2498–2504. doi: 10.1101/gr.1239303 14597658

92. Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21: 3448–3449. doi: 10.1093/bioinformatics/bti551 15972284

93. R Core Team (2014) R: A Language and Environment for Statistical Computing.

94. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B 57: pp. 289–300.

95. Wickham H (2009) ggplot2: elegant graphics for data analysis. Springer New York.

96. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, et al. (2009) Circos: an information aesthetic for comparative genomics. Genome Res 19: 1639–1645. doi: 10.1101/gr.092759.109 19541911

97. Schmittgen TD, Livak KJ (2008) Analyzing real-time PCR data by the comparative CT method. Nat Protoc 3: 1101–1108. doi: 10.1038/nprot.2008.73 18546601

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2015 Číslo 3
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

plice
INSIGHTS from European Respiratory Congress
nový kurz

Současné pohledy na riziko v parodontologii
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Svět praktické medicíny 3/2024 (znalostní test z časopisu)

Kardiologické projevy hypereozinofilií
Autoři: prof. MUDr. Petr Němec, Ph.D.

Střevní příprava před kolonoskopií
Autoři: MUDr. Klára Kmochová, Ph.D.

Všechny kurzy
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#