Yu Mu
Wei Su
Yingchun Mu
Li Jiang- 1School of Liquor and Food Engineering, Guizhou University, Guiyang, China
- 2Guizhou Provincial Key Laboratory of Agricultural and Animal Products Storage and Processing, Guizhou University, Guiyang, China
Panxian ham, a traditional Chinese dry-cured ham, is protected by national geographical indication. Similar to other fermented foods, the microbial population of dry-cured ham is pivotal to taste and flavor formation. This study aimed to establish the relationship between microorganisms and metabolites during the spontaneous fermentation of Panxian ham. Multivariate analysis based on metabolomics data revealed that continuous metabolic changes occurred during the entire fermentation process, with the most significant changes occurring in the initial stage of ripening. Thirty-one significantly different metabolites (SDMs) were identified as discriminant factor, and pathway analysis suggested that these metabolites were involved in 30 pathways, including alanine, aspartate, and glutamate metabolism; glycine, serine, and threonine metabolism; and arginine and proline metabolism. Microbial community analysis using the Illumina MiSeq platform indicated that the bacterial community was more complex than the fungal community, and their succession regulation differed during processing. At the genus level, 11 bacteria and five fungi were identified as core microbes, of which Staphylococcus was the dominant bacteria and Debaryomyces and Aspergillus were the dominant fungi. Further, statistical redundancy analysis (RDA) indicated that Staphylococcus, Debaryomyces, and Chromohalobacter promoted the production of amino and fatty acids; Cobetia and Aspergillus were associated with sugar metabolism, and Kushneria, Penicillium, and Yamadazyma were closely related with organic acids. These findings provide fundamental knowledge regarding the metabolically active microorganisms in Panxian ham, helping industrial processors to develop effective strategies for standardizing quality parameters.
Introduction
Dry-cured ham is a traditional fermented meat product with a long history of consumption that is very popular around the world for its unique sensory properties and high nutritive value. Previous studies have suggested that dry-cured ham may be included as a conventional part of a healthy diet (Fernández et al., 2007; Jimenez-Colmenero et al., 2010), because many nutritional and functional substances are accumulated during fermentation, such as the essential amino acids, taurine, cysteine, carnosine, creatinine, anserine, and glutathione (Marusic et al., 2013). In addition, dry-cured ham is also a valuable source of bioactive peptides. Recent reports indicate that regular consumption of dry-cured ham can provide protective health benefits, such as reducing the risk of cardiovascular disease and inhibiting the activation of platelets and monocytes (Martinez-Sanchez et al., 2017; Montoro-Garcia et al., 2017).
Panxian ham, one of the famous dry-cured hams in China, is mainly distributed in Guizhou and its neighboring provinces, and is protected as a national geographical indication product (following Jinhua and Xuanwei hams) due to its strong characteristic taste, flavor, aroma, and texture; which are attributable to the Panxian region. Similar to most dry-cured hams, the production of Panxian ham usually consists of raw material selection, salting, resting, and drying-ripening. The entire process is usually accompanied by a series of complicated biochemical and enzymatic reactions, such as protein degradation, lipid oxidation, glycolysis, and Maillard reactions (Zhou and Zhao, 2007; Petrova et al., 2015). It is well known that endogenous enzymes play essential in these reactions (Toldrá and Flores, 2000), but the enzymatic action of microorganisms cannot be discounted (Simoncini et al., 2007; Landeta et al., 2013). The diversity of fermented meat’s microbial community could be influenced by geographical location and production environment due to the traditional spontaneous fermentation that occurs during processing (Van Reckem et al., 2019). This is especially important because microbial diversity has the potential to affect the metabolome, and thus the quality and typicity of fermented foods (Montanari et al., 2018). However, previous studies on dry-cured ham have focused on the relationship between microorganisms and volatile flavors (Ge et al., 2017; Mu et al., 2019), while the effect of microorganisms on non-volatile metabolites has not yet been studied.
In recent years, multivariate analysis of microorganisms and metabolites has been applied to improving the quality of fermented foods (Park et al., 2019; Wurihan et al., 2019). Thus, exploration of microbial community succession, metabolic changes, and their interrelationship during processing, is warranted to improve Panxian ham quality. Previous studies, based on culture-dependent methods, have investigated changes in microbial diversity during ham production (Tu et al., 2010). In fact, the majority of microorganisms cannot be detected by culture-dependent methods due to their “viable but non-culturable” state (Fakruddin et al., 2013). High-throughput sequencing (HTS) technology allows a deeper and more precise evaluation of complex microbiota with lower cost and shorter analysis time (Tamang et al., 2016), and has successfully characterized microbial community succession during the spontaneous fermentation of sausage (Wang et al., 2019). Moreover, nuclear magnetic resonance (NMR)-based metabolomics analysis has revealed metabolic changes during dry-cured ham processing (Zhang et al., 2019), but has lower sensitivity for detecting low abundance metabolites compared with the MS-based technique (Sugimoto et al., 2017).
Indigenous microorganisms are known to play a pivotal role in the final taste, flavor, and texture of fermented foods. However, little information is available about the functional microorganisms in Panxian ham; this frequently leads to fluctuations in quality among different batches. Hence, a gas chromatography time of flight mass spectrometry (GC-TOF-MS)-based metabolomics approach and HTS technology were applied to monitor the evolution of the Panxian ham metabolite profile and microbial community structure during processing. Furthermore, multivariate redundancy analysis (RDA) was conducted to uncover the correlation between dominant microbes and significantly different metabolites (SDMs). It is hoped that this study will aid in the standardization of Panxian ham production throughout the region.
Materials and Methods
Processing of Dry-Cured Hams
Thirty fresh pork thighs weighing between 10–12 kg from a local slaughterhouse were selected and processed in a Panxian ham food factory (Panzhou, China) according to traditional methods, including aging, salting, resting, and ripening. The raw hams were aged for 24 h and then salted with sea salt in 4°C cold storage. The curing process occurred in three stages over 21 days; the first use of salt was 4% of the weight of the thigh, and the second and third uses were 2% each, for a total amount of 8% salt based on the weight of the thigh. After salting, the salt and blood remaining on the surface of the ham was rinsed. Then, the hams were transferred to a room for 30 days of rest at 10–15°C and 75–80% relative humidity (RH). Finally, the hams were hung in a fermentation room at 20–22°C and 65–70% RH until fully ripened.
Sampling
The biceps femoris muscles were collected at six different stages of ham processing: raw ham, post-salting (21 days), post-resting (51 days), and initial stage (180 days), middle stage (360 days), and final stage (540 days) of ripening. Five hams were sampled randomly at each stage during processing for parallel testing. After removal of subcutaneous fat and connective tissue, the samples were vacuum-packaged and stored at −80°C until analysis.
Metabolites Extraction and GC-TOF-MS Analysis
Ham sample (300 ± 1 mg) was placed in the 5 mL EP tubes with 3 mL extraction liquid (VMethanol: VChloroform = 3:1), then 5 μL of L-2-Chlorophenylalanine (1 mg/mL stock in dH2O) was added as internal standard followed by vortex-mixing for 30 s. The mixture was homogenized in a ball mill for 4 min at 45 Hz and treated with ultrasound for 5 min (incubated in ice water). Next, the centrifugal process was conducted at 12,000 rpm and 4°C for 15 min, and 50 μL of sample supernatant was transferred into a fresh 1.5 mL EP tubes. In addition, quality control (QC) samples were prepared by pooling 20 μL of supernatant from each sample. Then the supernatant was completely dried in a vacuum concentrator (without heating) and 80 μL of methoxyamination hydrochloride in pyridine solution (20 mg/mL) was added and incubated for 30 min at 80°C. Then derivatized using 40 μL of the BSTFA (Bistrifluoroacetamide) [1% TMCS (Templated mesoporous carbons), v/v] at 70°C for 1.5 h. Finally, 5 μL of FAMEs (Fatty acid methyl esters) in chloroform was added to the QC sample when cooling to the room temperature. All samples were analyzed by gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS).
GC-TOF-MS analysis was performed using an Agilent 7890 gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer. The system utilized a DB-5MS capillary column coated with 5% diphenyl cross-linked with 95% dimethylpolysiloxane (30 m × 250 μm inner diameter, 0.25 μm film thickness; J&W Scientific, Folsom, CA, United States). A 1 μL aliquot of the analyte was injected in splitless mode. Helium was used as the carrier gas, the front inlet purge flow was 3 mL/min, and the gas flow rate through the column was 1 mL/min. The initial temperature was kept at 50°C for 1 min, then raised to 310°C at a rate of 10°C/min, then kept for 8 min at 310°C. The injection, transfer line, and ion source temperatures were 280, 280, and 250°C, respectively. The energy was −70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with the m/z range of 50–500 at a rate of 12.5 spectra per second after a solvent delay of 6.33 min.
Metabolome Data Processing and Analysis
Chroma TOF 4.3X software of LECO Corporation and LECO-Fiehn Rtx5 database were used for raw peaks exacting, the data baselines filtering and calibration of the baseline, peak alignment, deconvolution analysis, peak identification and integration of the peak area (Kind et al., 2009). Both of mass spectrum match and retention index match were considered in metabolites identification. Remove peaks detected in <50% of QC samples or RSD > 30% in QC samples (Dunn et al., 2011). The data including the compound name, peak area, retention time, similarity to metabolites in the database, and mass (m/z) were ultimately imported into Microsoft Excel.
The compound name and peak area information were selected as variables for multivariate statistical analysis using SIMCA-P 14.1 (Umetrics, Umea, Sweden). The GC-TOF-MS data sets were visualized with unsupervised principal component analysis (PCA) and supervised orthogonal partial least squares discriminant analysis (OPLS-DA). Meanwhile, the cross-validation analysis of variance (CV-ANOVA) was performed to evaluate the OPLS-DA model. Based on OPLS-DA model and one-way analysis of variance (ANOVA), the variable importance in projection (VIP) values and p-value were calculated and the SDMs were determined. The peak area of SDMs was visualized via R software (v3.6.1) with the “pheatmap” package after log-transformation. Pathway and enrichment analyses were conducted by MetaboAnalyst 4.01, which is a free and web-based tool with a high-quality KEGG (Kyoto encyclopedia of genes and genomes) metabolic pathway.
DNA Extraction and PCR Amplification
Total bacterial DNA were extracted from samples using the Power Soil DNA Isolation Kit (MO BIO Laboratories) according to the manufacturer’s protocol. DNA quality and quantity were assessed by the ratios of 260/280 nm and 260/230 nm. Then DNA was stored at −80°C until further processing. The V3-V4 region of the bacteria 16S rRNA gene was amplified with primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). The ITS1 region of the fungi was amplified with the forward primer ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and the reverse primer ITS1R (5′-GCTGCGTTCTTCATCGATGC-3′).
PCR amplification was performed in a total volume of 50 μL, which contained 10 μL Buffer, 0.2 μL Q5 High-Fidelity DNA Polymerase, 10 μL High GC Enhancer, 1 μL dNTP, 10 μM of each primer and 60 ng genome DNA. Thermal cycling conditions were as follows: an initial denaturation at 95°C for 5 min, followed by 15 cycles at 95°C for 1 min, 50°C for 1 min and 72°C for 1 min, with a final extension at 72°C for 7 min. The PCR products from the first step PCR were purified through VAHTSTM DNA Clean Beads. A second round PCR was then performed in a 40 μL reaction which contained 20 μL 2 × Phusion HF MM, 8 μL ddH2O, 10 μM of each primer and 10 μL PCR products from the first step. Thermal cycling conditions were as follows: an initial denaturation at 98°C for 30 s, followed by 10 cycles at 98°C for 10 s, 65°C for 30 s and 72°C for 30 s, with a final extension at 72°C for 5 min. Finally, all PCR products were quantified by Nanodrop 2000 and pooled together. HTS analysis of bacterial rRNA and fungal ITS1 genes was performed on the purified, pooled sample using the Illumina Hiseq 2500 platform (2 × 250 paired ends).
Microbiome Data Processing and Bioinformatics Analysis
Raw fastq files were demultiplexed, quality-filtered by Trimmomatic and merged by FLASH (Magoc and Salzberg, 2011). Operational taxonomic units (OTUs) were clustered with 97% similarity cutoff using UPARSE (Bokulich et al., 2013). Then, the taxonomy of each sequence was analyzed by Ribosomal Database Project (RDP) Classifier against the Silva (SSU123) database using confidence threshold of 80% (Fadrosh et al., 2014). Finally, chimeric sequences were identified and removed using UCHIME. The alpha diversity indices, rarefaction curves, and the number of unique OTUs in each sample were evaluated using the Mothur software (Kozich et al., 2013). For beta diversity analysis, principal coordinate analysis (PCoA) and non-metric multidimensional scaling (NMDS) were performed using QIIME software (Caporaso et al., 2010). Prediction of gene function was performed using PICRUSt (phylogenetic investigation of communities by reconstruction of unobserved states) (Langille et al., 2013). To establish the co-occurrence network between core bacterial and fungal genus, we used R software with the “corrplot” package for calculating Pearson’s correlation coefficient and applied Cytoscape software for visualization.
Integrative Analysis of Metabolome and Microbiome Data
The relative abundance of core microbes (≥1%) and the relative content of SDMs (VIP > 1 and p < 0.05) were used for the integrative analysis. The multivariate RDA based on the “vegan” package was performed by R software.
Statistical Analysis
Statistical analysis was performed using SPSS 22.0 software. The data obtained were subjected to one-way analysis of variance (ANOVA), and the significance level in the analyses was considered at p < 0.05.
Results and Discussion
Metabolite Profile Comparison of Ham Samples at Different Stages
Typical GC-TOF-MS total ion chromatograms (TICs) of dry-cured ham samples from the six different stages are shown in Supplementary Figure S1. In total, 578 signal features were identified in all samples and 527 remained after applying an imputation model with an 80% limit for missing data. PCA of the entire set of measured analytes was used to investigate gross changes in metabolic physiology of Panxian ham during processing. The PCA score plot showed 40.5% and 15.6% variances by PC1 and PC2, respectively (Figure 1A). As fermentation progressed, a typical metabolite drift was observed. Specifically, metabolite profiles were distributed along PC1 and progressed from quadrants 2 to 4, which illustrated that continuous metabolic changes occurred during processing. In the PCA score plot, a clear separation between post-resting and initial stage of ripening was found, indicating that extremely rapid metabolite changes transpired in this period. To maximize the distinction between different samples, and to obtain a better understanding of the metabolites responsible for the separation, the OPLS-DA model was applied. Compared with the PCA model, the R2 and Q2 values for evaluation of model quality increased, and all samples in the OPLS-DA score plot were within the 95% confidence interval (Figure 1B). Meanwhile, the CV-ANOVA p-value for the OPLS-DA model generated from the original data matrix was 1.67 × 10–22, indicating this model’s improved predictive ability (Bradley and Robert, 2016).
Figure 1. Principal component analysis (PCA) (A) and orthogonal partial least squares discriminant analysis (OPLS-DA) (B) score plots derived from the GC-TOF-MS data from Panxian ham and Quality Control (QC) samples. A–F indicate the traditional spontaneous fermentation process stage of Panxian ham: (A) raw ham; (B) post-salting; (C) post-resting; (D) initial stage of ripening; (E) middle stage of ripening; (F) final stage of ripening. (QC) represents the quality control sample.
Identification of Significantly Different Metabolites and Analysis of Metabolic Pathways During Processing
A total of 226 metabolites were identified by mass spectrum matching and 80 reliable metabolites were determined based on spectral similarity values >700 (Hao et al., 2019). The differential variables among groups were selected based on VIP scores >1.0 from OPLS-DA and p < 0.05 from ANOVA. In total, 31 SDMs were identified, including 15 amino acids, 6 fatty acids, 3 organic acids, 2 sugars, 2 polyols, 2 nucleic acids, and 1 other metabolite (Supplementary Table S1). To identify the most relevant metabolic pathways during Panxian ham processing, all SDMs were evaluated using metabolic pathway and enrichment analysis with MetaboAnalyst 4.0. Thirty pathways were obtained and “alanine, aspartate, and glutamate metabolism,” “glycine, serine, and threonine metabolism,” and “arginine and proline metabolism” were determined as the most relevant metabolic pathways (Supplementary Figure S2). Figure 2 shows the integrated metabolic pathway adapted from the KEGG database2 and previous studies (Toldrá and Flores, 2000), including protein, fat, and glycogen metabolism.
Figure 2. The integrated metabolic pathway shows the dynamics of significantly different metabolites (SDMs) during Panxian ham processing. Beneath each metabolite, the color gradient and their values indicate the log2 (fold change) with respect to raw ham: red and blue represent up and downregulation of metabolites, respectively. Glycine, serine, and threonine metabolism, alanine, aspartate, and glutamate metabolism, and arginine and proline metabolism are shown in the blue, yellow, and green areas, respectively.
Amino acids, often the most important primary metabolites in dry-cured ham because of the high protein content of pork thigh, had a crucial role on the formation of taste in dry-cured hams (Zhang et al., 2018; Liu et al., 2019). Specifically, alanine, valine, and serine contributed to perceived sweetness, glutamate, pyroglutamate, and aspartate were responsible for the umami taste, while tyrosine and lysine were associated with the unique aged taste (Careri et al., 1993; Zheng et al., 2015). In addition, amino acids were the precursors of many volatile flavors and played an indispensable role in the flavor formation of dry-cured ham (Shi et al., 2019). In the present study, the levels of all amino acids increased rapidly until the initial stage of ripening. The concentrations of most amino acids such as serine, glutamate, alanine, L-allothreonine, and methionine then stabilized, whereas those of phenylalanine, ornithine, tyrosine, and lysine continued to increase until the final stage of the ripening. One possible interpretation of this result is that the most significant proteolysis usually occurred in the first 5 months of dry-cured ham maturation (Grossi et al., 2014). Rapid lipolysis of dry-cured ham also occurred during early maturation (Motilva et al., 1993). As shown in Figure 2, the levels of myristic acid, palmitoleic acid, heptadecanoic acid, arachidonic acid, and arachidic acid increased sharply during this period. As the second most abundant SDM in Panxian ham, fatty acids were also closely related to the flavor of dry-cured ham (Shi et al., 2019), and hexanal from the oxidation of n-6 fatty acids like arachidonic acid was an odor-active compound in Istrian and Iberian hams (Marušić et al., 2014).
Glucose, fructose, glycerol, and myo-inositol were identified as the major carbon compounds in Panxian ham. Among them, the concentrations of glucose and fructose declined continually throughout the process, while those of myo-inositol and glycerol increased steadily. Glucose and fructose are often directly associated with sweetness; however, some studies have reported that sugars had little contribution to the taste of dry-cured ham due to its low content in the final product (Zhang et al., 2018). Glycerol and myo-inositol also had a slight effect on perceived sweetness (Seo et al., 2016). Glycerol was a product of fat degradation. Therefore, the increase in glycerol content was linked to lipolytic enzyme activity (Koutina et al., 2012). Myo-inositol was a vitamin-like essential nutrient. Yang et al. (2018) found that higher concentrations of myo-inositol could improve digestion and growth.
Organic acids were also important taste-related compounds in foods. Succinic acid has been identified as a taste-active compound in boneless dry-cured ham (Zhang et al., 2019), and the combination of succinic, lactic, and citric acids can generate unique sour and umami tastes (Dashdorj et al., 2015). Malate was a mild acid and contributed to the red coloration and protein degradation of cured meat (Kim et al., 2019). Sugimoto et al. (2017) found that dry-cured ham with a longer ripening period had a much higher malate content, this corroborated with the present study (Figure 2). Creatine was linked to bitter taste and has been described as a taste-active compound in duck meat (Liu et al., 2013). In general, a decrease in creatine during the initial stage of dry-cured ham ripening was accompanied by an increase in creatinine (Mora et al., 2010). Therefore, the creatinine/creatine ratio was considered useful in estimating the minimum ripening time. However, Escudero et al. (2011) found that the hypoxanthine/inosine ratio was more suitable as a potential biomarker for the primary ripening of dry-cured ham. In addition, the high levels of hypoxanthine can not only enhance the cured meat taste (Ichimura et al., 2017), but also contribute to the taste of boneless dry-cured hams (Zhang et al., 2019).
Abundance and Diversity of Bacterial and Fungal Microorganisms
High-throughput sequencing generated a total of 1,991,084 V3-V4 16S rRNA regions and 2,158,470 ITS1 high-quality sequence reads from 30 Panxian ham samples. The numbers of effective sequences ranged from 62,788 to 71,736 for bacteria and from 68,114 to 75,068 for fungi (Supplementary Table S2). As shown in Figure 3A, the OTU numbers of bacteria and fungi decreased continually during Panxian ham processing, and the numbers of bacteria always higher than that of fungi. The rarefaction curve of each sample based on the OTU number was used to evaluate sequencing depth, which showed that sequencing data was sufficient for subsequent analyses (Figure 3B). Dynamic changes to α-diversity indices indicated that the fungal diversity of Panxian ham fluctuated greatly during processing, while bacterial diversity tended to continually decrease (Table 1). Similar trends have been reported during the Suan yu, a fermentation process for Chinese traditional fermented fish; Zang et al. (2018) suggested that chitin (a characteristic component of the fungal cell wall) may be responsible for these results, providing a certain protection against environmental stress. Microbiota structural changes were analyzed using NMDS and PCoA based on Binary-Jaccard distances to evaluate variation and similarity (Figures 3C,D). The results demonstrated that the succession of bacteria and fungi was not synchronous during Panxian ham processing.
Figure 3. Effects of spontaneous fermentation on microbial diversity and community structure of Panxian ham. (A) OTU number of bacteria and fungi at different processing stages, values are presented as mean ± standard error (n = 5), a–e different letters represent significant differences (p < 0.05). (B) rarefaction curves of bacteria and fungi for each sample. (C,D) NMDS and PCoA score plots of bacteria and fungi. A–F indicate the traditional spontaneous fermentation process stage of Panxian ham: (A) raw ham; (B) post-salting; (C) post-resting; (D) initial stage of ripening; (E) middle stage of ripening; (F) final stage of ripening.
Table 1. Richness and diversity of bacteria and fungi at different processing stages.
Bacterial and Fungal Community Succession During Panxian Ham Processing
The sequencing data of bacteria and fungi were classified at both the phylum and genus levels to better understand community succession. Ten bacterial phyla were identified; Firmicutes and Proteobacteria were predominant in all samples (Figure 4A). Proteobacteria was the most abundant phylum in raw and post-salting hams, while Firmicutes dominated the resting and ripening stages. With respect to fungi, four phyla were identified, and Ascomycota represented more than 92% of the sequences of each sample (Figure 4B). Furthermore, Cyanobacteria, Actinobacteria, Bacteroidetes, and Basidiomycota, primarily as a result of their presence in raw ham materials, were commonly identified in all early stage ham samples.
Figure 4. Relative abundance of bacteria (A) and fungi (B) at the phylum level and bacteria (C) and fungi (D) at the genus level in samples of Panxian ham. A–F indicate the traditional spontaneous fermentation process stage of Panxian ham: (A) raw ham; (B) post-salting; (C) post-resting; (D) initial stage of ripening; (E) middle stage of ripening; (F) final stage of ripening.
At the genus level, 154 bacterial genera were identified during Panxian ham processing and 11 core genera (whose relative abundance was larger than 1%) are shown in Figure 4C. The raw hams were characterized by Psychrobacter (16.28%), followed by uncultured_bacterium (12.70%), Acinetobacter (8.36%), Ochrobactrum (5.98%), Staphylococcus (5.44%), and Macrococcus (3.79%). After salting, the relative abundance of Psychrobacter was further increased to 31.57%, while the relative abundances of uncultured_bacterium and Ochrobactrum were remarkably reduced to 0.12% and 0.08%, respectively. Interestingly, the proportion of Cobetia, which was reportedly derived from salt (Jung et al., 2015), increased dramatically from an initial value of 1.38 to 27.76%. Staphylococcus then became the dominant genus during the resting and ripening stages, with a relative abundance ranging from 48.74 to 73.87%. Similarly, some researchers have reported that Staphylococcus was dominant in other fermented meat products (Wang et al., 2018; Gassem, 2019). Coagulase-negative staphylococci (CNS) has been regarded as the most important factor in the taste and flavor development of dry-cured ham, after endogenous enzymes (Landeta et al., 2007), owing to its strong enzyme activity, which includes nitrate reductase, catalase, and lipolytic and proteolytic enzymes (Landeta et al., 2013). In addition, CNS also contributed to the preservation of meat products through nitric oxide synthase (Sapp et al., 2014). For other core genera, the relative abundance of Chromohalobacter and Kushneria increased rapidly at the final stage of ripening. This variation might be related to the increase in salinity, because most members of these two genera were salt-tolerant and were often isolated from high salt foods such as cured meat (Zou and Wang, 2010) and salt-fermented fish (Yun et al., 2018).
Some 51 fungal genera were detected and only five taxa, Debaryomyces, Aspergillus, Yamadazyma, Candida, and Penicillium, were identified as core fungal genera (Figure 4D). Two dominant groups were present in all samples, Debaryomyces and Aspergillus, and exhibited obvious antagonism during the whole processing of Panxian ham, which was mainly due to wide nutritional competition between yeasts and molds (Adel Druvefors and Schnurer, 2005). Yeast could promote the sensory characteristics of dry-cured meat products due to exogenous enzyme activities (Dura et al., 2004). As one of the most commonly isolated yeasts from dry-cured ham, Debaryomyces often dominated the yeast population after salting (Núñez et al., 1996). Debaryomyces inoculation has been reported to promote the development of more volatile compounds related to the flavor of dry-cured meat products (Andrade et al., 2009). Thus, Debaryomyces has been widely used as a starter culture to enhance the flavor and aroma of dry-cured meat products (Laranjo et al., 2019). Aspergillus was one of the most abundant molds in dry-cured ham. High-quality molds can not only improve the typical features of fermented meat products, but also inhibit the growth of harmful microorganism (Bruna et al., 2003; Castellari et al., 2010). However, Wang et al. (2006) suggested that the production of some undesirable mycotoxins by molds such as Aspergillus and Penicillium may not contribute to Xuanwei ham quality. Rojas et al. (1991) and Comi and Iacumin (2013) analyzed the toxigenic potential of molds on the surface of Spanish and San Daniele dry-cured ham and concluded that their presence did not pose a health hazard. Meanwhile, the growth of autochthonous yeasts in dry-cured ham has been shown to inhibit the accumulation of ochratoxin A, a significant mycotoxin in processed meats (Simoncini et al., 2014).
Multivariate Statistical Analysis of Metabolites and Microorganisms
It is generally accepted that variance in the microbial community structure will affect the volatile and non-volatile metabolites present in fermented food that determine its sensory properties, which inevitably affects the final product’s industrial value. This is especially true for foods prepared by spontaneous fermentation such as dry-cured ham. Therefore, multivariate RDA was applied to reveal the correlations between Panxian ham samples, core microbial genera, and SDMs.
Statistical RDA triplot analysis (Figure 5) demonstrated that the changes in SDMs during Panxian ham processing were the result of interactions between bacterial and fungal communities, but seemed to be dominated by the bacterial population, particularly, Staphylococcus and Chromohalobacter. These two genera were significantly correlated with the production of all amino acids and fatty acids except stearic acid (Figure 5A). Moreover, Debaryomyces was also associated with the metabolism of six amino acids and two fatty acids (Figure 5B). This correlation could be attributed to microbial enzymes originated from these microorganisms. Of these, the proteolytic and lipolytic activities of Staphylococcus and Debaryomyces were well known (Dura et al., 2004; Landeta et al., 2013). The ability of Chromohalobacter to produce industrially important enzymes, including proteases and lipases, had also been previously reported (Kumar et al., 2012). Furthermore, the fermentation of ganjiang, a Korean traditional fermented soy sauce, was dominated by Chromohalobacter, which was conducive to the generation of amino and organic acids (Jung et al., 2015). Aspergillus and Cobetia had a significant and positive correlation coefficient with glucose and fructose. These findings concur with those reported by other authors who noticed that Aspergillus was primarily responsible for sugar metabolism in fermented foods (Lee et al., 2017; Tang et al., 2019). However, Jung et al. (2015) reported that Cobetia was not associated with the metabolites of ganjiang, although the early stage fermentation samples were represented by Cobetia. In addition, Kushneria, Yamadazyma, and Penicillium were closely related with organic acids, namely, succinic and malic acids. It has been reported that some Penicillium strains isolated on the surface of dry-cured sausages have malate dehydrogenase activity (Vila et al., 2019), which is crucial to the reversible conversion between malate and oxaloacetate. However, microbial community succession did not account for changes in inositol, creatine, creatinine, and inosine concentrations. Seo et al. (2016) suggested that myo-inositol was not technically a fermentation product and was mainly derived from Aspergillus fungus in koji. In addition, the conversion of creatine to creatinine during dry-cured ham processing was dependent on the processing time and initial creatine concentration (Mora et al., 2010).
Figure 5. Redundancy analysis (RDA) of Panxian ham samples, core microbial genera, and significantly different metabolites (SDMs) during Panxian ham processing. (A) Correlation between samples, core bacteria genera, and SDMs. (B) Correlation between samples, core fungi genera, and SDMs. A–F indicate the traditional spontaneous fermentation process stage of Panxian ham: (A) raw ham; (B) post-salting; (C) post-resting; (D) initial stage of ripening; (E) middle stage of ripening; (F) final stage of ripening. The core bacterial and fungal genera were consistent with Figures 4C,D; while SDMs are represented by M1–M31. M1, valine; M2, Isoleucine; M3, proline; M4, glycine; M5, serine; M6, L-allothreonine; M7, aspartic acid; M8, methionine; M9, pyroglutamate; M10, glutamic acid; M11, phenylalanine; M12, ornithine; M13, lysine; M14, tyrosine; M15, alanine; M16, myristic acid; M17, palmitoleic acid; M18, heptadecanoic acid; M19, arachidonic acid; M20, arachidic acid; M21, stearic acid; M22, succinic acid; M23, malic acid; M24, creatine; M25, glycerol; M26, myo-inositol; M27, fructose; M28, glucose; M29, hypoxanthine; M30, inosine; M31, creatinine.
It is worth noting that Acinetobacter, Psychrobacter, Brochothrix, Pseudomonas, Ochrobactrum, Macrococcus, Candida, and uncultured_bacterium showed significant and negative correlation coefficients with most SDMs (Figure 5); these microorganisms represented the microbial community in the early stages of Panxian ham processing (Figure 4). Usually, the incidence and/or prevalence of gram-negative bacteria in fermented meat products is undesirable. Acinetobacter, Psychrobacter, Brochothrix, and Pseudomonas were recognized as meat spoilage factors due to their production of off-flavors, and were a hazardous source of transmissible antibiotic resistance (Wang et al., 2018). Additionally, Pseudomonas has been proven to be related to the production of biogenic amines in fermented foods (Huang et al., 2018). Ochrobactrum is considered an opportunistic bacterium with antimicrobial resistance that has been associated with human wound infections (Furlan and Stehling, 2017). This genus was frequently isolated from soil samples and it was also recently detected in traditional Icelandic fermented fish called hakarl (Osimani et al., 2019). Macrococcus, a gram-positive bacterium closely related to Staphylococcus, was often found in animal skin samples and their meat products (Sae et al., 2010). Previous researchers concluded that Macrococcus has nothing to do with human or animal disease (Lory, 2014), but a strain of Macrococcus causing high mortality has been separated from commercial chickens in China (Li et al., 2018). Meanwhile, Macrococcus was also one of the main microorganisms responsible for chicken spoilage (Arthur et al., 2004). Candida was the most abundant yeast in the fresh ham (Núñez et al., 1996), but it was also recognized as a pathogenic yeast isolated from dry-cured hams (Levenson et al., 1991). Some pathogenic strains of Candida can cause both superficial and systemic infections (Carolus et al., 2019). Fortunately, Staphylococcus, Debaryomyces, Aspergillus, Kushneria, and Chromohalobacter showed strong antagonistic behavior toward these undesirable microorganisms (Supplementary Figure S3). As a result, the relative abundance of these potentially risky microorganisms was reduced and eventually disappeared through traditional spontaneous fermentation (Figure 4), indicating the safety of the final Panxian ham product. The results of gene function prediction from the PICRUSt also supported this viewpoint (Supplementary Figure S4). The abundance of predicted functional genes involved in human diseases decreased significantly during Panxian ham processing. Nevertheless, it is necessary to find measures to reduce the incidence and prevalence of these undesirable microorganisms in the early stages of processing.
Conclusion
Investigating the dynamic changes in microbial community and metabolite profiles provided information regarding the metabolically active microbiota that promote the production of primary metabolites during spontaneous fermentation of Panxian ham. The results showed that eight microbial genera, namely, Staphylococcus, Chromohalobacter, Cobetia, Kushneria, Debaryomyces, Aspergillus, Yamadazyma, and Penicillium, were significantly correlated with changes in SDMs. The prevalence of these microbial genera was crucial to the refinement of the sensory qualities and safety of Panxian ham by inhibiting the growth of some pathogen-related microorganisms. These findings play a vital role in understanding Panxian ham quality parameter development; allowing the production of safe, high quality products. Nevertheless, corresponding this type of microbial analysis with sensory parameters (meat color, pH, water loss, texture), volatiles flavor, and taste evaluation by a sensory panel would be highly relevant. Further, metagenomics and transcriptomics should be applied to elucidate the functions of these microorganisms.
Data Availability Statement
All sequencing data have been deposited at the Sequence Read Archive of the National Center for Biotechnology Information (NCBI), with SRA accession number: the accession number of bacteria and fungi are PRJNA579779 and PRJNA579784, respectively.
Author Contributions
YM and WS contributed to the experimental design. YM performed the statistical analysis and wrote the manuscript. YM, WS, YCM, and LJ contributed to manuscript revision, and read and approved the submitted version.
Funding
This research was financially supported by Guizhou University Doctoral Fund (2017-43).
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.
Acknowledgments
We thank Biotree Biotech Co., Ltd. (Shanghai, China) and Biomarker Technologies Co., Ltd. (Beijing, China) for technical assistance in carrying out the analyses. We would like to thank Editage (www.editage.cn) for English language editing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03012/full#supplementary-material
Footnotes
References
Adel Druvefors, U., and Schnurer, J. (2005). Mold-inhibitory activity of different yeast species during airtight storage of wheat grain. FEMS Yeast Res. 5, 373–378. doi: 10.1016/j.femsyr.2004.10.006
Andrade, M. J., Cordoba, J. J., Sanchez, B., Casado, E. M., and Rodriguez, M. (2009). Evaluation and selection of yeasts isolated from dry-cured Iberian ham by their volatile compound production. Food Chem. 113, 457–463. doi: 10.1016/j.foodchem.2008.07.080
Arthur, H., Cason, J. A., and Ingram, K. D. (2004). Tracking spoilage bacteria in commercial poultry processing and refrigerated storage of poultry carcasses. Int. J. Food Microbiol. 91, 155–165. doi: 10.1016/s0168-1605(03)00377-5
Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59. doi: 10.1038/Nmeth.2276
Bradley, W., and Robert, P. (2016). PCA as a practical indicator of OPLS-DA model reliability. Curr. Metab. 4, 97–103. doi: 10.2174/2213235X04666160613122429
Bruna, J. M., Hierro, E. M., Hoz, L. D. L., Mottram, D. S., Fernández, M., and Ordóñez, J. A. (2003). Changes in selected biochemical and sensory parameters as affected by the superficial inoculation of Penicillium camemberti on dry fermented sausages. Int. J. Food Microbiol. 85, 111–125. doi: 10.1016/S0168-1605(02)00505-6
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
Careri, M., Mangia, A., Barbieri, G., Bouoni, L., Virgili, R., and Parolari, G. (1993). Sensory property relationships to chemical data of Italian-type dry-cured ham. J. Food Sci. 58, 968–972. doi: 10.1111/j.1365-2621.1993.tb06090.x
Carolus, H., Dyck, K. V., and Dijck, P. V. (2019). Candida albicans and Staphylococcus species: a threatening twosome. Front. Microbiol. 10:2162. doi: 10.3389/fmicb.2019.02162
Castellari, C., Quadrelli, A. M., and Laich, F. (2010). Surface mycobiota on Argentinean dry fermented sausages. Int. J. Food Microbiol. 142, 149–155. doi: 10.1016/j.ijfoodmicro.2010.06.016
Comi, G., and Iacumin, L. (2013). Ecology of moulds during the pre-ripening and ripening of san daniele dry cured ham. Food Res. Int. 54, 1113–1119. doi: 10.1016/j.foodres.2013.01.031
Dashdorj, D., Amna, T., and Hwang, I. (2015). Influence of specific taste-active components on meat flavor as affected by intrinsic and extrinsic factors: an overview. Eur. Food Res. Technol. 241, 157–171. doi: 10.1007/s00217-015-2449-2443
Dunn, W. B., Broadhurst, D., Begley, P., Zelena, E., Francis-McIntyre, S., Anderson, N., et al. (2011). Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc. 6, 1060–1083. doi: 10.1038/nprot.2011.335
Dura, M. A., Flores, M., and Toldra, F. (2004). Effect of Debaryomyces spp. on the proteolysis of dry-fermented sausages. Meat Sci. 68, 319–328. doi: 10.1016/j.meatsci.2004.03.015
Escudero, E., Mora, L., Aristoy, M. C., and Toldra, F. (2011). Possible biological markers of the time of processing of dry-cured ham. Meat Sci. 89, 536–539. doi: 10.1016/j.meatsci.2011.05.002
Fadrosh, D. W., Ma, B., Gajer, P., Sengamalay, N., Ott, S., Brotman, R. M., et al. (2014). An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome 2, 2–6. doi: 10.1186/2049-2618-2-6
Fakruddin, M., Mannan, K. S., and Andrews, S. (2013). Viable but nonculturable bacteria: food safety and public health perspective. ISRN Microbiol. 2013, 1–6. doi: 10.1155/2013/703813
Fernández, M., Ordóñez, J. A., Cambero, I., Santos, C., Pin, C., and Hoz, L. D. L. (2007). Fatty acid compositions of selected varieties of Spanish dry ham related to their nutritional implications. Food Chem. 101, 107–112. doi: 10.1016/j.foodchem.2006.01.006
Furlan, J. P. R., and Stehling, E. G. (2017). High-level of resistance to β-lactam and presence of β-lactamases encoding genes in Ochrobactrum sp. and Achromobacter sp. isolated from soil. J. Glob. Antimicrob. Res. 11, 133–137. doi: 10.1016/j.jgar.2017.10.014
Gassem, M. A. (2019). Microbiological and chemical quality of a traditional salted-fermented fish (Hout-Kasef) product of Jazan Region, Saudi Arabia. Saudi J. Biol. Sci. 26, 137–140. doi: 10.1016/j.sjbs.2017.04.003
Ge, Q., Gu, Y., Zhang, W., Yin, Y., Yu, H., Wu, M., et al. (2017). Comparison of microbial communities from different Jinhua ham factories. AMB Express. 7:37. doi: 10.1186/s13568-017-0334-0
Grossi, A. B., do Nascimento, E. S. P., Cardoso, D. R., and Skibsted, L. H. (2014). Proteolysis involvement in zinc-protoporphyrin IX formation during Parma ham maturation. Food Res. Int. 56, 252–259. doi: 10.1016/j.foodres.2014.01.007
Hao, R., Du, X., Yang, C., Deng, Y., Zheng, Z., and Wang, Q. (2019). Integrated application of transcriptomics and metabolomics provides insights into unsynchronized growth in pearl oyster Pinctada fucata martensii. Sci. Total Environ. 666, 46–56. doi: 10.1016/j.scitotenv.2019
Huang, X. N., Yu, S. Z., Han, B. Z., and Chen, J. Y. (2018). Bacterial community succession and metabolite changes during sufu fermentation. LWT Food Sci. Technol. 97, 537–545. doi: 10.1016/j.lwt.2018.07.041
Ichimura, S., Nakamura, Y., Yoshida, Y., and Hattori, A. (2017). Hypoxanthine enhances the cured meat taste. Anim. Sci. J. 88, 379–385. doi: 10.1111/asj.12625
Jimenez-Colmenero, F., Ventanas, J., and Toldra, F. (2010). Nutritional composition of dry-cured ham and its role in a healthy diet. Meat Sci. 84, 585–593. doi: 10.1016/j.meatsci.2009.10.029
Jung, J. Y., Chun, B. H., and Jeon, C. O. (2015). Chromohalobacter is a causing agent for the production of organic acids and putrescine during fermentation of ganjang, a korean traditional soy sauce. J. Food Sci. 80, 2853–2859. doi: 10.1111/1750-3841.13114
Kim, T. K., Hwang, K. E., Lee, M. A., Paik, H. D., Kim, Y. B., and Choi, Y. S. (2019). Quality characteristics of pork loin cured with green nitrite source and some organic acids. Meat Sci. 152, 141–145. doi: 10.1016/j.meatsci.2019.02.015
Kind, T., Wohlgemuth, G., Lee, D. Y., Lu, Y., Palazoglu, M., Shahbaz, S., et al. (2009). Mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal. Chem. 81, 10038–10048. doi: 10.1021/ac9019522
Koutina, G., Jongberg, S., and Skibsted, L. H. (2012). Protein and lipid oxidation in Parma ham during production. J. Agri. Food Chem. 60, 9737–9745. doi: 10.1021/jf3026887
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/Aem.01043-13
Kumar, S., Karan, R., Kapoor, S., Singh, S. P., and Khare, S. K. (2012). Screening and isolation of halophilic bacteria producing industrially important enzymes. Braz. J. Microbiol. 43, 1595–1603. doi: 10.1590/S1517-83822012000400044
Landeta, G., Curiel, J. A., Carrascosa, A. V., Muñoz, R., and de las Rivas, B. (2013). Characterization of coagulase-negative staphylococci isolated from spanish dry cured meat products. Meat Sci. 93, 387–396. doi: 10.1016/j.meatsci.2012.09.019
Landeta, G., de las Rivas, B., Carrascosa, A. V., and Muñoz, R. (2007). Screening of biogenic amine production by coagulase-negative staphylococci isolated during industrial spanish dry-cured ham process. Meat Sci. 77, 556–561. doi: 10.1016/j.meatsci.2007.05.004
Langille, M. G. I., Jesse, Z., Caporaso, J. G., Daniel, M. D., Dan, K., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31:814. doi: 10.1038/nbt.2676
Laranjo, M., Potes, M. E., and Elias, M. (2019). Role of starter cultures on the safety of fermented meat products. Front. Microbiol. 10:853. doi: 10.3389/fmicb.2019.00853
Lee, S., Lee, S., Singh, D., Oh, J. Y., Jeon, E. J., Ryu, H. S., et al. (2017). Comparative evaluation of microbial diversity and metabolite profiles in doenjang, a fermented soybean paste, during the two different industrial manufacturing processes. Food Chem. 221, 1578–1586. doi: 10.1016/j.foodchem.2016.10.135
Levenson, D., Pfaller, M. A., Smith, M. A., Hollis, R., Gerarden, T., Tucci, C. B., et al. (1991). Candida zeylanoides: another opportunistic yeast. J. Clin. Microbiol. 29, 1689–1692.
Li, G., Du, X., Zhou, D., Li, C., Huang, L., Zheng, Q., et al. (2018). Emergence of pathogenic and multiple-antibiotic-resistant Macrococcus caseolyticus in commercial broiler chickens. Transbound. Emerg. Dis. 65, 1605–1614. doi: 10.1111/tbed.12912
Liu, C., Pan, D., Ye, Y., and Cao, J. (2013). 1H NMR and multivariate data analysis of the relationship between the age and quality of duck meat. Food Chem. 141, 1281–1286. doi: 10.1016/j.foodchem.2013.03.102
Liu, S. Y., Wang, G. Y., Xiao, Z. C., Pu, Y. H., Ge, C. R., and Liao, G. Z. (2019). 1H-NMR-based water-soluble low molecular weight compound characterization and free fatty acid composition of five kinds of Yunnan dry-cured hams. LWT Food Sci. Technol. 108, 174–182. doi: 10.1016/j.lwt.2019.03.043
Lory, S. (2014). “The family staphylococcaceae,” in The Prokaryotes: Firmicutes and Tenericutes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson, (Berlin: Springer), 363–366. doi: 10.1007/978-3-642-30120-9_350
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Martinez-Sanchez, S. M., Minguela, A., Prieto-Merino, D., Zafrilla-Rentero, M. P., Abellan-Aleman, J., and Montoro-Garcia, S. (2017). The Effect of regular intake of dry-cured ham rich in bioactive peptides on inflammation, platelet and monocyte activation markers in humans. Nutrients 9:321. doi: 10.3390/nu9040321
Marusic, N., Aristoy, M. C., and Toldra, F. (2013). Nutritional pork meat compounds as affected by ham dry-curing. Meat Sci. 93, 53–60. doi: 10.1016/j.meatsci.2012.07.014
Marušić, N., Vidaèek, S., Janèi, T., Petrak, T., and Mediæ, H. (2014). Determination of volatile compounds and quality parameters of traditional Istrian dry-cured ham. Meat Sci. 96, 1409–1416. doi: 10.1016/j.meatsci.2013.12.003
Montanari, C., Barbieri, F., Magnani, M., Grazia, L., Gardini, F., and Tabanelli, G. (2018). Phenotypic diversity of Lactobacillus sakei strains. Front. Microbiol. 9:2003. doi: 10.3389/fmicb.2018.02003
Montoro-Garcia, S., Zafrilla-Rentero, M. P., Celdran-de Haro, F., Pinero-de Armas, J. J., Toldra, F., Tejada-Portero, L., et al. (2017). Effects of dry-cured ham rich in bioactive peptides on cardiovascular health: a randomized controlled trial. J. Funct. Foods 38, 160–167. doi: 10.1016/j.jff.2017.09.012
Mora, L., Hernandez-Cazares, A. S., Sentandreu, M. A., and Toldra, F. (2010). Creatine and creatinine evolution during the processing of dry-cured ham. Meat Sci. 84, 384–389. doi: 10.1016/j.meatsci.2009.09.006
Motilva, M., Toldrfi, F., Nieto, P., and Flores, J. (1993). Muscle lipolysis phenomena in the processing of dry-cured ham. Food Chem. 48, 121–125. doi: 10.1016/0308-8146(93)90045-h
Mu, Y., Su, W., and Mu, Y. C. (2019). Analysis of microbial diversity and key volatile flavor compounds of Panxian dry-cured ham. Food Res. Dev. 40, 77–85. doi: 10.12161/j.issn.1005-6521.2019.15.013
Núñez, F., Rodríguez, M. M., Córdoba, J. J., Bermúdez, M. E., and Asensio, M. A. (1996). Yeast population during ripening of dry-cured Iberian ham. Int. J. Food Microbiol. 29, 271–280. doi: 10.1016/0168-1605(95)00037-2
Osimani, A., Ferrocino, I., Agnolucci, M., Cocolin, L., Giovannetti, M., Cristani, C., et al. (2019). Unveiling hakarl: a study of the microbiota of the traditional Icelandic fermented fish. Food Microbiol. 82, 560–572. doi: 10.1016/j.fm.2019.03.027
Park, S. E., Seo, S. H., Kim, E. J., Byun, S., Na, C. S., and Son, H. S. (2019). Changes of microbial community and metabolite in kimchi inoculated with different microbial community starters. Food Chem. 274, 558–565. doi: 10.1016/j.foodchem.2018.09.032
Petrova, I., Aasen, I. M., Rustad, T., and Eikevik, T. M. (2015). Manufacture of dry-cured ham: a review. Part 1. Biochemical changes during the technological process. Eur. Food Res. Technol. 241, 587–599. doi: 10.1007/s00217-015-2490-2
Rojas, F. J., Jodral, M., Gosalvez, F., and Pozo, R. (1991). Mycoflora and toxigenic aspergillus flavus in Spanish dry-cured ham. Int. J. Food Microbiol. 13, 249–255. doi: 10.1016/0168-1605(91)90082-z
Sae, T., Kyoko, K. A., Tadashi, B., and Keiichi, H. J. A. A. C. (2010). Staphylococcal cassette chromosome mec-like element in Macrococcus caseolyticus. Antimicrob. Agents Chemother. 54, 1469–1475. doi: 10.1128/AAC.00575-09
Sapp, A. M., Mogen, A. B., Almand, E. A., Rivera, F. E., Shaw, L. N., Richardson, A. R., et al. (2014). Contribution of the nos-pdt operon to virulence phenotypes in methicillin-sensitive Staphylococcus aureus. PLoS One 9:e108868. doi: 10.1371/journal.pone.0108868
Seo, S. H., Park, S. E., Yoo, S. A., Lee, K. I., Na, C. S., and Son, H. S. (2016). Metabolite profiling of Makgeolli for the understanding of yeast fermentation characteristics during fermentation and aging. Process. Biochem. 51, 1363–1373. doi: 10.1016/j.procbio.2016.08.005
Shi, Y., Li, X., and Huang, A. (2019). A metabolomics-based approach investigates volatile flavor formation and characteristic compounds of the Dahe black pig dry-cured ham. Meat Sci. 158:107904. doi: 10.1016/j.meatsci.2019.107904
Simoncini, N., Rotelli, D., Virgili, R., and Quintavalla, S. (2007). Dynamics and characterization of yeasts during ripening of typical italian dry-cured ham. Food Microbiol. 24, 577–584. doi: 10.1016/j.fm.2007.01.003
Simoncini, N., Virgili, R., Spadola, G., and Battilani, P. (2014). Autochthonous yeasts as potential biocontrol agents in dry-cured meat products. Food Control 46, 160–167. doi: 10.1016/j.foodcont.2014.04.030
Sugimoto, M., Obiya, S., Kaneko, M., Enomoto, A., Honma, M., Wakayama, M., et al. (2017). Metabolomic profiling as a possible reverse engineering tool for estimating processing conditions of dry-cured hams. J. Agric. Food Chem. 65, 402–410. doi: 10.1021/acs.jafc
Tamang, J. P., Koichi, W., and Holzapfel, W. H. (2016). Review: diversity of microorganisms in global fermented foods and beverages. Front. Microbiol. 7:377. doi: 10.3389/fmicb.2016.00377
Tang, Q., He, G., Huang, J., Wu, C., Jin, Y., and Zhou, R. (2019). Characterizing relationship of microbial diversity and metabolite in Sichuan Xiaoqu. Front. Microbiol. 10:696. doi: 10.3389/fmicb.2019.00696
Toldrá, F., and Flores, M. J. F. C. (2000). The use of muscle enzymes as predictors of pork meat quality. Food Chem. 69, 387–395. doi: 10.1016/S0308-8146(00)00052-2
Tu, R. J., Wu, H. Y., Lock, Y. S., and Chen, M. J. (2010). Evaluation of microbial dynamics during the ripening of a traditional Taiwanese naturally fermented ham. Food Microbiol. 27, 460–467. doi: 10.1016/j.fm.2009.12.011
Van Reckem, E., Geeraerts, W., Charmpi, C., Van der Veken, D., De Vuyst, L., and Leroy, F. (2019). Exploring the link between the geographical origin of European fermented foods and the diversity of their bacterial communities: the case of fermented meats. Front. Microbiol. 10:2302. doi: 10.3389/fmicb.2019.02302
Vila, G., Segura, J. A., Ludemann, V., and Pose, G. N. (2019). Surface mycobiota of home-made dry cured sausages from the main producing regions of Argentina and morphological and biochemical characterization of Penicillium nalgiovense populations. Int. J. Food Microbiol. 309:10. doi: 10.1016/j.ijfoodmicro.2019.108312
Wang, X., Ping, M., Jiang, D. F., Qian, P., and Yang, H. Y. (2006). The natural microflora of Xuanwei ham and the no-mouldy ham production. J. Food Eng. 77, 103–111. doi: 10.1016/j.jfoodeng.2005.06.047
Wang, X. H., Wang, S. H., and Zhao, H. (2019). Unraveling microbial community diversity and succession of Chinese Sichuan sausages during spontaneous fermentation by high-throughput sequencing. J. Food Sci. Technol. 56, 3254–3263. doi: 10.1007/s13197-019-03781-y
Wang, X. H., Zhang, Y. L., Ren, H. Y., and Zhan, Y. (2018). Comparison of bacterial diversity profiles and microbial safety assessment of salami, Chinese dry-cured sausage and Chinese smoked-cured sausage by high-throughput sequencing. LWT Food Sci. Technol. 90, 108–115. doi: 10.1016/j.lwt.2017.12.011
Wurihan, B. L. S., Hasigaowa, Bao, X. P., Dai, Y. J., and Jia, S. R. (2019). Bacterial community succession and metabolite changes during the fermentation of koumiss, a traditional Mongolian fermented beverage. Int. Dairy J. 98, 1–8. doi: 10.1016/j.idairyj.2019
Yang, C., Hao, R., Du, X., Wang, Q., Deng, Y., Sun, R., et al. (2018). GC-TOF/MS-based metabolomics studies on the effect of protein sources in formulated diet for pearl oyster Pinctada fucata martensii. Aquaculture 486, 139–147. doi: 10.1016/j.aquaculture.2017.12.020
Yun, J. H., Sung, H., Kim, H. S., Tak, E. J., Kang, W., Lee, J. Y., et al. (2018). Complete genome sequence of the halophile bacterium Kushneria konosiri X49T, isolated from salt-fermented Konosirus punctatus. Stand Genomic Sci. 13:19. doi: 10.1186/s40793-018-0324-0
Zang, J. H., Xu, Y. S., Xia, W. S., Yu, D. W., Gao, P., Jiang, Q. X., et al. (2018). Dynamics and diversity of microbial community succession during fermentation of Suan yu, a Chinese traditional fermented fish, determined by high throughput sequencing. Food Res. Int. 111, 565–573. doi: 10.1016/j.foodres.2018.05.076
Zhang, J., Ye, Y., Sun, Y., Pan, D., Ou, C., Dang, Y., et al. (2018). 1H NMR and multivariate data analysis of the differences of metabolites in five types of dry-cured hams. Food Res. Int. 113, 140–148. doi: 10.1016/j.foodres.2018.07.009
Zhang, J., Yi, Y., Pan, D. D., Zhou, G. H., Wang, Y., Dang, Y. L., et al. (2019). 1H NMR-based metabolomics profiling and taste of boneless dry-cured hams during processing. Food Res. Int. 122, 114–122. doi: 10.1016/j.foodres.2019.04.005
Zheng, J. Y., Tao, N. P., and Gong, J. (2015). Comparison of non-volatile taste-active compounds between the cooked meats of pre- and post-spawning Yangtze Coilia ectenes. Fish Sci. 81, 559–568. doi: 10.1007/s12562-015-0858-7
Zhou, G. H., and Zhao, G. M. (2007). Biochemical changes during processing of traditional Jinhua ham. Meat Sci. 77, 114–120. doi: 10.1016/j.meatsci.2007.03.028
Keywords: Panxian ham, metabolomics, high-throughput sequencing, non-volatile metabolites, metabolically active microorganism
Citation: Mu Y, Su W, Mu Y and Jiang L (2020) Combined Application of High-Throughput Sequencing and Metabolomics Reveals Metabolically Active Microorganisms During Panxian Ham Processing. Front. Microbiol. 10:3012. doi: 10.3389/fmicb.2019.03012
Received: 31 October 2019; Accepted: 16 December 2019;
Published: 10 January 2020.
Edited by:
Shao Quan Liu, National University of Singapore, SingaporeReviewed by:
Hong-Seok Son, Dongshin University, South KoreaKumaraswamy Jeyaram, Institute of Bioresources and Sustainable Development (IBSD), India
Copyright © 2020 Mu, Su, Mu and Jiang. 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: Wei Su, wsu@gzu.edu.cn