Next Article in Journal
Ethanol Extracts from Torreya grandis Seed Have Potential to Reduce Hyperuricemia in Mouse Models by Influencing Purine Metabolism
Previous Article in Journal
A Review of the Functional Characteristics and Applications of Aristotelia chilensis (Maqui Berry), in the Food Industry
Previous Article in Special Issue
Exploring the Antibacterial Potential and Underlying Mechanisms of Prunella vulgaris L. on Methicillin-Resistant Staphylococcus aureus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metataxonomic Identification of Microorganisms during the Coffee Fermentation Process in Colombian Farms (Cesar Department)

by
Carmenza E. Góngora
1,*,
Laura Holguín-Sterling
2,
Bertilda Pedraza-Claros
2,3,
Rosangela Pérez-Salinas
3,
Aristofeles Ortiz
4 and
Lucio Navarro-Escalante
1,†
1
Department of Entomology, National Coffee Research Center, Cenicafé, Manizales 170009, Colombia
2
Department of Post Harvesting, National Coffee Research Center, Cenicafé, Manizales 170009, Colombia
3
Cienciaudes Research Group, Faculty of Health Science, Universidad de Santander, Valledupar 200004, Colombia
4
Department of Plant Physiology, National Coffee Research Center, Cenicafé, Manizales 170009, Colombia
*
Author to whom correspondence should be addressed.
Current address: Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712, USA.
Foods 2024, 13(6), 839; https://doi.org/10.3390/foods13060839
Submission received: 11 January 2024 / Revised: 29 February 2024 / Accepted: 3 March 2024 / Published: 9 March 2024

Abstract

:
The metataxonomic diversity and microbial composition of microorganisms during the coffee fermentation process as well as their relationship with coffee quality were determined across 20 farms in the department of Cesar, Colombia, by sampling coffee fruits from Coffea arabica; Var. Castillo General®, Var. Colombia, and Var. Cenicafé 1. In each farm, the fruits were processed and the fermentation process took place between 10 and 42 h following this. Three samples of mucilage and washed coffee seeds were collected per farm during the fermentation process. The microorganisms present in the mucilage were identified using metataxonomic methods by amplifying the 16S rRNA gene for bacteria and ITS for fungi. The microorganisms’ morphotypes were isolated and identified. The analysis of bacteria allowed for the identification of the following genera: Gluconobacter, Leuconostoc, Acetobacter, Frateuria, Pantoea, Pseudomonas, Tatumella, and Weisella, as well as unclassified enterobacteria; the Lactobacillacea and Secundilactobacillus families were only identified in the Var. Cenicafé 1. For fungi, the top 11 genera and families found included Hanseniaspora, Candida, Meyerozyma, Wickerhamomyces, Pichia, f-Saccharomycodaceae, f-Nectriciae, unclassified fungi, and Saccharomycetaceae, which were only found in Cenicafé 1. A total of 92% of the coffee samples obtained scored between 80.1 and 84.9, indicating “Very Good” coffee (Specialty Coffee Association (SCA) scale). Farms with the longest fermentation times showed better coffee attributes related to acidity, fragrance, and aroma. During coffee fermentation, there is a central microbiome. The differences between the microorganisms’ genera could be influenced by the coffee variety, while the specific conditions of each farm (i.e., altitude and temperature) and its fermentation processes could determine the proportions of and interactions between the microbial groups that favor the sensory characteristics responsible for coffee cup quality.

1. Introduction

The coffee produced in Colombia is internationally recognized as a high quality, mild, washed arabica coffee due to the varieties cultivated, the climatic diversity in the cultivated areas, the production systems, and the harvesting and processing practices during post-harvest [1]. For this reason, since 2005, efforts have been made to establish an insignia, such as a Protected Designation of Origin, to guarantee consumers of the quality characteristics associated with the natural factors and traditional practices of the region in which the coffee is grown. Coffee consumers currently perceive the origin and type of processing as an additional attribute of the product that can directly influence purchase preferences [1], which is why it is essential to know the characteristics of the processes that influence the quality of coffee in certain places that are highly preferable by consumers [2].
The department of Cesar has an area of 22,905 km2 [3]; it is located in northeastern Colombia between 07°41′16′′ and 10°52′14′′ north latitude and 72°53′27′′ and 74°08′28′′ west longitude. There are two major mountain systems that are important for the country in terms of area and the ecosystem services they offer to the communities [3] where the coffee is grown: La Sierra Nevada de Santa Marta, which is located to the northwest of the department, and La Serranía del Perijá, which is a bordering mountain system between Colombia and Venezuela. The coffee growing process in Cesar plays a fundamental role in the agricultural and economic sectors due to the generation of both direct and indirect jobs. In this region, coffee growing is considered multi-ethnic, with the native populations Kankuamo, Kogui, Yupa, and Arhuaco being the main coffee growers [4].
The coffee from this area is characterized by a good sensory profile, high yield, and a Sierra Nevada denomination of origin with an organic seal [5]. For this reason, understanding the post-harvest processes, particularly fermentation, that the microorganisms are naturally involved in, the physical–chemical changes during the process, and the quality are important to generate added value and recognition of coffee growing in the region.
One of the most significant factors that affects the quality of the coffee drink in wet fermentation processing is the interactions that occur between pulped coffee and the microorganisms present during the fermentation process [6]. Various studies have been carried out worldwide to analyze microbial communities and their relation to different fermentation methods and their effects on the final quality of the coffee [7,8]. Coffee variety has also been considered [9,10], as well as fruit maturity and altitude [11,12], which are factors that affect microbial composition and, therefore, interfere with the biochemical changes associated with the metabolism of microorganisms.
The fermentation process plays a fundamental role in the degradation of the mucilage attached to coffee beans. The duration of the fermentation process in the farms will generally depend on the altitude and temperature of the farm [13]. Fermentation is performed by a diversity of native microorganisms present in the coffee fruit or by those who intervene in the process and are associated with the environment, water, air, machinery, containers, plant particles, or insects that are part of the diversity of microorganisms and produce significant changes due to the concentration of primary and secondary metabolites [6,14]. According to Puerta Quintero et al. [12], all the processes that occur in the fermentation stage are generated from the interaction of different factors; among the most relevant are the intrinsic quality of the coffee variety, the pH of the substrate, the microbial load, and the type of fermentation system. The quality of the coffee beverage is the result of both the intrinsic and extrinsic properties and factors from cultivation to post-harvest activities [15].
Additionally, in Colombia, recent interest in understanding the dynamics and diversity of microbial populations during coffee fermentation processes has resulted in the identification of predominant groups of bacteria and yeasts in most processes, as is the case for lactic and acetic acid bacteria (LAB-AAB) [16,17,18]. Due to their enzymatic potential, they degrade and/or ferment the substrates contained in the mucilage, favoring the production of primary, intermediate, and secondary metabolites, while consuming sugars, amino acids, lipids, citrate, lactate, pyruvate, and alcohol [19]. The enzymes involved, such as citrate permeases, citrate lyases, oxaloacetate-decarboxylases, pectin-lyases, polygalacturonases, and alcohol and aldehyde dehydrogenases, are those that allow for the formation of intermediate products in addition to the production of lactic, acetic, succinic, and malic acids, glycerol, butanediol, and acetoin [20], which provide differentiating properties to coffee quality. However, yeasts also play an important role in both the fermentation process and in the development of flavors and aromas in the final drinks [21,22,23] because they hydrolyze macromolecules, generating aroma precursors and reducing sugars, amino acids, and chlorogenic acids [24,25].
A previous study [10] characterized the microbiota of fermentation in three Colombian coffee varieties in Pueblo Bello, Cesar, showing differences in dominant microbial groups that influence the sensory characteristics of each variety. Other works have been carried out in the department of Nariño in the south of the country, in which De Oliverira Junqueira et al. [17], through molecular methods, identified the lactic acid bacterium Leuconostoc sp. as the most predominant during the entire fermentation process (0–48 h), and Pichia nakasei as the main yeast species during the process. To the north of the country, in a study by Cruz-O’Byrne et al. [16], Leuconostoc sp. was also found to be the main bacteria involved in the process, while the Saccharomycetaceae family was reported as the most abundant yeast, and Kazachstania humilis was identified for the first time in the coffee fermentation process. In the central zone of the country, Peñuela-Martínez et al. [18] also evaluated the fermentation in farms in the department of Quindío (Colombia), reporting 1349 bacterial genera, with enterobacteria being predominant at the beginning of fermentation; they mainly identified bacteria from the genus Leuconostoc, Acetobacter, Tatumella, and Lactobacillus. In fungi, genera of the Saccharomycedaceae family predominated.
The objective of this research was to characterize the microbiota present during the coffee fermentation processes in 20 farms in the department of Cesar with the usual post-harvest practices in each farm and to determine the relationships between the microorganisms, the quality control of the fermentation process, coffee variety, farm location, and coffee quality. For this, the microbial dynamics of bacteria, fungi, and yeasts were determined through culture-dependent and metataxonomic methods, and the physicochemical characteristics of the mucilage were also evaluated to finally correlate them with coffee quality and the physical and sensory characteristics of the beverage.

2. Materials and Methods

2.1. Location and Characterization of the Farms

A total of 20 farms (Fs) located in the municipalities of Agustín Codazzi, Pueblo Bello, La Paz, and La Jagua de Ibirico in the department of Cesar (Colombia) were visited. Sampling was performed in two coffee harvest periods: (1) between December 2021 and February 2022 and (2) between September and November 2022.
In each farm, the technical information of the crop and the traditional post-harvest practices carried out by the coffee grower were recorded. Information about the usual benefit process in each farm was obtained. The conditions of the fermentation process, as well as the material of the fermentation tank, the use of water during the process, and the times used for fermentation are recorded in Table S1 (Supplementary Materials).

2.2. Sampling, Characterization of the Harvest, and the Fermentation Process

In all the farms, the manual collection of the coffee fruits from the plants was carried out in plastic containers or fique sacks that were deposited on the ground while the work was finished. In each farm, once the coffee growers harvested the coffee fruits, a sample of the coffee collected from the container was taken at random to determine the quality of the harvest using the Mediverdes® tool (Agroinsumos del Café S. A, Bogotá, Colombia) and the Cromacafé® (Agroinsumos del café S. A) color chart [26,27].
These coffee fruits were transferred to the coffee mill, weighed, and then entered the wooden hopper directly from the sacks; in no case were the varieties or collections made on different days mixed. The fruits were pulped and incorporated into the containers to be fermented either with or without water for the time that each coffee grower determined on a regular basis.
On farms that used water during wet processing, a sample of the water was collected in a 250 mL Nalgene (Thermo Scientific, Waltham, MA, USA) sterile bottle. These water samples were stored at 4 °C until processing for microbiological identification.
Under the coffee processing conditions that each farmer usually performs on his farm, the day’s batch was pulped. At the end of the process, 1.2 kg of pulped coffee was collected and processed by means of experimental mucilage removal, as previously described by Holguin et al. [10]. At this point, the first samples of mucilage (MT1) and seed (ST1) contained in the equipment corresponding to T1 (zero hours of fermentation) were taken. The remaining pulped coffee was deposited in the tanks, as the coffee farmer traditionally has to initiate the spontaneous fermentation process. Accordingly, the time determined by each farmer based on their experience was recorded. In the middle of the fermentation process, the second coffee sample was taken, which was processed in the mucilage removal equipment, and the mucilage (MT2) and seed (ST2) samples were obtained. At the end of the fermentation process, the intake of mucilage (MT3) and seed (ST3) was repeated, as previously indicated. These times varied from 10 to 42 h (see Table S1).
Sampling was carried out under the guidelines of the Collection Framework Permit 01749 granted by the Environmental Licenses Authority-ANLA to the Universidad de Santander, Colombia.
In order to determine if a complete fermentation process would be carried out in terms of mucilage degradation, at the beginning of the fermentation process, the Fermaestro® tool (Agroinsumos del café S. A, Bogotá, Colombia) [28] was used, which was introduced into the total mass contained in the fermentation tank and was left there until the final fermentation time was indicated by the farmer. These data were also recorded.
In all the farms, the seed samples (ST1, ST2, and ST3) were dried in a parabolic dryer [29] until reaching a humidity between 10 and 12%; then, they were stored in Ziploc bags and taken for quality analysis at the Quality Laboratory of the Cesar-Guajira Coffee Growers Committee. All mucilage samples (MT1, MT2, and MT3) were collected in triplicate in sterile Nalgene (Thermo Scientific Waltham, MA, USA) flasks and Falcon tubes that were stored at 4 °C in a portable refrigerator (Klimber 27 L) and transferred for preservation at −40 °C and 4 °C at the University of Santander in Valledupar (Cesar).

2.3. Physicochemical Analysis

For the mucilage samples taken at three different times (T1, T2, and T3), the following variables were determined in the fermentation tanks: temperature with a digital thermometer (Brixco, Berlin, Germany), pH with a portable pH meter HI98167 (Hanna Instruments, Woonsocket, RI, USA), and Brix degrees with an HI96801 digital refractometer (Hanna Instruments). For each measurement, three repetitions were made at different points in the fermentation tanks and the average and standard deviation were calculated. The Brix degrees were also evaluated in the field once the mucilage samples were taken. For each measurement, three repetitions were made and the average and standard deviation were calculated.
The acidity was determined in the mucilage samples MT1, MT2, and MT3 with the potentiometry technique using a 0.1 N alkaline sodium hydroxide (NaOH) solution with a 50 mL burette. Each measurement was carried out in triplicate and the average and standard deviation were calculated.
The total acidity value was expressed in mg L−1 equivalent to calcium carbonate (CaCO3) [30].
For the sugar analysis (sucrose, glucose, and fructose), the mucilage samples were centrifuged at 12,600× g RCF for 5 min and 1 mL of supernatant was taken and filtered through a 0.22 µm PVDF polyvinylidene fluoride membrane and kept at 4 °C until analysis. The samples were analyzed by means of high-performance liquid chromatography (HPLC) in the same way as reported by Holguin et al. [10]. Each measurement was carried out in triplicate and the average and standard deviation were calculated.
Physicochemical data from all 20 farms were grouped into three groups according to the final fermentation time. The first group contained the farms with up to 16 h of fermentation. The group 2 farms had 18 to 20 h of fermentation and the group 3 farms had a fermentation time of between 36 and 40 h. The data average and standard deviation were calculated in order to see the differences between the groups.

2.4. Metataxonomic Analysis

For the coffee mucilage samples (MT1, MT2, and MT3), DNA extraction was performed using 2 mL of sample. After, centrifugation was carried out at 16,000× g for 5 min. On the sediment, the extraction of genomic DNA (gDNA) was carried out with the QIAGEN DNeasy Powerlyzer Powersoil Kit. DNA quantification was performed by measuring the light absorption at 260 nm using NanoDrop™ equipment (2000-Thermo Scientific™). The gDNA samples were frozen at −20 °C.
High-throughput sequencing:
For the samples of DNA, amplicons were generated and they were used to characterize the microbiome taxonomic diversity. For bacteria characterization, the 16S rRNA gene molecular marker in the variable regions V3 and V4 was PCR amplificated using the following primers: Bakt_341F: CCTACGGGNGGCWGCAG and Bakt_805R: GACTACHVGGGTATCTAATCC [31]. Fungal diversity was determined via the amplification of the ITS molecular marker using the following primers: ITS3F (GCATCGATGAAGAACGCAGC) and ITS4R (TCCTCCGCTTATTGATATGC) [32].
The amplicons from the 16S rRNA gene and ITS end (PARE) libraries were obtained using the Herculase II Fusion DNA Polymerase Nextera XT Index Kit V2. For library construction, the sequencing library was prepared using the random fragmentation of the DNA, followed by 5′ and 3′ adapter ligation. Alternatively, “tagmentation” combines the fragmentation and ligation reactions into a single step that greatly increases the efficiency of the library preparation process. Adapter-ligated fragments are then PCR amplified and gel purified.
High-throughput sequencing was performed on the Illumina MiSeq platform, generating paired-end reads (PEs) of 300 bases each. The reads were cleaned up to a Q30 quality threshold and singletons and sequences less than 200 bases in length were removed using the Cutadapt version 3.5 program [33]. Sequences were analyzed with the Mothur program version 1.44, following the standard protocol for Illumina MiSeq libraries (SOP) [34]. PEs were assembled with the Mothur make.contigs tool and then aligned to the 16SSILVA reference database (Silva.nr v138) [35]. Subsequently, the VSEARCH algorithm [36] was used to remove the chimeric sequences. Sequences from nonbacterial lineages (sequences of mitochondrial, chloroplast, archaeal, and eukaryotic origin) were removed. The dist.seqs routine was used to group the reads into operational taxonomic units (OTUs), considering a distance limit between sequences of 0.03. The data were normalized with the normalize.shared command. The phylogenetic classification of the OTUs for fungi and bacteria with taxonomic assignment at the family and genus level for both bacteria and fungi was performed with the 16S Silva database (Silva.nr v138) [35], with a threshold of 80 (boot threshold 80) using the RDP Classifier Algorithm [37]. Coverage analysis was performed to determine sequencing coverage using Mothur and the rarefaction curve was also analyzed.
With the OTU-abundance data from each sample, beta diversity was analyzed using the Bray–Curtis distance and their ordination were visualized through nonmetric multidimensional scaling (NMDS). Differences in the community structure (beta diversity) between the coffee varieties were assessed with the permutational multivariate analysis of variance PERMANOVA [38]. Alpha diversity was examined with the indices Chao1, Shannon, Simpson, and observed species. Differences in alpha diversity were assessed using Kruskal–Wallis chi-squared statistics. Beta and alpha diversity analyses were performed in the Phyloseq and Microbiome packages of the RStudio program (version 2022.07.0) (https://posit.co/download/rstudio-desktop/, accessed on 10 April 2023) [39]. Statistical tests and graphs were carried out with the same program.
Taxonomic assignment at the phylum, family, and genus levels was performed for both bacteria and fungi.

2.5. Cultivation Counting and the Isolation of Microorganisms from Mucilage and Water

The MT1, MT2, and MT3 mucilage samples were processed to determine bacteria, yeasts, and mycelial fungi under the criteria of the current Colombian regulations [40,41]. Serial dilutions were made from the coffee mucilage samples, from 10−1 to 10−4, and 100 µL was inoculated on the culture media for each microbial group; mesophilic aerobes were cultured in nutrient agar (AN) and agar plate count (APC) (OXOID, Basingstoke, UK), and enterobacteria were cultured in MacConkey agar. Incubation was carried out at 28 °C for 48 h. After incubation, the CFUs (colony-forming units) were counted in a 10−3 dilution, counting each one of the CFUs that grew in the plates.
For the determination of lactic acid bacteria (LAB), first, an enrichment of the sample was carried out by taking 5 mL of mucilage and adding it to 25 mL of Man Rogosa Sharpe-MRS broth (MERCK, Darmstadt, Germany), then incubating it at 29 +/− 1 °C for 18 to 24 h at 80 rpm in a Shaker 1000 (Heidolph, Schwabach, Germany). Subsequently, serial dilutions were made up to 10−5, then 100 µL of each dilution was plated on the MRS agar medium and incubated aerobically and anaerobically at 28–30 °C for 48 h. After incubations, the CFUs were counted in the same way as described previously.
For the isolation of AAB, enrichment broths II and III were used (glucose, yeast extract, peptone, ethanol, pH 4.5; additionally, broth II contained acetic acid and had a pH of 3.5), adding (in equal proportion) broth and mucilage to the test tubes and incubating at 30 °C for 72 h at 80 rpm in a Shaker 1000 Heidolph [42]. Subsequently, a 100 µL aliquot was taken and inoculated on GEY-CaCO3-selective medium (2% glucose, 20% ethanol, 1% yeast extract, 3% CaCO3, agar) with the addition of 0.3% spiramycin 3 MUI (Labinco). Incubation at 30 °C for 10 days was followed by biochemical tests for identification [43]. After the 10 days, the CFUs were counted in the same way as previously described.
In order to grow and count the group of yeast and mycelial fungi Potato dextrose agar (PDA) (OXOID, Basingstoke, UK), yeast extract peptone dextrose agar (YEPD) (MERCK, Darmstadt, Germany) and Sabouraud 4% dextrose Agar (SDA) (MERCK, Darmstadt, Germany) were used. The incubation was carried out at a temperature of 28 °C for 96 h; after that, the CFUs were counted.
The microbial count data from all 20 farms were grouped into three groups according to the final fermentation time. The first group contained the farms that had up to 16 h of fermentation. The group 2 farms had 18 to 20 h of fermentation and the group 3 farms had fermentation between 36 and 40 h. The average and standard deviation were calculated in order to see the differences between the groups.

2.6. Identification of Microbial Groups from Coffee Mucilage

The identification of enterobacteria was carried out through the phenotypic criteria of colony-forming units (CFUs) and microscopic characteristics from the isolated bacterial strains. Additionally, biochemical tests were used, such as catalase, oxidase, OF-glucose, motility, and growth on MacConkey agar. As a complementary step, 1 CFU was taken, diluted in 1 mL of 0.85% saline solution, and inoculated into the API 20E gallery (BIOMERIEUX, Craponne, France). A dilution process of pure CFUs in 0.85% saline solution was carried out to identify spore-forming, Gram-positive bacilli. This was followed by another dilution on a 2 McFarland scale from which twice the quantity was extracted and added to the API50CHB/E vial. Subsequently, the gallery was inoculated. Furthermore, 12 tests from the API20E gallery were performed, as follows: β-galactosidase (ONPG), arginine dehydrolase (ADH), lysine decarboxylase (LCD), ornithine decarboxylase (ODC), citrate utilization (CIT), hydrogen sulfide production (H2S), urease activity (URE), tryptophan deaminase (TDD), indole production (IND), Voges–Proskauer test (VP), gelatinase production (GEL), and nitrite reduction (NO2).
Catalase, oxidase, and coagulase tests were carried out to identify spore-forming, Gram-positive bacilli.
For the identification of the LAB group after incubation, the total count was performed and each morphotype was characterized via macroscopic observation and microscopy. Biochemical discard tests were also performed, as follows: catalase, oxidase, coagulase, nitrate reduction, gelatin, mobility, and indole. Finally, the API technique using 50CH and API 50CHL (Biomerieux®), which contain 50 biochemical tests that assess assimilation, oxidation, and fermentation, was performed. Identifications were made with APIWEBTM software.
For the AAB group after the incubation period, the characteristics of the morphotype were observed and a translucent halo formed around the CFUs. Smears and Gram staining were performed to confirm the morphology of the bacilli, as follows: medium and short Gram-negative, catalase tests, cytochrome oxidase, acetate oxidation, growth in 3% NaCl and D-mannitol broth, growth in glutamate agar, and growth in sugars such as L-arabinose, L-rhamnose, D-mannose, sucrose, and maltose [42,43]. Through these tests, the identification of the genera was achieved.
The identification of yeasts and mycelial fungi was carried out via the characterization of the isolated morphotypes by means of lactophenol blue staining and the observation of their microscopic characteristics. Then, the isolates were identified using the microbial DL-96II ID/AST system. Yeast colorimetry and turbidimetry were performed via the semi-quantitative analysis of antimicrobial MICs, following the manufacturer’s instructions (Zhuhai DL Biotech Co., Ltd., Zhuhai, China).
Mycelial morphotypes were identified from the development of colonies on PDA, SDA, and YEDP agar. The macroscopic recognition of the colonies (front and back), as well as their appearance (shape, size, color, and texture), was performed by reviewing the taxonomic keys [44].
Microscopic identification was carried out using the microculture technique via observing reproductive structures and using the taxonomic key of imperfect fungi [45].
Serial dilutions were carried out for the microbiological identification of the water samples. The most likely number (MPN) technique allows for the identification of total coliforms in Bright Green Broth according to the regulation. The number of mesophilic aerobes was also determined. The morphotypes were identified based on the macroscopic and microscopic characteristics and Gram staining [46,47]. Other bacteria were identified by plating in plate count media and eosin methylene blue or EMB agar using conventional biochemical tests that allowed for the isolation of different genera.

2.7. Physical and Sensory Analysis of Coffee and Coffee Cup Quality

Physical and sensory analyses were performed on the 60 seed samples (ST1, ST2, and ST3) from the 20 farms after they were dried and reached a humidity of between 10 and 12%.
The analysis of the physical quality of the green coffee beans was based on the determination of moisture, the percentage of loss, low quality, black and vinegar beans, beans infested with coffee borers, and the percentage of healthy beans according to what is established by the Colombian Institute of Technical Standards and Certification [48]. The sensory evaluation was carried out according to the Specialty Coffee Association (SCA) protocol with the participation of three Q-Grader-certified tasters by the CQI (Coffee Quality Institute), which belongs to the sensory panel of the Quality Laboratory of the Cesar-Guajira Coffee Growers Committee [49]. With this methodology, 10 coffee sensory attributes were recorded, as follows: fragrance/aroma, flavor, residual flavor, acidity, body, balance, uniformity, clean cup, sweetness, taster score, defects, and total, as well as sensorial, descriptors. The sensory quality expressed as the SCA total score was the response variable. The total score was the sum of each attribute score. The SCA total score and the coffee attributes data from all 20 farms were grouped into three groups according to the final fermentation time. The first group contained the farms that had up to 16 h of fermentation. The group 2 farms had 18 to 20 h of fermentation and the group 3 farms had fermentation between 36 and 40 h. The average and standard deviation were calculated in order to see the differences among the groups.
By using the Voyan Tools software [41] the coffee sensorial descriptors were analyzed in a word cloud.

3. Results

3.1. Characterization of the Harvest and Fermentation Processes

Table S1 in the Supplementary Materials shows the characteristics of the 20 sampled farms in the 4 municipalities of Cesar. Regarding the varieties of coffee planted on the farms, 55% corresponded to Var. General Castillo, 40% to Colombia, and one farm, F7, had Var. Cenicafé 1.
The farms were located between 1.116 and 1.938 m above sea level (m.a.s.l.). A total of 30% of the farms were located at an altitude between 1.116 and 1.390 m.a.s.l., 40% were at an altitude of 1.400–1.700 m.a.s.l., and the other 30% were between 1.740 and 1.938 m.a.s.l. A total of 50% of the farms evaluated did not use water for coffee processing, another 45% did use water during the process, and one farm (5%), farm F7, carried out submerged coffee fermentation, which consists of adding water to completely cover the pulped coffee mass. This measure was performed visually by the coffee farmer.
According to what was stated by the coffee growers, the processes in each farm are mainly established by the operational capacities of the mill and the coffee production on the farm; that is, the greater the amount of coffee collected, the shorter the fermentation times since they need to facilitate the spaces in the mill to process the next day’s collection, as was the case for farms F12 and F14, with 14 and 10 h of fermentation. In contrast, farms with a smaller amount of coffee collected prefer to leave their coffee fermenting for several days, as was observed in farms F6, F7, and F11, with fermentations between 36 and 42 h. Regarding the materials of the fermentation tanks, 50% of the tanks in which the fermentation took place were tiled, 35% were concrete tanks, and 15% were plastic tanks.
The characterization of the harvest with the Mediverdes® tool in the 20 farms is presented in Table S2 in the Supplementary Materials. According to the classification of the tool, 75% of the farms had an excellent harvest and 20% had a good harvest. In 5% (F18), an average harvest was obtained.
There were some farms with incomplete fermentation according to the indicator Fermaestro®. In this case, most of the farms correspond to the shortest fermentation time. Farms F12, F14, and F15 had fermentation between 10 and 16 h and two farms had 18 h of final fermentation (F18 and F20).

3.2. Physicochemical Analysis

3.2.1. Temperature in the Fermentation Tanks

The average temperature (Figure 1) in the fermentation tanks during the fermentation process ranged between 18.9 and 26 °C. At the end of the process, the highest temperature values were obtained in farms F11, F13, F1, and F9, ranging from 34 to 28.8 °C. A few hours after the start of the fermentation process, in some farms, a slight decrease in temperature was observed, mainly associated with environmental conditions, as these evaluation points were usually carried out at night. On average, 20.5 ± 1.93 °C was observed at 9 h of fermentation, 23 ± 3.12 °C at 18 h, and 26.3 ± 6.89 °C at 36 h of fermentation. Although the fermentation time varied depending on the farm, in general, an increase in temperature was observed in all fermentations between 1 and 10 °C. The farms where the temperature increased the most were farms F11, F17, and F1, with increments of 10.6, 8.4, and 6 °C, respectively. Furthermore, those that remained more stable were farms F2, F4, F7, and F8, with averages of 22.0 ± 0.01, 20.2 ± 0.01, 20.1 ± 1.1, and 20.8 ± 0.81, respectively.

3.2.2. pH, Brix, and Total Acidity

The changes in pH, Brix, and total acidity during fermentation in the 20 evaluated farms can be seen in Table S3 in the Supplementary Materials.
Additionally, Figure 2 shows the results of the average pH and total acidity in the farms, grouped into three groups according to the final fermentation time.
For pH, the three fermentation groups showed a significant decrease between T1 and T3 at the end of the process. At the end of the process, independent of the fermentation time, the final pH was similar and around 3.8. In the farms with the longest fermentation time, this value was reached around 18 h, as well as in the other two groups of farms.
The Brix degree measurements show that in 80% of the samples (16 farms), a decrease in the dissolved solids and sugars was evidenced at the end of fermentation. On average, the decrease was clearly observed in the three groups and the behavior was very similar to that observed with the pH.
For total acidity, in the three groups, an increase was observed and the highest acidity value was reached in the group with the longest fermentation time, up to 42 h. At the beginning of the process, the values recorded were between 4.1 and 333.4 mg L−1 CaCO3. Particularly, the most notable increases in acidity occurred in farms F1 and F20, with values at the beginning of the process of 248.4 and 10.1 mg L−1 CaCO3 and ending the fermentation with 741.7 and 314.7 mg L−1 CaCO3, respectively.

3.2.3. Sucrose and Reducing Sugars

The percentages of sucrose and reduced sugar in the mucilage are shown in Table S3 in the Supplementary Materials. Figure 3 shows the results of the average sugar data from all the 20 farms, grouped into the three groups according to the final fermentation time.
At the beginning of the fermentation process, the sucrose contents were between 0 and 0.75%, and, in 19 farms, the sucrose content decreased at the end of fermentation and the final total content was similar in the three groups. In F18, there was an increase from 0.04 to 0.23%.
The analysis results of reducing sugars composed of glucose and fructose monosaccharides in mucilage can be seen in Figure 3. The final total content was similar in the first and second farm groups, while the farms with longest fermentation time showed the lowest values and largest reduction. Particularly, in F18, an increase in reducing sugars was observed. On the other hand, the farms with the highest values of reducing sugars in mucilage both at the beginning and at the end of the fermentation process were F4 and F14.

3.3. Microbiologic Analysis

3.3.1. Microorganisms Count

Table S4 in the Supplementary Materials shows the counts in terms of CFUs obtained for mesophiles, LAB, coliforms, yeasts, and filamentous fungi in the mucilage.
Table 1 shows the microorganisms count data, grouped into three groups according to the farm’s final fermentation time.
The LAB counts, in general, were higher than those obtained for the other groups, with an average of 6.29 log10 CFUs/mL, followed by yeasts with an average of 5.75 log10 CFUs/mL, mesophiles at 5.47 log10 CFUs/mL, coliforms at 3.92 log10 CFUs/mL, and, to a lesser extent, filamentous fungi at 2.8 log10 CFUs/mL.
In general, regardless of the group of microorganisms, there were more farms in which the counts of all groups of bacteria and fungi decreased as the fermentation time elapsed. In the fermentation time of up to 16 h, in the five groups of microorganisms, a decrease is observed between T1 and T3. In the group of up to 20 h of fermentation, only for mesophiles is no decrease apparent. Moreover, in the group of farms with fermentation times between 36 h and 42 h, there is an increase in microorganisms belonging to the mesophile and yeast groups.
From time 1 (T1) to time 3 (T3), a decrease in mesophile and yeast numbers was observed in 12 farms (60%); in the LAB group, a decrease was observed in 15 farms (75%). In three farms (15%), filamentous fungi were not observed in the counts (F7, F8, and F13) and of the remaining 17 farms, only one showed a slight increase in this group from T1 to T3 (F11).
Mesophiles in F6 and F9 were the predominant group, with averages in each farm of 6.19 and 5.94 log10 CFUs/mL, respectively. In both farms, the population increased from the beginning to the end of fermentation.
On the other hand, in yeasts, a decrease in the population was observed in 11 farms, and it was the predominant microbial population in F1, F2, and F3, respectively. In terms of filamentous fungi, these populations remained the same or decreased in 14 farms. Additionally, in F7, F8, and F13, they were not detected at any stage of fermentation.
Individually, farm F7 had the highest LAB count at T1 with 9.51 log10 CFUs/mL; for yeast and mesophiles, farm F6 at T3 showed the highest counts with 6.41 log10 CFUs/mL and 6.60 log10 CFUs/mL, respectively.
The development of coliforms and mycelial fungi in 55% of the farms decreased as fermentation time elapsed.

3.3.2. Bacterial Metataxonomy Analysis

The analysis of sequencing coverage and rarefaction analysis with the 16S rRNA gene molecular marker presented percentages of greater than 96% (Table S5 in the Supplementary Materials), indicating that the sampling was sufficient and presented the appropriate characteristics to study the diversity in the samples. Likewise, it was possible to calculate the alpha diversity indices. Observed genus species or OTUs presented values between 172 and 1008, being higher in all farms at the beginning of fermentation. Farms F8 and F14 presented the highest values at the beginning of fermentation, with 1008 and 1002 OTUs, respectively, and at the end of the process, farms F10 and F17 had the highest OTUs with values of 436 and 457. No statistical differences for “Observed genus species” between coffee varieties (Kruskal–Wallis chi-squared = 1.5179; df = 2; p-value = 0.4681) were found.
The ACE richness estimator decreased from the beginning to the end of fermentation in the 15 farms, while an increase in wealth was observed on farms F8, F9, F10, F11, and F12. No statistical differences in the ACE index between coffee varieties (Kruskal–Wallis chi-squared = 1.2978; df = 2; p-value = 0.5226) were observed. For the Shannon diversity index, it was possible to show a decrease in all farms except F4; likewise, for the Simpson index, an increase was observed in farms F5, F15, and F16, while in the others, a decrease in diversity was identified from the beginning to the end of the process, particularly in F18. No statistical differences in the Shannon index between coffee varieties (Kruskal–Wallis chi-squared = 3.0791; df = 2; p-value = 0.2145) were observed. The Simpson index was the same at the beginning and at the end of the fermentation process, with a value of 0.90. No statistical differences in the Simpson index between coffee varieties (Kruskal–Wallis chi-squared = 1.8109, df = 2, p-value = 0.4044) were observed.
Beta diversity analysis was carried out using nonmetric multidimensional scaling (NMDS) (Figure 4), which allowed us to identify the relationship between the values obtained for alpha diversity and the coffee variety in each of the 20 farms evaluated. There is a significant difference in the genus microbial community structure (beta diversity) between coffee varieties (PERMANOVA p-value = 0.001). Few differences were observed, as shown by the dispersion of the points between the results of microbial diversity of the Var. Castillo general and Colombia, while a distant grouping was observed in the Var. Cenicafé 1.
The top 10 bacterial taxonomic assignments at the family and genus levels, with their relative abundances are shown in Figure 5 and Figure 6. The analysis at each farm was carried out at the initial, middle, and final points of the fermentation process, which allowed us to show the changes in the communities. At the family level (Figure 5), in all the farms, the bacteria of the Lactobacillaceae family were mainly identified; to a greater or lesser extent, the bacteria of the Acetobacteraceae, Enterobacterales, and Erwiniaceae families were also identified.
At the genus level (Figure 6), the metataxonomic analysis of bacteria allowed us to find the top 10 genera, where Gluconobacter and Leuconostoc stood out. Together, they were observed in a greater proportion in all the farms and evaluated times, followed by Acetobacter, Tatumella, and Weisella, and unclassified enterobacteria were present in most of the farms; however, the proportions were different in each farm. There is an assignment identified as “other”; in this, the genera from the lactobacillaceae family were present, as was observed in the 16S.trim.contigs.good.unique.good.fileter.uniqueprecluster.pick.nr_v132.wang.tax.summary.
In the wet processing of F1 and F11, a greater predominance of Acetobacter was observed in T3. On the other hand, in F2, Leuconostoc predominated during the process, while in F3, the bacteria of the genus Weissella predominated. Secundilactobacillus was predominant only in farm 7. In F4, bacteria of the genera Pantoea and Leuconostoc predominated. In F5, mainly the unclassified bacteria of the Enterobacterales and Pantoea genera were observed, and F6 had the highest predominance of Tatumella at T2 and T3; however, it could also be identified at T3 of F15, F16, and F18. In F7, with Var. Cenicafé 1, a decrease in the Acetobacter population was observed, while a progressive increase in the Secundilactobacillus population was observed. Farms F8, F9, F10, F12, F17, F19, and F20 showed an increase in the population of AAB Gluconobacter and Frateuria during fermentation. However, the abundance of the LAB Leuconostoc also increased during the fermentation process.

3.3.3. Bacteria Isolates

Table S6 in the Supplementary Materials describes the mesophilic aerobic bacteria isolated at the three mucilage sampling points. In total, 40 different isolates corresponding to mesophilic aerobic bacteria were identified, after been found in different farms at different times (Figure 7). Gram-negative strains, such as Stenotrophomonas maltophila, stood out and were observed in 45% of the farms. In 40% of the farms, it was possible to isolate Pantoea agglomerans, followed by the presence of Bacillus, distributed as follows: Bacillus firmus in 25% and B. subtilis in 20% of the farms. The other bacteria of this group were only isolated into pure culture in three or fewer farms. Of the Gram-positive bacteria in 25% of the farms, Micrococcus spp. were observed in greater abundance in T1.
Table 2 describes the LAB; eight species were identified. A total of 100% of the farms exhibited LAB, either in the form of Lactiplantibacillus plantarum, Leuconostoc mesenteroides, or Lactobacillus delbrueckii. The farms that had L. delbrueckii did not have L. plantarum. The strains that were most frequently observed in the fermentation process of the farms were L. plantarum, followed by L. mesenteroides, Levilactobacillus brevis, and Lactiplantibacillus pentosus in the three sampling times. Additionally, the L. plantarum strain, in addition to being the most abundant on the farms (85% presence), was also observed during the vast majority of sampling times (67%). This, in turn, was one of the LAB that showed a decrease from T1 to T3. For the three remaining genera, 45% of the farms recorded the presence of L. mesenteroides. A total of 40% of the farms showed the presence of L. brevis. Likewise, 30% of the farms showed the development of L. pentosus. Other strains, with respect to their frequency on the farms, were lower, such as Leuconostoc citreum, which was present in 15%, and L. delbrueckii in only 10% of the farms.
Regarding the group of acetic acid bacteria or AAB, in 11 farms, this type of bacteria was not isolated (F1, F2, F3, F4, F5, F6, F8, F10, F12, F14, and F20) (Table 3). In total, four genera were identified, as follows: Acetobacter spp. and Gluconobacter spp., which were found in 25% of the farms. Another identified genus was Saccharibacter, which was found in two farms (10%) at both the beginning and end of fermentation, and Acidomonas was found in one farm (F9).
The results of the microbiological analyses of the water samples used for fermentation on the 10 farms showed that 15 bacterial genera were isolated; a total of 13 were Gram-negative, corresponding to enterobacteria and Pseudomonas spp., and two were Gram-positive: Staphylococcus spp. and Micrococcus spp.
The evaluation of the thermotolerant coliform in water using the most likely number technique indicated that F5 obtained the highest concentration of 460 MPNCF/100 mL, and for the mesophilic count, F10 registered a concentration of 4.72 log10/mL from AM. The main genera found were Proteus mirabilis, P. vulgaris, Enterobacter spp., and Micrococcus spp. for the first, and Citrobacter freundii, C. gergoviae, C. murliniae, and Edwarsiella tarda for the second.

3.3.4. Yeast and Fungi Metataxonomy Analysis

The analysis of sequencing coverage with the ITS molecular marker for fungi presented percentages greater than 99% (Table S7 in the Supplementary Materials), indicating that the sampling was sufficient and included the appropriate characteristics to study the diversity in the samples. Likewise, it was possible to calculate the alpha diversity indices. Observed species or operational taxonomic units (OTUs) presented values between 61 and 629, unlike what was observed in the bacterial analysis. In the fungal communities, an increase in the species observed during the fermentation process was observed, with values from 285 to 319, 367 to 378, and 108 to 110 in farms F1, F3, and F11, respectively. In the other farms evaluated, a decrease in the species was evidenced from the initial point to the end of the fermentation process. No statistical differences for observed species between coffee varieties (Kruskal–Wallis chi-squared = 3.4776; df = 2; p-value = 0.1757) were found.
The ACE richness estimator varied between 105 and 2155, decreasing in 16 farms, whereas it increased in farms F1, F4, F5, and F7. No statistical differences in the ACE index between coffee varieties (Kruskal–Wallis chi-squared = 3.3426; df = 2; p-value = 0.188) were observed.
The results of the beta diversity analysis were carried out using nonmetric multidimensional scaling (NMDS) (Figure 8), which allowed us to identify a trend similar to that observed with the bacterial communities, since there were few differences between the diversity results of the genus coffee samples from farms with Var. Castillo general and Var. Colombia, while a distant grouping was observed in the samples of farms with Var. Cenicafé 1. There is a significant difference in the microbial community structure (beta diversity) between coffee varieties (PERMANOVA p-value = 0.001). These results indicate that there were marked differences in this group of microorganisms with respect to the other two varieties.
The top 10 metataxonomic assignments of fungi at the family and genus levels and their relative abundances are shown in Figure 9 and Figure 10. The analysis at each farm was carried out at the initial, middle, and final points of the fermentation process, which allowed for evidence of the changes in the fungal populations. At the family level (Figure 9), yeasts from the Saccharomycodaceae and Phiciaceae families were mainly identified; in F1, particularly at T1, unclassified fungi of the order Saccharomycetales predominated; however, at T1 and T2, the structure of the fungal community was of the fungi of the Saccharomycodaceae and Phiciaceae families, similar to what was observed in most of the processes. On the other hand, in the fermentation processes of F4, F5, F6, and F20, unclassified or “other” fungi were identified in the different stages, which may suggest the presence of new fungal species. At the genus level (Figure 10), a predominance was observed in the Saccharomycodaceae family in 19 of the 20 farms evaluated. In contrast, in F7, the presence of the Saccharomycetaceae and Nectriaceae families was evidenced.

3.3.5. Yeasts and Filamentous Fungi Isolates

The fungi (yeasts and mycelial) were isolated and identified during the three fermentation times, which are listed in Tables S8 and S9 in the Supplementary Materials.
Figure 11 shows the type of yeasts that were isolated in the natural fermentation process in the 20 farms.
A total of 100% of the farms showed Hanseniaspora and some Candida yeasts (Candida sp., C. albicans, C. guilliermondii, C. kruesi, and/or C. tropicalis). Finally, Cryptococcus neoformans was observed in 30% of the farms and the other yeasts appeared in less than three farms; among these yeasts are Rhodotorula sp., Pichia, and Hanseniaspora.
Table S8 shows the frequencies and times at which different mycelial fungi were observed in the farms: Penicillium sp. was found in 55%, Rhizopus spp., Aspergillus fumigatus, and Aspergillus flavus in 25%, and Cladosporium sp. in 20%. However, the highest number of Penicillium sp. strains was observed in the first and second fermentation times. A similar behavior was observed in Rhizopus spp. and Aspergillus flavus, where the presence of fungi was similar at times 1 and 2, but at the end of fermentation, they decreased.
Cladosporium spp. were present in 20% of the farms. Fungi, such as Mucor sp., Trichoderma sp., and Fusarium were seen in three or fewer farms. In three farms, F7 (36 h), F8, and F18 (18 h), no filamentous fungi were identified.
The mycelial fungi, although in low percentages, were present throughout the fermentation process of the mucilage at the three time points.

3.4. Quality Analysis

Table S10 in the Supplementary Materials shows the average humidity, total coffee cup SCA score, and sensory description for the 60 samples. The average humidity of the seed samples was 11.1%, so the samples were in the acceptable ranges of drying for the cup quality analysis.
The total SCA score data are grouped into three groups, according to the farm’s final fermentation time and are shown in Table 4. In this table, the data from farms F1 (T2), F1 (T3), and F 20 (T2), corresponding to defects associated with fermentation, mold, and contamination were not considered.
The average global score was 82.7. The three groups showed coffee considered “Very Good”, according to what is established by the Specialty Coffee Association [38], ranking between 82.3 and 83.7. A total of 3% of the samples achieved a score greater than 85, which is considered “Excellent Special Coffee”, 5% of the samples achieved values lower than 80 points, categorized as “Commercial Coffee”, and 92% of the samples achieved scores between 80.1 and 84.9, while the average highest scores (83.7) were observed in the farm group with the largest fermentation time (36–42 h) at the end of the fermentation. The best scores were observed in F2, F6, and F7.
In the total SCA score from all the sensory attributed evaluated factors, fragrance, flavor, and acidity, shown in Table 5, were the ones with the higher contribution to the total score. In this case, the highest cup acidity is observed in the samples with the longest fermentation time (T3), with a score of 7.73. These farms also were the ones that showed the highest total acidity in the mucilage physicochemical analysis (Figure 2).
With respect to the sensory attributes of fragrance and flavor, these are also related to the sensory descriptors. The coffee sensorial descriptors were analyzed in word cloud.
The Q-Grader certified tasters identified a particular sensory descriptor in the coffee cup samples from the farms with fermentation between 36 and 42 h at the end of the fermentation (T3). Figure 12 shows the word cloud that describes coffee sensory descriptors from those farms. The bigger the size of the word, the higher the frequency as a descriptor in the coffee cup of the farm. The sensory descriptors blackberries, fruity, citrus, and caramel are the most predominant. The others word clouds are not shown.

4. Discussion

The department of Cesar shows a large variability regarding the characteristics of the coffee farms, altitudes, temperature, fermentation times, and the material of the tanks where the wet fermentation takes place. However, the harvesting of the crop was appropriate based on the Mediverdes® analysis. This could be a first indicator of a good cup quality of the coffee obtained.
The temperature characterization showed that the lower temperature variations seem to be an indicator of good fermentation. The farms F2 with Var. Castillo General, F6 with Var. Colombia, and F7 with Var. Cenicafé 1, and a fermentation time of 36 h achieved the best sensory quality results, and these showed the lowest temperature variation. This could indicate that a constant and stable temperature may favor the presence of certain groups of microorganisms that positively affect fermentation. The evidence may be in the metataxonomic analysis, in which the microbial populations of both the bacteria and fungi in these farms were stable at the three evaluation times. The increase in temperature during fermentation may be due to the number of microorganisms, as well as their type and growth rate, coupled with the high acidity due to the metabolites generated.
Because good quality results were obtained at the end of fermentation in farms with a constant temperature, the temperature in farm F1 (36 h fermentation) stands out, where the highest increase was observed throughout the long fermentation process and in which the defects (ferment and mold) associated with this occurred; it is possible that this high temperature favors the production of acidic metabolites and the growth of unwanted microorganisms, such as the Citrobacter and Serratia plymuthica-type enterobacteria isolated only on this farm. The genus Serratia has been isolated from water, soil, plants, and air and is known to secrete many enzymes, which are capable of inhibiting the growth of other microorganisms beneficial to the process since they degrade compounds in the substrate that are important for the development of certain microbial groups [50,51]. Therefore, it would be recommended to establish temperature control points to avoid the fermentation defects associated with the presence of unwanted bacteria. Nonetheless, additional assays under controlled conditions would be necessary to address the impact of temperature variations on quality features.
Regarding the Brix parameter, a decrease in this indicator is generated in response to the degradation of the organic acids, and the oxidation and fermentation of other substrates caused by the microorganisms. This also causes a decrease in pH that influences the various types of microorganisms present and their metabolites released during this stage, such as lactic, citric, succinic, malic, and acetic acids, ethyl acetate, acetaldehyde, glycerol, and ethanol, which can add value to the sensory characteristics of the final product [19,52]. The group of farms with the longest fermentation time, independent of the coffee variety, were also the ones with the final highest total acidity measured in the mucilage. They also had a high coffee cup SCA score at the end of the fermentation process. This score was related to the sensory attributes, as follows: acidity, fragrance, and aroma and special descriptors such as blackberries and fruity flavors.
Similar to temperature, the farms that showed lower variations in total acidity (F2 and F7) were the ones that also showed better characteristics of the drink, and the highest variations in this variable were associated with fermentation defects in farms F1 and F20. Notably, in these two farms, an accelerated growth of AAB was observed, which may explain the variations in this characteristic. Particularly in the farms with longer fermentation processes, it was possible to show stabilization in relation to pH after 18 h of fermentation, because at the end of the fermentation in the three groups of farms, the pH was very similar.
The farms with the highest values of reducing sugars in mucilage both at the beginning and at the end of the fermentation process were F4 and F14, and they also showed incomplete fermentation, as indicated by Fermaestro®. The degradation of these carbohydrates may be associated with the mechanisms used by microorganisms for mucilage degradation. In general, with respect to the sucrose consumption, small variations were observed in all the farms. However, in respect to reducing sugar in the farm group with fermentation between 36 and 42 h, the highest consumption and lower values are observed at the end of the process; this can be associate to the long fermentation time.
The number of microorganisms of the different groups isolated on the farms showed a greater presence of LAB, which coincides with the results obtained by Pereira et al. [46]. Yeasts, mesophiles, AAB, and, to a lesser extent, fungi were found, independent of the fermentation time. In contrast, Peñuela-Martínez et al. [18] also evaluated fermentation in farms in the department of Quindío, Colombia, and reported enterobacteria as the predominant group at the beginning of fermentation, and their abundance decreased over the course of the process. This indicates that the abundance of these microorganisms depends on the geographical areas. Peñuela-Martínez et al. [13] reported similar mesophilic aerobic bacteria and LAB counts with averages of 8 log CFUs/mL and yeast counts of 7.7, which are values higher than those found in our study.
In general, regardless of the group of microorganisms present, there were more farms in which a tendency was observed for the counts of all groups of bacteria and fungi to decrease as the fermentation progressed. However, the largest decrease was observed in the LAB group in the farms with fermentation between 36 and 40 h. In this group, a slight increase in mesophiles and yeast was observed. In the mycelial fungi, there was very little variation in the counts. When comparing the fermentation processes with longer times (36 and 42 h) to the fermentation processes with shorter times (10, 14, and 18 h), it could be expected that the longer the time, the higher the counts, but this was not the case. In general, this microbial count seems to be independent of the fermentation time. No major changes in the counts were identified between the various groups evaluated on each farm, and these observations are similar to what was reported by Peñuela-Martínez et al. [13].
The metataxonomic analysis of bacteria allowed us to find the top 10 genera with greater abundance, among which Gluconobacter and Leuconostoc (LAB) stood out. These were observed together in the greatest proportion in all the farms and times evaluated, followed by Acetobacter, Tatumella, and Weisella, and the unclassified enterobacteria present on most of the farms; however, the proportions were different on each farm. In the case of Secundilactobacillus sp., it should be noted that it was only found in the Var. Cenicafé. Additionally, this farm had one of the best sensory quality scores at the end of the process, with 84.0 SCA points.
In this study, Lactiplantibacillus plantarum was the species most frequently isolated in the fermentation process, together with Leuconostoc. Regarding the presence of Leuconostoc, it was also one of the most predominant bacteria during the entire fermentation (0–48 h) process in the departments of Nariño and Magdalena [16,17]. The LAB are recognized together with yeasts as the main transformers and contributors of aromatic compounds for coffee during the fermentation process by metabolizing sugars into acidic compounds, esters, ketones, aldehydes, and alcohols, among others; likewise, LAB contribute to the degradation of complex compounds, such as pectins found in high concentrations in coffee [20]. LAB groups have a metabolic activity that provides high concentrations of organic acids, such as succinic acid and lactic acid, in addition to ethanol and CO2. This activity positively influences lactic fermentation, along with hetero- and homofermentative processes in coffee; this type of fermentation is related to the delicate floral, fruity, and body perception of the drink [53].
The main genera and species identified within the Enterobacteriaceae family were Pantoea agglomerans, Citrobacter freundii, Klebsiella pneumoniae, and Escherichia coli, and similar results were obtained by Vilela et al. [54]. These correspond to ubiquitous strains in plant environments that can reach the fermentation process through external factors but have the ability to degrade the components of the mucilage through a mixed acid metabolic pathway. In the case of Escherichia coli, its presence is related to the lack of good manufacturing practices. On the other hand, the following stand out: Stenotrophomonas maltophilia, Pasteurella pneumotropica, Hafnia alvei, and Chryseobacterium indolgenes; their presence may be due to environmental factors and to the lack of application of good practices since they are found in plants, animals, and insects, among others [55]. The Gram-positive bacteria Bacillus subtilis, Bacillus firmus, Bacillus cereus, and Bacillus megaterium were identified, as previously reported in coffee fermentation by Viela et al. [48]. Regarding AAB, the genera Acetobacter spp. and Gluconobacter spp. were identified and the concentration of acetic acid generated during the coffee fermentation process may be due to the enzymatic degradation of the sugars and the alcohol released by yeasts and LAB. These compounds are used as nutrients by acetic bacteria that are strictly aerobic [56].
Regarding the bacteria identified in the processing water used by some farms, Citrobacter spp., Enterobacter spp., and Klebsiella spp. are enterobacteria that ferment sugars via a mixed acid pathway, which indicates that they have the ability to exploit the mucilage components through degradation and fermentation with the generation of different organic acids, alcohols, and gases; however, these genera are sensitive to highly acidic conditions, the secondary metabolites generated by LAB, and temperature changes. Among the most frequent species in water and mucilage found in this study are Citrobacter and Enterobacter, previously reported in [57,58].
In the metataxonomic analysis, the greatest abundance of yeasts found at the family level corresponds to Saccharomycodaceae, Phiciaceae, and Saccharomycetaceae, and genera such as Hanseniaspora sp., Saccharomyces sp., Candida sp., Meyerozyma sp., Pichia sp., and Wickerhamomyces sp., together with sp. and another group of unidentified yeasts, suggest the presence of new fungal species. This finding demonstrates the diversity in the existing microbiota in the coffee of the department of Cesar. Various investigations have shown the presence of yeasts such as Candida parapsilosis, Saccharomyces cerevisiae, Pichia guilliermondii, and Hanseniaspora opuntiae at the beginning and end of the fermentation process of both washed and unwashed coffee in Brazil [22]. In particular, in this work, 20 yeasts grouped into six genera were found, and the Hanseniaspora sp. was the most abundant yeast species in coffee fermentation in the department of Cesar. Its presence within the consortium of microorganisms could guarantee a better quality in the drink due to its high pectinolytic activity, excellent fermentative capacity, and the production of desirable metabolisms in the final product, such as acids, alcohols, pyridines, aldehydes, and furans, which results in a favorable tasting cup of coffee [59,60,61]. Cruz-O’Byrne et al. [16] also indicated abundant yeasts of the Saccharomycetaceae family such as Kazachstania humilis in coffee fermentation processes in the department of Magdalena. While Peñuela-Martínez et al. [18] identified the yeasts of the Saccharomycedaceae family as the most abundant, with 33 to 52%, on average, in the fermentation process.
The presence of Hanseniaspora, Wickerhamomyces, and Pichia in this study is similar to that reported in the communities of microorganisms found during coffee fermentation in China, Ecuador, and Brazil [62,63,64]. The genus Pichia has been reported as a dominant yeast in coffee fermentation in different countries and it is associated with the production of high-quality coffees with sweet, floral, and fruity characteristics [58,59,65,66]. Because it has pectinolytic activity, it helps in the degradation of the pectin present in the pulp and mucilage of the coffee by producing metabolites that diffuse into the interior of the coffee beans, favoring the formation of the flavor of the final drink. According to Kurtzman [67], Meyerozyma, also identified in this study, is a species of yeast capable of fermenting glucose, sucrose, and raffinose and of maintaining it at a mesophilic temperature (20–45 °C) during fermentation [68]. This genus was also detected in coffee fermentation processes in Australia and Brazil [22,69,70].
Regarding the mycelial fungi found and belonging to the Nectriaceae family (Penicillium sp., Rhizopus spp., Aspergillus fumigatus, Aspergillus flavus, Cladosporium sp., and Fusarium), these fungi have also been reported by other authors in fermentation processes in Brazil, where genera such as Penicillium, Fusarium, Aspergillus, and Cladosporium represented 42.6% of the genera found in the fermentation process [54,71]. It is important to be aware of the presence of species such as Penicillum and Aspergillus because, under some conditions, they are producers of Ochratoxin A-OTA [72].
The diversity analysis of the bacterial and fungal communities according to the coffee variety showed differences in the group of bacteria present in Var. Cenicafé 1, with respect to Var. Colombia and Var. Castillo. There was a relative abundance of bacteria of the genus Secundilactobacillus in the fermentation process, whereas Acetobacter spp. decreased; in particular, a high level of participation of bacteria classified as “Others” (corresponding to the Lactobacillacea family) was observed. Similar results were found in the group of predominant yeasts, with the Saccharomycetaceae family being the most important in the coffee fermentation process of Var. Cenicafé 1, which differs from the other two varieties. In Var. Castillo General and Colombia, the dynamics were more variable between times and the yeasts of the Saccharomycodaceae family predominated.
The differences between the farms were analyzed according to the fermentation time that is also related to the farm altitude, with higher altitudes resulting in longer fermentation times. The results indicated that the abundance of the microorganisms during fermentation is not related to the altitude or fermentation times, but depends on other factors. The farm with the shortest fermentation time at a low altitude and the farm with the longest fermentation time at a high altitude both showed the presence of similar microorganisms including fungi. In this case, to reach a complete fermentation, allowing the mucilage to be degraded completely, seems to be the important point.
The sensory analyses of the 60 coffee seed samples with the SCA methodology showed that “Very Good” coffees were obtained from the three varieties, a total of 92% of the samples obtained scores between 80.1 and 84.9. The average global score was 82.7%. The samples from Cenicafé 1 (36 h) stood out, showing a clearly different pattern of microorganisms, as observed in the beta diversity analysis and metataxonomic analysis and 84.0 SCA. On the other hand, the other two varieties Var. Castillo General ® (F2—36 h) and Var. Colombia (F11—36 h), also show one of the highest scores between 83 and 84.25. The results also show that good practices such as a good quality of harvesting using the Mediverdes® tools and the determination of the coffee’s final fermentation using the Fermaestro® tool, allow us to achieve a very good cup.
In 50% of the farms, it was observed that the fermentation process allowed for an increase in the final sensory quality scores, and only in two farms (F1–F20) were the defects possibly associated with the increase in the population of bacteria of the genus Acetobacter and the interaction with other microorganisms during fermentation. In farm F1, at the fungal level, yeasts of the genus Candida also increased significantly, while in farm F20, unclassified genera were identified. It is possible that the interaction between acetic bacteria and certain fungal groups gives rise to undesirable metabolites for the development of sensory characteristics.
The best scores were obtained on farms F2 and F7 (with 36 h fermentation), which showed complete fermentation. In farm F2, in the metataxonomic observation of bacteria, the high abundance of Leucnostoc and Weisella stands out, and in farm F7, the abundance of the “Others” group and Secundilactobacillus stands out. Although there is also some presence of bacteria of the genus Acetobacter and Bacillus spp., unclassified Enterobacteria and yeasts in F2, the family Saccharomycodaceae-unclassified predominated, while in F7, f_Saccharomycetaceae_unclassified predominated.
On the other hand, in five farms, the fermentation was incomplete, so the quality of the coffee did not change or did not improve; therefore, it can be concluded that the time used for fermentation was not enough to guarantee the removal of the mucilage and improve its sensory quality. If the fermentation had reached its end, in at least 15 farms (75%), the wet fermentation process could have increased the final quality of the coffee. The use of the Fermaestro® tool, which indicates the final fermentation time, could have significantly improved the process on all farms. With these results, it is possible to show the potential of the region in the production of special coffee (scores > 80) SCA, with higher coffee attributes related to acidity, fragrance and aroma, and different coffee sensory descriptors associated with the longest fermentation times.

5. Conclusions

The natural microbial diversity that the Cesar region possesses is involved in the processes of spontaneous fermentation in coffee. Although the microbial diversity in coffee fermentation is specific to the geographical region, the specific conditions of each farm that combine altitude, temperature, and water quality, among others, could determine the proportions and interactions that occur between microbial groups and, therefore, the production of certain metabolites that favor (or do not) the sensory characteristics of coffee in a cup. The metataxonomic analyses and the isolated microorganism suggest that the central microbiome in coffee fermentation from Cesar is composed mainly of Enterobacter, LAB of the genera Leuconostoc and Lactiplantibacillus plantarum, AAB with Gluconobacter and yeasts of the genus Hanseniaspora, which are associated with the region. The differences between the microorganism genera could be influenced by the coffee variety, as has been shown in other investigations; this is because although most of the results correspond to the farms with Var. Colombia and Var. Castillo General, the only farm with Var. Cenicafé 1 showed differences between bacterial and fungal genera that need to be better understood. Finally, the farms with the longest fermentation time (36 to 42 h) showed the highest final total acidity measure in the mucilage and highest coffee cup sensory attributes, as follows: acidity, fragrance and aroma, and special descriptors such as blackberries and fruity flavors.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods13060839/s1, Table S1: Technical characteristics of the farms from the department of Cesar; Table S2: Characterization of the coffee harvest; Table S3: Physicochemical results of pH, Brix degrees, and total acidity of the mucilage; Table S4: Microbial counts during the fermentation process; Table S5: High throughput sequencing of the 16S rRNA gene with the MiSeq platform from Ilumina and richness and diversity indices; Table S6: Identification and frequency of mesophilic aerobic bacteria; Table S7: High-throughput sequencing of ITS with the MiSeq platform from Ilumina and richness and diversity indices; Table S8: Identification and frequency of yeasts isolated; Table S9: Identification and frequency of mycelial fungi into pure culture; Table S10: Total SCA score and sensory descriptors in farms from the department of Cesar.

Author Contributions

C.E.G., L.H.-S. and B.P.-C. contributed to the conceptualization and design of the study; C.E.G. supervised the work; L.H.-S., B.P.-C., R.P.-S. and A.O. were responsible for the methodology; C.E.G., L.H.-S. and B.P.-C. wrote the first draft of the manuscript; C.E.G., L.H.-S., B.P.-C. and R.P.-S. wrote sections of the manuscript; C.E.G., L.H.-S. and L.N.-E. performed the data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Contract No. 2019 02 1486 for the development of scientific and technological activities, signed between the Department of Cesar and the National Federation of Coffee Growers of Colombia for the purpose of the “Experimental Development for the improvement of the Competitiveness of the coffee sector of the Department of Cesar”, code BPIN 2017000100036, financed by the Fund for Science, Technology and Innovation-FCTeI, of the General System of Royalties—SGR and National Coffee Research Center (Cenicafé) project number POS101024.

Institutional Review Board Statement

Sampling was carried out under the guidelines of the Collection Framework Permit 01749 granted by the Environmental Licenses Authority-ANLA to the Universidad de Santander, Colombia.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repositories and accession number(s) can be found at https://www.ncbi.nlm.nih.gov/ (accessed on 30 May 2023) Bioproject: PRJNA977378 and PRJNA977712.

Acknowledgments

We thank all the growers that allowed sampling in the farms. We would also like to thank personnel from the Pueblo Bello station in Cesar, José Enrique Baute Belacazar, Pedro Hernandez, personnel of the Cesar-Guajira Committee of Colombian Coffee Growers Federation (FNC), Helvert Acuña, Liseth Javela, and field promotors for logistical support. Karen Arias Perez for laboratory support. Sensory panel of the Quality Laboratory of the Cesar-Guajira Coffee Growers Committee. We also thank the Multilab laboratory where the analysis of minor elements was carried out, Ruben Medina and Luz Fanny Echeverri for quality and word cloud analysis, Juan F Alzate and Katherine Bedoya from the National Center for Genomic Sequencing, University of Antioquia, and the FASPLAN Company for discussions about the metataxonomic analysis.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sabio, R.P.; Spers, E.E. Does Coffee Origin Matter? An Analysis of Consumer Behavior Based on Regional and National Origin. In Coffee Consumption and Industry Strategies in Brazil; Elsevier: Amsterdam, The Netherlands, 2020; pp. 297–320. ISBN 978-0-12-814721-4. [Google Scholar]
  2. Superintendencia de Industria y Comercio Resolución 4819. 2005. Available online: https://www.sic.gov.co/sites/default/files/files/Denominacion%20de%20Origen/Agro%20-%20Alimenticios/Caf%C3%A9%20de%20Colombia/cafe_de_colombia.pdf (accessed on 8 September 2021).
  3. Corpocesar Generalidades del Departamento de Cesar. Available online: https://www.corpocesar.gov.co/generalidades.html (accessed on 8 September 2021).
  4. Pertuz-Peralta, V.; Perez-Orozco, A.B. Modelo de Cultura Organizacional Innovadora En Caficultores Del Departamento Del Cesar, Colombia. Rev. Investig. Desarro. Innov. 2016, 6, 117–130. [Google Scholar] [CrossRef]
  5. Superintendencia de Industria y Comercio Resolución 2484. 2017. Available online: https://www.sic.gov.co/sites/default/files/files/Denominacion%20de%20Origen/Resoluci%C3%B3n%202484%20Caf%C3%A9%20de%20la%20Sierra%20Nevada.pdf (accessed on 10 September 2021).
  6. Louzada Pereira, L.; Rizzo Moreira, T. (Eds.) Quality Determinants in Coffee Production; Food Engineering Series; Springer International Publishing: Cham, Switzerland, 2021; ISBN 978-3-030-54436-2. [Google Scholar]
  7. Zhang, S.J.; De Bruyn, F.; Pothakos, V.; Contreras, G.F.; Cai, Z.; Moccand, C.; Weckx, S.; De Vuyst, L. Influence of Various Processing Parameters on the Microbial Community Dynamics, Metabolomic Profiles, and Cup Quality During Wet Coffee Processing. Front. Microbiol. 2019, 10, 2621. [Google Scholar] [CrossRef]
  8. De Bruyn, F.; Zhang, S.J.; Pothakos, V.; Torres, J.; Lambot, C.; Moroni, A.V.; Callanan, M.; Sybesma, W.; Weckx, S.; De Vuyst, L. Exploring the Impacts of Postharvest Processing on the Microbiota and Metabolite Profiles during Green Coffee Bean Production. Appl. Environ. Microbiol. 2017, 83, e02398-16. [Google Scholar] [CrossRef]
  9. Ribeiro, L.S.; Evangelista, S.R.; da Cruz Pedrozo Miguel, M.G.; van Mullem, J.; Silva, C.F.; Schwan, R.F. Microbiological and Chemical-Sensory Characteristics of Three Coffee Varieties Processed by Wet Fermentation. Ann. Microbiol. 2018, 68, 705–716. [Google Scholar] [CrossRef]
  10. Holguín-Sterling, L.; Pedraza-Claros, B.; Pérez-Salinas, R.; Ortiz, A.; Navarro-Escalante, L.; Góngora, C.E. Physical–Chemical and Metataxonomic Characterization of the Microbial Communities Present during the Fermentation of Three Varieties of Coffee from Colombia and Their Sensory Qualities. Agriculture 2023, 13, 1980. [Google Scholar] [CrossRef]
  11. Martinez, S.J.; Simão, J.B.P.; Pylro, V.S.; Schwan, R.F. The Altitude of Coffee Cultivation Causes Shifts in the Microbial Community Assembly and Biochemical Compounds in Natural Induced Anaerobic Fermentations. Front. Microbiol. 2021, 12, 671395. [Google Scholar] [CrossRef]
  12. Puerta Quintero, G.I.; Marín Mejía, J.; Osorio Betancur, G.A. Microbiología de La Fermentación Del Mucílago de Café Según Su Madurez y Selección. Rev. Cenicafé 2012, 63, 58–78. [Google Scholar]
  13. Peñuela-Martínez, A.E.; Moreno-Riascos, S.; Medina-Rivera, R. Influence of Temperature-Controlled Fermentation on the Quality of Mild Coffee (Coffea arabica L.) Cultivated at Different Elevations. Agriculture 2023, 13, 1132. [Google Scholar] [CrossRef]
  14. Elhalis, H.; Cox, J.; Zhao, J. Coffee Fermentation: Expedition from Traditional to Controlled Process and Perspectives for Industrialization. Appl. Food Res. 2023, 3, 100253. [Google Scholar] [CrossRef]
  15. Sunarharum, W.B.; Williams, D.J.; Smyth, H.E. Complexity of Coffee Flavor: A Compositional and Sensory Perspective. Food Res. Int. 2014, 62, 315–325. [Google Scholar] [CrossRef]
  16. Cruz-O’Byrne, R.; Piraneque-Gambasica, N.; Aguirre-Forero, S. Microbial Diversity Associated with Spontaneous Coffee Bean Fermentation Process and Specialty Coffee Production in Northern Colombia. Int. J. Food Microbiol. 2021, 354, 109282. [Google Scholar] [CrossRef]
  17. De Oliveira Junqueira, A.C.; de Melo Pereira, G.V.; Coral Medina, J.D.; Alvear, M.C.R.; Rosero, R.; de Carvalho Neto, D.P.; Enríquez, H.G.; Soccol, C.R. First Description of Bacterial and Fungal Communities in Colombian Coffee Beans Fermentation Analysed Using Illumina-Based Amplicon Sequencing. Sci. Rep. 2019, 9, 8794. [Google Scholar] [CrossRef] [PubMed]
  18. Peñuela-Martínez, A.E.; Velasquez-Emiliani, A.V.; Angel, C.A. Microbial Diversity Using a Metataxonomic Approach, Associated with Coffee Fermentation Processes in the Department of Quindío, Colombia. Fermentation 2023, 9, 343. [Google Scholar] [CrossRef]
  19. Zhang, S.J.; De Bruyn, F.; Pothakos, V.; Torres, J.; Falconi, C.; Moccand, C.; Weckx, S.; De Vuyst, L. Following Coffee Production from Cherries to Cup: Microbiological and Metabolomic Analysis of Wet Processing of Coffea arabica. Appl. Environ. Microbiol. 2019, 85, e02635-18. [Google Scholar] [CrossRef] [PubMed]
  20. De Melo Pereira, G.V.; da Silva Vale, A.; de Carvalho Neto, D.P.; Muynarsk, E.S.; Soccol, V.T.; Soccol, C.R. Lactic Acid Bacteria: What Coffee Industry Should Know? Curr. Opin. Food Sci. 2020, 31, 1–8. [Google Scholar] [CrossRef]
  21. Elhalis, H.; Cox, J.; Frank, D.; Zhao, J. The Crucial Role of Yeasts in the Wet Fermentation of Coffee Beans and Quality. Int. J. Food Microbiol. 2020, 333, 108796. [Google Scholar] [CrossRef] [PubMed]
  22. Evangelista, S.R.; da Cruz Pedroso Miguel, M.G.; Silva, C.F.; Pinheiro, A.C.M.; Schwan, R.F. Microbiological Diversity Associated with the Spontaneous Wet Method of Coffee Fermentation. Int. J. Food Microbiol. 2015, 210, 102–112. [Google Scholar] [CrossRef]
  23. Martins, P.M.M.; Ribeiro, L.S.; da Cruz Pedroso Miguel, M.G.; Evangelista, S.R.; Schwan, R.F. Production of Coffee (Coffea arabica) Inoculated with Yeasts: Impact on Quality. J. Sci. Food Agric. 2019, 99, 5638–5645. [Google Scholar] [CrossRef] [PubMed]
  24. Cordoba, N.; Fernandez-Alduenda, M.; Moreno, F.L.; Ruiz, Y. Coffee Extraction: A Review of Parameters and Their Influence on the Physicochemical Characteristics and Flavour of Coffee Brews. Trends Food Sci. Technol. 2020, 96, 45–60. [Google Scholar] [CrossRef]
  25. Lee, L.W.; Cheong, M.W.; Curran, P.; Yu, B.; Liu, S.Q. Coffee Fermentation and Flavor—An Intricate and Delicate Relationship. Food Chem. 2015, 185, 182–191. [Google Scholar] [CrossRef]
  26. Peñuela-Martínez, A.E.; Guerrero, Á.; Sanz-Uribe, J.R. Cromacafé® Herramienta Para Identificar Los Estados de Madurez de Las Variedades de Café de Fruto Rojo. Av. Técnicos Cenicafé 2022, 535, 1–8. [Google Scholar] [CrossRef]
  27. Guerrero, Á.; Sanz-Uribe, J.R.; Peñuela-Martínez, A.E.; Ramírez, C.A. Mediverdes®: Un Método Para Medir La Calidad de La Recolección Del Café En El Campo. Av. Técnicos Cenicafé 2022, 536, 1–8. [Google Scholar] [CrossRef]
  28. Peñuela, A.E.; Pabón, J.; Sanz-Uribe, J.R. Método Fermaestro: Para Determinar La Finalización de La Fermentación Del Mucílago de Café. Av. Técnicos Cenicafé 2013, 431, 1–8. [Google Scholar] [CrossRef]
  29. Oliveros-Tascón, C.E.; Ramírez, C.A.; Tibaduiza-Vianchá, C.A.; Sanz-Uribe, J.R. Construcción de Secadores Solares Tipo Túnel Con Nuevos Materiales. Av. Técnicos Cenicafé 2017, 482, 1–8. [Google Scholar]
  30. Tyl, C.; Sadler, G.D. pH and Titratable Acidity. In Food Analysis; Nielsen, S.S., Ed.; Food Science Text Series; Springer International Publishing: Cham, Switzerland, 2017; pp. 389–406. ISBN 978-3-319-45774-1. [Google Scholar]
  31. Mizrahi-Man, O.; Davenport, E.R.; Gilad, Y. Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs. PLoS ONE 2013, 8, e53608. [Google Scholar] [CrossRef]
  32. White, T.J.; Bruns, T.D.; Lee, S.B.; Taylor, J.W. Amplification and Direct Sequencing of Fungal Ribosomal RNA Genes for Phylogenetics. In PCR—Protocols and Applications—A Laboratory Manual; Academic Press: Cambridge, MA, USA, 1990; pp. 315–322. [Google Scholar]
  33. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  34. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing Mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [PubMed]
  35. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef] [PubMed]
  36. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef] [PubMed]
  37. Cole, J.R.; Wang, Q.; Fish, J.A.; Chai, B.; McGarrell, D.M.; Sun, Y.; Brown, C.T.; Porras-Alfaro, A.; Kuske, C.R.; Tiedje, J.M. Ribosomal Database Project: Data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014, 42, D633–D642. [Google Scholar] [CrossRef]
  38. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26, 32–46. [Google Scholar]
  39. Mcmurdie, P.J.; Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed]
  40. Instituto Colombiano de Normas Técnicas y Certificación NTC 4092. 2016. Available online: https://tienda.icontec.org/gp-microbiologia-de-alimentos-y-productos-para-alimentacion-animal-requisitos-generales-y-directrices-para-analisis-microbiologicos-ntc4092-2016.html (accessed on 20 July 2021).
  41. Instituto Colombiano de Normas Técnicas y Certificación NTC 5181. 2003. Available online: https://tienda.icontec.org/gp-buenas-practicas-de-manufactura-para-la-industria-del-cafe-ntc5181-2003.html (accessed on 20 July 2023).
  42. Gerard, L. Caracterización de Bacterias del ácido Acético Destinadas a la Producción de Vinagres de Frutas 2015. Available online: https://riunet.upv.es/bitstream/handle/10251/59401/GERARD%20-%20Caracterizaci%C3%B3n%20de%20bacterias%20del%20%C3%A1cido%20ac%C3%A9tico%20destinadas%20ª%20la%20producci%C3%B3n%20de%20vinagres%20de....pdf?sequence=1 (accessed on 13 March 2022).
  43. Seearunruangchai, A.; Tanasupawat, S.; Keeratipibul, S.; Thawai, C.; Itoh, T.; Yamada, Y. Identification of acetic acid bacteria isolated from fruits collected in Thailand. J. Gen. Appl. Microbiol. 2004, 50, 47–53. [Google Scholar] [CrossRef]
  44. Larone, D.H. Medically Important Fungi: A Guide to Identification, 6th ed.; Elsevier: Washington, DC, USA, 2018. [Google Scholar]
  45. Barnett, H. Illustrated Genera of Imperfect Fungi; The American Phytopathological Society: St. Paul, MN, USA, 1998. [Google Scholar]
  46. Ministerio de la Protección Social Ministerio de Ambiente, Vivienda y Desarrollo Territorial Resolución 2115. 2007. Available online: https://minvivienda.gov.co/sites/default/files/normativa/2115%20-%202007.pdf (accessed on 25 September 2021).
  47. Instituto Colombiano de Normas Técnicas y Certificación NTC 4939. 2001. Available online: https://tienda.icontec.org/gp-calidad-del-agua-enumeracion-de-coliformes-y-escherichia-coli-tecnica-con-tubos-de-fermentacion-y-tecnica-de-sustrato-enzimatico-ntc4939-2001.html (accessed on 15 June 2021).
  48. Instituto Colombiano de Normas Técnicas y Certificación NTC 2324. 2021. Available online: https://tienda.icontec.org/gp-cafe-verde-examen-olfativo-y-visual-y-determinacion-de-materia-extrana-y-defectos-ntc2324-2021.html (accessed on 26 July 2021).
  49. Protocols & Best Practices. Specialty Coffee Association Cupping Protocols; Protocols & Best Practices: Irvine, CA, USA, 2015. [Google Scholar]
  50. Sanchez Moreno, I. Dihidroxiacetona Quinasa de Citrobacter Freundii CECT 4626. Una Aproximación Integral a Un Nuevo Biocatalizador Para La Formación de Enlaces C-C 2009. Available online: https://digital.csic.es/bitstream/10261/17036/1/Tesis%20Israel%20Sh%C3%A1nchez%20Moreno.pdf (accessed on 23 July 2022).
  51. Horinouchi, S.; Ueda, K.; Nakayama, J.; Ikeda, T. Cell-to-Cell Communications among Microorganisms. In Comprehensive Natural Products II; Elsevier: Amsterdam, The Netherlands, 2010; pp. 283–337. ISBN 978-0-08-045382-8. [Google Scholar]
  52. Pereira, T.S.; Batista, N.N.; Santos Pimenta, L.P.; Martinez, S.J.; Ribeiro, L.S.; Oliveira Naves, J.A.; Schwan, R.F. Self-Induced Anaerobiosis Coffee Fermentation: Impact on Microbial Communities, Chemical Composition and Sensory Quality of Coffee. Food Microbiol. 2022, 103, 103962. [Google Scholar] [CrossRef]
  53. Abedi, E.; Hashemi, S.M.B. Lactic Acid Production—Producing Microorganisms and Substrates Sources-State of Art. Heliyon 2020, 6, e04974. [Google Scholar] [CrossRef] [PubMed]
  54. Vilela, D.M.; de, M. Pereira, G.V.; Silva, C.F.; Batista, L.R.; Schwan, R.F. Molecular Ecology and Polyphasic Characterization of the Microbiota Associated with Semi-Dry Processed Coffee (Coffea arabica L.). Food Microbiol. 2010, 27, 1128–1135. [Google Scholar] [CrossRef]
  55. Hameed, A.; Hussain, S.A.; Ijaz, M.U.; Ullah, S.; Pasha, I.; Suleria, H.A.R. Farm to Consumer: Factors Affecting the Organoleptic Characteristics of Coffee. II: Postharvest Processing Factors. Compr. Rev. Food Sci. Food Saf. 2018, 17, 1184–1237. [Google Scholar] [CrossRef]
  56. Saichana, N.; Matsushita, K.; Adachi, O.; Frébort, I.; Frebortova, J. Acetic Acid Bacteria: A Group of Bacteria with Versatile Biotechnological Applications. Biotechnol. Adv. 2015, 33, 1260–1271. [Google Scholar] [CrossRef]
  57. Avallone, S.; Guyot, B.; Brillouet, J.-M.; Olguin, E.; Guiraud, J.-P. Microbiological and Biochemical Study of Coffee Fermentation. Curr. Microbiol. 2001, 42, 252–256. [Google Scholar] [CrossRef] [PubMed]
  58. De Melo Pereira, G.; Soccol, V.T.; Brar, S.K.; Neto, E.; Soccol, C.R. Microbial Ecology and Starter Culture Technology in Coffee Processing. Crit. Rev. Food Sci. Nutr. 2017, 57, 2775–2788. [Google Scholar] [CrossRef]
  59. Silva, C.F.; Vilela, D.M.; de Souza Cordeiro, C.; Duarte, W.F.; Dias, D.R.; Schwan, R.F. Evaluation of a Potential Starter Culture for Enhance Quality of Coffee Fermentation. World J. Microbiol. Biotechnol. 2013, 29, 235–247. [Google Scholar] [CrossRef] [PubMed]
  60. Evangelista, S.R.; da Cruz Pedrozo Miguel, M.G.; de Souza Cordeiro, C.; Silva, C.F.; Marques Pinheiro, A.C.; Schwan, R.F. Inoculation of Starter Cultures in a Semi-Dry Coffee (Coffea arabica) Fermentation Process. Food Microbiol. 2014, 44, 87–95. [Google Scholar] [CrossRef]
  61. Jackels, S.C.; Jackels, C.F. Characterization of the Coffee Mucilage Fermentation Process Using Chemical Indicators: A Field Study in Nicaragua. J. Food Sci. 2005, 70, C321–C325. [Google Scholar] [CrossRef]
  62. Feng, X.; Dong, H.; Yang, P.; Yang, R.; Lu, J.; Lv, J.; Sheng, J. Culture-Dependent and -Independent Methods to Investigate the Predominant Microorganisms Associated with Wet Processed Coffee. Curr. Microbiol. 2016, 73, 190–195. [Google Scholar] [CrossRef]
  63. Martins, P.M.M.; Batista, N.N.; da Cruz Pedroso Miguel, M.G.; Simão, J.B.P.; Soares, J.R.; Schwan, R.F. Coffee Growing Altitude Influences the Microbiota, Chemical Compounds and the Quality of Fermented Coffees. Food Res. Int. 2020, 129, 108872. [Google Scholar] [CrossRef]
  64. Pothakos, V.; De Vuyst, L.; Zhang, S.J.; De Bruyn, F.; Verce, M.; Torres, J.; Callanan, M.; Moccand, C.; Weckx, S. Temporal Shotgun Metagenomics of an Ecuadorian Coffee Fermentation Process Highlights the Predominance of Lactic Acid Bacteria. Curr. Res. Biotechnol. 2020, 2, 1–15. [Google Scholar] [CrossRef]
  65. De Melo Pereira, G.V.; Soccol, V.T.; Pandey, A.; Medeiros, A.B.P.; Andrade Lara, J.M.R.; Gollo, A.L.; Soccol, C.R. Isolation, Selection and Evaluation of Yeasts for Use in Fermentation of Coffee Beans by the Wet Process. Int. J. Food Microbiol. 2014, 188, 60–66. [Google Scholar] [CrossRef] [PubMed]
  66. Evangelista, S.R.; Silva, C.F.; da Cruz Pedroso Miguel, M.G.; Cordeiro, C.d.S.; Pinheiro, A.C.M.; Duarte, W.F.; Schwan, R.F. Improvement of Coffee Beverage Quality by Using Selected Yeasts Strains during the Fermentation in Dry Process. Food Res. Int. 2014, 61, 183–195. [Google Scholar] [CrossRef]
  67. Kurtzman, C. Phylogenetic circumscription of Saccharomyces, Kluyveromyces and other members of the Saccharomycetaceae, and the proposal of the new genera Lachancea, Nakaseomyces, Naumovia, Vanderwaltozyma and Zygotorulaspora. FEMS Yeast Res. 2003, 4, 233–245. [Google Scholar] [CrossRef]
  68. Rodrigues, T.B.; Silva, A.E.T. (Eds.) Molecular Diversity of Environmental Prokaryotes; CRC Press: Boca Raton, FL, USA, 2016; ISBN 978-1-4822-3323-0. [Google Scholar]
  69. De Carvalho Neto, D.; de Melo Pereira, G.; Tanobe, V.; Thomaz Soccol, V.; da Silva, B.J.; Rodrigues, C.; Soccol, C. Yeast Diversity and Physicochemical Characteristics Associated with Coffee Bean Fermentation from the Brazilian Cerrado Mineiro Region. Fermentation 2017, 3, 11. [Google Scholar] [CrossRef]
  70. Elhalis, H.; Cox, J.; Zhao, J. Ecological Diversity, Evolution and Metabolism of Microbial Communities in the Wet Fermentation of Australian Coffee Beans. Int. J. Food Microbiol. 2020, 321, 108544. [Google Scholar] [CrossRef] [PubMed]
  71. Silva, C.F.; Schwan, R.F.; Sousa Dias, Ë.; Wheals, A.E. Microbial Diversity during Maturation and Natural Processing of Coffee Cherries of Coffea arabica in Brazil. Int. J. Food Microbiol. 2000, 60, 251–260. [Google Scholar] [CrossRef] [PubMed]
  72. Bessaire, T.; Perrin, I.; Tarres, A.; Bebius, A.; Reding, F.; Theurillat, V. Mycotoxins in Green Coffee: Occurrence and Risk Assessment. Food Control 2019, 96, 59–67. [Google Scholar] [CrossRef]
Figure 1. Temperature in the fermentation tanks during the fermentation process in the 20 farms identified: F1 to F20 in the department of Cesar.
Figure 1. Temperature in the fermentation tanks during the fermentation process in the 20 farms identified: F1 to F20 in the department of Cesar.
Foods 13 00839 g001
Figure 2. Physicochemical results of pH and total acidity of the mucilage during coffee fermentation in 20 farms. The data were grouped into three groups according to the farm’s final fermentation time. T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Figure 2. Physicochemical results of pH and total acidity of the mucilage during coffee fermentation in 20 farms. The data were grouped into three groups according to the farm’s final fermentation time. T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Foods 13 00839 g002
Figure 3. Physicochemical results of sucrose and reducing sugar of the mucilage during coffee fermentation in the 20 farms. The data were grouped into three groups according to the farm’s final fermentation time. T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Figure 3. Physicochemical results of sucrose and reducing sugar of the mucilage during coffee fermentation in the 20 farms. The data were grouped into three groups according to the farm’s final fermentation time. T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Foods 13 00839 g003
Figure 4. Nonmetric multidimensional scaling of the bacterial genus in mucilage from coffee samples during the fermentation process.
Figure 4. Nonmetric multidimensional scaling of the bacterial genus in mucilage from coffee samples during the fermentation process.
Foods 13 00839 g004
Figure 5. The relative abundance of sequences (%) corresponding to the bacterial taxonomic assignment to family category in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation.
Figure 5. The relative abundance of sequences (%) corresponding to the bacterial taxonomic assignment to family category in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation.
Foods 13 00839 g005
Figure 6. The relative abundance of sequences (%) corresponding to the bacterial taxonomic assignment to the genus category in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Figure 6. The relative abundance of sequences (%) corresponding to the bacterial taxonomic assignment to the genus category in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Foods 13 00839 g006
Figure 7. Mesophilic aerobic bacteria isolated during the coffee fermentation process at three different times in the 20 studied farms.
Figure 7. Mesophilic aerobic bacteria isolated during the coffee fermentation process at three different times in the 20 studied farms.
Foods 13 00839 g007
Figure 8. Nonmetric multidimensional scaling of fungi genera in mucilage from coffee samples during the fermentation process.
Figure 8. Nonmetric multidimensional scaling of fungi genera in mucilage from coffee samples during the fermentation process.
Foods 13 00839 g008
Figure 9. The relative abundance of sequences (%) corresponding to the fungi taxonomic assignment to the family level in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Figure 9. The relative abundance of sequences (%) corresponding to the fungi taxonomic assignment to the family level in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Foods 13 00839 g009
Figure 10. Relative abundance of sequences (%) corresponding to the fungi taxonomic assignment to the genus level in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Figure 10. Relative abundance of sequences (%) corresponding to the fungi taxonomic assignment to the genus level in coffee fermentation processes in 20 farms (Fs) in the department of Cesar. T01: Start of fermentation. T02: Midpoint during the fermentation and T03: end of fermentation. Detected OTUs below 0.5% are indicated as “Others”.
Foods 13 00839 g010
Figure 11. Yeast isolates during the coffee fermentation process at three different times in the 20 studied farms.
Figure 11. Yeast isolates during the coffee fermentation process at three different times in the 20 studied farms.
Foods 13 00839 g011
Figure 12. Word Cloud that shows the predominance of sensory descriptors in the coffee cup of farms with fermentation of 36–42 h.
Figure 12. Word Cloud that shows the predominance of sensory descriptors in the coffee cup of farms with fermentation of 36–42 h.
Foods 13 00839 g012
Table 1. Microbial counts during the fermentation process in 20 farms in the department of Cesar. Data are grouped into three groups according to the farm’s final fermentation time.
Table 1. Microbial counts during the fermentation process in 20 farms in the department of Cesar. Data are grouped into three groups according to the farm’s final fermentation time.
Fermentation Time (h)/Farms (Fs)Sample TimeMesophilesLactic Acid BacteriaColiforms YeastsFilamentous Fungi
(Log10 CFUs/mL)
X SDX SDX SDX SDX SD
10–16 h
(F14 − F12 − F13 − F15 − F16 − F17)
T15.48±0.116.57±0.654.62±0.455.75±0.122.40±1.86
T25.44±0.196.52±0.653.56±1.765.68±0.142.78±1.42
T35.35±0.406.44±0.603.80±1.875.73±0.342.30±1.79
18–20 h
(F3 − F4 − F5 − F8 − F9 − F10 − F18 − F19 − F20)
T15.50±0.3026.61±0.554.22±0.855.68±0.263.42±0.34
T25.49±0.4366.65±0.584.13±0.795.55±0.363.58±0.21
T35.62±0.4136.27±1.283.33±2.045.65±0.333.04±1.20
36–42 h
(F1 − F2 − F6 − F7 − F11)
T15.57±0.2025.80±2.353.99±0.855.75±0.212.93±1.70
T25.67±0.4825.18±2.423.53±0.645.85±0.262.94±1.70
T35.76±0.5394.98±2.453.66±0.935.82±0.402.83±1.69
T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Table 2. Identification and frequency of lactic acid bacteria (LAB) isolated into pure culture in coffee mucilage during the fermentation process at three different times in the 20 studied farms.
Table 2. Identification and frequency of lactic acid bacteria (LAB) isolated into pure culture in coffee mucilage during the fermentation process at three different times in the 20 studied farms.
MicroorganismTimeFarmsFrequency by Farms
F1F2F3F4F5F6F7F8F9F10F11F12F13F14F15F16F17F18F19F20#%
Fructilactobacillus fructivoransT1 15.0
T2 x
T3
Lactiplantibacillus plantarumT1 xxxxxxxxxxxxxx x 1785.0
T2 xxxxxxx xxxxxxxx
T3 xx xx xxx x xx
Lactiplantibacillus pentosusT1 x x x630.0
T2 x x
T3x x x
Lactobacillus delbrueckiiT1x 210.0
T2 x
T3
Lactococcus raffinolactisT1 x 15.0
T2
T3
Leuconostoc mesenteroidesT1 xx xx x x 945.0
T2 x xxx
T3 xxxx x x
Leuconostoc citreumT1 xx 315.0
T2 x x
T3
Levilactobacillus brevisT1 x x 840.0
T2 xx x
T3 x x x x
T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Table 3. Identification and frequency of AAB isolated into pure culture in coffee mucilage during the fermentation process at three different times in the 20 studied farms.
Table 3. Identification and frequency of AAB isolated into pure culture in coffee mucilage during the fermentation process at three different times in the 20 studied farms.
MicroorganismTimeFarmsFrequency by Farms
F1F2F3F4F5F6F7F8F9F10F11F12F13F14F15F16F17F18F19F20#%
Acetobacter spp.T1 x x x 525.0
T2 x
T3 x x xx
Acidomonas spp.T1 15.0
T2 x
T3
Gluconobacter spp.T1 x x 525.0
T2 x x x
T3
Sacharibacter spp.T1 x 210.0
T2
T3 x
T1 (zero hours of fermentation). T2 (middle of the fermentation time). T3 (end of the fermentation time).
Table 4. Total SCA score in the coffee cup in 20 farms in the department of Cesar. Data were grouped into three groups according to the farm’s final fermentation time.
Table 4. Total SCA score in the coffee cup in 20 farms in the department of Cesar. Data were grouped into three groups according to the farm’s final fermentation time.
Fermentation Time (h)/Farms (Fs) Sample Time
T1T2T3
X SDX SDX SD
10–16 h
(F12 − F13 − F14 − F15 − F16 − F17)
82.7±0.8682.7±0.4382.5±0.86
18–20 h
(F3 − F4 − F5 − F8 − F9 − F10 − F18 − F19 − F20)
82.3±1.4382.8±6.9682.5±0.71
36–42 h
(F1 − F2 − F6 − F7 − F11)
82.5±1.6082.7±8.9083.7±0.54
SCA score: 90–100 outstanding, 85.0–89.99 excellent, 80–84.99 very Good, <80 commercial coffee.
Table 5. Score of some of the sensory attributes identified in the coffee cup in 20 farms in the department of Cesar. Data were grouped into three groups according to the farm’s final fermentation time.
Table 5. Score of some of the sensory attributes identified in the coffee cup in 20 farms in the department of Cesar. Data were grouped into three groups according to the farm’s final fermentation time.
Fermentation Time (h)/Farms (Fs)Sample TimeFragance/AromaFlavor Acidity
X SDX SDX SD
10–16 h
(F12 − F13 − F14 − F15 − F16 − F17)
T17.55±0.157.61±0.147.55±0.10
T27.65±0.097.54±0.137.54±0.07
T37.50±0.117.57±0.137.52±0.14
18–20 h
(F3 − F4 − F5 − F8 − F9 − F10 − F18 − F19 − F20)
T17.54±0.227.53±0.247.52±0.25
T27.53±0.097.56±0.177.56±0.13
T37.55±0.097.54±0.127.50±0.11
36–42 h
(F1 − F2 − F6 − F7 − F11)
T17.65±0.277.58±0.177.55±0.19
T27.69±0.167.54±0.197.60±0.19
T37.69±0.137.85±0.127.73±0.21
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Góngora, C.E.; Holguín-Sterling, L.; Pedraza-Claros, B.; Pérez-Salinas, R.; Ortiz, A.; Navarro-Escalante, L. Metataxonomic Identification of Microorganisms during the Coffee Fermentation Process in Colombian Farms (Cesar Department). Foods 2024, 13, 839. https://doi.org/10.3390/foods13060839

AMA Style

Góngora CE, Holguín-Sterling L, Pedraza-Claros B, Pérez-Salinas R, Ortiz A, Navarro-Escalante L. Metataxonomic Identification of Microorganisms during the Coffee Fermentation Process in Colombian Farms (Cesar Department). Foods. 2024; 13(6):839. https://doi.org/10.3390/foods13060839

Chicago/Turabian Style

Góngora, Carmenza E., Laura Holguín-Sterling, Bertilda Pedraza-Claros, Rosangela Pérez-Salinas, Aristofeles Ortiz, and Lucio Navarro-Escalante. 2024. "Metataxonomic Identification of Microorganisms during the Coffee Fermentation Process in Colombian Farms (Cesar Department)" Foods 13, no. 6: 839. https://doi.org/10.3390/foods13060839

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop