茶树根际土壤细菌多样性及其主要环境影响因子研究

茶叶通讯　第48卷　第1期
Journal of Tea Communication Vol.48, No.1
投稿平台：
茶树根际土壤细菌多样性及其主要环境影响因子研究
贺屹成，王梦姣*，秦公伟
陕西理工大学 生物科学与工程学院，陕西 汉中 723000
摘　要：采用高通量扩增子测序法对秦巴山区镇巴、宁强和城固三地典型茶园茶树根际细菌多样性和茶树根际微环境进行分析。

茶树根际细菌多样性研究结果表明，采样地样品中所有土壤茶树根际细菌分布在41个门（其中13个门相对丰度大于1%）、1288个属（其中 33个属相对丰度大于1%）。

镇巴采样地2年生茶树根际细菌数量最少，多样性最差。

茶树根际微环境研究结果表明，所有采样地土壤的营养成分均能满足茶树生长需求，但城固采样地需要进行少量有机碳的补充；三个采样地土壤pH 值在4.5～5.5，含水量呈现出不同变化趋势。

茶树根际细菌多样性和根际微环境相关性分析结果表明，影响茶树根际细菌多样性的主要环境因子分别为土壤有机碳、有机氮及水解性氮含量。

关键词：茶树，秦巴山区，根际微生物，根际土壤元素含量，相关性
中图分类号：S571.1, S182 文献标识码：A 文章编号：1009-525X （2021）01-40-49
Study on Bacterial Diversity and Its Main Environmental Impact Factors in
Rhizosphere Soil of Tea plant
HE Yicheng ，WANG Mengjiao *，QIN Gongwei
School of Biological Science and Engineering, Shaanxi University of Technology, Hanzhong 723000, China
Abstract ：Rhizosphere bacteria diversity and rhizosphere microenvironment of tea plants in Zhenba, Ningqiang, and Chenggu in Qinba Mountain area were analyzed by high throughput amplino sequencing method. The results of bacteria diversity in the rhizosphere of tea plant showed that, the bacteria in the samples were distributed in 41 phyla (13 phyla with relative abundance ratio greater than 1%) and 1288 genera (33 genera with relative abundance ratio greater than 1%) were found. The number and diversity of rhizosphere bacteria of two years old tea plants in the sampling area of Zhenba were the least and the worst. The results of rhizosphere microenvironment of tea plants showed that the nutrients in the soil of all sampling sites could meet the growth demand of tea plants, but a small amount of organic carbon supplement was needed in the Chenggu solid sampling sites; The soil pH values in the three sampling sites ranged from 4.5 to 5.5, and soil water content showed different trends. The results of correlation analysis between rhizosphere bacterial diversity and rhizosphere microenvironment showed that the main environmental factors affecting rhizosphere bacterial diversity were soil organic carbon, organic nitrogen and hydrolytic nitrogen.
Key words ：Tea plant, Qinba mountain area, Rhizospheric microorganism, Rhizosphere soil element content, Correlation
收稿日期：2020-06-28 修订日期：2020-09-14
基金项目：陕西省教育厅专项（19JK0191），陕西省高校科协人才托举计划项目（20200206），国家自然科学基金（31870488、41571300），陕西省大学生创新创业训练项目（S201910720056），陕西省高层次人才特殊支持计划（2017）作者简介：贺屹成（1998 -），男，陕西西安人，在读本科生，研究方向：微生物学。

*通讯作者：王梦姣（1987 -），女，陕西宝鸡人，讲师，研究方向：茶树根际微环境。

E-mail:*****************
贺屹成, 王梦娇, 秦公伟, 等. 茶树根际土壤细菌多样性及其主要环境影响因子研究[J]. 茶叶通讯, 2021, 48(1): 40-49.
HE Yicheng, WANG Mengjiao, QIN Gongwei. Study on Bacterial Diversity and Its Main Environmental Impact Factors in Rhizosphere Soil of Tea plant[J]. Journal of Tea Communication, 2021, 48(1): 40-49.
41贺屹成等：茶树根际土壤细菌多样性及其主要环境影响因子研究
第1期
表1　采样信息表
Table 1　Sampling Information
采样地编号茶树品种茶树生长年份纬度（北纬）经度（东经）采样时间
镇巴
B1
福鼎大白茶20年以上32°12′30″108°0′30″
2019/6/20
B210年32°12′28″108°1′15″B32年32°12′32″108°1′50″
宁强N120年以上32°49′52″106°17′41″
2019/6/24 N210年32°50′46″106°17′25″
N32年32°51′10″106°18′11″
城固C120年以上33°7′40″107°14′41″
2019/6/30 C210年33°2′20″107°14′17″
C32年33°2′22″107°14′32″
1.2 茶树根际土壤理化指标测定
茶树根际土壤有机碳、总氮、总磷、全钾、水解性氨、有效磷、速效钾的含量测定分别参考森林土壤有机碳、氮、磷和钾的测定方法[12-15]。

土壤含水量及土壤pH值测定参考乔帅等[16]方法。

1.3 茶树根际土壤细菌多样性分析
采用高通量测序法测定分析茶树根际土壤细菌多样性，DNA 提取、PCR扩增、产物电泳检测及切胶纯化、构建DNA库、OTU数据的获取及测序过程以及OUT序列注释过程参考W ANG等[2]方法。

1.4 数据统计分析
1.4.1 系统进化树构建及多样性分析
选取总体相对丰度排在前50且有属分类信息的OTU 代表性序列，构建系统发育树，系统发育树构建参考Saitou[17]。

茶树[Camellia sinensis (L.) O. Kuntze]是我国重要经济作物，茶叶产业成为各主产茶区的支柱产业和部分贫困地区农民脱贫致富的重要产业。

但在茶叶产业发展过程中，也出现了一些严重掣肘问题，如盲目使用化学肥料导致的茶园土壤土质退化、高剂量化学农药施用导致茶叶农药残留日益严重等。

探索能显著降低化学肥料与农药施用量的新技术、新方法成为促进茶叶产业可持续发展的重要研究课题。

茶树根际土壤微环境的稳定性是茶树正常生长及茶叶品质的保证[1]，土壤微生物的多样性和稳定性还能促进恶劣土壤环境的生物修复[2]。

茶树根际土壤细菌，作为茶树根际土壤微环境中营养元素循环的重要驱动力，对土壤肥力和茶树正常生长具有显著促进作用[3]。

例如，某些根际促生菌能够分泌生长素类物质促进植物生长[4]；某些具有较高固氮酶活性的菌株可以进行生物固氮，提供植物生长必需的氮元素[5]和提高土壤氮素含量[6]；根际促生微生物的旺盛代谢可以促进土壤有机物的分解和营养元素的矿化，改善土壤营养环境[7]，土壤营养元素的变化又会显著影响土壤根际微生物的群落组成[8]，从而形成一个良性循环。

不仅如此，许多根际微生物通过分泌Phenazines及其衍生物、Tropolone和Pyocyunim等抑菌物质[9-10]，达到预防和降低某些植物病害的效果。

试验从根际微环境出发，采用高通量扩增子测序法对茶树根际细菌的结构和功能进行研究，并结合茶树根际土壤营养元素含量、pH值及含水量，分析土壤环境因素与根际细菌多样性的相互关系，以期为土壤改良和微生物菌肥开发提供理论基础和参考依据。

1 材料和方法
1.1 采样方法、采样地及茶树概况
采样点为秦巴山区的镇巴、宁强和城固三地的典型茶园（表1）。

选择距离根系5～15 cm的根周土壤，去除土壤表面异物后，取距离地表10～20 cm深度的土壤样品；每处茶园不同树龄样品均按“S”形补点，取10份土壤样品[11]。

收集土壤样品后迅速带回实验室进行充分混匀并过2 mm筛，分装、标记后分别保存在4 ℃冰箱、-80 ℃冰箱。

第48卷
茶叶通讯
42
使用QIIME（V1.9.1）软件中的alpha_diversity. py脚本对OTU序列进行分析，并计算α多样性指数（Chao1、Shannon以及Simpson）。

利用非度量多维尺度法（Nonmetric Multi-Dimensional Scaling, NMDS）计算样品的Weighted Unifrac 和Uweighted Unifrac距离，分析其β多样性。

1.4.2 NMDS（Non-metric Multidimensional scaling）分析
在beta多样性距离（Weighted Unifrac、Unweighted Unifrac）的计算基础上，基于Weighted Unifrac 和Unweighted Unifrac 距离矩阵，使用QIIME 软件和ggplot 2 软件包进行NMDS 分析并绘制出NMDS图。

1.4.3 茶树根际理化指标分析
茶树根际土壤元素含量、酸碱度和含水量取3次平行实验平均值，用SPSS和R软件，计算不同采样地之间物理化学指标的平均值、方差及其之间的差异性，p < 0.01表示组别之间差异性显著，用小写字母标明。

1.4.4 相关性分析
RDA/CCA 分析：使用R 软件对环境因子进行DCA（Detrended correspondence analysis）分析，根据梯度值确定模型类型，并分析绘图。

分别用SPSS（Statistical Product and Service Solutions）和R软件，进行环境因子与高通量测序结果中细菌α多样性指数及相对丰度大于1%的门的相关性分析。

2 结果与分析
2.1 茶树根际土壤细菌群落结构和多样性
2.1.1 茶树根际土壤细菌群落结构
高通量测序结果表明，三个茶园土样中细菌分布在41个门中，其中相对丰度大于1%的有13个门（图1A），变形菌门、酸杆菌门和拟杆菌门为优势菌门。

变形菌门所占百分比在镇巴及城固土样中随茶树树龄增长先下降后增高，在宁强土样中先升高后下降；酸杆菌门所占百分比在镇巴及城固土样中随茶树树龄增长先上升后下降，在宁强土样中则持续上升；拟杆菌门所占百分比在三个采样地均随茶树树龄增长而持续下降。

三个茶园茶树根际细菌分布于1288个属中，其中相对丰度占比大于1%的属有33个（图1B）。

地杆菌属（16.54%）为镇巴2年树龄茶树根际细菌的优势属,代尔夫特氏属（8.31%）为城固20年以上树龄茶树根际细菌的优势属。

选取所有采样地不同树龄根际土壤中总体相对丰度排在前50 且有属分类信息的OUT代表性序列构建系统发育树（图1C），并将其相对丰度进行可视化展示。

结果表明，Delftia OTU3是总体相对丰度最高的序列。

1 ~ 34分别代表相对丰度小于1%及未鉴定种等的总和、Bryobacter、Solibacter、RB41、嗜酸栖热菌属、Parafrigoribacterium、Paeniglutamicibacter、黄杆菌属、地杆菌属、芽孢八叠球菌属、芽单胞菌属、硝化螺菌属、伯克氏菌属、反硝化卡斯特兰尼氏属、代尔夫特氏属、Massilia、Chujaibacter、Rhodanobacter、寡养单胞菌属、Udaeobacter、Chthoniobacter、Granulicella、Koribacter、Mucilaginibacter、HSB_OF53-F07、慢生根瘤菌属、贪嗜菌属、Acidipila、Lysobacter、假单胞菌属、Haliangium、Pseudolabrys、Acidibacter、Terrimonas。

2.1.2 茶树根际土壤细菌多样性
土壤样品中细菌α多样性指数结果（图2A）表明，三个采样地不同树龄物种数量变化趋势不完全一致，镇巴采样地10年树龄茶树根际细菌数量最多；城固和宁强采样地不同树龄茶树根际细菌数量有一定变化，但是变化不显著。

香农指数和辛普森指数结果表明（图2B、图2C），三个采样地不同树龄茶树根际细菌多样性变化不一致，镇巴采样地2年树龄茶树根际细菌多样性最差，城固采样地10年树龄茶树根际细菌多样性最差，宁强采样地20年以上树龄茶树根际细菌多样性最差。

采用NMDS法对样品的Beta多样性进行分析，结果表明（图3），树龄越大，样品之间根际细菌群落多样性差异越小，不同采样地20年以上树龄茶树根际样品的细菌多样性差异最小。

第1期
贺屹成等：茶树根际土壤细菌多样性及其主要环境影响因子研究
43
茶叶通讯
第48卷44
图1　三地茶园茶树根际土壤细菌群落组成分析
Figure 1　Analysis of bacterial community composition in rhizosphere soil of tea trees in three tea gardens
45贺屹成等：茶树根际土壤细菌多样性及其主要环境影响因子研究
第1期
2.2 茶树根际土壤理化性质
茶树根际土壤元素含量分析结果（图4）表明，镇巴和城固两地茶树根际土壤的有机碳、总氮、水解性氮、有效磷以及速效磷等元素含量都随树龄增加而呈现下降的趋势，个别元素含量在20年以上树龄的茶树根际土壤中略有上升（图4A、图4B、图4E、图4F、图4G）。

宁强采样地10年树龄茶树根际土壤中有机碳、总氮、水解性氮及有效磷等元素含量均为最高（图4A、图4B、图4E、图4G）；2年树龄茶树根际土壤中总磷含量、有效磷含量及全钾含量均为最高（图4C、图4D、图4F）。

三地不同树龄土壤样品pH值无显著差异性（图4H）。

土壤含水量呈现出不同变化趋势，含水量最高的分别是镇巴采样地10年树龄茶树根际土壤、城固2年树龄茶树根际土壤及宁强
图2　三地茶园茶树根际土壤细菌Alpha多样性分析
Figure 2　Analysis of α-bacterial diversity in rhizosphere soil of tea trees in three tea gardens
其他样地数据于B1样地数据进行差异性分析，两组数据没有差异，则用相同小写字母标注；如果两组数据具有显著
性差异，则用不同小写字母标注，下同
图3　利用非度量多维尺度法三地茶园茶树根际土壤细菌Beta多样性分析
Figure 3　Analysis of Beta-bacterial diversity in rhizosphere soil of tea trees in three tea gardens by nonmetric multidimensional
scaling method
茶叶通讯
第48卷46
图 4　茶树根际土壤元素含量、pH值及含水量（数据差异性标注同图1和图2）Figure 4　Analysis of element content, pH value and water content in rhizosphere soil of tea trees in three tea gardens
47
贺屹成等：茶树根际土壤细菌多样性及其主要环境影响因子研究第1期20年以上树龄茶树根际土壤（图4I ）。

2.3 根际土壤理化性质与细菌多样性的相关性
采用RDA/CCA 方法对土壤理化性质与细菌多样性进行相关性分析，结果（图5）表明，对采样地根际微生物多样性影响最大的前四位环境因子分别为土壤总氮含量、pH 值、有机碳及水解性氮。

对环境因子与细菌多样性指数及13个主要细菌门相关性分析结果表明，有机碳、总氮和
水解氮与chao1和shaanon 指数呈极显著或显著负相关，但是对simpson 指数相关性不显著；土壤pH 值与simpson 指数显著相关（表2）。

有机碳、总氮、水解性氮、速效钾含量与各细菌门相关性显著或极显著，其余环境因子与细菌门的显著或极显著相关数量均远小于50%；醋酸杆菌门、放线菌门、Patescibacteria 的相对丰度受到环境因子影响大（表3）。

表2　环境因子与根际土壤细菌α多样性指数相关性
Table 2　Correlation between environmental factors and Alpha diversity index of rhizosphere soil bacteria
环境因子Chao1Shannon Simpson 有机碳含量-0.458*-0.401*-0.052总氮含量-0.572**-0.554**-0.271总磷含量 0.1100.3230.289全钾含量0.420*0.2520.114水解性氮含量-0.497**-0.438*-0.053有效磷含量-0.015-0.080-0.186速效钾含量-0.269-0.2200.024pH 值 0.1030.3150.391*土壤含水量
-0.009
-0.361
-0.269
注：当p < 0.05时，用*表示，当 p < 0.01时用**表示，下同
图5　RDA/CCA 法分析细菌多样性与环境因子关系
Figure 5　
Relationship between bacterial diversity and environmental factors analyzed by RDA/CCA method
第48卷
茶叶通讯
48
3 结论与讨论
本研究采用高通量测序法对秦巴山区镇巴、宁强和城固三地的典型茶园中不同树龄根际土壤细菌多样性进行分析，共检测到13个主要的细菌门（门的相对丰度大于1%），这些主要细菌门在每个采样地分布总丰度范围为80.25%～98.95%；变形菌门、酸杆菌门以及拟杆菌门为优势菌门。

这一结果与前人研究[18]相比，分离得到的门的种类多，优势菌门占比低，充分反映了样地茶树根际土壤细菌多样性较为丰度。

伯克氏菌、慢生根瘤菌、Chthoniobacter, Flavobacterium, Gemmatimonas, Massilia以及 RB41 都是常见的根际促生菌[19]，这些菌都能在本研究中找到，说明本研究样品采集地具有开发根际促生微生物的可能性。

不同树龄茶树根际土壤细菌多样性变化趋势表明，镇巴2年生树龄茶树根际土壤细菌丰富度最低，城固则为10年树龄、宁强20年树龄的茶树根际土壤细菌丰度最低。

值得注意的是，虽然不同采样地细菌丰富度指数变化趋势不一致，但是除去组中丰富度最低的树龄，本组中剩余两个年份树龄茶树根际细菌丰富度变化不显著。

因此，一方面可以得出，以上变化趋势可能和采样地细菌数量有一定关联；另一方面也说明，受到细菌数量的影响，在判断具体哪个采样地、哪个树龄茶树根际细菌多样性丰度最佳，还需继续深入研究。

对样品的Beta多样性分析结果表明，树龄越长，茶树根际土壤细菌结构组成越相似。

但是不同地区不同树龄的根际土壤细菌多样性变化不完全一致，这可能是不同地区根际土壤性质有一定差别，如营养元素含量、胶体结构、土壤性质、土壤颗粒度等，对细菌多样性有一定影响[20]。

因长期施肥和茶树根部有机酸的分泌，树龄较大的茶树根际土壤会出现明显的土壤酸化、板结、正常土壤团粒被破坏等现象，导致茶树根系延展困难，吸收功能降低，茶树长势逐渐减弱，分泌物逐渐降低，与根际微生物互作显著下降[21]。

细菌是土壤微生物，尤其是根际微生物中的主要种群，但是大多数细菌的最适酸碱度是中性或中性偏碱性，酸性土壤也不适合大部分细菌生存。

综合分析，本研究认为，在进行后续秦巴山区典型茶树根际土壤微环境与茶树互作，尤其是根际土壤细菌多样性与茶树的相互关系研究中，建议将10年树龄茶树根际土壤作为研究对象。

本研究对土壤环境因子与茶树根际微生物的相关性研究结果表明，土壤有机碳含量、总氮含量、水解性氮含量能够显著影响茶树根际土壤细菌多样性，pH值与速效钾与茶树根际土壤细菌多样性有一定影响，与O'Brien和Pasternak的结果较为一致[22-23]。

参考文献
[1] STUUTTER MI, Richards S. Relationships between soil
physicochemical, microbiological properties, and nutrient release in buffer soils compared to field soils[J]. Journal of
表3　环境因子与根际土壤细菌门的相关性
Table 3　Correlation between environmental factors and microbial phylum in rhizosphere soil
环境因子酸杆
菌门
放线
菌门
装甲
菌门
拟杆
菌门
绿弯
菌门
厚壁
菌门
芽单胞
菌门
硝化螺
旋菌门
Patesci-
bacteria
浮霉
菌门
变形
菌门
疣微球
菌门
WPS2
有机碳含量-0.622**0.485*-0.422*0.724**-0.0580.549**-0.252-0.453*-0.616**-0.878**0.477*-0.300-0.167总氮含量-0.701**0.429*-0.406*0.858**-0.0380.677**-0.202-0.311-0.527**-0.780**0.454*-0.435*-0.270总磷含量-0.594**0.694**0.1020.3060.480*0.1580.551**0.481*-0.2320.062-0.1020.067-0.532**全钾含量0.1160.007-0.301-0.024-0.102-0.0580.666**0.332-0.1900.346-0.315-0.1760.233水解性氮含量-0.602**0.214-0.3760.662**-0.1410.515**-0.343-0.387*-0.470*-0.885**0.632**-0.235-0.136有效磷含量0.176-0.083-0.258-0.060-0.439*0.0550.567**0.3540.1070.284-0.147-0.501**0.549**速效钾含量-0.683**0.719**-0.3060.736**0.1290.487*-0.129-0.338-0.624**-0.669**0.267-0.113-0.425* pH值-0.563**0.679**0.1330.3630.598**0.0510.1160.052-0.394*-0.171-0.0810.375-0.668**土壤含水量0.320-0.646**-0.3310.021-0.527**0.092-0.412*-0.386*-0.020-0.2120.233-0.1550.363
49贺屹成等：茶树根际土壤细菌多样性及其主要环境影响因子研究
第1期
Environmental Quality, 2012, 41(2): 400-409.
[2] WANG MJ, DENG BW, FU X, et al. Characterizations of
microbial diversity and machine oil degrading microbes in machine oil contaminated soil[J]. Environmental Pollution, 2019, 255: 113190.
[3] DENNIS PG, MILLER AJ, Hirsch PR. Are root exudates more
important than other sources of rhizodeposits in structuring rhizosphere bacterial communities? [J] Fems Microbiology Ecology, 2010, 72(3): 313-327.
[4] 王楠，李刚强，李云龙，等. 固氮类芽孢杆菌的分离鉴定
及其促生、抑菌活性的测定[J]. 中国农业科技导报，2019, 21(5): 95-103.
[5] BAUDOIN E, BENIZRI E, GUCKET A. Impact of artificial
root exudates on the bacterial community structure in bulk soil and maize rhizosphere[J]. Soil Biology and Biochemistry, 2003, 35(9), 1183-1192.
[6] GARCIA-SANCHEZ M, GARCIA-ROMERA I, CAJTHAML
T, et al. Changes in soil microbial community functionality and structure in a metal-polluted site: The effect of digestate and fly ash applications[J]. Journal of Environmental Management, 2015, 162, 63-73.
[7] QIN S, FENG WW, ZHANG YJ, et al. Diversity of Bacterial
Microbiota of Coastal Halophyte Limonium sinense and Amelioration of Salinity Stress Damage by Symbiotic Plant Growth-Promoting Actinobacterium Glutamicibacter halophytocola KLBMP 5180[J]. Applied and Environmental Microbiology, 2018, 84(19): 11-21.
[8] GUO A, DING L, TANG Z, et al. Microbial response to CaCO3
application in an acid soil in southern China[J]. Journal of Environmental Sciences-China, 2019, 79: 321-329.
[9] WELLER DM. Colonization of wheat roots by a fluorescent
pseudo-monad on suppressive to take-all[J].Phytopathology, 1983, (73): 1548-1553.
[10] HARIPRASAD P, CHANDRASHEKAR S, SINGH SB, et al.
Mechanisms of plant growth promotion and disease suppression by Pseudomonas aeruginosa strain 2apa[J].Journal of Basic Microbiology, 2014, 54(8): 792-801.
[11] 农业部全国技术推广服务中心. 土壤分析技术规范（第二
版）[M]. 北京: 农业出版社, 2009.
[12] LY/T 1237-1999 森林土壤有机碳的测定及碳氮比的计算
[S]. 北京: 中国出版社, 2000.
[12] LY/T 1228-2015, 森林土壤氮的测定[S]. 北京: 中国出版社,
2016.
[13] LY/T 1232-2015 森林土壤磷的测定[S]. 北京: 中国出版社,
2016.
[15] LY/T 1234-2015森林土壤钾的测定[S]. 北京: 中国出版社,
2016.
[16] 乔帅, 王梦姣, 邓百万, 等. 轮作区水稻根际土壤钙镁离子
含量、含水量和酸碱度变化趋势[J]. 江苏农业科学，2017, 45(5): 284-288.
[17] SAITOU N, NEI M. The neighbor-joining method: a new
method for reconstructing phylogenetic trees[J]. Mol Biol Evol, 1987, 4: 406-425.
[18] SAMET M, KARRAY F, MHIRI N, et al. Effect of
phosphogypsum addition in the composting process on the physico-chemical proprieties and the microbial diversity of the resulting compost tea[J]. Environmentak Science and Pollution Research, 2019, 26(21):21404-21415.
[19] AGUIRRE-VON-WOBESER, E., Rocha-Estrada, et
al. Enrichment of Verrucomicrobia, Actinobacteria and Burkholderiales drives selection of bacterial community fro m soil by maize roots in a traditional milpa agroecosystem[J].
PLoS One, 2018, 13(12), e0208852.
[20] 伍丽,余有本,周天山,等. 茶树根际土壤因子对根际微生
物数量的影响[J]. 西北农业学报, 2011, 20(4):159-16. [21] 殷佳丽. 茶园土壤团聚体氟、铝形态分布及其与土壤化学
性质的关系[D]. 成都: 四川农业大学, 2016.
[22] GREENLON A, CHANG PL, DAMTEW ZM, et al. Global-
level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria[J]. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116 (30): 15200-15209.
[23] LI Q, LI Y, LUO Y, et al. Shifts in diversity and function of the
bacterial community during the manufacture of Fu brick tea[J].
Food Microbiology, 2019, 80: 70-76.。