Bacterial and Fungal Communities of Table Grape Skins in Shanghai

: Grape skin serves as a reservoir for many bacteria and fungi, which can affect grape health, quality, and safety. However, grape skin microbiota and mycobiota of table grapes remain largely understudied. This study investigated bacterial and fungal communities residing on different table grapevine cultivars (‘Summer Black’ and ‘Kyoho’) as well as the impact of potential contributors thereby, via culture-dependent and culture-independent (analysis of 16S rRNA gene and internal transcribed spacer sequences) methods with different purposes. Microbiota of both grapevine cultivars were dominated by Cladosporium , Alternaria , Aspergillus , Thauera , and Pantoea . In addition, yeast strains belonging to Hanseniaspora opuntiae , Pichia terricola , Rhodotorula mucilaginosa , Candida stellimalicola , and Kodamaea ohmeri were enriched from the studied grapes, while some strains were considered as health-threatening pathogens. Differences in grapevine cultivars did not significantly affect their mycobiota and microbiota profiles. Nevertheless, their mycobiota exhibited significant variations across different grape-sampling sites in Shanghai, indicating the contribution of the grape-growing environment to grape skin mycobiota. Altogether, the current study demonstrated the contribution of the grape-growing environment to table grape skin mycobiota, and highlighted the importance of microbiota management in the production and consumption of table grapes.


Introduction
Grapevines (Vitis spp.) are one of the most widely cultivated fruits worldwide and are consumed fresh, dried, or processed into a variety of products, for instance wine, juice, and vinegar [1].Table grapes are one type of grapes that are often consumed fresh.Among popular table grapes, 'Summer Black' [2] and 'Kyoho' [3] are two of the most popular table grape varieties.They are known for their unique flavor, aroma, nutritional value, and health benefits, which, in part, can be attributed to the high phenolic content in their grape skins and their potential use to reduce the risk of cardiovascular disease and inflammation [4].
Grapes harbor a variety of bacteria and fungi that play an important role in grape health, quality, and food safety [5].Depending, in part, on the presence of beneficial or potentially harmful microbes on the grape, inhabiting bacteria and/or fungi could be beneficial and/or potentially harmful to grapes or consumers [5].The composition and dynamics of inhabiting bacteria and fungi can be affected by various factors, such as grapevine cultivars, geographic origin, climate, agricultural practices, and postharvest handling [6][7][8].Understanding these factors and their effects on grape skin bacteria and fungi is essential for developing effective strategies to manage grape quality and safety.
In recent years, advances in high-throughput sequencing technologies have enabled more comprehensive characterization of microbial communities in various environments, including grape skins [9].These methods, such as 16S rRNA gene sequencing and internal transcribed spacer (ITS) sequencing, allow for the identification and quantification of bacterial and fungal taxa present in complex samples [10,11].Nevertheless, to date, available studies investigating grape skin bacteria and fungi have mainly focused on wine-making grapes [5,9,12,13].Relatively little is known about the bacterial and fungal communities of table grapes, which are often harvested at a different stage of ripeness and sweetness compared to wine making grapes.
The current study aimed to characterize the bacterial and fungal communities of table grapes using 'Summer Black' and 'Kyoho' grapes as representatives.Based on cultureindependent methods, we studied grape skins' microbiota and mycobiota profiles, as well as the effect of grapevine cultivars, grape-sampling districts, and viticultural conditions thereof.Using a cultivation-dependent method, we isolated yeasts that could grow under given conditions.We hypothesized that grapevine cultivars, grape-sampling districts and viticultural conditions could contribute to variation in bacterial and fungal communities.

Study Area and Grape Sampling
In total, sixty bunches of undamaged table grapes were obtained from Jinshan, Pudong, and Jiading, which are key table grape production districts in Shanghai.Within each sampling district, grapes were collected from one open-air vineyard and one greenhouse vineyard.The distance from each sampling site to another is shown based on the projected coordinate system (WGS 1984 UTM zone 51N) in ArcGIS 10.8.2 (Figure 1).
including grape skins [9].These methods, such as 16S rRNA gene sequencing and internal transcribed spacer (ITS) sequencing, allow for the identification and quantification of bacterial and fungal taxa present in complex samples [10,11].Nevertheless, to date, available studies investigating grape skin bacteria and fungi have mainly focused on wine-making grapes [5,9,12,13].Relatively little is known about the bacterial and fungal communities of table grapes, which are often harvested at a different stage of ripeness and sweetness compared to wine making grapes.
The current study aimed to characterize the bacterial and fungal communities of table grapes using 'Summer Black' and 'Kyoho' grapes as representatives.Based on cultureindependent methods, we studied grape skins' microbiota and mycobiota profiles, as well as the effect of grapevine cultivars, grape-sampling districts, and viticultural conditions thereof.Using a cultivation-dependent method, we isolated yeasts that could grow under given conditions.We hypothesized that grapevine cultivars, grape-sampling districts and viticultural conditions could contribute to variation in bacterial and fungal communities.

Study Area and Grape Sampling
In total, sixty bunches of undamaged table grapes were obtained from Jinshan, Pudong, and Jiading, which are key table grape production districts in Shanghai.Within each sampling district, grapes were collected from one open-air vineyard and one greenhouse vineyard.The distance from each sampling site to another is shown based on the projected coordinate system (WGS 1984 UTM zone 51N) in ArcGIS 10.8.2 (Figure 1).Within each vineyard, two grapevine cultivars, namely 'Summer Black' (five bunches) and 'Kyoho' (five bunches) were collected.Obtained grapes were subjected to yeast cultivation and sequence analysis to investigate bacterial and fungal composition (Figure 2).Within each vineyard, two grapevine cultivars, namely 'Summer Black' (five bunches) and 'Kyoho' (five bunches) were collected.Obtained grapes were subjected to yeast cultivation and sequence analysis to investigate bacterial and fungal composition (Figure 2).

Yeast Cultivation
Differing from culturomics, in this study the yeast cultivation experiment was set to isolate yeasts that may act as new starters in the fermentation process.One kilogram of each cultivar of grapes from different growing districts and viticulture conditions were crushed into grape must and cultivated at 20 °C at 155 rpm for four consecutive days.Simultaneously, approximately three grams of grape skin were randomly extracted from five bunches of undamaged grapes.Extracted grape skins were then soaked in 80 mL yeast maintenance media (YMM, per liter composed of 3 g yeast extract, 3 g malt extract, 5 g peptone, and 10 g glucose, Sangon Biotech Co., Shanghai, China, pH 5.5) and cultivated at 30 °C for two consecutive days.Afterward, ten µL of each of the above cultures were evenly spread on YMM agar plates and cultivated overnight at 30 °C.Obtained colonies were selected based on colony shapes, colors, and surface features, as well as their morphologies, under light microscopy.Each colony was purified until the colony's morphology was consistent, and preserved at −80 °C with 30% glycerol.

DNA Extraction, Amplification, and Sequencing
Approximately ten grams of grape skin was randomly extracted from five bunches of undamaged grapes.Grape skins were washed with 20 mL PBS (Sangon Biotech, Shanghai, China) at 80 rpm/min for 30 min.The washing solution was then transferred into a new tube and centrifuged at 4 °C 10,000× g for 15 min.The obtained pellet was re-suspended with 1 mL PBS and centrifugated again at 4 °C 10,000× g for 5 min.The obtained pellet was stored at −20 °C for bacterial and fungal compositional analysis.
Total DNA was extracted using the Mabio DNA isolation kit (Mabio, Guangzhou, China) following the manufacturer's instructions.Obtained microbial DNA was used to study bacterial and fungal compositions.Specifically, for bacterial compositional analysis,

Yeast Cultivation
Differing from culturomics, in this study the yeast cultivation experiment was set to isolate yeasts that may act as new starters in the fermentation process.One kilogram of each cultivar of grapes from different growing districts and viticulture conditions were crushed into grape must and cultivated at 20 • C at 155 rpm for four consecutive days.Simultaneously, approximately three grams of grape skin were randomly extracted from five bunches of undamaged grapes.Extracted grape skins were then soaked in 80 mL yeast maintenance media (YMM, per liter composed of 3 g yeast extract, 3 g malt extract, 5 g peptone, and 10 g glucose, Sangon Biotech Co., Shanghai, China, pH 5.5) and cultivated at 30 • C for two consecutive days.Afterward, ten µL of each of the above cultures were evenly spread on YMM agar plates and cultivated overnight at 30 • C. Obtained colonies were selected based on colony shapes, colors, and surface features, as well as their morphologies, under light microscopy.Each colony was purified until the colony's morphology was consistent, and preserved at −80 • C with 30% glycerol.

DNA Extraction, Amplification, and Sequencing
Approximately ten grams of grape skin was randomly extracted from five bunches of undamaged grapes.Grape skins were washed with 20 mL PBS (Sangon Biotech, Shanghai, China) at 80 rpm/min for 30 min.The washing solution was then transferred into a new tube and centrifuged at 4 • C 10,000× g for 15 min.The obtained pellet was re-suspended with 1 mL PBS and centrifugated again at 4 • C 10,000× g for 5 min.The obtained pellet was stored at −20 • C for bacterial and fungal compositional analysis.

Statistical Analysis
The 16S rRNA gene sequence read counts and ITS read counts were normalized to relative abundance.Microbial diversity indices (Fisher alpha and observed species) were calculated based on ASVs.We used the Shapiro-Wilk test to test data normality, and used the Wilcoxon signed-rank test for comparisons of not-normally distributed data.Variation partitioning was assessed by fitting environmental factors (i.e., sampling district, viticultural conditions, and grapevine cultivars) to the composition of microbial populations using the RDA (Redundancy Analysis, RDA) function in combination with the Monte Carlo permutation in the vegan package.Obtained p values < 0.05 were considered to indicate significant differences.All statistics were conducted in R (R 4-4.2.3).

General Characteristics of Table Grape Microbiota
In total, we identified 252 fungi and 594 bacteria genera from collected grapes.The fungal population was on average dominated by Cladosporium (57.5%), followed by Aspergillus (16.6%),Alternaria (10.9%), and Zygosporium (4.9%).On average, Cladosporium accounted for 70.70% of the total fungal community of Kyoho, and 45.3% of that of Summer Black.Percentages of Cladosporium were 72.9% and 66.8% of the total fungal communities, respectively, of grapes collected from Jinshan and Pudong.The percentage of Cladosporium was only 41.6% in grapes collected from Jiading.
The bacterial population was on average dominated by 43.6% Enterobacteriaceae and 22.8% Rhodocyclaceae (Figure S1).At the genus level, unknown bacterial genera accounted for 32.7% of sequences; Thauera (22.8%) and Pantoea (17.3%) were the other dominant sequence types.Of note, Thauera and Pantoea accounted for more than 50% of the total identified bacteria on average.Nevertheless, this dominance was only observed in grapes collected from Jinshan and Jiading.Specifically, Thauera and Pantoea collectively accounted for 47.12% and 47.73% of the total bacterial communities, respectively, of grapes collected from Jiading and Jinshan.The percentage of Thauera and Pantoea was only 26.94% of the total bacterial community of grapes collected from Pudong, which was instead dominated by unknown bacterial genera.Nevertheless, despite differences in sampling districts, the Kyoho bacterial community was composed of 26.66% Thauera and 5.75% Pantoea, while Summer Black's was composed of 17.65% and 32.77%, respectively.

Potential Contributors to Grape Skin Mycobiota and Microbiota
The RDA investigation of potential contributions of differences in grapevine cultivars ('Summer Black' and 'Kyoho'), viticultural conditions (i.e., grapes grown in open air or in a greenhouse), and sampling locations (three different districts in Shanghai) to the variation of grape-skin microbiota and mycobiota found only a significant contribution from sampling locations, which explained 16.94% of the total mycobiota variation.
Permutational multivariate analysis of variance (PERMANOVA) revealed a significant difference only between the mycobiota of Jiading grapes and those Pudong grapes based on the Bray-Curtis distance matrix (considering the relative abundance of ASVs) (Figure 3A), but no significant difference based on the Jaccard distance matrix (considering only presence or absence of ASVs) (Figure 3B).Regarding specific fungal taxa, their mycobiota were all dominated by Cladosporium, which was highest in Pudong grapes and lowest in Jiading grapes (Figure 3A).In addition, the relative abundance of Aspergillus was higher in Jiading and Jinshan grapes compared to that of Pudong grapes.Moreover, the relative abundance of Alternaria was highest in Jiading grapes and lowest in Pudong grapes.Nevertheless, no significant difference was identified in the relative abundance of specific fungal genera between different districts.
for 47.12% and 47.73% of the total bacterial communities, respectively, of grapes collected from Jiading and Jinshan.The percentage of Thauera and Pantoea was only 26.94% of the total bacterial community of grapes collected from Pudong, which was instead dominated by unknown bacterial genera.Nevertheless, despite differences in sampling districts, the Kyoho bacterial community was composed of 26.66% Thauera and 5.75% Pantoea, while Summer Black's was composed of 17.65% and 32.77%, respectively.

Potential Contributors to Grape Skin Mycobiota and Microbiota
The RDA investigation of potential contributions of differences in grapevine cultivars ('Summer Black' and 'Kyoho'), viticultural conditions (i.e., grapes grown in open air or in a greenhouse), and sampling locations (three different districts in Shanghai) to the variation of grape-skin microbiota and mycobiota found only a significant contribution from sampling locations, which explained 16.94% of the total mycobiota variation.
Permutational multivariate analysis of variance (PERMANOVA) revealed a significant difference only between the mycobiota of Jiading grapes and those of Pudong grapes based on the Bray-Curtis distance matrix (considering the relative abundance of ASVs) (Figure 3A), but no significant difference based on the Jaccard distance matrix (considering only presence or absence of ASVs) (Figure 3B).Regarding specific fungal taxa, their mycobiota were all dominated by Cladosporium, which was highest in Pudong grapes and lowest in Jiading grapes (Figure 3A).In addition, the relative abundance of Aspergillus was higher in Jiading and Jinshan grapes compared to that of Pudong grapes.Moreover, the relative abundance of Alternaria was highest in Jiading grapes and lowest in Pudong grapes.Nevertheless, no significant difference was identified in the relative abundance of specific fungal genera between different districts.

Mycobiota and Microbiota of Different Table Grapes
Although no significant contribution from grapevine cultivars ('Summer Black' and 'Kyoho') was identified based on the RDA, it was generally believed that different grapevine cultivars enharbour distinct bacterial and fungal communities.As such, we further investigated the impact of grapevine cultivars on table grape microbiota and mycobiota.
PERMANOVA based on Bray-Curtis and Jaccard distance matrices revealed no significant difference between mycobiotas (Figure 4A,B) of 'Summer Black' and 'Kyoho' grapes.Nor was a significant difference identified in the fungal richness (observed ASVs) or diversity (inverse Simpson) between these two grape cultivars (Figure 4C,D).Similarly, no significant differences were detected in the bacterial profile (Figure 4E,F), richness, or diversity (Figure 4G,H) between these two grape varieties.

Mycobiota and Microbiota of Different Table Grapes
Although no significant contribution from grapevine cultivars ('Summer Black' and 'Kyoho') was identified based on the RDA, it was generally believed that different grapevine cultivars enharbour distinct bacterial and fungal communities.As such, we further investigated the impact of grapevine cultivars on table grape microbiota and mycobiota.
PERMANOVA based on Bray-Curtis and Jaccard distance matrices revealed no significant difference between mycobiotas (Figure 4A,B) of 'Summer Black' and 'Kyoho' grapes.Nor was a significant difference identified in the fungal richness (observed ASVs) or diversity (inverse Simpson) between these two grape cultivars (Figure 4C,D).Similarly, no significant differences were detected in the bacterial profile (Figure 4E,F), richness, or diversity (Figure 4G,H) between these two grape varieties.Regarding microbial functionality, predicted bacterial functionality profiles did not differ significantly between 'Kyoho' and 'Summer Black' grapes (Figure 5A), and redundancy analysis found no significant contribution from sampling districts and viticulture conditions, while studied grape skin microbiota enharboured pathways coding for aerobic respiration, fatty acid oxidation, and pathways for biosynthesis (Figure 5B).Regarding microbial functionality, predicted bacterial functionality profiles did not differ significantly between 'Kyoho' and 'Summer Black' grapes (Figure 5A), and redundancy analysis found no significant contribution from sampling districts and viticulture conditions, while studied grape skin microbiota enharboured pathways coding for aerobic respiration, fatty acid oxidation, and pathways for biosynthesis (Figure 5B).

Discussion
In the current study, we investigated the bacterial and fungal communities of two table grape types and explored the impact of potential contributors (grape-growing districts, viticultural conditions, and grapevine cultivars) thereof.Our hypothesis was that grapevine cultivars, grape-sampling districts, and viticultural conditions could contribute to variation in the bacterial and fungal communities of table grapes.We found that the microbiota of both grapevine cultivars was dominated by Cladosporium, Alternaria, Aspergillus, Thauera, and Pantoea.Of note, the impact of grapevine cultivars and viticultural conditions on grape skin microbiota and mycobiota variation was not detectable, while grapes' growing locations/districts significantly contributed to the variation thereof.Moreover, the enrichment and isolation of yeasts successfully identified Hanseniaspora opuntiae, Pichia terricola, Rhodotorula mucilaginosa, Candida stellimalicola, and Kodamaea ohmeri from the skins.Some of these are co-fermenters in fermentation, while others can act as spoilage yeast or even pathogens.Collectively, these findings document the potential contribution of grapegrowing environments to table grape skin microbiota and mycobiota, and highlight the importance of microbiota management in the production and consumption of these grapes.
In contrast to the thought that different grapevine cultivars select for distinct bacterial and fungal communities [7], in this study, cultivar differences did not clearly influence grape skin microbiota and mycobiota of 'Summer Black' and 'Kyoho' grapes.Grape-growing environments, i.e., sampling districts in Shanghai, affected their mycobiota, highlighting a regional effect on grape skin microbiota, as suggested for wine-making grapes [20].These differences in grape skin mycobiota among districts of Shanghai could, in part, be attribute to the environment or air fungal composition, which might be linked with existing biogeographical patterns and/or the presence of nearby industry sites.
Regarding specifical bacterial and fungal taxa, in line with earlier studies, this study found that grape skin mycobiota were dominated by Cladosporium and Alternaria [21].Nevertheless, in addition to these two genera, this study also identified that a large portion of table grapes' skin mycobiota was Aspergillus (16.6%).Aspergillus is present on mature berries and could be considered grape invaders that produce ochratoxin, i.e., a mycotoxin that might be carcinogenic to humans [22].In terms of grape skin microbiota, this study found that Thauera, Pantoea, and unknown bacterial genera dominate the skin.This finding, however, differed from that of wine grapes, which were dominated by Lactobacillus, Streptococcus, and Tatumella [23].This difference could therefore assist in the differentiation of table grape skin microbiota from those of wine grapes.Moreover, the presence and dominance of the identified Thauera and Pantoea could possibly affect the flavor of the grapes and their shelf-life post harvesting.
Complementary to high throughput sequence analysis, in the current study, the cultivation-based method successfully isolated a few rare yeasts.In addition to Hanseniaspora, Pichia, and Candida, which are commonly present on the surface of grapes [8], the current study also successfully cultivated and identified Rhodotorula mucilaginosa [24] and Kodamaea ohmeri [25].Notably, Rhodotorula mucilaginosa and Kodamaea might potentially be pathogenic to humans and lead to safety concerns regarding the studied grapes [26].Specifically, Rhodotorula mucilaginosa [24] could be an opportunistic infectious fungi leading to oral ulcers in AIDs patients [27], while Kodamaea ohmeri [25,26], also known as Pichia ohmeri or Yamadazyma ohmeri, could be a life-threatening pathogen in humans, especially immunocompromised patients [28,29].Furthermore, future studies are warranted to investigate characteristics of these isolated yeasts and uncover their role in grape production, as well as their potential application in food/beverage production.In addition to fungi identified after cultivation, the sequencing-based approach identified some plant pathogens, like Alternaria and Botrytis.The presence of these fungi could also have an impact on Vitis vinifera cultivars.Nevertheless, data on Vitis vinifera cultivars were not noted in the current study.Studies investigating the impact of these fungi on grapevines are still warranted.For instance, the presence of fungus Botrytis, as identified in the current study, could have an impact on grape flavors.Future studies could consider and evaluate the flavor profiles of grapes and co-existing bacterial and fungal communities.
In conclusion, instead of the type of grapevine cultivars ('Summer Black' and 'Kyoho'), the grape-growing environment, i.e., sampling districts in Shanghai, affected grapes' mycobiota significantly, while their mycobiota were dominated by Cladosporium, Alternaria, and Aspergillus.In terms of their microbiota, all were dominated by Thauera, Pantoea, and unknown bacterial genera.Of note, Hanseniaspora opuntiae, Pichia terricola, Rhodotorula mucilaginosa, Candida stellimalicola, and Kodamaea ohmeri were isolated and identified from grape skins.Collectively, the current study demonstrated the impact of the cultivation environment, i.e., sampling location, on grape skin mycobiota, as well as the increased need to monitor/manage table grape skin microbiota.This paper contributes to the growing body of knowledge regarding grape skin microbiota and its potential applications in grape production, food safety, and human health.

Figure 1 .
Figure 1.Sampling locations in Shanghai.Purple dots represented specific sampling locations in each district of Shanghai.

Figure 1 .
Figure 1.Sampling locations in Shanghai.Purple dots represented specific sampling locations in each district of Shanghai.

Figure 3 .
Figure 3. Impact of grape-growing district on overall mycobiota profile.PCoA of mycobiota based on (A) Bray-Curtis or (B) Jaccard distance matrices.(C) Relative abundance of most predominant fungal genus (top 20, ranked based on average relative abundance across entire dataset).Top 20 microbial genera are listed in the legend.Other genera are summarized as "Others".Each column represents a given type of sample from one district.

Figure 3 .
Figure 3. Impact of grape-growing district on overall mycobiota profile.PCoA of mycobiota based on (A) Bray-Curtis or (B) Jaccard distance matrices.(C) Relative abundance of most predominant fungal genus (top 20, ranked based on average relative abundance across entire dataset).Top 20 microbial genera are listed in the legend.Other genera are summarized as "Others".Each column represents a given type of sample from one district.

Figure 4 .
Figure 4. Impact of grape cultivars on mycobiota and microbiota of grape skins.PCoA of mycobiota based on (A) Bray-Curtis or (B) Jaccard distance matrices.(C) Mycobiota richness and (D) diversity.PCoA of microbiota based on (E) Bray-Curtis and (F) Jaccard distance matrices.(G) Microbial richness and (H) diversity.Significant differences between samples based on Bray-Curtis and Jaccard distance matrices (considering only the presence or absence of ASVs) at ASV level were tested by PERMANOVA.The percentage of variation in microbial composition explained by the two principal coordinates is shown at the axes.Purple triangles represent Kyoho grapes, and green circles represent 'Summer Black' grapes.Pairwise comparisons between different sample types were evaluated by the Wilcoxon signed-rank test.PCoA: principal coordinate analysis.InvSimpson: inverse Simpson.PERMANOVA: permutational multivariate analysis of variance.

Figure 4 .
Figure 4. Impact of grape cultivars on mycobiota and microbiota of grape skins.PCoA of mycobiota based on (A) Bray-Curtis or (B) Jaccard distance matrices.(C) Mycobiota richness and (D) diversity.PCoA of microbiota based on (E) Bray-Curtis and (F) Jaccard distance matrices.(G) Microbial richness and (H) diversity.Significant differences between samples based on Bray-Curtis and Jaccard distance matrices (considering only the presence or absence of ASVs) at ASV level were tested by PERMANOVA.The percentage of variation in microbial composition explained by the two principal coordinates is shown at the axes.Purple triangles represent Kyoho grapes, and green circles represent 'Summer Black' grapes.Pairwise comparisons between different sample types were evaluated by the Wilcoxon signed-rank test.PCoA: principal coordinate analysis.InvSimpson: inverse Simpson.PERMANOVA: permutational multivariate analysis of variance.

Horticulturae 2024 , 10 Figure 5 .
Figure 5. (A) PCoA of microbial functional pathways.(B) Relative abundance of most predominant bacterial functional pathways (top 10, ranked base on the averaged relative abundance across the entire dataset).Top 10 microbial functional pathways are listed on the left.

Figure 5 .
Figure 5. (A) PCoA of microbial functional pathways.(B) Relative abundance of most predominant bacterial functional pathways (top 10, ranked base on the averaged relative abundance across the entire dataset).Top 10 microbial functional pathways are listed on the left.

Figure 5 .
Figure 5. (A) PCoA of microbial functional pathways.(B) Relative abundance of most predominant bacterial functional pathways (top 10, ranked base on the averaged relative abundance across the entire dataset).Top 10 microbial functional pathways are listed on the left.

Figure 6 .
Figure 6.Phylogenetic tree constructed based on the sequence alignment of ITS1 and ITS2 regions.The phylogenetic tree was inferred from the "neighbour-joining method", and numbers on the branches represent percentages of supports from the bootstrap replication.The length of the branch represents bootstrap values; values smaller than 70 were collapsed.