玉米核不育基因 MS22 克隆及其同源基因生物信息學分析

高永鋼

摘要 利用水稻 ROXY1、ROXY2基因序列，在玉米基因組數據庫進行同源檢索，發現 GRMZM2G442791基因（ MS22）序列與 ROXY1、 ROXY2之間具有較高的相似度（73%～79%）。以玉米 B73自交系為實驗材料克隆 MS22啟動子區序列;用生物信息學方法對水稻 ROXY1、 ROXY2及 MS22基因核苷酸序列以及蛋白質結構、親水性和疏水性及進化樹等進行分析。結果表明，玉米雄性核不育相關基因 MS22編碼區長度為480 bp，編碼159個氨基酸; ROY1編碼區長度為 408 bp，編碼136個氨基酸; ROXY2編碼區長度為420 bp，編碼140個氨基酸。通過對蛋白功能結構域進行分析發現， MS22、 ROXY1、 ROXY2都為谷氧還蛋白家族成員，具有 GRX 家族基因相似的功能結構域和近似的三級結構。因此推測， MS22可能與 ROXY1、 ROXY2具有相似的功能，可能參與氧化還原調節、植物脅迫防御，在控制雄性配子發生中起主要調控作用。

關鍵詞 玉米;雄性不育; MS22基因;生物信息

中圖分類號 S513???? 文獻標識碼 A????? DOI：10.12008/j.issn.1009-2196.2022.03.009

Cloning of Maize Genic Male Sterility Gene MS22 and Bioinformatics Analysis of Its Homologous Gene

GAO Yonggang1，2

（1. Inner Mongolia University of Science and Technology， Baotou， Inner Mongolia 014010， China;2. College of Energy and Envi-ronment， Inner Mongolia University of Science and Technology， Baotou， Inner Mongolia 014010， China）

Abstract? Using rice ROXY1 and ROXY2 gene sequences， a homology search was conducted in the maize genome database. the searching results showed that the GRMZM2G442791（MS22） gene sequence had a high degree of similarity with ROXY1 and ROXY2（73%-79%）. Primers were designed by using Primer5.0， and MS22 was cloned by using the maize B73 inbred line. Bioinformatics methods were used to analyze the nucleotide sequences， protein structure， hydrophilia， hydrophobicity and evolution tree of rice ROXY1， ROXY2 and MS22. The results showed that the coding regions of the maize male sterility related gene MS22， ROY1 and ROXY2 are 480 bp， 408 bp and 420 bp in length， respectively and encode 159 amino acids， 136 amino acids and 140 amino acids， respectively. The analysis of protein functions and domains showed that MS22， ROXY1 and ROXY2 are members of the glutaredoxin （GRX） family and have the similar functional domains and similar tertiary structure to GRX family genes. Therefore， it is speculated that MS22 may have similar functions to ROXY1 and ROXY2， could participate in redox regulation， plant stress defense， and play a major role in controlling male gametogenesis.

Keywords? maize; male sterility; MS22; bioinformatics

玉米是我國最重要的糧食作物之一，其產量的提高主要依賴于雜種優勢的利用。玉米核雄性不育（ Genic male sterility）是由細胞核基因突變導致的雄性敗育現象，表現為植物有性生殖過程中不能產生正常的花藥、花粉或雄配子而雌蕊功能正常的遺傳現象[1]。玉米雄性不育基因的克隆不僅有助于解析玉米花藥發育的分子機理，同時也為玉米雜交種的生產提供理論支撐。玉米目前是國際上第一大糧食作物，玉米商業雜交種質量直接影響玉米的產量和品質。細胞質雄性不育系在抗病性、育性遺傳等方面存在諸多不確定性和干擾因素。近些年有關學者將目光再次集中于受環境影響較小、遺傳較為穩定且持續的細胞核雄性不育[2-4]。

迄今，國內外已經發現了200多個玉米雄性不育突變體，大部分為隱性核不育突變體，其中，部分完成了染色體的初步定位和細胞學分析，僅有少數基因完成了分離克隆。目前已經克隆的玉米隱性核不育基因有12個，包括 ms1、ms7、ms8、 ms9、 ms22/msca1、 ms23、 ms26、 ocl4、 mac1、 ms30、 ms32和 ms45[5-6]。

谷氧還蛋白（ GRX）是參與各種細胞過程的小型氧化還原酶，并在氧化應激反應中發揮重要作用。前人報道稱，2個 CC 型 GRX 基因 ROXY1和 ROXY2一起控制擬南芥中的花藥發育[7]。 ROXY1和 ROXY2是 CC 型花谷氧還蛋白，在擬南芥花藥發育中具有調控功能。 ROXY1在高等植物中編碼谷氧還蛋白亞基，這類蛋白調節植物花器官發育， ROXY1蛋白中保守的谷胱甘肽結合甘氨酸的誘變，表明 CC 型 GRX 需要與谷胱甘肽相互作用以催化必要的生物合成反應[8-9]。 roxy1和? roxy2的單突變體可育并產生正常的花藥;然而， roxy1、 roxy2雙突變體是不育的，不產生花粉。研究表明， ROXY1和 ROXY2在花藥發育過程中以復合模式表達， ROXY1和 ROXY2功能的缺乏在轉錄水平影響大量花藥基因，控制絨氈層發育。對 ROXY1和 ROXY2的分析表明，該類基因可能參與氧化還原調節、植物脅迫防御，在控制雄性配子發生中起主要作用[10-11]。

本研究在玉米基因組數據庫（ https ：//maizegdb.org/）中對水稻 ROXY1和 ROXY2基因序列進行同源比對、檢索，獲得玉米 GRMZM2G442791基因序列;經基因克隆、測序、序列分析后，運用生物信息學手段對該基因及 ROXY1和 ROXY2進行分析。

1 材料與方法

本實驗所用材料為玉米 B73 自交系，由實驗室盆栽種植，待玉米生長至 V3時期，取幼嫩葉片提取基因組 DNA，用于 MS22基因的擴增。

通過 NCBI 檢索 ROXY1和 ROXY2基因序列，在玉米基因組數據庫進行序列同源比對。利用expasy數據庫（http：//www.expasy.org /tools / protparam.html）對基因的理化性質、結構和功能域進行序列分析預測。通過InterPro（ http：//www.ebi.ac.uk/interpro/）對 ROXY1和 ROXY2全長 cDNA 序列編碼的氨基酸序列和玉米 MS22全長 cDNA 序列編碼的氨基酸序列進行結構分析。用 SOPMA（http：//npsa- pbil.ibcp.fr/cgi- bin/=/ NPSA/npsa_sopma.html ）對 MS22、 ROXY1 和 ROXY2基因編碼蛋白質產物的二級結構進行分析。利用Phyre（http：//www.sbg.bio.ic.ac.uk/ phyre/）程序對 MS22、 ROXY1和 ROXY2基因編碼蛋白質產物進行三級結構預測，并利用Rasmol軟件對三級結構進行圖形化分析。用Clustal X 對蛋白質多重序列進行比對。利用 Primer 5.0軟件進行引物設計，取玉米新鮮幼嫩的葉片，采用 CTAB 法提取玉米基因組 DNA;利用 Power? pol 2×PCR Mix 進行 PCR 擴增。膠回收試劑盒、大腸桿菌 DH5α、 DNA marker、 Simple -Blunt Peasy-T 載體等購自上海生物工程股份有限公司。所用引物由上海生工負責合成（表1）。

2 結果與分析

2.1 玉米 MS22啟動子克隆

玉米 MS22基因位于玉米第7 條染色體（ Chr7：5333364..5334245）短臂，通過 Primer 5.0進行引物設計，提取 B73材料幼嫩葉片的 DNA，進行 PCR 擴增，得到長度為1523 bp 的序列，見圖1。

通過對目的片段進行凝膠回收及連接轉化，挑取單菌落經菌落 PCR（圖2）鑒定正確后測序，對測序成功的序列進行比對分析。

2.2? ROXY1、ROXY2與 MS22基因功能結構及理化性質分析

以 MS22、ROXY1、ROXY2的氨基酸序列為模板，在 NCBI 網站中對其進行基因結構域分析，發現3個蛋白均屬于硫氧還蛋白超級家族，并具備該家族典型的保守結構域及相關功能位點（圖3）。

利用expasy數據庫 tools-protparam在線程序對 MS22、 ROXY1、 ROXY2的理化性質進行分析，結果發現， MS22和 ROXY2為疏水性蛋白，其親水指數分別為0.496和0.474;而 ROXY1為親水蛋白，其親水系數為0.604，3個蛋白的具體理化性質見表2。

2.3 蛋白序列及保守結構域分析

對 MS22、 ROXY1、 ROXY2進行氨基酸序列比對發現， MS22與 ROXY1、ROXY2在第50個氨基酸之后具有較高的匹配度。這一區段也正是 GSH 轉移酶（ GST， N 端結構域）活性位點結合區段，該區段有谷氧還蛋白家族（ GRXs ）特有的結構域，具有氧化還原活性 CXXC 的保守基序。

玉米的 MS22與其他植物的 GRX 類蛋白一樣具有谷氧還蛋白的特殊結構域。玉米的 MS22基因屬于谷氧還蛋白家族中Glutaredoxin–C9的成員。利用在線結構預測分析軟件 MEME3.0對玉米的 MS22、ROXY1、ROXY2基因的氨基酸序列進行分析發現，3個基因均具有 GRX 類蛋白典型的保守結構域（圖4）。

利用在線網站 http：//meme-suite.org/分析 Glutaredoxin-C9-like 蛋白的保守基序，在氨基酸序列的第49～97的位點處有一個含有49個氨基酸殘基的保守基序 motif。分析發現，該類蛋白的氨基酸序列包含 GRX 基因家族特有的活性位點結構域： CY LV C 或者 CMR K S 的保守基序，CXXC 與 CXXS 為該類群蛋白質的特征標志（圖5）。保守結構域分析發現，玉米 MS22屬于谷氧還蛋白亞類 III，結構域分析發現，該蛋白具有 GRX-subIII的CCx[C/S/GI]典型活性位點。

2.4? MS22及同源基因的蛋白結構預測

采用在線分析軟件 Swiss-Model （ http ：//swissmodel.expasy.org/）預測玉米 MS22基因的蛋白三級結構。 MS22蛋白的三級結構主要以α-螺旋、β-轉角、無規則卷曲及延伸結構為主。如圖6所示， a-螺旋為藍色，不規則卷曲為黃色，延伸鏈結構為紅色，β-轉角為綠色，玉米 MS22蛋白的三維結構組成以α-螺旋、不規則卷曲、延伸結構為主。應用 SOPMA 軟件預測 MS22蛋白的二級結構，其中α-螺旋占38.12%、β-轉角占14.37%、延伸鏈占21.25%、無規卷曲占26.25%。

利用Ensembl? Genomes 網站（http：//ense- mblgenomes.org/）檢索玉米的 zm00001d018802。結果發現，玉米 MS22有37個直系同源物種基因和19個旁系同源物種基因（圖7）。主要包括的農作物有小麥、大麥、水稻、大豆、高粱、谷子等，此外還有苜蓿、卷柏、樟子松、葡萄等。分析發現，這類基因主要存在于陸生禾本科開花植物中，在植物中主要對電子載體活性及蛋白質二硫氧化還原酶活性起調控作用。

3 討論與結論

GRXs 蛋白是一類小分子氧化還原酶，是一類約包含 100個氨基酸殘基及氧化還原活性 CXXC 結構的小分子蛋白;該類蛋白主要在植物中調節氧化還原代謝反應及相關通路信號的識別;該類基因在擬南芥、水稻、高粱、玉米、谷子及一些花卉植物、蔬菜中以不同的作用形式存在[12-15]。研究發現，玉米 MS22、水稻 ROXY1、 ROXY2蛋白均屬于硫氧還蛋白超級家族中谷氧還蛋白家族（ GRXs）成員，該類蛋白目前已經被證實參與植物不同的細胞代謝過程，例如植物的氧化脅迫、干旱、鹽害、高溫、冷害、重金屬等抗性，植物個體對一些真菌類病原菌的抗性水平方面[16-18]。

本研究對水稻 ROXY1、ROXY2及玉米 MS22 的 DNA 序列、氨基酸序列進行比對發現， MS22、 ROXY1、ROXY2具有較高的序列相似性，其在基因結構、保守基序、功能結構域上都具有 GRXs 蛋白家族典型的基因結構。說明 MS22在植物中不僅可導致玉米雄性不育發生，而且在不同植物、不同生育時期，其功能類似于同為禾本科作物水稻 ROXY1和 ROXY2等其他 GRXs 類蛋白家族成員，調節蛋白質的氧化還原狀態、抵御不同生物（非生物）逆境環境脅迫條件，參與植物不同器官的生長發育調控過程。本研究對玉米雄性核不育基因 MS22有進一步的認識，對其他作物中該類型基因的研究具有一定的理論參考價值。

參考文獻

[1]? 金危危， 賴錦盛， 張祖新， 等.玉米骨干自交系深度重測序及有利基因挖掘[J].科技創新導報， 2016， 13（1）：169.

[2]? 田有輝， 萬向元.玉米花藥發育的細胞生物學與分子遺傳學的研究方法[J].中國生物工程雜志， 2018， 38（1）：88-99.

[3]? 蘇愛國， 宋偉， 王帥帥， 等.玉米細胞質雄性不育及其育性恢復基因的研究進展[J].中國生物工程雜志 ， 2018， 38（1）：108-114.

[4]? 張勤 ， 金圣浩 ， 方芳 ， 等. 玉米光溫敏雄性不育系CB1208-82的生理生化代謝研究[J].玉米科學 ， 2019， 27（3）：48-53.

[5]? 田士可 ， 秦心兒 ， 張文亮 ， 等. 玉米雄性不育突變體mi-ms-3的遺傳分析及分子鑒定[J].作物學報.2020， 46（12）：1991-1996.

[6]? 馬萌.玉米細胞核雄性不育基因定位及克隆的研究進展[J].中國農業科學， 2020， 10（5）：332-341.

[7]? 黃淦， 王瀟， 金學鋒， 等.擬南芥谷氧還蛋白 GRXC9負調控葉片大小[J].植物學報， 2017， 52（5）：550-559.

[8]? 湯玲， 趙敏， 張軍， 等.蠟梅谷氧還蛋白基因 CpGRX1的克隆與表達分析[J].西南大學學報（自然科學版）， 2014， 36（2）：41-47.

[9]? 王雪， 阮孟斌， 王斌， 等.木薯干旱脅迫響應CC類谷氧還蛋白基因的篩選與體外表達分析[J].分子植物育種， 2017，15（4）：1204-1211.

[10] Ziemann M， Bhave M， Zachgo S. Origin and Diversificationof land plant CC-Type glutaredox -ins [J]. Genome Biology & Evolution， 2009， 1：265-277.

[11] Xing? S，? Zachgo? S. ROXY1 and ROXY2，? two? Arabidopsisglutaredox -in genes， are required for anther development[J].The Plant Journal， 2010， 53（5）：790-801.

[12] Pérez-Pérez M E， Lemaire S D， Crespo J L， et al. Reactiveoxygen species and autophagy in plants and algae1[J]. PlantPhysiology， 2012，160（1）：156-164.

[13] Yushuang Guo， Changjun Huang， Yan Xie， et al. A tomatoglutar-edoxin gene SlGRX1 regulates plant responses to oxi- dative， drou -ght and salt stresses[J]. Planta， 2010， 232（6）：1499-1509.

[14] Hu Y， Wu Q， Sprague S A， et al. Tomato expressing Arabi-dopsisglutaredoxin? gene? AtGRXS17 confers? tolerance? to? chilling stress via modulating cold responsive components[J]. Horticulture Research， 2015， 2：15051.

[15] Schieber M， Chandel N. ROS function in redox? signalingand? oxidative? stress[J]. Current? Biology， 2014， 24（10）：R453-R462.

[16] Zhang D， Wu S， An X， et al. Construction of a multicontrolsterility system for a maize malesterile line and hybrid seed production based on the ZmMs7 gene encoding a PHD-finger transcription? factor[J]. Plant? Biotechnology? Journal， 2018， 16（2）.

[17] TayebehKakeshpour，Tej Man Tamang，GergelyMotolai，et al.CGFS-type glutaredoxin mutations reduce tolerance to mu -ltiple abiotic stresses in tomato[J]. Ecophysiology， 2021， 3，173：1263-1279.

[18] Ma? B，? Suo? Y，? Zhang? J，? et? al. Glutaredoxin? like? protein（RtGRL1） regulates H2O2 and Na+ accumulation by main- taining? the? glutathione pool? during? abiotic? stress[J]. Plant Physiology and Biochemistry， 2020， 159（2SI）.

（責任編輯 林海妹）