神舟十號飛船下行菌株Microbacterium sp．LCT?H2的基因組測序及分析

周 宏，謝 瓊，劉長庭

（1.中國人民解放軍總醫院南樓呼吸科，北京100853；2.泰山醫學院，泰安271000；3.中國航天員訓練中心醫監醫保研究室，北京100094）

·基礎研究·

神舟十號飛船下行菌株Microbacterium sp．LCT?H2的基因組測序及分析

周 宏1，2，謝 瓊3，劉長庭1?

（1.中國人民解放軍總醫院南樓呼吸科，北京100853；2.泰山醫學院，泰安271000；3.中國航天員訓練中心醫監醫保研究室，北京100094）

神舟十號飛船返回艙內表面分離得到一株下行菌株LCT?H2，顯微觀察和革蘭氏染色結果為短桿狀的革蘭氏陽性菌，利用葡萄糖、纖維二糖等多種碳源，對該菌進行16S rRNA測序，鑒定其為微桿菌屬成員，并命名為Microbacterium sp．LCT?H2。對該下行菌進行了全基因組測序和基因功能注釋分析，其基因組大小為3.36 M，包括3235個編碼序列和53個RNA基因，COG數據庫比對共注釋到21組COG功能分類。可對研究密閉飛行器內微生物種類分布等提供一定幫助。

微桿菌；太空飛行；全基因組測序；功能分類

1 引言

隨著航天技術的發展，人類的外太空探索活動越來越頻繁。盡管載人飛船及內部物品等經過嚴格的真空消毒，載人航天器密閉環境條件仍然容易滋生細菌和真菌等微生物，這些微生物能在密閉艙室內的金屬、高分子復合材料等表面形成生物膜［1］，它們的生長繁殖和代謝會對材料產生腐蝕，嚴重威脅空間站的長期在軌運行安全，縮短空間站的服役時間［2］。研究密閉航天器內的微生物種類對于航天器內微生物的安全防護具有重要的意義。

早期人們對許多環境微生物的認識僅來自于對16S rRNA的研究［3］，一些重要的遺傳信息很難獲得。現在越來越多的新型培養方法和分子生物學方法，如基因組學、蛋白質組學、宏基因組學等，被用來研究新獲得微生物的潛在特征和功能等［4］。

本研究圍繞一株分離自神舟十號載人飛船返回艙內冷凝水的細菌Microbacterium sp．LCT?H2展開。通過16S rRNA測序鑒定該菌菌屬，完成菌株基因組的高通量測序、基因注釋分析和功能分類，研究該菌的表型及生理生化特征。

2 材料和方法

2.1 菌株與培養

菌株Microbacterium sp．LCT?H2為中國航天員訓練中心自神舟飛船的返回艙內冷凝水分離得到。菌株分別于37℃培養箱或搖床中進行靜置或振蕩培養，液體培養基為腦心浸出液培養基（BHI，Oxiod），配方為37 g／L BHI；固體培養基為MH平板。

2.2 基因組DNA的提取

Microbacterium sp．LCT?H2接種至BHI液體培養基，37℃，180 rpm過夜培養，6000 rpm離心5 min收集菌體，棄上清。利用Qiagen公司的基因組DNA提取試劑盒（Code No.51304，詳細步驟見產品說明書）提取收集的菌體的基因組DNA。

2.3 菌種鑒定及系統發育樹構建

以Microbacterium sp．LCT?H2的基因組DNA為模板，利用細菌16SrRNA通用測序引物27F（5’?AGAGTTTGATCCTGGCTCAG?3’）和1492R（5’?TACGGCTACCTTGTTACGACTT?3’）擴增其16S rRNA片段并測序，得到其堿基序列。用CLUSTALW［5］進行多序列比對，并用MEGA 5.0構建NJ系統發育樹［6］。

參照Chen等[21]的方法，采用透析測定樣品的釋放情況。取10 mL樣品于透析袋(分子截留量8 000～14 000 Da)中，密封好后浸沒于200 mL的pH 7.0的PBST緩沖溶液中(吐溫-20，0.5% v/v)。游離的丁香酚可透過透析袋進入到緩沖溶液中，而酪蛋白酸鈉穩定微乳則被截留在透析袋中。在預定時間內，取出適量透析液，用相同的緩沖溶液稀釋至適當濃度，280 nm下測定其吸光度并計算丁香酚釋放的量，繪制釋放曲線。

2.4 菌株性態及生理特性

2.4.1 菌體形態觀察

采用革蘭氏染色試劑盒（Code No.G1060，Solarbio）對LCT?H2菌體進行染色。BHI液體培養基過夜培養菌體，將菌液滴于載玻片中央，在酒精燈上略加溫，使之迅速干燥。結晶紫、沙黃染色液對干燥菌體一次染色脫色，詳細步驟見產品說明書。染色完成并晾干后，置于光學顯微鏡下觀察染色情況。

分別在透射電子顯微鏡（TEM）和掃描電子顯微鏡（SEM）下觀察菌體形態。1 ml培養至穩定前期的菌液，6000 rpm離心5 min收集菌體，棄上清；PBS緩沖液洗滌一次后將菌體懸浮至合適濃度。菌懸液等體積與4%多聚甲醛＋0.5%戊二醛溶液混合，滴在有支持膜的銅網上，15min后用濾紙吸取未吸附樣品；吸取磷鎢酸染液小心滴在銅網上，染色1 min后用濾紙吸取未吸附的樣品；樣品干燥后在Philips EM 400T型透射電子顯微鏡（TEM）下觀察菌體形態。首先用2.5%戊二醛、PBS緩沖液以及1%鋨酸清洗菌體，乙醇梯度脫水并乙酸異戊酯置換，樣品干燥后做噴金處理［7］，處理好的樣品置于FEIQuanta 200型掃描電鏡下觀察。

2.4.2 BIOLOG板分析

Biolog微孔板（Biolog GENIIIMicroPlate，Bi?olog）用來對LCT?H2菌株進行94種表型檢測，其中包括71種碳源利用率分析和23種化學敏感性分析。Microbacterium sp．LCT?H2接種至BHI液體培養基，37℃，180 rpm過夜培養，6000 rpm離心5 min收集菌體，棄上清。將菌體用IF?B接種液懸浮并調節至適當濃度（OD600約為0.8），分別取100μL接種液加入Biolog微孔板的96個孔中，將微孔板放置在37℃培養箱中培養24 h后，用酶標儀測定590 nm各孔的吸光度值（OD590）［8］。

2.5 基因組測序、組裝及注釋

基因組DNA樣品委托美吉生物有限公司進行全基因組序列測定、基因預測、基因功能注釋、非編碼RNA預測等工作，序列分析工作委托美吉生物有限公司和上海瀚宇生物科技有限公司共同完成。

基因組DNA提取后進行質量鑒定，在濃度和純度達到測序要求后進行全基因組測序。利用Covaris進行基因組DNA片段化，構建插入片段約300～500 bp的基因組測序文庫，通過橋式PCR擴增得到DNA簇，然后采用Illumina Hiseq 2000測序技術完成該菌株的基因組掃描測序，獲得約1 Gb的原始測序數據。Illumina HiSeq2000測序得到的原始圖像數據經過Base Calling轉化為序列數據，結果以FASTQ文件格式來存儲，包含測序read的序列信息以及測序質量信息［9］；原始數據經過預處理去除接頭、引物及低質量數據后，利用SOAPdenovo［10］拼接軟件對優化序列進行多個Kmer參數的拼接，得到最優的組裝結果，并運用Gap Closer軟件對組裝結果進行局部內洞填充和堿基校正；分別利用RNAmmer和tRNAscan?SE軟件對基因組中包含的rRNA和tRNA進行預測［11］；利用Glimmer 3.0軟件進行基因預測［12?14］；利用各種數據庫進行基因功能注釋和分類。

3 結果與分析

菌株LCT?H2可在MH固體平板上形成無色、表面光滑且邊緣整齊菌落。革蘭氏染色后在光鏡下觀察，發現菌體被染成紫色，為革蘭氏陽性菌特征（圖1（a））。在透射電鏡及掃描電鏡下觀察發現，菌體呈桿狀，分散排布，表面沒有鞭毛、菌毛等結構（圖1（b）、（c））。

3.2 Microbacterium sp．LCT?H2菌株鑒定

利用通用引物測序得到的16S rRNA基因序列，構建多株細菌的系統發育樹（圖2），顯示與Microbacterium屬多個菌株位于一個大的系統分支，且與Microbacterium paraoxydans進化關系最近。Biolog微孔板結果表明LCT?H2菌株可利用多種碳源，包括單糖（葡萄糖、果糖、甘露糖等）、二糖（纖維二糖、龍膽二糖、蔗糖等），多糖（果膠、糊精等）以及多種無機碳源等，對醋竹桃霉素（Troleandomycin）、硫酸四癸鈉（Niaproof 4）、二甲胺四環素（Minocycline）等敏感（圖3）。綜合該菌在生理生化特征、16S rRNA基因同源性和化學指標等方面的實驗結果，確定菌株LCT?H2是微桿菌屬，命名為Microbacterium sp．LCT?H2。

3.3 基因組測序與分析

3.3.1 基因組序列基本信息

通過Illumin Hiseq2000平臺測序，基于測序數據組裝得到LCT?H2基因組大小為3 356 231 bp，GC含量為70.5%，共30個scaffold，46個contig。基因組組分分析發現，LCT?H2預測出3235個編碼基因，總長度為3 084 376 bp，占基因組全長的91.9%。tRNA 50個，rRNA 3個。LCT?H2菌株的基因組序列具體信息詳見表1。

3.3.2 基因組序列基本信息

用Glimmer 3.0軟件對LCT?H2基因組的ORF進行預測，獲得3235個編碼基因，依據COG分類標準將1575個基因劃分為21類，以英文大寫字母表示每一類的代碼，每一類基因的數量和功能見表2：RNA形成及修飾（A），共1個基因；能量產生和轉化（C），共96個基因；細胞周期調控、細胞分裂、染色體分割（D），共15個基因；氨基酸運輸代謝（E），共204個基因；核苷酸運輸代謝（F），共57個基因；糖類運輸代謝（G），共135個基因；輔酶運輸代謝（H），共56個基因；脂質運輸代謝（I），共66個基因；翻譯、核糖體結構和生物合成（J），共121個基因；轉錄（K），共108個基因；復制、重組、修復（L），共80個基因；細胞壁、細胞膜合成（M），共59個基因；細胞運動（N），共11個基因；蛋白質修飾和折疊（O），共48個基因；無機離子運輸代謝（P），共134個基因；次生產物合成、轉運和分解（Q），共39個基因；一般功能預測（R），175個基因；功能未知（S），77個基因；信號轉導機制（T），共48個基因；細胞內轉運，分泌和膜泡運輸（U），共23個基因；防御機制（V），共22個基因；還有1660個未能找到對應的COG分類。

表1 M icrobacterium sp．LCT?H2菌株基因組的核苷酸含量和基因計數水平Table 1 Nucleotide content and gene count levels of the genome

4 討論

空間環境具有地面環境不存在的一些特殊環境因素，如微重力、高真空、離子輻射、電磁輻射和極度溫差等，這些因素除了可以誘導空間環境微生物發生基因水平的改變，進而影響微生物的生物學性狀和功能外，還會影響航天器材料、艙室環境等。

微桿菌屬是一類廣泛存在于土壤、昆蟲、海洋、病患、植物等多種環境中的微生物［15］，曾有報道微桿菌屬菌種可以引起人類、動物和植物疾病［16?18］。在神舟十號船返回艙內分離得的微桿菌屬Microbacterium sp．LCT?H2，參與氨基酸運輸代謝，糖類運輸代謝，翻譯、核糖體結構和生物合成等相關功能。

表2 LCT?H2基因組21個COG分類的基因數統計Table 2 Number of genes associated w ith the 21 general COG functional categories

5 結論

1）LCT?H2菌株經本研究鑒定為微桿菌屬，并命名為Microbacterium sp．LCT?H2，為短桿狀的革蘭氏陽性菌。

2）Microbacterium sp．LCT?H2基因組大小為3.36 M，包括3235個編碼序列和53個RNA基因，COG數據庫比對共注釋到21組COG功能分類。

（References）

［1］ Kim W，Tengra F K，Young Z，et al.Spaceflight promotes biofilm formation by Pseudomonasaeruginosa［J］.Plos One，2013，8（4）：e62437 1?8.

［2］ Novikova N D.Review of the knowledge ofmicrobial contami?nation of the Russianmanned spacecraft［J］.Microbial Ecol?ogy，2004，47（2）：127?132.

［3］ Handelsman J.Metagenomics：application of genomics to un?culturedmicroorganisms［J］.Microbiology andmolecular biol?ogy reviews，2004，68（4）：669?685.

［4］ Hemme C L，Deng Y，Gentry T J，et al.Metagenomic in?sights into evolution of a heavy metal?contaminated groundwa?ter microbial community［J］.The ISME journal，2010，4（5）：660?672.

［5］ CLUSTALW［EB／OL］.http：／／www.ebi.ac.uk／Tools／msa／clustalw2／.

［6］ Tamura K，Peterson D，Peterson N，et al.MEGA5：molecu?lar evolutionary genetics analysis using maximum likelihood，evolutionary distance，and maximum parsimony methods［J］. Molecular Biology and Evolution，2011，28（10）：2731?2739.

［7］ 謝家儀，董光軍，劉振英.掃描電鏡的微生物樣品制備方法［J］.電子顯微學報，2005，24（4）：440?440. Xie Jiangyi，Dong Guangjun，Liu Zhenying.Methodsof pre?paringmicrobial samplesfor scanning electron microscope［J］. Journal of Chinese Electron Microscopy Society，2005，24（4）：440?440.（in Chinese）

［8］ Li J，Liu F，Wang Q，et al.Genomic and transcriptomic a?nalysis of NDM?1 Klebsiella pneumoniae in spaceflight reveal mechanisms underlying environmental adaptability［J］.Scien?tific Reports，2014（4）：6216.

［9］ Dolan P C，Denver D R.TileQC：a system for tile?based quality control of Solexa data［J］.BMC bioinformatics，2008， 9（1）：1?10.

［10］ SOAPdenovo［EB／OL］.http：／／soap.genomics.org.cn／.

［11］ Li R，Zhu H，Ruan J，et al.De novo assembly of human ge?nomeswith massively parallel short read sequencing［J］.Ge?nome Research，2010，20（2）：265?272.

［12］ Glimmer 3.0［EB／OL］.http：／／www.cbcb.umd.edu／soft?ware／glimmer／.

［13］ Delcher A L，Bratke K A，Powers E C，et al.Identifying bacterial genes and endosymbiont DNA with Glimmer［J］. Bioinformatics，2007，23（6）：673?679.

［14］ Lukashin A V，Borodovsky M.GeneMark.hmm：new solu?tions for gene finding［J］.Nucleic Acids Research，1998，26（4）：1107?1115.

［15］ KumariP，Bandyopadhyay S，Das SK.Microbacterium oryzae sp.nov.，an actinobacterium isolated from rice field soil［J］. International Journal of Systematic and Evolutionary Microbiol?ogy，2013，63（7）：2442?2449.

［16］ Alonso?Echanove J，Shah SS，Valenti A J，etal.Nosocomial outbreak of Microbacterium species bacteremia among cancer patients［J］.Journal of Infectious Diseases，2001，184（6）：754?760.

［17］ Hodgkin J，Kuwabara PE，Corneliussen B.A novel bacterial pathogen，Microbacterium nematophilum，inducesmorphologi?cal change in the nematode C.elegans［J］.Current Biology，2000，10（24）：1615?1618.

［18］ Hisatoshi K，Subandiyah S，Ochiai H.Red stripe of rice is caused by a bacterium Microbacterium sp［J］.Journal of Gen?eral Plant Pathology，2000，66（2）：149?152.

Draft Genome Sequencing and Analysis of Strain M icrobacterium sp．LCT?H2 Isolated from ShenzhouⅩReturn Capsule

ZHOU Hong1，2，XIE Qiong3，LIU Changting1?
（1.Nanlou Respiratory Disease Department，The 301 Hospital，Beijing 100853，China；2.Taishan Medical College，Taian 271000，China；3.Department of Medical Monitoring and Medical Support，China Astronaut Research and Training Center，Beijing 100094，China）

Microbacterium sp．LCT?H2 strain was a short?rod，Gram?negative strain isolated from in?side of the ShenzhouⅩreturn capsule，which could utilizemany carbon sources including glucose and cellobiose.LCT?H2 was identified as amember of the Microbacterium according to the sequence of 16S rRNA and the draftgenome sequence.The total size of the genomewas3，356，231bp with 3，235 protein?coding and 50 RNA genes.21 COG functional categorieswere annotated based on COG database.The resultsmay be helpful to the study ofmicrobial species distribution inside the space?craft.

Microbacterium；space flight；whole?genome sequencing；functional category

R856

A

1674?5825（2017）01?0109?05

2015?12?02；

2016?12?21

國家重點基礎研究發展計劃（973）（2014CB744400）；載人航天預先研究項目（040203）；全軍醫學科研“十二五”課題重點項目（BWS12J046）

周宏，女，博士，講師，在站博士后，研究方向為微生物生化與分子生物學。E?mail：hongzhoujy＠hotmail.com

?通訊作者：劉長庭，男，碩士，教授，研究方向為空間生命科學。E?mail：changtingliu＠sohu.com