Molecular genetic analysis and regulation of aflatoxin biosynthesis

Aflatoxins, produced by some Aspergillus species, are toxic and extremely carcinogenic furanocoumarins. Recent investigations of the molecular mechanism of AFB biosynthesis showed that the genes required for biosynthesis are in a 70 kb gene cluster. They encode a DNA-binding protein functioning in aflatoxin pathway gene regulation, and other enzymes such as cytochrome p450-type monooxygenases, dehydrogenases, methyltransferases, and polyketide and fatty acid synthases. Information gained from these studies has led to a better understanding of aflatoxin biosynthesis by these fungi. The characterization of genes involved in aflatoxin formation affords the opportunity to examine the mechanism of molecular regulation of the aflatoxin biosynthetic pathway, particularly during the interaction between aflatoxin-producing fungi and plants.     

Cloning and characterization of the O-methyltransferase I gene (dmtA) from Aspergillus parasiticus associated with the conversions of demethylsterigmatocystin to sterigmatocystin and dihydrodemethylsterigmatocystin to dihydrosterigmatocystin in aflatoxin biosynthesis.

O-Methyltransferase I catalyzes both the conversion of demethylsterigmatocystin to sterigmatocystin and the conversion of dihydrodemethylsterigmatocystin to dihydrosterigmatocystin during aflatoxin biosynthesis. In this study, both genomic cloning and cDNA cloning of the gene encoding O-methyltransferase I were accomplished by using PCR strategies, such as conventional PCR based on the N-terminal amino acid sequence of the purified enzyme, 5' and 3' rapid amplification of cDNA ends PCR, and thermal asymmetric interlaced PCR (TAIL-PCR), and genes were sequenced by using Aspergillus parasiticus NIAH-26. A comparison of the genomic sequences with the cDNA of the dmtA region revealed that the coding region is interrupted by three short introns. The cDNA of the dmtA gene is 1,373 bp long and encodes a 386-amino-acid protein with a deduced molecular weight of 43,023, which is consistent with the molecular weight of the protein determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The C-terminal half of the deduced protein exhibits 76.3% identity with the coding region of the Aspergillus nidulans StcP protein, whereas the N-terminal half of dmtA exhibits 73.0% identity with the 5' flanking region of the stcP gene, suggesting that translation of the stcP gene may start at a site upstream from methionine that is different from the site that has been suggested previously. Also, an examination of the 5' and 3' flanking regions of the dmtA gene in which TAIL-PCR was used demonstrated that the dmtA gene is located in the aflatoxin biosynthesis cluster between (and in the same orientation as) the omtA and ord-2 genes. Northern blotting revealed that expression of the dmtA gene is influenced by both medium composition and culture temperature and that the pattern correlates with the patterns observed for other genes in the aflatoxin gene cluster. Furthermore, Southern blotting and PCR analyses of the dmtA gene showed that a dmtA homolog is present in Aspergillus oryzae SYS-2.     

黄曲霉毒素生物合成径调节基因aflR

黄曲霉毒素生物合成途径调节基因在黄曲霉毒素产生过程中发挥十分重要的作用,它为绝大多数黄曲霉毒素合成相关基因的表达所必需。黄曲霉毒素生物合成途径调节基因的启动子中,含有若干真菌转录因子同源物的假定结合位点。AflR蛋白是黄曲霉毒素生物合成途径中的主要正性转录因子,它调节大多数黄曲霉毒素合成相关基因,也包括其自身基因的表达。     

Cloning and Characterization of the O-Methyltransferase I Gene (dmtA) from Aspergillus parasiticus Associated with the Conversions of Demethylsterigmatocystin to Sterigmatocystin and Dihydrodemethylsterigmatocystin to Dihydrosterigmatocystin in Aflatoxin Biosynthesis

O-Methyltransferase I catalyzes both the conversion of demethylsterigmatocystin to sterigmatocystin and the conversion of dihydrodemethylsterigmatocystin to dihydrosterigmatocystin during aflatoxin biosynthesis. In this study, both genomic cloning and cDNA cloning of the gene encoding O-methyltransferase I were accomplished by using PCR strategies, such as conventional PCR based on the N-terminal amino acid sequence of the purified enzyme, 5′ and 3′ rapid amplification of cDNA ends PCR, and thermal asymmetric interlaced PCR (TAIL-PCR), and genes were sequenced by using Aspergillus parasiticus NIAH-26. A comparison of the genomic sequences with the cDNA of the dmtA region revealed that the coding region is interrupted by three short introns. The cDNA of the dmtA gene is 1,373 bp long and encodes a 386-amino-acid protein with a deduced molecular weight of 43,023, which is consistent with the molecular weight of the protein determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The C-terminal half of the deduced protein exhibits 76.3% identity with the coding region of the Aspergillus nidulans StcP protein, whereas the N-terminal half of dmtA exhibits 73.0% identity with the 5′ flanking region of the stcP gene, suggesting that translation of the stcP gene may start at a site upstream from methionine that is different from the site that has been suggested previously. Also, an examination of the 5′ and 3′ flanking regions of the dmtA gene in which TAIL-PCR was used demonstrated that the dmtA gene is located in the aflatoxin biosynthesis cluster between (and in the same orientation as) the omtA and ord-2 genes. Northern blotting revealed that expression of the dmtA gene is influenced by both medium composition and culture temperature and that the pattern correlates with the patterns observed for other genes in the aflatoxin gene cluster. Furthermore, Southern blotting and PCR analyses of the dmtA gene showed that a dmtA homolog is present in Aspergillus oryzae SYS-2.     

Effect of Dietary Aflatoxin Concentration on the Assessment of Genetic Variability of Response in a Randombred Population of Chickens

The effect of graded levels of dietary aflatoxin on the assessment of genetic variability of body weight and gain and plasma protein response was tested utilizing the Athens-Canadian randombred population of chickens. Dietary aflatoxin was administered at levels of either 0, 1.25, 2.50 or 5.0 µg/g of diet ad libitum from 7 to 21 days of age to progeny from 58 sire families. Twenty-one-day body weights, gain and plasma protein concentration were used to assess the variation in response.—The administration of increasing levels of aflatoxin resulted in a dose-related decrease of gains and plasma protein concentrations. Plasma protein concentrations were significantly higher among males than females within the control group; however, this difference was reversed as the severity of the aflatoxin challenge increased. Heritability estimates for all responses increased as the level of aflatoxin administered increased. This change was most notable for total plasma protein concentration. Phenotypic correlations for plasma protein concentration and growth measurements tended to diminish with increasing levels of aflatoxin. A similar trend was noted for the genetic correlations; however, a moderate correlation between growth responses and plasma protein response was detected in the 5.0-µg/g aflatoxin treatment group. Genetic correlations were calculated for the same characters between the different levels of aflatoxin. Regardless of which aflatoxin challenges were compared, a very high genetic correlation for 21-day body weight and 7- to 21-day gain was estimated. This variation in growth potential in the toxic environment paralleled that observed in the control environment but at a lower plane. Genetic correlations for plasma protein response across aflatoxin levels diminished as the difference between the levels of aflatoxin administered increased. Plasma protein concentration in the control environment was positively correlated with plasma protein response in groups fed a low level of aflatoxin, but negatively correlated when an aflatoxin challenge of 2.5 µg/g or more was given, suggesting that selection for aflatoxin resistance using plasma protein response as a selection criterion should be made under an aflatoxin stress environment.     

Aspergillus flavus expressed sequence tags and microarray as tools in understanding aflatoxin biosynthesis

Aflatoxins are the most toxic and carcinogenic naturally occurring mycotoxins. They are produced primarily byAspergillus flavus andA. parasiticus. In order to better understand the molecular mechanisms that control aflatoxin production, identification of genes usingA. flavus expressed sequence tags (ESTs) and microarrays is currently being performed. Sequencing and annotation ofA. flavus ESTs from a normalizedA. flavus cDNA library identified 7,218 unique EST sequences. Genes that are putatively involved in aflatoxin biosynthesis, regulation and signal transduction, fungal virulence or pathogenicity, stress response or antioxidation, and fungal development were identified from these ESTs. Microarrays containing over 5,000 uniqueA. flavus gene amplicons were constructed at The Institute for Genomic Research. Gene expression profiling under aflatoxin-producing and non-producing conditions using this microarray has identified hundreds of genes that are potentially involved in aflatoxin production. Further investigations on the functions of these genes by gene knockout experiments are underway. This research is expected to provide information for developing new strategies for controlling aflatoxin contamination of agricultural commodities.     

Participation in aflatoxin biosynthesis by a reductase enzyme encoded by vrdA gene outside the aflatoxin gene cluster

Three reactions from hydroxyversicolorone to versicolorone, from versiconal hemiacetal acetate to versiconol acetate, and from versiconal to versiconol are involved in a metabolic grid in aflatoxin biosynthesis. This work demonstrated that the same reductase of Aspergillus parasiticus catalyzes the three reactions. The gene (named vrdA) encoding the reductase was cloned, and its sequence did not show homology to any regions in aflatoxin gene cluster. Its cDNA encoding a 38,566 Da protein was separated by three introns in the genome. Deletion of the vrdA gene in A. parasiticus caused a significant decrease in enzyme activity, but did not affect aflatoxin productivity of the fungi. Although the vrdA gene was expressed in culture conditions conducive to aflatoxin production, it was expressed even in the aflR deletion mutant. These results suggest that the vrdA is not an aflatoxin biosynthesis gene, although it actually participates in aflatoxin biosynthesis in cells.     

Sequence comparison of aflR from different Aspergillus species provides evidence for variability in regulation of aflatoxin production

Aflatoxin contamination of foods and feeds is a world-wide agricultural problem. Aflatoxin production requires expression of the biosynthetic pathway regulatory gene, aflR, which encodes a Cys6Zn2-type DNA-binding protein. Homologs of aflR from Aspergillus nomius, bombycis, parasiticus, flavus, and pseudotamarii were compared to investigate the molecular basis for variation among aflatoxin-producing taxa in the regulation of aflatoxin production. Variability was found in putative promoter consensus elements and coding region motifs, including motifs involved in developmental regulation (AbaA, BrlA), regulation of nitrogen source utilization (AreA), and pH regulation (PacC), and in coding region PEST domains. Some of these elements may affect expression of aflJ, a gene divergently transcribed from aflR, that also is required for aflatoxin accumulation. Comparisons of phylogenetic trees obtained with either aligned aflR intergenic region sequence or coding region sequence and the observed divergence in regulatory features among the taxa provide evidence that regulatory signals for aflatoxin production evolved to respond to a variety of environmental stimuli under differential selective pressures. Phylogenetic analyses also suggest that isolates currently assigned to the A. flavus morphotype SBG represent a distinct species and that A. nomius is a diverse paraphyletic assemblage likely to contain several species.     

Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.     

Markov Chain modeling of pyelonephritis-associated pili expression in uropathogenic Escherichia coli

Pyelonephritis-associated pili (Pap) expression in uropathogenic Escherichia coli is regulated by a complex phase variation mechanism involving the competition between leucine-responsive regulatory protein (Lrp) and DNA adenine methylase (Dam). Population dynamics of pap gene expression has been studied extensively and the detailed molecular mechanism has been largely elucidated, providing sufficient information for mathematical modeling. Although the Gillespie algorithm is suited for modeling of stochastic systems such as the pap operon, it becomes computationally expensive when detailed molecular steps are explicitly modeled in a population. Here we developed a Markov Chain model to simplify the computation. Our model is analytically derived from the molecular mechanism. The model presented here is able to reproduce results presented using the Gillespie method, but since the regulatory information is incorporated before simulation, our model runs more efficiently and allows investigation of additional regulatory features. The model predictions are consistent with experimental data obtained in this work and in the literature. The results show that pap expression in uropathogenic E. coli is initial-state-dependent, as previously reported. However, without environment stimuli, the pap-expressing fraction in a population will reach an equilibrium level after approximately 50-100 generations. The transient time before reaching equilibrium is determined by PapI stability and Lrp and Dam copy numbers per cell. This work demonstrates that the Markov Chain model captures the essence of the complex molecular mechanism and greatly simplifies the computation.     

veA is required for toxin and sclerotial production in Aspergillus parasiticus

It was long been noted that secondary metabolism is associated with fungal development. In Aspergillus nidulans, conidiation and mycotoxin production are linked by a G protein signaling pathway. Also in A. nidulans, cleistothecial development and mycotoxin production are controlled by a gene called veA. Here we report the characterization of a veA ortholog in the aflatoxin-producing fungus A. parasiticus. Cleistothecia are not produced by Aspergillus parasiticus; instead, this fungus produces spherical structures called sclerotia that allow for survival under adverse conditions. Deletion of veA from A. parasiticus resulted in the blockage of sclerotial formation as well as a blockage in the production of aflatoxin intermediates. Our results indicate that A. parasiticus veA is required for the expression of aflR and aflJ, which regulate the activation of the aflatoxin gene cluster. In addition to these findings, we observed that deletion of veA reduced conidiation both on the culture medium and on peanut seed. The fact that veA is necessary for conidiation, production of resistant structures, and aflatoxin biosynthesis makes veA a good candidate gene to control aflatoxin biosynthesis or fungal development and in this way to greatly decrease its devastating impact on health and the economy.     

Cloning of the afl-2 gene involved in aflatoxin biosynthesis from Aspergillus flavus.

Aflatoxins are extremely potent carcinogens produced by Aspergillus flavus and Aspergillus parasiticus. Cloning of genes in the aflatoxin pathway provides a specific approach to understanding the regulation of aflatoxin biosynthesis and, subsequently, to the control of aflatoxin contamination of food and feed. This paper reports the isolation of a gene involved in aflatoxin biosynthesis by complementation of an aflatoxin-nonproducing mutant with a wild-type genomic cosmid library of A. flavus. Strain 650-33, blocked in aflatoxin biosynthesis at the afl-2 allele, was complemented by a 32-kb cosmid clone (B9), resulting in the production of aflatoxin. The onset and profile of aflatoxin accumulation was similar for the transformed strain and the wild-type strain (NRRL 3357) of the fungus, indicating that the integrated gene is under the same control as in wild-type strains. Complementation analyses with DNA fragments from B9 indicated that the gene resides within a 2.2-kb fragment. Because this gene complements the mutated afl-2 allele, it was designated afl-2. Genetic evidence obtained from a double mutant showed that afl-2 is involved in aflatoxin biosynthesis before the formation of norsolorinic acid, the first stable intermediate identified in the pathway. Further, metabolite feeding studies with the mutant, transformed, and wild-type cultures and enzymatic activity measurements in cell extracts of these cultures suggest that afl-2 regulates gene expression or the activity of other aflatoxin pathway enzymes. This is the first reported isolation of a gene for aflatoxin biosynthesis in A. flavus.