ORIGINAL RESEARCH article

Front. Microbiol., 16 February 2023

Sec. Microbe and Virus Interactions with Plants

Volume 14 - 2023 | https://doi.org/10.3389/fmicb.2023.1043430

Three new yeast species from flowers of Camellia sinensis var. assamica collected in Northern Thailand and their tannin tolerance characterization

  • 1. Division of Biochemistry and Biochemical Innovation, Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand

  • 2. Research Center for Multidisciplinary Approaches to Miang, Chiang Mai University, Chiang Mai, Thailand

  • 3. Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand

  • 4. Division of Food Science and Technology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand

  • 5. Division of Microbiology, Department of Biology, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand

  • 6. Research Center of Microbial Diversity and Sustainable Utilization, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand

  • 7. Department of Microbiology, Faculty of Science, Kasetsart University, Bangkok, Thailand

  • 8. Biodiversity Center Kasetsart University (BDCKU), Bangkok, Thailand

  • 9. Global Institute of Food Security and International Agriculture (GIFSIA), Department of Plant Sciences, North Dakota State University, Fargo, ND, United States

Abstract

Our recent research study focused on Miang fermentation revealed that tannin-tolerant yeasts and bacteria play vital roles in the Miang production process. A high proportion of yeast species are associated with plants, insects, or both, and nectar is one of the unexplored sources of yeast biodiversity. Therefore, this study aimed to isolate and identify yeasts of tea flowers of Camellia sinensis var. assamica and to investigate their tannin tolerance, which is a property essential to Miang production processes. A total of 82 yeasts were recovered from a total of 53 flower samples in Northern Thailand. It was found that two and eight yeast strains were distinct from all other known species within the genera Metschnikowia and Wickerhamiella, respectively. These yeast strains were described as three new species, namely, Metschnikowia lannaensis, Wickerhamiella camelliae, and W. thailandensis. The identification of these species was based on phenotypic (morphological, biochemical, and physiological characteristics) and phylogenetic analyses of a combination of the internal transcribed spacer (ITS) regions and the D1/D2 domains of the large subunit (LSU) ribosomal RNA gene. The yeast diversity in tea flowers acquired from Chiang Mai, Lampang, and Nan provinces had a positive correlation with those acquired from Phayao, Chiang Rai, and Phrae, respectively. Wickerhamiella azyma, Candida leandrae, and W. thailandensis were the species uniquely found in tea flowers collected from Nan and Phrae, Chiang Mai, and Lampang provinces, respectively. Some of the tannin-tolerant and/or tannase-producing yeasts were associated with yeasts in the commercial Miang process and those found during Miang production, i.e., C. tropicalis, Hyphopichia burtonii, Meyerozyma caribbica, Pichia manshurica, C. orthopsilosis, Cyberlindnera fabianii, Hanseniaspora uvarum, and Wickerhamomyces anomalus. In conclusion, these studies suggest that floral nectar could support the formation of yeast communities that are beneficial for Miang production.

1. Introduction

Camellia sinensis var. assamica is one of the important materials for processing various tea products, specifically oolong, black, and Pu-erh (Wang et al., 2013). It is suggested by some researchers that C. sinensis var. assamica is the variant of C. sinensis var. sinensis, another major variety grown for tea production to date (Zhang et al., 2019). Miang is a fermented food product traditionally produced from the leaves of C. sinensis var. assamica in the mountainous areas of Northern Thailand. It is different from other tea products, and it is known as eating tea or chewing tea. Miang also has long deep-rooted social and cultural integration and relevance with the people of Northern Thailand (Khanongnuch et al., 2017). In our previous studies, several synergistic interactions were observed among bacteria, yeasts, and filamentous fungi in Miang, which was associated with the release of several bioactive compounds particularly gallated catechins and gallated epicatechin (Unban et al., 2019, 2020a; Kodchasee et al., 2021). Most associated microorganisms, especially bacteria and yeasts, exhibited specific characteristics in terms of tannin tolerance and tannase production (Kanpiengjai et al., 2016; Chaikaew et al., 2017). Production of lignocellulose-degrading enzymes and tannase by tannin-tolerant microorganisms is believed to be essential for the release of useful therapeutic phenolic compounds, specifically catechins, tannins, and gallic acid (Unban et al., 2020b). To identify the deep-rooted microbial environment during tea leaf fermentation, bacterial and fungal diversities were determined (Unban et al., 2020a; Kodchasee et al., 2021). Based on the previous results, it was clearly demonstrated that tannin-tolerant ability of microorganisms is aligned with their taxonomy, specifically at family, genus, and species levels and is one of the key criteria for evaluating the evolution of microorganisms associated with Miang fermentation (Kanpiengjai et al., 2016; Chaikaew et al., 2017; Unban et al., 2020b). Plants offer a great variety of favorable microenvironments specifically for yeast growth. Unique yeast communities have been associated with the rhizosphere, the phylloplane, plant exudates, plant necrotic tissues, fruits, and flowers (De Vega et al., 2017). In a previous study, several yeasts, including Aureobasidium melanogenum, Kazachstania aquatica, Saturnispora diversa, Saturnispora sekii, Schwanniomyces pseudopolymorphus, Wickerhamomyces anomalus, Apiotrichum scarabaeorum, Curvibasidium pallidicorallinum, Papiliotrema laurentii, Rhodosporidiobolus ruineniae, Trichosporon asahii, and T. coremiiforme, were discovered from the soil of Assam tea (Kumla et al., 2020), but their taxonomic alignment is not associated with yeasts found in Miang. The current study has revealed that nectar-inhabiting microorganisms are likely to interact with plants and pollinators. Yeasts, commonly present in nectar, can modify a variety of nectar traits that are important for pollinator attraction (Schaeffer et al., 2017). Therefore, the yeast community in Assam tea flowers may act as a reservoir of yeasts that are associated with the Miang production process as well as unique Miang products from each production area. Therefore, this study aimed to investigate the yeast ecology of Assam tea flowers collected from six provinces in upper Northern Thailand by using the culture-dependent method and to characterize the tannin tolerance capability of the isolated yeasts. The study insights can explain the association between flower yeasts and Miang fermentation. To our knowledge, this is the first report to describe the isolation and molecular-based identification of yeasts associated with tea flowers.

2. Materials and methods

2.1. Chemical and media

Tannin (gallotannin) was purchased from LOBA Chemie (Mumbai, India), and methyl gallate was purchased from Sigma (Steinheim, Germany). Ingredients such as peptone, yeast extract, malt extract, and carbon sources such as glucose, galactose, sucrose, maltose, lactose, raffinose, and trehalose were purchased from HiMedia (Nashik, India). The genomic DNA of yeast was extracted using a Wizard® Genomic DNA purification kit (Promega Corp., Madison, WI, United States). Nucleic acid amplifications were performed using TaKaRa Ex Taq® (Shiga, Japan). Other reagents and solvents were of analytical grade.

2.2. Sampling site and sample collection

Assam tea flower samples were directly collected from 53 sites that were located within 19 tea plantations used for Miang production from December 2015 to May 2016. These locations covered all regions of the Miang production area, which include six provinces of upper Northern Thailand, namely, Chiang Mai, Chiang Rai, Phrae, Nan, Lampang, and Phayao (Figure 1), where they were recorded to be the Miang production area in previous studies (Khanongnuch et al., 2017). The main tea plantations for Miang production are in Chiang Mai (28 sites in eight locations), Phrae (10 sites in three locations), Nan (seven sites in three locations), Chiang Rai (four sites in two locations), Phayao (two sites in two locations), and Lampang (one location). Each tea flower sample was aseptically collected from each sampling site, put in a sterile plastic bag, date recorded, stored on ice, and delivered to the laboratory within 24 h prior to the sample preparation and further yeast isolation. The sampling site, location, and date of the collection were recorded and given code numbers (Supplementary Table S1).

Figure 1

2.3. Isolation and purification of yeasts

Approximately 1 g of the tea flower sample was homogenized in 10 mL of sterile 0.85% NaCl as a diluent using a Masticator® homogenizer blender (Barcelona, Spain) for 5 min. After that, the homogenate was serially diluted to 100–105 with the diluent, and 0.1 mL of each dilution was plated on yeast–malt extract agar (YMA) (1% yeast extract, 2% peptone, 2% glucose, and 1.5% agar) supplemented with 100 mg/L of chloramphenicol and then incubated at 30°C for 48 h. Yeast colonies with distinct morphologies were isolated and streaked on the same medium without the supplementation of chloramphenicol for purification (Khunnamwong et al., 2022). For storage, a single colony of each yeast isolate was inoculated in yeast–malt extract broth (YMB) (1% yeast extract, 2% peptone, and 2% glucose) and incubated at 30°C on a 150-rpm rotary shaker for 24 h. The culture was mixed with 30% (v/v) glycerol to make a final concentration of 15% before being stored at −80°C.

2.4. Yeast identification

2.4.1. Molecular identification

Amplification of the D1/D2 domain of the large subunit (LSU) rRNA gene or the internal transcribed spacer (ITS) region was performed by a standard polymerase chain reaction (PCR) using genomic DNA of each yeast as a template with forward primer NL1 (5′-GCATATCAATAAGCGGGGAAAAG-3′) and reverse primer NL4 (5′-GGTCCGTGTTTCAAGACGG-3′), and forward primer ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′) and reverse primer ITS4 (5′-TCCTCCGCTTATTGATATGC-3′), respectively. Amplicons obtained from PCR were purified and sequenced by a sequencing service provider (1st BASE Ptd Ltd., Singapore).

Analysis of the ITS and D1/D2 sequences was conducted by similarity searches using the BLAST program available at NCBI.1 The yeast strains showing nucleotide substitutions of greater than 1% in the D1/D2 domain of the LSU rRNA gene are usually considered different species (Kurtzman and Robnett, 1997) and were selected for further phylogenetic analysis and phenotypic characterization. The holotypes of new species were preserved in a metabolically inactive state in the culture collection of the Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand. The ex-type cultures were deposited in a metabolically inactive state at the culture collection of the Sustainable Development of Biological Resources, Faculty of Science, Chiang Mai University (SDBR-CMU), Chiang Mai, Thailand.

2.4.2. Phylogenetic analysis of new yeast species

The sequences from this study and those obtained from previous studies together with sequences downloaded from the nucleotide GenBank database were selected and provided for phylogenetic analysis. Multiple sequence alignment was performed with MUSCLE (Edgar, 2004) and manually checked using BioEdit v. 6.0.7 (Hall, 1999). Phylogenetic analysis was performed based on the combined dataset of ITS and D1/D2 sequences. A phylogenetic tree was constructed using maximum likelihood (ML) and Bayesian inference (BI) methods. The best substitution models for ML and BI analyses were estimated by the Akaike information criterion (AIC) in jModelTest 2.1.10 (Darriba et al., 2012). ML analysis was carried out on RAxML v7.0.3 under the GTRCAT model with 25 categories and 1,000 bootstrap (BS) replications (Felsenstein, 1985; Stamatakis, 2006) via the online portal CIPRES Science Gateway v. 3.3 (Miller et al., 2010). BI analysis was performed by MrBayes v3.2.6 (Ronquist and Huelsenbeck, 2003). For the BI analysis, six simultaneous Markov chains were run for one million generations with random initial trees, wherein every 1,000 generations were sampled. A burn-in phase was employed to discard the first 2,000 of the trees, while the remaining trees were used to construct the 50% majority-rule consensus phylogram with calculated Bayesian posterior probabilities (PPs). Tree topologies were visualized in FigTree v1.4.0 (Rambaut, 2016).

2.4.3. Morphology characterization of new yeast species

The morphological characteristics of expected new yeast species were determined according to the standard methods described by Kurtzman et al. (2011), De Vega et al. (2017), and Shibayama et al. (2020). Colony formation was observed on YMA after incubation at 25°C for 3 days. The formation of pseudohyphae and true hyphae was investigated by cultivation on potato dextrose agar (PDA) (20% potato infusion, 2% glucose, and 1.5% agar) in slide culture at 25°C for up to 4 weeks. Ascospore formation was investigated for individual strains and strain pairs on PDA, corn meal agar, 5% malt extract agar (5% malt extract and 1.5% agar), yeast carbon base agar supplemented with 0.01% yeast extract (YCBY), Fowell’s acetate agar, yeast extract-peptone dextrose (YPD) agar (1% yeast extract, 2% peptone, 2% glucose, and 1.5% agar), glucose yeast (GY) agar (1% glucose, 1% yeast extract, 2% agar), and YMA at 25°C for 4 weeks. Macromorphological characteristics were examined under a light microscope (Nikon Eclipse Ni U, Tokyo, Japan). Size data of the anatomical structure such as cells, pseudohyphae, asci, and ascospores were based on at least 50 measurements of each structure using the Tarosoft® Image Frame Work program.

2.4.4. Biochemical and physiological characterization of new yeast species

All strains of the new yeast species were characterized biochemically and physiologically, according to the standard methods described by Kurtzman et al. (2011). Fermentation of carbohydrates was carried out in a liquid medium using Durham fermentation tubes. Carbon source and nitrogen source assimilation tests were carried out in a liquid medium, and starved inocula were used in the nitrogen assimilation test. Cycloheximide resistance was also performed in a liquid medium, while urea hydrolysis was conducted on an agar slant. Acid production and the diazonium blue B (DBB) reaction were investigated on the solid media in Petri dishes. The effect of temperature on yeast growth was investigated using YMA at various temperatures ranging from 15 to 37°C. The ability to grow in media of high osmotic pressure was performed using 50 and 60% glucose agar, and 10 and 16% sodium chloride plus 5% glucose medium.

2.5. Determination of yeast diversity isolated from tea flowers

In the yeast isolation step, each tea flower sample was collected from each sampling site. The number of isolated yeasts from each sample was recorded, calculated, and further expressed in terms of the number of yeasts per sample. After the yeast identification step, the number of samples found for each species was recorded. The frequency of each species was calculated by the number of samples that contained certain species divided by the total number of tea flower samples. Principal component analyses (PCA) based on the relative abundance of yeast species in tea flower samples and provinces were performed using GraphPad Prism9 software version 9.4.1 (GraphPad Software, San Diego, CA, United States).

2.6. Tannin-tolerant test and detection of tannase

Tannin-tolerant test was evaluated by inoculating a single yeast colony on the YMA plate with and without the supplementation of 10, 20, 30, 40, and 50 g/L tannic acid according to the method described in a previous study (Kanpiengjai et al., 2016). Growth of yeasts was observed after incubation at 30°C for 24 to 72 h. Principal component analyses (PCA) based on the relative abundance of tannin-tolerant yeast species in tea flower samples and provinces were performed.

Detection of tannase was also performed by a visual reading method (Osawa and Walsh, 1993). Briefly, fresh cultures on YMA containing 10 g/L tannic acid were harvested and suspended in 750 μL of substrate medium containing 33 mM NaH2PO4 and 20 mM methyl gallate for the purposes of preparing a dense suspension. The substrate medium was incubated aerobically at 30°C for 24 h. The mixture was later alkalinized by adding an equal amount of a saturated NaHCO3 solution (pH 8.6) and then incubating it at 25°C for 1 h. This suspension was centrifuged to remove any cells, and the absorbance of the supernatant was measured at 440 nm. The supernatant changed color from green to brown. An absorbance value of higher than 0.5 was considered a positive indication of tannase activity.

3. Results

3.1. Yeast isolation and identification

Based on colony and shape, a total of 82 yeast isolates were generated from 47 out of 53 tea flower samples. The number of yeast strains was dependent on the number of sampling sites derived from each province; thus, Chiang Mai was the major source of the isolated yeasts. The number of yeast isolates per sample was between 0 and 6 isolates per sample (Supplementary Table S1). Two samples possessed the highest number of yeast isolates with 4–6 isolates per sample. They included tea flowers from Na Rai Luang subdistrict, Nan (FLA31), and Longkhod subdistrict, Chiang Mai (FLA11), while others were recorded to have the ratio of yeast species and site in a range of 0 to 3 isolates per site.

Initially, all yeasts were identified using sequence analysis of the D1/D2 domain of the LSU rRNA gene. The sequence of the D1/D2 domain was aligned with the related type strains deposited in NCBI. For each species, the sequence of the D1/D2 domain that shared more than 98% with those of the closest species is shown in Table 1. In total, 72 yeast strains were identified within 39 known yeast species, while another 10 yeast strains were representative of three undescribed species in the Wickerhamiella (strains FLA5.4, FLA13.1, FLA23.2, FLA46.4, FLA4.2, FLA24.1, FLA27.3, and FLA30.2) and Metschnikowia (FLA8.2 and FLA27.1) clade (Supplementary S2). In addition, each strain of the three undescribed species was isolated from a different location in the tea plantation area.

Table 1

IsolateAccession No.CBS type speciesIdentification/Similarity (%)LocationSites (subdistrict/district)ProvinceNo. of sample found for each species% Frequency of occurrence
FLA4.1MW602315Candida leandrae CBS 9735C. leandrae/99.8019°07′00.3″N 98°43′35.0″ESite 3: Papae/ Mae TaengChiang Mai910.6
FLA5.1MW602316Candida leandrae CBS 9735C. leandrae/99.8019°44′01.9″N 98°54′04.1″ESite 1: Mae Na/ Chiang DaoChiang Mai
FLA7.1MW602319Candida leandrae CBS 9735C. leandrae/99.8018°50′29.3″N 98°54′15.9″ESite 2: Suthep/ MuangChiang Mai
FLA8.1MW602320Candida leandrae CBS 9735C. leandrae/99.8018°50′29.3″N 98°54′15.9″ESite 3: Suthep/ MuangChiang Mai
FLA9.1MW602322Candida leandrae CBS 9735C. leandrae/99.8018°50′29.3″N 98°54′15.9″ESite 4: Suthep/ MuangChiang Mai
FLA15.1MW602331Candida leandrae CBS 9735C. leandrae/99.8018°57′41.2″N 99°14′06.9″ESite 1: Pa Miang/ Doi SaketChiang Mai
FLA16.1MW602332Candida leandrae CBS 9735C. leandrae/99.8018°57′46.2″N 99°21′22.3″ESite 1: Thepsadet/ Doi SaketChiang Mai
FLA19.1MW602335Candida leandrae CBS 9735C. leandrae/99.8018°57′46.2″N 99°21′22.3″ESite 3: Thepsadet/ Doi SaketChiang Mai
FLA28.2MW602340Candida leandrae CBS 9735C. leandrae/99.8019°35′19.1″N 100°22′30.0″ESite 1: Phu Sang/ Chiang KhamPhayao
FLA44.2MW602372Candida orthopsilosis CBS 10744C. orthopsilosis/100.0017°59′59.2″N 100°15′19.7″ESite 1: Cho Hae/ MuangPhrae22.4
FLA46.2MW602376Candida orthopsilosis CBS 10744C. orthopsilosis/100.0018°07′33.4″N 100°18′57.5″ESite 2: Suan Khuen/ MuangPhrae
FLA40.1MW602364Clavispora pseudohaemulonis CBS 10004Cl. pseudohaemulonis/99.2118°07′02.0″N 100°17′50.9″ESite 3: Pa Daeng/ MuangPhrae11.2
FLA39.2MW602345Cyberlindnera fabianii CBS 5640Cy. fabianii/100.0018°07′02.0″N 100°17′50.9″ESite 2: Pa Daeng/ MuangPhrae11.2
FLA32.2MW602349Debaryomyces fabryi CBS 789D. fabryi/100.0019°14′41.0″N 100°59′52.0″ESite 1: Sakad/ PuaNan11.2
FLA31.1MW602342Debaryomyces nepalensis CBS 7761D. natalensis/100.0019°23′06.9″N 100°41′46.7″ESite 1: NaRai Luang/ Song KhwaeNan22.4
FLA46.3MW602377Debaryomyces nepalensis CBS 7761D. nepalensis/100.0018°07′33.4″N 100°18′57.5″ESite 2: Suan Khuen/ MuangPhrae
FLA34MW602352Debaryomyces subglobosus CBS 792D. subglobosus/100.0019°14′41.0″N 100°59′52.0″ESite 3: Sakad/ PuaNan11.2
FLA33MW602351Diutina catenulata CBS 565Di. catenulata/100.0019°14′41.0″N 100°59′52.0″ESite 2: Sakad/ PuaNan11.2
FLA9.2MW602323Hanseniaspora uvarum CBS 104H. uvarum/99.6518°50′29.3″N 98°54′15.9″ESite 4: Suthep/ MuangChiang Mai11.2
FLA31.6MW602346Hyphopichia buttonii CBS 4571Hy. burtonii/100.0019°23′06.9″N 100°41′46.7″ESite 1: NaRai Luang/ Song KhwaeNan11.2
FLA31.4MW602344Kodamaea ohmeri CBS 5367K. ohmeri/100.0019°23′06.9″N 100°41′46.7″ESite 1: NaRai Luang/ Song KhwaeNan22.4
FLA46.1MW602375Kodamaea ohmeri CBS 5367K. ohmeri/100.0018°07′33.4″N 100°18′57.5″ESite 2: Suan Khuen/ MuangPhrae
FLA11.1MW602325Kodamaea restingae CBS 8492K. restingae/98.5819°05′43.3″N 99°10′54.5″ESite 1: Longkhod/ PhraoChiang Mai11.2
FLA11.3MW602327Kurtzmaniella quercitrusa CBS 4412Ku. quercitrusa/100.0019°05′43.3″N 99°10′54.5″ESite 1: Longkhod/ PhraoChiang Mai33.5
FLA12.1MW602330Kurtzmaniella quercitrusa CBS 4412Ku. quercitrusa/99.8219°05′43.3″N 99°10′54.5″ESite 1: Mae Ai/ Mai AiChiang Mai
FLA29.1MW602341Kurtzmaniella quercitrusa CBS 4412Ku. quercitrusa/100.0019°14′18.4″N 100°19′50.1″ESite 1: Oi/ PongPhayao
FLA47.1MW602378Meryerozyma caribbica CBS 9966Me. caribbica/100.0018°55′24.7″N 99°12′03.5″ESite 4: Pa Miang/ Doi SaketChiang Mai11.2
FLA22.1MW602337Metschnikowia chrysomelidarum CBS 9904M. chrysomelidarum/99.5819°05′16.8″N 99°22′29.6″ESite 1: MaeChedi Mai/ Wiang PapaoChiang Rai11.2
FLA6.1MW602318Metschnikowia hawaiiana CBS 9146M. hawaiiana/99.8018°50′29.3″N 98°54′15.9″ESite 1: Suthep/ MuangChiang Mai33.5
FLA32.1MW602348Metschnikowia hawaiiana CBS 9146M. hawaiiana/100.0019°14′41.0″N 100°59′52.0″ESite 1: Sakad/ PuaNan
FLA43.1MW602370Metschnikowia hawaiiana CBS 9146M. hawaiiana/100.0018°07′02.0″N 100°17′50.9″ESite 6: Pa Daeng/ MuangPhrae
FLA1.2MW602309Metschnikowia koreensis CBS 8854M. koreensis/100.0018°47′08.5″N 100°39′37.5″ESite 1: Ruang/ MuangNan22.4
FLA41.1MW602367Metschnikowia koreensis CBS 8854M. koreensis/100.0018°07′02.0″N 100°17′50.9″ESite 4: Pa Daeng/ MuangPhrae
FLA31.2MW602343Metschnikowia laotica CBS 12961M. laotica/99.2019°23′06.9″N 100°41′46.7″ESite 1: Na Rai Luang/ Song KhwaeNan11.2
FLA2.3MW602312Metschnikowia rancensis CBS 8174M. rancensis/100.0019°07′16.8″N 98°42′28.9″ESite 1: Papae/ Mae TaengChiang Mai44.7
FLA18.1MW602334Metschnikowia rancensis CBS 8174M. rancensis/99.8018°57′42.3″N 99°14′04.9″ESite 2: Pa Miang/ Doi SaketChiang Mai
FLA22.2MW602338Metschnikowia rancensis CBS 8174M. rancensis/99.8019°05′16.8″N 99°22′29.6″ESite 1: MaeChedi Mai/ Wiang PapaoChiang Rai
FLA38.1MW602359Metschnikowia rancensis CBS 8174M. rancensis/99.8018°47′08.5″N 100°39′37.5″ESite 3: Ruang/ MuangNan
FLA11.4MW602328Nakazawaea ishiwadae CBS 6022N. ishiwadae/100.0019°05′43.3″N 99°10′54.5″ESite 1: Longkhod/ PhraoChiang Mai11.2
FLA8.2KY640629Metschinikowia hawiiana CBS 9146Metschinikowia sp./89.3518°50′29.3″N 98°54′15.9″ESite 3: Suthep/ MuangChiang Mai22.4
FLA27.1MW542582Metschinikowia hawiiana CBS 9146Metschinikowia sp./89.1618°49′43.4″N 99°23′14.6″ESite 2: Chae Son/ Muang PanLampang
FLA5.4KY411895Wickerhamiella natalensis CBS 14161Wickerhamiella sp./96.9519°44′01.9″N 98°54′04.1″ESite 1: Mae Na/ Chiang DaoChiang Mai44.7
FLA13.1MW542579Wickerhamiella natalensis CBS 14161Wickerhamiella sp./96.5920°05′09.9″N 99°17′11.0″ESite 2: Mae Ai/ Mae AiChiang Mai
FLA23.2MW542580Wickerhamiella natalensis CBS 14161Wickerhamiella sp./96.9519°05′16.8″N 99°22′29.6″ESite 2: MaeChedi Mai/ Wiang PapaoChiang Rai
FLA46.4MW542585Wickerhamiella natalensis CBS 14161Wickerhamiella sp./96.9518°07′33.4″N 100°18′57.5″ESite 2: Suan Khuen/ MuangPhrae
FLA4.2KY411893Wickerhamiella musiphila CBS10697Wickerhamiella sp./94.3819°07′00.3″N 98°43′35.0″ESite 3: Papae/ Mae TaengChiang Mai44.7
FLA24.1MW542581Wickerhamiella musiphila CBS 10697Wickerhamiella sp./94.3818°49′43.4″N 99°23′14.6″ESite 1: Chae Son/ Muang PanLampang
FLA27.3KY640634Wickerhamiella musiphila CBS 10697Wickerhamiella sp./94.5718°49′43.4″N 99°23′14.6″ESite 2: Chae Son/ Muang PanLampang
FLA30.2MW542584Wickerhamiella musiphila CBS 10697Wickerhamiella sp./94.5719°34′00.7″N 100°02′35.8″ESite 1: Mae Loi/ ThoengChiang Rai
FLA39.3MW602363Pichia manshurica CBS 209P. manshurica/100.0018°07′02.0″N 100°17′50.9″ESite 2: Pa Daeng/ MuangPhrae33.5
FLA40.3MW602365Pichia manshurica CBS 209P. manshurica/100.0018°07′02.0″N 100°17′50.9″ESite 3: Pa Daeng/ MuangPhrae
FLA41.3MW602368Pichia manshurica CBS 209P. manshurica/100.0018°07′02.0″N 100°17′50.9″ESite 4: Pa Daeng/ MuangPhrae
FLA48.1MW602739Priceomyces melissophilus CBS 6344Pr. melissophilus/98.8018°55′24.7″N 99°12′03.5″ESite 5: Pa Miang/ Doi SaketChiang Mai11.2
FLA42.1MW602369Saccharomyces cerevisiae CBS 5493S. cerevisiae/100.0018°07′02.0″N 100°17′50.9″ESite 5: Pa Daeng/ MuangPhrae11.2
FLA17.1MW602333Saccharomycopsis fodiens CBS 8332Sc. fodiens/100.0018°55′15.8″N 99°19′49.0″ESite 2: Thepsadet/ Doi SaketChiang Mai11.2
FLA5.3KY411894Suhomyces xylopsoci CBS 6037Su. xylopsoci/98.4019°44′01.9″N 98°54′04.1″ESite 1: Mae Na/ Chiang DaoChiang Mai33.5
FLA44.1MW602371Suhomyces xylopsoci CBS 6037Su. xylopsoci/98.5717°59′59.2″N 100°15′19.7″ESite 1: Cho Hae/ MuangPhrae
FLA45.1MW602373Suhomyces xylopsoci CBS 6037Su. xylopsoci/98.5717°59′59.2″N 100°15′19.7″ESite 2: Cho Hae/ MuangPhrae
FLA1.3MW602310Wickerhamiella azyma CBS 6826W. azyma/99.2818°47′08.5″N 100°39′37.5″ESite 1: Ruang/ MuangNan89.4
FLA8.3MW602321Wickerhamiella azyma CBS 6826W. azyma/99.2818°50′29.3″N 98°54′15.9″ESite 3: Suthep/ MuangChiang Mai
FLA20.1MW602336Wickerhamiella azyma CBS 6826W. azyma/99.2818°55′18.7″N 99°19′54.4″ESite 4: Thepsadet/ Doi SaketChiang Mai
FLA35.2MW602354Wickerhamiella azyma CBS 6826W. azyma/99.2818°07′33.4″N 100°18′57.5″ESite 1: Suan Khuen/ MuangPhrae
FLA36.3MW602357Wickerhamiella azyma CBS 6826W. azyma/99.2818°07′02.0″N 100°17′50.9″ESite 1: Pa Daeng/ MuangPhrae
FLA37.1MW602358Wickerhamiella azyma CBS 6826W. azyma/99.6418°47′08.5″N 100°39′37.5″ESite 2: Ruang/ MuangNan
FLA38.2MW602360Wickerhamiella azyma CBS 6826W. azyma/99.2818°47′08.5″N 100°39′37.5″ESite 3: Ruang/ MuangNan
FLA40.4MW602366Wickerhamiella azyma CBS 6826W. azyma/99.2718°07′02.0″N 100°17′50.9″ESite 3: Pa Daeng/ MuangPhrae
FLA39.1MW602361Wickerhamomyces anomalus CBS 5759Wc. anomalus/100.0018°07′02.0″N 100°17′50.9″ESite 2: Pa Daeng/ MuangPhrae11.2
FLA11.2MW602326Wickerhamomyces ciferrii CBS 111Wc. ciferrii/100.0019°05′43.3″N 99°10′54.5″ESite 1: Longkhod/ PhraoChiang Mai11.2
FLA45.4MW602374Yamadazyma dushanensis CBS 13914Y. dushanensis/98.6817°59′59.2″N 100°15′19.7″ESite 2: Cho Hae/ MuangPhrae11.2
FLA2.4MW602313Hannaella luteola CBS 943Ha. luteola/98.9919°07′16.8″N 98°42′28.9″ESite 1: Papae/ Mae TaengChiang Mai11.2
FLA26.1MW602339Hannaella pagnoccae CBS 11142Ha. pagnoccae/100.0019°11′58.0″N 99°31′00.5″ESite 3: MaeChedi Mai/ Wiang PapaoChiang Rai11.2
FLA2.1MW602311Hannaella sinensis CBS 7225Ha. sinensis/100.0019°07′16.8″N 98°42′28.9″ESite 1: Papae/ Mae TaengChiang Mai11.2
FLA50MW602380Moesziomyces antarcticus CBS 5955Mo. antarcticus/100.0018°57′46.2″N 99°21′22.3″ESite 7: Thepsadet/ Doi SaketChiang Mai22.4
FLA53MW602382Moesziomyces antarcticus CBS 5955Mo. antarcticus/100.0018°47′20.4″N 99°15′09.5″ESite 2: On Nuea/ Mae OnChiang Mai
FLA31.5MW602345Papiliotrema flavoscens CBS 6475Pa. flavescens/100.0019°23′06.9″N 100°41′46.7″ESite 1: Na Rai Luang/ Song KhwaeNan22.4
FLA51.1MW602381Papiliotrema flavoscens CBS 7140Pa. flavoscens/100.0018°57′46.2″N 99°21′22.3″ESite 8: Thepsadet/ Doi SaketChiang Mai
FLA3.1MW602314Pseudozyma hubeiensis CBS 10077Ps. hubeiensis/99.6719°07′00.3″N 98°43′35.0″ESite 2: Papae/ Mae TaengChiang Mai33.5
FLA32.3MW602350Pseudozyma hubeiensis CBS 10077Ps. hubeiensis/100.0019°14′41.0″N 100°59′52.0″ESite 1: Sakad/ PuaNan
FLA35.1MW602353Pseudozyma hubeiensis CBS 10077Ps. hubeiensis/99.6718°07′33.4″N 100°18′57.5″ESite 1: Suan Khuen/ MuangPhrae
FLA36.1MW602355Sporidobolus pararoseus CBS 491Sd. pararoseus/99.3018°07′02.0″N 100°17′50.9″ESite 1: Pa Daeng/ MuangPhrae11.2
FLA31.7MW602347Sporobolomyces bannaensis CBS 9204Sb. bannaensis/100.0019°23′06.9″N 100°41′46.7″ESite 1: Na Rai Luang/ Song KhwaeNan22.4
FLA36.2MW602356Sporobolomyces bannaensis CBS 9204Sb. bannaensis/99.1318°07′02.0″N 100°17′50.9″ESite 1: Pa Daeng/ MuangPhrae

Identification of yeast strains isolated from tea flower and their GenBank accession numbers of the D1/D2 domain of LSU rRNA gene.

3.2. Phylogenetic results of new yeast species

The sequences obtained in this study were deposited in the GenBank database. Phylogenetic analyses of a combined ITS and D1/D2 dataset of the genera Metschnikowia and Wickerhamiella were presented in this study. The best substitution models were GTR + I + G for ITS and nrLSU in both analyses. Our phylogenetic analyses revealed that phylograms of the ML and BI analyses in each genus were similar in terms of topology (data not shown). Therefore, the phylogram obtained from the ML analysis was selected and presented in this study.

Phylogenetic analysis of Metschnikowia included 77 sequences of Metschnikowia and two sequences (Eremothecium cymbalariae CBS 1171 and Saccharomyces cerevisiae CBS 270.75) of outgroup (Supplementary S3). The aligned dataset was comprised of 1,093 characters including gaps (ITS: 1–664 and D1/D2: 665–1,093). ML analysis of the combined dataset yielded the best scoring tree with a final ML optimization likelihood value of −13,839.8771. A phylogram clearly separated the two Metschnikowia strains (SDBR-CMU426 and SDBR-CMU427) (introduced as M. lannaensis) obtained in this study from the previously known Metschnikowia species with high support values (BS = 100% and PP = 1.0) (Figure 2). Metschnikowia lannaensis formed a sister taxon to M. hawaiiana and M. miensis with 85% BS and 0.95 PP support values.

Figure 2

Phylogenetic analysis of Wickerhamiella used 52 sequences of Wickerhamiella and two sequences of the outgroup (E. cymbalariae CBS 1171 and S. cerevisiae CBS 270.75) (Supplementary S4). The aligned dataset was comprised of 1,249 characters including gaps (ITS: 1–656 and D1/D2: 657–1,249). ML analysis of the combined dataset yielded the best scoring tree with a final ML optimization likelihood value of −14,850.8649. A phylogram clearly separated eight yeast strains into two monophyletic clades with high support values (BS = 100% and PP = 1.0) (Figure 3). The results indicated that four yeast strains, such as SDBR-CMU428, SDBR-CMU429, SDBR-CMU430, and SDBR-CMU431 (introduced as W. camelliae), were clearly distinguished from the previously known species of Wickerhamiella. It was found that W. camelliae formed a sister taxon to W. azyma, W. azymoides, W. nectarea, and W. parazyma with 95% BS and 1.0 PP support values. Moreover, four yeast strains in this study, such as SDBR-CMU432, SDBR-CMU433, SDBR-CMU434, and SDBR-CMU435 (described here as W. thailandensis) formed a sister clade to W. musiphila with high support (BS = 100% and PP = 1.0).

Figure 3

3.3. Taxonomic description of novel yeast species

3.3.1. Metschnikowia lannaensis

A. Kanpiengjai, P. Kodchasee, K. Unban, J. Kumla, S. Lumyong, P. Khunnamwong, Sarkar, K. Shetty and C. Khanongnuch sp. nov. (Figure 4).

Figure 4

MycoBank No.: MB 845148.

Etymology:lannaensis (lan.na.en’sis) N.L. fem. adj. lanna, referring to an ancient area of Northern Thailand, where the type strain was isolated.

Holotype: FLA8.2T is the holotype of Metschnikowia lannaensis. It was isolated from the tea flower of Camellia sinensis var. assamica collected from the Suthep subdistrict, Muang district, Chiang Mai province, Thailand. It has been preserved in a metabolically inactive state in the culture collection of the Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand. The culture ex-type has been permanently deposited in a metabolically inactive state at the Sustainable Development of Biological Resources, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand (SDBR) as SDBR-CMU426.

Description: After 3 days at 25°C on YM agar, cells are ovoid (2–3.5 × 5.5–7 μm) and occur singly or in pairs. Budding is multilateral. The colony is a circular form (2–3 mm), white to cream in color, convex, smooth, soft to butyrous, and has an entire margin. Pseudohyphae are formed in slide culture on YMA within 4 weeks at 25°C. Ascospore formation is not observed on PDA, corn meal agar, 5% malt extract agar, Fowell’s acetate agar, YCBY agar, YPD agar, and YM agar at 15°C and 25°C for 4 weeks. Fermentation of glucose is negative. D-Glucose, D-galactose, L-sorbose, sucrose, salicin, raffinose, glycerol (or latent positive), ribitol, D-glucitol, D-mannitol, D-glucono-1,5-lactone, 5-ketogluconic acid, succinate, ethanol, and xylitol (or slow positive) are assimilated, but N-acetyl glucosamine, ribose, xylose, L-arabinose, D-arabinose, L-rhamnose, maltose, trehalose, methyl-α-D-glucoside, cellobiose, melibiose, lactose, melizitose, inulin, soluble starch, erythritol, galactitol, myo-inositol, 2-ketogluconic acid, D-gluconate, D-glucuronate, D-galacturonic acid, lactate, citrate, and methanol are not assimilated. For the assimilation of nitrogen compounds, growth on ammonium sulfate, ethylamine HCl, and L-lysine is positive, but growth on nitrate, nitrite, and cadaverine is negative.

Growth in the vitamin-free medium is negative. Growth with 0.01% cycloheximide and 0.1% cycloheximide is negative. No growth occurs on media containing 10% NaCl and 16% NaCl. Growth is observed in the presence of 50% glucose, 60% glucose, and at 20, 25, and 30°C, but not at 35 and 37°C. Soluble starch-like extracellular carbohydrates are not produced. The acid formation is negative. Hydrolysis of urea and DBB reaction is negative.

Additional strains examined: Thailand, Lampang province, Muang Pan district, Chae Son, isolated from flowers of Assam tea (C. sinensis var. assamica), January 2016, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch, SDBR-CMU427.

GenBank accession number: holotype FLA8.2T (D1/D2: KY640629, ITS: OM213034); additional strain SDBR-CMU427 (D1/D2: MW542582, ITS: MZ857157).

Note: Phylogenetically, M. lannaensis formed a sister taxon to M. hawaiiana and M. miensis. However, W. lannaensis could be clearly distinguished from M. miensis and W. hawaiina by its ability to assimilate D-mannitol in addition to its inability to assimilate N-acetyl glucosamine (Shibayama et al., 2020). A pairwise nucleotide comparison of ITS data also indicated that M. lannaensis differed from M. hawaiiana and M. miensis in 19.3% (72/373 bp) and 25.7% (95/369), respectively. In addition, D1/D2 data of M. lannaensis differed from M. hawaiiana in 11.7% (62/529 bp) and M. miensis in 14.5% (78/535 bp).

3.3.2. Wickerhamiella camelliae

A. Kanpiengjai, P. Kodchasee, K. Unban, J. Kumla, S. Lumyong, P. Khunnamwong, D. Sarkar, K. Shetty and C. Khanongnuch sp. nov. (Figure 5).

Figure 5

MycoBank No.: MB 845147.

Etymology:camelliae (ca.mel.li’ae) L. fem. adj. camelliae, from the flower of Camellia sinensis, referring to the substrate from which the species was isolated.

Holotype: FLA5.4T is the holotype of Wickerhamiella camelliae. It was isolated from the tea flower of Camellia sinensis var. assamica collected from the Mae Na subdistrict, Chiang Dao district, Chiang Mai province, Thailand. It has been preserved in a metabolically inactive state in the culture collection of the Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand. The culture ex-type has been permanently deposited in a metabolically inactive state at the Sustainable Development of Biological Resources, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand (SDBR) as SDBR-CMU428.

Description: After 3 days at 25°C on YM agar, cells are ovoid (2–3 × 5–7 μm) and occur singly or in pairs. Budding is multilateral. The colony is a circular form (2–3 mm), white to cream in color, convex, smooth, and has an entire margin. Pseudohyphae or true hyphae are not formed in slide culture on YMA within 4 weeks at 25°C. Ascospore formation is not observed on PDA, corn meal agar, 5% malt extract agar, Fowell’s acetate agar, YCBY agar, YPD agar, and YM agar at 15°C and 25°C for 4 weeks. Fermentation of glucose is negative. D-Glucose, D-galactose, L-sorbose, xylose, L-arabinose, sucrose, maltose, trehalose, methyl-α-D-glucoside, melizitose, glycerol (or slow positive), D-glucitol, D-mannitol, galactitol, D-glucono-1,5-lactone, 2-ketogluconic acid, succinate, and xylitol are assimilated, but N-acetyl glucosamine, D-ribose, L-rhamnose, cellobiose, salicin, melibiose, lactose, raffinose, inulin, soluble starch, erythritol, ribitol, myo-inositol, 5-ketogluconic acid, D-glucuronate, D-galacturonic acid, and methanol are not assimilated. Assimilation is variable for DL-lactate, citrate, ethanol, D-arabinose, and D-gluconate. Nitrogen assimilation is positive for ammonium sulfate, ethylamine HCl, cadaverine, and L-lysine but is negative for nitrate and nitrite.

Growth in the vitamin-free medium is negative. Growth with 0.01% cycloheximide and 0.1% cycloheximide is negative. No growth is observed in the presence of 10% NaCl and 16% NaCl. Growth is observed in the presence of 50% glucose, 60% glucose and at 20, 25, and 30°C but not at 35 and 37°C. Soluble starch-like carbohydrates are not produced. The acid formation is negative. Hydrolysis of urea and DBB reaction is negative.

Additional strains examined: Thailand, Mae Ai subdistrict, Mae Ai, isolated from the flower of Assam tea (C. sinensis var. assamica), December 2015, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch, SDBR-CMU429; Chiang Rai province, Wiang Pa Pao district, Mae Chedi, isolated from the flower of Assam tea (C. sinensis var. assamica), December 2015, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch SDBR-CMU430; Phrae province, Muang district, Suan Khuen, isolated from the flower of Assam tea (C. sinensis var. assamica), March 2016, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch, SDBR-CMU431.

GenBank accession number: holotype FLA5.4T (D1/D2: KY411895, ITS: MZ857161); additional strains SDBR-CMU429 (D1/D2: MW542579, ITS: MZ857162), SDBR-CMU430 (D1/D2: MW542580, ITS: 857163), SDBR-CMU431 (D1/D2: MW542585, ITS: MZ857167).

Note: Based on phylogenetic analysis, W. camelliae formed a sister taxon to W. azyma, W. azymoides, W. nectarea, and W. paraazyma. However, W. camelliae could be clearly distinguished from W. azyma, W. azymoides, W. nectarea, and W. paraazyma by its inability to assimilate cadaverine. A pairwise nucleotide comparison of ITS data also indicated that W. camelliae differed from W. azyma, W. azymoides, and W. parazyma in 17.2% (70/406 bp), 16.2% (66/407 bp), and 19.6% (82/417 bp), respectively. Furthermore, D1/D2 of W. camelliae differed from W. azyma, W. azymoides, W. nectarea, and W. parazyma in 8.3% (50/556 bp), 9.1% (51/556 bp), 8.1% (45/557 bp), and 12.9% (72/557 bp), respectively.

3.3.3. Wickerhamiella thailandensis

A. Kanpiengjai, P. Kodchasee, K. Unban, J. Kumla, S. Lumyong, P. Khunnamwong, D. Sarkar, K. Shetty and C. Khanongnuch sp. nov. (Figure 6).

Figure 6

MycoBank No.: MB 845146.

Etymology:thailandensis thai.land.en’sis. N.L. fem. adj. thailand, refers to Thailand, where the type strain was isolated.

Holotype: FLA4.2T is the holotype of Wickerhamiella thailandensis. It was isolated from the tea flower of Camellia sinensis var. assamica collected from the Pa Pae subdistrict, Mae Taeng district, Chiang Mai province, Thailand. It has been preserved in a metabolically inactive state in the culture collection of the Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, Chiang Mai, Thailand. The culture ex-type has been permanently deposited in a metabolically inactive state at the Sustainable Development of Biological Resources, Faculty of Science, Chiang Mai University, Chiang Mai, Thailand (SDBR) as SDBR-CMU432.

Description: After 3 days at 25°C on YM agar, cells are subglobose to ovoid (2–3.5 × 5.5–7 μm) and occur singly or in pairs. Budding is multilateral. The colony is a circular form (2–3 mm), white to cream in color, convex, smooth, and has an entire margin. Pseudohyphae or true hyphae are not formed in slide culture on YMA within 4 weeks at 25°C. Ascospore formation is not observed on PDA, corn meal agar, 5% malt extract agar, Fowell’s acetate agar, YCBY agar, YPD agar, and YM agar at 15 and 25°C for 4 weeks. Fermentation of glucose is negative. D-Glucose, D-galactose, L-sorbose, sucrose, raffinose, glycerol (or slow positive), ribitol, D-glucitol, D-mannitol, D-glucono-1,5-lactone, 5-ketogluconic acid (or weak positive), DL-lactate (or slow positive), succinate, ethanol (or weak positive), xylitol (or weak positive), and D-gluconate are assimilated, but N-acetyl glucosamine, ribose, xylose, L-arabinose, D-arabinose, L-rhamnose, maltose, trehalose, methyl-α-D-glucoside, cellobiose, salicin, melibiose, lactose, melizitose, inulin, soluble starch, erythritol, galactitol, myo-inositol, 2-ketogluconic acid, D-glucuronate, D-galacturonic acid, citrate, and methanol are not assimilated.

Assimilation is variable for D-gluconate. Nitrogen assimilation is positive for ammonium sulfate, ethylamine HCl, lysine, and cadaverine, but is negative for nitrate and nitrite.

Growth in the vitamin-free medium is negative. Growth with 0.01% cycloheximide and 0.1% cycloheximide is negative. Growth in the presence of 50% glucose and 60% glucose is positive, but growth in the presence of 10% NaCl and 16% NaCl is negative. Growth is observed at 20, 25, and 30°C, but not at 35 and 37°C. Soluble starch-like carbohydrates are not produced. The acid formation is negative. Hydrolysis of urea and DBB reaction is negative.

Additional strains examined: Thailand, Lampang province, Muang Pan district, Chae Son, isolated from the flower of Assam tea (C. sinensis var. assamica), January 2016, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch, SDBR-CMU433, SDBR-CMU434; Chiang Rai province, Thoeng district, Mae Loi, isolated from the flower of Assam tea (C. sinensis var. assamica), January 2016, A. Kanpiengjai, P. Kodchasee, K. Unban, and C. Khanongnuch, SDBR-CMU435.

GenBank accession numbers: holotype FLA4.2T (D1/D2: KY411893, ITS: MZ857160); additional strains SDBR-CMU433 (D1/D2: MW542581, ITS: MZ857164), SDBR-CMU434 (D1/D2: KY640634, ITS: MZ857165), and SDBR-CMU435 (D1/D2: MW542584, ITS: MZ857166).

Note: Phylogenetically, W. thailandensis formed a sister taxon to W. musiphila. However, W. thailandensis could be clearly distinguished from W. musiphila by its ability to assimilate D-mannitol, D-xylose, and glycerol and its inability to assimilate cadaverine (Wang et al., 2008). Furthermore, ITS and D1/D2 data of W. thailandensis differed from W. musiphila in ITS 12.3% (49/397 bp) and in D1/D2 5.6% (31/552 bp).

3.4. Yeast diversity in tea flowers acquired from the Miang production area of Northern Thailand

Based on family identification, all 82 yeast strains could be assigned into 12 families. The frequency aligned each species is as follows: Debaryomycetaceae (34.1%), Trichomonascaceae (19.5%), Metschnikowiaceae (17.1%), Ustilaginaceae (6.1%), Pichiaceae (4.9%), Bulleribasdiaceae (3.7%), Phaffomycetaceae (3.7%), Sporidiobolaceae (3.7%), Rhynchogastremataceae (2.4%), Saccharomycetaceae (2.4%), Saccharomycetales (1.2%), and Saccharomycopsidaceae (1.2%). It was observed that in Debaryomycetaceae, 17.1% frequency of occurrence belonged to Kodamaea/Candida (11.0%), Kurtzmaniella/Candida (3.7%), and Candida/Lodderomyces (2.4%) clades. According to the results, Debaryomyces spp., Candida spp., and Wickerhamiella spp. were recorded as the main yeast population found in tea flowers, while Metschnikowia spp. was considered the second most abundant yeasts. Candida leandrae (10.6%) was the predominant species that was mainly found in samples collected from Chiang Mai. W. azyma was the second most abundant species (9.4%) which was acquired from tea flowers collected from Chiang Mai, Phrae, and Nan provinces. Metschnikowia rancensis, W. thailandensis, and W. camelliae were partially found in tea flowers with an equivalent percentage of occurrence (4.7%). The PCA plot reveals three different correlations of yeast species in tea flowers acquired from six provinces of Northern Thailand (Figure 7). The species diversity of yeasts in tea flowers collected from Chiang Mai, Lampang, and Nan provinces showed a high positive correlation with those acquired from Phayao, Chiang Rai, and Phrae, respectively. Candida leandrae and Ku. quercitrusa were the unique yeast species found in Chiang Mai and Phayao provinces and W. thailandensis was mainly found in tea flowers collected from Lampang province. There is a high opportunity to recover W. azyma from tea flowers acquired from Phrae and Nan provinces. Other yeast species were generally distributed within tea flowers of six provinces of Northern Thailand. These yeast species could be simply divided into two main groups according to the PCA plot, the first group was yeast species that were widely distributed in tea flowers among Chiang Mai, Phayao, Lampang, and Chiang Rai provinces, the second group was yeast species that were markedly found in tea flowers collected from Phrae and Nan provinces.

Figure 7

3.5. Tannin-tolerant yeasts and distribution in Northern Thailand

All yeast strains were assessed for tannin tolerance using YMA supplemented with different concentrations of tannic acid ranging from 10 to 50 g/L. The results revealed that approximately 78% of the tested yeast strains (64 out of 82 yeast strains) were confirmed for their tannin tolerance at a concentration of 50 g/L tannic acid (Figure 8). These yeasts belong to three main families including Saccharomycetaceae, Metschnikowiaceae, and Trichomonascaceae. A total of 12 out of 64 yeast strains exhibited positive results for tannase activity. These yeasts were Cy. (Cyberlindnera) fabianii FLA39.2, D. nepalensis FLA46.3, D. subglobosus FLA34, M. rancensis FLA38.1, M. koreensis FLA41.1, W. thailandensis strains FLA4.2., FLA24.1, FLA27.3, and FLA30.2, Wc. (Wickerhamomyces) anomalus FLA39.1, Pr. (Priceomyces) melissophilus FLA48, and Y. dushanensis FLA45.4. It was observed that tannin-tolerant ability had species variability. The PCA plot reveals similar correlations of yeast species in tea flowers to those previously described (Figure 9). Interestingly, Phrae and Nan provinces were considered the most diverse areas in upper Northern Thailand for the isolation of tannin-tolerant yeast although the tea plantation area and sampling site were in lower numbers than those in Chiang Mai. Wickerhamiell azyma was identified as tannin-tolerant yeast that was uniquely distributed in tea flowers from these areas. The tannin-tolerant yeasts isolated from tea flowers collected from the Phayao province showed poor representation due to the low number of tannin-tolerant yeasts, but the yeast species seem to closely relate with those obtained from samples collected from Phrae and Nan provinces. It was observed that sampling sites of tea flowers from Phayao province align within the eastern part of Northern Thailand where they are in the same location as Phrae and Nan provinces. Here, the distribution of dominant tannin-tolerant yeasts was found within tea flowers collected from Phare, Nan, and Phayao provinces.

Figure 8

Figure 9

4. Discussion

The tea plants that are cultivated for Miang production are usually scattered in the forested areas of the mountainous region of Northern Thailand. Generally, the Miang production area is located within the same area of the tea plantation where it belongs to the Miang entrepreneur. Depending on the area of production, two different production processes have been established which include the filamentous fungi growth-based process (FFP) and the non-filamentous fungi growth-based process (NFP; Khanongnuch et al., 2017). In previous studies on the microbial community of Miang fermentation, lactic acid bacteria, Bacillus spp., yeasts, and filamentous fungi contributed to physico-chemical changes in Miang directly affected flavor and aroma, health-relevant bioactive compounds, and antioxidant activity of the Miang (Unban et al., 2019; Abdullahi et al., 2021). The evolution of microorganisms associated with Miang fermentation has not yet been elucidated through in-depth studies. Flowers offer favorable microenvironments for yeast growth rather than bacteria (Hausmann et al., 2017) and are recognized as rich sources of novel yeast (De Vega et al., 2017). Unique yeast communities occur to be associated with the rhizosphere, the phylloplane, plant extrudates, necrotic tissues, fruits, and flowers. Like other wild flowers, tea flowers remain a comparatively unexplored reservoir of yeast biodiversity.

The highest number of sampling sites was in Chiang Mai province, followed by Phrae and Nan provinces. However, flowers collected from Phrae and Nan provinces showed more diverse yeast species than those in Chiang Mai. The samples FLA11 and FLA31 were acquired from unexplored locations for Miang production due to difficulty in accessing the tea plantations which are located in high mountains. The incidence of specialist nectar-borne yeast species decreases with increasing pollution index due to the limitation of pollinator foraging behavior; thus, urbanization may restrict the movement of nectar-specialized yeasts (Wehner et al., 2017). Nectar serves as a substrate for the proliferation of yeasts especially Ascomycetes with regard to the various nutrient compositions including lipids, proteins, and amino acids (De Vega et al., 2017). This implication was in accordance with the finding of this study in which most of the yeasts found in tea flowers are specifically ascomycete yeasts. Moreover, three new yeast species, namely, M. lannaensis, W. camelliae, and W. thailandensis isolated from the flower of Assam tea in Northern Thailand were proposed based on identification through molecular phylogenetic and phenotypic (morphological, biochemical, and physiological characteristics) analyses. Based on nucleotide comparisons, all new yeasts described in this study exhibited more than a 1.5% difference in ITS and D1/D2 domain regions with their sister taxa. The yeast strain showing nucleotide substitutions of greater than 1% in the D1/D2 domain of the LSU rRNA gene are usually considered different species (Kurtzman and Robnett, 1997). Therefore, M. lannaensis, W. camelliae, and W. thailandensis can be considered different species. The findings that most yeast strains including three new yeast species were isolated from tea flowers agree with the results from the previous studies which stated that the most frequent yeasts that were recovered from floral nectar and insects belong to the genera Metschnikowia, Starmerella, Kodamaea, and Wickerhamiella, but their role on plant–pollinator interactions has not yet been studied (De Vega et al., 2017, 2018). Other yeast genera commonly found in nectar and flora included Aureobasidium, Candida, Clavispora, Cryptococcus, Debaryomyces, Hanseniaspora, Kodamaea, Papiliotrema, Rhodotorula, Starmerella, Sporobolomyces, and Wickerhamiella (Klaps et al., 2020).

The diversity of yeast species in tea flowers was influenced by the sampling site and location. Nan and Phrae provinces were unique areas for the isolation of W. azyma, while C. leandrae was predominantly found in Chiang Mai. Species diversity of yeasts in flowers is correlated with many factors such as location, the density of flowers (Belisle et al., 2012), pollinators, and other flower-visiting animals (Hausmann et al., 2017; Chappell and Fukami, 2018). Similar to nectar-inhabiting yeast communities between Phrae and Nan provinces may be dependent on geographically distant locations. Therefore, the yeast communities in these provinces were different from those in Chiang Mai.

To evaluate yeasts that may be associated with Miang and the evolution of Miang fermentation, all yeast strains were assessed for tannin-tolerant ability. Based on the results, approximately 78% of total yeast strains were able to tolerate high concentrations at up to 50 g/L tannic acid except for C. leandrae, the most abundant yeast found in tea flowers and all strains of Su. (Suhomyces) xylosoci. Some yeasts were identical species reported from both commercial Miang and Miang fermentation processes, i.e., C. tropicalis, Hy. burtonii, Me. caribbica, and P. manshurica but the later species including C. orthopsilosis, Cy. fabianii, H. uvarum, Wc. anomalus were found only in the Miang fermentation process (Kanpiengjai et al., 2016; Unban et al., 2020a; Kodchasee et al., 2021) and Wc. anomalus, a part of the tannin-tolerant yeasts, is capable of producing tannase. Further tannase-producing yeasts found from tea flowers were Cy. fabianii, D. nepalensis, D. subglobosus, M. rancensis, W. thailandensis, Pc. melissophilus, and Y. dushanensis. As described in the previous studies, tannin tolerance and tannase-producing ability are good attributes for microorganisms associated with Miang fermentation (Kanpiengjai et al., 2016; Chaikaew et al., 2017; Unban et al., 2020b). The tannin-tolerant yeasts isolated from Miang, Cy. rhodanensis, D. hansenii, and S. ruineniae have been confirmed for their ability to produce cell-associated tannase that played a crucial role in the biotransformation of bioactive phenolic compounds, namely, epigallocatechin gallate (EGCG) and epicatechin gallate (ECG) (Leangnim et al., 2021). For traditional fermented tea leaves such as Pu-erh tea, Fu brick tea, and Laphet-so, it was proposed that properties of yeasts such as S. cerevisiae, Blastobotrys adeninivorans, Cyberlindnera sp., and Candida sp. may improve the quality of Fu brick tea via their enzyme activities toward the degradation of phenolic compounds and producing the sweet substance xylitol and other flavor substances (Li et al., 2017; Zhao et al., 2018). We found coding regions of feruloyl esterase B (FaeB), which is a member of the tannase family, and tannase in whole genome sequences of the reported yeast detected during tea fermentation such as B. adeninivorans (accession number CAT07374) and Cy. fabianii (accession number ONH67742). Furthermore, our current research revealed FaeB and tannase in the transcriptome of tannin-tolerant and tannase-producing yeasts isolated from Miang, namely, Cy. rhodanensis and Sd. ruineniae (data not shown). Considering the catalytic region of all above-mentioned tannases and FaeBs, GCSXGG sequence at the catalytic region and CS-D-HC motif as the catalytic triads of the enzymes are always conserved for fungal and yeast tannases that have high specificity toward catechins and condensed tannins from the plant cell wall (Suzuki et al., 2014; Dong et al., 2021). It is believed that the coding regions of these enzymes may either be conserved within species or varied among different yeast genera. Although the main yeasts in Miang and Miang fermentation processes were found to be C. ethanolica and P. manshurica, it is assumed that a part of tannin-tolerant yeasts isolated from tea flowers which are conserved within Miang fermentation may be alternatively associated with the fermentation by playing roles in the production of Miang fermentation such as promoting an appropriate environment for tannase-producing yeasts and other microorganisms associated with developing steamed tea leaves into Miang, and enhancement of flavor and aroma to serve a unique Miang product of each Miang production area.

In summary, we present the yeast ecology of tea flowers and suggest that it is partially associated with Miang fermentation based on their tannin tolerance and tannase-producing ability. The yeast ecology is mainly restricted to the flower habitat, and this strengthens the relevance of plant–pollinator interactions. It can be further concluded that tea flowers are a potential source of tannin-tolerant and tannase-producing yeasts, which will shed more light as the important microbial source for the development of Miang fermentation and the production of bioactive compounds from Miang.

Funding

This study was supported financially by the National Research Council of Thailand (NRCT) and also funding from the Royal Golden Jubilee Ph.D. (RGJ) Programme (grant no. PHD/0137/2560). We are also grateful to Chiang Mai University for financial support via a postdoctoral fellowship.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

AK, CK, JK, and SL were responsible for conceptualization and funding acquisition. AK, CK, KU, and PKo collected the species, performed the morphological analysis, and molecular phylogenetic analyses. PKh, AK, and PKo performed biochemical and physiological analyses. AK wrote the first draft. AK, CK, KU, JK, and PKh wrote, reviewed, and edited the manuscript. DS and KS gave some suggestions for English improvement and the flow of concepts. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors would like to acknowledge the Division of Biotechnology, Faculty of Agro-Industry, Chiang Mai University, and Research Center of Microbial Diversity and Sustainable Utilization, Faculty of Sciences, Chiang Mai University for the research facilities.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1043430/full#supplementary-material

References

  • 1

    AbdullahiA. D.KodchaseeP.UnbanK.PattananandechaT.SaenjumC.KanpiengjaiA.et al. (2021). Comparison of phenolic contents and scavenging activities of Miang extracts derived from filamentous and non-filamentous fungi-based fermentation processes. Antioxidants10:1144. doi: 10.3390/antiox10071144

  • 2

    BelisleM.PeayK. G.FukamiT. (2012). Flowers as islands: spatial distribution of nectar-inhabiting microfungi among plants of Mimulus aurantiacus, a hummingbird-pollinated shrub. Microb. Ecol.63, 711718. doi: 10.1007/s00248-011-9975-8

  • 3

    ChaikaewS.BaipongS.SoneT.KanpiengjaiA.Chui-ChaiN.AsanoK.et al. (2017). Diversity of lactic acid bacteria from Miang, a traditional fermented tea leaf in northern Thailand and their tannin-tolerant ability in tea extract. J. Microbiol.55, 720729. doi: 10.1007/s12275-017-7195-8

  • 4

    ChappellC. R.FukamiT. (2018). Nectar yeasts: a natural microcosm for ecology. Yeast35, 417423. doi: 10.1002/yea.3311

  • 5

    DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods9:772. doi: 10.1038/nmeth.2109

  • 6

    De VegaC.AlbaladejoR. G.GuzmánB.SteenhuisenS. L.JohnsonS. D.HerreraC. M.et al. (2017). Flowers as a reservoir of yeast diversity: description of Wickerhamiella nectarea f.a. sp. nov., and Wickerhamiella natalensis f.a. sp. nov. from South African flowers and pollinators, and transfer of related Candida species to the genus Wickerhamiella as new combinations, fox054. FEMS Yeast Res.17. doi: 10.1093/femsyr/fox054

  • 7

    De VegaC.AlbaladejoR. G.LachanceM.-A. (2018). Metschnikowia maroccana fa, sp. nov., a new yeast species associated with floral nectar from Morocco. Int. J. Syst. Evol. Microbiol.68, 20282035. doi: 10.1099/ijsem.0.002784

  • 8

    DongL.MckinstryW. J.PanL.NewmanJ.RenB. (2021). Crystal structure of fungal tannase from aspergillus Niger. Acta Crystallogr. D Struct. Biol.77, 267277. doi: 10.1107/s2059798320016484

  • 9

    EdgarR. C. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics5:113. doi: 10.1186/1471-2105-5-113

  • 10

    FelsensteinJ. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution39, 783791. doi: 10.2307/2408678

  • 11

    HallT. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program forWindows 95/98/NT. Nucleic Acids Symp. Ser.41, 9598.

  • 12

    HausmannS. L.TietjenB.RilligM. C. (2017). Solving the puzzle of yeast survival in ephemeral nectar systems: exponential growth is not enough. FEMS Microbiol. Ecol.93:fix150. doi: 10.1093/femsec/fix150

  • 13

    KanpiengjaiA.Chui-ChaiN.ChaikaewS.KhanongnuchC. (2016). Distribution of tannin-tolerant yeasts isolated from Miang, a traditional fermented tea leaf (Camellia sinensis var. assamica) in northern Thailand. Int. J. Food Microbiol.238, 121131. doi: 10.1016/j.ijfoodmicro.2016.08.044

  • 14

    KhanongnuchC.UnbanK.KanpiengjaiA.SaenjumC. (2017). Recent research advances and ethno-botanical history of Miang, a traditional fermented tea (Camellia sinensis var. assamica) of northern Thailand. J. Ethn. Foods4, 135144. doi: 10.1016/j.jef.2017.08.006

  • 15

    KhunnamwongP.KingphadungK.LomthongT.KanpiengjaiA.KhanongnuchC.LimtongS. (2022). Wickerhamiella nakhonpathomensis f.a. sp. nov., a novel ascomycetous yeast species isolated from a mushroom and a flower in Thailand. Int. J. Syst. Evol. Microbiol.72:005191. doi: 10.1099/ijsem.0.005191

  • 16

    KlapsJ.LievensB.Álvarez-PérezS. (2020). Towards a better understanding of the role of nectar-inhabiting yeasts in plant–animal interactions. Fungal Biol. Biotechnol.7:1. doi: 10.1186/s40694-019-0091-8

  • 17

    KodchaseeP.NainK.AbdullahiA. D.UnbanK.SaenjumC.ShettyK.et al. (2021). Microbial dynamics-linked properties and functional metabolites during Miang fermentation using the filamentous fungi growth-based process. Food Biosci.41:100998. doi: 10.1016/j.fbio.2021.100998

  • 18

    KumlaJ.NundaengS.SuwannarachN.LumyongS. (2020). Evaluation of multifarious plant growth promoting trials of yeast isolated from the soil of Assam tea (Camellia sinensis var. assamica) plantations in Northern Thailand. Microorganisms8:1168. doi: 10.3390/microorganisms8081168

  • 19

    KurtzmanC. P.FellJ. W.BoekhoutT.RobertV. (2011). “Chapter 7 – Methods for isolation, phenotypic characterization and maintenance of yeasts” in The Yeasts. eds. KurtzmanC. P.FellJ. W.BoekhoutT.. 5th ed (London: Elsevier), 87110.

  • 20

    KurtzmanC. P.RobnettC. J. (1997). Identification of clinically important ascomycetous yeasts based on nucleotide divergence in the 5′ end of the large-subunit (26S) ribosomal DNA gene. J. Clin. Microbiol.35, 12161223. doi: 10.1128/jcm.35.5.1216-1223.1997

  • 21

    LeangnimN.AisaraJ.UnbanK.KhanongnuchC.KanpiengjaiA. (2021). Acid stable yeast cell-associated tannase with high capability in gallated catechin biotransformation. Microorganisms9:1418. doi: 10.3390/microorganisms9071418

  • 22

    LiQ.HuangJ.LiY.ZhangY.LuoY.ChenY.et al. (2017). Fungal community succession and major components change during manufacturing process of Fu brick tea. Sci. Rep.7:6947. doi: 10.1038/s41598-017-07098-8

  • 23

    MillerM. A.PfeifferW.SchwartzT. (2010). Creating the cipres science gateway for inference of large phylogenetic trees. Proceedings of the 2010 Gateway Computing Environments Workshop (GCE); 2010 November 14; New Orleans, LA, USA. IEEE: Manhattan, NY.

  • 24

    OsawaR.WalshT. (1993). Visual reading method for detection of bacterial tannase. Appl. Environ. Microbiol.59, 12511252. doi: 10.1128/aem.59.4.1251-1252.1993

  • 25

    RambautA. (2016). FigTree Tree Figure Drawing Tool Version 131. Institute of Evolutionary Biology, University of Edinburgh. Available at: http://treebioedacuk/software/figtree/ (Accessed October 2, 2022).

  • 26

    RonquistF.HuelsenbeckJ. P. (2003). MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics19, 15721574. doi: 10.1093/bioinformatics/btg180

  • 27

    SchaefferR. N.MeiY. Z.AndicoecheaJ.MansonJ. S.IrwinR. E. (2017). Consequences of a nectar yeast for pollinator preference and performance. Funct. Ecol.31, 613621. doi: 10.1111/1365-2435.12762

  • 28

    ShibayamaK.OtoguroM.NakashimaC.YanagidaF. (2020). Metschnikowia miensis f.a., sp. nov., isolated from flowers in Mie prefecture, Japan. Antonie van Leeuwenhoek113, 321329. doi: 10.1007/s10482-019-01340-1

  • 29

    StamatakisA. (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics22, 26882690. doi: 10.1093/bioinformatics/btl446%JBioinformatics

  • 30

    SuzukiK.HoriA.KawamotoK.ThanguduR. R.IshidaT.IgarashiK.et al. (2014). Crystal structure of a feruloyl esterase belonging to the tannase family: a disulfide bond near a catalytic triad. Proteins82, 28572867. doi: 10.1002/prot.24649

  • 31

    UnbanK.KhatthongngamN.PattananandechaT.SaenjumC.ShettyK.KhanongnuchC. (2020a). Microbial community dynamics during the non-filamentous fungi growth-based fermentation process of Miang, a traditional fermented tea of North Thailand and their product characterizations. Front. Microbiol.11:1515. doi: 10.3389/fmicb.2020.01515

  • 32

    UnbanK.KhatthongngamN.ShettyK.KhanongnuchC. (2019). Nutritional biotransformation in traditional fermented tea (Miang) from North Thailand and its impact on antioxidant and antimicrobial activities. J. Food Sci. Technol.56, 26872699. doi: 10.1007/s13197-019-03758-x

  • 33

    UnbanK.KodchaseeP.ShettyK.KhanongnuchC. (2020b). Tannin-tolerant and extracellular tannase producing Bacillus isolated from traditional fermented tea leaves and their probiotic functional properties. Foods9:490. doi: 10.3390/foods9040490

  • 34

    WangS.-A.JiaJ.-H.BaiF.-Y. (2008). Candida alocasiicola sp. nov., Candida hainanensis sp. nov., Candida heveicola sp. nov. and Candida musiphila sp. nov., novel anamorphic, ascomycetous yeast species isolated from plants. Antonie Van Leeuwenhoek94, 257265. doi: 10.1007/s10482-008-9238-y

  • 35

    WangD.MengJ.GaoH.XuK.XiaoR.ZhongY.et al. (2013). Evaluation of reproductive and developmental toxicities of Pu-erh black tea (Camellia sinensis var. assamica) extract in Sprague Dawley rats. J. Ethnopharmacol.148, 190198. doi: 10.1016/j.jep.2013.04.009

  • 36

    WehnerJ.MittelbachM.RilligM. C.VerbruggenE. (2017). Specialist nectar-yeasts decline with urbanization in Berlin. Sci. Rep.7, 18. doi: 10.1038/srep45315

  • 37

    ZhangL.HoC.-T.ZhouJ.SantosJ. S.ArmstrongL.GranatoD. (2019). Chemistry and biological activities of processed Camellia sinensis teas: a comprehensive review. Compr. Rev. Food Sci. Food Saf.18, 14741495. doi: 10.1111/1541-4337.12479

  • 38

    ZhaoZ.-J.SuiY.-C.WuH.-W.ZhouC.-B.HuX.-C.ZhangJ. (2018). Flavour chemical dynamics during fermentation of kombucha tea. Emirates J. Food Agric.30, 732741. doi: 10.9755/ejfa.2018.v30.i9.1794

Summary

Keywords

Miang, tannin-tolerance, tannase, taxonomy, tea flower, yeast

Citation

Kanpiengjai A, Kodchasee P, Unban K, Kumla J, Lumyong S, Khunnamwong P, Sarkar D, Shetty K and Khanongnuch C (2023) Three new yeast species from flowers of Camellia sinensis var. assamica collected in Northern Thailand and their tannin tolerance characterization. Front. Microbiol. 14:1043430. doi: 10.3389/fmicb.2023.1043430

Received

13 September 2022

Accepted

19 January 2023

Published

16 February 2023

Volume

14 - 2023

Edited by

Samantha Chandranath Karunarathna, Qujing Normal University, China

Reviewed by

Zhipeng Wang, Qingdao Agricultural University, China; Nalin Nilusha Wijayawardene, Qujing Normal University, China; Kasun Thambugala, University of Sri Jayewardenepura, Sri Lanka

Updates

Copyright

*Correspondence: Chartchai Khanongnuch,

This article was submitted to Microbe and Virus Interactions with Plants, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics