Next Article in Journal
Botulinum Neurotoxins (BoNTs)—Antibody and Vaccine
Next Article in Special Issue
In Vitro and in Field Response of Different Fungicides against Aspergillus flavus and Fusarium Species Causing Ear Rot Disease of Maize
Previous Article in Journal
Fungal Endophytes Control Fusarium graminearum and Reduce Trichothecenes and Zearalenone in Maize
Previous Article in Special Issue
Phylogeny and Mycotoxin Characterization of Alternaria Species Isolated from Wheat Grown in Tuscany, Italy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 10
Issue 12
10.3390/toxins10120494
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:
Article

Analysis of the Fungal Community in Ziziphi Spinosae Semen through High-Throughput Sequencing

by
Mengyue Guo
,
Wenjun Jiang
,
Jiaoyang Luo
,
Meihua Yang
and
Xiaohui Pang
*
Key Lab of Chinese Medicine Resources Conservation, State Administration of Traditional Chinese Medicine of the People’s Republic of China, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100193, China
*
Author to whom correspondence should be addressed.
Toxins 2018, 10(12), 494; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120494
Submission received: 27 October 2018 / Revised: 14 November 2018 / Accepted: 23 November 2018 / Published: 25 November 2018
(This article belongs to the Special Issue 2nd MycoKey International Conference - Integrated Solutions for Mycotoxin Management)
Browse Figures
Review Reports Versions Notes

Abstract

:
Ziziphi Spinosae Semen (ZSS) has been widely used in traditional Chinese medicine system for decades. Under proper humidity and temperature, ZSS is easily contaminated by fungi and mycotoxins during harvest, storage, and transport, thereby posing a considerable threat to consumer health. In this study, we first used the Illumina MiSeq PE250 platform and targeted the internal transcribed spacer 2 sequences to investigate the presence of fungi in moldy and normal ZSS samples collected from five producing areas in China. Results showed that all 14 samples tested were contaminated by fungi. Ascomycota was the dominant fungus at the phylum level, accounting for 64.36–99.74% of the fungal reads. At the genus level, Aspergillus, Candida, and Wallemia were the most predominant genera, with the relative abundances of 13.52–87.87%, 0.42–64.56%, and 0.06–34.31%, respectively. Meanwhile, 70 fungal taxa were identified at the species level. Among these taxa, three potential mycotoxin-producing fungi, namely, Aspergillus flavus, A. fumigatus, and Penicillium citrinum that account for 0.30–36.29%, 0.04–7.37%, and 0.01–0.80% of the fungal reads, respectively, were detected in all ZSS samples. Moreover, significant differences in fungal communities were observed in the moldy and normal ZSS samples. In conclusion, our results indicated that amplicon sequencing is feasible for the detection and analysis of the fungal community in the ZSS samples. This study used a new approach to survey the fungal contamination in herbal materials. This new approach can provide early warning for mycotoxin contamination in herbal materials, thereby ensuring drug efficacy and safety.
Keywords:
Ziziphi Spinosae Semen; fungal contamination; amplicon sequencing; toxigenic fungi
Key Contribution: We utilized the amplicon sequencing platform and targeted the internal transcribed spacer 2 (ITS2) sequences to investigate the fungal contamination in herbal materials for the first time. A total of 70 fungal species were accurately identified in Ziziphi Spinosae Semen, among which three (i.e., Aspergillus flavus, A. fumigatus, and Penicillium citrinum) were potential toxigenic fungi.

1. Introduction

Herbal medicines, which are commonly used to prevent, diagnose, and cure diseases, have played an important role in health care since ancient times [1]. However, many studies have described the occurrence of mycotoxins in medicinal plants and herbal medicines from various countries [2,3,4], thereby attracting considerable attention worldwide due to drug efficacy and safety. For instance, 58 (8.29%) and 17 (2.43%) of the 700 herbal medicine samples in South Korea are aflatoxin B1 (AFB1)- and total aflatoxin (AF)-positive, respectively, and the AFB1 (up to 73.27 mg/kg) and total aflatoxin contents (up to 108.42 mg/kg) in some samples exceeded the legal limits (10 mg/kg) [2]. Both AFs and ochratoxin A (OTA) are detected in all Glycyrrhiza uralensis samples (six moldy and nine normal samples) that were collected from different areas in China [5]. Moreover, the transfer rate of OTA and AFs were investigated in decoctions of herbal medicines. The transfer rate of OTA from the three herbal medicines, namely Trichosanthis Semen, Eucommiae Cortex, and Rubi Fructus, to decoctions is 12.72–61.33% [6]. Considerable transfer rates of AFB1 (4.37–26.37%), AFB2 (9.64–47.68%), and AFG2 (7.26–115.36%) are also observed from the five herbal medicines (i.e., Lilii Bulbus, Hordei Fructus Germinatus, Nelumbinis Semen, Polygalae Radix, and Bombyx Batryticatus) to decoctions [7]. Chronic toxicity is the most common form of mycotoxicoses and is caused by low-dose exposure over an extended period, leading to cancers and other irreversible effects [8,9]. The presence of mycotoxins in the decoctions poses a considerable direct threat to consumer safety, because the ingestion of even extremely small amounts of mycotoxins can lead to many diseases or even death in humans and animals [10]. Among the approximately 400 recognized mycotoxins [11], AFs, OTA, fumonisins (FBs), zearalenone (ZEA), deoxynivalenol (DON), and patulin are important contaminants on humans in terms of health perspective [12,13]. AFB1 has been classified as a group 1 carcinogen by the International Agency for Research on Cancer (IARC) [14] and is risk factors for human hepatocellular carcinoma [15]. Meanwhile, OTA has been classified in Group 2B as possibly carcinogenic to humans by IARC [14].
Some species from the Aspergillus and Penicillium genera are AFs and OTA producers [3,16]. The detection of toxigenic fungi in herbal medicines has received considerable attention. Singh et al. [17] found that Aspergillus is the most prevailing genus that infects the raw materials of six medicinal plants from the herbal markets of Varanasi, and 13 out of 32 isolates of Aspergillus flavus were AFB1-positive. Another study showed that 90% of the 30 medicinal plant samples from Pakistan are contaminated with molds. The frequently isolated fungi predominantly consist of A. flavus, A. niger, A. parasiticus, and Penicillium spp., and 31% of the 47 isolates tested are toxigenic [3]. Fungal contamination in herbal medicines occurs during all the processes, including growing, harvesting, cleaning, transporting, and storage [18,19,20,21]. Under favorable temperature and humidity conditions, toxigenic fungi belonging to Aspergillus, Penicillium, and Fusarium produce mycotoxins [13]. Among these fungi, 21 species from Flavi, Ochraceorosei, and Nidulantes in Aspergillus produce AFs, while various Aspergillus and Penicillium species produce OTA [22,23]. The occurrence of toxigenic fungi in herbal medicines has the potential to produce mycotoxins. Therefore, the simultaneous detection of fungi, especially toxigenic fungi, in herbal medicines is important for early warning for mycotoxin contamination.
Ziziphi Spinosae Semen (ZSS, suanzaoren in Chinese), which is derived from the dried ripe seed of Ziziphus jujuba Mill. var. spinosa (Bunge) Hu ex H. F. Chou, is widely used to treat insomnia and palpitation in traditional Chinese medicine [24]. Pharmacological studies have demonstrated that ZSS regulates immune function, ameliorates learning and memory, and has a hypnotic-sedative effect [25,26,27]. Ziziphi Spinosae Semen is mainly produced in the Henan, Shanxi, Liaoning, and Hebei provinces of China and sold throughout the country. ZSS samples are prone to fungal and mycotoxin contamination if improperly harvested, processed, and stored. At present, regulation has been set for the maximum limits of AFB1 (5 μg/kg) and the sum of AFB1, AFB2, AFG1, and AFG2 in ZSS (10 μg/kg) in Chinese Pharmacopeia, respectively [24]. The accurate and rapid differentiation of fungi in ZSS is important to guarantee the safe use of ZSS.
At present, the identification of fungi in herbal medicines is mainly based on fungi isolation and culture. The morphological and microscopic features of pure isolates are generally combined with DNA sequences for a comprehensive analysis [4,28,29]. However, fungi isolation and culture is a complicated and time-consuming procedure with the risk of missing some strains, thereby leading to an imprecise characterization of the fungal diversities. Therefore, an accurate and rapid method to detect fungi in herbal medicines is urgently needed. High-throughput sequencing (HTS) is a growth-independent method that can provide mass data of the composition of mixed microbial communities in low abundances, such as soil, sediment, and air filter samples. HTS has been widely used in fungal ecology studies [30,31,32], thereby providing a new prospect for the detection of fungal diversity in herbal medicines.
Here, we first investigated the presence of fungi in ZSS by using the Illumina MiSeq PE250 platform and demonstrated the diversity of fungal contamination in the ZSS samples to provide early warning for mycotoxin contamination in ZSS.

2. Results

2.1. Analyses of the Diversity of Fungal Communities in the ZSS Samples

All 14 ZSS samples were successfully amplified by PCR for the internal transcribed spacer 2 (ITS2) sequences. A total of 912,941 ITS2 sequences that were longer than 200 bp were produced after excluding the chimeric sequences. The optimized sequences were divided into 210 operational taxonomic units (OTUs) after cluster analysis.
Table S1 shows the calculated Chao 1, Shannon, and Good’s coverage to measure the α-diversity. High Chao 1 value demonstrated a large variation of species in each sample. High Shannon value indicated high community diversity in each sample. Meanwhile, the results of Good’s coverage, which is an estimator of sampling completeness, indicated good overall sampling with levels of >99.8%. Rarefaction curve analysis showed that all samples were almost parallel to the x-axis, thereby indicating that the obtained reads were sufficient to represent the overall fungal diversity (Figure S1). The fungal community in the WM group also had higher Chao 1 estimate and Shannon index than those of the FM group (Figure 1a). The number of unique and common OTUs for the two groups are shown in a Venn diagram (Figure 1b). The results showed that the WM group possessed more unique OTUs than the FM group.
For β-diversity, principal coordinate analysis (PCoA) showed that samples were clustered according to the presence of mold. This result was consistent with the results of hierarchical clustering analysis (Figure 1c–d). Significant differences were reported between the FM and WM groups (analysis of similarity, ANOSIM, R = 0.524 and p = 0.002; Figure S2a). By contrast, the samples from the five producing areas showed an insignificant difference (Figure S2b).

2.2. Fungal Community Composition in the ZSS Samples

Three fungal phyla, namely, Ascomycota, Basidiomycota, and Mucoromycota, were identified in all ZSS samples. Among the three phyla, Ascomycota was the dominant fungus, accounting for 64.36–99.74% of the fungal reads. Meanwhile, other fungal phyla were detected with extremely low relative abundances (0–0.33%, Figure 2a).
Further taxonomical classification detected 61 genera, and the genera with high relative abundance are shown in Figure 2b. Among these genera, Aspergillus, Candida, and Wallemia were the most common, with the relative abundances of 13.52–87.87%, 0.42–64.56%, and 0.06–34.31%, respectively. The differences in the fungal community compositions between the FM and WM groups were significant. The relative abundances of the Aspergillus species in the FM group (42.66–87.87%) were higher than those of the Aspergillus species in the WM group (13.52 to 49.67%). Meanwhile, the relative abundances of the Candida species in the FM group (0.42–23.95%) were relatively lower than those of the Candida species in the WM group (3.13–64.56%). The linear discriminant analysis effect size (LEfSe) algorithm, which was used to identify the different relative abundances of the fungal taxa in the two groups, showed the same results (Figure 2c). Compared with the WM group, Aspergillus was the most dominant genus in the FM group. Meanwhile, the WM group possessed a much higher proportion of Candida, Meyerozyma, and Botryosphaeria species than those in the FM group.
A total of 162 fungal taxa were identified in the ZSS samples. Among these fungal taxa, 70 can be identified at the species level, while the remaining 92 can be resolved at the genus level or higher. Among the 70 accurately differentiated species, potential mycotoxin-producing fungi, namely, A. flavus, A. fumigatus, and Penicillium citrinum that account for 0.30–36.29%, 0.04–7.37%, and 0.01–0.80% of the fungal reads, respectively, were detected in all ZSS samples (Figure 2d). A higher percentage of potential mycotoxin-producing fungi was also detected in the WM group (1.14–36.37%) than that in the FM group (0.61–5.03%).

3. Discussion

3.1. Prevalence of Fungal Contamination in the Commercial ZSS Samples

Fungi are common natural contaminants of herbal medicines due to their widespread distribution in the air. Fungal contamination in medicinal herbs has received considerable concern from the public. An investigation [33] on the occurrence of fungi in medicinal herbs and spices from India showed that 92% samples are contaminated by fungi, and 47% of the contaminated samples exceed the permissible limits set by the World Health Organization. Kong et al. [4] also found that 14 functional foods and 10 spices from Chinese markets are infected by molds. Zheng et al. [34] also showed the widespread fungal contamination in 15 medicinal herbs collected from China. Meanwhile, Su et al. [35] found that 48 samples of eight root herbs from Chinese markets are all contaminated by fungi. In all these studies, Aspergillus and Penicillium are identified as the most common contaminants. In the present study, the results indicated that Aspergillus was the most abundant genus in the ZSS samples, followed by Candida, accounting for 13.52–87.87% and 0.42–64.56% of the total reads, respectively. In contrast to previous studies, Penicillium species were detected at relatively low abundances (<1%) in all ZSS samples, except for the LY01 sample (7.77%). The relationship between fungal species and herbal materials cannot be fully explained owing to the complicated contamination reasons [34,36]. Aspergillus and Penicillium species were the principal storage fungi [36]. The substrate composition and storage conditions, such as moisture content, aeration, and temperature, might be the factors to influence the fungal community [34,36].
Meanwhile, Aspergillus, Penicillium, Fusarium, and Alternaria are the most common contaminants in food and herbal medicines, which contain the majority of mycotoxin producers [3,16]. Potential toxigenic fungi, such as A. fumigatus, A. flavus, A. tubingensis, and A. aculeatus, were previously isolated in different herbal medicines [3,34,37]. In the present study, three potential mycotoxin-producing fungi, namely, A. flavus, A. fumigatus, and P. citrinum that are potential producers of AFB1, OTA, and citrinin, respectively, were identified in the ZSS samples. The high relative abundances of the potential mycotoxin-producing fungi ranging from 0.61–36.37% in ZSS poses a considerable threat to public health. The relative abundances of the potential toxigenic fungi in normal samples were much higher than those in the moldy samples. Therefore, evaluating the safety of herbal medicines according to the presence or absence of molds is unadvisable.
The differences in the fungal community structure between the FM and WM groups were significant. Candida species were dominant in the normal ZSS samples, whereas Aspergillus was the predominant genus in the presence of molds. Meanwhile, all samples were clustered according to the initial grouping in PCoA, thereby indicating the meaningfulness of grouping samples according to the presence of microscopic molds. However, the differences in the fungal community structures among the SX, HN, LN, HB, and SD groups were insignificant, thereby demonstrating that the presence or absence of molds affected the fungal community structure in ZSS samples to a greater extent than the producing areas.

3.2. DNA Marker Selection to Analyze the Fungal Community in ZSS

An effective molecular marker is the precondition for the application of DNA barcoding technique. Schoch et al. [38] evaluated six candidate DNA regions, namely, ITS, LSU, SSU, RPB1, RPB2, and MCM7, for fungal identification and proposed ITS as the primary barcode marker. The length of the ITS region differs substantially among the different fungal genera and species, which distort the community description due to the preferential amplification of shorter sequences, thereby resulting in a biased quantification of the taxon relative abundances in the HTS research [39,40]. ITS1 and ITS2, which were shorter than the complete ITS, were of appropriate length for HTS and have become the most commonly used markers for fungal diversity analyses. Considering the lower phylogenetic richness and fewer OTUs generated by ITS1 region than those of the ITS2 region, ITS2 is recommended for metabarcoding [41]. In the present study, the results showed that the ITS2 region exhibited good identification ability for fungi in the ZSS samples. A total of 162 fungal taxa were detected in the ZSS samples. Among these fungal taxa, 92 cannot be resolved at the species level mainly due to the following reasons. First, some information on species is lacking on the UNITE database. Second, the ITS2 region cannot provide an adequate variation information in distinguishing some sibling species in Aspergillus, Fusarium, Penicillium, and Candida [42,43,44,45].

3.3. Prospects of Applying Amplicon Sequencing for the Analysis of Fungal Diversity in Herbal Materials

The application of HTS platform in analyzing the abundance and richness of the microbial species with low abundances in environmental samples overcomes the isolation limitations of the culture-based approach. Most microorganisms cannot be cultivated using traditional cultivation techniques [46]. Xia et al. [47] identified 55 fungal genera from Chinese Cordyceps through Illumina Miseq sequencing that were not observed using culture-dependent methods. The amplicon sequencing platform is widely used for the analysis of microbial diversity in soil, outdoor air, food, and other environmental samples [32,48,49,50,51]. The application of HTS platform in analyzing the fungal diversity in herbal materials has not been reported to date. In the present study, we utilized the amplicon sequencing platform and targeted the ITS2 sequences to investigate the fungal contamination in herbal materials for the first time. Fungal contamination in the ZSS samples is extremely common. Hence, all samples were infected. HTS platform is a good tool for investigating the occurrence of fungi, especially the potential toxigenic fungi, in herbal materials. The application of amplicon sequencing provides a new approach to analyze the fungal diversity in herbal materials, thereby providing an early warning for mycotoxin contamination in herbal materials.

4. Materials and Methods

4.1. Sampling

A total of 14 commercial ZSS samples were collected from Henan, Shanxi, Hebei, Liaoning, and Shandong, which are the five main ZSS producing provinces in China. The samples were divided into five groups, namely, SX, HN, LN, HB, and SD, by production areas. Among these commercial samples, six samples were affected with mildew due to inappropriate storage. Samples were also divided into two groups, namely, FM and WM, according to the presence or absence of macroscopic molds. The detailed information of the samples is listed in Table 1.

4.2. DNA Extraction

Approximately 3 g ZSS samples were transferred into a 15 mL sterilized centrifuge tube with 10 mL of sterilized water and shaken with a vortex mixer for 20 min. Then, the mixture was filtered through a single layer of sterile gauze, and the microorganisms were collected from the filtrate by centrifugation at 7830 rpm for 15 min (Centrifuge 5430 R, Eppendorf AG, Hamburg, Germany) [52]. The total DNA was extracted using the EZNA® soil DNA kit (Omega Bio-tek., Inc., Norcross, GA, USA) according to the manufacturer’s instructions.

4.3. PCR Amplification and HTS

ITS2 sequences were amplified with the primer pairs ITS3 (5′-GCATCGATGAAGAACGCAGC-3′) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) [53]. PCR was performed with the following conditions: Initial denaturation at 95 °C for 5 min; 35 cycles of denaturation at 95 °C for 45 s, annealing at 58 °C for 50 s, and elongation at 72 °C for 45 s; and a final extension at 72 °C for 10 min. Amplifications were conducted for each sample in triplicate, and the PCR products were pooled to minimize the PCR bias. The PCR products were detected using electrophoresis on a 2% agarose gel and purified using the DNA gel extraction kit (Axygen, Union City, CA, USA). Subsequently, the PCR products were further identified on a 2% agarose gel and quantified using QuantiFluorTM-ST (Promega, Madison, WI, USA). Then, amplicons were sequenced using the Illumina MiSeq PE250 platform (Illumina, San Diego, CA, USA) by AuwiGene Technology Co., Ltd. (Beijing, China). Raw sequences are available in the Sequence Read Archive of the NCBI under the accession numbers SAMN10275058–SAMN10275071.

4.4. Sequence Analysis

Raw FASTQ files were demultiplexed, and their quality was filtered using the Quantitative Insights into Microbial Ecology (QIIME, version 1.8, http://qiime.org) software [54]. The reads were truncated at any site, thereby receiving an average quality score of <20 over a 50 bp sliding window. The primers were exactly matched, thereby allowing two nucleotide mismatches, and the reads containing ambiguous bases were removed. Only sequences that overlap longer than 10 bp were merged according to their overlap sequences. The sequences were clustered into OTUs at 97% sequence similarity by using UPARSE (version 7.1, http://drive5.com/uparse/) [55], and chimeras were removed using USEARCH (version 8.1.1861, Http://Www.Drive5.Com/Usearch/) [56]. OTUs were denominated at the kingdom, phylum, class, order, family, genus, and species levels according to the UNITE database [57]. To estimate the α-diversity, we calculated the three metrics, including Chao 1, Shannon, and Good’s coverage, after the reads were normalized to the minimum reads (35,924 reads) in each sample. β-Diversity was measured using the weighted UniFrac and unweighted UniFrac distance matrices. PCoA was performed based on the weighted UniFrac distance matrices. Samples were also hierarchically clustered based on the unweighted Unifrac distances by using Unweighted Pair Group Method with Arithmetic Mean (UPGMA). LEfSe was applied to identify the differentially abundant taxa between the FM and WM groups [58]. Statistical differences among the sample groups were tested using ANOSIM (available through QIIME). Venn diagram and heat map were drawn using R tools.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/10/12/494/s1, Figure S1: Rarefaction curves of OTUs for the ZSS samples, Figure S2: ANOSIM analysis of the ZSS samples. (a) The differences in the fungal community structures between the FM and WM groups were significant. (b) The differences in the fungal community structures among samples from the five producing areas were insignificant. Table S1: α-Diversity of the fungal community in the ZSS samples.

Author Contributions

Conceptualization, X.P.; Data curation, M.G. and W.J.; Formal analysis, M.G., W.J. and X.P.; Funding acquisition, M.Y. and X.P.; Methodology, M.G. and X.P.; Resources, J.L. and M.Y.; Validation, X.P.; Writing—original draft, M.G.; Writing—review and editing, X.P.

Funding

This work was supported by a grant from CAMS Innovation Fund for Medical Sciences (CIFMS) (No. 2017-I2M-1-013).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moreira, D.D.L.; Teixeira, S.S.; Monteiro, M.H.D.; De-Oliveira, A.C.A.; Paumgartten, F.J. Traditional use and safety of herbal medicines. Rev. Bras. Farmacogn. 2014, 24, 248–257. [Google Scholar] [CrossRef]
  2. Shim, W.B.; Kim, K.; Ofori, J.A.; Chung, Y.C.; Chung, D.H. Occurrence of aflatoxins in herbal medicine distributed in South Korea. J. Food Prot. 2012, 75, 1991–1999. [Google Scholar] [CrossRef] [PubMed]
  3. Ahmad, B.; Ashiq, S.; Hussain, A.; Bashir, S.; Hussain, M. Evaluation of mycotoxins, mycobiota, and toxigenic fungi in selected medicinal plants of Khyber Pakhtunkhwa, Pakistan. Fungal Biol. 2014, 118, 776–784. [Google Scholar] [CrossRef] [PubMed]
  4. Kong, W.; Wei, R.; Logrieco, A.F.; Wei, J.; Wen, J.; Xiao, X.; Yang, M. Occurrence of toxigenic fungi and determination of mycotoxins by HPLC-FLD in functional foods and spices in China markets. Food Chem. 2014, 146, 320–326. [Google Scholar] [CrossRef] [PubMed]
  5. Wei, R.; Qiu, F.; Kong, W.; Wei, J.; Yang, M.; Luo, Z.; Qin, J.; Ma, X. Co-occurrence of aflatoxin B1, B2, G1, G2 and ochrotoxin A in Glycyrrhiza uralensis analyzed by HPLC-MS/MS. Food Control 2013, 32, 216–221. [Google Scholar] [CrossRef]
  6. Shim, W.B.; Ha, K.S.; Kim, M.G.; Kim, J.S.; Chung, D.H. Evaluation of the transfer rate of ochratoxin A to decoctions of herbal medicines. Food Sci. Biotechnol. 2014, 23, 2103–2108. [Google Scholar] [CrossRef]
  7. Nian, Y.; Wang, H.; Ying, G.; Yang, M.; Wang, Z.; Kong, W.; Yang, S. Transfer rates of aflatoxins from herbal medicines to decoctions determined by an optimized high-performance liquid chromatography with fluorescence detection method. J. Pharm. Pharmacol. 2018, 70, 278–288. [Google Scholar] [CrossRef] [PubMed]
  8. Roberts, S.M.; James, R.C.; Williams, P.L. Principles of Toxicology: Environmental and Industrial Applications; John Wiley & Sons: Hoboken, NJ, USA, 2015. [Google Scholar]
  9. Marin, S.; Ramos, A.J.; Cano-Sancho, G.; Sanchis, V. Mycotoxins: Occurrence, toxicology, and exposure assessment. Food Chem. Toxicol. 2013, 60, 218–237. [Google Scholar] [CrossRef] [PubMed]
  10. Prouillac, C.; Koraichi, F.; Videmann, B.; Mazallon, M.; Rodriguez, F.; Baltas, M.; Lecoeur, S. In vitro toxicological effects of estrogenic mycotoxins on human placental cells: Structure activity relationships. Toxicol. Appl. Pharmacol. 2012, 259, 366–375. [Google Scholar] [CrossRef] [PubMed]
  11. Zinedine, A.; Mañes, J. Occurrence and legislation of mycotoxins in food and feed from Morocco. Food Control 2009, 20, 334–344. [Google Scholar] [CrossRef]
  12. Ashiq, S. Natural occurrence of mycotoxins in food and feed: Pakistan perspective. Compr. Rev. Food Sci. Food Saf. 2015, 14, 159–175. [Google Scholar] [CrossRef]
  13. El Darra, N.; Gambacorta, L.; Solfrizzo, M. Multimycotoxins occurrence in spices and herbs commercialized in Lebanon. Food Control 2019, 95, 63–70. [Google Scholar] [CrossRef]
  14. IARC. Some naturally occurring substances: Food items and constituents, heterocyclic aromatic amines and mycotoxins. In IARC Monographs on the Evaluation of Carcinogenic Risks to Humans; International Agency for Research on Cancer: Lyon, France, 1993; Volume 56. [Google Scholar]
  15. Chawanthayatham, S.; Valentine, C.C.; Fedeles, B.I.; Fox, E.J.; Loeb, L.A.; Levine, S.S.; Slocum, S.L.; Wogan, G.N.; Croy, R.G.; Essigmann, J.M. Mutational spectra of aflatoxin B1 in vivo establish biomarkers of exposure for human hepatocellular carcinoma. Proc. Natl. Acad. Sci. USA 2017, 114, E3101–E3109. [Google Scholar] [CrossRef] [PubMed]
  16. Niessen, L. PCR-based diagnosis and quantification of mycotoxin producing fungi. Int. J. Food Microbiol. 2007, 119, 38–46. [Google Scholar] [CrossRef] [PubMed]
  17. Singh, P.; Srivastava, B.; Kumar, A.; Dubey, N.K. Fungal contamination of raw materials of some herbal drugs and recommendation of Cinnamomum camphora oil as herbal fungitoxicant. Microb. Ecol. 2008, 56, 555–560. [Google Scholar] [CrossRef] [PubMed]
  18. De Curtis, F.; De Felice, D.V.; Ianiri, G.; De Cicco, V.; Castoria, R. Environmental factors affect the activity of biocontrol agents against ochratoxigenic Aspergillus carbonarius on wine grape. Int. J. Food Microbiol. 2012, 159, 17–24. [Google Scholar] [CrossRef] [PubMed]
  19. Jalili, M.; Jinap, S. Natural occurrence of aflatoxins and ochratoxin A in commercial dried chili. Food Control 2012, 24, 160–164. [Google Scholar] [CrossRef]
  20. Romagnoli, B.; Menna, V.; Gruppion, N.; Bergamini, C. Afatoxins in spices, aromatic herbs, herb-teas and medicinal plants marketed in Italy. Food Control 2007, 18, 697–701. [Google Scholar] [CrossRef]
  21. Stevic, T.; Pavlovic, S.; Stankovic, S.; Savikin, K. Pathogenic microorganisms of medicinal herbal drugs. Arch. Biol. Sci. 2012, 64, 49–58. [Google Scholar] [CrossRef]
  22. Niessen, L.; Bechtner, J.; Fodil, S.; Taniwaki, M.H.; Vogel, R.F. LAMP-based group specific detection of aflatoxin producers within Aspergillus section Flavi in food raw materials, spices, and dried fruit using neutral red for visible-light signal detection. Int. J. Food Microbiol. 2018, 266, 241–250. [Google Scholar] [CrossRef] [PubMed]
  23. Duarte, S.C.; Pena, A.; Lino, C.M. A review on ochratoxin A occurrence and effects of processing of cereal and cereal derived food products. Food Microbiol. 2010, 27, 187–198. [Google Scholar] [CrossRef] [PubMed]
  24. Chinese Pharmacopeia Commission. Pharmacopoeia of People’s Republic of China (Part 1); Press of Chemical Industry: Beijing, China, 2015; pp. 366–367. [Google Scholar]
  25. Fang, X.S.; Hao, J.F.; Zhou, H.Y.; Zhu, L.X.; Wang, J.H.; Song, F.Q. Pharmacological studies on the sedative-hypnotic effect of Semen Ziziphi spinosae (Suanzaoren) and Radix et Rhizoma Salviae miltiorrhizae (Danshen) extracts and the synergistic effect of their combinations. Phytomedicine 2010, 17, 75–80. [Google Scholar] [CrossRef] [PubMed]
  26. Xie, J.; Guo, L.; Pang, G.; Wu, X.; Zhang, M. Modulation effect of Semen Ziziphi Spinosae extracts on IL-1β, IL-4, IL-6, IL-10, TNF-α and IFN-γ in mouse serum. Nat. Prod. Res. 2011, 25, 464–467. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, Y.; Qiao, L.; Song, M.; Wang, L.; Xie, J.; Feng, H. Hplc-ESI-MS/MS analysis of the water-soluble extract from Ziziphi spinosae semen and its ameliorating effect of learning and memory performance in mice. Pharmacogn. Mag. 2014, 10, 509. [Google Scholar] [PubMed] [Green Version]
  28. Chen, A.J.; Huang, L.F.; Wang, L.Z.; Tang, D.; Cai, F.; Gao, W.W. Occurrence of toxigenic fungi in ochratoxin A contaminated liquorice root. Food Addit Contam A 2011, 28, 1091–1097. [Google Scholar] [CrossRef] [PubMed]
  29. Al-Hindi, R.R.; Aly, S.E.; Hathout, A.S.; Alharbi, M.G.; Al-Masaudi, S.; Al-Jaouni, S.K.; Harakeh, S.M. Isolation and molecular characterization of mycotoxigenic fungi in agarwood. Saudi J. Biol. Sci. 2017. [Google Scholar] [CrossRef]
  30. Tedersoo, L.; Bahram, M.; Põlme, S.; Kõljalg, U.; Yorou, N.S.; Wijesundera, R.; Ruiz, L.; Vasco-Palacios, A.; Thu, P.; Suija, A.; et al. Global diversity and geography of soil fungi. Science 2014, 346, 1256688. [Google Scholar] [CrossRef] [PubMed]
  31. Bittinger, K.; Charlson, E.S.; Loy, E.; Shirley, D.J.; Haas, A.R.; Laughlin, A.; Yi, Y.; Wu, G.; Lewis, J.; Frank, I.; et al. Improved characterization of medically relevant fungi in the human respiratory tract using next-generation sequencing. Genome Biol. 2014, 15, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lee, S.; An, C.; Xu, S.; Yamamoto, N. High-throughput sequencing reveals unprecedented diversities of Aspergillus species in outdoor air. Lett. Appl. Microbiol. 2016, 63, 165–171. [Google Scholar] [CrossRef] [PubMed]
  33. Aiko, V.; Mehta, A. Prevalence of toxigenic fungi in common medicinal herbs and spices in India. 3 Biotech. 2016, 6, 159. [Google Scholar] [CrossRef] [PubMed]
  34. Zheng, R.S.; Wang, W.L.; Tan, J.; Xu, H.; Zhan, R.T.; Chen, W.W. An investigation of fungal contamination on the surface of medicinal herbs in China. Chin. Med. 2017, 12, 2. [Google Scholar] [CrossRef] [PubMed]
  35. Su, C.; Hu, Y.; Gao, D.; Luo, Y.I.; Chen, A.J.; Jiao, X.; Gao, W. Occurrence of Toxigenic Fungi and Mycotoxins on Root Herbs from Chinese Markets. J. Food Prot. 2018, 81, 754–761. [Google Scholar] [CrossRef] [PubMed]
  36. Mandeel, Q.A. Fungal contamination of some imported spices. Mycopathologia 2005, 159, 291–298. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, A.J.; Jiao, X.; Hu, Y.; Lu, X.; Gao, W. Mycobiota and mycotoxins in traditional medicinal seeds from China. Toxins 2015, 7, 3858–3875. [Google Scholar] [CrossRef] [PubMed]
  38. Schoch, C.L.; Seifert, K.A.; Huhndorf, S.; Robert, V.; Spouge, J.L.; Levesque, C.A.; Chen, W.; Fungal Barcoding Consortium. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc. Natl. Acad. Sci. USA 2012, 109, 6241–6246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bellemain, E.; Carlsen, T.; Brochmann, C.; Coissac, E.; Taberlet, P.; Kauserud, H. ITS as an environmental DNA barcode for fungi: An in silico approach reveals potential PCR biases. BMC Microbiol. 2010, 10, 189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. De Filippis, F.; Laiola, M.; Blaiotta, G.; Ercolini, D. Different amplicon targets for sequencing-based studies of fungal diversity. Appl. Environ. Microbiol. 2017, 83, pii: e00905-17. [Google Scholar] [CrossRef]
  41. Tedersoo, L.; Anslan, S.; Bahram, M.; Põlme, S.; Riit, T.; Liiv, I.; Urmas, K.; Veljo, K.; Henrik, N.; Falk, H.; et al. Shotgun metagenomes and multiple primer pair-barcode combinations of amplicons reveal biases in metabarcoding analyses of fungi. MycoKeys 2015, 10, 1–43. [Google Scholar] [CrossRef] [Green Version]
  42. Samson, R.A.; Visagie, C.M.; Houbraken, J.; Hong, S.B.; Hubka, V.; Klaassen, C.H.; Perrone, G.; Seifert, K.A.; Susca, A.; Tanney, J.B.; et al. Phylogeny, identification and nomenclature of the genus Aspergillus. Stud. Mycol. 2014, 78, 141–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bakker, M.G. A fungal mock community control for amplicon sequencing experiments. Mol. Ecol. Resour. 2018, 18, 541–556. [Google Scholar] [CrossRef] [PubMed]
  44. Visagie, C.M.; Houbraken, J.; Frisvad, J.C.; Hong, S.B.; Klaassen, C.H.W.; Perrone, G.; Seifert, K.A.; Varga, J.; Yaguchi, T.; Samson, R.A. Identification and nomenclature of the genus Penicillium. Stud. Mycol. 2014, 78, 343–371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. De Llanos Frutos, R.; Fernández-Espinar, M.T.; Querol, A. Identification of species of the genus Candida by analysis of the 5.8 S rRNA gene and the two ribosomal internal transcribed spacers. Anton. Leeuw. 2004, 85, 175–185. [Google Scholar] [CrossRef] [PubMed]
  46. Daniel, R. The soil metagenome–a rich resource for the discovery of novel natural products. Curr. Opin. Biotechnol. 2004, 15, 199–204. [Google Scholar] [CrossRef] [PubMed]
  47. Xia, F.; Chen, X.; Guo, M.Y.; Bai, X.H.; Liu, Y.; Shen, G.R.; Li, Y.; Lin, J.; Zhou, X.W. High-throughput sequencing-based analysis of endogenetic fungal communities inhabiting the Chinese Cordyceps reveals unexpectedly high fungal diversity. Sci. Rep. 2016, 6, 33437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Fierer, N.; Leff, J.W.; Adams, B.J.; Nielsen, U.N.; Bates, S.T.; Lauber, C.L.; Owens, S.; Gilbert, J.; Wall, D.; Caporaso, J.G. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc. Natl. Acad. Sci. USA. 2012, 109, 21390–21395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Žifčáková, L.; Větrovský, T.; Howe, A.; Baldrian, P. Microbial activity in forest soil reflects the changes in ecosystem properties between summer and winter. Environ. Microbiol. 2016, 18, 288–301. [Google Scholar] [CrossRef] [PubMed]
  50. Ercolini, D. High-throughput sequencing and metagenomics: Moving forward in the culture-independent analysis of food microbial ecology. Appl. Environ. Microbiol. 2013, 79, 3148–3155. [Google Scholar] [CrossRef] [PubMed]
  51. Marsh, A.J.; O’Sullivan, O.; Hill, C.; Ross, R.P.; Cotter, P.D. Sequence-based analysis of the bacterial and fungal compositions of multiple kombucha (tea fungus) samples. Food Microbiol. 2014, 38, 171–178. [Google Scholar] [CrossRef] [PubMed]
  52. Zhao, M.; Wang, B.; Li, F.; Qiu, L.; Li, F.; Wang, S.; Cui, J. Analysis of bacterial communities on aging flue-cured tobacco leaves by 16S rDNA PCR–DGGE technology. Appl. Microbiol. Biotechnol. 2007, 73, 1435–1440. [Google Scholar] [CrossRef] [PubMed]
  53. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols a Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: San Diego, CA, USA, 1990; pp. 315–322. [Google Scholar]
  54. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.; Goodrich, J.; Gordon, J.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Med. 2010, 7, 335–336. [Google Scholar] [Green Version]
  55. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996. [Google Scholar] [CrossRef] [PubMed]
  56. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kõljalg, U.; Nilsson, R.H.; Abarenkov, K.; Tedersoo, L.; Taylor, A.F.; Bahram, M.; Bates, S.; Bruns, T.; Bengtsson-Palme, J.; Callaghan, T.; et al. Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 2013, 22, 5271–5277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed]
Toxins 10 00494 g001 550
Figure 1. Diversity analyses of the fungal community in the Ziziphi Spinosae Semen (ZSS) samples. (a) α-Diversity plots of taxa richness (Chao 1 estimate) and diversity (Shannon index) in the FM and WM groups. (b) Venn diagram of operational taxonomic units (OTUs) in the FM and WM groups. (c) principal coordinate analysis (PCoA) plots on the basis on weighted UniFrac distance matrices. (d) Unweighted Pair Group Method with Arithmetic Mean (UPGMA) clustering on the basis of unweighted UniFrac distance analysis.
Figure 1. Diversity analyses of the fungal community in the Ziziphi Spinosae Semen (ZSS) samples. (a) α-Diversity plots of taxa richness (Chao 1 estimate) and diversity (Shannon index) in the FM and WM groups. (b) Venn diagram of operational taxonomic units (OTUs) in the FM and WM groups. (c) principal coordinate analysis (PCoA) plots on the basis on weighted UniFrac distance matrices. (d) Unweighted Pair Group Method with Arithmetic Mean (UPGMA) clustering on the basis of unweighted UniFrac distance analysis.
Toxins 10 00494 g001
Toxins 10 00494 g002 550
Figure 2. Composition analyses of the fungal communities in the ZSS samples. (a) Fungal composition at the phylum level. (b) Fungal composition at the genus level. (c) Fungal taxa that were differentially abundant in the FM and WM groups visualized using linear discriminant analysis effect size (LEfSe) analysis. (d) Relative abundances of the potential toxigenic fungi detected in the ZSS samples.
Figure 2. Composition analyses of the fungal communities in the ZSS samples. (a) Fungal composition at the phylum level. (b) Fungal composition at the genus level. (c) Fungal taxa that were differentially abundant in the FM and WM groups visualized using linear discriminant analysis effect size (LEfSe) analysis. (d) Relative abundances of the potential toxigenic fungi detected in the ZSS samples.
Toxins 10 00494 g002
Table
Table 1. Voucher information and GenBank accession numbers for the ZSS samples in this study.
Table 1. Voucher information and GenBank accession numbers for the ZSS samples in this study.
NameVoucher No.SourcesGroupMildewyGroupGenBank Accession No.
Ziziphi Spinosae SemenRY01Ruyang, HenanHNNoWMSAMN10275060
Ziziphi Spinosae SemenRY02Ruyang, HenanHNNoWMSAMN10275061
Ziziphi Spinosae SemenHX01Hui, HenanHNNoWMSAMN10275062
Ziziphi Spinosae SemenSJ01Shijiazhuang, HebeiHBNoWMSAMN10275069
Ziziphi Spinosae SemenQH01Qinhuangdao, HebeiHBNoWMSAMN10275067
Ziziphi Spinosae SemenSL01Shengli, LiaoningLNNoWMSAMN10275065
Ziziphi Spinosae SemenSL02Shengli, LiaoningLNNoWMSAMN10275066
Ziziphi Spinosae SemenDP01Dongping, ShandongSDNoWMSAMN10275070
Ziziphi Spinosae SemenLL01Lanling, ShandongSDYesFMSAMN10275071
Ziziphi Spinosae SemenQH02Qinhuangdao, HebeiHBYesFMSAMN10275068
Ziziphi Spinosae SemenYC01Yuncheng, ShanxiSXYesFMSAMN10275058
Ziziphi Spinosae SemenYC02Yuncheng, ShanxiSXYesFMSAMN10275059
Ziziphi Spinosae SemenLY01Lingyuan, LiaoningLNYesFMSAMN10275063
Ziziphi Spinosae SemenLY02Lingyuan, LiaoningLNYesFMSAMN10275064

Share and Cite

MDPI and ACS Style

Guo, M.; Jiang, W.; Luo, J.; Yang, M.; Pang, X. Analysis of the Fungal Community in Ziziphi Spinosae Semen through High-Throughput Sequencing. Toxins 2018, 10, 494. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120494

AMA Style

Guo M, Jiang W, Luo J, Yang M, Pang X. Analysis of the Fungal Community in Ziziphi Spinosae Semen through High-Throughput Sequencing. Toxins. 2018; 10(12):494. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120494

Chicago/Turabian Style

Guo, Mengyue, Wenjun Jiang, Jiaoyang Luo, Meihua Yang, and Xiaohui Pang. 2018. "Analysis of the Fungal Community in Ziziphi Spinosae Semen through High-Throughput Sequencing" Toxins 10, no. 12: 494. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120494

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