Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum

While studying myoblast methylomes and transcriptomes, we found that CDH15 had a remarkable preference for expression in both myoblasts and cerebellum. To understand how widespread such a relationship was and its epigenetic and biological correlates, we systematically looked for genes with similar transcription profiles and analyzed their DNA methylation and chromatin state and accessibility profiles in many different cell populations. Twenty genes were expressed preferentially in myoblasts and cerebellum (Myob/Cbl genes). Some shared DNA hypo- or hypermethylated regions in myoblasts and cerebellum. Particularly striking was ZNF556, whose promoter is hypomethylated in expressing cells but highly methylated in the many cell populations that do not express the gene. In reporter gene assays, we demonstrated that its promoter's activity is methylation sensitive. The atypical epigenetics of ZNF556 may have originated from its promoter's hypomethylation and selective activation in sperm progenitors and oocytes. Five of the Myob/Cbl genes (KCNJ12, ST8SIA5, ZIC1, VAX2, and EN2) have much higher RNA levels in cerebellum than in myoblasts and displayed myoblast-specific hypermethylation upstream and/or downstream of their promoters that may downmodulate expression. Differential DNA methylation was associated with alternative promoter usage for Myob/Cbl genes MCF2L, DOK7, CNPY1, and ANK1. Myob/Cbl genes PAX3, LBX1, ZNF556, ZIC1, EN2, and VAX2 encode sequence-specific transcription factors, which likely help drive the myoblast and cerebellum specificity of other Myob/Cbl genes. This study extends our understanding of epigenetic/transcription associations related to differentiation and may help elucidate relationships between epigenetic signatures and muscular dystrophies or cerebellar-linked neuropathologies.

Epigenetics of Mitochondria-Associated Genes in Striated Muscle

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3' end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson's disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene's four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

Epigenetics and expression of key genes associated with cardiac fibrosis: NLRP3, MMP2, MMP9, CCN2/CTGF and AGT

Aims:

Excessive inflammatory signaling and pathological remodeling of the extracellular matrix drive cardiac fibrosis and require changes in gene expression.

Materials and methods:

Using bioinformatics, both tissue-specific expression profiles and epigenomic profiles of some genes critical for cardiac fibrosis were examined, namely, NLRP3, MMP2, MMP9, CCN2/CTGF, AGT (encodes angiotensin II precursors) and hsa-mir-223 (post-transcriptionally regulates NLRP3).

Results:

In monocytes, neutrophils, fibroblasts, venous cells, liver and brain, enhancers or super-enhancers were found that correlate with high expression of these genes. One enhancer extended into a silent gene neighbor. These enhancers harbored tissue-specific foci of DNA hypomethylation, open chromatin and transcription factor binding.

Conclusions:

This study identified previously undescribed enhancers containing hypomethylated transcription factor binding subregions that are predicted to regulate expression of these cardiac fibrosis-inducing genes.

Epigenetics of Muscle- and Brain-Specific Expression of KLHL Family Genes

KLHL and the related KBTBD genes encode components of the Cullin-E3 ubiquitin ligase complex and typically target tissue-specific proteins for degradation, thereby affecting differentiation, homeostasis, metabolism, cell signaling, and the oxidative stress response. Despite their importance in cell function and disease (especially, KLHL40, KLHL41, KBTBD13, KEAP1, and ENC1), previous studies of epigenetic factors that affect transcription were predominantly limited to promoter DNA methylation. Using diverse tissue and cell culture whole-genome profiles, we examined 17 KLHL or KBTBD genes preferentially expressed in skeletal muscle or brain to identify tissue-specific enhancer and promoter chromatin, open chromatin (DNaseI hypersensitivity), and DNA hypomethylation. Sixteen of the 17 genes displayed muscle- or brain-specific enhancer chromatin in their gene bodies, and most exhibited specific intergenic enhancer chromatin as well. Seven genes were embedded in super-enhancers (particularly strong, tissue-specific clusters of enhancers). The enhancer chromatin regions typically displayed foci of DNA hypomethylation at peaks of open chromatin. In addition, we found evidence for an intragenic enhancer in one gene upregulating expression of its neighboring gene, specifically for KLHL40/HHATL and KLHL38/FBXO32 gene pairs. Many KLHL/KBTBD genes had tissue-specific promoter chromatin at their 5' ends, but surprisingly, two (KBTBD11 and KLHL31) had constitutively unmethylated promoter chromatin in their 3' exons that overlaps a retrotransposed KLHL gene. Our findings demonstrate the importance of expanding epigenetic analyses beyond the 5' ends of genes in studies of normal and abnormal gene regulation.

Mycotoxins in Grain Dust: Method for Analysis of Aflatoxins, Ochratoxin A, Zearalenone,Vomitoxin, and Secalonic Acid D

A multimycotoxin method is presented to quantitate aflatoxins, ochratoxin A, zearalenone, secalonic acid D, and vomitoxin in grain dust. Dust spiked with these mycotoxins was extracted sequentially with methylene chloride followed by acetonitrile-water (86 + 14). Vomitoxin was recovered in the latter extract and all other mycotoxins were recovered in the methylene chloride. Aflatoxins and ochratoxin were quantitated by fluorescence measurement on silica thin layer chromatographic plates. The other mycotoxins were quantitated after cleanup by reverse phase liquid chromatography and ultraviolet detection. Recoveries from dust spiked in the parts per billion (ng/g) range were approximately 80% (SD = 15-29%) for all mycotoxins. Minimum detectable amounts ranged from less than 0.5 ng/g for aflatoxins to 20 ng/g for zearalenone.

Epigenetics of Skeletal Muscle-Associated Genes in the ASB, LRRC, TMEM, and OSBPL Gene Families

Much remains to be discovered about the intersection of tissue-specific transcription control and the epigenetics of skeletal muscle (SkM), a very complex and dynamic organ. From four gene families, Leucine-Rich Repeat Containing (LRRC), Oxysterol Binding Protein Like (OSBPL), Ankyrin Repeat and Socs Box (ASB), and Transmembrane Protein (TMEM), we chose 21 genes that are preferentially expressed in human SkM relative to 52 other tissue types and analyzed relationships between their tissue-specific epigenetics and expression. We also compared their genetics, proteomics, and descriptions in the literature. For this study, we identified genes with little or no previous descriptions of SkM functionality (ASB4, ASB8, ASB10, ASB12, ASB16, LRRC14B, LRRC20, LRRC30, TMEM52, TMEM233, OSBPL6/ORP6, and OSBPL11/ORP11) and included genes whose SkM functions had been previously addressed (ASB2, ASB5, ASB11, ASB15, LRRC2, LRRC38, LRRC39, TMEM38A/TRIC-A, and TMEM38B/TRIC-B). Some of these genes have associations with SkM or heart disease, cancer, bone disease, or other diseases. Among the transcription-related SkM epigenetic features that we identified were: super-enhancers, promoter DNA hypomethylation, lengthening of constitutive low-methylated promoter regions, and SkM-related enhancers for one gene embedded in a neighboring gene (e.g., ASB8-PFKM, LRRC39-DBT, and LRRC14B-PLEKHG4B gene-pairs). In addition, highly or lowly co-expressed long non-coding RNA (lncRNA) genes probably regulate several of these genes. Our findings give insights into tissue-specific epigenetic patterns and functionality of related genes in a gene family and can elucidate normal and disease-related regulation of gene expression in SkM.

Osteoporosis- and obesity-risk interrelationships: an epigenetic analysis of GWAS-derived SNPs at the developmental gene TBX15

A major challenge in translating findings from genome-wide association studies (GWAS) to biological mechanisms is pinpointing functional variants because only a very small percentage of variants associated with a given trait actually impact the trait. We used an extensive epigenetics, transcriptomics, and genetics analysis of the TBX15/WARS2 neighbourhood to prioritize this region's best-candidate causal variants for the genetic risk of osteoporosis (estimated bone density, eBMD) and obesity (waist-hip ratio or waist circumference adjusted for body mass index). TBX15 encodes a transcription factor that is important in bone development and adipose biology. Manual curation of 692 GWAS-derived variants gave eight strong candidates for causal SNPs that modulate TBX15 transcription in subcutaneous adipose tissue (SAT) or osteoblasts, which highly and specifically express this gene. None of these SNPs were prioritized by Bayesian fine-mapping. The eight regulatory causal SNPs were in enhancer or promoter chromatin seen preferentially in SAT or osteoblasts at TBX15 intron-1 or upstream. They overlap strongly predicted, allele-specific transcription factor binding sites. Our analysis suggests that these SNPs act independently of two missense SNPs in TBX15. Remarkably, five of the regulatory SNPs were associated with eBMD and obesity and had the same trait-increasing allele for both. We found that WARS2 obesity-related SNPs can be ascribed to high linkage disequilibrium with TBX15 intron-1 SNPs. Our findings from GWAS index, proxy, and imputed SNPs suggest that a few SNPs, including three in a 0.7-kb cluster, act as causal regulatory variants to fine-tune TBX15 expression and, thereby, affect both obesity and osteoporosis risk.

Data showing atherosclerosis-associated differentially methylated regions are often at enhancers

Atherosclerosis involves phenotypic modulation and transdifferentiation of vascular smooth muscle cells (SMCs). Data are given in tabular or figure format that illustrate genome-wide DNA methylation alterations in atherosclerotic vs. control aorta (athero DMRs). Data based upon publicly available chromatin state profiles are also shown for normal aorta, monocyte, and skeletal muscle tissue-specific DMRs and for aorta-specific chromatin features (enhancer chromatin, promoter chromatin, repressed chromatin, actively transcribed chromatin). Athero hypomethylated and hypermethylated DMRs as well as epigenetic and transcription profiles are described for the following genes: ACTA2, MYH10, MYH11 (SMC-associated genes); SMAD3 (a signaling gene for SMCs and other cell types); CD79B and SH3BP2 (leukocyte-associated genes); and TBX20 and genes in the HOXA, HOXB, HOXC, and HOXD clusters (T-box and homeobox developmental genes). The data reveal strong correlations between athero hypermethylated DMRs and regions of enhancer chromatin in aorta, which are discussed in the linked research article “Atherosclerosis-associated differentially methylated regions can reflect the disease phenotype and are often at enhancers” (M. Lacey et al., 2019).

Tissue-specific epigenetics of atherosclerosis-related ANGPT and ANGPTL genes

Aim:

To understand tissue-specific regulation of angiopoietin/angiopoietin-like (ANGPT/ANGPTL) genes (especially the five genes embedded in introns of host genes) and their association with atherosclerosis.

Methods:

Transcription and epigenomic databases from various normal tissues were examined in the vicinity of ANGPT1, ANGPT2, ANGPTL1, ANGPTL2, ANGPTL3, ANGPTL4 and ANGPTL8.

Results:

We identified tissue-specific enhancer chromatin regions that are likely to regulate transcription of ANGPT/ANGPTL genes and were intragenic, intergenic or host gene-linked. In addition, we found atherosclerosis-linked differentially methylated regions associated with ANGPT2 and with sequences encoding miR-145, a microRNA that targets ANGPT2 mRNA in cancers.

Conclusion:

Our findings implicate enhancers as major contributors to tissue-specific expression of ANGPT/ANGPTL genes, which play critical roles in angiogenesis, atherosclerosis, cancer, and inflammatory and metabolic diseases.

Atherosclerosis-associated differentially methylated regions can reflect the disease phenotype and are often at enhancers

BACKGROUND AND AIMS

Atherosclerosis is a widespread and complicated disease involving phenotypic modulation and transdifferentiation of vascular smooth muscle cells (SMCs), the predominant cells in aorta, as well as changes in endothelial cells and infiltrating monocytes. Alterations in DNA methylation are likely to play central roles in these phenotypic changes, just as they do in normal differentiation and cancer.

METHODS

We examined genome-wide DNA methylation changes in atherosclerotic aorta using more stringent criteria for differentially methylated regions (DMRs) than in previous studies and compared these DMRs to tissue-specific epigenetic features.

RESULTS

We found that disease-linked hypermethylated DMRs account for 85% of the total atherosclerosis-associated DMRs and often overlap aorta-associated enhancer chromatin. These hypermethylated DMRs were associated with functionally different sets of genes compared to atherosclerosis-linked hypomethylated DMRs. The extent and nature of the DMRs could not be explained as direct effects of monocyte/macrophage infiltration. Among the known atherosclerosis- and contractile SMC-related genes that exhibited hypermethylated DMRs at aorta enhancer chromatin were ACTA2 (aorta α2 smooth muscle actin), ELN (elastin), MYOCD (myocardin), C9orf3 (miR-23b and miR-27b host gene), and MYH11 (smooth muscle myosin). Our analyses also suggest a role in atherosclerosis for developmental transcription factor genes having little or no previous association with atherosclerosis, such as NR2F2 (COUP-TFII) and TBX18.

CONCLUSIONS

We provide evidence for atherosclerosis-linked DNA methylation changes in aorta SMCs that might help minimize or reverse the standard contractile character of many of these cells by down-modulating aorta SMC-related enhancers, thereby facilitating pro-atherosclerotic phenotypic changes in many SMCs.

Developmentally linked human DNA hypermethylation is associated with down-modulation, repression, and upregulation of transcription

DNA methylation can affect tissue-specific gene transcription in ways that are difficult to discern from studies focused on genome-wide analyses of differentially methylated regions (DMRs). To elucidate the variety of associations between differentiation-related DNA hypermethylation and transcription, we used available epigenomic and transcriptomic profiles from 38 human cell/tissue types to focus on such relationships in 94 genes linked to hypermethylated DMRs in myoblasts (Mb). For 19 of the genes, promoter-region hypermethylation in Mb (and often a few heterologous cell types) was associated with gene repression but, importantly, DNA hypermethylation was absent in many other repressed samples. In another 24 genes, DNA hypermethylation overlapped cryptic enhancers or super-enhancers and correlated with down-modulated, but not silenced, gene expression. However, such methylation was absent, surprisingly, in both non-expressing samples and highly expressing samples. This suggests that some genes need DMR hypermethylation to help repress cryptic enhancer chromatin only when they are actively transcribed. For another 11 genes, we found an association between intergenic hypermethylated DMRs and positive expression of the gene in Mb. DNA hypermethylation/transcription correlations similar to those of Mb were evident sometimes in diverse tissues, such as aorta and brain. Our findings have implications for the possible involvement of methylated DNA in Duchenne's muscular dystrophy, congenital heart malformations, and cancer. This epigenomic analysis suggests that DNA methylation is not simply the inevitable consequence of changes in gene expression but, instead, is often an active agent for fine-tuning transcription in association with development.

Association of 5-hydroxymethylation and 5-methylation of DNA cytosine with tissue-specific gene expression

Differentially methylated or hydroxymethylated regions (DMRs) in mammalian DNA are often associated with tissue-specific gene expression but the functional relationships are still being unraveled. To elucidate these relationships, we studied 16 human genes containing myogenic DMRs by analyzing profiles of their epigenetics and transcription and quantitatively assaying 5-hydroxymethylcytosine (5hmC) and 5-methylcytosine (5mC) at specific sites in these genes in skeletal muscle (SkM), myoblasts, heart, brain, and diverse other samples. Although most human promoters have little or no methylation regardless of expression, more than half of the genes that we chose to study-owing to their myogenic DMRs-overlapped tissue-specific alternative or cryptic promoters displaying corresponding tissue-specific differences in histone modifications. The 5mC levels in myoblast DMRs were significantly associated with 5hmC levels in SkM at the same site. Hypermethylated myogenic DMRs within CDH15, a muscle- and cerebellum-specific cell adhesion gene, and PITX3, a homeobox gene, were used for transfection in reporter gene constructs. These intragenic DMRs had bidirectional tissue-specific promoter activity that was silenced by in vivo-like methylation. The CDH15 DMR, which was previously associated with an imprinted maternal germline DMR in mice, had especially strong promoter activity in myogenic host cells. These findings are consistent with the controversial hypothesis that intragenic DNA methylation can facilitate transcription and is not just a passive consequence of it. Our results support varied roles for tissue-specific 5mC- or 5hmC-enrichment in suppressing inappropriate gene expression from cryptic or alternative promoters and in increasing the plasticity of gene expression required for development and rapid responses to tissue stress or damage.

DNA Hypomethylation in Intragenic and Intergenic Enhancer Chromatin of Muscle-Specific Genes Usually Correlates with their Expression

Tissue-specific enhancers are critical for gene regulation. In this study, we help elucidate the contribution of muscle-associated differential DNA methylation to the enhancer activity of highly muscle-specific genes. By bioinformatic analysis of 44 muscle-associated genes, we show that preferential gene expression in skeletal muscle (SkM) correlates with SkM-specific intragenic and intergenic enhancer chromatin and overlapping foci of DNA hypomethylation. Some genes, e.g., CASQ1 and FBXO32, displayed broad regions of both SkM- and heart-specific enhancer chromatin but exhibited focal SkM-specific DNA hypomethylation. Half of the genes had SkM-specific super-enhancers. In contrast to simple enhancer/gene-expression correlations, a super-enhancer was associated with the myogenic MYOD1 gene in both SkM and myoblasts even though SkM has < 1 percent as much MYOD1 expression. Local chromatin differences in this super-enhancer probably contribute to the SkM/myoblast differential expression. Transfection assays confirmed the tissue-specificity of the 0.3-kb core enhancer within MYOD1's super-enhancer and demonstrated its repression by methylation of its three CG dinucleotides. Our study suggests that DNA hypomethylation increases enhancer tissue-specificity and that SkM super-enhancers sometimes are poised for physiologically important, rapid up-regulation.

The Aspergillus flavus fluP-associated metabolite promotes sclerotial production

Aspergillus flavus is able to synthesize a variety of polyketide-derived secondary metabolites including the hepatocarcinogen, aflatoxin B1. The fungus reproduces and disseminates predominantly by production of conidia. It also produces hardened mycelial aggregates called sclerotia that are used to cope with unfavourable growth environments. In the present study, we examined the role of A. flavus fluP, the backbone polyketide synthase gene of secondary metabolite gene cluster 41, on fungal development. The A. flavus CA14 fluP deletion mutant (AfΔfluP) grew and accumulated aflatoxin normally but produced a lower amount of sclerotia than the parental strain. This was also true for the Aspergillus parasiticus BN9 fluP deletion mutant (ApΔfluP). The A. flavus fluP gene was positively regulated by developmental regulators of VeA and VelB but not by the global regulator of secondary metabolism, LaeA. Overexpression of fluP in AfΔfluP (OEfluP) elevated its ability to produce sclerotia compared to that of the parental strain. Coculture of OEfluP with CA14, AfΔfluP, ApΔfluP, or an A. flavus pptA deletion mutant incapable of producing functional polyketide synthases also allowed increased sclerotial production of the respective strains at edges where colonies made contact. Acetone extracts of OEfluP but not of AfΔfluP exhibited the same effect in promoting sclerotial production of AfΔfluP. These results suggest that FluP polyketide synthase is involved in the synthesis of a diffusible metabolite that could serve as a signal molecule to regulate sclerotiogenesis.

Comparison of expression of secondary metabolite biosynthesis cluster genes in Aspergillus flavus, A. parasiticus, and A. oryzae

Fifty six secondary metabolite biosynthesis gene clusters are predicted to be in the Aspergillus flavus genome. In spite of this, the biosyntheses of only seven metabolites, including the aflatoxins, kojic acid, cyclopiazonic acid and aflatrem, have been assigned to a particular gene cluster. We used RNA-seq to compare expression of secondary metabolite genes in gene clusters for the closely related fungi A. parasiticus, A. oryzae, and A. flavus S and L sclerotial morphotypes. The data help to refine the identification of probable functional gene clusters within these species. Our results suggest that A. flavus, a prevalent contaminant of maize, cottonseed, peanuts and tree nuts, is capable of producing metabolites which, besides aflatoxin, could be an underappreciated contributor to its toxicity.

Non-aflatoxigenic Aspergillus flavus to prevent aflatoxin contamination in crops: advantages and limitations

Aspergillus flavus is a diverse assemblage of strains that include aflatoxin-producing and non-toxigenic strains with cosmopolitan distribution. The most promising strategy currently being used to reduce preharvest contamination of crops with aflatoxin is to introduce non-aflatoxin (biocontrol) A. flavus into the crop environment. Whether or not introduction of biocontrol strains into agricultural fields is enough to reduce aflatoxin contamination to levels required for acceptance of the contaminated food as fit for consumption is still unknown. There is no question that biocontrol strains are able to reduce the size of the populations of aflatoxin-producing strains but the available data suggests that at most only a four- to five-fold reduction in aflatoxin contamination is achieved. There are many challenges facing this strategy that are both short term and long term. First, the population biology of A. flavus is not well understood due in part to A. flavus's diversity, its ability to form heterokaryotic reproductive forms, and its unknown ability to survive for prolonged periods after application. Second, biocontrol strains must be selected that are suitable for the environment, the type of crop, and the soil into which they will be introduced. Third, there is a need to guard against inadvertent introduction of A. flavus strains that could impose an additional burden on food safety and food quality, and fourth, with global warming and resultant changes in the soil nutrients and concomitant microbiome populations, the biocontrol strategy must be sufficiently flexible to adapt to such changes. Understanding genetic variation within strains of A. flavus is important for developing a robust biocontrol strategy and it is unlikely that a "one size fits all" strategy will work for preharvest aflatoxin reduction.

Identification of novel metabolites from Aspergillus flavus by high resolution and multiple stage mass spectrometry

The filamentous fungus Aspergillus flavus is one of the most important species in the Aspergillus genus and is distributed worldwide as a prevalent aflatoxin-producing food and feed contaminant. A. flavus contains more than 55 gene clusters that are predicted to encode proteins involved in secondary metabolite production. One of these, cluster 27, contains a polyketide synthase (pks27) gene that encodes a protein that is highly homologous to the aflatoxin cluster PKS. Comparative metabolomics, using ultra-high performance liquid chromatography (UHPLC) coupled to high resolution Orbitrap mass spectrometry (MS) was used to detect metabolites differentially expressed in the A. flavus wild-type and ∆pks27 mutant strains. Metabolite profiling was aided by a statistical differential analysis of MS data using SIEVE software. This differential analysis combined with accurate mass data from the Orbitrap and ion trap multiple stage MS allowed four metabolites to be identified that were produced only by the wild-type culture. These included asparasone A (358 Da), an anthraquinone pigment, and related anthraquinones with masses of 316, 340 and 374 Da. These latter three compounds had similar fragmentation patterns to that of asparasone A. The 316 Da anthraquinone is particularly interesting because it is most likely formed by incorporation of seven malonyl-CoA units rather than the eight units required for the formation of asparasone A. The 340 and 374 Da metabolites are the dehydration and an oxy-derivative of asparasone A, respectively. Asparasone A was also identified in extracts from several other Aspergillus species.

Aspergillus flavus VelB acts distinctly from VeA in conidiation and may coordinate with FluG to modulate sclerotial production

The proteins VeA, VelB and LaeA of Aspergillus nidulans form a heterotrimeric complex (the velvet complex) in the dark to coordinate sexual development and production of some secondary metabolites. VeA and VelB of A. nidulans and Aspergillus fumigatus also are repressors of conidiation, but VeA of Aspergillus flavus in studied strains acts positively on conidiation. In the present study, we show via yeast-two hybrid assays that interactions among A. flavus VeA, VelB, and LaeA are conserved as in the A. nidulans velvet complex. We found that FluG, which is required for conidiophore formation in A. nidulans but whose deletion in A. flavus delays onset of conidiation, was probably an interacting partner of VelB. Deletion of velB in A. flavus CA14 severely impaired conidiation in the dark although to a lesser extent than deletion of veA. In both mutants fluG deletion resulted in further decreased conidiation even in the light. Deletion of fluG in the ΔlaeA strain, however, did not affect conidiation. All mutant types were unable to produce aflatoxin and sclerotia. Cross-complementation of the ΔvelB strain with gpdA::veA restored conidiation but not aflatoxin production although aflR, the aflatoxin pathway regulatory gene, was expressed at a normal level. Cross-complementation of the ΔveA strain with gpdA::velB failed to restore conidiation and aflatoxin production. The ΔvelB strain complemented with or a wild type transformed by gpdA::velB had elevated sclerotial production as the ΔfluG strain. Concerted interactions of A. flavus VeA and VelB with LaeA are critical for conidiation and aflatoxin biosynthesis. VelB may have a dual role and likely coordinates with FluG to modulate sclerotial production.

Association with AflR in endosomes reveals new functions for AflJ in aflatoxin biosynthesis

Aflatoxins are the most potent naturally occurring carcinogens of fungal origin. Biosynthesis of aflatoxin involves the coordinated expression of more than 25 genes. The function of one gene in the aflatoxin gene cluster, aflJ, is not entirely understood but, because previous studies demonstrated a physical interaction between the Zn₂Cys₆ transcription factor AflR and AflJ, AflJ was proposed to act as a transcriptional co-activator. Image analysis revealed that, in the absence of aflJ in A. parasiticus, endosomes cluster within cells and near septa. AflJ fused to yellow fluorescent protein complemented the mutation in A. parasiticus ΔaflJ and localized mainly in endosomes. We found that AflJ co-localizes with AflR both in endosomes and in nuclei. Chromatin immunoprecipitation did not detect AflJ binding at known AflR DNA recognition sites suggesting that AflJ either does not bind to these sites or binds to them transiently. Based on these data, we hypothesize that AflJ assists in AflR transport to or from the nucleus, thus controlling the availability of AflR for transcriptional activation of aflatoxin biosynthesis cluster genes. AflJ may also assist in directing endosomes to the cytoplasmic membrane for aflatoxin export.

Fragmentation of an aflatoxin-like gene cluster in a forest pathogen

Plant pathogens use a complex arsenal of weapons, such as toxic secondary metabolites, to invade and destroy their hosts. Knowledge of how secondary metabolite pathways evolved is central to understanding the evolution of host specificity. The secondary metabolite dothistromin is structurally similar to aflatoxins and is produced by the fungal pine pathogen Dothistroma septosporum. Our study focused on dothistromin genes, which are widely dispersed across one chromosome, to determine whether this unusual distributed arrangement evolved from an ancestral cluster. We combined comparative genomics and population genetics approaches to elucidate the origins of the dispersed arrangement of dothistromin genes over a broad evolutionary time-scale at the phylum, class and species levels. Orthologs of dothistromin genes were found in two major classes of fungi. Their organization is consistent with clustering of core pathway genes in a common ancestor, but with intermediate cluster fragmentation states in the Dothideomycetes fungi. Recombination hotspots in a D. septosporum population matched sites of gene acquisition and cluster fragmentation at higher evolutionary levels. The results suggest that fragmentation of a larger ancestral cluster gave rise to the arrangement seen in D. septosporum. We propose that cluster fragmentation may facilitate metabolic retooling and subsequent host adaptation of plant pathogens.

Functional and phylogenetic analysis of the Aspergillus ochraceoroseus aflQ (ordA) gene ortholog

Within the Aspergillus parasiticus and A. flavus aflatoxin (AF) biosynthetic gene cluster the aflQ (ordA) and aflP (omtA) genes encode respectively an oxidoreductase and methyltransferase. These genes are required for the final steps in the conversion of sterigmatocystin (ST) to aflatoxin B(1) (AFB(1)). Aspergillus nidulans harbors a gene cluster that produces ST, as the aflQ and aflP orthologs are either non-functional or absent in the genome. Aspergillus ochraceoroseus produces both AF and ST, and it harbors an AF/ST biosynthetic gene cluster that is organized much like the A. nidulans ST cluster. The A. ochraceoroseus cluster also does not contain aflQ or aflP orthologs. However the ability of A. ochraceoroseus to produce AF would indicate that functional aflQ and aflP orthologs are present within the genome. Utilizing degenerate primers based on conserved regions of the A. flavus aflQ gene and an A. nidulans gene demonstrating the highest level of homology to aflQ, a putative aflQ ortholog was PCR amplified from A. ochraceoroseus genomic DNA. The A. ochraceoroseus aflQ ortholog demonstrated 57% amino acid identity to A. flavus AflQ. Transformation of an O-methylsterigmatocystin (OMST)-accumulating A. parasiticus aflQ mutant with the putative A. ochraceoroseus aflQ gene restored AF production. Although the aflQ gene does not reside in the AF/ST cluster it appears to be regulated in a manner similar to other A. ochraceoroseus AF/ST cluster genes. Phylogenetic analysis of AflQ and AflQ-like proteins from a number of ST- and AF-producing Aspergilli indicates that A. ochraceoroseus might be ancestral to A. nidulans and A. flavus.

Effects of laeA deletion on Aspergillus flavus conidial development and hydrophobicity may contribute to loss of aflatoxin production

LaeA of Aspergillus nidulans is a putative methyltransferase and a component of the velvet complex; it is thought to mainly affect expression of genes required for the production of secondary metabolites. We found that although Aspergillus flavus CA14 laeA deletion mutants showed no aflatoxin production, expression of some of the early genes involved in aflatoxin formation, but not the later genes, could still be detected. The mutants grown in minimal medium supplemented with simple sugars and on some complex media exhibited altered conidial development. On potato dextrose agar (PDA) medium the deletion mutants showed reduced conidial chain elongation, increased production of conidiophores, and decreased colony hydrophobicity when compared to the parental strain. The loss of hydrophobicity and the other developmental changes in the laeA deletion mutants could affect the ability of the fungus to produce aflatoxins.

Expression profiling of non-aflatoxigenic Aspergillus parasiticus mutants obtained by 5-azacytosine treatment or serial mycelial transfer

Aflatoxins are carcinogenic secondary metabolites produced by the fungi Aspergillus flavus and Aspergillus parasiticus. Previous studies found that repeated serial mycelial transfer or treatment of A. parasiticus with 5-azacytidine produced colonies with a fluffy phenotype and inability to produce aflatoxins. To understand how these treatments affect expression of genes involved in aflatoxin production and development, we carried out expressed sequence tag (EST)-based microarray assays to identify genes in treated clones that are differentially expressed compared to the wild-type. Expression of 183 genes was significantly dysregulated. Of these, 38 had at least two-fold or lower expression compared to the untreated control and only two had two-fold or higher expression. The most frequent change was downregulation of genes predicted to encode membrane-bound proteins. Based on this result we hypothesize that the treatments cause changes in the structure of cellular and organelle membranes that prevent normal development and aflatoxin biosynthesis.

Conversion of 11-hydroxy-O-methylsterigmatocystin to aflatoxin G1 in Aspergillus parasiticus

In aflatoxin biosynthesis, aflatoxins G(1) (AFG(1)) and B(1) (AFB(1)) are independently produced from a common precursor, O-methylsterigmatocystin (OMST). Recently, 11-hydroxy-O-methylsterigmatocystin (HOMST) was suggested to be a later precursor involved in the conversion of OMST to AFB(1), and conversion of HOMST to AFB(1) was catalyzed by OrdA enzyme. However, the involvement of HOMST in AFG(1) formation has not been determined. In this work, HOMST was prepared by incubating OrdA-expressing yeast with OMST. Feeding Aspergillus parasiticus with HOMST allowed production of AFG(1) as well as AFB(1). In cell-free systems, HOMST was converted to AFG(1) when the microsomal fraction, the cytosolic fraction from A. parasiticus, and yeast expressing A. parasiticus OrdA were added. These results demonstrated (1) HOMST is produced from OMST by OrdA, (2) HOMST is a precursor of AFG(1) as well as AFB(1), and (3) three enzymes, OrdA, CypA, and NadA, and possibly other unknown enzymes are involved in conversion of HOMST to AFG(1).

Absence of the aflatoxin biosynthesis gene, norA, allows accumulation of deoxyaflatoxin B1 in Aspergillus flavus cultures

Biosynthesis of the highly toxic and carcinogenic aflatoxins in select Aspergillus species from the common intermediate O-methylsterigmatocystin has been postulated to require only the cytochrome P450 monooxygenase, OrdA (AflQ). We now provide evidence that the aryl alcohol dehydrogenase NorA (AflE) encoded by the aflatoxin biosynthetic gene cluster in Aspergillus flavus affects the accumulation of aflatoxins in the final steps of aflatoxin biosynthesis. Mutants with inactive norA produced reduced quantities of aflatoxin B(1) (AFB(1)), but elevated quantities of a new metabolite, deoxyAFB(1). To explain this result, we suggest that, in the absence of NorA, the AFB(1) reduction product, aflatoxicol, is produced and is readily dehydrated to deoxyAFB(1) in the acidic medium, enabling us to observe this otherwise minor toxin produced in wild-type A. flavus.

Cyclopiazonic acid biosynthesis of Aspergillus flavus and Aspergillus oryzae

Cyclopiazonic acid (CPA) is an indole-tetramic acid neurotoxin produced by some of the same strains of A. flavus that produce aflatoxins and by some Aspergillus oryzae strains. Despite its discovery 40 years ago, few reviews of its toxicity and biosynthesis have been reported. This review examines what is currently known about the toxicity of CPA to animals and humans, both by itself or in combination with other mycotoxins. The review also discusses CPA biosynthesis and the genetic diversity of CPA production in A. flavus/oryzae populations.

Predicted roles of the uncharacterized clustered genes in aflatoxin biosynthesis

Biosynthesis of the toxic and carcinogenic aflatoxins (AFs) requires the activity of more than 27 enzymes. The roles in biosynthesis of newly described enzymes are discussed in this review. We suggest that HypC catalyzes the oxidation of norsolorinic acid anthrone; AvfA (AflI), the ring-closure step in formation of hydroxyversicolorone; HypB, the second oxidation step in conversion of O-methylsterigmatocystin to AF; and HypE and NorA (AflE), the final two steps in AFB(1) formation. HypD, an integral membrane protein, affects fungal development and lowers AF production while AflJ (AflS), has a partial methyltransferase domain that may be important in its function as a transcriptional co-activator.

Characterization of the Aspergillus ochraceoroseus aflatoxin/sterigmatocystin biosynthetic gene cluster

Production of carcinogenic aflatoxins has been reported from members of Aspergillus section Flavi, Aspergillus section Nidulantes and a newly proposed Aspergillus section Ochraceorosei that consists of Aspergillus ochraceoroseus and A. rambellii. Unlike members of section Flavi, A. ochraceoroseus and A. rambellii have been shown to accumulate both aflatoxin (AF) and the aflatoxin precursor sterigmatocystin (ST). Alhough morphologically distinct from A. nidulans, molecular characterization of A. ochraceoroseus AF/ST genes and physiological characteristics of AF/ST production indicated that A. ochraceoroseus is more closely related to A. nidulans than to A. flavus. Knowing that the A. nidulans ST gene cluster is organized differently from the A. flavus AF gene cluster, we determined the genetic organization of the AF/ST biosynthetic cluster in A. ochraceoroseus. Sequencing of overlapping lambda clones and genomic PCR fragments obtained by gene-walking techniques demonstrated that the A. ochraceoroseus AF/ST gene cluster is organized much like the A. nidulans ST gene cluster except that the region from aflN to aflW is located directly upstream of aflC and in reverse orientation such that aflW represents the distal end and aflY the proximal end of the cluster. The A. ochraceoroseus cluster genes demonstrated 62-76% nucleotide identity to their A. nidulans ST cluster gene homologs. Transformation of an A. nidulans aflR mutant with the A. ochraceoroseus aflR restored ST production in A. nidulans transformants. PCR amplification of A. rambellii genomic DNA demonstrated that the AF/ST gene cluster is organized in the same manner as that of A. ochraceoroseus.

Are the genes nadA and norB involved in formation of aflatoxin G(1)?

Aflatoxins, the most toxic and carcinogenic family of fungal secondary metabolites, are frequent contaminants of foods intended for human consumption. Previous studies showed that formation of G-group aflatoxins (AFs) from O-methylsterigmatocystin (OMST) by certain Aspergillus species involves oxidation by the cytochrome P450 monooxygenases, OrdA (AflQ) and CypA (AflU). However, some of the steps in the conversion have not yet been fully defined. Extracts of Aspergillus parasiticus disruption mutants of the OYE-FMN binding domain reductase-encoding gene nadA (aflY) contained a 386 Da AFG(1) precursor. A compound with this mass was predicted as the product of sequential OrdA and CypA oxidation of OMST. Increased amounts of a 362 Da alcohol, the presumptive product of NadA reduction, accumulate in extracts of fungi with disrupted aryl alcohol dehydrogenase-encoding gene norB. These results show that biosynthesis of AFG(1) involves NadA reduction and NorB oxidation.

Use of pyrosequencing to quantify incidence of a specific Aspergillus flavus strain within complex fungal communities associated with commercial cotton crops

Atoxigenic strains of Aspergillus flavus have been used as aflatoxin management tools on over 50,000 hectares of commercial crops since 2000. To assess treatment efficacy, atoxigenic strain incidence is routinely monitored by vegetative compatibility analyses (VCA) that require culturing, generation of auxotrophs, and complementation with tester mutants. Two pyrosequencing assays (PA) that require no culturing were developed for monitoring incidences of atoxigenic strains on ginned cottonseed. The assays, which quantify frequencies of characteristic single nucleotide polymorphisms (SNPs) in the aflR and pksA genes, were validated against standard VCA on cottonseed collected from commercial gins in South Texas, Arizona, and Southern California where the atoxigenic strain AF36 is used to manage aflatoxin contamination. Cottonseed washings were subjected to both VCA and PA. PA was performed directly on DNA isolated from particulates pelleted from the wash water by centrifugation. Addition of CaCl(2) and diatomaceous earth prior to pelleting increased the amount of DNA isolated. Accuracy and reproducibility of the PA were contrasted with those for the VCA that has been used for over a decade. Correlation coefficients between VCA and PA indicated good correspondence between the results from the two assays (r = 0.91 for aflR assay and r = 0.80 for pksA assay). PAs were highly variable for samples with low incidences of A. flavus due to variability in the initial polymerase chain reaction step. This held for both DNA isolated from cottonseed washes and for mixtures of purified DNA. For samples yielding low quantities of A. flavus DNA, averaging of results from 4 to 5 replicates was required to achieve acceptable correlations with VCA. Pyrosequencing has the potential to become a powerful tool for monitoring atoxigenic strains within complex A. flavus communities without limitations imposed by traditional culturing methods.

Starter unit specificity directs genome mining of polyketide synthase pathways in fungi

Search of the protein database with the aflatoxin pathway polyketide synthase (PKS) revealed putative PKSs in the pathogenic fungi Coccidioides immitis and Coccidioides posadasii that could require partnerships with a pair of fatty acid synthase (FAS) subunits for the biosynthesis of fatty acid-polyketide hybrid metabolites. A starter unit:acyl-carrier protein transacylase (SAT) domain was discovered in the nonreducing PKS. This domain is thought to accept the fatty acid product from the FAS to initiate polyketide synthesis. We expressed the C. immitis SAT domain in Escherichia coli and showed that this domain, unlike that from the aflatoxin pathway PKS, transferred octanoyl-CoA four times faster than hexanoyl-CoA. The SAT domain also formed a covalent octanoyl intermediate and transferred this group to a free-standing ACP domain. Our results suggest that C. immitis/posadasii, both human fungal pathogens, contain a FAS/PKS cluster with functional similarity to the aflatoxin cluster found in Aspergillus species. Dissection of the PKS and determination of in vitro SAT domain specificity provides a tool to uncover the growing number of similar sequenced pathways in fungi, and to guide elucidation of the fatty acid-polyketide hybrid metabolites that they produce.

Aflatoxin-producing Aspergillus species from Thailand

Aflatoxin-producing Aspergillus species were isolated from soil samples from ten different regions within Thailand. Aspergillus flavus was present in all of the soil samples. Unlike previous studies, we found no A. parasiticus or A. flavus capable of both B- and G-type aflatoxin production in any of the samples. A. pseudotamarii, which had not been previously reported from Thailand, was found in four soil samples. In two of the samples A. nomius was determined to be the most abundant aflatoxin-producing species. Based on sequence alignments for three DNA regions (Taka-amylase A (taa), the rRNA internal transcribed spacer (ITS), and the intergenic region for the aflatoxin biosynthesis genes aflJ and aflR) the A. nomius isolates separated into three well-supported clades. Isolates from one of the A. nomius clades had morphological properties similar to those found for S-type isolates capable of B and G aflatoxin production and could easily be mistaken for these isolates. Our results suggest that such unusual A. nomius isolates could be a previously unrecognized agent for aflatoxin contamination in Thailand.

Analysis of single nucleotide polymorphisms in three genes shows evidence for genetic isolation of certain Aspergillus flavus vegetative compatibility groups

Genetic exchange by asexual filamentous fungi is presumed to be limited to isolates in the same vegetative compatibility group (VCG). To evaluate genetic isolation of Aspergillus flavus due to vegetative incompatibility, three gene regions were chosen that contained closely spaced nucleotides that were polymorphic among some of the six VCGs examined. A member of each VCG was collected from five regions across the southern United States. Isolates belonging to the same VCG had similar sets of single nucleotide polymorphisms regardless of isolate origin. The six VCGs formed four genetically distinct groups. Although recombination between certain pairs of VCGs could not be excluded, none was found for YV36, the VCG that includes the atoxigenic A. flavus isolate currently used to mitigate aflatoxin contamination in cotton in Arizona.

Aflatoxigenicity in Aspergillus: molecular genetics, phylogenetic relationships and evolutionary implications

Aflatoxins (AFs) are toxic and carcinogenic secondary metabolites produced by isolates of Aspergillus section Flavi as well as a number of Aspergillus isolates that are classified outside of section Flavi. Characterization of the AF and sterigmatocystin (ST) gene clusters and analysis of factors governing regulation of their biosynthesis has resulted in these two mycotoxins being the most extensively studied of fungal secondary metabolites. This wealth of information has allowed the determination of the molecular basis for non-production of AF in natural isolates of A. flavus and domesticated strains of A. oryzae. This review provides an overview of the molecular analysis of the AF and ST gene clusters as well as new information on an AF gene cluster identified in the non-section Flavi isolate, Aspergillus ochraceoroseus. Additionally, molecular phylogenetic analysis using AF biosynthetic gene sequences as well as ribosomal DNA internal transcribed spacer (ITS) sequences between various section Flavi and non-section Flavi species has enabled determination of the probable evolutionary history of the AF and ST gene clusters. A model for the evolution of the AF and ST gene clusters as well as possible biological roles for AF are discussed.

The aflatoxin biosynthesis cluster gene, aflX, encodes an oxidoreductase involved in conversion of versicolorin A to demethylsterigmatocystin

Biosynthesis of the toxic and carcinogenic aflatoxins by the fungus Aspergillus flavus is a complicated process involving more that 27 enzymes and regulatory factors encoded by a clustered group of genes. Previous studies found that three enzymes, encoded by verA, ver-1, and aflY, are required for conversion of versicolorin A (VA), to demethylsterigmatocystin. We now show that a fourth enzyme, encoded by the previously uncharacterized gene, aflX (ordB), is also required for this conversion. A homolog of this gene, stcQ, is present in the A. nidulans sterigmatocystin (ST) biosynthesis cluster. Disruption of aflX in Aspergillus flavus gave transformants that accumulated approximately 4-fold more VA and fourfold less aflatoxin than the untransformed strain. Southern and Northern blot analyses confirmed that aflX was the only gene disrupted in these transformants. Feeding ST or O-methylsterigmatocystin, but not VA or earlier precursor metabolites, restored normal levels of AF production. The protein encoded by aflX is predicted to have domains typical of an NADH-dependent oxidoreductase. It has 27% amino acid identity to a protein encoded by the aflatoxin cluster gene, aflO (avfA). Some of domains in the protein are similar to those of epoxide hydrolases.

Aflatoxin biosynthesis gene clusters and flanking regions

AIMS

To compare the biosynthetic gene cluster sequences of the main aflatoxin (AF)-producing Aspergillus species.

METHODS AND RESULTS

Sequencing was on fosmid clones selected by homology to Aspergillus parasiticus sequence. Alignments revealed that gene order is conserved among AF gene clusters of Aspergillus nomius, A. parasiticus, two sclerotial morphotypes of Aspergillus flavus, and an unnamed Aspergillus sp. Phylogenetic relationships were established using the maximum likelihood method implemented in PAUP. Based on the Eurotiomycete/Sordariomycete divergence time, the A. flavus-type cluster has been maintained for at least 25 million years. Such conservation of the genes and gene order reflects strong selective constraints on rearrangement. Phylogenetic comparison of individual genes in the cluster indicated that ver-1, which has homology to a melanin biosynthesis gene, experienced selective forces distinct from the other pathway genes. Sequences upstream of the polyketide synthase-encoding gene vary among the species, but a four-gene sugar utilization cluster at the distal end is conserved, indicating a functional relationship between the two adjacent clusters.

CONCLUSIONS

The high conservation of cluster components needed for AF production suggests there is an adaptive value for AFs in character-shaping niches important to those taxa.

SIGNIFICANCE AND IMPACT OF THE STUDY

This is the first comparison of the complete nucleotide sequences of gene clusters harbouring the AF biosynthesis genes of the main AF-producing species. Such a comparison will aid in understanding how AF biosynthesis is regulated in experimental and natural environments.

An aflatoxin biosynthesis cluster gene encodes a novel oxidase required for conversion of versicolorin a to sterigmatocystin

Disruption of the aflatoxin biosynthesis cluster gene aflY (hypA) gave Aspergillus parasiticus transformants that accumulated versicolorin A. This gene is predicted to encode the Baeyer-Villiger oxidase necessary for formation of the xanthone ring of the aflatoxin precursor demethylsterigmatocystin.

Cladal relatedness among Aspergillus oryzae isolates and Aspergillus flavus S and L morphotype isolates

Aspergillus flavus is the main etiological agent for aflatoxin contamination of crops. Its close relative, A. oryzae, does not produce aflatoxins and has been widely used to produce fermented foods. We compared the phylogeny of A. oryzae isolates and L- and S-type sclerotial isolates of A. flavus using single nucleotide polymorphisms in the omtA gene in the aflatoxin biosynthesis gene cluster and deletions in and distal to the norB-cypA intergenic region as phylogenetic signals. Aflatoxin-producing ability and sclerotial size also were weighted in the analysis. Like A. flavus, the A. oryzae isolates form a polyphyletic assemblage. A. oryzae isolates in one clade strikingly resemble an A. flavus subgroup of atoxigenic L-type isolates. All toxigenic S-type isolates closely resemble another subgroup of atoxigenic L-type isolates. Because atoxigenic S-type isolates are extremely rare, we hypothesize that loss of aflatoxin production in S-type isolates may occur concomitantly with a change to L-type sclerotia. All toxigenic L-type isolates, unlike A. oryzae, have a 1.0 kb deletion in the norB-cypA region. Although A. oryzae isolates, like S-type, have a 1.5 kb deletion in the norB-cypA region, none were cladally related to S-type A. flavus isolates. Our results show that A. flavus populations are genetically diverse. A. oryzae isolates may descend from certain atoxigenic L-type A. flavus isolates.

Divergent regulation of aflatoxin production at acidic pH by two Aspergillus strains

Production of aflatoxins (AF) by Aspergillus flavus and A. parasiticus is known to occur only at acidic pH. Although typical A. flavus isolates produced more AF as the external pH became increasingly acidic, an atypical strain from West Africa produced less. The lower AF production was not well correlated with decreases in expression of the aflatoxin pathway regulatory gene, aflR, or of two other biosynthesis genes.