Next Article in Journal
Phytotoxic Metabolites Isolated from Neufusicoccum batangarum, the Causal Agent of the Scabby Canker of Cactus Pear (Opuntia ficus-indica L.)
Next Article in Special Issue
Signatures of TRI5, TRI8 and TRI11 Protein Sequences of Fusarium incarnatum-equiseti Species Complex (FIESC) Indicate Differential Trichothecene Analogue Production
Previous Article in Journal
2-Phenylethyl Isothiocyanate Exerts Antifungal Activity against Alternaria alternata by Affecting Membrane Integrity and Mycotoxin Production
Previous Article in Special Issue
The In Vivo and In Vitro Toxicokinetics of Citreoviridin Extracted from Penicillium citreonigrum
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 12
Issue 2
10.3390/toxins12020125
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Whole Genome Analysis Revealed the Genes Responsible for Citreoviridin Biosynthesis in Penicillium citreonigrum

by
Takumi Okano
1,†,
Naoki Kobayashi
1,*,†,
Kazuki Izawa
2,
Tomoya Yoshinari
3 and
Yoshiko Sugita-Konishi
1
1
Graduate School of Life and Environmental Sciences, Azabu University, Kanagawa 252-5201, Japan
2
Department of Computer Science, Tokyo Institute of Technology, Tokyo 152-8550, Japan
3
Division of Microbiology, National Institute of Health and Sciences, Kanagawa 210-9501, Japan
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Toxins 2020, 12(2), 125; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12020125
Submission received: 23 January 2020 / Revised: 12 February 2020 / Accepted: 13 February 2020 / Published: 15 February 2020
(This article belongs to the Special Issue Mycotoxins and Mycotoxin Producing Fungi: Genetic and Toxicity Aspects)
Browse Figures
Versions Notes

Abstract

:
Citreoviridin (CTV) is a mycotoxin that is produced by Aspergillus terreus, Eupenicillium ochrosalmoneum and Penicillium citreonigrum, and CTV has been detected in a wide range of cereal grains throughout the world. Furthermore, it is especially a serious problem in regions where rice is consumed as a staple food. Moreover, CTV is a well-known yellow rice toxin, and outbreaks of beriberi have occurred due to consumption of rice that is contaminated by CTV even in the recent years. Although CTV biosynthetic genes of A. terreus have been described, those of P. citreonigrum remain unclear, which is concerning since P. citreonigrum is the main cause of CTV contamination in rice. In the present study, we determined the draft genome of the P. citreonigrum strain IMI92228 and revealed the presence of all four genes that form a gene cluster and that are homologous to the CTV biosynthesis genes of A. terreus. The expression of these four homologous genes were highly correlated with CTV production, suggesting that they may play an important role in CTV biosynthesis in P. citreonigrum. We concluded that the gene cluster is a CTV biosynthesis cluster of P. citreonigrum. The findings will contribute to the understanding of the biosynthetic pathway of CTV and will ultimately lead to improvements in the CTV management of agricultural products.
Keywords:
biosynthesis gene cluster; citreoviridin; draft genome; Penicillium citreonigrum
Key Contribution: This is the first report on the draft genome of P. citreonigrum and we revealed that P. citreonigrum has a gene cluster that is homologous to the CTV biosynthesis gene cluster described for A. terreus.

1. Introduction

Citreoviridin (CTV) is one of the mycotoxins that contaminates a wide range of agricultural products, including rice. It is produced by Aspergillus terreus, Eupenicillium ochrosalmoneum and Penicillium citreonigrum [1,2,3]. Toxicological studies indicate that CTV could induce acute neurotoxic effects, such as respiratory failure, circulatory failure, paralysis and convulsions, in humans and experimental animals [4,5,6]. Previous studies demonstrate the relationship between CTV and several human diseases. Several decades ago, shoshin-kakke (cardiac beriberi) occurred in Japan, and it is related to the consumption of molded yellow rice, which is contaminated with CTV produced by P. citreonigrum [7]. Recently, outbreaks of beriberi have occurred in Brazil, and P. citreonigrum and CTV were detected in rice samples of related outbreaks [8,9]. CTV has also been implicated in Keshan disease, which is an endemic cardiomyopathy condition of China and North Korea, through an oxidative stress mechanism [10]. Additionally, it is possible that CTV is a risk factor for the development of atherosclerosis [11].
CTV is a highly reduced polyketide product, and its chemical structure was determined in the 1960s by Sakabe [12]. The structure is similar to aurovertins, which are potent inhibitors of mitochondrial respiration, and in vitro studies demonstrate that CTV inhibits triphosphate adenosine and thiamine diphosphate, which may reveal a potential mechanism for how CTV causes cardiac beriberi [13,14,15].
Recently, the gene cluster for CTV biosynthesis was identified in A. terreus ver. aureus [16]. Lin et al. revealed the cluster by a resistance-gene-driven genome mining method. To achieve self-resistance against its self-produced metabolites, fungi may occasionally harbor duplicated resistant targets within the biosynthetic gene cluster [17,18]. As such, an extra gene copy of the F1-ATPase β-subunit, a well-known target of CTV, was found and named ctvE. Furthermore, ctvE was located next to putative enzyme genes that are involved in the synthesis of CTV. A highly reducing polyketide synthase (HR-PKS) gene (ctvA), and the following genes for other putative enzymes are required to synthesize the tetrahydrofuran ring and methylated α-pyrone in CTV: methyltransferase (ctvB), flavin-dependent mono-oxygenase (ctvC) and hydrolase (ctvD). Heterologous expression in A. nidulans revealed that these four genes together are sufficient for CTV formation [16]. Nevertheless, the gene cluster for CTV biosynthesis in P. citreonigrum, which is a species that is believed to be responsible for yellow rice toxicity and recent outbreaks of beriberi, has still been unclear. In the present study, we hypothesized that homologous genes to the CTV biosynthesis genes in A. terreus are involved in the biosynthesis of CTV in P. citreonigrum, and as such, these genes are searched for in the genome.

2. Results

2.1. Whole Genome Sequence

To investigate whether P. citreonigrum has homologous genes to the CTV biosynthesis genes in A. terreus, we conducted a whole genome shotgun analysis using the P. citreonigrum strain IMI92228. In total, 76,116,858 paired-end raw reads were obtained. After quality filtering, 37,808,090 forward reads and 32,889,633 reverse reads remained. 32,731,394 reads had pairs and 5,076,696 reads were single. After the assemble and scaffolding, we obtained 79 scaffolds that were more than 500 bp each and approximately 27 Mbp in total length (Table 1).
From the open reading frame (ORF) prediction, 9805 ORFs were predicted, and analysis of single copy genes conserved in Eurotiomycetes showed that 97.3% (3936 genes out of 4046 genes) of genes were found in the obtained scaffolds. The estimated complete genome size was approximately 28 Mbp (27,997,905 bp). All the predicted ORFs were subjected to alignment using protein basic local alignment search tool (BLASTP) with the following five genes that are involved in the CTV biosynthesis of A. terreus and which form a gene cluster in A. terreus genome: ctvA, ctvB, ctvC, ctvD and ctvE (Figure 1a). As a result, we found that the predicted ORFs in the P. citreonigrum strain IMI92228 genome showed high homology to all five CTV biosynthesis genes in A. terreus (Table 2). Among these, the ORFs that showed high homology to ctvA, ctvB, ctvC and ctvD were located on scaffold 16 (2,344,991 bp) (Table 2 and Figure 1). We considered that these specific predicted ORFs were homologous to ctvA to ctvD in A. terreus and were named ci-ctvA to ci-ctvD (Accession number: LC517105 and LC517107 to LC517109). These predicted genes were also found to be arranged in the same order and direction as that of A. terreus, although there is a gene (g1457, LC517106) that cannot be found in A. terreus, between ci-ctvC and ci-ctvD (Figure 1b). The gene g1457 is similar to HC-toxin efflux carrier TOXA of Talaromyces islandicus. Meanwhile, the ORF that showed high homology to ctvE was separated from the others and located on scaffold 19 (2,964,612 bp) (Table 2 and Figure 1). There was only one ORF (g2666, LC517110) that showed high homology to ctvE in the P. citreonigrum genome and the ORF showed higher identity to one of the duplicated genes of the F1-ATPase β-subunit in A. terreus, other than ctvE (ATEG_07609, identity = 90.0%).

2.2. Citreoviridin Production Profile

P. citreonigrum was cultured on YES agar and YES agar with 5% NaCl, and CTV production was monitored for 2 weeks. CTV production on the YES agar peaked on day 7 (110.3 ± 11.7 ng/sample) and decreased on day 10 (Figure 2). CTV production was suppressed at a value lower than 1.1 ± 0.4 ng/sample on YES agar with 5% NaCl.

2.3. Gene Expression

To validate whether the homologous genes to the CTV biosynthesis genes of A. terreus are involved in CTV biosynthesis in P. citreonigrum, a gene expression analysis of these genes was carried out during 14-day cultures on YES agar and YES with 5% NaCl agar. The expression of ci-ctvA was highest at day 4 and then gradually decreased on YES agar (Figure 3a). The expressions of ci-ctvB, ci-ctvC and ci-ctvD peaked on day 7 (Figure 3b–d). On YES agar with 5% NaCl, the expression of ci-ctvA, ci-ctvB, ci-ctvC and ci-ctvD tended to decrease compared to those of YES agar.

3. Discussion

The genes responsible for CTV biosynthesis were previously identified in A. terreus ver. aureus [16], while those of P. citreonigrum, which is an important species that causes most of the world’s CTV contamination of rice, are veiled. Thus, in the present study, we aimed to identify the genes for CTV biosynthesis in the P. citreonigrum genome.
Biosynthesis of fungal secondary metabolites, including mycotoxins, are generally produced through numerous enzymatic steps [19,20]. The genes that encode the enzymes to produce the secondary metabolites are generally located near each other and form a biosynthetic gene cluster on the chromosome. For example, biosynthesis of aflatoxin, which is the most studied and well-known mycotoxin as a highly toxic and carcinogenic secondary metabolite, requires 28 genes to be produced, and these genes are clustered on chromosome 3 in A. flavus [21,22,23]. The genes for CTV biosynthesis in A. terreus also form a cluster (Figure 1a) [16]. In the present study, we determined the draft genome of P. citreonigrum strain IMI92228 and searched for genes that are homologous to the biosynthesis genes of CTV in A. terreus.
Next-generation sequencing of the P. citreonigrum strain IMI92228 yielded a total of 7611 Mbp, and the final assembly was found to contain 27 Mbp in 79 scaffolds. The total length of scaffolds and estimated genome size were relatively small compared to previously reported genomes of the Penicillium species [24,25,26,27], although these scaffolds code 97.3% of single copy genes that were conserved in Eurotiomycetes, which is enough to search for the homologous gene cluster that corresponds to the CTV biosynthetic cluster of A. terreus. The estimated genome size may be underestimated because only short inset sizes were estimated, and short read analysis may be insufficient to evaluate repeat regions and transposable elements of the genome. Further analysis on the genome structure of P. citreonigrum strain IMI92228 will provide more accurate information on genome size. Furthermore, we found genes that show high homology to all four CTV biosynthesis genes in A. terreus-ctvA to ctvD, and which are sufficient for CTV formation (Table 2). We named these genes ci-ctvA to ci-ctvD. These four predicted genes were located near each other on scaffold 16 and form a gene cluster, such as that in the A. terreus genome (Figure 1). Interestingly, the direction and order were also similar to those in the A. terreus genome except for the presence of an extra gene, g1457, between ci-ctvC and ci-ctvD, and which is not found in the A. terreus genome. Generally, in the biosynthetic gene clusters for the same secondary metabolite, the order of the homologous genes differs among species, which suggests that the precise gene order is not crucial for function. For example, the gene order markedly differs between the aflatoxin biosynthetic gene cluster in A. flavus and the sterigmatocystin cluster in A. nidulans, although A. flavus also produces sterigmatocystin as an intermediate of aflatoxin [28,29]. The gene order of the CTV biosynthesis gene cluster was tremendously similar when comparing A. terreus and P. citreonigrum (Figure 1), suggesting that the origin of them may be same. This is evolutionally interesting and the function of the genes may be conserved.
On the other hand, the gene that shows high homology to ctvE was located on scaffold 19, which is different from the scaffold that ci-ctvA to ci-ctvD were located on (Figure 1, g2666). Scaffold 16 and 19 can be combined if more precise genome sequencing emerges. However, these genes that are found in the middle of each scaffold imply that ci-ctvA to ci-ctvD and g2666 do not form a gene cluster in P. citreonigrum. ctvE is thought to be a F1-ATPase β-subunit, which is a known target of citreoviridin [16]. Lin et al. [16] demonstrated that A. terreus has an extra copy of the F1-ATPase β-subunit, known as ctvE, near the CTV biosynthesis cluster, and they predicted that this copy ensures self-resistance against CTV. The genome analysis of the P. citreonigrum strain IMI9228 that was performed in this study revealed that there was only one F1-ATPase β-subunit gene (g2666), which is more similar to ATEG_07609, also known as the F1-ATPase β-subunit, rather than ctvE. This indicates that the F1-ATPase β-subunit does not duplicate in the P. citreonigrum genome. Therefore, P. citreonigrum should protect itself from CTV by a different method than by harboring an extra F1-ATPase β-subunit gene. The genome analysis that was performed in this study predicted a gene that was not found in the A. terreus genome, specifically between ci-ctvC and ci-ctvD. This gene is similar to the HC-toxin efflux carrier TOXA and has an MFS_Azr1_MDR_like domain, which is shared among multidrug resistance transporters of the major facilitator superfamily [30]. Efflux-mediated toxin resistance to trichothecenes is well-regarded in certain trichothecene-producing Fusarium species [31,32]. P. citreonigrum may achieve self-resistance against CTV by possessing an MFS transporter in the biosynthesis cluster of CTV. These results strongly suggest that ci-ctvA to ci-ctvD are involved in the CTV biosynthesis in P. citreonigrum.
To validate the responsibility of the genes for CTV biosynthesis in P. citreonigrum, expression analysis was performed (Figure 3). The expression level of ci-ctvA was highest at day 4, that is immediately before CTV production peak and then gradually decreased on YES agar (Figure 2 and Figure 3a). Since ci-ctvA is believed to be a HR-PKS gene, in which the encoding enzyme accepts acetyl-CoA as the starter unit for CTV and catalyzes iterations of both the malonyl-CoA extension and SAM-dependent methylation at the beginning of CTV biosynthesis [16], the expression profile is correlated with CTV production. On the other hand, the expressions of ci-ctvB, ci-ctvC and ci-ctvD peaked on day 7, which is correlated with CTV production. These three genes operate in the steps of CTV biosynthesis after the ci-ctvA step, and their expression profiles correlate with CTV production. Moreover, expressions of all four genes decreased on YES agar with 5% NaCl, whereby CTV production was suppressed. These results ascertain the potential of ci-ctvA to ci-ctvD as key genes for CTV biosynthesis in P. citreonigrum, although further analysis, such as deficient mutant analyses, must be conducted to confirm responsibility for CTV production.

4. Conclusions

The draft genome of the P. citreonigrum strain IMI92228 was first determined and the genome analysis revealed the presence of a gene cluster quite similar to the CTV biosynthetic gene cluster of A. terreus. In the cluster, there were all four genes that were believed to be sufficient for CTV formation in A. terreus, and the gene expression profiles were highly correlated with CTV production, which indicate involvement in CTV production in P. citreonigrum. Furthermore, it was found that the cluster contains a putative transporter, which may be involved in self-resistance against CTV. Lastly, we conclude that the gene cluster found is a CTV biosynthesis cluster in P. citreonigrum. This finding of a gene cluster for CTV in P. citreonigrum will expedite efforts towards better understanding of the biosynthetic pathway of CTV, which ultimately leads to the management of CTV contamination in agricultural products, especially rice.

5. Materials and Methods

5.1. Strain

Penicillium citreonigrum strain IMI92228 was obtained from the culture collection of Centre for Agriculture and Bioscience International (CABI) in the United Kingdom and stored on potato dextrose agar (PDA; Eiken, Tokyo, Japan) slant medium.

5.2. Whole Genome Sequencing

The P. citreonigrum strain IMI92228 was cultured on Potato Dextrose (Difco Laboratories, Franklin Lakes, NJ, USA) liquid medium for 2 days at 25 °C. Genomic DNA was extracted using DNeasy (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions and quantified using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, MA, USA). Library construction and sequencing for the Illumina HiSeq 2500 was provided as a custom service of Eurofins Genomics K. K. (Tokyo, Japan). Genomic DNA samples were sheared to 350 bp fragments by sonication with a sonicator (Covaris, MA, USA). The resulting DNA fragments were processed for adaptor ligation and underwent amplification to generate DNA libraries. Prepared libraries were subjected to paired-end 2 × 100 bp sequencing with the HiSeq 2500 platform (Illumina, San Diego, CA, USA). Obtained raw reads were trimmed and filtered by quality. Each read was trimmed from both ends until all bases in the read surpassed 20 in a Phred quality score. De novo assembly was performed using Platanus 1.2.4 [33] under default settings. The scaffolds that are shorter than 500 bp were discarded, and 80 scaffolds were subjected to the following analysis. For ORF prediction, a web service version of AUGUSTUS (http://bioinf.uni-greifswald.de/webaugustus/prediction/create) [34] was used. The settings of A. fumigatus were used as the reference for ORF prediction and the predicted ORFs were subjected to BLAST. The predicted ORFs that aligned to the genes that were involved in the CTV biosynthesis of A. terreus were manually checked and identified as the CTV biosynthesis genes in P. citreonigrum. The genome completeness of the P. citreonigrum strain IMI92228 were estimated by a single copy gene set of Eurotiomycetes using the BUSCO pipeline [35].

5.3. Citreoviridin Production

P. citreonigrum strain IMI92228 was cultured on either YES agar or YES–5% NaCl agar for 14 days at 25 °C. The YES agar had the following formulation: 4 g/L of yeast extract (Difco Laboratories), 20 g/L of sucrose (Mitsui Sugar Co. Ltd., Tokyo, Japan) and 15 g/L of agar (Difco Laboratories) in distilled water. CTV was extracted from the plug with 1 mL of chromophore. After sonication for 20 min, the chloroform layer was collected and dried up under N2. The extract was resolved with 1 mL of acetonitrile. Detection and quantification were performed using liquid chromatography–tandem mass spectrometry (LC-MS/MS).

5.4. LC–MS/MS Conditions

LC–MS/MS analysis was performed with a Triple Quad 4500 system (AB Sciex, Foster City, CA) equipped with an ESI source and an LC-20A series high-performance LC system (Shimadzu Corp., Kyoto, Japan). The column used was a 150 mm × 2.1 mm i.d., 3 μm, InertSustain C18 (GL Sciences Inc., Tokyo, Japan). The mobile phase was a binary gradient of solvent A (water containing 0.1% formic acid) and solvent B (acetonitrile containing 0.1% formic acid) programmed as follows: at 0 min, 30% B; at 8 min, 80% B; and at 10 min, 80% B. The total run time was 15 min, which included 5 min of equilibration. The flow rate was set at 0.2 mL/min. Chromatographic separation was achieved at 40 °C. The injection volume was 2 μL. The ESI source was operated at 300 °C in the positive ionization mode. Other MS parameters were as follows: curtain gas, 40 psi; ion spray voltage, 5500 V; nebulizer gas (GS1), 30 psi; turbo heater gas (GS2), 80 psi; and collision-activated dissociation gas, 7 (arbitrary units). The multiple reaction monitoring transitions of CTV were 403 [M+H]+ to 315 (collision energy, 9 eV) and 403 [M+H]+ to 139 (collision energy, 31 eV). The retention time of the CTV in the sample solution (8.8 min) was consistent with that of the analytical standard. The standard of CTV was purchased from FERMENTEK Ltd. (Jerusalem, Israel).

5.5. Gene Expression

The P. citreonigrum strain IMI92228 was cultured on YES agar or YES–5% NaCl agar for 14 days at 25 °C. The mycelia were stored in RNAlaterTM Stabilization Solution (Thermo Fisher Scientific, Waltham, MA USA) at −80 °C. Mycelia were frozen by liquid nitrogen and clashed using an SK mill (Tokken, inc., Chiba, Japan). Total RNA were isolated using the Maxwell® RSC instrument and SimplyRNA tissue Kit (Promega, Madison, WI, USA), according to the manufacturer’s instructions. RNA quantity and quality were assessed using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific). cDNA synthesis was performed using PrimeScript™ RT reagent Kit with gDNA Eraser (Perfect Real Time) (TaKaRa, Shiga, Japan) using 250 ng of total RNA, according to the manufacturer’s instructions.
qPCR was performed with a Thermal Cycler Dice® Real Time System III (TaKaRa) using TB Green® Premix DimerEraser™ (Perfect Real Time) (TaKaRa). Each 20 μl of reaction mixture contained 2 μL of cDNA sample, 10 μL of 2 × TB Green Fast qPCR Mix and 100 nM of each primer. The amplification condition was one cycle for 2 min at 50 °C and 95 °C each, followed by 45 cycles at both 95 °C for 15 s and at 60 °C for 30 s. Melting curve analysis was performed to ensure the specificity of the PCR products. The PCR for each target gene was performed in duplicate for every sample. The data from qPCR were normalized to β-tubulin amplification using the comparative threshold cycle (Ct) method [36]. The primers used for amplification are listed in Table 3. The experiments were repeated with three independent cultures. A Student’s t-test was used to determine the differences in the expression levels compared to those of a 4-day culture sample on YES agar.

5.6. Statistical Analyses

The mean and SD of each value were calculated. Differences between cultures on YES agar and those on YES agar with 5% NaCl were analyzed using a Student’s t-test. Statistical p values of <0.05 were considered statistically significant.

Author Contributions

Conceptualization, T.O., N.K., K.I., T.Y. and Y.S.-K.; methodology, T.O., N.K., K.I., T.Y. and Y.S.-K.; investigation, T.O., N.K. K.I. and T.Y.; formal analysis, T.O. and N.K.; writing—original draft preparation, T.O., N.K. and Y.S.-K.; writing—review and editing, all authors contributed to writing and reviewing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We would like to thank Editage (www.editage.com) for English language editing.

Conflicts of Interest

The authors declare no conflict of interest in association with this study.

References

  1. Shiratori, N.; Kobayashi, N.; Tulayakul, P.; Sugiura, Y.; Takino, M.; Endo, O.; Sugita-Konishi, Y. Occurrence of Penicillium brocae and Penicillium citreonigrum, which produce a mutagenic metabolite and a mycotoxin citreoviridin, respectively, in selected commercially available rice grains in Thailand. Toxins 2017, 9, 194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lima, H.C.A.V.; Porto, E.A.S.; Marins, J.R.P.; Alves, R.M.; Machado, R.R.; Braga, K.N.L.; de Paiva, F.B.; Carmo, G.M.I.; Faria Silva e Santelli, A.C.; Sobel, J. Outbreak of beriberi in the state of Maranhão, Brazil: Revisiting the mycotoxin aetiologic hypothesis. Trop. Doct. 2010, 40, 95–97. [Google Scholar] [CrossRef] [PubMed]
  3. Wicklow, D.T.; Cole, R.J. Citreoviridin in standing corn infested by Eupenicillium ochrosalmoneum. Mycologia 1984, 76, 959–961. [Google Scholar] [CrossRef]
  4. Morrissey, R.E.; Vesonder, R.F. Teratogenic potential of the mycotoxin, citreoviridin, in rats. Food Chem. Toxicol. 1986, 24, 1315–1320. [Google Scholar] [CrossRef]
  5. Datta, S.C.; Ghosh, J.J. Effect of citreoviridin, a mycotoxin from Penicillium citreoviride, on kinetic constants of acetylcholinesterase and ATPase in synaptosomes and microsomes from rat brain. Toxicon 1981, 19, 555–562. [Google Scholar] [CrossRef]
  6. Ueno, Y. Temperature-dependent production of citreoviridin, a neurotoxin of Penicillium citreo-viride Biourge. Jpn. J. Exp. Med. 1972, 42, 107–114. [Google Scholar] [PubMed]
  7. Kushiro, M. Historical review of researches on yellow rice and mycotoxigenic fungi adherent to rice in Japan. JSM Mycotoxins 2015, 65, 19–23. [Google Scholar] [CrossRef] [Green Version]
  8. Padilha, E.M.; Fujimori, E.; Borges, A.L.V.; Sato, A.P.S.; Gomes, M.N.; dos Remédios Freitas Carvalho Branco, M.; dos Santos, H.J.; Junior, N.L. Perfil epidemiológico do beribéri notificado de 2006 a 2008 no Estado do Maranhão, Brasil. Cad. Saude Publica 2011, 27, 449–459. [Google Scholar] [CrossRef]
  9. Da Rocha, M.W.; Resck, I.S.; Caldas, E.D. Purification and full characterisation of citreoviridin produced by Penicillium citreonigrum in yeast extract sucrose (YES) medium. Food Addit. Contam. Part A Chem. Anal. Control. Exp. Risk Assess. 2015, 32, 584–595. [Google Scholar] [CrossRef]
  10. Sun, S. Chronic exposure to cereal mycotoxin likely citreoviridin may be a trigger for Keshan disease mainly through oxidative stress mechanism. Med. Hypotheses 2010, 74, 841–842. [Google Scholar] [CrossRef]
  11. Hou, H.; Zhou, R.; Jia, Q.; Li, Q.; Kang, L.; Jiao, P.; Li, D.; Jiang, B. Citreoviridin enhances tumor necrosis factor-α-induced adhesion of human umbilical vein endothelial cells. Toxicol. Ind. Health 2015, 31, 193–201. [Google Scholar] [CrossRef] [PubMed]
  12. Sakabe, N.; Goto, T.; Hirata, Y. The structure of citreoviridin, a toxic compound produced by P. citreoviride molded on rice. Tetrahedron Lett. 1964, 5, 1825–1830. [Google Scholar] [CrossRef]
  13. Linnett, P.E.; Mitchell, A.D.; Osselton, M.D.; Mulheirn, L.J.; Beechey, R.B. Citreoviridin, a specific inhibitor of the mitochondrial adenosine triphosphatase. Biochem. J. 1978, 170, 503–510. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sayood, S.F.; Suh, H.; Wilcox, C.S.; Schuster, S.M. Effect of citreoviridin and isocitreoviridin on beef heart mitochondrial ATPase. Arch. Biochem. Biophys. 1989, 270, 714–721. [Google Scholar] [CrossRef]
  15. Cataldi de Flombaum, M.A.; Stoppani, A.O.M. Influence of efrapeptin, aurovertin and citreoviridin on the mitochondrial adenosine triphosphatase from Trypanosoma cruzi. Mol. Biochem. Parasitol. 1981, 3, 143–155. [Google Scholar] [CrossRef]
  16. Lin, T.S.; Chiang, Y.M.; Wang, C.C.C. Biosynthetic pathway of the reduced polyketide product citreoviridin in Aspergillus terreus var. aureus revealed by heterologous expression in Aspergillus nidulans. Org. Lett. 2016, 18, 1366–1369. [Google Scholar] [CrossRef]
  17. Tang, X.; Li, J.; Millán-Aguiñaga, N.; Zhang, J.J.; O’Neill, E.C.; Ugalde, J.A.; Jensen, P.R.; Mantovani, S.M.; Moore, B.S. Identification of thiotetronic acid antibiotic biosynthetic pathways by target-directed genome mining. ACS Chem. Biol. 2015, 10, 2841–2849. [Google Scholar] [CrossRef]
  18. Regueira, T.B.; Kildegaard, K.R.; Hansen, B.G.; Mortensen, U.H.; Hertweck, C.; Nielsen, J. Molecular basis for mycophenolic acid biosynthesis in Penicillium brevicompactum. Appl. Environ. Microbiol. 2011, 77, 3035–3043. [Google Scholar] [CrossRef] [Green Version]
  19. Osbourn, A. Secondary metabolic gene clusters: Evolutionary toolkits for chemical innovation. Trends Genet. 2010, 26, 449–457. [Google Scholar] [CrossRef]
  20. Keller, N.P.; Turner, G.; Bennett, J.W. Fungal secondary metabolism—From biochemistry to genomics. Nat. Rev. Microbiol. 2005, 3, 937–947. [Google Scholar] [CrossRef]
  21. Yu, J.; Chang, P.K.; Ehrlich, K.C.; Cary, J.W.; Bhatnagar, D.; Cleveland, T.E.; Payne, G.A.; Linz, J.E.; Woloshuk, C.P.; Bennett, J.W. Clustered pathway genes in aflatoxin biosynthesis. Appl. Environ. Microbiol. 2004, 70, 1253–1262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Payne, G.A.; Nierman, W.C.; Wortman, J.R.; Pritchard, B.L.; Brown, D.; Dean, R.A.; Bhatnagar, D.; Cleveland, T.E.; Machida, M.; Yu, J. Whole genome comparison of Aspergillus flavus and A. oryzae. Med. Mycol. 2006, 44, S9–S11. [Google Scholar] [CrossRef] [PubMed]
  23. Bennett, J.W.; Klich, M. Mycotoxins. Clin. Microbiol. Rev. 2003, 16, 497–516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yu, J.; Jurick, W.M.; Cao, H.; Yin, Y.; Gaskins, V.L.; Losada, L.; Zafar, N.; Kim, M.; Bennett, J.W.; Nierman, W.C. Draft Genome Sequence of Penicillium expansum Strain R19, Which Causes Postharvest Decay of Apple Fruit. Genome Announc. 2014, 2, e00635. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yu, J.; Wu, G.; Jurick, W.M.; Gaskins, V.L.; Yin, Y.; Yin, G.; Bennett, J.W.; Shelton, D.R. Genome sequence of Penicillium solitum RS1, which causes postharvest apple decay. Genome Announc. 2016, 4, e00363. [Google Scholar] [CrossRef] [Green Version]
  26. Van Den Berg, M.A.; Albang, R.; Albermann, K.; Badger, J.H.; Daran, J.M.; Driessen, A.J.; Garcia-Estrada, C.; Fedorova, N.D.; Harris, D.M.; Heijne, W.H.M.; et al. Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum. Nat. Biotechnol. 2008, 26, 1161–1168. [Google Scholar] [CrossRef] [Green Version]
  27. Mardones, W.; Di Genova, A.; Cortés, M.P.; Travisany, D.; Maass, A.; Eyzaguirre, J. The genome sequence of the soft-rot fungus Penicillium purpurogenum reveals a high gene dosage for lignocellulolytic enzymes. Mycology 2018, 9, 59–69. [Google Scholar] [CrossRef] [Green Version]
  28. Carbone, I.; Ramirez-Prado, J.H.; Jakobek, J.L.; Horn, B.W. Gene duplication, modularity and adaptation in the evolution of the aflatoxin gene cluster. BMC Evol. Biol. 2007, 7, 111. [Google Scholar] [CrossRef] [Green Version]
  29. Brown, D.W.; Yu, J.H.; Kelkar, H.S.; Fernandes, M.; Nesbitt, T.C.; Keller, N.P.; Adams, T.H.; Leonard, T.J. Twenty-five coregulated transcripts define a sterigmatocystin gene cluster in Aspergillus nidulans. Proc. Natl. Acad. Sci. USA 1996, 93, 1418–1422. [Google Scholar] [CrossRef] [Green Version]
  30. Quistgaard, E.M.; Löw, C.; Guettou, F.; Nordlund, P. Understanding transport by the major facilitator superfamily (MFS): Structures pave the way. Nat. Rev. Mol. Cell Biol. 2016, 17, 123–132. [Google Scholar] [CrossRef]
  31. Wuchiyama, J.; Kimura, M.; Yamaguchi, I. A Trichothecene Efflux Pump Encoded by Tri102 in the Biosynthetic Gene Cluster of Fusarium graminearum. J. Antibiot. 2000, 53, 196–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Alexander, N.J.; McCormick, S.P.; Hohn, T.M. TRI12, a trichothecene efflux pump from Fusarium sporotrichioides: Gene isolation and expression in yeast. Mol. Gen. Genet. 1999, 261, 977–984. [Google Scholar] [CrossRef] [PubMed]
  33. Kajitani, R.; Toshimoto, K.; Noguchi, H.; Toyoda, A.; Ogura, Y.; Okuno, M.; Yabana, M.; Harada, M.; Nagayasu, E.; Maruyama, H.; et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 2014, 24, 1384–1395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Stanke, M.; Morgenstern, B. AUGUSTUS: A web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res. 2005, 33, W465–W467. [Google Scholar] [CrossRef] [Green Version]
  35. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  37. Yin, H.B.; Chen, C.H.; Kollanoor-Johny, A.; Darre, M.J.; Venkitanarayanan, K. Controlling Aspergillus flavus and Aspergillus parasiticus growth and aflatoxin production in poultry feed using carvacrol and trans-cinnamaldehyde. Poult. Sci. 2015, 94, 2183–2190. [Google Scholar] [CrossRef]
Toxins 12 00125 g001 550
Figure 1. The gene clusters for citreoviridin biosynthesis in Aspergillus terreus (a) and Penicillium citreonigrum (b). Black arrows indicate the genes of the enzymes involved in citreoviridin biosynthesis. Striped arrows indicate the gene for the F1-ATPase β-subunit. Gray arrow indicates a gene for the putative transporter. Scale bar indicates 1 kb.
Figure 1. The gene clusters for citreoviridin biosynthesis in Aspergillus terreus (a) and Penicillium citreonigrum (b). Black arrows indicate the genes of the enzymes involved in citreoviridin biosynthesis. Striped arrows indicate the gene for the F1-ATPase β-subunit. Gray arrow indicates a gene for the putative transporter. Scale bar indicates 1 kb.
Toxins 12 00125 g001
Toxins 12 00125 g002 550
Figure 2. Citreoviridin production by the Penicillium citreonigrum, strain IMI92228. Penicillium citreonigrum was cultivated on YES agar (black circle) and YES agar with 5% NaCl (gray triangle) at 25 °C. Error bars represent the standard error of three replicates. Different letters indicate statistically significant differences (p < 0.05).
Figure 2. Citreoviridin production by the Penicillium citreonigrum, strain IMI92228. Penicillium citreonigrum was cultivated on YES agar (black circle) and YES agar with 5% NaCl (gray triangle) at 25 °C. Error bars represent the standard error of three replicates. Different letters indicate statistically significant differences (p < 0.05).
Toxins 12 00125 g002
Toxins 12 00125 g003 550
Figure 3. Gene expression profile of the citreoviridin biosynthesis genes of Penicillium citreonigrum during citreoviridin production: (a) ci-ctvA, (b) ci-ctvB, (c) ci-ctvC, (d) ci-ctvD. The x-axis denotes culture time (day) and the y-axis indicates the relative expression level of each gene. Black bars represent the expression of the culture on YES agar and gray bars represents that on YES agar with 5% NaCl. Error bars represent the standard error of three replicates. Asterisks indicate statistically significant differences between cultures on YES agar and YES agar with 5% NaCl (p < 0.05).
Figure 3. Gene expression profile of the citreoviridin biosynthesis genes of Penicillium citreonigrum during citreoviridin production: (a) ci-ctvA, (b) ci-ctvB, (c) ci-ctvC, (d) ci-ctvD. The x-axis denotes culture time (day) and the y-axis indicates the relative expression level of each gene. Black bars represent the expression of the culture on YES agar and gray bars represents that on YES agar with 5% NaCl. Error bars represent the standard error of three replicates. Asterisks indicate statistically significant differences between cultures on YES agar and YES agar with 5% NaCl (p < 0.05).
Toxins 12 00125 g003
Table
Table 1. The assemble status of Penicillium citreonigrum strain IMI92228.
Table 1. The assemble status of Penicillium citreonigrum strain IMI92228.
Number of ScaffoldsMinimum Length of ScaffoldMaximum Length of ScaffoldN50Sum of Length
79503 bp2,963,310 bp1,448,320 bp27,997,905 bp
Table
Table 2. The result of homology search for citreoviridin biosynthesis genes by protein basic local alignment search tool.
Table 2. The result of homology search for citreoviridin biosynthesis genes by protein basic local alignment search tool.
Query Gene 1Hit Open Reading Frame (ORF)E-valueIdentity
(%)		Coverage
(%)
NameSize
(aa)		NumberSize
(aa)		ScaffoldDirectionPosition in Scaffold
ctvA2436g1460251316+519,757..527,298070.499.9
ctvB228g145923316-518,431..519,1323E-14379.9100
ctvC473g145647216-512,588..514,173077.497.9
ctvD354g145834116-517,137..518,2298E-14659.999.7
ctvE468g266651919+939,425..941,478079.495.9
1 The citreoviridin biosynthesis genes in Aspergillus terreus were used as a query.
Table
Table 3. Primers used for q PCR.
Table 3. Primers used for q PCR.
GenePrimer NameSequence (5′ to 3′)Reference
ctvActvA_1FAGCGTGGCATGATTACACCAAACCthis study
ctvA_1RCAACGTCGGCCATTGAAGAACCTC
ctvBctvB_1FTGTCTGAGAAAGGCTGCCAATCGTGthis study
ctvB_1RCAGCACGTACATCAGCGAGATGGA
ctvCctvC_1FCCAATACCGCCATGGAAGCAGthis study
ctvC_1RGCAAGCGCTCGTTCAATCGTATC
ctvDctvD_1FGGCGAATCTCTTGGCAGACATCthis study
ctvD_1RCCACAGCAAGAAACCACTCATCC
ctvEctvE_1FGTGACTCCAAGGTGTCTCTGthis study
ctvE_1RCAATGAAGAGCAGCACATCC
β-tubulinβ-tublinFCGTGTCGGCGACCAGTTC[37]
β-tublinRCCTCACCAGTGTACCAATGCA

Share and Cite

MDPI and ACS Style

Okano, T.; Kobayashi, N.; Izawa, K.; Yoshinari, T.; Sugita-Konishi, Y. Whole Genome Analysis Revealed the Genes Responsible for Citreoviridin Biosynthesis in Penicillium citreonigrum. Toxins 2020, 12, 125. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12020125

AMA Style

Okano T, Kobayashi N, Izawa K, Yoshinari T, Sugita-Konishi Y. Whole Genome Analysis Revealed the Genes Responsible for Citreoviridin Biosynthesis in Penicillium citreonigrum. Toxins. 2020; 12(2):125. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12020125

Chicago/Turabian Style

Okano, Takumi, Naoki Kobayashi, Kazuki Izawa, Tomoya Yoshinari, and Yoshiko Sugita-Konishi. 2020. "Whole Genome Analysis Revealed the Genes Responsible for Citreoviridin Biosynthesis in Penicillium citreonigrum" Toxins 12, no. 2: 125. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12020125

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop