Skip to main content

Draft genome sequence of Aspergillus flavus isolate TERIBR1, a highly tolerant fungus to chromium stress



Aspergillus flavus isolate TERIBR1 was isolated from tannery sludge highly contaminated with chromium. During characterization process, it exhibited capability to adapt and grow in fungal growth media amended with chromium concentration as high as 250 mg/l. In order to understand the genetic underpinnings of the chromium tolerance trait, whole genome sequencing of the TERIBR1 genome was carried out. Information from the current genome will facilitate an understanding of the mechanisms underlying fungal adaptation to heavy metal stress and also heavy metal bioremediation.

Data description

Here, we report the draft genome sequence along with the assembly and annotation methods used for genome sequence of the A. flavus isolate TERIBR1. The draft genome assembly size is estimated at 37.7 Mb coding for 13,587 genes and has high similarity to the reference genome of A. flavus strain NRRL3357.


Several species of filamentous fungi have been identified for their bioaccumulation or biosorption potentials [1,2,3,4]. Reduced cost and environmental-toxicity through microbial bioremediation approach makes it favorable over the conventional methods [5]. The genome of several A. flavus strain have been reported previously The ability of the A. flavus isolate TERIBR1 to adapt and grow in tannery sludge highly contaminated with chromium inspired us to carry out its whole genome sequencing. The genome sequence reported here was utilized for comparative genomics study to understand the putative influence of the abundantly present non-synonymous SNP in TERIBR1 on the function of candidate genes involved in chromium tolerance [6].

Data description

Pure culture of A. flavus isolate TERIBR1 was recovered through an enrichment culture technique from tannery sludge [containing very high concentration of of Cr(III)] and molecularly characterized by the universal fungal primer set for Ascomycetes (ITS1: 5′ TCCGTAGGTGAACCTGCGG, 3′ (Eurofins India, Cat. No. 24-1023-5/6) and ITS4A: 5′ CGCCGTTACTGGGGCAATCCCTG 3′ (Eurofins India, Cat. No. 24-2002-1/6). Genomic DNA was extracted using the DNeasy plant maxi kit (QIAGEN, USA; cat. No. 68163). Using a whole-genome shotgun approach, two TruSeq paired-end (PE) libraries (insert sizes 180 bp and 500 bp) and a mate pair (MP) library (insert size ~ 5 Kb) was generated. An Illumina (HiSeq 2000) machine at a commercial facility (MOgene LC, USA) was used for sequencing. DNA libraries were loaded into Illumina flow-cells at concentrations of 1.4–1.75 pM. Cluster generation was performed in a cBOT automated cluster generation system. Real Time Analysis (RTA) software (rta_1–13) was used to process the image analysis and base calling. Sequencing of the DNA libraries yielded 5.4 Gb of PE reads and 2.6 Gb of MP reads. The raw reads were trimmed using Trimmomatic V 0.36 [7]. Quality-passed reads were assembled using the de novo genome assembler ALLPATHS-LG. PE reads with overlaps were first combined to form contigs. MP reads were used for gap filling in order to get sequences with minimal N’s and the longest length. Table 1 presents webpage links for genome assembly and annotation data files. The resulting 3,77,32,467 bp (100 X coverage) draft genome assembly [10] comprises of 322 contigs greater than 900 bp and has an N50 of 1,536,000 bp and an L50 of 9 contigs (Additional file 1). The GC content of the assembled genome is 48.30%. 225 out of 248 ultra-conserved eukaryotic genes were identified in the assembly through CEGMA ([8], Additional file 2). The MAKER v2.31.9 [9] genome annotation and curation pipeline predicted 13,587 protein coding genes as compared to 13,659 in NRRL3357. Using blastp search in the NCBI NR database, significant matches were identified for 11,120 protein-coding genes. An InterProScan analysis was also performed in order to further annotate the predicted genes with protein functional domains. 2551 proteins with InterProScan domains were identified (Additional file 3); major protein families included, Major facilitator superfamily (n = 334), fungal specific transcription factor domain (n = 190), Cytochrome P450 (n = 140), sugar (and other) transporters (n = 127), Protein kinase domain (n = 112), short chain dehydrogenase (n = 112) and fungal Zn(2)-Cys(6) binuclear cluster domain (n = 94) (Additional file 4). Genes were also annotated by using Blast2GO V5 basic [10] based on the term “biological function” in Gene Ontology (GO) (Additional file 5).

Table 1 Overview of data files


Illumina sequencing reads generated in this study were de novo assembled and annotated to understand the gene/protein repertoires in the chromium tolerant isolate of A. flavus. Since the whole genome sequencing project involved use of both PE and MP libraries for scaffold development, a high quality assembly with 100 X coverage could be generated. Therefore, we did not notice any serious limitations of the data.

Availability of data and materials

Genomic assembly of A. flavus TERIBR1; The additional files (Additional File 1, Additional File 2, Additional File 3, Additional File 4, Additional File 5) can be accessed openly on Figshare (



non-synonymous single nucleotide polymorphism




mate pair


  1. 1.

    Lacina C, Germain G, Spiros NA. Utilization of fungi for bio treatment of raw wastewaters. Afr J Biotechnol. 2003;2(12):620–30.

    Article  Google Scholar 

  2. 2.

    Volesky B, Holan ZR. Biosorption of heavy metals. Biotechnol Progr. 1995;11:235–50.

    CAS  Article  Google Scholar 

  3. 3.

    Huang CP, Huang CP. Application of Aspergillus oryzae and Rhizopus oryzae. Water Res. 1996;30:1985–90.

    CAS  Article  Google Scholar 

  4. 4.

    Zafar S, Aqil F, Ahmad I. Metal tolerance and biosorption potential of filamentous fungi isolated from metal contaminated agriculture soil. Biores Technol. 2007;98:2557–61.

    CAS  Article  Google Scholar 

  5. 5.

    Hansda A, Kumar V, Anshumali. A comparative review towards potential of microbial cells for heavy metal removal with emphasis on biosorption and bioaccumulation. World J Microbiol Biotechnol. 2016;32(10):170.

    Article  Google Scholar 

  6. 6.

    Jaiswar A, Varshney D, Adholeya A, Prasad P. Do environmentally induced DNA variations mediate adaptation in Aspergillus flavus exposed to chromium stress in tannery sludge? BMC Genomics. 2018;19(1):868.

    CAS  Article  Google Scholar 

  7. 7.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.

    CAS  Article  Google Scholar 

  9. 9.

    Cantarel BL, Korf I, Robb SM, Parra G, Ross E, Moore B, Holt C, Sánchez Alvarado A, Yandell M. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008;18(1):188–96.

    CAS  Article  Google Scholar 

  10. 10.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Bioproject: Aspergillus flavus strain: TERIBR1 Genome sequencing and assembly; 2019.

  11. 11.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Additional file 1. figshare; 2019.

  12. 12.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Additional file 2. figshare; 2019.

  13. 13.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Additional file 3. figshare; 2019.

  14. 14.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Additional file 4. figshare; 2019.

  15. 15.

    Prasad P, Jaiswar A, Srivastava D, Adholeya A. Additional file 5. figshare; 2019.

Download references


We acknowledge Mr. Aditya Gaur for providing support on symmetric multi-processing (SMP) platform with Unix-based operating system for bioinformatics analysis.

Previous data citation

The genomic assembly has been previously cited: Jaiswar A, Varshney D, Adholeya A, Prasad P. Do environmentally induced DNA variations mediate adaptation in Aspergillus flavus exposed to chromium stress in tannery sludge? BMC Genomics. 2018 Dec 4;19(1):868. PubMed PMID: 30509176; PubMed Central PMCID: PMC6278149.


Not applicable

Author information




PPS was involved in conceptualization of the study, genome study, manuscript writing and editing; AJ was involved in manuscript writing and data curation, DS was involved in data curation, AA was involved in conceptualization of the study, supervised the work towards purification of the microbial isolate and chromium tolerance activity assessment of TERIBR1. All authors read and approved the final manuscript

Corresponding author

Correspondence to Pushplata Prasad Singh.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Singh, P.P., Jaiswar, A., Srivastava, D. et al. Draft genome sequence of Aspergillus flavus isolate TERIBR1, a highly tolerant fungus to chromium stress. BMC Res Notes 12, 443 (2019).

Download citation


  • A. flavus isolate TERIBR1
  • Whole genome sequencing
  • de novo assembly