Fungal community succession and major components change during manufacturing process of Fu brick tea AbstractFu brick tea is a unique post-fermented tea product which is fermented with microorganism during the manufacturing process. Metabolic analysis showed that most metabolites content were decreased during the manufacturing process of Fu brick tea, except GA (gallic acid). Illumina MiSeq sequencing of ITS gene amplicons was applied to analyze the fungal community succession. The genera Aspergillus, Cyberlindnera and Candida were predominant at the early stage of manufacturing process (from 鈥減rimary dark tea鈥?to 鈥渇ermentation for 3 days鈥?, but after the stage of 鈥渇ermentation for 3 days鈥?only Aspergillus was still dominated, and maintain a relatively constant until to the end of manufacturing process. The effects of metabolites on the structure of the fungal community were analyzed by redundancy analysis (RDA) and variation partitioning analysis (VPA). The results indicated that GCG (gallocatechin gallate), EGCG (epigallocatechin gallate) and GA as well as the interactions among them were the most probably ones to influence, or be influenced by the fungal communities during the fermentation process of Fu brick tea. This study revealed fungal succession, metabolite changes and their relationships, provided new insights into the mechanisms for manufacturing process of Fu brick tea. IntroductionIn general, Chinese teas are divided into six categories according to the processing technology and oxidation degree: green tea (no oxidization), white tea (slightly oxidized), yellow tea (lightly oxidized), oolong tea (partially oxidized), black tea (fully oxidized), dark tea (post-fermented)1. These six types of tea have different flavor quality and health benefits because of their characteristic chemical components. The manufacturing process of Fu brick tea involves steaming, piling, pressing, fermentation (microbial growth) and drying. A complex biochemical changes take place during the manufacturing process producing its special flavor and health benefits of Fu brick tea, such as anti-hyperlipidemia, anti-obesity, anti-hyperglycemia, and anti-dysentery activities2,3,4,5.Since the quality of Fu brick tea is closely related to the microbial activities during the manufacturing process, microbial community have been intensive studied in the past few years. Microbial counting and identification revealed that genera Aspergillus, Penicillium, and Eurotium were the dominating fungus during the manufacturing process of Fu brick tea6, 7. However, traditional microbiological methods, such as cell cultures and colony counting, are of limited value for exploring the variation and structure of microbial communities. Moreover, isolation media may be suitable for only some types of microorganisms, as strains of microorganisms cannot be accurately discriminated based on appearance8. Several molecular biology techniques, such as terminal restriction fragment length polymorphism (T-RFLP), polymerase chain reaction denaturing gradient gel electrophoresis (PCR-DGGE), and gene library construction, were used for determining microbial communities and structures9. A recent studies using PCR-DGGE analysis revealed that the microorganisms found in Fu brick tea were from or closely related to the genera Aspergillus, Beauveria, Debaryomyces, Eurotium, Pestalotiopsis, Pichia, Rhizomucor, and Verticillium 10, while other studies found that Aspergillus niger, Blastobotrys adeninivorans and Bacillus, Enterobacteriaceae were the major fungal and bacterial communities involved in Pu-erh tea, another kind of post-fermentation tea11, 12. However, these techniques provide limited information on community information, because only a few sequences can be separated and analyzed. Recently, a high-throughput sequencing technique, also called next generation sequencing technology, was used to investigate microbial composition and diversity in various fermented foods, such as doubanjiang-meju, wine grape and Pu-erh tea13,14,15. Compared to previous culture-dependent or low resolution molecular method (e.g. PCR-DEGG), this technology produces a very large number of reads in a run of many different samples and provides a powerful tool for profiling microbial structure and diversity at genus or even species levels16, 17.The aim of this study was to investigate the fungal community succession, metabolite changes, and assess relationships between fungal community and metabolites during the manufacturing process of Fu brick tea. This knowledge could give new insight into the manufacturing process and provide valuable knowledge to improve the quality of Fu brick tea.ResultsMetabolite changes during the manufacturing process of Fu brick teaThe metabolites, including water extract (WE), soluble sugar (SS), flavonoid (FLA), organic acid (OA), tea polyphenol (TP), amino acid (AA), gallic acid (GA), caffeine (CAF), epigallocatechin gallate (EGCG), epicatechin gallate (ECG), gallocatechin gallate (GCG), epigallocatechin (EGC), epicatechin (EC), catechins (C), were analyzed during the manufacturing process of Fu brick tea. As shown in Fig.聽1, significant decreasing of WE and SS were beginning at S5 and S6, then continued decreased to the end of manufacturing process of Fu brick tea (P鈥?lt;鈥?.05). A 16.68% decrease in WE and a 26.87% decrease in SS were observed at S10 compared with S1. The level of FLA and OA ranged from 0.86鈥壜扁€?.01% to 1.05鈥壜扁€?.03% and 1.34鈥壜扁€?.20% to 1.65鈥壜扁€?.02%, but no significant difference were observed during the manufacturing process of Fu brick tea (Fig.聽1A). The concentration of TP was significantly decreased by 30.61% from S1 to S10 during the manufacturing process of Fu brick tea (P鈥?lt;鈥?.05). But AA ranged from 1.10鈥壜扁€?.01% to 1.35鈥壜扁€?.02%, which slightly decreased but no significantly difference were observed during the process of Fu brick tea (P鈥?gt;鈥?.05) (Fig.聽1B). The level of CAF ranged from 0.91鈥壜扁€?.01% to 1.17鈥壜扁€?.01%, which was slightly decreased beginning at S2 and continued to the end of the process (P鈥?gt;鈥?.05). The level of GA was relatively stable from S1 to S5 and sharply increased from S5 to S7, then decreased until the end of the process, but the concentration also increased by 58.48% in S10, compared with S1 (P鈥?lt;鈥?.05) (Fig.聽1C). The level of EGC was slightly decreased from S1 to S8 and dramatically decreased from S9 to S10 and a 58.83% decrease was observed at S10 compared with S1 (P鈥?lt;鈥?.05). But the concentration of C and EC ranged from 0.25鈥壜扁€?.01% to 0.34鈥壜扁€?.01% and 0.16鈥壜扁€?.01% to 0.24鈥壜扁€?.01%, which were relatively stable during the manufacturing process. The level of EGCG, GCG and ECG were slightly decreased from S1 to S7 (P鈥?gt;鈥?.05), and significantly decreased from S8 to the end of the manufacturing process of Fu brick tea (P鈥?lt;鈥?.05) (Fig.聽1D).Figure 1Metabolite changes during manufacturing process of Fu brick tea. water extract (WE), soluble sugar (SS), flavonoid (FLA), organic acid (OA), tea polyphenol (TP), amino acid (AA), gallic acid (GA), caffeine (CAF), epigallocatechin gallate (EGCG), epigallocatechin (EGC), epicatechin gallate (ECG), gallocatechin gallate (GCG), epicatechin (EC), catechins (C). *Indicated that the significance 0.05, compared with S1; **indicated that the significance 0.01, compared with S1.Full size image Diversity of fungal communities during the manufacturing process of Fu brick teaAfter quality filtering and chimera removal, a total of 285,184 of valid sequences were generated from the samples, with an average sequence number of 28,518 for each sample (range from 20846 to 37278). At 97% sequence identity, 216 OTUs (operational taxonomic units) were identified (Supplementary Information Table聽S1), ranged from 3 to 130 through the manufacturing process. The value of Good鈥檚 coverage were 99% for all sequences from ten groups, indicating sufficient sequencing depth for microbial communities in Fu brick tea. The Shannon and Simpson index, the Chao1 and ACE estimator demonstrated that both the diversity and richness of fungus were no significant change from S1 to S4 sample, but dramatically decreased in S5 and maintained relatively constant until the end of the manufacturing process of Fu brick tea (Supplementary Information Table聽S2). These changes were related to the dominance of specific taxonomic group during the manufacturing process. The rarefaction curves and Shannon-Wiener curves for all samples almost reached the saturation phase, suggesting that there few new microbes would be identified by increasing the sequence depth, and the majority of fungal microbes in the samples had already been captured in the analysis (Supplementary Figure聽S1A and B). Rank-abundance curves indicated that a majority of the reads belonged to rare organisms represented by only a few sequences, particularly for S5~S10 samples (Supplementary Figure聽S2). To evaluate the similarities of fungal communities in different process stage of Fu brick tea samples, NMDS analysis based on Bray-Curtis distance was performed. The result showed that the samples could be divided into four groups: the first group, including S1; the second group, including S2; the third group, including S3, S4; the fourth group, including S5~S10 (Fig.聽2A). The Venn diagram analysis was used to show the shared OTUs among different stages of manufacturing process after singleton sequences being removed. The samples were divided into four groups based on the NMDS analysis. This analysis showed that only 5 of the identified OTUs were shared among all the groups, while 38 of 174 OTUs, 54 of 229 OTUs and 6 of 133 OTUs were shared between the groups in group I and group II, group II and group III, group III and group IV, respectively (Fig.聽2B). The detail information of Venn diagram analysis was provided in Supplementary Information Table聽S3.Figure 2NMDS (A) and Venn diagram analysis (B) of fungal communities in samples during manufacturing process of Fu brick tea.Full size image Composition of fungal communities during the process of Fu brick teaTo identify the taxonomic composition of fungal communities among three groups, the RDP classifier was used to assign the sequence tags to different taxonomic levels at 70% threshold. All the identified OTUs could be assigned into 4 different phyla, 16 classes, 40 orders, 69 families or 102 genera. As shown in Fig.聽3, Ascomycota was the predominant phylum in all samples, accounting for 95.05%~99.99% of total effective sequences. Basidiomycota and some other phylum fungi were also identified in the early stage of process samples, but almost not detected after S5 of manufacturing process (Fig.聽3A). The order Eurotiales, Saccharomycetales and Hypocreales predominant from S1 to S4 samples, but quickly replaced by Eurotiales after S5 and remained stable until the end of manufacturing process in Fu brick tea (Fig.聽3B).Figure 3Fungal taxonomic compositions showing the fungal successions at phylum (A), order (B) and genus (C) level during manufacturing process of Fu brick tea. The taxonomic abundance 1% were classified into 鈥渙thers鈥? Phylum level (A), order level (B), genus level (C).Full size image At the genus level, Aspergillus, Cyberlindnera and Candida comprised 80% of all sequence. Aspergillus dominated the whole manufacturing process of Fu brick tea, which represented the largest fraction, first increased from 65.13% in S1 sample to 93.58% in S2 sample and significantly decreased to 43.63% in S3, then dominant until the end of process, with relative abundance clearly increased to 99.99% in S6. Cyberlindnera account for 23.46% of all sequences in S1 and was dramatically decreased to 1.47% in S2, then sharply increased to 25.72% in S3, and almost not detected from S6 sample to the end of manufacturing process. Candida increased from the initial 3.25% in S1 to 27.46% in S2 and decreased to 13.23% in S3, but also not detected from S6 to the end of manufacturing process. Variations in some small proportions in the groups Wallemia (0.81% to 0.57%), Penicillium (2.49% to 0.01%) were also observed during manufacturing process of Fu brick tea (Fig.聽3C).Relationships between fungal community and metabolitesIn order to evaluate the effects of metabolite factors on the structure of fungal community during manufacturing process of Fu brick tea, the multivariate analysis were performed. The gradient lengths were 0.959, 0.652, 0.186 and 0.183 from the first to the fourth axis, as retrieved by detrended correspondence analysis (DCA). The result of DCA indicated that the gradient length of the first axis was less than 3SD. Therefore, the linear model with RDA was considered to be the appropriate ordination method for direct gradient analysis. As shown in Fig.聽4, the results of forward selection of RDA showed that the most variation of the fungal community structure could be well explained by GCG, EGCG and GA during the manufacturing process. The model statistically explained up to 90.30% of the variation (P鈥?鈥?.004) (Fig.聽4A). Furthermore, the result of VPA indicated that GCG, EGCG and GA explained 38.10% (P鈥?鈥?.002), 23.80% (P鈥?鈥?.004) and 7.10% (P鈥?鈥?.026) of the variation in fugal community structure, respectively. The variation shared by GCG, EGCG and GA catechins was 21.30% (Fig.聽4B). Those factors as well as the interactions among them were the most probably ones to influence, or be influenced by the fungal communities during the fermentation process of Fu brick tea.Figure 4RDA and VPA showing the correlations between fungal structure and metabolites concentrations during manufacturing process of Fu brick tea. Ordination diagram of RDA (A). Variance decomposition of VPA (B). Gallic acid (GA), epigallocatechin gallate (EGCG), gallocatechin gallate (GCG). Red arrow indicated the metabolite variance significant effect on the fungal community structure (P鈥?lt;鈥?.05).Full size image As shown in Table聽1, Pearson鈥檚 correlation between fungal community and major metabolites were also analyzed. Genus of Aspergillus was significantly (P鈥?lt;鈥?.05) correlated with WE, TP, AA, CAF, EGCG, GCG and ECG. Cyberlindnera was significantly (P鈥?lt;鈥?.05) correlated with WE, TP, CAF, GCG and ECG. Candida was significantly (P鈥?lt;鈥?.05) correlated with AA, GCG and ECG. Wallemia and Fungi_unclassified were significantly (P鈥?lt;鈥?.05) correlated with WE, AA, GCG and ECG. Penicillium was significantly (P鈥?lt;鈥?.05) correlated with WE, TP, CAF, EGC, EGCG, GCG and ECG. Ascomycota_unclassified was significantly (P鈥?lt;鈥?.05) correlated with TP and CAF. Eucasphaeria was significantly (P鈥?lt;鈥?.05) correlated with TP and CAF. However, no fungal genera were significantly (P鈥?gt;鈥?.05) correlated with FLA, OA, EC and C.Table 1 Pearson鈥檚 correlation between the relative abundances of important genera (relative abundance 1%) and metabolite variablesa.Full size table DiscussionFu brick tea is a unique post-fermented tea product which is fermented with microorganism during the manufacturing process. To the best of our knowledge, this is the first report on the fungal community succession and metabolite changes during the manufacturing process of Fu brick tea by high throughput Illumina MiSeq sequencing, spectrophotometry and HPLC. The results, from dramatic changes in fungal community structure and metabolites as well as correlation analysis, indicated strong relationships between metabolites and fungal community during the manufacturing process of Fu brick tea.Our metabolic analysis showed that WE, SS, TP and catehins were significantly decreased, but GA was significantly increased during the manufacturing process of Fu brick tea. These changes probably resulted from moisture and heat reaction and microbial metabolism18. Several enzymes secreted from genus of Aspergillus, such as cellulases, hemicellulases, proteases and 伪-amylases, could catalyze the major metabolites in Pu-erh tea during fermentation process11, 19. The decrease of SS may be related to microbial growth that SS could be used as the carbon source of microbial growth during the manufacturing process20. Our previous study indicated that the degradation of TP and catechins in the early stage of process might be resulted from moisture and heat reaction, but in the later stage of process was mainly related to the microbial metabolism. And this degradation was good for the formation of Fu brick tea鈥檚 taste and its special aroma20. The aldehyde compounds of stale aroma and terpene alcohols of flower aroma which specially increased after the fermentation process were from decarboxylation and oxidative deamination of amino acids21. The increase of GA may be also the result of hydrolysis of galloylated catechins and microbial metabolism, which might be helpful on some bioactivity of Fu brick tea10. Interesting, these decreased components were the major functional components of tea, but Fu brick tea still has a strong health benefits. Mo et al. reported that the activities increased with the course of microbial fermentation process (MFP) of Fu brick tea, which indicated that this process produced new varieties of functional metabolites2. Furthermore, new metabolites such as three triterpenoids22, one norisoprenoid23, seven new B ring fission catechins (flavan-3-ols) derivatives24,25,26 were reported from a few phytochemical investigations on Fu brick tea. Especially these B ring fission catechins have not been identified in other tea products, which suggested that the biochemical profiles of Fu brick tea were influenced significantly by MFP.Fungal community succession showed that the genera Aspergillus, Cyberlindnera and Candida were predominant from S0 to S4, but only the genus of Aspergillus was still dominated in Fu brick tea from S5 and maintained constant until the end of manufacturing process. The genus Aspergillus is a group of filamentous fungi consist of more than 250 species, which is the most economically important of the fungal genera27. Many species of Aspergillus are used in biotechnology for the production of various metabolites, such as antibiotics, organic acids, medicines or enzymes, or as agents in many food fermentations28. The fungal genus Eurotium, which is the teleomorph of Aspergillus, has been proved to be a rich source of novel bioactive metabolites29, 30. Since the recent synonymization of the teleomorph-based genus Eurotium with Aspergillus by the International Commission on Aspergillus (ICPA, 2012), which adopted the newly established principle 鈥渙ne fungus, one name鈥?(Norvell, 2011), species formerly included in the genus Eurotium are displayed with their Aspergillus name27. Many previous reports showed that Aspergillus cristatum (Eurotium cristatum) was the dominant fungus during the manufacturing process of Fu brick tea31, 32, and it is considered to be safe under low- and high-osmolarity conditions by genomic, transcriptomic and HPLC analysis33. Several indole alkaloids and indole diketopiperazine alkaloid were identified from the culture extract of Aspergillus cristatum and these compounds were proved to possess brine shrimp lethality, antibacterial activity against E. coli, radical scavenging activity against DPPH radicals, and exhibited marginal attenuation of 3T3L1 pre-adipocytes34,35,36. Both emodin and physcion were also identified in Fu brick tea fermented by Aspergillus cristatum 10. Non-mycotoxigenic strains of Aspergillus repens and Aspergillus rubrum are used as starter cultures in the manufacture of the traditional fermented food katsuobushi, from bonito (Katsuwonuspelamis), in Japan37. During katsuobushi fermentation these xerophilic fungi produce effective antioxidants, which participate in the suppression of lipid oxidation such that the fermented katsuobushi is resistant to lipid oxidation and has a long shelf-life38.The genus Lindnera was established to accommodate some Pichia species in the Pichia anomala clade39,40,41. But Lindnera is a later homonym of a valid published plant genus, the genus Cyberlindnera was introduced as a replacement name for Lindnera, with 21 new combinations42. The Cyberlindnera clade contains 23 recognized teleomorphic species and 12 Candida species are also related to this clade43. The yeast Cyberlindner ajadinii, previously referred to as Pichia jadinii, Hansenul ajadinii or Candida guilliermondii, was reported to be the teleomorph of Candida utilis (also referred to as Torula yeast)44, which was used for almost one hundred years for various biotechnological applications45. Because of its capability to efficiently assimilate pentose, including xylose and arabinose, and its robust fermentation characteristics, Cyberlindner ajadinii was exploited to grow large amounts on waste hardwood hydrolysates from the pulp industry. Moreover, some species of Cyberlindnera and Candida, which had a high capability of tannin tolerance, were also isolated from Miang (a fermented food product prepared from tea leaves)46. The genus Candida has become one of the largest in species number, present in almost every environment. Yeasts of this genus are abundantly distributed in nature on land and sea, associated with animals or plants and inanimate objects47. This genus is distributed across the ascomycetous yeast domain, overlapping with other genera, according to phylogenetic analysis using ribosomal genes. Some species of Candida genus were used in food fermentation, such as Candida etchellsii, Candida milleri, Candida rugosa and Candida tropicalis 28. Some species of Candida has high phenol degradation ability48,49,50,51, bioconversion of corn fiber, sugarcane bagasse and xylose into the sweetener xylitol52, 53, enhance of clove and smoked favors and biotransformation of maltose into ethanol in soy sauce process54. These properties of yeast Cyberlindnera and Candida genera indicate that their enzyme activities improve the quality of Fu brick tea by degradation of phenol, producing the sweet substance xylitol, and other favors flavour substance. The yeast genera may play an indispensable role in the initial stage of Fu brick tea fermentation.In conclusion, predominant fungus and major metabolites were analyzed using Illumina MiSeq sequencing, spectrophotometry and HPLC methods, respectively. The metabolic analysis showed that most metabolites content were decreased during the manufacturing process of Fu brick tea, except GA. Illumina MiSeq sequencing analysis showed that the genera Aspergillus, Cyberlindnera and Candida were dominant at the early stage of manufacturing process, but after the stage of 鈥渇ermentation for 3 days鈥?only the genus Aspergillus was still dominated and maintain a relatively constant until to the end of manufacturing process of Fu brick tea. In addition, the effects of metabolites on the structure of the fungal community were analyzed by RDA and VPA. The results indicated that GCG, EGCG and GA as well as the interactions among them were the most probably ones to influence, or be influenced by the fungal communities during the fermentation process of Fu brick tea.Materials and MethodsTea leaf samples and process characterizationSampling of Fu brick-tea was carried out in a major Fu brick tea production factory, Yiyang Fu Cha Industry Development Co. Ltd. (Yiyang, Hunan Province, China). Primary dark tea, used as raw materials for the process of Fu brick tea, was manufactured from fresh tea leaves as described by DB43/T660-2011 (Supplement Material聽1)55. And the manufacturing process of Fu brick tea was descripted as follow: a 100鈥塳g primary dark tea was moistened by steaming then piled up to 1.5鈥塵 in height and subjected to temperatures as high as 80鈥壜癈 overnight. After pile fermentation, the tea materials were mixed with 15鈥塋 of tap water, then partitioned and pressed into desired sizes of brick tea. And the brick teas were placed in the fermentation room for 20~22 days. During the fermentation process, the temperature and relative humidity of fermentation room was 30鈥壜扁€?鈥壜癈 and 75鈥壜扁€?% from 0 day to 15 day. Then the temperature of fermentation room was increased by 2鈥壜癈 per day until the end of the manufacturing process during the drying process. Finally, the tea products were packaged and stored. Three independent batches of Fu brick tea were prepared simultaneously by the same manufacturing process. The samples were collected from three independent batches of Fu brick tea at different stage of manufacturing process: primary dark tea (S1), piling-fermentation tea materials (S2), fermenting brick-tea in the fermentation room for 0 days (S3), 3 days (S4), 6 days (S5), 9 days (S6), 12 days (S7), 15 days (S8), 18 days (S9) and newly made tea products (S10) at day 22. The samples were packaged in sterile polyethylene bags, transported to the laboratory and stored at 鈭?0鈥壜癈 until required.Chemical analysisThe content of WE and OA was analyzed determined as described by GB/T 8305-2013 and GB/T 12456-2008, respectively56, 57. The concentration of TP, and AA in tea leaves were determined using the spectrophotometric method based on Folin-Phenol and the ninhydrin assay as described by GB/T 8313-2008 and GB/T 8314-2013, respectively58, 59. The content of GA, CAF, EGCG, ECG, GCG, EGC, CEC and C were determined, as described previously60. The concentration of FLA and SS were determined by a colorimetric method, as described previously61, 62. All analysis was performed in three times. The detailed approaches are described in Supplementary Material聽2.DNA extraction and Illumina MiSeq sequencingEach tea sample (25鈥塯) was mixed with sterile water (125鈥塵L), stirred thoroughly, filtered through three layers of coarse sterile gauze to remove large particles, and centrifuged at18,000鈥塯 for 10鈥塵in at 4鈥壜癈. The microbial genomic DNA was extracted by the E.Z.N.A.Soil DNA Kit (Omega Bio-tek, Norcross, GA, US), and the procedure was according to the manufacturer鈥檚 protocols.The ITS1 region of fungal was amplified by PCR for MiSeq sequencing. The ITS1 region of fungi was amplified with the forward primer (5鈥?CTTGGTCATTTAGAGGAAGTAA-3鈥? and the reverse primer ITS2 (5鈥?GCTGCGTTCTTCATCGATGC-3鈥?63. Different barcode sequences were added at the 5鈥?end of the forward primer for multiplexed pyrosequencing. PCR were carried out in a 20鈥壩糒 reaction volumes containing 10 ng DNA template, 2鈥壩糒 dNTPs (2.5鈥塵M), 0.4鈥壩糒 of each primer (5 渭螠) and0.4 U渭L FastPfu Polymerase (Applied Biosystems) in the appropriate 56 FastPfu Buffer (4鈥壩糒) and de-ionized ultrapure water (to 20鈥壩糒). The protocol was optimized with low cycles for better accuracy and reliability of the subsequent data analysis. The PCR condition were initial denaturation at 95鈥壜癈 for 2鈥塵in, followed by 25 cycles of denaturation at 95鈥壜癈 for 30鈥塻, annealing at 55鈥壜癈 for 30 sand extension at 72鈥壜癈 for 45鈥塻, with a final extension phase at 72鈥壜癈 for 10鈥塵in. The PCR products were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, US) according to the protocol and quantified using QuantiFluor鈩?ST (Promega, U.S.). Samples were then pooled at equal concentrations. Parallel tagged sequencing was performed using Illumina Miseq Sequencing in Majorbio Bio-Pharm Technology Co., Ltd., Shanghai, China.Data processing and statistical analysesRaw reads were demultiplexed and quality-filtered using Trimmomatic software (version 0.36, http://www.usadellab.org/cms/?page=trimmomatic) with the following criteria64: (a) The 250鈥塨p reads were truncated at any site receiving an average quality score less than 20 over a 10鈥塨p sliding window, discarding the truncated reads that were shorter than 50鈥塨p; (b) Exact barcode matching, two nucleotide mismatch in primer matching, reads containing ambiguous characters were removed; (c) Only sequences that overlap longer than 10鈥塨p were assembled according to their overlap region. Reads that could not be assembled were also discarded. These sequences were clustered to OTUs (operational taxonomic units) at 97% sequence identity by using Mothur (version 1.25.0, http://www.mothur.org)65. The OTUs were used for alpha diversity (Shannon and Simpson), richness (ACE and Chao1), Good鈥檚 coverage, Rarefaction curve, Shannon-Wiener curve and Venn diagram were performed by Mothur65. In addition, Nonmetric Multidimensional Scaling (NMDS) diagram and rank_abundance curves were generated by using R package vegan (version 3.4.0, https://mirrors.tuna.tsinghua.edu.cn/CRAN/). CANOCO 4.5 software was used to conduct multivariate analysis using detrended correspondence analysis (DCA), redundancy analysis (RDA) and variation partitioning analysis (VPA)66. One-way ANOVA test, Pearson鈥檚 correlation coefficients and P values were calculated with IBM SPSS 19.0 (https://www.ibm.com/analytics/us/en/technology/spss/). All results are presented as the mean value (卤SE). Differences between groups were declared significant at P鈥?lt;鈥?.05. References1.Wu, Q. J. et al. Discrimination of Chinese teas with different fermentation degrees by stepwise linear discriminant analysis (S-LDA) of the chemical compounds. Journal of Agricultural and Food Chemistry 62, 9336鈥?344 (2014).CAS聽 Article聽 PubMed聽Google Scholar聽 2.Mo, H., Zhu, Y. Chen, Z. Microbial fermented tea鈥揳 potential source of natural food preservatives. Trends in Food Science Technology 19, 124鈥?30 (2008).CAS聽 Article聽Google Scholar聽 3.Fu, D. et al. Fermented Camellia sinensis, Fu Zhuan Tea, regulates hyperlipidemia and transcription factors involved in lipid catabolism. Food research international 44, 2999鈥?005 (2011).CAS聽 Article聽Google Scholar聽 4.Li, Q. et al. Anti-obesity and hypolipidemic effects of Fuzhuan brick tea water extract in high-fat diet-induced obese rats. Journal of the Science of Food and Agriculture 93, 1310鈥?316 (2013).ADS聽 CAS聽 Article聽 PubMed聽Google Scholar聽 5.Zhang, L., Zhang, Z., Zhou, Y., Ling, T. Wan, X. Chinese dark teas: Postfermentation, chemistry and biological activities. Food Research International 53, 600鈥?07 (2013).CAS聽 Article聽Google Scholar聽 6.Wen, Q. Liu, S. Evolutionary regulation of dominant fungi in Fuzhuan brick tea during the fungus growing process. Journal of tea science 11, 56鈥?2 (1991). Google Scholar聽 7.Mo, H., Zhang, H., Li, Y. Zhu, Y. Antimicrobial activity of the indigenously microbial fermented Fuzhuan brick-tea. Journal of Biotechnology 136, S722鈥揝722 (2008).Article聽Google Scholar聽 8.Amann, R. Ludwig, W. Ribosomal RNA-targeted nucleic acid probes for studies in microbial ecology. FEMS Microbiology Reviews 24, 555鈥?65 (2000).CAS聽 Article聽 PubMed聽Google Scholar聽 9.Siqueira, J. F. Jr., Sakamoto, M. Rosado, A. S. Microbial community profiling using terminal restriction fragment length polymorphism (T-RFLP) and denaturing gradient gel electrophoresis (DGGE). Methods in molecular biology (Clifton, N.J.) 666, 71鈥?5 (2010).CAS聽 Article聽Google Scholar聽 10.Xu, A. et al. Fungal community associated with fermentation and storage of Fuzhuan brick-tea. International Journal of Food Microbiology 146, 14鈥?2 (2011).CAS聽 Article聽 PubMed聽Google Scholar聽 11.Abe, M. et al. Characteristic fungi observed in the fermentation process for Puer tea. International Journal of Food Microbiology 124, 199鈥?03 (2008).CAS聽 Article聽 PubMed聽Google Scholar聽 12.Zhao, M. et al. Structure and dynamics of the bacterial communities in fermentation of the traditional Chinese post-fermented pu-erh tea revealed by 16S rRNA gene clone library. World journal of microbiology biotechnology 29, 1877鈥?884 (2013).CAS聽 Article聽Google Scholar聽 13.Li, Z. et al. Bacterial community succession and metabolite changes during doubanjiang-meju fermentation, a Chinese traditional fermented broad bean (Vicia faba L.) paste. Food Chemistry 218, 534鈥?42 (2017).CAS聽 Article聽 PubMed聽Google Scholar聽 14.Bokulich, N. A., Thorngate, J. H., Richardson, P. M. Mills, D. A. Microbial biogeography of wine grapes is conditioned by cultivar, vintage, and climate. Proceedings of the National Academy of Sciences of the United States of America 111, E139鈥揈148 (2014).ADS聽 CAS聽 Article聽 PubMed聽Google Scholar聽 15.Zhao, M. et al. An Integrated Metagenomics/Metaproteomics Investigation of the Microbial Communities and Enzymes in Solid-state Fermentation of Pu-erh tea. Scientific Reports 5, 10117 (2015).ADS聽 CAS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 16.Miller, S. R., Strong, A. L., Jones, K. L. Ungerer, M. C. Bar-coded pyrosequencing reveals shared bacterial community properties along the temperature gradients of two alkaline hot springs in Yellowstone National Park. Applied and environmental microbiology 75, 4565鈥?572 (2009).CAS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 17.Delgado, S. et al. Diversity of thermophilic bacteria in raw, pasteurized and selectively-cultured milk, as assessed by culturing, PCR-DGGE and pyrosequencing. Food Microbiology 36, 103鈥?11 (2013).CAS聽 Article聽 PubMed聽Google Scholar聽 18.Wang, Z., Zhang, Y., Tong, X. Liu, Z. Dynamics of polyphenils and carbohydrates during the primary processing of dark green tea. Journal of tea science 11, 23鈥?8 (1991).CAS聽Google Scholar聽 19.Ward, O. P., Qin, W. M., Dhanjoon, J., Ye, J. Singh, A. In Advances in Applied Microbiology Vol. 58 (eds Joan W. Bennett Geoffrey M. Gadd Allen I. Laskin Sariaslani Sima) 1鈥?75 (Academic Press, 2005).20.Fu, D. et al. Variations of components of Fuzhuan tea during processing. Food Science 29, 64鈥?7 (2008).CAS聽Google Scholar聽 21.Xu, X. Q., Mo, H. Z., Yan, M. C. Zhu, Y. Analysis of characteristic aroma of fungal fermented Fuzhuan brick-tea by gas chromatography/mass spectrophotometry. Journal of the Science of Food and Agriculture 87, 1502鈥?504 (2007).CAS聽 Article聽Google Scholar聽 22.Ling, T. et al. New Triterpenoids and Other Constituents from a Special Microbial-Fermented Tea-Fuzhuan Brick Tea. Journal of Agricultural and Food Chemistry 58, 4945鈥?950 (2010).CAS聽 Article聽 PubMed聽Google Scholar聽 23.Luo, Z. et al. A New Norisoprenoid and Other Compounds from Fuzhuan Brick Tea. Molecules 17, 3539 (2012).CAS聽 Article聽 PubMed聽Google Scholar聽 24.Jiang, H. et al. A new catechin oxidation product and polymeric polyphenols of post-fermented tea. Food Chemistry 129, 830鈥?36 (2011).CAS聽 Article聽 PubMed聽Google Scholar聽 25.Luo, Z. et al. Fuzhuanins A and B: The B-ring Fission Lactones of Flavan-3-ols from Fuzhuan Brick-Tea. Journal of Agricultural and Food Chemistry 61, 6982鈥?990 (2013).CAS聽 Article聽 PubMed聽Google Scholar聽 26.Zhu, Y. et al. Changes of major tea polyphenols and production of four new B-ring fission metabolites of catechins from post-fermented Jing-Wei Fu brick tea. Food Chemistry 170, 110鈥?17 (2015).CAS聽 Article聽 PubMed聽Google Scholar聽 27.Samson, R. A. et al. Phylogeny, identification and nomenclature of the genus Aspergillus. Studies in Mycology 78, 141鈥?73 (2014).CAS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 28.Bourdichon, F. et al. Food fermentations: microorganisms with technological beneficial use. International Journal of Food Microbiology 154, 87鈥?7 (2012).CAS聽 Article聽 PubMed聽Google Scholar聽 29.Ishikawa, Y., Morimoto, K. Hamasaki, T. Flavoglaucin, a metabolite of Eurotium chevalieri, its antioxidation and synergism with tocopherol. Journal of the American Oil Chemists鈥?Society 61, 1864鈥?868 (1984).CAS聽 Article聽Google Scholar聽 30.Slack, G. J., Puniani, E., Frisvad, J. C., Samson, R. A. Miller, J. D. Secondary metabolites from Eurotium species, Aspergillus calidoustus and A. insuetus common in Canadian homes with a review of their chemistry and biological activities. Mycological Research 113, 480鈥?90 (2009).CAS聽 Article聽 PubMed聽Google Scholar聽 31.Qi, Z. Sun, Z. Identification of predominant species in brick tea. Acta Mycologica Sinica 9, 176鈥?79 (1990). Google Scholar聽 32.Mo, H., Zhu, Y. Chen, Z. Microbial fermented tea 鈥?a potential source of natural food preservatives. Trends in Food Science Technology 19, 124鈥?30 (2008).CAS聽 Article聽Google Scholar聽 33.Ge, Y. et al. Comparative genomic and transcriptomic analyses of the Fuzhuan brick tea-fermentation fungus Aspergillus cristatus. BMC Genomics 17, 428 (2016).Article聽 PubMed聽 PubMed Central聽Google Scholar聽 34.Du, F. Y., Li, X. M., Li, C. S., Shang, Z. Wang, B. G. Cristatumins A-D, new indole alkaloids from the marine-derived endophytic fungus Eurotium cristatum EN-220. Bioorganic medicinal chemistry letters 22, 4650鈥?653 (2012).CAS聽 Article聽Google Scholar聽 35.Zou, X. et al. A new prenylated indole diketopiperazine alkaloid from Eurotium cristatum. Molecules 19, 17839鈥?7847 (2014).Article聽 PubMed聽Google Scholar聽 36.Du, F. Y., Li, X., Li, X. M., Zhu, L. W. Wang, B. G. Indolediketopiperazine Alkaloids from Eurotium cristatum EN-220, an Endophytic Fungus Isolated from the Marine Alga Sargassum thunbergii. Marine drugs 15, 1鈥?0 (2017).Article聽Google Scholar聽 37.Manabe, M. Fermented foods and mycotoxins. JSM Mycotoxins 51, 25鈥?9 (2001).CAS聽 Article聽Google Scholar聽 38.Kaminishi, Y., Egusa, J. M, K. Antioxidant production from several xerophilous fungi used in 鈥淜atsuobushi鈥?molding. Journal of National Fisheries University 47, 113鈥?20 (1999).CAS聽Google Scholar聽 39.Kurtzman, C. P., Robnett, C. J. Basehoar-Powers, E. Phylogenetic relationships among species of Pichia, Issatchenkia and Williopsis determined from multigene sequence analysis, and the proposal of Barnettozyma gen. nov., Lindnera gen. nov. and Wickerhamomyces gen. nov. FEMS yeast research 8, 939鈥?54 (2008).CAS聽 Article聽 PubMed聽Google Scholar聽 40.Kurtzman, C. P., Fell, J. W. Boekhout, T. In The Yeasts (Fifth Edition) 137鈥?44 (Elsevier, 2011).41.Kurtzman, C. P., Fell, J. W., Boekhout, T. Robert, V. In The Yeasts (Fifth Edition) 87鈥?10 (Elsevier, 2011).42.Minter, D. Cyberlindnera, a replacement name for Lindnera Kurtzman et al., nom. illegit. Mycotaxon 110, 473鈥?76 (2009).Article聽Google Scholar聽 43.Chang, C. F. et al. Five novel species of the anamorphic genus Candida in the Cyberlindnera clade isolated from natural substrates in Taiwan. Antonie Van Leeuwenhoek 102, 9鈥?1 (2012).Article聽 PubMed聽Google Scholar聽 44.Kurtzman, C. P., Johnson, C. J. Smiley, M. J. Determination of conspecificity of Candida utilis and Hansenula jadinii through DNA reassociation. Mycologia 11, 844鈥?47 (1979).CAS聽 Article聽 PubMed聽Google Scholar聽 45.Bekatorou, A., Psarianos, C. Koutinas, A. A. Production of Food Grade Yeasts. Food Technology Biotechnology 44, 407鈥?15 (2006). Google Scholar聽 46.Kanpiengjai, A., Chui-Chai, N., Chaikaew, S. Khanongnuch, C. Distribution of tannin-鈥榯olerant yeasts isolated from Miang, a traditional fermented tea leaf (Camellia sinensis var. assamica) in northern Thailand. International Journal of Food Microbiology 238, 121鈥?31 (2016).CAS聽 Article聽 PubMed聽Google Scholar聽 47.Meyer, S. A., Payne, R. W. Yarrow, D. In The Yeasts (Fourth Edition) (ed Jack W. Fell) 454鈥?73 (Elsevier, 1998).48.Chen, K., Lin, Y., Chen, W. Liu, Y. Degradation of phenol by PAA-immobilized Candida tropicalis. Enzyme and Microbial Technology 31, 490鈥?97 (2002).CAS聽 Article聽Google Scholar聽 49.Yan, J., Jianping, W., Hongmei, L., Suliang, Y. Zongding, H. The biodegradation of phenol at high initial concentration by the yeast Candida tropicalis. Biochemical Engineering Journal 24, 243鈥?47 (2005).Article聽Google Scholar聽 50.Adav, S. S., Chen, M. Y., Lee, D. J. Ren, N. Q. Degradation of phenol by aerobic granules and isolated yeast Candida tropicalis. Biotechnology and bioengineering 96, 844鈥?52 (2007).CAS聽 Article聽 PubMed聽Google Scholar聽 51.Varma, R. J. Gaikwad, B. G. Rapid and high biodegradation of phenols catalyzed by Candida tropicalis NCIM 3556 cells. Enzyme and Microbial Technology 43, 431鈥?35 (2008).CAS聽 Article聽Google Scholar聽 52.Horitsu, H. et al. Production of xylitol from D-xylose by Candida tropicalis: optimization of production rate. Biotechnology and bioengineering 40, 1085鈥?091 (1992).CAS聽 Article聽 PubMed聽Google Scholar聽 53.Garc铆a Mart铆n, J. F., Cuevas, M., Bravo, V. S谩nchez, S. Ethanol production from olive prunings by autohydrolysis and fermentation with Candida tropicalis. Renewable Energy 35, 1602鈥?608 (2010).Article聽Google Scholar聽 54.Feng, J., Zhan, X. B., Wang, D., Zhang, L. M. Lin, C. C. Identification and analysis of the metabolic functions of a high-salt-tolerant halophilic aromatic yeast Candida etchellsii for soy sauce production. World journal of microbiology biotechnology 28, 1451鈥?458 (2012).CAS聽 Article聽Google Scholar聽 55.Liu, Z. et al. In The technique specification for raw Anhua dark tea processing. Vol. DB43/T, 660鈥?011 (China Agriculture Press, Changsha, 2011). Google Scholar聽 56.Zhou, W., Xu, J., Lu, X., Wang, J. Sha, H. Tea-Determination of water extract content. Vol. GB/T 8305鈥?013 (China Agriculture Press, Beijing, 2013). Google Scholar聽 57.Gong, L. Xu, Q. Determination of total acid in foods . Vol. GB/T 12456鈥?008 (China Agriculture Press, Beijing, 2008). Google Scholar聽 58.Zhou, W., Xu, J. Xu, L. Determination of total polyphenols and catechins content in tea. Vol. GB/T 8313鈥?008 (China Agriculture Press, Beijing, 2008). Google Scholar聽 59.Xu, J., Zhou, W., Lu, X., Wang, J. Sha, H. Tea-Determination of free amino acids content. Vol. GB/T 8314鈥?013 (China Agriculture Press, Beijing, 2013). Google Scholar聽 60.Li, Y., Li, J., Gong, X. Liu, Z. Simultaneous determination of eight catechins, three purine alkaloids and gallic acid in tea by high-performance liquid chromatography. Food Science 32, 214鈥?17 (2011).CAS聽Google Scholar聽 61.Petry, R. D., Ortega, G. G. Silva, W. B. Flavonoid content assay: influence of the reagent concentration and reaction time on the spectrophotometric behavior of the aluminium chloride鈥揻lavonoid complex. Die Pharmazie 56, 465鈥?70 (2001).CAS聽 PubMed聽Google Scholar聽 62.Mengxi, L., Zhigang, X., Yang, Y. Yijie, F. Effects of different spectral lights on Oncidium PLBs induction, proliferation, and plant regeneration. Plant Cell, Tissue and Organ Culture (PCTOC) 106, 1鈥?0 (2011).Article聽Google Scholar聽 63.Mukherjee, P. K. et al. Oral mycobiome analysis of HIV-infected patients: identification of Pichia as an antagonist of opportunistic fungi. PLoS pathogens 10, e1003996 (2014).Article聽 PubMed聽 PubMed Central聽Google Scholar聽 64.Bolger, A. M., Lohse, M. Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114鈥?120 (2014).CAS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 65.Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and environmental microbiology 75, 7537鈥?541 (2009).CAS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 66.Braak, C. J. F. t. Smilauer, P. CANOCO Reference Manual and CanoDraw for Windows User鈥檚 Guide: Software for Canonical Community Ordination (version 4.5). (www.canoco.com, Ithaca NY, USA, 2002).Download referencesAcknowledgementsThe work was financially supported by the National Natural Science Foundation of China (31471706), Natural Science Foundation of Hunan Province (13JJ4067), Scientific Research Fund of Hunan Agricultural University (13YJ13) and the 鈥?515 Talent Project鈥?of Hunan Agricultural University, Open Fund Research Project of Hunan Provicial Key Laboratory of Crop Germplasm Innovation and Utilization (16KFXM11).Author informationAffiliationsKey Laboratory of Tea Science of Ministry of Education, Hunan Agricultural University, Changsha, Hunan, 410128, P.R. ChinaQin Li,聽Jianan Huang,聽Yongdi Li,聽Yiyang Zhang聽 聽Kunbo WangHunan Provincial Key Laboratory of Crop Germplasm Innovation and Utilization, Hunan Agricultural University, Changsha, Hunan, 410128, P.R. ChinaQin Li,聽Jianan Huang聽 聽Zhonghua LiuNational Research Center of Engineering Technology for Utilization of Functional Ingredients from Botanicals, Hunan Agricultural University, Changsha, Hunan, 410128, P.R. ChinaQin Li,聽Haiyan Lin,聽Kunbo Wang聽 聽Zhonghua LiuCollaborative Innovation Centre of Utilization of Functional Ingredients from Botanicals, Hunan Agricultural University, Changsha, Hunan, 410128, P.R. ChinaQin Li,聽Jianan Huang,聽Kunbo Wang聽 聽Zhonghua LiuInstitute of Soil and Water Resources and Environmental Sciences, Zhejiang University, Hangzhou, Zhejiang, 3100058, P.R. ChinaYu LuoCollege of Plant Protection, Hunan Agricultural University, Changsha, Hunan, 410128, P.R. ChinaYuan ChenAuthorsQin LiView author publicationsYou can also search for this author in PubMed聽Google ScholarJianan HuangView author publicationsYou can also search for this author in PubMed聽Google ScholarYongdi LiView author publicationsYou can also search for this author in PubMed聽Google ScholarYiyang ZhangView author publicationsYou can also search for this author in PubMed聽Google ScholarYu LuoView author publicationsYou can also search for this author in PubMed聽Google ScholarYuan ChenView author publicationsYou can also search for this author in PubMed聽Google ScholarHaiyan LinView author publicationsYou can also search for this author in PubMed聽Google ScholarKunbo WangView author publicationsYou can also search for this author in PubMed聽Google ScholarZhonghua LiuView author publicationsYou can also search for this author in PubMed聽Google ScholarContributionsQ.L., Z.H.L. and J.A.H. designed the research. Q.L. conducted the experiments, analyzed the data, and drafted the manuscript. Y.D.L. and Y.Y.Z. helped carry out the metabolic analysis. Y.L. and Y.C. helped analyze the data and revise the manuscript. K.B.W. and H.Y.L. helped prepare and collect the tea samples.Corresponding authorCorrespondence to Zhonghua Liu.Ethics declarations 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.Electronic supplementary material Supplementary infromationRights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article鈥檚 Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article鈥檚 Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and PermissionsAbout this articleCite this articleLi, Q., Huang, J., Li, Y. et al. Fungal community succession and major components change during manufacturing process of Fu brick tea. Sci Rep 7, 6947 (2017). https://doi.org/10.1038/s41598-017-07098-8Download citationReceived: 10 April 2017Accepted: 22 June 2017Published: 31 July 2017DOI: https://doi.org/10.1038/s41598-017-07098-8 Wenping Kong, Chengshuang Huang, Jie Shi, Yu Li, Xinxin Jiang, Quwen Duan, Yong Huang, Yanwen Duan Xiangcheng Zhu Microbial Cell Factories (2019) CommentsBy submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Sign up for the Nature Briefing newsletter 鈥?what matters in science, free to your inbox daily.