Next Article in Journal
Smart Project Management: Interactive Platform Using Natural Language Processing Technology
Next Article in Special Issue
Effects of Different LED Light Recipes and NPK Fertilizers on Basil Cultivation for Automated and Integrated Horticulture Methods
Previous Article in Journal
Efficient Sulfate Adsorption on Modified Adsorbents Prepared from Zea mays Stems
Previous Article in Special Issue
Optimization of the Acetification Stage in the Production of Wine Vinegar by Use of Two Serial Bioreactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kombucha Tea as a Reservoir of Cellulose Producing Bacteria: Assessing Diversity among Komagataeibacter Isolates

Department of Life Sciences, University of Modena and Reggio Emilia, 42122 Reggio Emilia, Italy
*
Author to whom correspondence should be addressed.
Submission received: 31 December 2020 / Revised: 3 February 2021 / Accepted: 6 February 2021 / Published: 10 February 2021

Abstract

:
Bacterial cellulose (BC) is receiving a great deal of attention due to its unique properties such as high purity, water retention capacity, high mechanical strength, and biocompatibility. However, the production of BC has been limited because of the associated high costs and low productivity. In light of this, the isolation of new BC producing bacteria and the selection of highly productive strains has become a prominent issue. Kombucha tea is a fermented beverage in which the bacteria fraction of the microbial community is composed mostly of strains belonging to the genus Komagataeibacter. In this study, Kombucha tea production trials were performed starting from a previous batch, and bacterial isolation was conducted along cultivation time. From the whole microbial pool, 46 isolates were tested for their ability to produce BC. The obtained BC yield ranged from 0.59 g/L, for the isolate K2G36, to 23 g/L for K2G30—which used as the reference strain. The genetic intraspecific diversity of the 46 isolates was investigated using two repetitive-sequence-based PCR typing methods: the enterobacterial repetitive intergenic consensus (ERIC) elements and the (GTG)5 sequences, respectively. The results obtained using the two different approaches revealed the suitability of the fingerprint techniques, showing a discrimination power, calculated as the D index, of 0.94 for (GTG)5 rep-PCR and 0.95 for ERIC rep-PCR. In order to improve the sensitivity of the applied method, a combined model for the two genotyping experiments was performed, allowing for the ability to discriminate among strains.

1. Introduction

The attraction of bacterial cellulose (BC) as a natural biopolymer, produced by microorganisms, arises from its main functional properties and subsequent applications, especially in the biomedical, food, and engineering fields [1,2,3]. Bacteria able to synthesize cellulose can be considered ubiquitous since they have been found to inhabit different ecosystems. They can be isolated from environmental sources (such as soil and plants), from insects, humans, and from food sources [4].
Within the family Acetobacteraceae, the genus Komagataeibacter includes species, such as K. xylinus, K. europaeus and K. hansenii, that have been described as cellulose producing organisms [5]. The main studied strains belong to the K. xylinus species, which is recognized as the model organism for studying the mechanism of BC synthesis [6,7].
One of the most important reservoirs in which it is possible to detect strains belonging to the Komagataeibacter genus is kombucha tea, a fermented beverage traditionally produced in Asia, probably originating from the northeast of China [8]. Kombucha tea is characterized by a microbial community in which different microbial groups, mainly yeasts and bacteria, live in a symbiotic lifestyle [9]. Among bacteria groups, acetic acid bacteria are the main functional organisms, although lactic acid bacteria are also found; moreover, among yeasts, strains of the genera Saccharomyces and Zygosaccharomyces have been described [10]. Kombucha tea production consists of two fermentation reactions operated by yeasts and acetic acid bacteria [11]. The symbiotic relation between yeast and bacteria is established during the fermentation process, in which the yeasts are capable of hydrolysing sucrose, the main substrate for kombucha tea production, into glucose and fructose, and then ethanol. Acetic acid bacteria use glucose, fructose and ethanol to produce gluconic acid, glucuronic acid, acetic acid, and BC [12,13].
Kombucha tea can be considered as a dynamic environment, rapidly changing its composition, acting as an ecophysiological reservoir of very specialized organisms. Considering the occurrence and function of acetic acid bacteria within kombucha, previous studies have described K. xylinus as the main species.
The ability of K. xylinus strains to produce BC is highly variable, making the selection of the best ones a crucial step in order to obtain candidate strains for large scale BC production. In fact, the availability of selected organisms in producing BC is a requirement for exploiting them as industrial machinery. The production of BC is mainly influenced by two factors. The first one is represented by the carbon sources used in the medium [14,15]; pure carbon sources could increase the production cost, making the process unsustainable. The second issue is related to the ability of the organism to grow and produce BC in the selected conditions [16]. Some studies have described the abilities of strains belonging to the Komagataeibacter genus for producing BC at different yields [5].
In this work, acetic acid bacteria were collected from liquid and pellicle fraction of kombucha tea samples. In total, 46 isolates were screened for their ability to produce BC and for the main phenotypic characteristics. Moreover, they were typed by applying a fingerprint analysis based on (GTG)5 repetitive elements and Enterobacterial Repetitive Intragenic Consensus (ERIC) elements in order to assay the intraspecific genetic diversity.

2. Materials and Methods

2.1. Kombucha Preparation, Sampling and Physico-Chemical Determinations

Two native local kombucha precultures (P1 and P2) were provided by the National Institute of Applied Sciences and Technology (INSAT, Tunisia), previously cultivated on two different substrates: black tea (P1) and green tea (P2). Each preculture was composed of kombucha cellulose pellicle with approximately 200 mL of starter culture. Duplicates of kombucha, referred to green tea kombucha (GTK) and black tea kombucha (BTK), were prepared and inoculated according to Jayabalan et al. [17]. Briefly, sucrose 10% (w/v) was added to demineralized water and allowed to boil. Afterwards, 1.2% (w/v) black or green tea (Les Jardins du thé, Office Tunisien du Commerce, Tunisia) was added and allowed to infuse for about 5 min. After filtration of tea leaves through a sterile sieve, 200 mL of tea was poured into 500 mL glass jars that had been previously sterilized at +121 °C for 20 min. Once the sucrose–tea solution had cooled at room temperature, 10% (v/v) of the fermented tea (of the previous fermentation brew from kombucha with the same origin, in this case black or green tea as substrate, corresponding to the aforementioned starter culture) was added. Then, 3% (w/v) of representative pellicle fragments cut from the preculture were placed in the culture light side up. The jars were covered with a sterile gauze and fixed with an elastic band. Fermentations were carried out in the dark for 12 days. All experiments were carried out aseptically.
Cell morphology and direct observation of samples were carried out using optical microscopy at 100× of magnification, using C. Zeiss microscope apparatus (Axiolab).
Titratable acidity was determined by neutralizing samples at pH 7.2 with 0.1 N of NaOH (it was assumed that all media acidity was due to acetic acid) and expressed as % (wt/wt); pH was measured by CRISON, MicropH, 2002 pHmeter. Ethanol, expressed as % (v/v), was checked by densitometry measure using a hydrostatic balance after distillation.

2.2. Plating and Isolation

Plating was performed immediately after sampling. Therefore, 10 mL of 0.1% (w/v) peptone water was added to 1 g of black or green kombucha pellicle in a sterile stomacher bag that was vigorously shaken for 5 min in a laboratory blender STOMACHER® 400 (Seward, England) to obtain a uniform homogenate. Samples (1 mL each) of the homogenate and liquid phase were 10-fold serially diluted in 0.1% (w/v) peptone water, from which aliquots (0.1 mL) were plated on GYC medium (50 g/L glucose, 10 g/L yeast extract, 15 g/L CaCO3, 9 g/L bacteriological agar) and ACB agar medium (30 g/L yeast extract, 1 mL of a 22 g/L bromocresol green solution, 20 g/L bacteriological agar). The media were previously autoclaved at 121 °C for 15 min and after 20 mL/L of filtered ethanol 95% (v/v) was added aseptically to the ACB medium. A total of 0.2 mL of cycloheximide hydroalcoholic solution at 25% m/v was added directly to petri dishes.
Plates were incubated at +28 °C for 72–120 h. Colonies were picked up from a suitable dilution of each sample on GYC and ACB agar media, grown for 3–5 days at +28 °C and purified through subculturing and plating. Pure kombucha isolates were named and strains were cultivated on the corresponding isolation medium for 3–5 days at +28 °C [18]. Long-term preservation methods of the strains were performed according to the standard procedures of the Microbial Resource Research Infrastructure—Italian Joint Research Unit (MIRRI-IT) [19]. Specifically, a seed lot of the strains was stored at −80 °C, adding glycerol 50% (v/v) to the liquid media. The strains with significant cellulose production were deposited in the University of Modena and Reggio Emilia (UNIMORE) Microbial Culture Collection (UMCC, www.umcc.unimore.it, accessed on 3 February 2020).

2.3. Phenotypic Characterization

Strains were phenotypically characterized by determination of cell morphology, gram staining, catalase activity, potassium hydroxide (KOH) test, consumption of calcium carbonate and oxidation of ethanol and acetic acid on ACB agar medium, as previously reported in Gullo et al. [20].
The BC production test was carried out by collecting the pellicles and boiling them in 4 mL of 5.0% NaOH solution for 2 h. Cellulose was confirmed when the pellicle did not dissolve after boiling [21]. DSMZ (German Collection of Microorganisms and Cell Cultures GmbH) strains K. xylinus (DSMZ 2004) and G. oxydans (DSMZ 3503) were used as positive and negative controls, respectively. Tests were conducted in triplicate.
The amount of BC was estimated, as reported by Gullo et al. [22]. Briefly, pellicles were washed with distilled water and then treated with 1% NaOH at 90 °C for 30 min. Treated BC was washed twice with distilled water, the pH of the residual water was determined; the washing step was repeated until reaching neutral pH. Drying was conducted at 80 °C. BC weight was measured by analytical balance (Gibertini E42S) [23]. The yield of BC produced was expressed in grams of dried BC per liter of broth.

2.4. Genomic DNA Extraction, RFLP and Amplification of (GTG)5/REP-PCR and ERIC Elements

Genomic DNA extraction was performed on liquid cultures, which were previously incubated for 72 h at +28 °C. Flasks were vigorously shaken to remove the embedded cells from BC, the liquid medium was transferred to a sterile tube and centrifuged at 12,000× g for 5 min. gDNA was extracted using the method previously proposed by Gullo et al. [24]. gDNA was checked by 3% gel agarose in 1X TBE buffer and quantified by spectrophotometric measure (NanoDrop ND-1000). Band sizes were determined using a 100 bp DNA ladder (Invitrogen, Carlsbad, CA, USA).
RFLP analysis on the full-length 16S rRNA gene was performed using RsaI and AluI endonucleases (Fermentas, Hanover, ND, USA), incubating at +37 °C for 2 h. The restriction reaction was stopped by incubation at 80 °C for 20 min. The fragments were checked on 3% gel agarose in 1X TBE buffer and sizes were determined using 100 bp DNA ladder (Invitrogen).
(GTG)5 rep-PCR fingerprinting was carried out according to the method of Versalovic et al. [25], with some minor modifications. Isolates were subjected to REP-PCR with a single oligonucleotide GTG5 (5′- GTGGTGGTGGTGGTG-3′). Samples were incubated for 5 min at 94 °C and then cycled 35 times at 94 °C for 30 s, 40 °C for 1 min, and 72 °C for 4 min. The samples were incubated for 7 min at 72 °C for final extension and kept at 4 °C. ERIC elements were initially amplified according to Versalovic et al. [26] using the primers ERIC 1R (5′-ATGTAAGCTCCTGGGGATTCAC-3′) and ERIC 2 (5′-AAGTAAGTGGGGTGAGCG-3′). Samples were incubated for 5 min at 94 °C and then cycled 30 times at 94 °C for 30 s, 57 °C for 30 s, and 65 °C for 4 min. The final extension was performed at 65 °C at 8 min and kept at 4 °C until tested. The reproducibility of ERIC/PCR was also tested by amplifying gDNA from randomly chosen strains several times. For both rep-PCR, pattern band lengths were determined by comparison against a 100 bp plus DNA ladder for the smallest bands and by 1 Kbp DNA ladder for the largest bands (Takara Bio, Inc., Otsu, Shiga, Japan).

2.5. Statistics and Clustering Analysis of (GTG)5 and ERIC Rep-PCR Patterns

Patterns obtained from (GTG)5 and ERIC rep-PCR runs were imported into BioNumerics tool, package version 8. The cluster analysis was performed calculating Dice coefficients using a tolerance of 1% and optimization of 1%. Clustering based on Dice coefficients was performed using the unpaired group method with arithmetic average (UPGMA) method. For each fingerprint assay, the cluster cut-off analysis was applied to define the most reliable clusters. The cophenetic correlation index was used as a statistical tool to evaluate the quality of branches. For each of the genotyping experiments, the discriminatory index (D), expressed as the probability that two strains consecutively taken from a sample would be placed into different clades, was calculated as described by Hunter and Gaston [27].
Correlation analysis was performed considering the most abundant clades (Clade 5, 6 and 7), in order to obtain statistically significant results. The Spearman correlation index was calculated using the Hmisc v4.4–2 R package [28] implemented in R v 4.0.3 [29]. The correlation plot was obtained using the Corrplot v 0.84 package [30].

3. Results

3.1. Kombucha Characteristics and Isolated Strains

Within 3 days of cultivation in 600 mL beakers, the two kombucha tea samples, GTK and BTK, demonstrated a thin exopolysaccharide layer that became thicker with time (from 2–3 mm to 10 mm at the end of the cultivation period). Optical microscopy observations showed a high number of free bacterial and yeast cells, as well as aggregates of cells within the matrix, making cell counting uncertain and not enlightening even from the first day of fermentation.
Titratable acidity reached a maximum of 12 g/L at the end of fermentation in the black kombucha trial (GTK) and 6 g/L in the green one (BTK). pH dropped from approximately 3.7 to 2.75 for both samples, as a result of acid formation. Although no inhibition compounds were determined in this study, the difference in the final amount of acetic acid of the GTK and BTK samples could be due to the occurrence of more antibacterial compounds in GTK [31]. Ethanol was nearly 0 at the beginning of the cultivation time and reached 0.28% (v/v) and 0.30% (v/v) at the final time, for samples GTK and BTK, respectively. The low value of observed ethanol is in agreement with the gradual increase in acetic acid during cultivation time.
The isolation was performed from preculture (P1 and P2) and samples at 0, 6, and 12 days of cultivation. A total of 122 isolates were selected based on the colony and/or cellular morphologies. Microscopically, presumptive acetic acid bacteria cells appeared in single or grouped pairs or short chains, Gram-negative, catalase-positive, and KOH positive. A total of 46 isolates, which oxidized both ethanol and acetic acid, as confirmed by the ethanol chalk–ethanol test on GYC and the colorimetric assay on ACB, and produced a well-defined cellulose layer and/or formed both a cellulose layer and cellulose particles dispersed in the surrounding liquid, were considered in this study (Table 1).
Qualitative cellulose tests confirmed their ability to produce BC. Significant macroscopic diversity in the cellulosic material was observed, coupled with high variability in terms of amount of BC, as reported in Figure 1. The strains with a BC production level of at least higher than 4 g/L were deposited in the UMCC culture collection with the code reported in Table 1.

3.2. Diversity of Acetic Acid Bacteria from Kombucha Tea

The RFLP analysis on 16S rRNA, using RsaI and AluI restriction endonucleases, was conducted as a preliminary test to investigate the diversity of isolates recovered from kombucha tea samples. The mapping analysis based on the size of the fragments obtained by RsaI and AluI showed high homogeneous results, indicating the isolates as members of the same species. This evidence was in agreement with a previous study [16] which grouped Komagataeibacter strains according to RsaI and AluI restriction enzymes and observed similarities between the pattern of K. xylinus DSMZ 2004 with the Komagataeibacter strains investigated in this study. The low variability of the fingerprinting pattern is also consistent with previous studies on the high similarity of the 16S rRNA gene among the Komagataeibacter genus and especially within the K. xylinus species [32,33].
From the genotyping of strains performed by using (GTG)5 and ERIC fingerprinting techniques, two dendrograms were generated (Figure 2). Based on the cluster cut-off analysis, profiles obtained using the amplification of interspersed tandem repeats GTG (Figure 2) were capable of discriminating between the strains considering a minimum percentage of similarity of 27.9%. A total of five major clades were created with a discrimination power of 0.42, calculated using the D index (Simpson index) [27]. In order to improve the discrimination power of the analysis, we considered a minimum percentage of similarity of 70%. At this percentage of similarity, the resulted discrimination power was 0.94, defining 20 different biotypes. Among the detected clades, eight of them included just one isolate while seven were formed by two isolates. It was possible to distinguish two major clades, of which, one included nine isolates (clade 18) and another included five isolates (clade 19). The two reference strains, represented by K1G4 (=UMCC 2947) and K2G30 (=UMCC 2756) clustered into two different clades, represented by clades 19 and 7, respectively. The similarity percentage among the two references was 40.2%.
The phylogenetic tree obtained from the digitized pattern profile using ERIC rep-PCR is represented in Figure 3. Based on the clustering cut-off analysis, ERIC was able to discriminate using a minimum similarity percentage of 40%, identifying five major clades, as in the case of (GTG)5. The clustering cut-off analysis discriminated a total of 19 biotypes, with a discrimination power of 0.95, slightly higher compared to (GTG)5. A total of six clades were represented by just one isolate, while the remaining isolates were well distributed among the detected clades. Three clades were represented by five isolates (clade 7, 16, and 17) and two clusters by four isolates (clade 8 and 9). The remaining clades were represented by three and two isolates. The number of detected clades was not so different compared to (GTG)5 and, also in this case, the reference strains (K1G4 and K2G30) were found to be separated into two different groups (clade 15 for K1G4 and clade 12 for K2G30). The similarity of Dice indexes based on the patterns obtained from ERIC rep-PCR was 56.7%; quite similar to that obtained using (GTG)5 rep-PCR (40.2%).

3.3. Improving the Fingerprinting of Komagataeibacter Strains

In order to improve the sensitivity of the intraspecific discrimination, the data obtained from both rep-PCR assays were combined. The clustering analysis is reported in Figure 4. Based on the cluster cut-off analysis, the combination of (GTG)5 and ERIC rep-PCR allowed for discriminating the isolates, grouping them into 22 clades, considering a similarity percentage threshold of 94%. The discrimination power, calculated as the Simpson index, was 0.94, mining a high intraspecies diversity. Most of the isolates were grouped into three clusters (clade 5, 6, and 7) representing 43.5% of the isolates (n = 20). The remaining isolates (n = 26) were dispersed in 19 clades, of which, most of them consisted of just one isolate (n = 12). The reference strains were clustered into two different clades, as observed in both the case of (GTG)5 and ERIC rep-PCR. Clade 10, in which K1G4 was clustered with K1G6, and clade 19, where K2G30 clustered alone, had a similarity percentage of Dice index of 84.9%, higher considering the separated results from (GTG)5 and ERIC rep-PCR. This result is in agreement with the meaning of the combined analysis, since the two reference strains belong to the species of K. xylinus, as previously stated [14,34].

3.4. Correlation between Strains and Cellulose Yield

In order to understand putative correlation of the BC yield obtained from isolates with time and environmental factors, such as sampling fraction (pellicle or liquid) and clustering methods applied in this study, a Spearman correlation index was calculated. At the beginning all clades that resulted from the combined fingerprinting approach (Figure 4) were compared with all variables (Figure 5a). All factors seemed to strongly influence the BC yields. Only the time of sampling seemed to slightly influence BC production, but with a low value (R2 = −0.38, p = 0.009). In order to understand how the sampling time could affect BC yield, the three most abundant clades, represented by clade 5, 6 and 7 were evaluated by calculating the Spearman coefficient, considering both time and sampling fraction. Regarding clade 5 (n = 6) (Figure 5b), the time seems to correlate positively (R2 = 0.84, p = 0.03). Indeed, the amount of BC (Figure 1) was higher for isolates sampled at day 12 (K2A46, K2A47, K2G46 and K1A34) and lower for isolates sampled at day 0 and from the preculture (K1G5 and K1G3 respectively). In clade 6 (n = 8) (Figure 5c), BC yield seems to be not statistically correlated by the sampling time (R2 = −0.70, p > 0.05). Regarding clade 7 (n = 6) (Figure 5d), the sampling time negatively affects the BC amount (R2 = −0.92, p = 0.009).

4. Discussion

Bacterial cellulose is receiving a great deal of attention due to several features that make it suitable for many biotechnological applications, such as biomedical devices, food packaging components and engineering materials [1,2,35].
It is well known from the literature that kombucha tea represents one of the most abundant reservoirs of BC producer bacteria, in particular, strains belonging to the Komagataeibacter genus [10,12].
In this study, kombucha tea was chosen as a selective source for recovering cellulose producing bacteria. A total of 46 isolates were selected as the main components of the liquid and pellicle fractions of GTK and BTK kombucha tea samples and their phenotypic and genotypic variability was evaluated. RFLP of 16S rRNA gene and an analysis of the interspecific variability of isolates was conducted, and the results were discussed, considering the available genome sequences of strains K1G4 and K2G30.
The data obtained by RFLP analysis showed a very low variability in terms of the number of different species isolated from kombucha tea. This result is in agreement with previous studies reflecting both the low species diversity in kombucha tea [10], as well as the difficulty of taxonomic assignment of Acetobacteraceae by using 16S rRNA gene, due to its similarity in high phylogenetic related organisms, as extensively reported in the literature [36].
Considering the analysis of the intraspecific variability within bacteria, usually, genotyping is performed using the amplification of highly conserved tandem repeat regions that are widely dispersed in the genome. The amplification of repeated elements of genomic DNA, such as (GTG)5 and ERIC sequences, offers some advantages in terms of costs and rapidity, and is well developed for several Gram-negative and Gram-positive bacterial species, such as Lactobacilli [37] and Enterococci [38] and acetic acid bacteria [39].
In this study, the clustering analysis using the digitized patterns of (GTG)5 elements revealed the clustering of the 46 isolates into five large clusters considering a cluster cut-off value of 50% in similarity. The discrimination power of this analysis was 0.94 considering a similarity threshold of 70%. At this similarity value, a total of 20 different biotypes were detected. Similar results were obtained using the amplification of ERIC repeated elements. Considering the cluster cut-off analysis, the comparison of the digitized patterns of ERIC rep-PCR showed a clustering of the 46 isolates with a similarity percentage of 50%, the same result obtained using (GTG)5 rep-PCR. The discrimination power of ERIC fingerprint was 0.95, very similar to the (GTG)5 rep-PCR, considering a similarity percentage of 70. Based on these results, a total of 19 different clades were identified, compared to the 20 obtained using (GTG)5 rep-PCR. In both the fingerprint approaches, K1G4 and K2G30, previously identified as K. xylinus [14,34], were considered as reference strains and they were clustered into two different clades. Among the two reference strains, the percentage of similarity was 44.2% in the case of (GTG)5 rep-PCR and 56.7% for ERIC rep-PCR. These results support the hypothesis that most of the strains analyzed belong to the K. xylinus species and that the discrimination was performed at the strain level.
In order to improve the genotyping analysis, a combined model of the two fingerprint assays was performed (Figure 3). The resulting analysis showed a clusterization of the isolates into seven clades based on the cluster cut-off value of 88% of pattern similarity. The discrimination power of this analysis was 0.94 at 94% of similarity. The number of biotypes defined by the combined model was 22; meaning that an improvement of the discrimination of the isolates was achieved compared to the single results obtained by the two approaches. Additionally, concerning the reference strains (K1G4 and K2G30) in this case, we observed a grouping into two different clades, with a similarity percentage of 84.9%, meaning that all the isolates with equal or major than that value possibly belong to the species of K. xylinus. Based on this consideration, all of the isolates could be assigned to the K. xylinus except one isolate, the K2G44, clustering out of the group, with a similarity percentage of 78.8%. The data were in agreement with previous studies of the fingerprinting analysis of acetic acid bacteria [39,40], in which the differences of strains among the same species (as Komagataeibacter xylinus) were about 80%. In this case, combining the two methods, the detection of different strains was improved.
The correlation analysis deriving from the combined strain clustering allowed us to link the BC yield with isolation time and environmental variables. A possible explanation of the obtained results may lie in the genetic plasticity features of acetic acid bacteria. It is well known in the literature that acetic acid bacteria are able to adapt themselves to environmental stressors, such as temperature and high concentration of organic acids, by acquiring genetic material from the environment or losing it. Indeed, acetic acid bacteria possess a large number of transposons genes and plasmids [41]. One example is provided by Acetobacter species which strains adapt themselves to environmental conditions along cultivation and fermentation cycles [42,43]. In Komagataeibacter species the genetic plasticity has not yet been evaluated. However, we believe this aspect should be elucidated in the light to better understand the variability of BC yield within strains of Komagataeibacter genus.
The kombucha tea samples analyzed in this study allowed us to recover a pool of acetic acid bacteria strains, which have been characterized and deposited into the UMCC culture collection as promising candidates for BC production.

Author Contributions

M.G. conceptualized the research, provided resources and supervised the work; D.M. performed experiments; S.L.C. performed genotyping analysis and wrote part of the original draft; L.D.V. reviewed and edited the manuscript; K.A. wrote part of the original draft; M.B. participated in the investigation, reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Part of this research was granted by The Research Doctorate School in Food and Agricultural Science, Technology and Biotechnology (STEBA).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The National Institute of Applied Sciences and Technology (INSAT, Tunisia), is acknowledged for providing Kombucha samples. A special thanks to The Research Doctorate School in Food and Agricultural Science, Technology and Biotechnology (STEBA) is acknowledged for assigning the “Michela Stanca” award. The JRU MIRRI-IT is greatly acknowledged for technical and scientific support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Robotti, F.; Bottan, S.; Fraschetti, F.; Mallone, A.; Pellegrini, G.; Lindenblatt, N.; Starck, C.; Falk, V.; Poulikakos, D.; Ferrari, A. A micron-scale surface topography design reducing cell adhesion to implanted materials. Sci. Rep. 2018. [Google Scholar] [CrossRef]
  2. Haghighi, H.; Gullo, M.; La China, S.; Pfeifer, F.; Siesler, H.W.; Licciardello, F.; Pulvirenti, A. Characterization of bio-nanocomposite films based on gelatin/polyvinyl alcohol blend reinforced with bacterial cellulose nanowhiskers for food packaging applications. Food Hydrocoll. 2020, 106454. [Google Scholar] [CrossRef]
  3. Barbi, S.; Taurino, C.; La China, S.; Anguluri, K.; Gullo, M.; Montorsi, M. Mechanical and structural properties of environmental green composites based on functionalized bacterial cellulose. Cellulose 2021, 1–7. [Google Scholar] [CrossRef]
  4. Valera, M.J.; Torija, M.J.; Mas, A. Detrimental effects of acetic acid bacteria in foods. In Acetic Acid Bacteria: Fundamentals and Food Applications, 1st ed.; Sengun, I., Ed.; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: Boca Raton, FL, USA, 2017; Volume 10, pp. 299–320. [Google Scholar]
  5. Gullo, M.; La China, S.; Falcone, P.M.; Giudici, P. Biotechnological production of cellulose by acetic acid bacteria: Current state and perspectives. Appl. Microbiol. Biotechnol. 2018, 1–14. [Google Scholar] [CrossRef] [PubMed]
  6. Son, H.J.; Heo, M.S.; Kim, Y.G.; Lee, S.J. Optimization of fermentation conditions for the production of bacterial cellulose by a newly isolated Acetobacter sp. A9 in shaking cultures. Biotechnol. Appl. Biochem. 2001, 33, 1. [Google Scholar] [CrossRef] [PubMed]
  7. La China, S.; Zanichelli, G.; De Vero, L.; Gullo, M. Oxidative fermentations and exopolysaccharides production by acetic acid bacteria: A mini review. Biotechnol. Lett. 2018, 40, 1289–1302. [Google Scholar] [CrossRef]
  8. Coton, M.; Pawtowski, A.; Taminiau, B.; Burgaud, G.; Deniel, F.; Coulloumme-Labarthe, L.; Fall, A.; Daube, G.; Coton, E. Unraveling microbial ecology of industrial-scale Kombucha fermentations by metabarcoding and culture-based methods. FEMS Microbiol. Ecol. 2017, 93. [Google Scholar] [CrossRef]
  9. Laureys, D.; Britton, S.J.; De Clippeleer, J. Kombucha tea fermentation: A review. J. Am. Soc. Brew. Chem. 2020, 78, 165–174. [Google Scholar] [CrossRef]
  10. Sievers, M.; Lanini, C.; Weber, A.; Schuler-Schmid, U.; Teuber, M. Microbiology and fermentation balance in a Kombucha beverage obtained from a tea fungus fermentation. Syst. Appl. Microbiol. 1995, 18, 590–594. [Google Scholar] [CrossRef]
  11. Arıkan, M.; Mitchell, A.L.; Finn, R.D.; Gürel, F. Microbial composition of Kombucha determined using amplicon sequencing and shotgun metagenomics. J. Food Sci. 2020, 85, 455–464. [Google Scholar] [CrossRef]
  12. Villarreal-Soto, S.A.; Beaufort, S.; Bouajila, J.; Souchard, J.P.; Taillandier, P. Understanding Kombucha tea fermentation: A review. J. Food Sci. 2018, 83, 580–588. [Google Scholar] [CrossRef]
  13. Gaggìa, F.; Baffoni, L.; Galiano, M.; Nielsen, D.; Jakobsen, R.; Castro-Mejía, J.; Bosi, S.; Truzzi, F.; Musumeci, F.; Dinelli, G.; et al. Kombucha beverage from green, black and rooibos teas: A comparative study looking at microbiology, chemistry and antioxidant activity. Nutrients 2018, 11, 1. [Google Scholar] [CrossRef] [Green Version]
  14. Gullo, M.; La China, S.; Petroni, G.; Di Gregorio, S.; Giudici, P. Exploring K2G30 genome: A high bacterial cellulose producing strain in glucose and mannitol based media. Front. Microbiol. 2019, 10, 1–12. [Google Scholar] [CrossRef]
  15. Singhsa, P.; Narain, R.; Manuspiya, H. Physical structure variations of bacterial cellulose produced by different Komagataeibacter xylinus strains and carbon sources in static and agitated conditions. Cellulose 2018, 25, 1571–1581. [Google Scholar] [CrossRef]
  16. Mamlouk, D.; Gullo, M. Acetic acid bacteria: Physiology and carbon sources oxidation. Indian J. Microbiol. 2013, 53, 377–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Jayabalan, R.; Marimuthu, S.; Swaminathan, K. Changes in content of organic acids and tea polyphenols during kombucha tea fermentation. Food Chem. 2007, 102, 392–398. [Google Scholar] [CrossRef]
  18. Mamlouk, D. Insight Into Physiology and Functionality of Acetic Acid Bacteria Through a Multiphasic Approach. Ph.D. Thesis Dissertation, University of Modena and Reggio Emilia, Reggio Emilia, Italy, 2012. [Google Scholar]
  19. De Vero, L.; Boniotti, M.B.; Budroni, M.; Buzzini, P.; Cassanelli, S.; Comunian, R.; Gullo, M.; Logrieco, A.F.; Mannazzu, I.; Musumeci, R.; et al. Preservation, characterization and exploitation of microbial biodiversity: The perspective of the Italian Network of Culture Collections. Microorganisms 2019, 7, 685. [Google Scholar] [CrossRef] [Green Version]
  20. Gullo, M.; Mamlouk, D.; De Vero, L.; Giudici, P. Acetobacter pasteurianus strain AB0220: Cultivability and phenotypic stability over 9 years of preservation. Curr. Microbiol. 2012, 64, 576–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Navarro, R.; Uchimura, T.; Komagata, K. Taxonomic heterogeneity of strains comprising Gluconacetobacter hansenii. J. Gen. Appl. Microbiol. 1999, 45, 295–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Gullo, M.; Sola, A.; Zanichelli, G.; Montorsi, M.; Messori, M.; Giudici, P. Increased production of bacterial cellulose as starting point for scaled-up applications. Appl. Microbiol. Biotechnol. 2017, 101, 8115–8127. [Google Scholar] [CrossRef] [PubMed]
  23. Hwang, J.W.; Yang, Y.K.; Hwang, J.K.; Pyun, Y.R.; Kim, Y.S. Effects of pH and dissolved oxygen on cellulose production by Acetobacter xylinum BRC5 in agitated culture. J. Biosci. Bioeng. 1999, 88, 183–188. [Google Scholar] [CrossRef]
  24. Gullo, M.; De Vero, L.; Giudici, P. Succession of selected strains of Acetobacter pasteurianus and other acetic acid bacteria in traditional balsamic vinegar. Appl. Environ. Microbiol. 2009, 75, 2585–2589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Versalovic, J.; Schneider, M.; De Bruijn, F.J.; Lupski, J.R. Genomic fingerprinting of bacteria using repetitive sequence-based polymerase chain reaction. Methods Mol. Cell. Biol. 1994, 5, 25–40. [Google Scholar]
  26. Versalovic, J.; Koeuth, T.; Lupski, R. Distribution of repetitive DNA sequences in eubacteria and application to fingerprinting of bacterial genomes. Nucleic Acids Res. 1991, 19, 6823–6831. [Google Scholar] [CrossRef] [PubMed]
  27. Hunter, P.R.; Gaston, M.A. Numerical index of the discriminatory ability of typing systems: An application of Simpson’s index of diversity. J. Clin. Microbiol. 1988, 26, 2465–2466. [Google Scholar] [CrossRef] [Green Version]
  28. Harrell, F.E., Jr.; Dupont, M.C. R Package Hmisc; CRAN; R Foundation for Statistical Computing: Vienna, Austria, 2014; Available online: https://CRAN.R-project.org/package=Hmisc (accessed on 3 February 2020).
  29. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 3 February 2020).
  30. Taiyun, W.; Viliam, S. R package “corrplot”: Visualization of a Correlation Matrix: Version 0.84. 2017. Available online: https://github.com/taiyun/corrplot (accessed on 3 February 2020).
  31. Battikh, H.; Chaieb, K.; Bakhrouf, A.; Ammar, E. Antibacterial and antifungal activities of black and green kombucha teas. J. Food Biochem. 2013. [Google Scholar] [CrossRef]
  32. Sievers, M.; Sellmer, S.; Teuber, M. Acetobacter europaeus sp. nov., a main component of industrial vinegar fermenters in central Europe. Syst. Appl. Microbiol. 1992, 15, 386–392. [Google Scholar] [CrossRef]
  33. Ludwig, W.; Klenk, H.P. Overview: A phylogenetic backbone and taxonomic framework for procaryotic systematics. In Bergey’s Manual of Systematic Bacteriology, The Proteobacteria, Part A, Introductory Essays, 2nd ed.; Brenner, D.J., Krieg, N.R., Staley, J.T., Garrity, G.M., Eds.; Springer: New York, NY, USA, 2005; Volume 2, pp. 49–65. [Google Scholar]
  34. La China, S.; Bezzecchi, A.; Moya, F.; Petroni, G.; Di Gregorio, S.; Gullo, M. Genome sequencing and phylogenetic analysis of K1G4: A new Komagataeibacter strain producing bacterial cellulose from different carbon sources. Biotechnol. Lett. 2020, 42, 807–818. [Google Scholar] [CrossRef]
  35. Jakmuangpak, S.; Prada, T.; Mongkolthanaruk, W.; Harnchana, V.; Pinitsoontorn, S. Engineering bacterial cellulose films by nanocomposite approach and surface modification for biocompatible triboelectric nanogenerator. ACS Appl. Electron. Mater. 2020, 2, 2498–2506. [Google Scholar] [CrossRef]
  36. Cleenwerck, I.; Vandemeulebroecke, K.; Janssens, D.; Swings, J. Re-examination of the genus Acetobacter, with descriptions of Acetobacter cerevisiae sp. nov. and Acetobacter malorum sp. nov. Int. J. Syst. Evol. Microbiol. 2002, 52, 1551–1558. [Google Scholar] [CrossRef] [Green Version]
  37. Gevers, D.; Huys, G.; Swings, J. Applicability of rep-PCR fingerprinting for identification of Lactobacillus species. FEMS Microbiol. Lett. 2001, 205, 31–36. [Google Scholar] [CrossRef] [Green Version]
  38. Švec, P.; Vancanneyt, M.; Seman, M.; Snauwaert, C.; Lefebvre, K.; Sedláček, I.; Swings, J. Evaluation of (GTG)5-PCR for identification of Enterococcus spp. FEMS Microbiol. Lett. 2005, 247, 59–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. De Vuyst, L.; Camu, N.; De Winter, T.; Vandemeulebroecke, K.; Van De Perre, V.; Vancanneyt, M.; De Vos, P.; Cleenwerck, I. Validation of the (GTG)5-rep-PCR fingerprinting technique for rapid classification and identification of acetic acid bacteria, with a focus on isolates from Ghanaian fermented cocoa beans. Int. J. food Microbiol. 2007. [Google Scholar] [CrossRef]
  40. Cleenwerck, I.; De Wachter, M.; González, Á.; De Vuyst, L.; De Vos, P. Differentiation of species of the family Acetobacteraceae by AFLP DNA fingerprinting: Gluconacetobacter kombuchae is a later heterotypic synonym of Gluconacetobacter hansenii. Int. J. Syst. Evol. Microbiol. 2009, 59, 1771–1786. [Google Scholar] [CrossRef] [PubMed]
  41. Azuma, Y.; Hosoyama, A.; Matsutani, M.; Furuya, N.; Horikawa, H.; Harada, T.; Hirakawa, H.; Kuhara, S.; Matsushita, K.; Fujita, N.; et al. Whole-genome analyses reveal genetic instability of Acetobacter pasteurianus. Nucleic Acids Res. 2009, 37, 5768–5783. [Google Scholar] [CrossRef] [Green Version]
  42. Matsutani, M.; Nishikura, M.; Saichana, N.; Hatano, T.; Masud-Tippayasak, U.; Theergool, G.; Yakushi, T.; Matsushita, K. Adaptive mutation of Acetobacter pasteurianus SKU1108 enhances acetic acid fermentation ability at high temperature. J. Biotechnol. 2013, 165, 109–119. [Google Scholar] [CrossRef] [PubMed]
  43. Pothimon, R.; Gullo, M.; La China, S.; Thompson, A.K.; Krusong, W. Conducting High acetic acid and temperature acetification processes by Acetobacter pasteurianus UMCC 2951. Process Biochem. 2020, 98, 41–50. [Google Scholar] [CrossRef]
Figure 1. Bacterial cellulose yields of selected strains isolated from Kombucha tea. Each value is the mean of three parallel replicates ± standard deviation.
Figure 1. Bacterial cellulose yields of selected strains isolated from Kombucha tea. Each value is the mean of three parallel replicates ± standard deviation.
Applsci 11 01595 g001
Figure 2. Dendrograms obtained from UPGMA analysis, using Dice’s coefficient, of the digitized patterns obtained from (GTG)5 rep-PCR. The similarity threshold for biotypes discrimination was 70%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality.
Figure 2. Dendrograms obtained from UPGMA analysis, using Dice’s coefficient, of the digitized patterns obtained from (GTG)5 rep-PCR. The similarity threshold for biotypes discrimination was 70%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality.
Applsci 11 01595 g002
Figure 3. Dendrograms obtained from UPGMA analysis, using Dice’s coefficient, of the digitized patterns obtained from ERIC rep-PCR. The similarity threshold for biotypes discrimination was 70%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality.
Figure 3. Dendrograms obtained from UPGMA analysis, using Dice’s coefficient, of the digitized patterns obtained from ERIC rep-PCR. The similarity threshold for biotypes discrimination was 70%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality.
Applsci 11 01595 g003
Figure 4. Dendrogram obtained combining the digitalized patterns from both (GTG)5 rep-PCR and ERIC rep-PCR. The dendrogram was drawn from UPGMA analysis using Dice’s coefficient. The discrimination of biotypes was performed considering a similarity threshold of 94%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality. Cluster colors were defined in agreement with the Figure 1.
Figure 4. Dendrogram obtained combining the digitalized patterns from both (GTG)5 rep-PCR and ERIC rep-PCR. The dendrogram was drawn from UPGMA analysis using Dice’s coefficient. The discrimination of biotypes was performed considering a similarity threshold of 94%. The cophenetic coefficient is represented by numbers and dots colored as red–orange–yellow–green, based on the branch quality. Cluster colors were defined in agreement with the Figure 1.
Applsci 11 01595 g004
Figure 5. Correlation plot of the most abundant clades. Spearman correlation among time, environmental variables and obtained yields. (a) All retrieved clades were compared with time and environmental variables. Yield obtained within (b) clade 5, (c) clade 6 were compared with time and fraction variables. For clade 7 (d), includes only isolates sampled from liquid fraction.
Figure 5. Correlation plot of the most abundant clades. Spearman correlation among time, environmental variables and obtained yields. (a) All retrieved clades were compared with time and environmental variables. Yield obtained within (b) clade 5, (c) clade 6 were compared with time and fraction variables. For clade 7 (d), includes only isolates sampled from liquid fraction.
Applsci 11 01595 g005
Table 1. Features of 46 isolates from Kombucha tea used in this study. Site of isolation and RFLP analysis. In bold are reported the strains deposited in UNIMORE Microbial Culture Collection (UMCC) collection; within the brackets is the assigned accession number.
Table 1. Features of 46 isolates from Kombucha tea used in this study. Site of isolation and RFLP analysis. In bold are reported the strains deposited in UNIMORE Microbial Culture Collection (UMCC) collection; within the brackets is the assigned accession number.
* StrainSampleTimeBC (g/L) Colony MorphologyLiquid GrowthAluI (pb)GroupRsaI (bp)GroupReference
K1A18liquid64.0448 ± 0.0008 green cellulosiccellulose layer--120, 400, 425, 550AThis study
K1A34liquid125.6214 ± 0.0010 green with light cellulosic edgecellulose layer150, 450, 8001120, 400, 425, 550AThis study
K1A6liquid03.0024 ± 0.0009 cellulosic smallcellulose layer150, 450, 8001120, 400, 425, 550AThis study
K1A7liquid01.5886 ± 0.0004 cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550AThis study
K1A8liquid02.5886 ± 0.0002 cellulosiccellulose layer150, 450, 8001120, 400, 425, 550AThis study
K1A9liquid03.4014 ± 0.0004 green cellulosiccellulose layer----This study
K1G2
(=UMCC 2964)
liquidP14.7895 ± 0.0014 white cellulosic edge cellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K1G22liquid120.7057 ± 0.0003 cellulose lightcellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K1G23liquid124.8122 ± 0.0003dark cellulosecellulose layer150, 450, 8001120, 400, 425, 550A[20]
K1G24liquid124.1210 ± 0.0001white celluloseonly cellulose deposit--120, 400, 425, 550A[20]
K1G3liquidP12.1857 ± 0.0052 light beige cellulose layer150, 450, 8001120, 400, 425, 550A[20]
K1G4liquid01.1738 ± 0.0032 cellulosic cellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K1G5liquid04.9590 ± 0.0027 cellulosic cellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K1G6liquid05.3676 ± 0.0025 white cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2A10
(=UMCC 2965)
liquid05.0405 ± 0.0008 cellulosic green smallcellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A28liquid65.8962 ± 0.0007 green cellulosiccellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A32pellicle65.6143 ± 0.0001 green cellulosiccellulose layer----This study
K2A33liquid60.6414 ± 0.0005 green cellulosiccellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A44liquid125.1971 ± 0.0002 small green cellulose layer/deposit150, 450, 8001120, 400, 425, 550AThis study
K2A45pellicle125.2133 ± 0.0004 light edgecellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A46liquid125.7833 ± 0.0001 dark center light edgecellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A47pellicle125.8019 ± 0.0001 blue medium displayedcellulose layer150, 450, 8001120, 400, 425, 550AThis study
K2A7liquid04.6076 ± 0.0007 green cellulosic smallcellulose layer150, 450, 8001400, 450, 550BThis study
K2A8liquid04.4781 ± 0.0004 green cellulosic smallcellulose layer----This study
K2G1liquidP26.8843 ± 0.0004 cellulosic cellulose layer/deposit150, 450, 8001120, 400, 425, 550AThis study
K2G10liquid04.9752 ± 0.0016 cellulosic smallcellulose layer/deposit150, 450, 8001400, 450, 550B[20]
K2G11liquid01.0638 ± 0.0003 white cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G12liquid06.2805 ± 0.0001 small cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G14liquid07.7414 ± 0.0005 white cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G15liquid07.5448 ± 0.0001 creamycellulose layer/deposit150, 450, 8001400, 450, 550A[20]
K2G2pellicleP26.2548 ± 0.0004 cellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G29
(=UMCC 2967)
pellicle67.5848 ± 0.0002 cellulosiccellulose layer150, 450, 8001400, 450, 550B[20]
K2G30pellicle623.1829 ± 0.0001 cellulosiccellulose layer/deposit 150, 450, 8001400, 450, 550B[20]
K2G31pellicle61.1267 ± 0.0007 cellulosiccellulose layer/deposit----[20]
K2G32
(=UMCC 2968)
liquid69.9976 ± 0.0003 cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G34liquid66.6614 ± 0.0001 cellulosiccellulose layer----[20]
K2G35liquid61.9314 ± 0.0005 cellulosiccellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G36liquid60.5929 ± 0.0004 yellow cellulosic cellulose layer150, 450, 8001400, 450, 550B[20]
K2G38
(=UMCC 2969)
liquid69.6957 ± 0.0009 cellulosiccellulose layer150, 450, 8001120, 400, 425, 550A[20]
K2G39
(=UMCC 2970)
liquid610.3957 ± 0.0001 cellulosiccellulose 150, 450, 8001400, 450, 550B[20]
K2G41
(=UMCC 2971)
liquid69.1443 ± 0.0003 cellulosiccellulose layer150, 450, 8001400, 450, 550B[20]
K2G44
(=UMCC 2972)
pellicle125.9500 ± 0.0005 dark brown cellulosecellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G46liquid129.1195 ± 0.0007 brown cellulosecellulose layer/deposit----[20]
K2G5
(=UMCC 2966)
liquidP28.0605 ± 0.0003 white cellulosiccellulose layer--120, 400, 425, 550A[20]
K2G6liquidP27.7533 ± 0.0003 white cellulosic cellulose layer/deposit150, 450, 8001120, 400, 425, 550A[20]
K2G8liquid06.8343 ± 0.0002 white cellulosiccellulose layer150, 450, 8001120, 400, 425, 550A[20]
(-) not detected; * Strain labeling: K (Kombucha tea); 1 (BTK); 2 (GTK); A (ACB medium); G (GYC medium); Number: progressive number. BC: bacterial cellulose; BTK: black tea kombucha; GTK: green tea kombucha.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

La China, S.; De Vero, L.; Anguluri, K.; Brugnoli, M.; Mamlouk, D.; Gullo, M. Kombucha Tea as a Reservoir of Cellulose Producing Bacteria: Assessing Diversity among Komagataeibacter Isolates. Appl. Sci. 2021, 11, 1595. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041595

AMA Style

La China S, De Vero L, Anguluri K, Brugnoli M, Mamlouk D, Gullo M. Kombucha Tea as a Reservoir of Cellulose Producing Bacteria: Assessing Diversity among Komagataeibacter Isolates. Applied Sciences. 2021; 11(4):1595. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041595

Chicago/Turabian Style

La China, Salvatore, Luciana De Vero, Kavitha Anguluri, Marcello Brugnoli, Dhouha Mamlouk, and Maria Gullo. 2021. "Kombucha Tea as a Reservoir of Cellulose Producing Bacteria: Assessing Diversity among Komagataeibacter Isolates" Applied Sciences 11, no. 4: 1595. https://0-doi-org.brum.beds.ac.uk/10.3390/app11041595

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