Meichun Xiang1,2*
Jun Chu1
Wenjiao Cai1
Haikun Ma3
Weijing Zhu1Xiaoling Zhang1
Jinwei Ren1
Lizheng Xiao4
Dongpo Liu4Xingzhong Liu3*- 1State Key Laboratory of Mycology, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3State Key Laboratory of Medicinal Chemical Biology, Key Laboratory of Molecular Microbiology and Technology of the Ministry of Education, Department of Microbiology, College of Life Sciences, Nankai University, Tianjin, China
- 4Department of Tea Science, College of Horticulture, Hunan Agricultural University, Changsha, China
Fu Brick tea is a very popular post-fermented tea that is known for its “golden flower fungus,” Aspergillus cristatus, which becomes the dominant microbe during the maturation process. This study used both culture-dependent methods and high-throughput sequencing to track microbial succession and interactions during the development of the golden flower fungus, a crucial component of the manufacturing process of Fu Brick tea. Among the bacterial communities, Klebsiella and Lactobacillus were consistently cultured from both fresh tea leaves and in post-fermentation Fu Brick tea. Methylobacterium, Pelomonas, and Sphingomonas were dominant genera in fresh tea leaves but declined once fermentation started, while Bacillus, Kluyvera, and Paenibacillus became dominant after piling fermentation. The abundance of A. cristatus increased during the manufacturing process, accounting for over 98% of all fungi present after the golden flower bloom in the Fu Brick tea product. Despite their consistent presence during culture work, network analysis showed Lactobacillus and Klebsiella to be negatively correlated with A. cristatus. Bacillus spp., as expected from culture work, positively correlated with the presence of golden flower fungus. This study provides complete insights about the succession of microbial communities and highlights the importance of co-occurrence microbes with A. cristatus during the manufacturing process of Fu Brick tea.
Introduction
Fu Brick tea is made from the leaves of Camellia sinensis var. sinensis after undergoing indigenous microbial fermentation. It is traditionally produced in the Hunan Province of China, and it dates back to the Ming Dynasty (Zhang et al., 2013). Fu Brick tea is consumed widely, especially as a necessary component of milk tea or buttered tea for the nomadic groups in the border regions of southern and western China (Chen, 2000). The health benefits of drinking Fu Brick tea include suppression of fatty acid synthase expression (Chiang et al., 2006), inhibition of lipid, and non-lipid oxidative damage (Jie et al., 2006), as well as anti-obesity (Li et al., 2013; Chen et al., 2018), anti-hyperlipidemia (Fu et al., 2011), and anti-hyperglycemia (Wu et al., 2010). Unlike the other types of fermentation that depend on the enzyme oxidation induced by tea leaves, the post-fermentation process of Fu Brick tea depends on the auto-oxidation and non-enzyme auto-oxidation induced by the indigenous microorganisms on tea leaves (Zheng et al., 2015). The unique hallmark for the quality of Fu Brick tea is the “golden flower fungus,” which earned this name because of the golden cleistothecia produced by the fungus Aspergillus cristatus (sexual morph is Eurotium-like) (Xu et al., 2011; Chen G. et al., 2017). The formation of small golden spheres throughout the interior of Fu Brick tea during fermentation is crucial for achieving high-quality aroma and taste (Liu et al., 2019).
The manufacture of Fu Brick tea generally involves two separate processes. The first process converts fresh tea leaves to raw dark tea through withering, steaming, rolling, pile fermentation, and drying. The second process (likely occurring in a separate facility) subjects the raw dark tea to blending, steaming, short-term piling, compression into bricks, and fungal fermentation (e.g., golden flower blooming) (Zheng et al., 2015). Pile fermentation and golden flower blooming are two key steps involved in the microbial fermentation process (Liu et al., 2019). Pile fermentation is normally conducted under high humidity (over 86% RH) with temperature above 25°C for up to 40 h to enhance the microbial activities in tea leaves (Zheng et al., 2015). Golden flower blooming occurs with relative humidity under 70% and temperature above 25°C for around 21 days (Wen and Liu, 1991). Several investigations have investigated the dynamics of microbial communities during the Fu Brick manufacturing process that focused solely on either fungi (Li et al., 2017), bacteria (Li et al., 2019), or the pile fermentation process (Li et al., 2018). Microbes other than A. cristatus have shown prevalence during the pile fermentation stage, such as fungal species of Cyberlindnera and Candida, or bacterial species from the Methylobacterium, Klebsiella, and Lactobacillus genera (Li et al., 2018). This indicates that lactic acid bacteria and yeasts are involved in the fermentation of Fu Brick tea (Zhang et al., 2013). Thus, it is necessary to know the relationships between these prevalent microbes during the manufacturing processes for Fu Brick tea.
In this study, tea samples were collected from different factories and at six stages of processing. By combining culture-dependent methods with high-throughput sequencing, we assessed the dynamics and succession of bacterial and fungal communities in Fu Bricks during the manufacturing process. Specifically, the objectives were to (1) assess the dominant microbes present and (2) investigate microbial interactions at each step of the Fu Brick manufacturing process. Our findings could advance the understanding about microbial fermentation of Fu Brick tea and may offer means for improving the manufacturing process.
Materials and Methods
Sample Collection and Preparation
A total of 90 tea samples were collected at different processing stages from four companies in Anhua, Hunan Province, and one processing facility in Meitan, Guizhou Province, China. The sampled fermentation stages were (in chronological order) as follows: fresh tea (9 samples), pile processing (14 samples), drying (14 samples), early stage golden flower bloom (2–12 h) (22 samples), full golden flower bloom (1–20 days) (23 samples), and post-fermentation (8 samples; Figure 1). The first three stages cover the processing of fresh tea leaves to raw dark tea, and the last three stages cover the processing of raw dark tea to Fu Bricks (Zheng et al., 2015). The pile fermentation and golden flower stages were each sampled at different time points, within each stage, to more accurately assess changes to the microbial community, which were averaged to obtain an overall effect during each stage. The samples were placed in sterile polyethylene bags and transported on ice to the laboratory. All the samples were portioned immediately, with one stored portion at 4°C for culturing of microorganisms and the other stored at −80°C for DNA extraction.
Figure 1. The manufacturing stages of Fu Brick tea. (A) The processing of fresh tea leaves to make raw dark tea. (B) The processing of raw dark tea to make Fu Bricks.
Culture-Independent Method: DNA Extraction and Sequencing
For each sample, 2.5 g of DNA was extracted with the MoBio PowerSoil® DNA Isolation Kit (MO BIO Laboratories, Inc.) by following the manufacturer's protocols. The extracted DNA was stored at −20°C for further PCR amplification and barcoded amplicon sequencing. Primer sets ITS1-F (5'-CTTGGTCATTTAGAGGAAGTAA-3') and ITS2-R (5'-GCTGCGTTC TTCATCGATGC-3') were used for the amplification of fungal DNA, and primer sets V4-F (5'-GTGCCAGCMGCCGCGGTAA-3') and V4-R (5'-GGACTACHVGGGTWTCTAAT-3') were used for the amplification of bacterial DNA. Both the forward and reverse primers had 6-bp barcodes unique to each sample, which were used to permit multiplexing of samples. The DNA concentrations were adjusted to 20 ng/μl for each sample in the PCR. Each 50 μl PCR solution consisted of 25 μl 2x Taq PCR MasterMix (TIANGEN Company, Beijing), 1.5 μl sense primer (10 μM), 1.5 μl anti-sense primer (10 μM), 3 μl DNA template, and 19 μl ddH2O. PCR amplifications were performed in triplicate reactions on a Veriti thermal cycler (Applied Biosystems). PCR amplifications were performed under the following conditions: initial denaturing at 94°C for 4 min, followed by 35 cycles of 94°C for 40 s, 54°C for 45 s (for ITS), or 52°C for 45 s (for V4), 72°C for 50 s, and a final extension step at 72°C for 10 min. Negative controls with no template DNA were also included in PCR amplifications. The products of the triplicate PCRs of each soil sample were pooled and subjected to gel electrophoresis, from which purified amplicons were harvested using the EasyPure Quick Gel Extraction Kit (TransGen, Beijing). The purified amplicons and products encoding different barcodes were combined into a single library for automated cluster generation and paired-end sequencing on the MiSeq system.
Primary data analysis (image analysis, base calling) was performed on a MiSeq instrument. The sequences were analyzed using the Quantitative Insights Into Microbial Ecology (QIIME) toolkit with default parameters for each step (Caporaso et al., 2010). Fungal and bacterial sequences with lengths shorter than 220 nucleotides were removed. De novo and reference-based chimera checking were performed, and sequences that were characterized as chimeric were also removed. Bacterial reads were binned into OTUs at the ≥97% sequence similarity level using the open-reference OTU picking protocol in UPARSE-pipeline (Edgar, 2013), and the most abundant sequences from each OTU were taken as representative sequences for the respective OTU. Taxonomic configuration of bacterial OTUs was performed by searching the Basic Local Alignment Search Tool (BLAST) for each representative sequence against a subset of the Silva database (Quast et al., 2013). The fungal sequences were subjected to ≥97% sequence similarity single-linkage clustering using GenBank (Benson et al., 2005) and UNITE (Abarenkov et al., 2012) as reference taxonomic databases.
Culture-Dependent Method
For each sample, 10 g was homogenized with sterile distilled water using a blender at 12,000 rpm using three 30 s pulses. The homogenate was serially diluted in sterile distilled water and spread onto each of four types of media: MEA (2% malt extract, 0.1% peptone, 2% glucose, 1.8% agar; containing 100 mg/L chloramphenicol and 100 mg/L streptomycin), LA medium (BD Co.), YDP (1% yeast extract, 1% peptone, 2% glucose, 1.8% agar; containing 100 mg/L chloramphenicol), and Gauze No.1 (AOBOX BIOTECHNOLOGY; 25 mg/L nalidixic acid and 50 mg/L cycloheximide) (Liu et al., 1996; Rutiaga-Quiñones et al., 2012). Plates were prepared for each sample and medium in triplicate. The plates were incubated at 25°C for 2 days, and colonies were transferred to new plates to purify and preserve.
For the identification of microbes from pure cultures, morphological examinations were conducted with microscopy. These identifications were then followed by the implementation of DNA-based methods. The genomic DNA was extracted with a simple and rapid “thermolysis” method (Zhang et al., 2010). Briefly, a sterile toothpick was used to transfer mycelia or yeast cells from the colony into 100 μl pure water. The mixture was vortexed thoroughly and centrifuged at 10,000 g for 1 min. The supernatant was removed, and then, 100 μl of lysis solution [EasySelectTM Pichia Expression Kit (Invitrogen, USA)] was added into the microcentrifuge tube. The mixture was then incubated at 85°C for 20–30 min. The extract contained genomic DNA and was stored at −20°C for future use. Multiple primers are used for confirmation of our morphological identifications; they are listed in Supplementary Table 1, and they included the ITS region of fungi (White et al., 1990), the D1/D2 region of yeasts (Boekhout et al., 1995), the 16S rRNA gene of bacteria with degenerate bases (DeLong et al., 2006), and the 16S rRNA gene of actinomycetes without degenerate bases (Lane, 1991). Special primers were used for better differentiating yeast and actinomycetes species based on literature. Each 25 μl PCR solution consisted of 12.5 μl 2x Taq PCR MasterMix (TIANGEN Company, Beijing), 1 μl sense primer (10 μM), 1 μl anti-sense primer (10 μM), 2 μl DNA template, and 8.5 μl ddH2O. PCR amplifications were performed on a Veriti Thermal Cycler (Applied Biosystems) under the following conditions: initial denaturing at 94°C for 3 min, followed by 35 cycles at 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min, with a final extension step at 72°C for 10 min (fungi); initial denaturing at 94°C for 3 min, followed by 35 cycles at 94°C for 40 s, 55°C for 30 s, and 72°C for 80 s, with a final extension step at 72°C for 10 min (bacteria); initial denaturing at 94°C for 3 min, followed by 35 cycles at 94°C for 30 s, 55°C for 30 s and 72°C for 1 min, with a final extension step at 72°C for 10 min (yeast); initial denaturing at 94°C for 3 min, followed by 35 cycles at 94°C for 30 s, 55°C for 30 s, and 72°C for 90 s, with a final extension step at 72°C for 10 min (actinomycetes). The PCR products were sent to Omega Genetics Co. Ltd., China, for downstream sequencing, and the obtained sequences were uploaded to GenBank database to search for homologous sequences by BLAST.
Statistical Analysis
Sequencing data were normalized using total sum scaling (Weiss et al., 2017). Culture-dependent data were shown as counts (CFU/g). All data were checked for homogeneity of variance, and normality was confirmed by inspection of residuals. The overall effects of manufacturing stages on the composition of bacterial and fungal communities were examined using permutational multivariate ANOVA (PERMANOVA) based on a Bray-Curtis dissimilarity matrix (Oksanen et al., 2006). Non-metric multidimensional scaling (NMDS) plots, also based on Bray-Curtis, were used to visualize the separations between manufacturing stages. Linear mixed models were used to examine the effects of the manufacturing process on the Simpson diversities of bacteria and fungi, and the relative abundances of bacterial phyla, fungal classes, A. cristatus, yeast, and lactic acid bacteria. In the model, processing stages were used as the fixed factors, and different sampling facilities were used as random factors. Then, post-hoc Tukey's tests were used for the pairwise comparisons between all processing stages. The linear mixed model and post-hoc Tukey's test described earlier were also used to examine the effects of processing stages on the relative abundances of core microbes (those that dominated in 80% of all samples or had relative abundances >1%). Network analysis based on Pearson's correlations (P < 0.05) was used to explore only the relationships between core microorganisms because too few sample sizes would lead to false-positive correlations. The network was constructed using Gephi-0.9.2. Given the likelihood that not all microbes present in our samples would be culturable during the culture-dependent experiments and would result in missing values in the culture-dependent dataset, statistical analysis may not be possible for the culture-dependent dataset. All analyses were performed in R (version 3.0.1, R Development Core Team, 2013).
Results
Culture-Independent Analyses Reveal Significant Fluctuations in Microbial Composition
Fermentation stage one explained 13% of the variances in the bacterial community and 17% in the fungal community (Table 1). Fermentation stage two did not significantly influence the microbial communities. The overall effects of fermentation stages explained 16% of the variances in the bacterial community and 15% in the fungal community (Table 1 and Figure 2). Fresh tea leaves contained more distinct fungal and bacterial communities compared with those present during pile fermentation (Figures 2A,B). Manufacturing stages significantly influenced the Simpson diversity of bacterial and fungal communities (Figures 2C,D). With regard to bacterial diversity, pile fermentation had the lowest level, while early stage golden flower bloom had the highest (Figure 2C). For fungi, the drying stage had the highest amount of diversity, while the Fu Brick tea product (i.e., the final stage of processing) had the lowest (Figure 2D).
Table 1. The effects of fermentation stages on the composition of bacterial and fungal communities.
Figure 2. Assessing microbial diversity at different stages of the Fu Brick manufacturing process. Non-metric Multidimensional Scaling (NMDS) plots with Bray-Curtis dissimilarity distances for bacterial (A) and fungal (B) communities, and bar plots showing statistical results following linear mixed models for bacterial (C) and fungal (D) diversity. Standard errors from the mean for each stage are shown with each bar, and letters above bars indicate significant differences between manufacturing stages. “GF bloom (h)” is the early stage golden flower bloom and “GF bloom (d)” is the full golden flower bloom.
Across all stages of the manufacturing process, the dominant phyla comprising bacterial communities included the Proteobacteria, Firmicutes, Bacteroidetes, Actinobacteria, and Acidobacteria (Figure 3A and Supplementary Table 2). Fresh tea leaves had the highest relative abundances of Proteobacteria and Acidobacteria and the lowest relative abundance of Firmicutes (Supplementary Table 2). Throughout both fermentation stages, however, the relative abundances of Proteobacteria and Acidobacteria decreased, while the relative abundance of Firmicutes increased (Supplementary Table 2). The final Fu Brick product contained the highest relative abundance of Firmicutes and the lowest relative abundances of Proteobacteria and Acidobacteria (Supplementary Table 2).
Figure 3. The relative abundances of bacterial phyla (A) and fungal classes (B) in each examined stage of the Fu Brick manufacturing process. “GF bloom (h)” is the early stage golden flower bloom, and “GF bloom (d)” is the full golden flower bloom.
At the genus level, the bacterial communities on fresh tea leaves were dominated by Klebsiella, followed by Methylobacterium. The pile fermentation and drying stages increased the relative abundances of Klebsiella, Sphingomonas, Paenibacillus, and Streptomyces, which then decreased at the golden flower bloom stages (Table 2). Throughout the fermentation stages, the relative abundance of Methylobacterium decreased significantly and reached its lowest level in the Fu Brick product (Table 2). Conversely, the abundances of Bacillus and Pelomonas species were lowest in fresh tea leaves, but they significantly increased during pile fermentation and reached their highest levels in the Fu Brick product (Table 2).
Table 2. The relative abundances of core microbes at each stage of the Fu Brick manufacturing process.
The dominant classes within fungal communities throughout the manufacturing process consisted of Eurotiomycetes, Sordariomycetes, Saccharomycetes, and Dothideomycetes (Figure 3B and Supplementary Table 2). Fresh tea leaves contained the lowest relative abundance of Eurotiomycetes and the highest relative abundances of Dothideomycetes and Sordariomycetes (Supplementary Table 3). As subsequent processing stages progressed, the relative abundance of Eurotiomycetes increased and the relative abundances of Dothideomycetes and Sordariomycetes decreased, with Fu Bricks containing the highest relative abundance of Eurotiomycetes and the lowest relative abundances of Dothideomycetes and Sordariomycetes (Supplementary Table 3). The relative abundance of Saccharomycetes first increased, and became highest in the pile fermentation stage, and then decreased, with Fu Bricks containing the lowest relative abundance of Saccharomycetes (Supplementary Table 3).
So far, Aspergillus, especially A. cristatus, was the dominant fungal genus, which increased significantly across all processing stages from fresh tea leaves (low levels) to comprise 98% of the Fu Brick product (Table 2). Other fungal genera predominated in fresh tea leaves, such as Penicillium, Purpureocillium, and Cladosporium. As fermentation progressed, their relative abundances fluctuated. By the drying stage, Penicillium and Purpureocillium species had the lowest relative abundances (Table 2). Both fermentation stages caused consistent decreases in the relative abundance of Cladosporium fungi with the Fu Brick product containing the lowest relative abundance of this genus (Table 2). Despite being considered a dominant fungal class, no genus within the Saccharomycetes satisfied the two standards of core microbes for this study.
Culture-Dependent Analyses Support the Fluctuation of Dominant Microbes During the Manufacturing Process
In this study, we classified the microbes as A. cristatus, yeasts, or lactic acid bacteria based on their colony forming units (cfu) (Figure 4). Overall, we found that the culture-based data at each processing stage supported the sequencing data: fresh tea leaves contained the lowest abundance of A. cristatus cfu, which increased along the manufacturing process to attain the greatest abundance in the Fu Brick product (Figures 4A,E). Yeasts were abundant in fresh tea leaves, but during fermentation, their numbers decreased. This fate continued through the final Fu Brick stage where they had the lowest abundance (Figures 4B,F). The culture-dependent data showed that all isolated bacteria cfu (including the lactic acid bacteria) were highest at the drying stage and lowest in the Fu Brick product (Figure 4G), while the sequencing data suggested that the relative abundances of lactic acid bacteria did not significantly fluctuate during the manufacturing process (Figure 4C).
Figure 4. Bar plots showing the relative abundances (A–C) and colony forming units (D–F) of A. cristatus golden flower fungus, yeasts, and lactic acid bacteria, respectively, based on culture-dependent analysis at each examined stage of Fu Brick manufacturing. Standard errors from the mean for each stage are shown with each bar. Statistical results following linear mixed models (A-C) are presented above each panel. “n.s.” indicates no significant effect. Letters above bars indicate significant differences between manufacturing stages. “GF bloom (h)” is the early stage golden flower bloom and “GF bloom (d)” is the full golden flower bloom.
Network Analyses Reveal Important Relationships Among Core Microbial Groups
Network analysis was conducted for the microbes that presented over 80% of all samples (Figure 5). Klebsiella, Staphylococcus, and Acinetobacter were in the core of the network, Klebsiella negatively correlated with most other microbes, while Staphylococcus and Acinetobacter had both positive and negative relationships with other microbes. Bacillus, Penicillium, and Streptomyces positively correlated with A. cristatus, but Lactobacillus, Klebsiella, and Cyberlindnera negatively correlated with this fungus.
Figure 5. Network graph showing relationships between core microbes (bacteria and fungi) presented in Table 2. Larger nodes indicate higher connections with other genera. Red connections indicate negative correlations, while green connections indicate positive correlations. Pink nodes indicate bacterial genera and green nodes indicate fungal genera.
Discussion
Many beneficial qualities of Fu Brick tea exist because of fermentation microbes and their metabolites (Zheng et al., 2015), especially the golden flower fungus, which is inherently present on fresh tea leaves and becomes the dominant organism solely through indigenous microbial succession during the manufacturing process. In this study, both the high-throughput sequencing and culture-dependent data revealed drastic changes to microbial community structure that continued to fluctuate throughout each stage of the manufacturing process. Many bacterial genera that were present in fresh tea leaves were lost during the manufacturing process. However, the relative abundance of Bacillus bacteria increased as each stage of manufacturing progressed. In fact, this bacterial genus was observed to be positively correlated with the increase in abundance of A. cristatus. The positive correlation between Bacillus and A. cristatus indicated that these bacteria might facilitate and enhance the development of the golden flower fungus. Although both Penicillium and Aspergillus were the dominant genera in fresh tea leaves, only the abundance of A. cristatus continued to increase at each stage of the manufacturing process until it became the dominant microbe, an indicator of high-quality Fu Brick tea.
The succession in bacterial and fungal communities of Fu Brick tea has been extensively investigated (Xu et al., 2011; Li et al., 2017, 2018, 2019). However, those studies mainly correlated microbiota with the biochemical components of Fu Brick tea. In this study, we focused not only on the golden flower fungus but also on the co-occurrence of fungi and bacteria, thus providing a comprehensive picture of the fluctuation of the microbiome during the manufacturing processes. Moreover, our study explored microbial community changes at different stages of the Fu Brick manufacturing process. In other studies, samples were collected at a single time point, such as in the pile fermentation stage, which lasts for 24 h, or at different time points but only in the golden flower blooming stage, which takes around 20 days (Xu et al., 2011; Zheng et al., 2015; Li et al., 2019). Not including microbial communities from the first fermentation stage (piling), in addition to the second fermentation stage (golden flower bloom), leads to an unbalanced representation of the manufacturing process where the pile fermentation microbes would be overlooked, and the golden flower fermentation microbes would be over-represented (as if they are present at all stages). In this study, we offer a more accurate representation of the dynamic microbial changes that occur over the entire processing of Fu Brick tea by sampling several time points in both stages of fermentation. Due to our thorough sampling, we illustrated the variable succession of microbial communities throughout the Fu Brick manufacturing process. More importantly, we showed that the golden flower blooming stage had the greatest influence on the composition of bacterial communities and that abundance of one bacterial genus, in particular, Bacillus, was directly correlated with the increase of A. cristatus, the golden flower fungus.
Previous investigations associating the Fu Brick microbial community with biochemical components of the tea provided essential information to support our observations (Li et al., 2017, 2018, 2019). For example, Klebsiella had the highest relative abundance on both fresh tea leaves and the final Fu Brick tea product despite decreases in abundance during the other processing stages. Previous studies have shown the predominant role of Klebsiella during the Fu Brick manufacturing process (Liu et al., 2019), associating this bacterium with the metabolism of biochemical components in Fu Bricks, such as polyphenol, caffeine, and gallocatechin gallate (Li et al., 2019), as well as the process of saccharification by producing ethanol (Doran et al., 1994; Zhou and Ingram, 2001; Golias et al., 2002) and butanediol (Cao et al., 1997; Sun et al., 2009; Moon et al., 2018), especially in co-culture with fungi (Yu et al., 1985; Cao et al., 1997). Other bacterial genera (Methylobacterium, Pelomonas, and Sphingomonas) had high relative abundances on fresh tea leaves, which greatly decreased during the Fu Brick manufacturing process. Methylobacterium is known for the production of poly-β-hydroxybutyrate (Bourque et al., 1992; Yellore et al., 2000; Ghatnekar et al., 2002) and was identified as one of the core bacteria during the pile fermentation processing stage (Liu et al., 2019). It is also one of the predominant genera of raw Pu'er tea, which is another type of fermented tea derived from the leaves of C. sinensis var. sinensis (Ma et al., 2017). Although Pelomonas was only recently discovered in Fu Brick tea (Li et al., 2018, 2019), it has been found in other fermented foods, such as kefir grains (Gao et al., 2013), Shaoxing liquor (Chen et al., 2020), and fermented grape, where it was associated with the production of lactic and acetic acids (del Carmen Portillo and Mas, 2016). Thus, the presence of Pelomonas indicates the flavor of the food products. Three bacterial genera (i.e., Bacillus, Kluyvera, and Paenibacillus) increased in abundance to become dominant during the manufacturing process of Fu Brick tea. Bacillus can produce pectinase (Rehman et al., 2015; Uzuner and Cekmecelioglu, 2015) and oxidase polyphenol (Mohammad and Alireza, 2007), which are important for oxidating the major catechins into theaflavin (Zhu et al., 2015). Bacillus has been linked to the flavor profile of Fu Brick (Zheng et al., 2015; Liu et al., 2019). Moreover, a recent study has found that Bacillus correlated with 44 metabolites of Fu Brick tea, highlighting the importance of Bacillus in manufacturing Fu Brick tea (Xia et al., 2022). Although Pseudomonas were not significantly changed between different manufacturing stages in our study, they were revealed by metagenomic methods as the most important bacteria associated with 35 probiotics in Fu Brick tea (Wang et al., 2021). Other dominant bacteria, such as Paenibacillus, were reported to be one of the key microbes in the manufacture of other traditionally fermented foods, including Pu'er tea (Oh et al., 2008; Kim et al., 2009; Zhao et al., 2015).
With regard to fungal genera, A. cristatus accounted for about 50% of all fungi in fresh tea leaves and increased up to 98% in the Fu Brick tea product. This pattern was further confirmed by the culture-dependent data, whereby A. cristatus was the dominant fungus. Some species of Aspergillus are commonly used in the traditional fermentation of foods where they reportedly degrade plant cell walls, produce peroxiredoxins, and contribute to the oxidation of catechins, such as during the manufacture of Pu'er tea (Chen A. J. et al., 2017). As the uniquely dominant fungus in Fu Brick tea, A. cristatus was reported to be responsible for the production of polyphenol oxidase, cellulase, and pectinase (Li et al., 2018), as well as serving as an important probiotic that reduces obesity in humans (Kang et al., 2019). Purpureocillium and Cyberlindnera were also dominant fungi in the fresh tea leaves used for Fu Brick tea, accounting for 28 and 6% of all fungal genera observed, respectively. The presence of Purpureocillium fungi in Fu Brick tea has not been previously reported; however, they have been observed in other fermented foods, such as beans (Serra et al., 2019), Chinese Paocai (Xiao et al., 2018), and liquor (Song et al., 2017), as well as being associated with the production of amino acids (Xiao et al., 2018). Cyberlindnera has been reported as one of the core microbes observed during the piling fermentation process to create raw dark tea (Li et al., 2017; Liu et al., 2019). As a genus of yeast, Cyberlindnera may facilitate the bioconversion of sugars into ethanol, CO2, and other metabolites during the piling fermentation process of Fu Brick manufacture (Tofalo et al., 2019).
Generalized relationships between bacteria and fungi in Fu Brick tea have been reported (Rui et al., 2019). Our study provided further information about the dominant bacterial genera that are highly correlated with fermentation fungi (i.e., A. cristatus) throughout the entire manufacturing process. As the golden flower fungus has been extensively credited for improving the flavor and quality of Fu Brick tea in the fermentation process (Angayarkanni et al., 2002; Wulandari et al., 2011; Yao et al., 2017; Jiang et al., 2018), it is important to identify any microbes that may facilitate its growth and development. However, the mechanisms underscoring these relationships between A. cristatus and dominant Fu Brick bacteria remain unknown. Understanding these relationships throughout the entire manufacturing process should help improve tea quality and streamline the Fu Brick manufacturing process.
Conclusion
The succession of Fu Brick core microbes over time indicated that certain bacterial and yeast genera do not survive the fermentation stages, especially in the presence of the golden flower fungus, A. cristatus, that greatly increases its infiltration of the tea product during fermentation. However, some bacterial and yeast genera do survive and increase in direct proportion with A. cristatus. Although the golden flower fungus has long been considered the most important microbial component in Fu Brick tea, other microbes (e.g., Bacillus) that positively correlated with A. cristatus growth also deserve further investigation to explore the causal relationships between fermentation-thriving organisms in Fu Bricks. Associating microbial succession and interactions that affect the biochemical components of Fu Brick tea could improve forecasting of tea quality and improve the manufacturing process.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://osf.io/2cbq4, https://osf.io/2cbq4.
Author Contributions
MX and XL designed the experiment. MX, JC, and WC took the tea samples. WZ and XZ isolated microbes from the samples. WC and HM analyzed the data. XL and MX lead the writing of the manuscript. All authors revised and give comments to the manuscripts. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.892437/full#supplementary-material
References
Abarenkov, K., Henrik Nilsson, R., Larsson, K., Alexander, I., Eberhardt, U., Erland, S., et al. (2012). The UNITE database for molecular identification of fungi–recent updates and future perspectives. N. Phytol. 186, 281–285. doi: 10.1111/j.1469-8137.2009.03160.x
Angayarkanni, J., Palaniswamy, M., Murugesan, S., and Swaminathan, K. (2002). Improvement of tea leaves fermentation with Aspergillus spp. pectinase. J. Biosci. Bioeng. 94, 299–303. doi: 10.1263/jbb.94.299
Benson, D. A., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J., and Wheeler, D. L. (2005). GenBank. Nucleic Acids Res. 36, 25–30. doi: 10.1093/nar/gkj157
Boekhout, T., Fell, J. W., and O’Donnell, K. (1995). Molecular systematics of some yeast-like anamorphs belonging to the Ustilaginales and Tilletiales. Stud. Mycol. 38, 175–183.
Bourque, D., Ouellette, B., Andre, G., and Groleau, D. (1992). Production of poly-β-hydroxybutyrate from methanol: characterization of a new isolate of Methylobacterium extorquens. Appl. Microbiol. Biotechnol. 37, 7–12.
Cao, N., Xia, Y., Gong, C., and Tsao, G. T. (1997). “Production of 2, 3-butanediol from pretreated corn cob by Klebsiella oxytoca in the presence of fungal cellulase,” in Biotechnology for Fuels and Chemicals, ed J. Sachs (Humana Press), 129–139. doi: 10.1007/BF02920419
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chen, A. J., Hubka, V., Frisvad, J. C., Visagie, C. M., Houbraken, J., Meijer, M., et al. (2017). Polyphasic taxonomy of Aspergillus section Aspergillus (formerly Eurotium), and its occurrence in indoor environments and food. Stud. Mycol. 88, 37–135. doi: 10.1016/j.simyco.2017.07.001
Chen, C., Liu, Y., Tian, H., Ai, L., and Yu, H. (2020). Metagenomic analysis reveals the impact of JIUYAO microbial diversity on fermentation and the volatile profile of Shaoxing-jiu. Food Microbiol. 86, 103326. doi: 10.1016/j.fm.2019.103326
Chen, G., Chen, C., and Lei, Z. (2017). Meta-omics insights in the microbial community profiling and functional characterization of fermented foods. Trends Food Sci. Technol. 65, 23–31. doi: 10.1016/j.tifs.2017.05.002
Chen, G., Xie, M., Dai, Z., Wan, P., Ye, H., Zeng, X., et al. (2018). Kudingcha and fuzhuan brick tea prevent obesity and modulate gut microbiota in high-fat diet fed mice. Mol. Nutr. Food Res. 62, 1700485. doi: 10.1002/mnfr.201700485
Chiang, C. T., Weng, M. S., Lin-Shiau, S. Y., Kuo, K. L., Tsai, Y. J., and Lin, J. K. (2006). Pu-erh tea supplementation suppresses fatty acid synthase expression in the rat liver through downregulating Akt and JNK signalings as demonstrated in human hepatoma HepG2 cells. Oncol. Res. 16, 119–128. doi: 10.3727/000000006783981143
del Carmen Portillo, M., and Mas, A. (2016). Analysis of microbial diversity and dynamics during wine fermentation of Grenache grape variety by high-throughput barcoding sequencing. LWT Food Sci. Technol. 72, 317–321. doi: 10.1016/j.lwt.2016.05.009
DeLong, E. F., Preston, C. M., Mincer, T., Rich, V., Hallam, S. J., Frigaard, N., et al. (2006). Community genomics among stratified microbial assemblages in the ocean's interior. Science 311, 496–503. doi: 10.1126/science.1120250
Doran, J. B., Aldrich, H. C., and Ingram, L. O. (1994). Saccharification and fermentation of sugar cane bagasse by Klebsiella oxytoca P2 containing chromosomally integrated genes encoding the Zymomonas mobilis ethanol pathway. Biotechnol. Bioeng. 44, 240–247. doi: 10.1002/bit.260440213
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Fu, D., Ryan, E. P., Huang, J., Liu, Z., Weir, T. L., Snook, R. L., et al. (2011). Fermented Camellia sinensis, Fu Zhuan Tea, regulates hyperlipidemia and transcription factors involved in lipid catabolism. Food Res. Int. 44, 2999–3005. doi: 10.1016/j.foodres.2011.07.008
Gao, J., Gu, F., He, J., Xiao, J., Chen, Q., Ruan, H., et al. (2013). Metagenome analysis of bacterial diversity in Tibetan kefir grains. Eur. Food Res. Technol. 236, 549–556. doi: 10.1007/s00217-013-1912-2
Ghatnekar, M. S., Pai, J. S., and Ganesh, M. (2002). Production and recovery of poly-3-hydroxybutyrate from Methylobacterium sp V49. J. Chem. Technol. Biol. 77, 444–448. doi: 10.1002/jctb.570
Golias, H., Dumsday, G. J., Stanley, G. A., and Pamment, N. B. (2002). Evaluation of a recombinant Klebsiella oxytoca strain for ethanol production from cellulose by simultaneous saccharification and fermentation: comparison with native cellobiose-utilising yeast strains and performance in co-culture with thermotolerant yeast and Zymomonas mobilis. J. Biotechnol. 96, 155–168. doi: 10.1016/s0168-1656(02)00026-3
Jiang, C., Zeng, Z., Huang, Y., and Zhang, X. (2018). Chemical compositions of Pu'er tea fermented by Eurotium Cristatum and their lipid-lowering activity. LWT Food Sci. Technol. 98, 204–211. doi: 10.1016/j.lwt.2018.08.007
Jie, G. L., Lin, Z., Zhang, L., Lv, H., He, P., and Zhao, B. (2006). Free radical scavenging effect of Pu-erh tea extracts and their protective effect on oxidative damage in human fibroblast cells. J. Agric. Food Chem. 54, 8058–8064. doi: 10.1021/jf061663o
Kang, D., Su, M., Duan, Y., and Huang, Y. (2019). Eurotium cristatum, a new fungi probiotic from Fuzhuan brick tea, alleviated obesity in mice by modulating gut microbiota. Food Funct. 10, 5032–5045. doi: 10.1039/c9fo00604d
Kim, B. C., Jeong, W. J., Oh, H. W., Kim, H., Park, D., Park, H., et al. (2009). Paenibacillus pueri sp. nov., isolated from Pu'er tea. Int. J. Syst. Evol. Microbiol. 59, 1002–1006. doi: 10.1099/ijs.0.002352-0
Lane, D. J. (1991). “V4/23S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt and M. Goodfellow (Chichester: John Wiley and Sons), 115–175.
Li, Q., Chai, S., Li, Y., Huang, J., Luo, Y., Xiao, L., et al. (2018). Biochemical components associated with microbial community shift during the pile-fermentation of primary dark tea. Front. Microbiol. 9, 1509. doi: 10.3389/fmicb.2018.01509
Li, Q., Huang, J., Li, Y., Zhang, Y., Luo, Y., Chen, Y., et al. (2017). Fungal community succession and major components change during manufactural process of Fu Brick tea. Sci. Rep. 7, 1–9. doi: 10.1038/s41598-017-07098-8
Li, Q., Li, Y., Luo, Y., Zhang, Y., Chen, Y., Lin, H., et al. (2019). Shifts in diversity and function of the bacterial community during the manufacture of Fu Brick tea. Food Microbiol. 80, 70–76. doi: 10.1016/j.fm.2019.01.001
Li, Q., Liu, Z., Huang, J., Luo, G., Liang, Q., Wang, D., et al. (2013). Anti-obesity and hypolipidemic effects of Fuzhuan brick tea water extract in high-fat diet-induced obese rats. J. Sci. Food Agric. 93, 1310–1316. doi: 10.1002/jsfa.5887
Liu, C., Hsu, W. H., Lee, F. L., and Liao, C. C. (1996). The isolation and identification of microbes from a fermented tea beverage, Haipao, and their interactions during Haipao fermentation. Food Microbiol. 13, 407–415. doi: 10.1006/fmic.1996.0047
Liu, Z., Gao, L., Chen, Z., Zeng, X., Huang, J., Gong, Y., et al. (2019). Leading progress on genomics, health benefits and utilization of tea resources in China. Nature 566, 7742. Available online at: https://www.nature.com/articles/d42473-019-00032-8
Ma, Y., Duan, S., Zhang, D., Su, X., Zhang, D., Lv, C., et al. (2017). Microbial succession and the dynamics of chemical compounds during the solid-state fermentation of pu-erh tea. Appl. Sci. 7, 166. doi: 10.3390/app7020166
Mohammad, F. G., and Alireza, T. (2007). Isolation and characterization of polyphenol oxidase-and peroxidase-producing Bacillus strains from fully fermented tea (Camellia sinensis). World J. Microb. Biotechnol. 23, 1327–1332. doi: 10.1007/S11274-007-9363-1
Moon, S., Kim, K., Park, J. M., Min, J., and Song, H. (2018). Development of a semi-continuous two-stage simultaneous saccharification and fermentation process for enhanced 2,3-butanediol production by Klebsiella oxytoca. Lett. Appl. Microbiol. 66, 300–305. doi: 10.1111/lam.12845
Oh, H., Kim, B., Lee, K. H., Kim, D. Y., Park, D., Park, H., et al. (2008). Paenibacillus camelliae sp. nov., isolated from fermented leaves of Camellia sinensis. J. Microbiol. 46, 530–534. doi: 10.1007/s12275-008-0233-9
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'hara, R. B., et al. (2006). Vegan: Community Ecology Package. R Package Version 2. 3–3.
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596. doi: 10.1093/nar/gks, 1219
R Development Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rehman, H. U., Siddique, N. N., Aman, A., Nawaz, M. A., Baloch, A. H., and Qader, S. A. U. (2015). Morphological and molecular based identification of pectinase producing Bacillus licheniformis from rotten vegetable. J. Genet. Eng. Biotechnol. 13, 139–144. doi: 10.1016/j.jgeb.2015.07.004
Rui, Y., Wan, P., Chen, G., Xie, M., Sun, Y., Zeng, X., et al. (2019). Analysis of bacterial and fungal communities by Illumina MiSeq platforms and characterization of Aspergillus cristatus in Fuzhuan brick tea. LWT Food Sci. Technol. 110, 168–174. doi: 10.1016/j.lwt.2019.04.092
Rutiaga-Quiñones, O. M., Córdova, É., Martell-Nevárez, M. A., Guillamón, J. M., Rozes, N., and Páez, J. (2012). Volatile compound production in Agave duranguensis juice fermentations using four native yeasts and NH4Cl supplementation. Eur. Food Res. Technol. 235, 29–35. doi: 10.1007/s00217-012-1729-4
Serra, J. L., Moura, F. G., de Melo Pereira, G. V., Soccol, C. R., Rogez, H., and Darnet, S. (2019). Determination of the microbial community in Amazonian cocoa bean fermentation by Illumina-based metagenomic sequencing. LWT Food Sci. Technol. 106, 229–239. doi: 10.1016/j.lwt.2019.02.038
Song, Z., Du, H., Zhang, Y., and Xu, Y. (2017). Unraveling core functional microbiota in traditional solid-state fermentation by high-throughput amplicons and metatranscriptomics sequencing. Front. Microbiol. 8, 1294. doi: 10.3389/fmicb.2017.01294
Sun, L. H., Wang, X. D., Dai, J. Y., and Xiu, Z. (2009). Microbial production of 2, 3-butanediol from Jerusalem artichoke tubers by Klebsiella pneumoniae. Appl. Microbiol. Biotechnol. 82, 847–852. doi: 10.1007/s00253-008-1823-5
Tofalo, R., Fusco, V., Böhnlein, C., Kabisch, J., Logrieco, A. F., Habermann, D., et al. (2019). The life and times of yeasts in traditional food fermentations. Crit. Rev. Food Sci. Nutr. 60, 3103–3132. doi: 10.1080/10408398.2019.1677553
Uzuner, S., and Cekmecelioglu, D. (2015). Enhanced pectinase production by optimizing fermentation conditions of Bacillus subtilis growing on hazelnut shell hydrolyzate. J. Mol. Catal. 113, 62–67. doi: 10.1016/j.molcatb.2015.01.003
Wang, X., Du, G., Chen, H., Zeng, X., Liu, B., Guo, C., et al. (2021). Comparative metagenomics reveals microbial communities and their associated functions in two types of Fuzhuan brick tea. Front. Microbiol. 12, 705681. doi: 10.3389/fmicb.2021.705681
Weiss, S., Xu, Z. Z., Peddada, S., Amir, A., Bittnger, K., Gonzalez, A., et al. (2017). Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome 5, 27. doi: 10.1186/s40168-017-0237-y
Wen, Q., and Liu, S. (1991). Evolutionary regulation of dominant fungi in fuzhuan brick tea during the fungus growing process. J. Tea Sci. 11, 56–62.
White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols: A Guide to Methods and Applications, eds M. A., Gelfand, D. H., Sninsky, J. J. Innis, and T. J. White (San Diego, CA: Academic Press), 315–322.
Wu, Y. Y., Ding, L., Xia, H. L., and Tu, Y. Y. (2010). Analysis of the major chemical compositions in Fuzhuan brick-tea and its effect on activities of pancreatic enzymes in vitro. Afr. J. Biotechnol. 9, 6748–6754. doi: 10.5897/AJB10.947
Wulandari, R. A., Amano, M., Yanagita, T., Tanaka, T., Kouno, I., Kawamura, D., et al. (2011). New phenolic compounds from Camellia sinensis L. leaves fermented with Aspergillus sp. J. Nat. Med. 65, 594–597. doi: 10.1007/s11418-011-0515-0
Xia, F., Hu, S., Zheng, X., Wang, M. W., Zhang, C. C., Wu, Z. N., et al. (2022). New insights into metabolomics profile generation in fermented tea: the relevance of bacteria and metabolites in Fuzhuan brick tea. J. Sci. Food Agric. 102, 350–359. doi: 10.1002/jsfa.11365
Xiao, Y., Xiong, T., Peng, Z., Liu, C., Huang, T., Yu, H., et al. (2018). Correlation between microbiota and flavours in fermentation of Chinese Sichuan Paocai. Food Res. Int. 114, 123–132. doi: 10.1016/j.foodres.2018.06.051
Xu, A., Wang, Y., Wen, J., Liu, P., Liu, Z., and Li, Z. (2011). Fungal community associated with fermentation and storage of Fuzhuan brick-tea. Int. J. Food Microbiol. 146, 14–22. doi: 10.1016/j.ijfoodmicro.2011.01.024
Yao, Y., Wu, M., Huang, Y., Li, C., Pan, X., Zhu, W., et al. (2017). Appropriately raising fermentation temperature beneficial to the increase of antioxidant activity and gallic acid content in Eurotium cristatum-fermented loose tea. LWT Food Sci. Technol. 82, 248–254. doi: 10.1016/j.lwt.2017.04.032
Yellore, V. S., Ghatnekar, M. S., Pai, J. S., and Desai, A. J. (2000). Effect of total nutrient feed on production of poly-3-hydroxybutyrate by Methylobacterium sp. ZP24 grown on sugars. J. Ind. Microbiol. Biotechnol. 25, 276–279. doi: 10.1038/sj.jim.7000073
Yu, E. K., Deschatelets, L., Louis-Seize, G., and Saddler, J. N. (1985). Butanediol production from cellulose and hemicellulose by Klebsiella pneumoniae grown in sequential coculture with Trichoderma harzianum. Appl. Environ. Microbiol. 50, 924–929. doi: 10.1128/aem.50.4.924-929.1985
Zhang, L., Zhang, Z., Zhou, Y., Ling, T., and Wan, X. (2013). Chinese dark teas: postfermentation, chemistry and biological activities. Food Res. Int. 53, 600–607. doi: 10.1016/j.foodres.2013.01.016
Zhang, Y. J., Zhang, S., Liu, X. Z., Wen, H. A., and Wang, M. (2010). A simple method of genomic DNA extraction suitable for analysis of bulk fungal strains. Lett. Appl. Microbiol. 51, 114–118. doi: 10.1111/j.1472-765X.2010.02867.x
Zhao, M., Zhang, D., Su, X., Duan, S., Wan, J., Yuan, W., et al. (2015). An integrated metagenomics/metaproteomics investigation of the microbial communities and enzymes in solid-state fermentation of Pu-erh tea. Sci. Rep. 5, 10117. doi: 10.1038/srep10117
Zheng, W., Wan, X., and Bao, G. (2015). Brick dark tea: a review of the manufacture, chemical constituents and bioconversion of the major chemical components during fermentation. Phytochem. Rev. 14, 499–523. doi: 10.1007/s11101-015-9402-8
Zhou, S., and Ingram, L. O. (2001). Simultaneous saccharification and fermentation of amorphous cellulose to ethanol by recombinant Klebsiella oxytoca SZ21 without supplemental cellulase. Biotechnol. Lett. 23, 1455–1462. doi: 10.1023/A:1011623509335
Keywords: Aspergillus cristatus, Fu Brick tea, pile fermentation, golden flower fungus, microbiota
Citation: Xiang M, Chu J, Cai W, Ma H, Zhu W, Zhang X, Ren J, Xiao L, Liu D and Liu X (2022) Microbial Succession and Interactions During the Manufacture of Fu Brick Tea. Front. Microbiol. 13:892437. doi: 10.3389/fmicb.2022.892437
Received: 09 March 2022; Accepted: 23 May 2022;
Published: 23 June 2022.
Edited by:
Abd El-Latif Hesham, Assiut University, EgyptReviewed by:
Ajar Nath Yadav, Eternal University, IndiaBeatrix Stessl, University of Veterinary Medicine Vienna, Austria
Geromy G. Moore, Southern Regional Research Center (USDA-ARS), United States
Copyright © 2022 Xiang, Chu, Cai, Ma, Zhu, Zhang, Ren, Xiao, Liu and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Meichun Xiang, xiangmc@im.ac.cn; Xingzhong Liu, liuxz@nankai.edu.cn