## 3 Results
### VOC Profile and Correlation Analysis  
The detected VOCs throughout storage at 0, 4 and 8 \({}^{\circ}\)C mainly consisted of esters, alcohols, ketones and aldehydes, while the presence of other compounds (e.g., limonene, ethane, 1,2-diethoxy and n-nonane) was also noticeable (Table S1). A total of 92, 60 and 78 VOCs were detected at temperatures of 0, 4 and 8 \({}^{\circ}\)C, respectively (without taking into account terpen, aromatic compounds and hydrocarbons that were found in trace amounts). Several--56, 26 and 41 VOCs at 0, 4 and 8 \({}^{\circ}\)C, respectively--were produced in inoculated models only (e.g., methanethiol, 2-methyl-1-butanol, isopropyl alcohol, dimethyl sulfide, dimethyl disulfide, 2-methylbutyratealdehyde, 3-methylbutyratealdehyde, 2-pentanone, 2-heptanone, n-nonanone, 2-undecanone, 2-nonanone, 2-nonanone, 3-pentanone, butyl butyrate, isobutyl, ethyl-2 methylbutyrate, ethyl isovalerate, ethyl isovalerate, ethyl isobutyrate, ethyl lostanoate, ethyl propionate, ethyl lgiate, lauric aldehyde, etc.). The rest of the VOCs (36, 34 and 37 at 0, 4 and 8 \({}^{\circ}\)C, respectively) were found in both inoculated and non-inoculated models (e.g., 2-ethyl-1-hexanol, 1-penten-3-ol, 2-hexen-1-ol (trans), nonanal, n-decanal, heptanal, octanal, hexenal, 6-methyl-5-hepten-2-one, acetophenone, etc. ), respectively (Table S1). At the end of storage, the VOC profiles differed to some extent among the three studied temperatures, while in some cases, it was strain dependent. According to HCA, a group of VOCs exhibited a stable profile--always detectable and closely related with either initial or final stages of storage--depending on the strain and storage temperature (Figure 1). Another group of VOCs exhibited a quite unstable tendency, since many of them were detectable sporadically at intervals. In the non-inoculated FJA in general, the majority of VOCs that were found initially (D 0) were also closely linked to this treatment throughout the storage period, with specific exceptions (e.g., 1-decanol, n-valeraldehyde and some other alcohols and ketones). Furthermore, there was no compound that was produced solely in the non-inoculated models during the latest stages of storage in any of the studied temperatures. According to the heatmap plot at 0 \({}^{\circ}\)C (Figure 1a), the models could be divided into two main clusters, clearly separated based on storage time. More specifically, the first cluster included both single and co-cultured FJA models at the initial and middle stages of storage. The non-inoculated models, including all storage days, were also a part of this group. The first cluster was mainly linked with the production of alcohols (e.g., ethyl alcohol, isopropyl alcohol, 1-penten-3-ol, 2-hexe-1-ol (trans), etc. ), aldehydes (e.g., nonanal, hexenal, heptanal, octanal, hexanal, etc.) and ketones (e.g., 2-octanone, 4-methyl-2-pentanone, etc.). The second clustered group included samples after the middle of storage (D 16 and D 20). The samples belonging to the second group exhibited a different VOC profile, since they were closely related with esters such as ethyl isovalerate, ethyl isobutyrate, isobutyl? isobutyrate, ethyl-2-methylbutyrate, dimethyl sulfide, butyl acetate and ketones such as 2- heptanone and 2-nonanone. Apart from those VOCs well-shared between the different treatments, it is crucial to highlight a unique and/or strong correlation profile between specific VOC and strains at D 20. For instance, lauric aldehyde, isobutyraldehyde and hexyl acetate were exclusively linked with strain 55 (Bacterium SBF-B18-dsp: _P. fragi_), while 2-methylbutyrateade and 4-methyl-3-penten were solely linked and strongly correlated, respectively, with strain 83 (_Pseudomonas_ sp. GSB-b11). Additionally, strain 35 (Bacterium SBF-B5-d0: _P. fluorescens_) was strongly linked with 3-octanone and 2-nonanol, while the mixed culture was the only case linked with 1-octen-3-ol, trans 2-octenal and 1-hexanol at the end of the storage. Figure 1: _Cont_.