Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Transcriptome analysis reveals key genes associated with root-lesion nematode Pratylenchus thornei resistance in chickpea.

The root-lesion nematode, Pratylenchus thornei, is one of the major plant-parasitic nematode species causing significant yield losses in chickpea (Cicer arietinum). In order to identify the underlying mechanisms of resistance to P. thornei, the transcriptomes of control and inoculated roots of three chickpea genotypes viz. D05253 > F3TMWR2AB001 (resistant advanced breeding line), PBA HatTrick (moderately resistant cultivar), and Kyabra (susceptible cultivar) were studied at 20 and 50 days post inoculation using the RNA-seq approach. On analyzing the 633.3 million reads generated, 962 differentially expressed genes (DEGs) were identified. Comparative analysis revealed that the majority of DEGs upregulated in the resistant genotype were downregulated in the moderately resistant and susceptible genotypes. Transcription factor families WRKY and bZIP were uniquely expressed in the resistant genotype. The genes Cysteine-rich receptor-like protein kinase 10, Protein lifeguard-like, Protein detoxification, Bidirectional sugar transporter Sugars Will Eventually be Exported Transporters1 (SWEET1), and Subtilisin-like protease were found to play cross-functional roles in the resistant chickpea genotype against P. thornei. The identified candidate genes for resistance to P. thornei in chickpea can be explored further to develop markers and accelerate the introgression of P. thornei resistance into elite chickpea cultivars.

Genome-Wide SNP Discovery and Mapping QTLs for Seed Iron and Zinc Concentrations in Chickpea (Cicer arietinum L.).

Biofortification through plant breeding is a cost-effective and sustainable approach towards addressing micronutrient malnutrition prevailing across the globe. Screening cultivars for micronutrient content and identification of quantitative trait loci (QTLs)/genes and markers help in the development of biofortified varieties in chickpea (Cicer arietinum L.). With the aim of identifying the genomic regions controlling seed Fe and Zn concentrations, the F2:3 population derived from a cross between MNK-1 and Annigeri 1 was genotyped using genotyping by sequencing approach and evaluated for Fe and Zn concentration. An intraspecific genetic linkage map comprising 839 single nucleotide polymorphisms (SNPs) spanning a total distance of 1,088.04 cM with an average marker density of 1.30 cM was constructed. By integrating the linkage map data with the phenotypic data of the F2:3 population, a total of 11 QTLs were detected for seed Fe concentration on CaLG03, CaLG04, and CaLG05, with phenotypic variation explained ranging from 7.2% (CaqFe3.4) to 13.4% (CaqFe4.2). For seed Zn concentration, eight QTLs were identified on CaLG04, CaLG05, and CaLG08. The QTLs individually explained phenotypic variations ranging between 5.7% (CaqZn8.1) and 13.7% (CaqZn4.3). Three QTLs for seed Fe and Zn concentrations (CaqFe4.4, CaqFe4.5, and CaqZn4.1) were colocated in the "QTL-hotspot" region on CaLG04 that harbors several drought tolerance-related QTLs. We identified genes in the QTL regions that encode iron-sulfur metabolism and zinc-dependent alcohol dehydrogenase activity on CaLG03, iron ion binding oxidoreductase on CaLG04, and zinc-induced facilitator-like protein and ZIP zinc/iron transport family protein on CaLG05. These genomic regions and the associated markers can be used in marker-assisted selection to increase seed Fe and Zn concentrations in agronomically superior chickpea varieties.

Carbohydrate, glutathione, and polyamine metabolism are central to Aspergillus flavus oxidative stress responses over time.

BACKGROUND: The primary and secondary metabolites of fungi are critical for adaptation to environmental stresses, host pathogenicity, competition with other microbes, and reproductive fitness. Drought-derived reactive oxygen species (ROS) have been shown to stimulate aflatoxin production and regulate in Aspergillus flavus, and may function in signaling with host plants. Here, we have performed global, untargeted metabolomics to better understand the role of aflatoxin production in oxidative stress responses, and also explore isolate-specific oxidative stress responses over time. RESULTS: Two field isolates of A. flavus, AF13 and NRRL3357, possessing high and moderate aflatoxin production, respectively, were cultured in medium with and without supplementation with 15 mM H2O2, and mycelia were collected following 4 and 7 days in culture for global metabolomics. Overall, 389 compounds were described in the analysis which encompassed 9 biological super-pathways and 47 sub-pathways. These metabolites were examined for differential accumulation. Significant differences were observed in both isolates in response to oxidative stress and when comparing sampling time points. CONCLUSIONS: The moderately high aflatoxin-producing isolate, NRRL3357, showed extensive stimulation of antioxidant mechanisms and pathways including polyamines metabolism, glutathione metabolism, TCA cycle, and lipid metabolism while the highly aflatoxigenic isolate, AF13, showed a less vigorous response to stress. Carbohydrate pathway levels also imply that carbohydrate repression and starvation may influence metabolite accumulation at the later timepoint. Higher conidial oxidative stress tolerance and antioxidant capacity in AF13 compared to NRRL3357, inferred from their metabolomic profiles and growth curves over time, may be connected to aflatoxin production capability and aflatoxin-related antioxidant accumulation. The coincidence of several of the detected metabolites in H2O2-stressed A. flavus and drought-stressed hosts also suggests their potential role in the interaction between these organisms and their use as markers/targets to enhance host resistance through biomarker selection or genetic engineering.

The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication.

High oil and protein content make tetraploid peanut a leading oil and food legume. Here we report a high-quality peanut genome sequence, comprising 2.54 Gb with 20 pseudomolecules and 83,709 protein-coding gene models. We characterize gene functional groups implicated in seed size evolution, seed oil content, disease resistance and symbiotic nitrogen fixation. The peanut B subgenome has more genes and general expression dominance, temporally associated with long-terminal-repeat expansion in the A subgenome that also raises questions about the A-genome progenitor. The polyploid genome provided insights into the evolution of Arachis hypogaea and other legume chromosomes. Resequencing of 52 accessions suggests that independent domestications formed peanut ecotypes. Whereas 0.42-0.47 million years ago (Ma) polyploidy constrained genetic variation, the peanut genome sequence aids mapping and candidate-gene discovery for traits such as seed size and color, foliar disease resistance and others, also providing a cornerstone for functional genomics and peanut improvement.

Responses of Aspergillus flavus to Oxidative Stress Are Related to Fungal Development Regulator, Antioxidant Enzyme, and Secondary Metabolite Biosynthetic Gene Expression.

The infection of maize and peanut with Aspergillus flavus and subsequent contamination with aflatoxin pose a threat to global food safety and human health, and is exacerbated by drought stress. Drought stress-responding compounds such as reactive oxygen species (ROS) are associated with fungal stress responsive signaling and secondary metabolite production, and can stimulate the production of aflatoxin by A. flavus in vitro. These secondary metabolites have been shown to possess diverse functions in soil-borne fungi including antibiosis, competitive inhibition of other microbes, and abiotic stress alleviation. Previously, we observed that isolates of A. flavus showed differences in oxidative stress tolerance which correlated with their aflatoxin production capabilities. In order to better understand these isolate-specific oxidative stress responses, we examined the transcriptional responses of field isolates of A. flavus with varying levels of aflatoxin production (NRRL3357, AF13, and Tox4) to H2O2-induced oxidative stress using an RNA sequencing approach. These isolates were cultured in an aflatoxin-production conducive medium amended with various levels of H2O2. Whole transcriptomes were sequenced using an Illumina HiSeq platform with an average of 40.43 million filtered paired-end reads generated for each sample. The obtained transcriptomes were then used for differential expression, gene ontology, pathway, and co-expression analyses. Isolates which produced higher levels of aflatoxin tended to exhibit fewer differentially expressed genes than isolates with lower levels of production. Genes found to be differentially expressed in response to increasing oxidative stress included antioxidant enzymes, primary metabolism components, antibiosis-related genes, and secondary metabolite biosynthetic components specifically for aflatoxin, aflatrem, and kojic acid. The expression of fungal development-related genes including aminobenzoate degradation genes and conidiation regulators were found to be regulated in response to increasing stress. Aflatoxin biosynthetic genes and antioxidant enzyme genes were also found to be co-expressed and highly correlated with fungal biomass under stress. This suggests that these secondary metabolites may be produced as part of coordinated oxidative stress responses in A. flavus along with antioxidant enzyme gene expression and developmental regulation.

The binding modes and binding affinities of artemisinin derivatives with Plasmodium falciparum Ca2+-ATPase (PfATP6).

Noncompetitive inhibitors of sarco/endoplasmic reticulum Ca(2+)-ATPase (SERCA) orthologue (PfATP6) of P. falciparum have important therapeutic value in the treatment of malaria. Artemisinin and its analogues are one such class of inhibitors which bind to a hydrophobic pocket located in the transmembrane region of PfATP6 near the biomembrane surface and interfere with calcium transport. The 3D structure of PfATP6 was modeled by homology modeling. A library consisting of 150 artemisinin analogues has been designed. Their molecular interactions and binding affinities with modeled PfATP6 protein have been studied using the docking, molecular mechanics based on generalized Born/surface area (MM-GBSA) solvation model and multi-ligand bimolecular association with energetics (eMBrAcE). Structure activity relationship models were developed between the antimalarial activity (log RA) and molecular descriptors like docking score and binding free energy. For both the cases the r(2) was in the range of 0.538-0.688 indicating good data fit and r(2)(cv) was in the range of 0.525-0.679 indicating that the predictive capabilities of the models were acceptable. Besides, a scheme similar to linear response was used to develop free energy of binding (FEB) relationship based on electrostatic (∆G(ele)), van der Waal (∆G(vdW)) and surface accessible surface area (SASA), which can express the activity of these artemisinin derivatives. It has been seen that ∆G(vdW) has most significant correlation to the activity (log RA) and electrostatic energy (∆G(ele)) has less significant correlation. It indicates that the binding of these artemisinin derivatives to PfATP6 is almost hydrophobic. Low levels of root mean square error for the majority of inhibitors establish the docking, Prime/MM-GBSA and eMBrAcE based prediction model is an efficient tool for generating more potent and specific inhibitors of PfATP6 by testing rationally designed lead compound based on artemisinin derivatization.