Deciphering the Mechanism by Which Carbon Dioxide Extends the Shelf Life of Raw Milk: A Microbiomics- and Metabolomics-Based Approach

1. Introduction

Raw milk is commonly produced intensively due to its collection from dispersed farms and subsequent transportation to processing plants for uniform treatment. Throughout the stages of collection, transportation, and waiting for processing, raw milk undergoes a significant duration (several hours), imposing heightened demands on its quality [1,2,3]. Microorganisms play a pivotal role in influencing the quality of raw milk. Given its nutrient-rich composition, raw milk provides an optimal environment and nutrient matrix that facilitate rapid microbial reproduction [4]. Consequently, this can result in undesirable outcomes, such as off-flavors, sourness, discoloration, and a loss of the characteristic freshness associated with milk [1]. Additionally, microorganisms break down and consume essential nutrients present in milk, including fats, proteins and vitamins. This process may lead to a reduction in the nutrient content of raw milk, thereby impacting the nutritional value of the final dairy product [5]. To ensure the provision of high-quality and safe dairy products while prioritizing consumer health and satisfaction, it is imperative to implement effective control measures aimed at mitigating or eliminating microbial hazards inherent in raw milk [2,4].

The sources of microorganisms in raw milk are much more varied; microorganisms can be present in the skin of the cow's udder and teats, and are passed on to the milk when skin microorganisms come into contact with the milk during the milking process [4,5]. In addition, the environment, insects and rodents, milking equipment, and collection devices are all potential sources of microorganisms in raw milk [6]. Notably, these sources of microbial contamination highlight the need to inhibit microbial colonisation through effective technological means [7]. Researchers have explored the use of high pressure to inactivate microorganisms in raw milk while preserving its nutritional quality. However, due to the complexity of high-pressure equipment and limited treatment efficiency, it is not suitable for large-scale implementation [8]. Non-thermal plasma treatment offers an alternative approach to generating reactive oxygen and nitrogen species that can penetrate microorganisms and inhibit or kill them. However, these oxides can also cause oxidation of raw milk, leading to a decline in its quality as a complex emulsion system [9]. Conventional sterilization methods may destabilize raw milk by affecting its composition and structure. To overcome these challenges, alternative technologies, such as gas conditioning or non-thermal treatments, are being investigated to minimize adverse effects on raw milk while preserving its nutritional value and sensory properties [10].

CO₂ has been used for fresh milk storage and preservation because it does not affect the nutritional content or flavor compared to conventional heat treatments. Also, CO₂ is a natural compound that does not leave harmful residues in the product or the environment [11,12]. CO₂ has been used for food preservation. When CO₂ is dissolved in milk, it forms carbonic acid, which lowers the pH of the milk. Many bacteria, including those that cause spoilage, cannot survive or grow effectively in this more acidic environment [13,14]. In addition, by replacing oxygen with CO₂ in the storage container, packaging or food systems (a technique known as Modified Atmosphere Packaging), the growth of aerobic bacteria and fungi is inhibited. These microorganisms require oxygen to grow, thus removing or reducing oxygen levels can help extend the shelf life of raw milk [15].

In view of the complexity of the raw milk system and microbial community succession, the metabolic pattern of CO₂ treatment on raw milk and microorganisms was analyzed by using macro-genome and metabolome, and then the molecular mechanism by which CO₂ prolongs the shelf life of raw milk was clarified. The results of the study provide a theoretical basis for the application of CO₂ in raw milk.

2. Materials and methods

2.1. Sample collection

Fresh raw milk (the somatic cells < 2×10⁵/mL, total bacteria count < 10³ CFU/mL）was collected aseptically from milk tanks at Helan Mountain Dairy Farm, Ningxia, Yinchuan, China in June 2022, refrigerated at 4 ℃ and transported to the laboratory within one hour. Then, a known mass of solid CO₂ (dry ice, food grade) was added to 500 mL PET bottles (Servicebio, Wuhan, China) containing 300 mL of raw milk to achieve a CO₂ concentration of 2000 ppm in a sterile room immediately [16]. Raw milk without CO₂ treatment was used as a control group. Each group of samples was divided into 100 ml aliquots in sterile sealed PET bottles and stored at 4 ℃. One aliquot was used for analysis each day. The day of sample collection was defined as day 0.

2.2. Bacterial growth studies

The total bacterial count (TBC) was measured by the spread plate method. Briefly, 1 mL sample was taken and serially diluted 10 times with 0.85% sterile saline. Then the 200 uL dilutions were spread onto Plate Count Agar (PCA) and cultured at 37 ℃ for 48 h. Analysis of a sample was terminated when it was spoiled, i.e. the TBC reached a threshold of 6 lg(cfu/mL) [17].

2.3. pH measurement

Raw milk pH was determined using a pH acidimeter (PHS-3C, Shanghai Yidian Scientific Instrument Co. Ltd., Shanghai, China) with a three-point calibration standard prior to testing.

2.4. Metagenomics analysis of raw milk

2.4.1. Raw milk DNA extraction and processing

After centrifugation of raw cow's milk, all microbial DNA in raw cow's milk was extracted by DNeasy PowerFood Microbial DNA Isolation Kit (QIAGEN, Dusseldorf, Germany). DNA integrity was determined by 1.2% agarose gel and gel Doc EZ imager (Bio-Rad, California, USA), and DNA concentration was measured by Qubit 4.0 Fluorescence Quantification Instrument (Invitrogen, Carlsbad, California, US). Qubit™ dsDNA HS Assay Kit (Invitrogen, Carlsbad, California, US) was used to accurately quantify the concentration of the genome and determine the total amount of DNA required for library construction. The DNA was fragmented using Covaris 220 ultrasonic crusher (Woburn, MA., US), concentrated and recovered, and then purified and its concentration determined using Qubit 4.0.

2.4.2. Library construction and sequencing

Library construction was carried out according to the instructions of the library construction kit (ONT, Oxford, UK) for Hieff NGS® high-throughput sequencing platform. Briefly, the purified DNA was subjected to end modification using NEBNext FFPE DNA Repair kit (New England Biolabs, Ipswich, MA, USA), junction ligation and magnetic bead sorting to purify the ligated product, and the purified product was subjected to PCR amplification and enrichment, and the PCR amplification product was purified, and the size of library was detected by 2% agarose gel electrophoresis, and the concentration of library was measured by Qubit 4.0, and the qualified library was used. After the libraries were qualified, the Illumina NovaSeq platform was used to analyze the libraries by end-pair sequencing.

2.4.3. Data analysis

In order to ensure the quality of information, the raw data were filtered, and the clean data were obtained by removing the sequences with junctions and low-quality sequences. After removing the junctions in the reads, global cropping and window quality cropping, and the host contamination, the quality of the sequenced raw data were evaluated by Fastp-0.36, and the validated reads were analyzed statistically after the data filtering.

2.5. Metabolomics analysis of raw milk

2.5.1. Metabolite extraction

Raw cow's milk was vortexed at a low speed, 100 μL of sample was pipetted into a centrifuge tube, then 400 μL of pre-cooled pure methanol solution was added, vortexed and mixed, then sonicated in an ice bath for 20 min, and then allowed to stand at -20℃ for 1 h. Frozen samples were centrifuged at 16,000 × g, 4 ℃ for 20 min, and the supernatant was dried in a high-speed vacuum concentrator centrifuge. The supernatant was dried in a high-speed vacuum concentration centrifuge. The dried sample was raw milk metabolites, and 100 μL of methanol-water solution (1:1, v/v) was added to dissolve the metabolites for spectrometry, and then centrifuged at 20,000 × g, at 4 ℃ for 15 min, and then the supernatant was taken as the sample for analysis.

2.5.2. LC-MS/MS analysis

2.5.2.1. Chromatographic separation

The samples were analyzed in an autosampler at 4°C. The samples were separated on a SHIMADZU-LC30 ultra-high-performance liquid chromatography (UHPLC) system using an ACQUITY UPLC® HSS T3 column (2.1×100 mm, 1.8 µm) (Waters, Milford, MA, USA). The flow rate was 0.3 mL/min at 40 ℃ with an injection volume of 4 μL. The chromatographic mobile phases were A: 0.1% formic acid aqueous solution, and B: acetonitrile. The chromatographic gradient elution program was as follows: from 0 to 2 min, 100% A; from 2 to 6 min, A varied linearly from 100% to 52%; from 6 to 10 min, A varied linearly from 52% to 0%; from 10 to 12 min, A was maintained at 0%; from 12 to 12.1 min, A was maintained at 0%; from 12 to 12.1 min, A was maintained at 0%; and from 12 to 12.1 min, A was maintained at 0%. 0%; from 12 to 12.1min, A varied linearly from 0% to 100%, and from 12.1 to 15min, A was maintained at 100%.

2.5.2.2. Mass spectrometry acquisition

Each sample was collected by electrospray ionization (ESI) in positive (+) and negative (-) modes. The samples were separated by UPLC and analyzed by mass spectrometry using a QE Plus mass spectrometer (Thermo Scientific). The ionization conditions were as follows: Spray Voltage: 3.8 kv (+) and 3.2 kv (-); Capillary Temperature: 320 (±); Spray Voltage: 3.8kv (+) and 3.2kv (-); Capillary Temperature: 320 (±); Sheath Gas: 30 (±); Aux Gas: 5 (±); Probe Heater Temp: 350 (±); S-Lens RF Level: 50.

2.5.3. Data pre-processing

The raw data were subjected to peak alignment, retention time correction and peak area extraction using MSDIAL software. The structural identification of metabolites was performed by exact mass number matching (mass tolerance < 10 ppm) and secondary spectra matching (mass tolerance < 0.01 Da), and the public databases, such as HMDB, MassBank, GNPS, etc., as well as the self-constructed metabolite standard library of Dubai Spectrum (BP-DB), were searched. Ion peaks with > 50% missing values were removed from the extracted data and were not included in the subsequent statistical analysis.

2.6. Statistical analysis

To identify the perturbed biological pathways, the differential metabolite data were subjected to KEGG pathway analysis using KEGG database (http://www.kegg.jp). KEGG enrichment analyses were carried out with the Fisher’s exact test, and an FDR correction was applied to account for multiple testing. Enriched KEGG pathways were considered statistically significant at the P < 0.05 level. All physicochemical data were analysed by one way ANOVA and t-test using SPSS 25 program (SPSS Inc. NY, USA) and significant difference was considered as P < 0.05. Plotting was done using Origin 2023 (Origin Lab., Northampton, USA).

3. Results

3.1. Physicochemical characteristics of raw milk

Total bacterial counts were utilized to assess the storage period of raw milk (Figure 1). The control group was identified as spoilage on day 6, with the TBC greater than 6 lg(cfu/mL). In contrast, the TBC of the CO₂-treated group was only 4.03 lg(cfu/mL) on day 6, and did not exceed 6 lg(cfu/mL) until day 16, indicating that the CO₂ treatment prolonged the storage time of raw milk by significantly inhibiting the growth of the bacteria. It is noteworthy that the TBC of the control raw milk increased steadily during the first 3 days, and showed a logarithmic increase on day 4, suggesting that day 4 of refrigeration may be a critical time point for changes in the quality of raw milk. Therefore, this study selected samples from days 0, 4 and 6 for metagenomic and metabolomic analysis to explore the mechanism of CO₂ prolonging the cold storage period of raw milk. In addition, CO₂ at 2000 ppm of significantly reduced the pH of raw milk (P < 0.05).

3.2. Microbial diversity of raw milk

Analysis of the composition of the biome in the raw milk using the Illumina sequencing platform showed that there were 221 ASVs in the raw milk. The control group contained a total of 189 ASVs, 136 of which were unique to that group. Of these, group A0 had 62 ASVs, 27 of which were unique to that group; group A4 had 67 ASVs, 3 of which were unique to that group; group A6 had 176 ASVs, 114 of which were unique to that group; and the CO₂-treated group contained a total of 85 ASVs, 5 of which were unique to that group. Of these, the day 4 group had 41 ASVs, 1 of which was unique to that group; and the day 6 group had 48 ASVs, 4 of which were unique to that group (Figure 2A). These findings suggested that CO₂ treatment reduced the microbial diversity compared to the control group. The community richness of raw milk microbiota (Simpson diversity index and Chao1 richness estimator) also showed that CO₂ treatment reduced raw milk microbial α-diversity (Figure 2B and D). Principal Component Analysis (PCA) showed that the samples from day 0 clustered in the third quadrant. With the extension of refrigeration time, the control group samples continued to move towards the fourth quadrant along the x-axis direction, while the CO₂ treated samples continued to move towards the second quadrant along the y-axis direction. There was a significant separation between the two groups of samples, confirming a significant difference in microbial β-diversity between the CO₂ treated raw milk and the control group (Figure 2C).

3.3. Microbiological composition of raw milk

The differences in microbial community composition between the two groups were analyzed by nonparametric Kruskal-Wallis and Wilcoxon rank-sum tests, combined with linear discriminant analysis (LDA) and effect size. A total of 100 microbial markers with LDA > 2 and P < 0.05 were identified at different levels. Notably, Figure 2A and 3A showed consistent differences in α- and β-diversity indices, especially in the A6 microbial community, the number of biomarkers was much higher than the D6 group, suggesting that the CO₂ treatment inhibited microbial reproduction. The LDA scores for microbial markers were shown in Figure 3B, and the initial microbial communities significantly enriched in raw milk included s_Stenotrophomonas_maltophilia, g_Stenotrophomonas, f_Xanthomonadaceae, s_Lactobacillus_helveticus, g_Lactobacillus_helveticus and g_Lactobacillus. When raw milk spoiled, the microorganisms significantly enriched in group A6 were mainly k_Bacteria, p_Proteobacteria, c_Gammaproteobacteria, f_Pseudomonadaceae, o_Pseudomonas and g_Pseudomonas, suggesting that these may be the key microorganisms responsible for the deterioration of raw milk quality. The significantly enriched in the D6 group treated with CO₂ were p_Firmicutes, c_Bacilli，o_Lactobacillales, f_Lactobacillaceae, g_Leuconostoc and s_Leuconostoc_mesenteroides. The sequences were annotated at the genus level for classification, and the results were shown in Figure 3C. Stenotrophomonas accounted for 70.2% of the initial flora of raw milk (A0/D0), and the abundance gradually decreased with the extension of the cold storage time, but still accounted for the highest proportion in groups A4, D4 and D6 (45.7%, 62.1% and 57.9%, respectively). CO₂ treatment had a significant effect on the microbial species at the genus level, with Pseudomonas, Actinobacteria and Serratia significantly increasing to 49.2%, 20.9% and 14.6% in the control group at day 6, but only to 0.09%, 3.0% and 0.03% in the D6 group.

3.4. Metabolites

3.4.1. Data quality control

Microorganisms use raw milk as a substrate for reproduction and, in addition to changes in microbial numbers and species, alter the raw milk metabolite system through their own metabolic activities. The study aimed to assess how CO₂ treatment affects raw milk metabolite evolution by influencing microbial succession and, consequently, raw milk quality. LC-MS-based untargeted metabolomics identified a total of 1234 metabolites in six batches of raw milk. PCA was initially used as a way to explore the structure of the data and assess its distribution and metabolite variability in PC space. Six raw milk samples and four QC samples were used as the dataset. First, according to Figure 4 A and B, the distributions of the four QC samples in the PCA model overlapped almost completely, and the correlation coefficients were all greater than 0.998, suggesting that mass spectrometry detection was reliable. As shown in Figure 4A, A0 and D0 sample metabolites had good similarity, indicating that carbon dioxide did not affect the quality of raw milk. A4, A6, D4 and D6 were distinguished from the 0-day raw milk group. However, the A4, D4 and D6 metabolites partially overlaped in the whole, indicating that it is not feasible to use PCA unsupervised modeling alone to distinguish the raw milk metabolite groups.

3.4.2. Overview of differences between samples

Orthogonal partial least squares data analysis (OPLS-DA) was further applied to develop discriminative models that could differentiate the metabolites of the control group and CO₂-treated group with different storage times. The results of the OPLS-DA models developed showed that R2X = 32.8%, R2Y = 99.7%, and Q2 = 95.2 in the control group and R2X = 35%, R2Y = 98.3%, and Q2 = 89.7% in the CO₂-treated group, and the Q2 was greater than 0.5 based on the results of 500 tests of cross-validated alignments (Figure 5) with Q2s of 0.8988 and 0.8735 both being greater than 0.5, indicating that these models are reliable, and different treatments and storage time have significant effects on raw milk metabolites.

3.4.3. Differential metabolites (DMs)

OPLS-DA analysis was used to obtain the variable importance for the projection (VIP) of metabolites to measure the strength of the influence of metabolite expression patterns on the categorical discrimination of each group of samples and the ability to explain, so as to assist in the screening of the marker metabolites (VIP > 1.0 was used as the screening criteria). DMs were indicated by VIP >1 and p < 0.05, as shown in Figure 6 A and C, 402 and 377 DMs were screened in the control and CO₂-treated groups, respectively, where red color in the clustering heatmap indicates up-regulation of metabolites and blue color indicates down-regulation.

The structure and function of the metabolites and the differential metabolites in each comparison group were categorized and counted by the Class and Ontology substance classification method of HMDB. As shown in Figure 6B and D, the main DMs in the raw milk were carboxylic acids and derivatives (A-14.32%, D-14.75%), fatty Acyls (A-12.24%, D-7.77%), benzene and substituted derivatives (A-9.11%, D-12.06%), organooxygen compounds (A-5.47%, D-6.7%), and so on. Notably, the Class classification of metabolites in the control and CO₂-treated groups showed significant differences, mainly in the increase of 10 benzene and substituted derivatives and the decrease of 18 fatty acyl groups in the CO₂-treated group.

3.5. Combined microbiomics and metabolomics analysis

Correlation analysis and elemental screening were performed to investigate the relationship between key microbial community structures and metabolite level parameters. As shown in Figure 7 A and B, some microbes with metabolite changes were identified. Notably, the control microbial community showed a stronger positive correlation with metabolites, indicating that CO₂ treatment was effective in reducing the metabolic activities of the raw milk system. The relationship between microbes and metabolites was shown by correlation hierarchical clustering heatmaps (Figure 7 C and D). In the control group, the top 28 microorganisms (Mycoplasmopsis, Aequoribacter, Aeromonas, Carnobacterium, Citrobacter, Cupriavidus, etc.) and metabolites (delta-tridecalactone, cytidine, alpha-aminoorcein, D-glucosamine, 3-hydroxycapric acid, etc.) were correlated. Only five microorganisms, Stenotrophomonas, Anaerobutyricum, Blautia, Leuconostoc and Propionibacterium were correlated with the metabolites in the CO₂-treated group compared to the control group.

3.6. Metabolic pathway of microbiomics and metabolomics

Based on Illumina NovaSeq/HiSeq high-throughput sequencing platform, Whole Genome Shotgun (WGS) strategy was used to extract the total DNA of the macro-genome of the obtained colonies, or the cDNA double strand of the macro-transcriptome synthesized by using mRNA as template. The coding regions of prokaryotic microbial and macrogenomic gene sequences were predicted to obtain the corresponding gene sequence files, protein sequence files for protein sequence species and functional annotation. As shown in Figure 8A and Table S1, CO₂ treatment significantly reduced metabolic processes, such as metabolism (limonene and pinene degradation, energy metabolism, xenobiotics biodegradation and metabolism, C5-branched dibasic acid metabolism, β-alanine metabolism, monobactam biosynthesis, caprolactam and styrene degradation ), cellular processes (bacterial chemotaxis and biofilm formation - pseudomonas aeruginosa), environmental information processing (ABC transporters), and genetic information processing (sulfur relay system).

The metabolic activities of environmental information processing (ABC transporters), metabolism (biosynthesis of cofactors, pyrimidine metabolism, arginine and proline metabolism, biosynthesis of amino acids, carbon metabolism, phenylalanine metabolism, 2-oxocarboxylic acid metabolism, and pantothenate and coa biosynthesis), and organismal systems (protein digestion and absorption) were significantly stronger in the control group than in the CO₂-treated group, and microbial metabolism and cellular processes altered raw material metabolites (Figure 8B and Table S1). As shown in Table 1 of Figure 8C, shared metabolic pathways identified by the macrogenome and metabolome, mainly included Cellular Processes, Environmental Information Processing, Genetic Information Processing, Metabolism, and Organismal Systems. A total of 61 shared metabolic pathways involving 53 DMs were screened. Among them, ABC transporters, pyrimidine metabolism, arginine and proline metabolism, and phenylalanine metabolism metabolic processes had the most significant effects on metabolites, involving metabolites, such as pyruvic acid (C00022), glutamic acid (C00025), 5'-CMP (C00055), arginine (C00062), 2-propenoic acid (C00074), and phenylalanine(C00079).

4. Discussion

The diversity of microbial communities in raw milk is usually considered as one of the key factors affecting milk quality, because udder surfaces, milking instruments, containers and environment are potential sources of microorganisms, which leads to microorganisms being inevitable in raw milk [18]. Milk being a nutrient rich liquid is rich in water, proteins, fats, lactose, etc. which makes it a good medium for microbial growth and multiplication [19]. These microorganisms include bacteria, molds and yeasts, among others. Firstly, it was found that Stenotrophomonas maltophilia, Lactococcus lactis and Chryseobacterium sp. NEB161 were the dominant strains of raw milk, with a high percentage of 81.7%. The cold storage process altered the microbial colony diversity of raw milk, and the dominant strains at 6 days were Acinetobacter guillouiae, Pseudomonas fluorescens, Serratia liquefaciens, Pseudomonas simiae and Leuconostoc mesenteroides, with a percentage of 56.2%. These results also verified that the microbial community is a key factor in the spoilage of raw milk. Therefore, it is necessary to use appropriate techniques to inhibit the growth and multiplication of microorganisms in raw milk. Several researchers have shown that adding CO₂ to food reduced the rate of food spoilage and the growth of disease-causing microorganisms [19]. We found that CO₂ treatment significantly inhibited the propagation of the dominant colonies Acinetobacter guillouiae, Pseudomonas fluorescens, Serratia liquefaciens and Pseudomonas simiae in the CO₂-treated group as compared with the control group. The percentage in the CO₂ treated group was only 0.2%. CO₂ treatment significantly improved the storage period of raw milk (Figure 1), indicating that CO₂ improved the quality of raw milk by suppressing microbial levels and diversity.

Through microbiomics analysis, CO₂ was found to prolong the storage period of raw milk by inhibiting microbial reproduction, but the exact mechanism of its antimicrobial effect is still unclear [12]. CO₂ is highly soluble in water and lipids, and can form carbonic acid in aqueous solutions, which can lead to changes in pH in the suspension medium, and a decrease in pH inhibits the growth and metabolism of certain microorganisms [20]. CO₂ significantly lowered the pH of raw milk (Figure 1). Another view is that oxygen is a necessary electron acceptor in the metabolism of aerobic bacteria and is a specific electron acceptor for parthenogenetic anaerobes. Supplementation of raw milk with CO₂ reduced the oxygen content of the system, which in turn inhibited bacterial metabolism and reduced growth rates [21]. CO₂ treatment significantly inhibited the levels of Acinetobacter guillouiae, Pseudomonas fluorescens, Serratia liquefaciens and Pseudomonas simiae (Figure 3), which inhibited Gram-negative aerobic bacteria, suggesting that CO₂ inhibited the growth and reproduction of aerobic microorganisms and prolonged the shelf life of raw milk by reducing the oxygen content of the system.

Microorganisms in raw milk systems can utilize the nutrients in milk for growth and reproduction while breaking down substances, such as proteins, fats and lactose. For example, the pH of milk decreases when lactose is broken down to lactic acid [22,23]. Proteins are broken down, producing toxic substances, such as indole, sulfur and fecal deodorant, as well as foul odors. To be precise, the growth and reproduction of microorganisms is the metabolic process, which leads to changes in the metabolites of the raw milk system through a variety of metabolic processes, reducing the stability of the raw milk system and even spoilage [23,24]. It is worth noting that microbial effects on metabolites are mainly metabolic, including pyrimidine metabolism, arginine and proline metabolism, phenylalanine metabolism, aminoacyl tRNA biosynthesis, glycine, serine and threonine metabolism, lysine degradation, pantothenate and coenzyme A biosynthesis, and pyruvate metabolism (Table 1).

Pyrimidine metabolism refers to the process of pyrimidine nucleotide synthesis and degradation in the cell and is an important component of nucleic acid synthesis. Through pyrimidine metabolism, cells are able to synthesize and maintain sufficient levels of pyrimidine nucleotides to support DNA and RNA synthesis and the stability of genetic material [25,26]. Microbial reproduction requires replication of genetic material, and reduced levels of metabolites, such as cytidylic acid, orotate and uridine, inhibit microbial growth and reproduction (Table S2). In addition, pyrimidine nucleotide degradation produces uracil, which in turn is converted to β-aminobutyric acid, which is involved in amino acid and energy metabolic pathways. The inhibition of energy production results from reduced levels of orotate [27]. The accumulation of N-carbamoyl-L-aspartate, a key compound in the metabolism of alanine, aspartic acid and glutamic acid, suggested that the processes involved in the synthesis of amino acids and the metabolism of other nitrogen compounds were inhibited.

Arginine and proline metabolism can generate nitric oxide (NO) and polyamine compounds, which play important roles in physiological processes, such as cell signaling, immunomodulation and blood flow regulation [28]. In addition, proline metabolism produces polyamines and cofactors required for post-translational modification of proteins, which play key roles in cell growth and development [29]. Phenylalanine metabolism produces tyrosine, and other important bioactive substances, which play important roles in neurotransmitter, hormone and pigment synthesis. Glycine, serine, and threonine metabolism involves a variety of enzymes and intermediates that have multiple effects on metabolites. For example, glycine metabolism produces alanine, which is involved in glucose and fatty acid metabolic pathways. Serine metabolism produces tryptophan. Threonine metabolism produces threonine and methionine, the compounds that play important roles in methionine metabolism and thionine metabolism [30]. Amino acid metabolic processes are attenuated, but accumulation of betaine occurs in the glycerophospholipid metabolic pathway of glycine, serine, and threonine metabolism, which may be due to enhanced upstream glycerophospholipid metabolic activity. Lysine degradation generates intermediates, such as malondialdehyde, acetaldehyde and aminobutyric acid, which can be further involved in oxidation that can alter the stability of the raw milk system [31]. Reduction in glutarate and N6, N6, N6-trimethyl-L-lysine levels were found in the lysine degradation pathway. Generally speaking, it was found that CO₂ treatment significantly inhibited amino acid metabolism and reduced microbial biological activity and acid accumulation. In addition, metabolic processes, such as carbohydrate metabolism, lipid metabolism, cell growth and death, and energy metabolism were also inhibited. Overall, CO₂ treatment effectively reduced the metabolism of microorganisms and thus inhibited their growth and reproduction.

5. Conclusions

In this study, it was found that CO₂ treatment extended the shelf life of raw milk from 6 days to 16 days. This is valuable for the storage and freshness preservation and quality control of raw milk. Carbon dioxide mainly inhibited Gram-negative aerobic bacteria, such as Acinetobacter guillouiae, Pseudomonas fluorescens, Serratia liquefaciens and Pseudomonas simiae, and their growth and multiplication accelerated ABC transporters, pyrimidine metabolism, amino acid metabolism, carbohydrate and lipid metabolism, and metabolism-produced metabolites can lead to raw milk quality deterioration.