The pH-responsive PacC transcription factor of Aspergillus fumigatus governs epithelial entry and tissue invasion during pulmonary aspergillosis.

Destruction of the pulmonary epithelium is a major feature of lung diseases caused by the mould pathogen Aspergillus fumigatus. Although it is widely postulated that tissue invasion is governed by fungal proteases, A. fumigatus mutants lacking individual or multiple enzymes remain fully invasive, suggesting a concomitant requirement for other pathogenic activities during host invasion. In this study we discovered, and exploited, a novel, tissue non-invasive, phenotype in A. fumigatus mutants lacking the pH-responsive transcription factor PacC. Our study revealed a novel mode of epithelial entry, occurring in a cell wall-dependent manner prior to protease production, and via the Dectin-1 ß-glucan receptor. ΔpacC mutants are defective in both contact-mediated epithelial entry and protease expression, and significantly attenuated for pathogenicity in leukopenic mice. We combined murine infection modelling, in vivo transcriptomics, and in vitro infections of human alveolar epithelia, to delineate two major, and sequentially acting, PacC-dependent processes impacting epithelial integrity in vitro and tissue invasion in the whole animal. We demonstrate that A. fumigatus spores and germlings are internalised by epithelial cells in a contact-, actin-, cell wall- and Dectin-1 dependent manner and ΔpacC mutants, which aberrantly remodel the cell wall during germinative growth, are unable to gain entry into epithelial cells, both in vitro and in vivo. We further show that PacC acts as a global transcriptional regulator of secreted molecules during growth in the leukopenic mammalian lung, and profile the full cohort of secreted gene products expressed during invasive infection. Our study reveals a combinatorial mode of tissue entry dependent upon sequential, and mechanistically distinct, perturbations of the pulmonary epithelium and demonstrates, for the first time a protective role for Dectin-1 blockade in epithelial defences. Infecting ΔpacC mutants are hypersensitive to cell wall-active antifungal agents highlighting the value of PacC signalling as a target for antifungal therapy.

Clonality despite sex: the evolution of host-associated sexual neighborhoods in the pathogenic fungus Penicillium marneffei.

Molecular genetic approaches typically detect recombination in microbes regardless of assumed asexuality. However, genetic data have shown the AIDS-associated pathogen Penicillium marneffei to have extensive spatial genetic structure at local and regional scales, and although there has been some genetic evidence that a sexual cycle is possible, this haploid fungus is thought to be genetically, as well as morphologically, asexual in nature because of its highly clonal population structure. Here we use comparative genomics, experimental mixed-genotype infections, and population genetic data to elucidate the role of recombination in natural populations of P. marneffei. Genome wide comparisons reveal that all the genes required for meiosis are present in P. marneffei, mating type genes are arranged in a similar manner to that found in other heterothallic fungi, and there is evidence of a putatively meiosis-specific mutational process. Experiments suggest that recombination between isolates of compatible mating types may occur during mammal infection. Population genetic data from 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam show that recombination is most likely to occur across spatially and genetically limited distances in natural populations resulting in highly clonal population structure yet sexually reproducing populations. Predicted distributions of three different spatial genetic clusters within P. marneffei overlap with three different bamboo rat host distributions suggesting that recombination within hosts may act to maintain population barriers within P. marneffei.

Evidence for genetic differentiation and variable recombination rates among Dutch populations of the opportunistic human pathogen Aspergillus fumigatus.

As the frequency of antifungal drug resistance continues to increase, understanding the genetic structure of fungal populations, where resistant isolates have emerged and spread, is of major importance. Aspergillus fumigatus is an ubiquitously distributed fungus and the primary causative agent of invasive aspergillosis (IA), a potentially lethal infection in immunocompromised individuals. In the last few years, an increasing number of A. fumigatus isolates has evolved resistance to triazoles, the primary drugs for treating IA infections. In most isolates, this multiple-triazole-resistance (MTR) phenotype is caused by mutations in the cyp51A gene, which encodes the protein targeted by the triazoles. We investigated the genetic differentiation and reproductive mode of A. fumigatus in the Netherlands, the country where the MTR phenotype probably originated, to determine their role in facilitating the emergence and distribution of resistance genotypes. Using 20 genome-wide neutral markers, we genotyped 255 Dutch isolates including 25 isolates with the MTR phenotype. In contrast to previous reports, our results show that Dutch A. fumigatus genotypes are genetically differentiated into five distinct populations. Four of the five populations show significant linkage disequilibrium, indicative of an asexual reproductive mode, whereas the fifth population is in linkage equilibrium, indicative of a sexual reproductive mode. Notably, the observed genetic differentiation among Dutch isolates does not correlate with geography, although all isolates with the MTR phenotype nest within a single, predominantly asexual, population. These results suggest that both reproductive mode and genetic differentiation contribute to the structure of Dutch A. fumigatus populations and are probably shaping the evolutionary dynamics of drug resistance in this potentially deadly pathogen.

Development stage-specific proteomic profiling uncovers small, lineage specific proteins most abundant in the Aspergillus Fumigatus conidial proteome.

BACKGROUND: The pathogenic mold Aspergillus fumigatus is the most frequent infectious cause of death in severely immunocompromised individuals such as leukemia and bone marrow transplant patients. Germination of inhaled conidia (asexual spores) in the host is critical for the initiation of infection, but little is known about the underlying mechanisms of this process. RESULTS: To gain insights into early germination events and facilitate the identification of potential stage-specific biomarkers and vaccine candidates, we have used quantitative shotgun proteomics to elucidate patterns of protein abundance changes during early fungal development. Four different stages were examined: dormant conidia, isotropically expanding conidia, hyphae in which germ tube emergence has just begun, and pre-septation hyphae. To enrich for glycan-linked cell wall proteins we used an alkaline cell extraction method. Shotgun proteomic resulted in the identification of 375 unique gene products with high confidence, with no evidence for enrichment of cell wall-immobilized and secreted proteins. The most interesting discovery was the identification of 52 proteins enriched in dormant conidia including 28 proteins that have never been detected in the A. fumigatus conidial proteome such as signaling protein Pil1, chaperones BipA and calnexin, and transcription factor HapB. Additionally we found many small, Aspergillus specific proteins of unknown function including 17 hypothetical proteins. Thus, the most abundant protein, Grg1 (AFUA_5G14210), was also one of the smallest proteins detected in this study (M.W. 7,367). Among previously characterized proteins were melanin pigment and pseurotin A biosynthesis enzymes, histones H3 and H4.1, and other proteins involved in conidiation and response to oxidative or hypoxic stress. In contrast, expanding conidia, hyphae with early germ tubes, and pre-septation hyphae samples were enriched for proteins responsible for housekeeping functions, particularly translation, respiratory metabolism, amino acid and carbohydrate biosynthesis, and the tricarboxylic acid cycle. CONCLUSIONS: The observed temporal expression patterns suggest that the A. fumigatus conidia are dominated by small, lineage-specific proteins. Some of them may play key roles in host-pathogen interactions, signal transduction during conidial germination, or survival in hostile environments.

Genome sequencing and analysis of Aspergillus oryzae.

The genome of Aspergillus oryzae, a fungus important for the production of traditional fermented foods and beverages in Japan, has been sequenced. The ability to secrete large amounts of proteins and the development of a transformation system have facilitated the use of A. oryzae in modern biotechnology. Although both A. oryzae and Aspergillus flavus belong to the section Flavi of the subgenus Circumdati of Aspergillus, A. oryzae, unlike A. flavus, does not produce aflatoxin, and its long history of use in the food industry has proved its safety. Here we show that the 37-megabase (Mb) genome of A. oryzae contains 12,074 genes and is expanded by 7-9 Mb in comparison with the genomes of Aspergillus nidulans and Aspergillus fumigatus. Comparison of the three aspergilli species revealed the presence of syntenic blocks and A. oryzae-specific blocks (lacking synteny with A. nidulans and A. fumigatus) in a mosaic manner throughout the genome of A. oryzae. The blocks of A. oryzae-specific sequence are enriched for genes involved in metabolism, particularly those for the synthesis of secondary metabolites. Specific expansion of genes for secretory hydrolytic enzymes, amino acid metabolism and amino acid/sugar uptake transporters supports the idea that A. oryzae is an ideal microorganism for fermentation.

Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus.

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility (het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer (HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated "gene dumps" and, perhaps, simultaneously, as "gene factories".

SMURF: Genomic mapping of fungal secondary metabolite clusters.

Fungi produce an impressive array of secondary metabolites (SMs) including mycotoxins, antibiotics and pharmaceuticals. The genes responsible for their biosynthesis, export, and transcriptional regulation are often found in contiguous gene clusters. To facilitate annotation of these clusters in sequenced fungal genomes, we developed the web-based software SMURF (www.jcvi.org/smurf/) to systematically predict clustered SM genes based on their genomic context and domain content. We applied SMURF to catalog putative clusters in 27 publicly available fungal genomes. Comparison with genetically characterized clusters from six fungal species showed that SMURF accurately recovered all clusters and detected additional potential clusters. Subsequent comparative analysis revealed the striking biosynthetic capacity and variability of the fungal SM pathways and the correlation between unicellularity and the absence of SMs. Further genetics studies are needed to experimentally confirm these clusters.

Sub-telomere directed gene expression during initiation of invasive aspergillosis.

Aspergillus fumigatus is a common mould whose spores are a component of the normal airborne flora. Immune dysfunction permits developmental growth of inhaled spores in the human lung causing aspergillosis, a significant threat to human health in the form of allergic, and life-threatening invasive infections. The success of A. fumigatus as a pathogen is unique among close phylogenetic relatives and is poorly characterised at the molecular level. Recent genome sequencing of several Aspergillus species provides an exceptional opportunity to analyse fungal virulence attributes within a genomic and evolutionary context. To identify genes preferentially expressed during adaptation to the mammalian host niche, we generated multiple gene expression profiles from minute samplings of A. fumigatus germlings during initiation of murine infection. They reveal a highly co-ordinated A. fumigatus gene expression programme, governing metabolic and physiological adaptation, which allows the organism to prosper within the mammalian niche. As functions of phylogenetic conservation and genetic locus, 28% and 30%, respectively, of the A. fumigatus subtelomeric and lineage-specific gene repertoires are induced relative to laboratory culture, and physically clustered genes including loci directing pseurotin, gliotoxin and siderophore biosyntheses are a prominent feature. Locationally biased A. fumigatus gene expression is not prompted by in vitro iron limitation, acid, alkaline, anaerobic or oxidative stress. However, subtelomeric gene expression is favoured following ex vivo neutrophil exposure and in comparative analyses of richly and poorly nourished laboratory cultured germlings. We found remarkable concordance between the A. fumigatus host-adaptation transcriptome and those resulting from in vitro iron depletion, alkaline shift, nitrogen starvation and loss of the methyltransferase LaeA. This first transcriptional snapshot of a fungal genome during initiation of mammalian infection provides the global perspective required to direct much-needed diagnostic and therapeutic strategies and reveals genome organisation and subtelomeric diversity as potential driving forces in the evolution of pathogenicity in the genus Aspergillus.

Transcriptional regulation of chemical diversity in Aspergillus fumigatus by LaeA.

Secondary metabolites, including toxins and melanins, have been implicated as virulence attributes in invasive aspergillosis. Although not definitively proved, this supposition is supported by the decreased virulence of an Aspergillus fumigatus strain, DeltalaeA, that is crippled in the production of numerous secondary metabolites. However, loss of a single LaeA-regulated toxin, gliotoxin, did not recapitulate the hypovirulent DeltalaeA pathotype, thus implicating other toxins whose production is governed by LaeA. Toward this end, a whole-genome comparison of the transcriptional profile of wild-type, DeltalaeA, and complemented control strains showed that genes in 13 of 22 secondary metabolite gene clusters, including several A. fumigatus-specific mycotoxin clusters, were expressed at significantly lower levels in the DeltalaeA mutant. LaeA influences the expression of at least 9.5% of the genome (943 of 9,626 genes in A. fumigatus) but positively controls expression of 20% to 40% of major classes of secondary metabolite biosynthesis genes such as nonribosomal peptide synthetases (NRPSs), polyketide synthases, and P450 monooxygenases. Tight regulation of NRPS-encoding genes was highlighted by quantitative real-time reverse-transcription PCR analysis. In addition, expression of a putative siderophore biosynthesis NRPS (NRPS2/sidE) was greatly reduced in the DeltalaeA mutant in comparison to controls under inducing iron-deficient conditions. Comparative genomic analysis showed that A. fumigatus secondary metabolite gene clusters constitute evolutionarily diverse regions that may be important for niche adaptation and virulence attributes. Our findings suggest that LaeA is a novel target for comprehensive modification of chemical diversity and pathogenicity.

Using aCGH to study intraspecific genetic variability in two pathogenic molds, Aspergillus fumigatus and Aspergillus flavus.

We have examined the feasibility of using array comparative genomic hybridization (aCGH) to explore intraspecific genetic variability at the genomic level in two pathogenic molds, Aspergillus fumigatus and Aspergillus flavus. Our analysis showed that strain-specific genes may comprise up to 2% of their genomes in comparison to isolates from different vegetative (heterokaryon) compatibility groups (VCGs). In contrast, isolates with the same VCG affiliations have almost identical gene content. Most isolate-specific genes are annotated as 'hypothetical' and located in a few large subtelomeric indels. The list includes highly polymorphic loci that contain putative het (heterokaryon compatibility) loci, which determine the individual's VCG during parasexual crossing. Incidentally, VCGs in both species seem to be significantly associated with either alpha or HMG mating type (Chi-square test, P=0.05). In conclusion CGH can be used to effectively to identify isolate-specific genes in Aspergillus species. Preliminary evidence suggests that gene flow in both species is largely constrained by VCG boundaries, although further VCG sampling is required to confirm this observation.

CDD: a curated Entrez database of conserved domain alignments.

The Conserved Domain Database (CDD) is now indexed as a separate database within the Entrez system and linked to other Entrez databases such as MEDLINE(R). This allows users to search for domain types by name, for example, or to view the domain architecture of any protein in Entrez's sequence database. CDD can be accessed on the WorldWideWeb at http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=cdd. Users may also employ the CD-Search service to identify conserved domains in new sequences, at http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi. CD-Search results, and pre-computed links from Entrez's protein database, are calculated using the RPS-BLAST algorithm and Position Specific Score Matrices (PSSMs) derived from CDD alignments. CD-Searches are also run by default for protein-protein queries submitted to BLAST(R) at http://www.ncbi.nlm.nih.gov/BLAST. CDD mirrors the publicly available domain alignment collections SMART and PFAM, and now also contains alignment models curated at NCBI. Structure information is used to identify the core substructure likely to be present in all family members, and to produce sequence alignments consistent with structure conservation. This alignment model allows NCBI curators to annotate 'columns' corresponding to functional sites conserved among family members.

MMDB: Entrez's 3D-structure database.

Three-dimensional structures are now known within most protein families and it is likely, when searching a sequence database, that one will identify a homolog of known structure. The goal of Entrez's 3D-structure database is to make structure information and the functional annotation it can provide easily accessible to molecular biologists. To this end, Entrez's search engine provides several powerful features: (i) links between databases, for example between a protein's sequence and structure; (ii) pre-computed sequence and structure neighbors; and (iii) structure and sequence/structure alignment visualization. Here, we focus on a new feature of Entrez's Molecular Modeling Database (MMDB): Graphical summaries of the biological annotation available for each 3D structure, based on the results of automated comparative analysis. MMDB is available at: http://www.ncbi.nlm.nih.gov/Entrez/structure.html.

Comparative analysis of programmed cell death pathways in filamentous fungi.

BACKGROUND: Fungi can undergo autophagic- or apoptotic-type programmed cell death (PCD) on exposure to antifungal agents, developmental signals, and stress factors. Filamentous fungi can also exhibit a form of cell death called heterokaryon incompatibility (HI) triggered by fusion between two genetically incompatible individuals. With the availability of recently sequenced genomes of Aspergillus fumigatus and several related species, we were able to define putative components of fungi-specific death pathways and the ancestral core apoptotic machinery shared by all fungi and metazoa. RESULTS: Phylogenetic profiling of HI-associated proteins from four Aspergilli and seven other fungal species revealed lineage-specific protein families, orphan genes, and core genes conserved across all fungi and metazoa. The Aspergilli-specific domain architectures include NACHT family NTPases, which may function as key integrators of stress and nutrient availability signals. They are often found fused to putative effector domains such as Pfs, SesB/LipA, and a newly identified domain, HET-s/LopB. Many putative HI inducers and mediators are specific to filamentous fungi and not found in unicellular yeasts. In addition to their role in HI, several of them appear to be involved in regulation of cell cycle, development and sexual differentiation. Finally, the Aspergilli possess many putative downstream components of the mammalian apoptotic machinery including several proteins not found in the model yeast, Saccharomyces cerevisiae. CONCLUSION: Our analysis identified more than 100 putative PCD associated genes in the Aspergilli, which may help expand the range of currently available treatments for aspergillosis and other invasive fungal diseases. The list includes species-specific protein families as well as conserved core components of the ancestral PCD machinery shared by fungi and metazoa.

Genomics of Aspergillus fumigatus.

Aspergillus fumigatus is a filamentous fungal saprophyte that is ubiquitous in the environment. It is also a human pathogen and induces allergenic response, negatively impacting health care and associated costs significantly around the world. Much of the basic biology of this organism is only poorly understood, but the recent completion and publication of its genome sequence provides an excellent tool for researchers to gain insight into these processes. In this review we will summarize some of the more salient features revealed by analysis of the genome, including the search for candidate pathogenicity genes and the switch to a pathogenic lifestyle, allergen proteins, DNA repair, secondary metabolite gene clusters that produce compounds both useful and toxic, a theoretical capability of this asexual organism to reproduce sexually, signalling, and transcription. A. fumigatus was compared with the food biotechnology fungus Aspergillus oryzae and sexual fungus Aspergillus nidulans, as well as other fungi, in an attempt to discern key differences between these organisms.

The COG database: an updated version includes eukaryotes.

BACKGROUND: The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies. RESULTS: We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or approximately 54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of approximately 20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (approximately 1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes. CONCLUSION: The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.

Mobile elements and mitochondrial genome expansion in the soil fungus and potato pathogen Rhizoctonia solani AG-3.

The soil fungus Rhizoctonia solani is an economically important pathogen of agricultural and forestry crops. Here, we present the complete sequence and analysis of the mitochondrial genome of R. solani, field isolate Rhs1AP. The genome (235 849 bp) is the largest mitochondrial genome of a filamentous fungus sequenced to date and exhibits a rich accumulation of introns, novel repeat sequences, homing endonuclease genes, and hypothetical genes. Stable secondary structures exhibited by repeat sequences suggest that they comprise functional, possibly catalytic RNA elements. RNA-Seq expression profiling confirmed that the majority of homing endonuclease genes and hypothetical genes are transcriptionally active. Comparative analysis suggests that the mitochondrial genome of R. solani is an example of a dynamic history of expansion in filamentous fungi.

Bioinformatics approaches and software for detection of secondary metabolic gene clusters.

The accelerating pace of microbial genomics is sparking a renaissance in the field of natural products research. Researchers can now get a preview of the organism's secondary metabolome by analyzing its genomic sequence. Combined with other -omics data, this approach may provide a cost-effective alternative to industrial high-throughput screening in drug discovery. In the last few years, several computational tools have been developed to facilitate this process by identifying genes involved in secondary metabolite biosynthesis in bacterial and fungal genomes. Here, we review seven software programs that are available for this purpose, with an emphasis on antibiotics & Secondary Metabolite Analysis SHell (antiSMASH) and Secondary Metabolite Unknown Regions Finder (SMURF), the only tools that can comprehensively detect complete secondary metabolite biosynthesis gene clusters. We also discuss five related software packages-CLUster SEquence ANalyzer (CLUSEAN), ClustScan, Structure Based Sequence Analysis of Polyketide Synthases (SBSPKS), NRPSPredictor, and Natural Product searcher (NP.searcher)-that identify secondary metabolite backbone biosynthesis genes. This chapter offers detailed protocols, suggestions, and caveats to assist researchers in using these tools most effectively.

Triallelic SNP-mediated genotyping of regenerated protoplasts of the heterokaryotic fungus Rhizoctonia solani.

The aneuploid and heterokaryotic nuclear condition of the soil fungus Rhizoctonia solani have provided challenges in obtaining a complete genome sequence. To better aid in the assembly and annotation process, a protoplast and single nucleotide polymorphism (SNP)-based method was developed to identify regenerated protoplasts with a reduced nuclear genome. Protocol optimization experiments showed that enzymatic digestion of mycelium from a 24 h culture of R. solani increased the proportion of protoplasts with a diameter of ≤7.5 µm and 1-4 nuclei. To determine whether strains regenerated from protoplasts with a reduced number of nuclei were genetically different from the parental strain, triallelic SNPs identified from variance records of the genomic DNA sequence reads of R. solani were used in PCR-based genotyping assays. Results from 16 of the 24 SNP-based PCR assays provided evidence that one of the three alleles was missing in the 11 regenerated protoplast strains, suggesting that these strains represent a reduced genomic complement of the parental strain. The protoplast and triallelic SNP-based method used in this study may be useful in strain development and analysis of other basidiomycete fungi with complex nuclear genomes.

Tight control of mycotoxin biosynthesis gene expression in Aspergillus flavus by temperature as revealed by RNA-Seq.

To better understand the effect of temperature on mycotoxin biosynthesis, RNA-Seq technology was used to profile the Aspergillus flavus transcriptome under different temperature conditions. This approach allowed us to quantify transcript abundance for over 80% of fungal genes including 1153 genes that were differentially expressed at 30 and 37 °C. Eleven of the 55 secondary metabolite clusters were upregulated at the lower temperature, including aflatoxin biosynthesis genes, which were among the most highly upexpressed genes. On average, transcript abundance for the 30 aflatoxin biosynthesis genes was 3300 times greater at 30 °C as compared with 37 °C. The results are consistent with the view that high temperature negatively affects aflatoxin production by turning down transcription of the two key transcriptional regulators, aflR and aflS. Subtle changes in the expression levels of aflS to aflR appear to control transcription activation of the aflatoxin cluster.