扩展功能
文章信息
- 陈文峰
- CHEN Wen-Feng
- 根瘤菌系统学研究进展与展望
- Progress and perspective of systematics of rhizobia
- 微生物学通报, 2016, 43(5): 1095-1100
- Microbiology China, 2016, 43(5): 1095-1100
- DOI: 10.13344/j.microbiol.china.150994
-
文章历史
- 收稿日期: 2015-12-01
- 接受日期: 2015-12-31
- 优先数字出版日期(www.cnki.net): 2016-01-04
根瘤菌是一类革兰氏阴性细菌,其主要特征是能与豆科植物共生固氮,但近年来的研究发现,许多非共生的细菌根据系统发育分析也归属到根瘤菌的范围内,使得根瘤菌既包括了共生固氮的细菌,也包括了非共生固氮的细菌。2001年,Moulin等[1]发现分离自含羞草的根瘤菌的分类地位与传统观念上认为的归属于α-变形杆菌纲的根瘤菌属(Rhizobium)、慢生根瘤菌属(Bradyrhizobium)、中华根瘤菌属(Sinorhizobium)、中慢生根瘤菌属(Mesorhizobium)等不同,而是归属到β-变形杆菌纲的伯克霍尔德氏菌属(Burkholderia )和贪铜菌属(Cupriavidus)。近几年,随着DNA测序技术的进步,多基因序列分析和全基因组序列测定已广泛应用于包括根瘤菌在内的原核生物系统学研究当中。之前曾发现一株根瘤菌——Sinorhizobium sp. NGR234能与120个属的豆科植物共生;而最新的研究表明,一种药用豆科植物——苦参(Sophora flavescens)能与α-纲和β-纲的6个属超过30个种的根瘤菌共生,彻底改变了之前人们认为根瘤菌与豆科植物共生专一性认识。本文对近些年根瘤菌的系统学和分类学方面的研究做一综合分析和讨论,并展望未来的发展趋势。
1 根瘤菌的系统发育地位根瘤菌归属于原核生物细菌域、变形杆菌门下的α-变形杆菌纲和β-变形杆菌纲。在α-纲和β-纲内,分别包括了1个目,总共17个属(表 1)。查找某个属的分类纲要,可通过NCBI Taxonomy网站(http://www.ncbi.nlm.nih.gov/Taxonomy/)进行查找。在该网站上,输入关键词,如“Rhizobium”,就能查到它的分类谱系(Lineage)为:root (根)→cellular organisms (细胞生物)→Bacteria (细菌域)→Proteobacteria (变形杆菌门)→ Alphaproteobacteria (α-变形杆菌纲) →Rhizobiales (根瘤菌目)→Rhizobiaceae (根瘤菌科)→ Rhizobium/Agrobacterium group (根瘤菌属/土壤杆菌属类群)。
纲 Class | 目 Order | 科 Family | 属 Genus |
Alphaproteobacteria | Rhizobiales | Bradyrhizobiaceae | Bradyrhizobium |
Brucellaceae | Ochrobactrum | ||
Hyphomicrobiaceae | Devosia | ||
Methylobacteriaceae | Methylobacterium,Microvirga | ||
Phyllobacteriaceae | Mesorhizobium,Phyllobacterium | ||
Rhizobiaceae | Rhizobium/Agrobacterium,Shinella,Sinorhizobium/Ensifer,Neorhizobium,Pararhizobium | ||
Xanthobacteraceae | Azorhizobium | ||
Betaproteobacteria | Burkholderiales | Burkholderiaceae | Burkholderia,Cupriavidus |
α-变形杆菌纲的根瘤菌在世界范围内都有分布,且是与多属种豆科植物共生固氮的优势根瘤菌。归属于β-变形杆菌纲的伯克霍尔德氏菌属(Burkholderia)和贪铜菌属(Cupriavidus,之前称为Ralstonia)的根瘤菌,最初于2001年由Moulin等[1]及Chen等[2]提出,分离自热带地区的含羞草属植物,后来在菜豆属植物也发现有β-变形杆菌纲根瘤菌存在[3]。在γ-变形纲内曾有人报道Pseudomonas能结瘤固氮[4-5],但还未被广泛公认。2014以来,新提出的两个属为:新根瘤菌属(Neorhizobium)[6]和副根瘤菌属(Pararhizobium)[7],均属于α-变形杆菌纲。之前被归入到根瘤菌属的土壤杆菌属(Agrobacterium)和其它根瘤菌属(Allorhizobium)得以恢复[6]。
当前要确认某个菌株究竟属于哪个纲、哪个属的根瘤菌的分类依据16S rRNA基因的系统发育。获得菌株的16S rRNA基因序列后,在NCBI网站(http://blast.ncbi.nlm.nih.gov/)上进行序列比对,如果与某已知属的菌株的16S rRNA基因相似性大于及等于97%以上,可初步确定供试菌株与已知菌种为同一个属;如果低于97%,则可能代表一个新属[8]。
2 根瘤菌结瘤基因的有无及根瘤菌的范畴传统观念认为,根瘤菌要与豆科结瘤,必需具有结瘤基因,如nodABC基因等,以及能够合成结瘤因子[9]。但2007年Giraud等的研究却发现,有两株光合慢生根瘤菌BTAi1和ORS278,其基因组上没有常规的nodABC基因,也不能合成结瘤因子,但却能在一种热带豆科植物合萌属(Aeschynomene)上结瘤固氮[10]。这一现象目前没有在更多的根瘤菌及豆科植物中发现,而且大多数根瘤菌仍然是依赖于结瘤基因才能与豆科植物建立共生关系[11]。
没有结瘤基因但能结瘤的根瘤菌归属于根瘤菌的范畴是可以被理解的,而目前发现的一些细菌没有结瘤基因且不能结瘤,但它们在系统发育地位上,仍然可归属于根瘤菌的范畴。这是因为根瘤菌的分类学和系统学不以结瘤固氮基因作为定属、定种的标准,而只是作为共生基因型。属种的分类地位依据是16S rRNA基因及持家基因分析。没有结瘤基因的根瘤菌,多数是从土壤或根际土壤内分离得到的,如土壤根瘤菌(Rhizobium soli)[12]、恒河平原硫自养中慢生根瘤菌(Mesorhizobium thiogangeticum)[13]等;也有从根瘤内分离到的,但不能结瘤的根瘤菌,如中国农大根瘤菌(Rhizobium cauense)[14]等。
3 Ensifer和Sinorhizobium的属名之争及使用建议Ensifer和Sinorhizobium的中文译名分别为“剑菌属”和“中华根瘤菌属”[15]。Ensifer最初是从土壤内分离到的能捕食其它细菌的细菌,以Ensifer adhaerens为模式种,不能与豆科植物结瘤固氮,于1982年提出[16];Sinorhizobium由Chen等在1988年提出[17],分离自大豆的根瘤,是快生型的大豆根瘤菌,与之前研究者们发现的大豆慢生型根瘤菌不同,以Sinorhizobium fredii为模式种,能与大豆有效共生固氮,不能捕食其它细菌。2003年时,Willems等的研究发现,Sinorhizobium和Ensifer的不同菌种之间的16S rRNA基因序列相似性范围为97.9%−99.9%,因此,建议合并为一个属[18]。根据细菌命名法规,Ensifer发表的早,理应用此属名,但作者又认为:Sinorhizobium包含多个种,Ensifer只包含一个种;Sinorhizobium从字面上能体现出根瘤菌(rhizobium)的范畴;Sinorhizobium meliloti已广泛使用,是与苜蓿结瘤的高效固氮根瘤菌,据此,建议使用Sinorhizobium,而不用Ensifer[18]。但新西兰的学者J. M. Young不同意,从而引发Ensifer和Sinorhizobium两个属名之争。最终由国际细菌分类委员会裁定,认为Ensifer为合法属名,Sinorhizobium为同物异名,但Ensifer有优先权[19, 20]。目前在官方的分类学杂志——“International Journal of Systematic and Evolutionary Microbiology” (IJSEM)上,发表归属于这两个属的系统分育分支的新种时,必须使用Ensifer作为属名,如Ensifer soyae sp. nov.[21];而在其它非官方的非分类学杂志上,则根据自己的喜好,选择使用Ensifer或Sinorhizobium,但以使用Sinorhizobium的居多数[22, 23]。本文作者认为,如果所研究的菌株确实与Ensifer adhaerens关系十分接近,且不能结瘤固氮,那么就应该使用Ensifer这个属名;如果分离自根瘤内,且能与豆科植物结瘤固氮,又与Sinorhizobium内的一些种接近,那么最好使用Sinorhizobium这个属名。并且,为避免混淆,在研究与结瘤固氮相关的工作时,发表的非IJSEM杂志的论文时,则最好要使用Sinorhizobium,以使读者更容易理解论文内容,不产生歧义[18]。
4 ANI值取代DNA杂交70%的DNA同源性值曾是确定细菌之间是否为同一种的“黄金标准”之一。大于及等于70%的为同一个种,低于70%的则为不同种。传统上的DNA同源性分析常采用液相复性速率法、膜杂交法等。由于这些方法可重复性低,不同实验室的结果没法直接比较等弊端,现在已被快速发展的全基因组测序技术所取代。ANI是指Average nucleotide identity,即平均核苷酸一致性。两株细菌的全基因组的ANI比较值以94%为阈值,它相当于DNA同源性的70%这个值[24]。当ANI值大于等于95%时,为同一个菌种;反之,低于95%,则为不同的菌 种[25]。目前使用JSpecies软件计算ANI值是一个优选的方法[26],推荐使用。
虽然全基因组测定方便快速,但还是需要一定的时间分析和注释。对于一个未知菌株,想快速判断其是否是一个新种,可以采用只测定一些持家基因的方法。如Zhang等采用多基因序列(如SMc00019-truA-thrA)分析的ANI值比较,以94%或96%作为区分不同种的阈值标准,可以很好地将不同根瘤菌物种区分开来[27],此方法与其它持家基因及全基因组分析结果有很好地一致性。结合16S rRNA基因序列分析确定属的系统发育地位,以持家基因ANI比较,确定种间关系,就能确定未知根瘤菌的分类地位了。
5 表型和化学成分分析的快速化根瘤菌的表型分析,包括对不同碳源和氮源等的利用情况、抗生素抗性、不同盐浓度的适应性、不同生长温度耐受性,以及氧化酶、硝酸还原酶、BTB产酸或产碱判定等[28]。目前Biolog系统广泛应用于根瘤菌的生理生化表型测定,可以快速获得根瘤菌对不同碳源的利用情况,再补充其它的表型分析,即可得到根瘤菌的表型数据。虽然有学者指出,表型数据不稳定,重复性不是很好,建议用基因组数据代替表型分析,但就当前的情况看,根瘤菌的表型分析仍然是根瘤菌分类鉴定及发表新种的必需条件[29]。本文作者认为,表型数据对根瘤菌的分类鉴定仍然很重要,加上利用Biolog系统,可以短期内获得菌株的生理生化特征,要比测序、注释、比较等更快速。
根瘤菌的化学成分分析包括脂肪酸的测定、醌的分析和极性脂分析等[30]。它们的测定都有相应的快速分析方法,如气相色谱法、HPLC法、TCL法。在测定这些化学成分时,要保证供比较的多个菌株采用同样的培养方法、培养时间、培养基等,以使数据结果具有可比性。由于根瘤菌的醌在属及其以上水平有一致性,因此,醌的分析有时并非必需指标。
6 根瘤菌与宿主的共生专一性与混杂性根瘤菌与豆科宿主的共生关系通常分为以下两类:(1) 一种根瘤菌对应于一种或一个互接种族的植物,如Sinorhizobium meliloti只能与苜蓿、草木樨、葫芦巴结瘤固氮,不能与其它植物结瘤固 氮[31, 32],这称为共生专一性;(2) 一种根瘤菌可以与多属种豆科植物结瘤固氮,如Sinorhizobium sp. NGR234,它可与120个属的豆科植物结瘤固氮[33],这称为根瘤菌的宿主混杂性。
最新的研究表明[34],自然条件下,与药用豆科植物苦参(Sophora flavescens)共生的根瘤菌可归属到5个属中;实验室条件下的交叉结瘤研究进一步发现,苦参能与来自α-变形杆菌纲和β-变形杆菌纲的超过30个种的根瘤菌结瘤固氮。特别是相互之间不能在对方的宿主上交叉结瘤的根瘤菌也都能与苦参结瘤固氮,如来自苜蓿、大豆、花生、紫云英、菜豆、豌豆的根瘤菌,除个别菌种外,都能与苦参结瘤固氮[34]。这一研究结果说明,苦参是极其混杂性的豆科宿主植物,能与很多根瘤菌共生固氮,我们称之为豆科植物的根瘤菌混杂性。这一发现彻底打破了根瘤菌与豆科植物之间的共生专一性认识,这项研究将促进人们重新考虑并研究豆科植物与根瘤菌之间的分子识别与信号交流的关系。
7 总结与展望本文对2001年β-变形杆菌纲的根瘤菌发现以来的根瘤菌系统学、根瘤菌分类学上的争议,以及根瘤菌与豆科植物关系等方面的研究做了总结和讨论。基因组测序技术和根瘤菌表型特征、化学特征的快速测定,加速了根瘤菌的分类鉴定和系统学研究。
未来根瘤菌的系统学研究中,将向更多的豆科宿主结瘤调查发展。国内的根瘤菌系统学研究中,分离自花生的根瘤菌新种将不断涌现(个人通讯)。而国外的根瘤菌系统学研究中,则偏重于从没有调查过的豆科植物中发现新的根瘤菌物种,特别是从非洲地区的豆科植物根瘤中发现新的根瘤菌物种。大量的新种仍有待描述,而这仍会导致建立一些新属,或在原本不包括根瘤菌的细菌属内分离到一些共生固氮的种或菌株。对于苦参这种混杂性宿主,还有可能从不同的生态区发现新的根瘤菌物种。目前,国内报道的分离自苦参的根瘤菌归属于根瘤菌属[35]和叶杆菌属[36],国外报道了分离自与苦参同属的槐属(Sophora)中的其它树种根瘤分离到的根瘤菌新种[37],槐属植物的根瘤菌还可能有更多的新种出现。
除此之外,自土壤中分离到的不能在豆科植物上结瘤固氮的,但系统发育地位在根瘤菌范畴之内的根瘤菌新种将会不断出现,这将丰富根瘤菌的系统学研究,也有助于人们更好地认识共生和非共生的根瘤菌系统发育、物种进化及与豆科宿主的关系。再者,苦参这种混杂性极高的植物,将成为豆科植物中的“NGR234”,激发人们去研究为什么苦参这种植物对根瘤菌选择性这么低,是不是有其独特的结瘤因子受体结构?还有,一个正在越来越热的研究领域——根瘤菌的试验进化[38, 39],结合前面提到的全基因组测序比较、表型快速分析技术,以及最近发展起来的甲基化组测定等手段,将有助于人们认识根瘤菌是如何进化的,进化的方向是什么,结瘤固氮基因是如何起源及进化的等科学问题。再结合大数据分析,开源协作分析系统(如R语言[40]及BioConductor[41]等)的应用,根瘤菌的系统学和进化学正向着一个崭新的方向前进!
[1] | Moulin L, Munive A, Dreyfus B, et al. Nodulation of legumes by members of the β-subclass of Proteobacteria[J]. Nature, 2001, 411 (6840): 948–950. |
[2] | Chen WM, Laevens S, Lee TM, et al. Ralstonia taiwanensis sp. nov., isolated from root nodules of Mimosa species and sputum of a cystic fibrosis patient[J]. International Journal of Systematic and Evolutionary Microbiology, 2001, 51 (5): 1729–1735. |
[3] | Talbi C, Delgado MJ, Girard L, et al. Burkholderia phymatum strains capable of nodulating Phaseolus vulgaris are present in Moroccan soils[J]. Applied and Environmental Microbiology, 2010, 76 (13): 4587–4591. |
[4] | Benhizia Y, Benhizia H, Benguedouar A, et al. Gamma proteobacteria can nodulate legumes of the genus Hedysarum[J]. Systematic and Applied Microbiology, 2004, 27 (4): 462–468. |
[5] | Shiraishi A, Matsushita N, Hougetsu T. Nodulation in black locust by the Gammaproteobacteria Pseudomonas sp. and the BetaproteobacteriaBurkholderia sp.[J]. Systematic and Applied Microbiology, 2010, 33 (5): 269–274. |
[6] | Mousavi SA, Österman J, Wahlberg N, et al. Phylogeny of the Rhizobium-Allorhizobium-Agrobacterium clade supports the delineation of Neorhizobium gen. nov.[J]. Systematic and Applied Microbiology, 2014, 37 (3): 208–215. |
[7] | Mousavi SA, Willems A, Nesme X, et al. Revised phylogeny of Rhizobiaceae: proposal of the delineation of Pararhizobium gen. nov., and 13 new species combinations[J]. Systematic and Applied Microbiology, 2015, 38 (2): 84–90. |
[8] | Stackebrandt E, Goebel BM. Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology[J]. International Journal of Systematic Bacteriology, 1994, 44 (4): 846–849. |
[9] | Relić B, Perret X, Estrada-García MT, et al. Nod factors of Rhizobium are a key to the legume door[J]. Molecular Microbiology, 1994, 13 (1): 171–179. |
[10] | Giraud E, Moulin L, Vallenet D, et al. Legumes symbioses: absence of nod genes in photosynthetic bradyrhizobia[J]. Science, 2007, 316 (5829): 1307–1312. |
[11] | D’Haeze W, Holsters M. Nod factor structures, responses, and perception during initiation of nodule development[J]. Glycobiology, 2002, 12 (6): 79R–105R. |
[12] | Yoon JH, Kang SJ, Yi HS, et al. Rhizobium soli sp. nov., isolated from soil[J]. International Journal of Systematic and Evolutionary Microbiology, 2010, 60 (Pt 6): 1387–1393. |
[13] | Ghosh W, Roy P. Mesorhizobium thiogangeticum sp. nov., a novel sulfur-oxidizing chemolithoautotroph from rhizosphere soil of an Indian tropical leguminous plant[J]. International Journal of Systematic and Evolutionary Microbiology, 2006, 56 (1): 91–97. |
[14] | Liu TY, Li Jr Y, Liu XX, et al. Rhizobium cauense sp. nov., isolated from root nodules of the herbaceous legume Kummerowia stipulacea grown in campus lawn soil[J]. Systematic and Applied Microbiology, 2012, 35 (7): 415–420. |
[15] | Yang RF, Tao TS, et al. Dictionary of Bacterial Names with English Explanation and Chinese Translation[M]. Beijing: Chemical Industry Press, 2011 . (in chinese) 杨瑞馥, 陶天申, 等. 细菌名称双解及分类词典[M]. 北京: 化学工业出版社, 2011 . |
[16] | Casida Jr LE. Ensifer adhaerens gen. nov., sp. nov.: a bacterial predator of bacteria in soil[J]. International Journal of Systematic Bacteriology, 1982, 32 (3): 339–345. |
[17] | Chen WX, Yan GH, Li JL. Numerical taxonomic study of fast-growing soybean rhizobia and a proposal that Rhizobium fredii be assigned to Sinorhizobium gen. nov[J]. International Journal of Systematic Bacteriology, 1988, 38 (4): 392–397. |
[18] | Willems A, Fernández-López M, Muoz-Adelantado E, et al. Description of new Ensifer strains from nodules and proposal to transfer Ensifer adhaerens Casida 1982 to Sinorhizobium as Sinorhizobium adhaerens comb. nov. Request for an opinion[J]. International Journal of Systematic and Evolutionary Microbiology, 2003, 53 (4): 1207–1217. |
[19] | Young JM. Sinorhizobium versus Ensifer: may a taxonomy subcommittee of the ICSP contradict the Judicial Commission[J]. International Journal of Systematic and Evolutionary Microbiology, 2010, 60 (7): 1711–1713. |
[20] | Judicial Commission of the International Committee on Systematics of Prokaryotes. The genus name Sinorhizobium Chen, et al. 1988 is a later synonym of Ensifer Casida 1982 and is not conserved over the latter genus name, and the species name ‘Sinorhizobium adhaerens’ is not validly published. Opinion 84[J]. International Journal of Systematic and Evolutionary Microbiology, 2008, 58 (8): 1973. |
[21] | Li QQ, Wang ET, Chang YL, et al. Ensifer sojae sp. nov., isolated from root nodules of Glycine max grown in saline-alkaline soils[J]. International Journal of Systematic and Evolutionary Microbiology, 2011, 61 (8): 1981–1988. |
[22] | Blanquet P, Silva L, Catrice O, et al. Sinorhizobium meliloti controls nitric oxide-mediated post-translational modification of a Medicago truncatula nodule protein[J]. Molecular Plant-Microbe Interactions, 2015, 28 (12): 1353–1363. |
[23] | Ghnaya T, Mnassri M, Ghabriche R, et al. Nodulation by Sinorhizobium meliloti originated from a mining soil alleviates Cd toxicity and increases Cd-phytoextraction in Medicago sativa L.[J]. Front Plant Science, 2015, 6 : 863. |
[24] | Goris J, Konstantinidis KT, Klappenbach JA, et al. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities[J]. International Journal of Systematic and Evolutionary Microbiology, 2007, 57 : 81–91. |
[25] | Arahal DR. Whole-genome analyses: average nucleotide identity[J]. Methods in Microbiology, 2014, 41 : 103–122. |
[26] | Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106 (45): 19126–19131. |
[27] | Zhang YM, Tian CF, Sui XH, et al. Robust markers reflecting phylogeny and taxonomy of rhizobia[J]. PLoS One, 2012, 7 (9): e44936. |
[28] | Gao JL, Sun JG, Li Y, et al. Numerical taxonomy and DNA relatedness of tropical rhizobia isolated from Hainan Province, China[J]. International Journal of Systematic Bacteriology, 1994, 44 (1): 151–158. |
[29] | Ormeo-Orrillo E, Martínez-Romero E. Phenotypic tests in Rhizobium species description: an opinion and (a sympatric speciation) hypothesis[J]. Systematic and Applied Microbiology, 2013, 36 (3): 145–147. |
[30] | Zhang JJ, Liu TY, Chen WF, et al. Mesorhizobium muleiense sp. nov., nodulating with Cicer arietinum L.[J]. International Journal of Systematic and Evolutionary Microbiology, 2012, 62 (Pt 11): 2737–2742. |
[31] | Yan AM, Wang ET, Kan FL, et al. Sinorhizobium meliloti associated with Medicago sativa and Melilotus spp. in arid saline soils in Xinjiang, China[J]. International Journal of Systematic and Evolutionary Microbiology, 2000, 50 (5): 1887–1891. |
[32] | He YR, Wang JY, Wang ET, et al. Trigonella arcuata-associated rhizobia-an Ensifer (Sinorhizobium) meliloti population adapted to a desert environment[J]. Plant and Soil, 2011, 345 (1/2): 89–102. |
[33] | Pueppke SG, Broughton WJ. Rhizobium sp. strain NGR234 and R. fredii USDA257 share exceptionally broad, nested host ranges[J]. Molecular Plant-Microbe Interactions, 1999, 12 (4): 293–318. |
[34] | Jiao YS, Liu YH, Yan H, et al. Rhizobial diversity and nodulation characteristics of the extremely promiscuous legume Sophora flavescens[J]. Molecular Plant-Microbe Interactions, 2015, 28 (12): 1338–1352. |
[35] | Jiao YS, Yan H, Ji ZJ, et al. Rhizobium sophorae sp. nov. and Rhizobium sophoriradicis sp. nov., nitrogen-fixing rhizobial symbionts of the medicinal legume Sophora flavescens[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65 (Pt 2): 497–503. |
[36] | Jiao YS, Yan H, Ji ZJ, et al. Phyllobacterium sophorae sp. nov., a symbiotic bacterium isolated from root nodules of Sophora flavescens[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65 (Pt 2): 399–406. |
[37] | de Meyer SE, Tan HW, Heenan PB, et al. Mesorhizobium waimense sp. nov. isolated from Sophora longicarinata root nodules and Mesorhizobium cantuariense sp. nov. isolated from Sophora microphylla root nodules[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65 (10): 3419–3426. |
[38] | Marchetti M, Capela D, Glew M, et al. Experimental evolution of a plant pathogen into a legume symbiont[J]. PLoS Biology, 2010, 8 (1): e1000280. |
[39] | Sachs JL, Russell JE, Hollowell AC. Evolutionary instability of symbiotic function in Bradyrhizobium japonicum[J]. PloS One, 2011, 6 (11): e26370. |
[40] | Ihaka R, Gentleman R. R: a language for data analysis and graphics[J]. Journal of Computational and Graphical Statistics, 1996, 5 (3): 299–314. |
[41] | Gentleman R, Carey VJ, Huber W, et al. Bioinformatics and Computational Biology Solutions Using R and Bioconductor[M]. New York: Springer, 2005 . |