Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Transcriptome Sequencing and <italic>De Novo</italic> Analysis of Cytoplasmic Male Sterility and Maintenance in JA-CMS Cotton</title><author><name sortKey="Yang, Peng" sort="Yang, Peng" uniqKey="Yang P" first="Peng" last="Yang">Peng Yang</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Agronomy, Henan Agricultural University, Zhengzhou, Henan, China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Department of Rural Development, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Han, Jinfeng" sort="Han, Jinfeng" uniqKey="Han J" first="Jinfeng" last="Han">Jinfeng Han</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Agronomy, Henan Agricultural University, Zhengzhou, Henan, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Huang, Jinling" sort="Huang, Jinling" uniqKey="Huang J" first="Jinling" last="Huang">Jinling Huang</name><affiliation><nlm:aff id="aff3"><addr-line>Department of Agronomy, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">25372034</idno><idno type="pmc">4221291</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4221291</idno><idno type="RBID">PMC:4221291</idno><idno type="doi">10.1371/journal.pone.0112320</idno><date when="2014">2014</date><idno type="wicri:Area/Pmc/Corpus">000216</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Transcriptome Sequencing and <italic>De Novo</italic> Analysis of Cytoplasmic Male Sterility and Maintenance in JA-CMS Cotton</title><author><name sortKey="Yang, Peng" sort="Yang, Peng" uniqKey="Yang P" first="Peng" last="Yang">Peng Yang</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Agronomy, Henan Agricultural University, Zhengzhou, Henan, China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Department of Rural Development, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Han, Jinfeng" sort="Han, Jinfeng" uniqKey="Han J" first="Jinfeng" last="Han">Jinfeng Han</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Agronomy, Henan Agricultural University, Zhengzhou, Henan, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Huang, Jinling" sort="Huang, Jinling" uniqKey="Huang J" first="Jinling" last="Huang">Jinling Huang</name><affiliation><nlm:aff id="aff3"><addr-line>Department of Agronomy, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2014">2014</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Cytoplasmic male sterility (CMS) is the failure to produce functional pollen, which is inherited maternally. And it is known that anther development is modulated through complicated interactions between nuclear and mitochondrial genes in sporophytic and gametophytic tissues. However, an unbiased transcriptome sequencing analysis of CMS in cotton is currently lacking in the literature. This study compared differentially expressed (DE) genes of floral buds at the sporogenous cells stage (SS) and microsporocyte stage (MS) (the two most important stages for pollen abortion in JA-CMS) between JA-CMS and its fertile maintainer line JB cotton plants, using the Illumina HiSeq 2000 sequencing platform. A total of 709 (1.8%) DE genes including 293 up-regulated and 416 down-regulated genes were identified in JA-CMS line comparing with its maintainer line at the SS stage, and 644 (1.6%) DE genes with 263 up-regulated and 381 down-regulated genes were detected at the MS stage. By comparing the two stages in the same material, there were 8 up-regulated and 9 down-regulated DE genes in JA-CMS line and 29 up-regulated and 9 down-regulated DE genes in JB maintainer line at the MS stage. Quantitative RT-PCR was used to validate 7 randomly selected DE genes. Bioinformatics analysis revealed that genes involved in reduction-oxidation reactions and alpha-linolenic acid metabolism were down-regulated, while genes pertaining to photosynthesis and flavonoid biosynthesis were up-regulated in JA-CMS floral buds compared with their JB counterparts at the SS and/or MS stages. All these four biological processes play important roles in reactive oxygen species (ROS) homeostasis, which may be an important factor contributing to the sterile trait of JA-CMS. Further experiments are warranted to elucidate molecular mechanisms of these genes that lead to CMS.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Young, Eg" uniqKey="Young E">EG Young</name></author><author><name sortKey="Hanson, Mr" uniqKey="Hanson M">MR Hanson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chase, Cd" uniqKey="Chase C">CD Chase</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Linke, B" uniqKey="Linke B">B Linke</name></author><author><name sortKey="Borner, T" uniqKey="Borner T">T Borner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Unseld, M" uniqKey="Unseld M">M Unseld</name></author><author><name sortKey="Marienfeld, Jr" uniqKey="Marienfeld J">JR Marienfeld</name></author><author><name sortKey="Brandt, P" uniqKey="Brandt P">P Brandt</name></author><author><name sortKey="Brennicke, A" uniqKey="Brennicke A">A Brennicke</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kubo, T" uniqKey="Kubo T">T Kubo</name></author><author><name sortKey="Nishizawa, S" uniqKey="Nishizawa S">S Nishizawa</name></author><author><name sortKey="Sugawara, A" uniqKey="Sugawara A">A Sugawara</name></author><author><name sortKey="Itchoda, N" uniqKey="Itchoda N">N Itchoda</name></author><author><name sortKey="Estiati, A" uniqKey="Estiati A">A Estiati</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Notsu, Y" uniqKey="Notsu Y">Y Notsu</name></author><author><name sortKey="Masood, S" uniqKey="Masood S">S Masood</name></author><author><name sortKey="Nishikawa, T" uniqKey="Nishikawa T">T Nishikawa</name></author><author><name sortKey="Kubo, N" uniqKey="Kubo N">N Kubo</name></author><author><name sortKey="Akiduki, G" uniqKey="Akiduki G">G Akiduki</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sugiyama, Y" uniqKey="Sugiyama Y">Y Sugiyama</name></author><author><name sortKey="Watase, Y" uniqKey="Watase Y">Y Watase</name></author><author><name sortKey="Nagase, M" uniqKey="Nagase M">M Nagase</name></author><author><name sortKey="Makita, N" uniqKey="Makita N">N Makita</name></author><author><name sortKey="Yagura, S" uniqKey="Yagura S">S Yagura</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rhoads, Dm" uniqKey="Rhoads D">DM Rhoads</name></author><author><name sortKey="Subbaiah, Cc" uniqKey="Subbaiah C">CC Subbaiah</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Suzuki, H" uniqKey="Suzuki H">H Suzuki</name></author><author><name sortKey="Rodriguez Uribe, L" uniqKey="Rodriguez Uribe L">L Rodriguez-Uribe</name></author><author><name sortKey="Xu, J" uniqKey="Xu J">J Xu</name></author><author><name sortKey="Zhang, J" uniqKey="Zhang J">J Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schuster, Sc" uniqKey="Schuster S">SC Schuster</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ansorge, Wj" uniqKey="Ansorge W">WJ Ansorge</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Z" uniqKey="Wang Z">Z Wang</name></author><author><name sortKey="Gerstein, M" uniqKey="Gerstein M">M Gerstein</name></author><author><name sortKey="Snyder, M" uniqKey="Snyder M">M Snyder</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, Bb" uniqKey="Zheng B">BB Zheng</name></author><author><name sortKey="Wu, Xm" uniqKey="Wu X">XM Wu</name></author><author><name sortKey="Ge, Xx" uniqKey="Ge X">XX Ge</name></author><author><name sortKey="Deng, Xx" uniqKey="Deng X">XX Deng</name></author><author><name sortKey="Grosser, Jw" uniqKey="Grosser J">JW Grosser</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yan, X" uniqKey="Yan X">X Yan</name></author><author><name sortKey="Dong, C" uniqKey="Dong C">C Dong</name></author><author><name sortKey="Yu, J" uniqKey="Yu J">J Yu</name></author><author><name sortKey="Liu, W" uniqKey="Liu W">W Liu</name></author><author><name sortKey="Jiang, C" uniqKey="Jiang C">C Jiang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, C" uniqKey="Liu C">C Liu</name></author><author><name sortKey="Ma, N" uniqKey="Ma N">N Ma</name></author><author><name sortKey="Wang, Py" uniqKey="Wang P">PY Wang</name></author><author><name sortKey="Fu, N" uniqKey="Fu N">N Fu</name></author><author><name sortKey="Shen, Hl" uniqKey="Shen H">HL Shen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huang, J" uniqKey="Huang J">J Huang</name></author><author><name sortKey="Yang, P" uniqKey="Yang P">P Yang</name></author><author><name sortKey="Li, B" uniqKey="Li B">B Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jia, Zc" uniqKey="Jia Z">ZC Jia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, H" uniqKey="Zhao H">H Zhao</name></author><author><name sortKey="Huang, J" uniqKey="Huang J">J Huang</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Grabherr, Mg" uniqKey="Grabherr M">MG Grabherr</name></author><author><name sortKey="Haas, Bj" uniqKey="Haas B">BJ Haas</name></author><author><name sortKey="Yassour, M" uniqKey="Yassour M">M Yassour</name></author><author><name sortKey="Levin, Jz" uniqKey="Levin J">JZ Levin</name></author><author><name sortKey="Thompson, Da" uniqKey="Thompson D">DA Thompson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, B" uniqKey="Li B">B Li</name></author><author><name sortKey="Dewey, Cn" uniqKey="Dewey C">CN Dewey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robinson, Md" uniqKey="Robinson M">MD Robinson</name></author><author><name sortKey="Oshlack, A" uniqKey="Oshlack A">A Oshlack</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, L" uniqKey="Wang L">L Wang</name></author><author><name sortKey="Feng, Z" uniqKey="Feng Z">Z Feng</name></author><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Young, Md" uniqKey="Young M">MD Young</name></author><author><name sortKey="Wakefield, Mj" uniqKey="Wakefield M">MJ Wakefield</name></author><author><name sortKey="Smyth, Gk" uniqKey="Smyth G">GK Smyth</name></author><author><name sortKey="Oshlack, A" uniqKey="Oshlack A">A Oshlack</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xie, C" uniqKey="Xie C">C Xie</name></author><author><name sortKey="Mao, X" uniqKey="Mao X">X Mao</name></author><author><name sortKey="Huang, J" uniqKey="Huang J">J Huang</name></author><author><name sortKey="Ding, Y" uniqKey="Ding Y">Y Ding</name></author><author><name sortKey="Wu, J" uniqKey="Wu J">J Wu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Livak, Kj" uniqKey="Livak K">KJ Livak</name></author><author><name sortKey="Schmittgen, Td" uniqKey="Schmittgen T">TD Schmittgen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Subramanian, A" uniqKey="Subramanian A">A Subramanian</name></author><author><name sortKey="Tamayo, P" uniqKey="Tamayo P">P Tamayo</name></author><author><name sortKey="Mootha, Vk" uniqKey="Mootha V">VK Mootha</name></author><author><name sortKey="Mukherjee, S" uniqKey="Mukherjee S">S Mukherjee</name></author><author><name sortKey="Ebert, Bl" uniqKey="Ebert B">BL Ebert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Balk, J" uniqKey="Balk J">J Balk</name></author><author><name sortKey="Leaver, Cj" uniqKey="Leaver C">CJ Leaver</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Papini, A" uniqKey="Papini A">A Papini</name></author><author><name sortKey="Mosti, S" uniqKey="Mosti S">S Mosti</name></author><author><name sortKey="Brighigna, L" uniqKey="Brighigna L">L Brighigna</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Doyle, Sm" uniqKey="Doyle S">SM Doyle</name></author><author><name sortKey="Diamond, M" uniqKey="Diamond M">M Diamond</name></author><author><name sortKey="Mccabe, Pf" uniqKey="Mccabe P">PF McCabe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jiang, Pd" uniqKey="Jiang P">PD Jiang</name></author><author><name sortKey="Zhang, Xq" uniqKey="Zhang X">XQ Zhang</name></author><author><name sortKey="Zhu, Yg" uniqKey="Zhu Y">YG Zhu</name></author><author><name sortKey="Zhu, W" uniqKey="Zhu W">W Zhu</name></author><author><name sortKey="Xie, Hy" uniqKey="Xie H">HY Xie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wan, Cx" uniqKey="Wan C">CX Wan</name></author><author><name sortKey="Li, Sq" uniqKey="Li S">SQ Li</name></author><author><name sortKey="Wen, L" uniqKey="Wen L">L Wen</name></author><author><name sortKey="Kong, J" uniqKey="Kong J">J Kong</name></author><author><name sortKey="Wang, K" uniqKey="Wang K">K Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gechev, Ts" uniqKey="Gechev T">TS Gechev</name></author><author><name sortKey="Van Breusegem, F" uniqKey="Van Breusegem F">F Van Breusegem</name></author><author><name sortKey="Stone, Jm" uniqKey="Stone J">JM Stone</name></author><author><name sortKey="Denev, I" uniqKey="Denev I">I Denev</name></author><author><name sortKey="Laloi, C" uniqKey="Laloi C">C Laloi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Moller, Im" uniqKey="Moller I">IM Moller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sharma, P" uniqKey="Sharma P">P Sharma</name></author><author><name sortKey="Jha, Ab" uniqKey="Jha A">AB Jha</name></author><author><name sortKey="Dubey, Rs" uniqKey="Dubey R">RS Dubey</name></author><author><name sortKey="Pessarakli, M" uniqKey="Pessarakli M">M Pessarakli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dos Santos, Cv" uniqKey="Dos Santos C">CV Dos Santos</name></author><author><name sortKey="Rey, P" uniqKey="Rey P">P Rey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Leon, Ip" uniqKey="De Leon I">IP De Leon</name></author><author><name sortKey="Sanz, A" uniqKey="Sanz A">A Sanz</name></author><author><name sortKey="Hamberg, M" uniqKey="Hamberg M">M Hamberg</name></author><author><name sortKey="Castresana, C" uniqKey="Castresana C">C Castresana</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mcconn, M" uniqKey="Mcconn M">M McConn</name></author><author><name sortKey="Browse, J" uniqKey="Browse J">J Browse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Park, Jh" uniqKey="Park J">JH Park</name></author><author><name sortKey="Halitschke, R" uniqKey="Halitschke R">R Halitschke</name></author><author><name sortKey="Kim, Hb" uniqKey="Kim H">HB Kim</name></author><author><name sortKey="Baldwin, It" uniqKey="Baldwin I">IT Baldwin</name></author><author><name sortKey="Feldmann, Ka" uniqKey="Feldmann K">KA Feldmann</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Yang, S" uniqKey="Yang S">S Yang</name></author><author><name sortKey="Terachi, T" uniqKey="Terachi T">T Terachi</name></author><author><name sortKey="Yamagishi, H" uniqKey="Yamagishi H">H Yamagishi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wei, M" uniqKey="Wei M">M Wei</name></author><author><name sortKey="Song, M" uniqKey="Song M">M Song</name></author><author><name sortKey="Fan, S" uniqKey="Fan S">S Fan</name></author><author><name sortKey="Yu, S" uniqKey="Yu S">S Yu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fini, A" uniqKey="Fini A">A Fini</name></author><author><name sortKey="Brunetti, C" uniqKey="Brunetti C">C Brunetti</name></author><author><name sortKey="Di Ferdinando, M" uniqKey="Di Ferdinando M">M Di Ferdinando</name></author><author><name sortKey="Ferrini, F" uniqKey="Ferrini F">F Ferrini</name></author><author><name sortKey="Tattini, M" uniqKey="Tattini M">M Tattini</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dong, X" uniqKey="Dong X">X Dong</name></author><author><name sortKey="Feng, H" uniqKey="Feng H">H Feng</name></author><author><name sortKey="Xu, M" uniqKey="Xu M">M Xu</name></author><author><name sortKey="Lee, J" uniqKey="Lee J">J Lee</name></author><author><name sortKey="Kim, Yk" uniqKey="Kim Y">YK Kim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Heazlewood, Jl" uniqKey="Heazlewood J">JL Heazlewood</name></author><author><name sortKey="Whelan, J" uniqKey="Whelan J">J Whelan</name></author><author><name sortKey="Millar, Ah" uniqKey="Millar A">AH Millar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kempken, F" uniqKey="Kempken F">F Kempken</name></author><author><name sortKey="Howard, W" uniqKey="Howard W">W Howard</name></author><author><name sortKey="Pring, Dr" uniqKey="Pring D">DR Pring</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De La Canal, L" uniqKey="De La Canal L">L de la Canal</name></author><author><name sortKey="Crouzillat, D" uniqKey="Crouzillat D">D Crouzillat</name></author><author><name sortKey="Quetier, F" uniqKey="Quetier F">F Quetier</name></author><author><name sortKey="Ledoigt, G" uniqKey="Ledoigt G">G Ledoigt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kurek, I" uniqKey="Kurek I">I Kurek</name></author><author><name sortKey="Ezra, D" uniqKey="Ezra D">D Ezra</name></author><author><name sortKey="Begu, D" uniqKey="Begu D">D Begu</name></author><author><name sortKey="Erel, N" uniqKey="Erel N">N Erel</name></author><author><name sortKey="Litvak, S" uniqKey="Litvak S">S Litvak</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Szklarczyk, M" uniqKey="Szklarczyk M">M Szklarczyk</name></author><author><name sortKey="Oczkowski, M" uniqKey="Oczkowski M">M Oczkowski</name></author><author><name sortKey="Augustyniak, H" uniqKey="Augustyniak H">H Augustyniak</name></author><author><name sortKey="Borner, T" uniqKey="Borner T">T Borner</name></author><author><name sortKey="Linke, B" uniqKey="Linke B">B Linke</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hanson, Mr" uniqKey="Hanson M">MR Hanson</name></author><author><name sortKey="Bentolila, S" uniqKey="Bentolila S">S Bentolila</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kim, Dh" uniqKey="Kim D">DH Kim</name></author><author><name sortKey="Kang, Jg" uniqKey="Kang J">JG Kang</name></author><author><name sortKey="Kim, Bd" uniqKey="Kim B">BD Kim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yamamoto, Mp" uniqKey="Yamamoto M">MP Yamamoto</name></author><author><name sortKey="Kubo, T" uniqKey="Kubo T">T Kubo</name></author><author><name sortKey="Mikami, T" uniqKey="Mikami T">T Mikami</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, H" uniqKey="Zhang H">H Zhang</name></author><author><name sortKey="Li, S" uniqKey="Li S">S Li</name></author><author><name sortKey="Yi, P" uniqKey="Yi P">P Yi</name></author><author><name sortKey="Wan, C" uniqKey="Wan C">C Wan</name></author><author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, F" uniqKey="Wang F">F Wang</name></author><author><name sortKey="Feng, Cd" uniqKey="Feng C">CD Feng</name></author><author><name sortKey="O Connell, Ma" uniqKey="O Connell M">MA O'Connell</name></author><author><name sortKey="Stewart, Jm" uniqKey="Stewart J">JM Stewart</name></author><author><name sortKey="Zhang, Jf" uniqKey="Zhang J">JF Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Jy" uniqKey="Wu J">JY Wu</name></author><author><name sortKey="Gong, Yc" uniqKey="Gong Y">YC Gong</name></author><author><name sortKey="Cui, Mh" uniqKey="Cui M">MH Cui</name></author><author><name sortKey="Qi, Tx" uniqKey="Qi T">TX Qi</name></author><author><name sortKey="Guo, Lp" uniqKey="Guo L">LP Guo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Suzuki, H" uniqKey="Suzuki H">H Suzuki</name></author><author><name sortKey="Yu, Jw" uniqKey="Yu J">JW Yu</name></author><author><name sortKey="Wang, F" uniqKey="Wang F">F Wang</name></author><author><name sortKey="Zhang, Jf" uniqKey="Zhang J">JF Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, J" uniqKey="Li J">J Li</name></author><author><name sortKey="Pandeya, D" uniqKey="Pandeya D">D Pandeya</name></author><author><name sortKey="Jo, Yd" uniqKey="Jo Y">YD Jo</name></author><author><name sortKey="Liu, Wy" uniqKey="Liu W">WY Liu</name></author><author><name sortKey="Kang, Bc" uniqKey="Kang B">BC Kang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yuan, J" uniqKey="Yuan J">J Yuan</name></author><author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name></author><author><name sortKey="Liu, X" uniqKey="Liu X">X Liu</name></author><author><name sortKey="Hao, R" uniqKey="Hao R">R Hao</name></author><author><name sortKey="Wang, H" uniqKey="Wang H">H Wang</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Ma, W" uniqKey="Ma W">W Ma</name></author><author><name sortKey="Ren, A" uniqKey="Ren A">A Ren</name></author><author><name sortKey="Cui, M" uniqKey="Cui M">M Cui</name></author><author><name sortKey="Han, Q" uniqKey="Han Q">Q Han</name></author><author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">25372034</article-id><article-id pub-id-type="pmc">4221291</article-id><article-id pub-id-type="publisher-id">PONE-D-14-15411</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0112320</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology and Life Sciences</subject><subj-group><subject>Agriculture</subject><subj-group><subject>Agronomy</subject><subj-group><subject>Plant Products</subject><subj-group><subject>Cotton</subject></subj-group></subj-group></subj-group><subj-group><subject>Crop Science</subject><subj-group><subject>Crops</subject><subj-group><subject>Fiber Crops</subject></subj-group></subj-group></subj-group></subj-group><subj-group><subject>Genetics</subject><subj-group><subject>Plant Genetics</subject><subj-group><subject>Crop Genetics</subject></subj-group></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Transcriptome Sequencing and <italic>De Novo</italic> Analysis of Cytoplasmic Male Sterility and Maintenance in JA-CMS Cotton</article-title><alt-title alt-title-type="running-head">Transcriptomic Analysis of Cytoplasmic Male Sterile Cotton</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Yang</surname><given-names>Peng</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Han</surname><given-names>Jinfeng</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Huang</surname><given-names>Jinling</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Agronomy, Henan Agricultural University, Zhengzhou, Henan, China</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Department of Rural Development, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></aff><aff id="aff3"><label>3</label><addr-line>Department of Agronomy, Shanxi Agricultural University, Taigu, Shanxi, China</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Zhang</surname><given-names>Jinfa</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>New Mexico State University, United States of America</addr-line></aff><author-notes><corresp id="cor1">* E-mail: <email>huangjl@sxau.edu.cn</email></corresp><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con"><p>Conceived and designed the experiments: PY J. Han J. Huang. Performed the experiments: PY J. Huang. Analyzed the data: PY J. Huang. Contributed reagents/materials/analysis tools: J. Han J. Huang. Contributed to the writing of the manuscript: PY J. Han J. Huang.</p></fn></author-notes><pub-date pub-type="collection"><year>2014</year></pub-date><pub-date pub-type="epub"><day>5</day><month>11</month><year>2014</year></pub-date><volume>9</volume><issue>11</issue><elocation-id>e112320</elocation-id><history><date date-type="received"><day>6</day><month>4</month><year>2014</year></date><date date-type="accepted"><day>8</day><month>10</month><year>2014</year></date></history><permissions><copyright-year>2014</copyright-year><copyright-holder>Yang et al</copyright-holder><license><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><abstract><p>Cytoplasmic male sterility (CMS) is the failure to produce functional pollen, which is inherited maternally. And it is known that anther development is modulated through complicated interactions between nuclear and mitochondrial genes in sporophytic and gametophytic tissues. However, an unbiased transcriptome sequencing analysis of CMS in cotton is currently lacking in the literature. This study compared differentially expressed (DE) genes of floral buds at the sporogenous cells stage (SS) and microsporocyte stage (MS) (the two most important stages for pollen abortion in JA-CMS) between JA-CMS and its fertile maintainer line JB cotton plants, using the Illumina HiSeq 2000 sequencing platform. A total of 709 (1.8%) DE genes including 293 up-regulated and 416 down-regulated genes were identified in JA-CMS line comparing with its maintainer line at the SS stage, and 644 (1.6%) DE genes with 263 up-regulated and 381 down-regulated genes were detected at the MS stage. By comparing the two stages in the same material, there were 8 up-regulated and 9 down-regulated DE genes in JA-CMS line and 29 up-regulated and 9 down-regulated DE genes in JB maintainer line at the MS stage. Quantitative RT-PCR was used to validate 7 randomly selected DE genes. Bioinformatics analysis revealed that genes involved in reduction-oxidation reactions and alpha-linolenic acid metabolism were down-regulated, while genes pertaining to photosynthesis and flavonoid biosynthesis were up-regulated in JA-CMS floral buds compared with their JB counterparts at the SS and/or MS stages. All these four biological processes play important roles in reactive oxygen species (ROS) homeostasis, which may be an important factor contributing to the sterile trait of JA-CMS. Further experiments are warranted to elucidate molecular mechanisms of these genes that lead to CMS.</p></abstract><funding-group><funding-statement>This project was supported by Shanxi Scientific and Technological Project Grant 20130311004–3 to JH and Shanxi Agricultural University Ph.D. fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="13"></page-count></counts><custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name><meta-value>The authors confirm that all data underlying the findings are fully available without restriction. The transcriptome sequencing data from this study have been deposited in the NCBI SRA database and are accessible through accession numbers SRX547770, SRX547777, SRX547779, and SRX547781 (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link>).</meta-value></custom-meta></custom-meta-group></article-meta><notes><title>Data Availability</title><p>The authors confirm that all data underlying the findings are fully available without restriction. The transcriptome sequencing data from this study have been deposited in the NCBI SRA database and are accessible through accession numbers SRX547770, SRX547777, SRX547779, and SRX547781 (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link>).</p></notes></front><body><sec id="s1"><title>Introduction</title><p>Cytoplasmic male sterility (CMS) is a maternally inherited trait in higher plants incapable of producing functional pollen <xref rid="pone.0112320-Young1" ref-type="bibr">[1]</xref>. However, CMS plants exhibit normal vegetative growth and female fertility, and this trait can be restored by nuclear genes known as <italic>restorer-of-fertility</italic> (<italic>Rf</italic>) genes <xref rid="pone.0112320-Chase1" ref-type="bibr">[2]</xref>. On one hand, it has been widely accepted that CMS is closely related to mitochondrial genome rearrangement, which creates chimeric genes that disturb normal pollen development <xref rid="pone.0112320-Linke1" ref-type="bibr">[3]</xref>. More than 50 mitochondrial genes have been identified as CMS-related in various plants <xref rid="pone.0112320-Unseld1" ref-type="bibr">[4]</xref>–<xref rid="pone.0112320-Sugiyama1" ref-type="bibr">[7]</xref>, which are valuable in producing F1 hybrid cultivars with heterosis. On the other hand, changes in mitochondrial function trigger altered nuclear gene expression by a mysterious process called mitochondrial retrograde regulation (MRR) <xref rid="pone.0112320-Rhoads1" ref-type="bibr">[8]</xref>. Much is still unknown about signalling pathways involved in this process, MRR targets, and how cells switch to programmed cell death (PCD) mode in the case of CMS.</p><p>Transcriptomic analysis using DNA microarray or RNA-seq technology has been implemented to investigate molecular markers or differentially expressed (DE) genes implicated in a variety of traits. Suzuki et al. <xref rid="pone.0112320-Suzuki1" ref-type="bibr">[9]</xref> presented the first genome-wide transcriptome analysis of CMS and restoration in cotton using Affymetrix GeneChips Cotton Genome Array. They compared DE genes of floral buds at the meiosis stage between CMS-D8 and its restorer line, and found 458 (1.9%) DE genes including 127 up-regulated and 331 down-regulated ones. The most frequent DE gene group was involved in cell wall expansion <xref rid="pone.0112320-Suzuki1" ref-type="bibr">[9]</xref>.</p><p>In recent years, the emergence and advancement of next generation sequencing technology has made a breakthrough in life sciences, by offering unprecedented speed and cost efficiency to study genomic and transcriptomic data <xref rid="pone.0112320-Schuster1" ref-type="bibr">[10]</xref>, <xref rid="pone.0112320-Ansorge1" ref-type="bibr">[11]</xref>. Comparing with DNA microarray, RNA-seq has very low, if any, background signal, and has also been shown to be highly accurate for quantifying expression levels <xref rid="pone.0112320-Wang1" ref-type="bibr">[12]</xref>. Researchers have used RNA-seq to study CMS in sweet orange (<italic>Citrus sinensis</italic>) <xref rid="pone.0112320-Zheng1" ref-type="bibr">[13]</xref>, rape (<italic>Brassica napus</italic>) <xref rid="pone.0112320-Yan1" ref-type="bibr">[14]</xref>, and chili pepper (<italic>Capsicum annuum</italic> L.) <xref rid="pone.0112320-Liu1" ref-type="bibr">[15]</xref>, to give a few examples.</p><p>Zheng et al. <xref rid="pone.0112320-Zheng1" ref-type="bibr">[13]</xref> compared nuclear gene expression profiles of male sterile cybrid and fertile pummelo floral buds by RNA-seq analysis. <italic>APETALA3</italic> and <italic>PISTILLATA</italic> transcripts, which encode key transcription factors for stamen identification and are restricted to normal floral whorls, were repressed in the sterile line. As the authors stated, these citrus class-B MADS-box genes are likely to be targets for CMS retrograde signalling <xref rid="pone.0112320-Zheng1" ref-type="bibr">[13]</xref>. In <italic>Brassica napus</italic>, 3,231 genes of <italic>Brassica rapa</italic> and 3,371 genes of <italic>Brassica oleracea</italic> were detected with significantly different expression levels from young floral buds of sterile and fertile plants, which were derived from self-pollinated offspring of the F<sub>1</sub> hybrid from novel restorer line NR1 and <italic>Nsa</italic> CMS line <xref rid="pone.0112320-Yan1" ref-type="bibr">[14]</xref>. Altered expressions of genes involved in carbon metabolism, tricarboxylic acid cycle (TCA cycle), oxidative phosphorylation, oxidoreductase activity and pentatricopeptide repeat (PPR) proteins were observed by Yan et al. <xref rid="pone.0112320-Yan1" ref-type="bibr">[14]</xref>. Liu et al. <xref rid="pone.0112320-Liu1" ref-type="bibr">[15]</xref>’s study profiled anther transcriptomes of a chili pepper CMS line 121A and its restorer line 121C, and found the top three pathways covering the most DE unigenes were starch and sucrose metabolism, oxidative phosphorylation and plant-pathogen interaction <xref rid="pone.0112320-Liu1" ref-type="bibr">[15]</xref>.</p><p>However, unbiased transcriptome sequencing analysis of CMS in cotton is lacking. To narrow this gap, we isolated, sequenced and quantified transcriptomes of floral buds at the sporogenous cells stage (SS) and microsporocyte stage (MS), which were demonstrated to be the two most important stages for pollen sterility <xref rid="pone.0112320-Huang1" ref-type="bibr">[16]</xref> in JA-CMS, a CMS line cultivated by Cotton Breeding Lab of Shanxi Agricultural University (Please see <xref ref-type="sec" rid="s2">materials and methods</xref> for details about JA-CMS). The obtained transcriptomic sequences were then assembled, annotated, and analysed to discover DE genes between the two lines at the two specific stages, to find pathways that were enriched with these DE genes, and to further determine potential key genes functional in CMS.</p></sec><sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title><sec id="s2a"><title>Plant materials and RNA extraction</title><p>JA-CMS and its maintainer line JB cotton plants were grown in the fields of Shanxi Agricultural University, Taigu, Shanxi, China. The two lines were isolated by repeated backcrossing of (<italic>Gossypium hirsutum</italic>×<italic>G. thurberi</italic>)×(<italic>G. arboreum</italic>×<italic>G. hirsutum</italic>) for more than 10 generations using MB177, one kind of <italic>G. hirsutum</italic>, as the recurrent parent. With 100% sterile plants, JA-CMS is unique in the following three ways: it has the genetic background of both <italic>G. thurberi</italic> and <italic>G. arboreum</italic>; it matures around 20 days earlier than 104-7A <xref rid="pone.0112320-Jia1" ref-type="bibr">[17]</xref>, another CMS line cultivated by Chinese scientists through backcrossing of <italic>G. hirsutum</italic> and <italic>G. barbadense</italic>; it is conditioned by <italic>G. hirsutum</italic> cytoplasm, which is different from CMS lines of <italic>G. harknessii</italic>, 104-7A and <italic>G. trilobum</italic>, CMS-D8, based on origins and phenotypes. Please refer to the supplementary note for a detailed comparison of JA-CMS, 104-7A, CMS-D2 and CMS-D8. Our group has studied morphological and cellular characteristics of JA-CMS, determined major development stages for pollen abortion, and conducted biochemical and molecular biological research on CMS. In collaboration with the Institute of Cotton Study, Shanxi Agricultural Academy, a restorer line of JA-CMS named JB was also cultivated <xref rid="pone.0112320-Huang1" ref-type="bibr">[16]</xref>, <xref rid="pone.0112320-Zhao1" ref-type="bibr">[18]</xref>.</p><p>Because of the stable correlation between development stages of pollen mother cells and floral bud morphology <xref rid="pone.0112320-Hou1" ref-type="bibr">[19]</xref>, following the method described by Zhao and Huang <xref rid="pone.0112320-Zhao1" ref-type="bibr">[18]</xref>, three floral buds, with diameters less than 2 mm or between 2 and 3 mm, were collected from five sterile and fertile plants at the same time, respectively, to confirm the correlation in our samples using cytological methods. As the results suggested, for JA-CMS and its maintainer line, a floral bud is at the SS stage when its diameter is between 1.5–2.2 mm, and at the MS stage when its diameter reaches 2.2–2.6 mm.</p><p>Floral buds were harvested from at least 50 plants of each line at apporoximately the same time of the day; 2 g of floral buds were mixed as a pool based on different stages and cotton lines. Pooled materials were snap-frozen in liquid nitrogen and stored at −80°C for RNA preparation. Total RNA was isolated using an RN37-EASYspin RNA extraction kit (Aidlab Biotechnologies, China) according to the manufacturer’s protocol. The integrity of total RNA was checked by 1% agarose gel electrophoresis; the concentration and purity were determined using NanoDrop (Thermo Scientific, USA) and Agilent 2100 Bioanalyzer (Agilent, USA). A total amount of 3 µg RNA per sample was used as input material for the RNA sample preparations. All four samples had RIN (RNA Integrity Number) values above 8.</p></sec><sec id="s2b"><title>Library preparation and sequencing</title><p>Sequencing libraries were generated using Illumina TruSeq RNA Sample Preparation Kit (Illumina, USA) following manufacturer’s recommendations and four index codes were added to attribute sequences to each sample. Briefly, after integrity, concentration and purity check, poly-(T) oligo-attached magnetic beads were utilized for mRNA enrichment. Then fragmentation buffer was added to break mRNA into short fragments, which served as templates for the first strand cDNA synthesis with random hexamer-primer. Buffer, dNTPs, RNase H and DNA polymerase I were added afterwards to synthesize the second strand. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities and enzymes were removed. The double-stranded cDNA was purified with AMPure XP beads, followed by end repair, poly-(A) addition and sequencing adaptor ligation. Fragments preferentially 200 bp in length were enriched using Illumina PCR Primer Cocktail in a 10 cycle PCR amplification to obtain the cDNA library. To ensure quality of the library, Qubit 2.0 (Life Technologies, USA), Agilent 2100 (Agilent, USA) and quantitative PCR were used for initial quantification, insert size check and effective concentration determination (library effective concentration>2 nM), respectively.</p><p>The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina, USA) according to the manufacturer’s instructions. After cluster generation, the library was sequenced on Illumina HiSeq 2000 platform; 100 bp paired-end reads were generated. Low-quality reads containing ambiguous nucleotides or adaptor sequences were removed from the raw reads to gain clean reads, which were then assembled using Trinity to construct unique consensus sequences without a reference genome <xref rid="pone.0112320-Grabherr1" ref-type="bibr">[20]</xref>.</p></sec><sec id="s2c"><title>Bioinformatics analysis</title><p>Seven databases were used for unigene annotation, including the NCBI protein database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/protein">http://www.ncbi.nlm.nih.gov/protein</ext-link>), the NCBI nucleotide database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/nucleotide">http://www.ncbi.nlm.nih.gov/nucleotide</ext-link>), the Pfam database (<ext-link ext-link-type="uri" xlink:href="http://pfam.xfam.org/">http://pfam.xfam.org/</ext-link>), the euKaryotic Orthologous Groups (KOG; <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/COG/">http://www.ncbi.nlm.nih.gov/COG/</ext-link>), the UniProt Knowledgebase (UniProtKB)/Swiss-Prot (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/uniprot">http://www.ebi.ac.uk/uniprot</ext-link>), the KEGG Orthology database (KO; <ext-link ext-link-type="uri" xlink:href="http://www.genome.jp/kegg/ko.html">http://www.genome.jp/kegg/ko.html</ext-link>), and the Gene Ontology (GO; <ext-link ext-link-type="uri" xlink:href="http://www.geneontology.org/">http://www.geneontology.org/</ext-link>). We performed these annotations by using a combination of BLAST (<ext-link ext-link-type="uri" xlink:href="http://blast.ncbi.nlm.nih.gov/Blast.cgi">http://blast.ncbi.nlm.nih.gov/Blast.cgi</ext-link>), HMMER (<ext-link ext-link-type="uri" xlink:href="http://hmmer.janelia.org/">http://hmmer.janelia.org/</ext-link>), Blast2GO (<ext-link ext-link-type="uri" xlink:href="http://www.blast2go.com/b2ghome">http://www.blast2go.com/b2ghome</ext-link>), and KEGG Automatic Annotation Server (KAAS; <ext-link ext-link-type="uri" xlink:href="http://www.genome.jp/kegg/kaas/">http://www.genome.jp/kegg/kaas/</ext-link>).</p><p>Using the transcriptome assembled by Trinity <xref rid="pone.0112320-Grabherr1" ref-type="bibr">[20]</xref> as reference, clean reads of each sample were mapped to it by RSEM (RNA-Seq by Expectation-Maximization) <xref rid="pone.0112320-Li1" ref-type="bibr">[21]</xref>. For biologically non-duplicating samples, mapped read counts were normalized through applying trimmed mean of M values (TMM) method <xref rid="pone.0112320-Robinson1" ref-type="bibr">[22]</xref>. Then normalized counts were analysed and unigenes showing significant expression differences were determined using DEGseq <xref rid="pone.0112320-Wang2" ref-type="bibr">[23]</xref>, with FDR q-value threshold set to 0.005 and absolute value of expression fold change set to 2.</p><p>After selection of unigenes with significant transcript abundance difference (DE genes), we conducted enrichment analyses, i.e. GO and KEGG analysis, to elucidate distributions of these unigenes based on their functions and attached biological pathways. GO analysis was performed using Goseq <xref rid="pone.0112320-Young2" ref-type="bibr">[24]</xref>, which is based on Wallenius non-central hypergeometric distribution. KEGG analysis utilizes hypergeometric test to search DE genes that are significantly enriched in certain KEGG pathways compared with other genes in the genome. KOBAS 2.0 <xref rid="pone.0112320-Xie1" ref-type="bibr">[25]</xref> was used for KEGG analysis. GO categories and KEGG pathways with FDR q-value ≤ 0.05 were selected as significantly enriched.</p><p>KEGG analysis results are presented graphically using scatter plots, with three parameters for enrichment level assessment, which are rich factor, the ratio between counts of DE genes and all annotated genes enriched in a certain pathway (higher the ratio, more significant the pathway); qvalue, P value after multiple hypothesis testing correction with a range between 0 and 1 (closer to 0, more significant the pathway); and number of genes enriched. Twenty most significant pathways were plotted, when more than twenty pathways were identified for each of the four comparisons.</p></sec><sec id="s2d"><title>Quantitative RT-PCR</title><p>Quantitative RT-PCR (qRT-PCR) was performed to validate transcript abundance of unigenes observed by RNA-seq. Primer-BLAST (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/tools/primer-blast/">http://www.ncbi.nlm.nih.gov/tools/primer-blast/</ext-link>) was used to design gene-specific primers for seven randomly selected DE genes plus elongation factor-1α gene as the internal control (<xref ref-type="supplementary-material" rid="pone.0112320.s014">Table S1</xref>). Reactions were carried out using the SYBR <italic>Premix Ex Taq</italic> II (Tli RnaseH Plus) Kit (Takara, USA) in the iQ5 Multicolour Real-Time PCR Detection System (Bio-Rad, USA), according to the manufacturer’s instructions.</p><p>The first-strand cDNA was synthesized using the PrimeScript RT Master Mix Kit (Takara, USA) under the following conditions: 37°C for 15 min (reverse transcription) and 85°C for 5 s (denaturing reverse transcriptase). Then amplification reactions of 20 µL volume consisting of 10 µL SYBR <italic>Premix Ex Taq</italic> II, 0.8 µL of 10 µM forward and reverse primers each, 1 µL cDNA template and sterile water, were conducted, with cycling parameters as follows: 95°C for 30 s, 40 cycles of 95°C for 5 s and 60°C for 30 s. Relative expression levels of unigenes from different samples were calculated using 2<sup>−ΔΔCT</sup> method <xref rid="pone.0112320-Livak1" ref-type="bibr">[26]</xref>. The qRT-PCR was conducted with 3 replicates for each sample, and data were indicated as means ± standard errors (n = 3) in <xref ref-type="supplementary-material" rid="pone.0112320.s001">Figure S1</xref>.</p></sec><sec id="s2e"><title>Data access</title><p>The transcriptome sequencing data from this study have been deposited in the NCBI SRA database and are accessible through accession numbers SRX547770, SRX547777, SRX547779, and SRX547781 (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link>).</p></sec></sec><sec id="s3"><title>Results</title><sec id="s3a"><title>RNA sequencing and gene annotations</title><p>A total of 2.8×10<sup>8</sup> raw reads were obtained from the Illumina sequencing platform. Each stage of either sterile or fertile line had around 6.8 ± 0.6 giga base pairs (Gb) of sequencing data. After data quality control as described in materials and methods, 97.54% ± 0.26% of the raw reads were clean reads for each of the four samples (<xref ref-type="table" rid="pone-0112320-t001">Table 1</xref>). Trinity <xref rid="pone.0112320-Grabherr1" ref-type="bibr">[20]</xref> reconstructed all the clean reads into a <italic>de novo</italic> reference transcriptome; the longest transcript for each gene was recognized as the unigene. Following these steps, 206,496 transcripts and 86,093 unigenes were identified, length distributions of which are presented in <xref ref-type="supplementary-material" rid="pone.0112320.s002">Figure S2</xref> and <xref ref-type="supplementary-material" rid="pone.0112320.s003">S3</xref>. The transcripts and unigenes are between 201 bp and 17,233 bp in length; the mean lengths of the transcripts and unigenes are 1,184 bp and 755 bp, respectively. Around 36.3% of the transcripts and 59.8% of the unigenes are within the range of 200–500 bp.</p><table-wrap id="pone-0112320-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.t001</object-id><label>Table 1</label><caption><title>Transcriptome sequencing data quality.</title></caption><alternatives><graphic id="pone-0112320-t001-1" xlink:href="pone.0112320.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sample</td><td align="left" rowspan="1" colspan="1">Raw reads</td><td align="left" rowspan="1" colspan="1">Clean reads</td><td align="left" rowspan="1" colspan="1">Clean bases</td><td align="left" rowspan="1" colspan="1">Error (%)</td><td align="left" rowspan="1" colspan="1">Q20 (%)</td><td align="left" rowspan="1" colspan="1">Q30 (%)</td><td align="left" rowspan="1" colspan="1">GC (%)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">JA-CMS at the SS stage</td><td align="left" rowspan="1" colspan="1">66695222</td><td align="left" rowspan="1" colspan="1">64995444</td><td align="left" rowspan="1" colspan="1">6.50 Gb</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">98.05</td><td align="left" rowspan="1" colspan="1">93.60</td><td align="left" rowspan="1" colspan="1">43.79</td></tr><tr><td align="left" rowspan="1" colspan="1">JA-CMS at the MS stage</td><td align="left" rowspan="1" colspan="1">77998708</td><td align="left" rowspan="1" colspan="1">75844362</td><td align="left" rowspan="1" colspan="1">7.58 Gb</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">98.02</td><td align="left" rowspan="1" colspan="1">93.57</td><td align="left" rowspan="1" colspan="1">43.65</td></tr><tr><td align="left" rowspan="1" colspan="1">JB at the SS stage</td><td align="left" rowspan="1" colspan="1">71264156</td><td align="left" rowspan="1" colspan="1">69566718</td><td align="left" rowspan="1" colspan="1">6.96 Gb</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">98.01</td><td align="left" rowspan="1" colspan="1">93.45</td><td align="left" rowspan="1" colspan="1">43.74</td></tr><tr><td align="left" rowspan="1" colspan="1">JB at the MS stage</td><td align="left" rowspan="1" colspan="1">62779022</td><td align="left" rowspan="1" colspan="1">61431768</td><td align="left" rowspan="1" colspan="1">6.14 Gb</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">98.11</td><td align="left" rowspan="1" colspan="1">93.71</td><td align="left" rowspan="1" colspan="1">43.75</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>Notes: Gb = Giga base pair; Q20 represents an error rate of 1 in 100, with a corresponding call accuracy of 99%; Q30 represents an error rate of 1 in 1000, with a corresponding call accuracy of 99.9%.</p></fn></table-wrap-foot></table-wrap><p>The percentage of successfully annotated unigenes differs among the seven databases specified in the <xref ref-type="sec" rid="s2">Materials and Methods</xref> part, with the NCBI protein database having the highest, 42.05% and KO providing the lowest, 11.76%. Only 5.15% of the unigenes are annotated in all the databases, while 49.47% are annotated in at least one database. Annotations are provided whenever available.</p></sec><sec id="s3b"><title>Genome-wide analysis of DE genes</title><p>Gene expression levels at the SS or MS stage for JA-CMS and JB floral buds were quantified and compared using RSEM <xref rid="pone.0112320-Li1" ref-type="bibr">[21]</xref>, TMM <xref rid="pone.0112320-Robinson1" ref-type="bibr">[22]</xref>, and DEGseq <xref rid="pone.0112320-Wang2" ref-type="bibr">[23]</xref>. At the SS stage, there were 11,111 genes uniquely expressed in JA-CMS plants and 7,811 in JB plants; 40,318 genes were shared between the two lines at this stage. At the MS stage, 8,502 and 8,636 genes were specifically expressed in JA-CMS and JB lines, and the number of commonly expressed genes was 39,967. Regarding gene expression level comparison between the two stages of the same plant material, JA-CMS plants had 9,334 genes uniquely expressed at the SS stage and 6,374 at the MS stage; 42,905 genes were expressed at both stages; and there were 7,366 and 7,840 genes expressed at the SS and MS stages of JB plants, respectively, while 40,763 genes were commonly found at both stages (<xref ref-type="fig" rid="pone-0112320-g001">Figure 1</xref>).</p><fig id="pone-0112320-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.g001</object-id><label>Figure 1</label><caption><title>Venn diagram of DE gene counts at the SS and MS stages for JA-CMS and JB floral buds.</title><p>B2 = SS stage of JA-CMS, K2 = SS stage of JB, K3 = MS stage of JB, B3 = MS stage of JA-CMS (from right to left).</p></caption><graphic xlink:href="pone.0112320.g001"></graphic></fig><p>Due to the large number of genes, we used FDR q value<0.005 combined with |log<sub>2</sub>(fold change)|>1 to select DE genes between the two materials at the same stage, or between different stages of the same material. By comparing with the JB line, JA-CMS plants had 709 genes with significantly different expression levels at the SS stage, among which 293 were up-regulated and 416 were down-regulated (<xref ref-type="supplementary-material" rid="pone.0112320.s015">Table S2</xref>). At the MS stage, 644 genes passed the significance threshold, with 263 up-regulated and 381 down-regulated in JA-CMS plants (<xref ref-type="supplementary-material" rid="pone.0112320.s016">Table S3</xref>). Fewer genes showed significant expression difference between the two stages of the same material. Seventeen genes, with 8 up-regulated and 9 down-regulated, had significantly different expression levels at the MS stage of JA-CMS comparing with the SS stage (<xref ref-type="supplementary-material" rid="pone.0112320.s017">Table S4</xref>). There were 38 DE genes showing significance in JB line between the two stages, among which 29 were up-regulated and 9 were down-regulated at the MS stage (<xref ref-type="supplementary-material" rid="pone.0112320.s018">Table S5</xref>).</p><p>All the 854 DE genes that showed significance in at least one of the four comparisons have 80.76 ± 158.10 and 78.07 ± 150.42 RPKM (Reads per kilobase per million) values at the SS and MS stages of JA-CMS, respectively. The RPKM values for JB plants have larger means and standard deviations than JA-CMS at both stages: 91.40 ± 301.55 (SS stage) and 96.25 ± 304.16 (MS stage). These values indicate reasonable coverage for the DE genes. Please refer to <xref ref-type="supplementary-material" rid="pone.0112320.s004">Figure S4</xref> for the density distributions of log10(RPKM) for the four samples.</p></sec><sec id="s3c"><title>Enrichment analysis of DE genes</title><p>Enrichment analysis was performed for the purpose of identifying major functions or pathways associated with CMS when genes function co-ordinately under certain biological mechanisms instead of individually to cause specific phenotype, which is the case for CMS. This type of analysis prioritizes DE genes for further investigation, if many DE genes are identified, and significantly increases statistical power to find predisposing genes, if few are available <xref rid="pone.0112320-Subramanian1" ref-type="bibr">[27]</xref>.</p><sec id="s3c1"><title>GO analysis</title><p>GO covers three domains: biological process (BP), cellular component (CC) and molecular function (MF). DE genes between the two materials at the SS stage were significantly enriched in 27 GO categories. Processes, components and functions related to reduction-oxidation (redox) reactions, photosynthesis, and transcription factors showed up in two out of the three GO domains. Oxidation-reduction process (GO:0055114) encompassed 100 DE genes with the smallest q value of all the 27 categories, which equals to 7.9×10<sup>−4</sup> (<xref ref-type="table" rid="pone-0112320-t002">Table 2</xref>). We also conducted GO analysis on up-regulated and down-regulated DE genes at the SS stage, separately (<xref ref-type="supplementary-material" rid="pone.0112320.s005">Figure S5</xref> and <xref ref-type="supplementary-material" rid="pone.0112320.s006">S6</xref>). For the up-regulated genes, their product properties were mostly enriched in the biological process of photosynthesis and cellular components of photosynthetic membranes with both q values equal to 6.31×10<sup>−6</sup>; besides these two categories, other processes, components and functions related to photosynthesis were also prominent, such as photosystem (q = 7.48×10<sup>−6</sup>), thylakoid part and thylakoid (q = 1.50×10<sup>−5</sup>; <xref ref-type="supplementary-material" rid="pone.0112320.s005">Figure S5</xref>). On the other hand, down-regulated gene product properties were over-represented in only four categories, which are all pertinent to redox reactions (<xref ref-type="supplementary-material" rid="pone.0112320.s006">Figure S6</xref>).</p><table-wrap id="pone-0112320-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.t002</object-id><label>Table 2</label><caption><title>GO analysis results of DE genes at the SS and MS stages between JA-CMS and JB.</title></caption><alternatives><graphic id="pone-0112320-t002-2" xlink:href="pone.0112320.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">GOaccessionID</td><td align="left" rowspan="1" colspan="1">Description</td><td align="left" rowspan="1" colspan="1">Domain</td><td colspan="3" align="left" rowspan="1">SS stage</td><td colspan="5" align="left" rowspan="1">MS stage</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Pvalue</td><td align="left" rowspan="1" colspan="1">Qvalue</td><td align="left" rowspan="1" colspan="1">DEgene #</td><td colspan="2" align="left" rowspan="1">Pvalue</td><td colspan="2" align="left" rowspan="1">Qvalue</td><td align="left" rowspan="1" colspan="1">DEgene #</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GO:0055114</td><td align="left" rowspan="1" colspan="1">Oxidation-reduction process</td><td align="left" rowspan="1" colspan="1">BP</td><td align="left" rowspan="1" colspan="1">1.53E-07</td><td align="left" rowspan="1" colspan="1">7.90E-04</td><td align="left" rowspan="1" colspan="1">100</td><td colspan="2" align="left" rowspan="1">6.79E-08</td><td colspan="2" align="left" rowspan="1">2.43E-04</td><td align="left" rowspan="1" colspan="1">92</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016705</td><td align="left" rowspan="1" colspan="1">Oxidoreductase activity, acting on paired donors, withincorporation or reduction of molecular oxygen</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">6.46E-07</td><td align="left" rowspan="1" colspan="1">1.66E-03</td><td align="left" rowspan="1" colspan="1">35</td><td colspan="2" align="left" rowspan="1">1.02E-07</td><td colspan="2" align="left" rowspan="1">2.43E-04</td><td align="left" rowspan="1" colspan="1">34</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0001071</td><td align="left" rowspan="1" colspan="1">Nucleic acid binding transcription factor activity</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">1.58E-06</td><td align="left" rowspan="1" colspan="1">2.04E-03</td><td align="left" rowspan="1" colspan="1">48</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0003700</td><td align="left" rowspan="1" colspan="1">Sequence-specific DNA binding transcription factor activity</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">1.58E-06</td><td align="left" rowspan="1" colspan="1">2.04E-03</td><td align="left" rowspan="1" colspan="1">48</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016491</td><td align="left" rowspan="1" colspan="1">Oxidoreductase activity</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">2.73E-06</td><td align="left" rowspan="1" colspan="1">2.44E-03</td><td align="left" rowspan="1" colspan="1">96</td><td colspan="2" align="left" rowspan="1">4.37E-07</td><td colspan="2" align="left" rowspan="1">4.51E-04</td><td align="left" rowspan="1" colspan="1">89</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005667</td><td align="left" rowspan="1" colspan="1">Transcription factor complex</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">2.83E-06</td><td align="left" rowspan="1" colspan="1">2.44E-03</td><td align="left" rowspan="1" colspan="1">55</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0009521</td><td align="left" rowspan="1" colspan="1">Photosystem</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">4.40E-06</td><td align="left" rowspan="1" colspan="1">3.24E-03</td><td align="left" rowspan="1" colspan="1">17</td><td colspan="2" align="left" rowspan="1">6.62E-05</td><td colspan="2" align="left" rowspan="1">3.10E-02</td><td align="left" rowspan="1" colspan="1">14</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0034357</td><td align="left" rowspan="1" colspan="1">Photosynthetic membrane</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">5.60E-06</td><td align="left" rowspan="1" colspan="1">3.61E-03</td><td align="left" rowspan="1" colspan="1">18</td><td colspan="2" align="left" rowspan="1">6.58E-05</td><td colspan="2" align="left" rowspan="1">3.10E-02</td><td align="left" rowspan="1" colspan="1">15</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0046872</td><td align="left" rowspan="1" colspan="1">Metal ion binding</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">8.76E-06</td><td align="left" rowspan="1" colspan="1">4.86E-03</td><td align="left" rowspan="1" colspan="1">116</td><td colspan="2" align="left" rowspan="1">1.89E-07</td><td colspan="2" align="left" rowspan="1">2.43E-04</td><td align="left" rowspan="1" colspan="1">110</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0043169</td><td align="left" rowspan="1" colspan="1">Cation binding</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">9.43E-06</td><td align="left" rowspan="1" colspan="1">4.86E-03</td><td align="left" rowspan="1" colspan="1">118</td><td colspan="2" align="left" rowspan="1">1.83E-07</td><td colspan="2" align="left" rowspan="1">2.43E-04</td><td align="left" rowspan="1" colspan="1">112</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005654</td><td align="left" rowspan="1" colspan="1">Nucleoplasm</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.20E-05</td><td align="left" rowspan="1" colspan="1">5.13E-03</td><td align="left" rowspan="1" colspan="1">58</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044451</td><td align="left" rowspan="1" colspan="1">Nucleoplasm part</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.20E-05</td><td align="left" rowspan="1" colspan="1">5.13E-03</td><td align="left" rowspan="1" colspan="1">58</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0015979</td><td align="left" rowspan="1" colspan="1">Photosynthesis</td><td align="left" rowspan="1" colspan="1">BP</td><td align="left" rowspan="1" colspan="1">1.40E-05</td><td align="left" rowspan="1" colspan="1">5.54E-03</td><td align="left" rowspan="1" colspan="1">26</td><td colspan="2" align="left" rowspan="1">1.07E-04</td><td align="left" rowspan="1" colspan="1">4.59E-02</td><td colspan="2" align="left" rowspan="1">22</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0008152</td><td align="left" rowspan="1" colspan="1">Metabolic process</td><td align="left" rowspan="1" colspan="1">BP</td><td align="left" rowspan="1" colspan="1">1.57E-05</td><td align="left" rowspan="1" colspan="1">5.77E-03</td><td align="left" rowspan="1" colspan="1">394</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044436</td><td align="left" rowspan="1" colspan="1">Thylakoid part</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">2.34E-05</td><td align="left" rowspan="1" colspan="1">8.05E-03</td><td align="left" rowspan="1" colspan="1">18</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0009579</td><td align="left" rowspan="1" colspan="1">Thylakoid</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">2.81E-05</td><td align="left" rowspan="1" colspan="1">8.58E-03</td><td align="left" rowspan="1" colspan="1">18</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016706</td><td align="left" rowspan="1" colspan="1">Oxidoreductase activity, acting on paired donors, with incorporationor reduction of molecular oxygen, 2-oxoglutarate as one donor,and incorporation of one atom each of oxygen into both donors</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">2.83E-05</td><td align="left" rowspan="1" colspan="1">8.58E-03</td><td align="left" rowspan="1" colspan="1">14</td><td colspan="2" align="left" rowspan="1">7.40E-06</td><td align="left" rowspan="1" colspan="1">5.45E-03</td><td colspan="2" align="left" rowspan="1">14</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0031981</td><td align="left" rowspan="1" colspan="1">Nuclear lumen</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">7.38E-05</td><td align="left" rowspan="1" colspan="1">2.11E-02</td><td align="left" rowspan="1" colspan="1">63</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0032787</td><td align="left" rowspan="1" colspan="1">Monocarboxylic acid metabolic process</td><td align="left" rowspan="1" colspan="1">BP</td><td align="left" rowspan="1" colspan="1">8.31E-05</td><td align="left" rowspan="1" colspan="1">2.25E-02</td><td align="left" rowspan="1" colspan="1">24</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0009538</td><td align="left" rowspan="1" colspan="1">Photosystem I reaction center</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">9.21E-05</td><td align="left" rowspan="1" colspan="1">2.37E-02</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">5.10E-05</td><td colspan="2" align="left" rowspan="1">2.92E-02</td><td colspan="2" align="left" rowspan="1">5</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044428</td><td align="left" rowspan="1" colspan="1">Nuclear part</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.07E-04</td><td align="left" rowspan="1" colspan="1">2.62E-02</td><td align="left" rowspan="1" colspan="1">75</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0043233</td><td align="left" rowspan="1" colspan="1">Organelle lumen</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.21E-04</td><td align="left" rowspan="1" colspan="1">2.71E-02</td><td align="left" rowspan="1" colspan="1">63</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0070013</td><td align="left" rowspan="1" colspan="1">Intracellular organelle lumen</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.21E-04</td><td align="left" rowspan="1" colspan="1">2.71E-02</td><td align="left" rowspan="1" colspan="1">63</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044710</td><td align="left" rowspan="1" colspan="1">Single-organism metabolic process</td><td align="left" rowspan="1" colspan="1">BP</td><td align="left" rowspan="1" colspan="1">1.47E-04</td><td align="left" rowspan="1" colspan="1">3.15E-02</td><td align="left" rowspan="1" colspan="1">184</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0031974</td><td align="left" rowspan="1" colspan="1">Membrane-enclosed lumen</td><td align="left" rowspan="1" colspan="1">CC</td><td align="left" rowspan="1" colspan="1">1.55E-04</td><td align="left" rowspan="1" colspan="1">3.19E-02</td><td align="left" rowspan="1" colspan="1">63</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0051213</td><td align="left" rowspan="1" colspan="1">Dioxygenase activity</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">1.90E-04</td><td align="left" rowspan="1" colspan="1">3.77E-02</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">3.09E-06</td><td colspan="2" align="left" rowspan="1">2.65E-03</td><td colspan="2" align="left" rowspan="1">17</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0048037</td><td align="left" rowspan="1" colspan="1">Cofactor binding</td><td align="left" rowspan="1" colspan="1">MF</td><td align="left" rowspan="1" colspan="1">2.17E-04</td><td align="left" rowspan="1" colspan="1">4.13E-02</td><td align="left" rowspan="1" colspan="1">38</td><td colspan="5" align="left" rowspan="1">Notsignificant</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0046914</td><td align="left" rowspan="1" colspan="1">Transition metal ion binding</td><td align="left" rowspan="1" colspan="1">MF</td><td colspan="3" align="left" rowspan="1">Notsignificant</td><td align="left" rowspan="1" colspan="1">2.03E-05</td><td colspan="2" align="left" rowspan="1">1.31E-02</td><td colspan="2" align="left" rowspan="1">83</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Notes: 1) BP = Biological process, CC = Cellular component, MF = Molecular function; Q value = FDR corrected p value. 2) Only GO categories with FDR q value ≤ 0.05 are presented here. 3) GO categories are arranged based on q values from smallest to largest at the SS stage.</p></fn></table-wrap-foot></table-wrap><p>For the MS stage, product properties of DE genes between JA-CMS and JB lines were separated into 12 categories (<xref ref-type="table" rid="pone-0112320-t002">Table 2</xref>). The most significant categories (q = 2.43×10<sup>−4</sup>) were oxidation-reduction process, oxidoreductase activity (acting on paired donors, with incorporation or reduction of molecular oxygen), cation binding and metal ion binding. Same as the SS stage, redox reactions and photosynthesis covered most of the GO categories. Additional GO analysis on up-regulated DE genes showed that, besides photosynthesis related categories, 7 out of the 43 significantly enriched classes were pertinent to transcription (<xref ref-type="supplementary-material" rid="pone.0112320.s007">Figure S7</xref>). Consistent with the SS stage analysis results, every category identified in GO analysis on down-regulated DE genes at the MS stage was involved in redox reactions (<xref ref-type="supplementary-material" rid="pone.0112320.s008">Figure S8</xref>). By comparing JA-CMS and JB lines at the two stages, it became obvious that genes implicated in redox reactions were persistently down-regulated, while genes related to photosynthesis were continuously up-regulated.</p><p>GO analysis did not find any significantly enriched category concerning up-regulated DE genes at the MS stage comparing with the SS stage of sterile plants; however, 24 categories were over-represented by down-regulated DE gene products, all of which were related to ATP biosynthesis (<xref ref-type="supplementary-material" rid="pone.0112320.s009">Figure S9</xref>). Combined analysis of up-regulated and down-regulated DE genes revealed the same category distribution as the down-regulated analysis. No GO enrichment was detected for DE genes between the two development stages of fertile plants.</p></sec><sec id="s3c2"><title>KEGG analysis</title><p>There were 134 pathways determined for DE genes between JA-CMS and JB plants at the SS stage, among which the most significant ones were “photosynthesis” (gene number [N] = 15), “flavonoid biosynthesis” (N = 12), “metabolic pathways” (N = 112), “DNA replication” (N = 9), and “alpha-linolenic acid metabolism” (N = 9; <xref ref-type="fig" rid="pone-0112320-g002">Figure 2A</xref>). Consistent with GO analysis, “photosynthesis” was significantly enriched with up-regulated DE genes, as well as “flavonoid biosynthesis” and “DNA replication” (<xref ref-type="supplementary-material" rid="pone.0112320.s010">Figure S10</xref>). Down-regulated DE genes were over-represented in “alpha-linolenic acid metabolism” (<xref ref-type="supplementary-material" rid="pone.0112320.s011">Figure S11</xref>). Four out of the five pathways identified by KEGG for the MS stage were the same as the SS stage, which are “photosynthesis” (N = 14), “flavonoid biosynthesis” (N = 10), “DNA replication” (N = 9), and “metabolic pathways” (N = 99). Another pathway showing significance was “drug metabolism” (N = 7) among the total 140 pathways identified at the MS stage (<xref ref-type="fig" rid="pone-0112320-g002">Figure 2B</xref>). The three pathways significantly enriched with up-regulated genes at this stage were the same as the SS stage (<xref ref-type="supplementary-material" rid="pone.0112320.s012">Figure S12</xref>); “Metabolism of xenobiotics by cytochrome P450” and “drug metabolism” were enriched with down-regulated genes (<xref ref-type="supplementary-material" rid="pone.0112320.s013">Figure S13</xref>).</p><fig id="pone-0112320-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.g002</object-id><label>Figure 2</label><caption><title>KEGG analysis results of DE genes at the SS and MS stages comparing JA-CMS and JB lines.</title><p>Rich factor is the ratio between counts of DE genes and all annotated genes enriched in a certain pathway; qvalue is P value after multiple hypothesis testing correction with a range between 0 and 1. Twenty most significant pathways were plotted, when more than 20 pathways were identified. A) Results for the SS stage. B) Results for the MS stage.</p></caption><graphic xlink:href="pone.0112320.g002"></graphic></fig><p>Because fewer DE genes were determined between the two stages of the same material, no significantly enriched pathway was detected.</p></sec></sec><sec id="s3d"><title>qRT-PCR verification on gene expression patterns</title><p>Although quantitative differences existed between qRT-PCR analysis results and sequencing data for the selected seven DE genes (<xref ref-type="supplementary-material" rid="pone.0112320.s014">Table S1</xref> and <xref ref-type="supplementary-material" rid="pone.0112320.s001">Figure S1</xref>), overall expression tendencies were the same (<xref ref-type="fig" rid="pone-0112320-g003">Figure 3</xref>). For example, comp54402_c0 is 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) gene included in amino acid metabolism pathway, expression of which was down-regulated in sterile plants. Sequencing data showed 16-fold difference of transcript abundance between JA-CMS and JB lines at the SS stage, while qRT-PCR analysis yielded around 14-fold difference. Chlorophyll a/b-binding protein gene, comp70789_c0, was up-regulated in sterile plants compared with fertile counterparts at the MS stage based on sequencing data (15-fold), which was confirmed by the 11-fold increase shown in the corresponding qRT-PCR analysis result.</p><fig id="pone-0112320-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.g003</object-id><label>Figure 3</label><caption><title>Quantitative RT-PCR validations for expression profiles of seven selected DE genes.</title><p>Comp54402_c0 is 1-aminocyclopropane-1-carboxylate synthase gene, comp77990_c0 is alpha-expansin 11 precursor gene, comp71873_c0 is a predicted protein gene, comp68006_c0 is a conserved hypothetical protein gene, comp70789_c0 is chloroplast chlorophyll A–B binding protein gene, comp88938_c1 is transcription factor LHY gene, comp86790_c1 is <italic>G. hirsutum</italic> cab gene for chlorophyll A–B binding protein. Ratio stands for log2 ratio of relative expression levels for qRT-PCR and RPKM (Reads per kilobase per million) values for RNA sequencing, respectively. B2 = SS stage of JA-CMS, K2 = SS stage of JB, B3 = MS stage of JA-CMS, K3 = MS stage of JB.</p></caption><graphic xlink:href="pone.0112320.g003"></graphic></fig></sec></sec><sec id="s4"><title>Discussion</title><p>As a result of our unbiased transcriptome sequencing analysis of JA-CMS and its maintainer line, four candidate functional categories or pathways related to CMS have been identified, which include two down-regulated gene groups in JA-CMS: genes related to redox reactions and alpha-linolenic acid metabolism, and two up-regulated gene groups: genes pertaining to photosynthesis and flavonoid biosynthesis. <xref ref-type="sec" rid="s4">Discussion</xref> of these genes would not be possible without mentioning the preceding study results of JA-CMS.</p><p>According to previous cytological studies, premature microspore abortion and tapetum degradation at the SS stage is one prominent characteristic of JA-CMS, and enzymatic activities of peroxidase, superoxide dismutase, cytochrome oxidase and succinate dehydrogenase were found significantly lower in JA-CMS comparing with its fertile counterpart <xref rid="pone.0112320-Huang1" ref-type="bibr">[16]</xref>. Tapetum PCD plays an important role in pollen development <xref rid="pone.0112320-Balk1" ref-type="bibr">[28]</xref>–<xref rid="pone.0112320-Papini1" ref-type="bibr">[30]</xref>, and abnormal reactive oxygen species (ROS) level is believed to be a signal promoting PCD in plants <xref rid="pone.0112320-Doyle1" ref-type="bibr">[31]</xref>. Disrupted ROS homeostasis was found in several CMS systems, e.g., excessive accumulation of O<sub>2</sub><sup>−</sup>, H<sub>2</sub>O<sub>2</sub> and malondialdehyde (MDA) existed in CMS cotton at the abortion peak <xref rid="pone.0112320-Jiang1" ref-type="bibr">[32]</xref>, and mitochondria of a rice CMS line suffered from serious oxidative stress, which was induced by abnormal increased ROS at the meiosis stage <xref rid="pone.0112320-Wan1" ref-type="bibr">[33]</xref>.</p><p>In plants, chloroplasts are a major site of ROS generation, where photosynthetic electron transfer chains produce O<sub>2</sub><sup>−</sup><xref rid="pone.0112320-Gechev1" ref-type="bibr">[34]</xref>. Because ROS can cause damage to proteins, lipids, and DNA, its production and removal must be strictly controlled <xref rid="pone.0112320-Moller1" ref-type="bibr">[35]</xref>. Plants possess complex antioxidative defense system comprising of nonenzymatic and enzymatic components to scavenge ROS. Nonenzymatic components include the major cellular redox buffers ascorbate (AsA) and glutathione (γ-glutamyl-cysteinyl-glycine, GSH), which act as cofactors for enzymes of AsA-GSH cycle ascorbate peroxidase (APX), monodehydroascorbate reductase (MDHAR), dehydroascorbate reductase (DHAR), glutathione reductase (GR), glutathione peroxidase (GPX), and glutathione S transferase (GST) <xref rid="pone.0112320-Sharma1" ref-type="bibr">[36]</xref>.</p><p>Based on both GO and KEGG analysis results, up-regulated genes in JA-CMS were significantly enriched in photosynthetic categories and pathways, while oxidation-reduction process and oxidoreductase activity categories were significantly over-represented with down-regulated DE genes at the SS and MS stages in sterile plants. There are 14 DE genes up-regulated in JA-CMS belonged to the KEGG pathway of “photosynthesis”, among which Photosystem I subunit PsaO encoding gene (comp65579_c0) showed the maximum fold change of expression at both SS (14-fold) and MS (12-fold) stages (<xref ref-type="table" rid="pone-0112320-t003">Table 3</xref>).</p><table-wrap id="pone-0112320-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.t003</object-id><label>Table 3</label><caption><title>Up-regulated photosynthesis and flavonoid biosynthesis related genes in JA-CMS at the SS and/or MS stages.</title></caption><alternatives><graphic id="pone-0112320-t003-3" xlink:href="pone.0112320.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Function</td><td align="left" rowspan="1" colspan="1">Gene ID</td><td colspan="3" align="left" rowspan="1">SS stage</td><td colspan="3" align="left" rowspan="1">MS stage</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Fold</td><td align="left" rowspan="1" colspan="1">P value</td><td align="left" rowspan="1" colspan="1">Q value</td><td align="left" rowspan="1" colspan="1">Fold</td><td align="left" rowspan="1" colspan="1">P value</td><td align="left" rowspan="1" colspan="1">Q value</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Photosynthesis</td><td align="left" rowspan="1" colspan="1">Photosystem I subunit X</td><td align="left" rowspan="1" colspan="1">comp53760_c0</td><td align="left" rowspan="1" colspan="1">2.53</td><td align="left" rowspan="1" colspan="1">2.23E-08</td><td align="left" rowspan="1" colspan="1">2.32E-06</td><td align="left" rowspan="1" colspan="1">2.56</td><td align="left" rowspan="1" colspan="1">1.23E-07</td><td align="left" rowspan="1" colspan="1">1.23E-05</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Plastocyanin</td><td align="left" rowspan="1" colspan="1">comp65386_c0</td><td align="left" rowspan="1" colspan="1">3.02</td><td align="left" rowspan="1" colspan="1">1.48E-30</td><td align="left" rowspan="1" colspan="1">7.37E-28</td><td align="left" rowspan="1" colspan="1">2.61</td><td align="left" rowspan="1" colspan="1">6.19E-24</td><td align="left" rowspan="1" colspan="1">2.34E-21</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem I subunit PsaO</td><td align="left" rowspan="1" colspan="1">comp65579_c0</td><td align="left" rowspan="1" colspan="1">3.77</td><td align="left" rowspan="1" colspan="1">8.58E-17</td><td align="left" rowspan="1" colspan="1">1.99E-14</td><td align="left" rowspan="1" colspan="1">3.61</td><td align="left" rowspan="1" colspan="1">4.29E-14</td><td align="left" rowspan="1" colspan="1">8.67E-12</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II PsbW protein</td><td align="left" rowspan="1" colspan="1">comp65639_c0</td><td align="left" rowspan="1" colspan="1">2.34</td><td align="left" rowspan="1" colspan="1">5.34E-07</td><td align="left" rowspan="1" colspan="1">4.77E-05</td><td align="left" rowspan="1" colspan="1">2.35</td><td align="left" rowspan="1" colspan="1">1.28E-06</td><td align="left" rowspan="1" colspan="1">1.14E-04</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II 10kDa protein</td><td align="left" rowspan="1" colspan="1">comp68020_c0</td><td align="left" rowspan="1" colspan="1">1.17</td><td align="left" rowspan="1" colspan="1">3.48E-10</td><td align="left" rowspan="1" colspan="1">4.54E-08</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II PsbY protein</td><td align="left" rowspan="1" colspan="1">comp68793_c0</td><td align="left" rowspan="1" colspan="1">1.72</td><td align="left" rowspan="1" colspan="1">4.38E-07</td><td align="left" rowspan="1" colspan="1">3.96E-05</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem I subunit XI</td><td align="left" rowspan="1" colspan="1">comp70643_c0</td><td align="left" rowspan="1" colspan="1">1.83</td><td align="left" rowspan="1" colspan="1">3.25E-15</td><td align="left" rowspan="1" colspan="1">6.79E-13</td><td align="left" rowspan="1" colspan="1">1.60</td><td align="left" rowspan="1" colspan="1">9.63E-11</td><td align="left" rowspan="1" colspan="1">1.44E-08</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Ferredoxin</td><td align="left" rowspan="1" colspan="1">comp71277_c0</td><td align="left" rowspan="1" colspan="1">3.11</td><td align="left" rowspan="1" colspan="1">4.21E-18</td><td align="left" rowspan="1" colspan="1">1.07E-15</td><td align="left" rowspan="1" colspan="1">2.68</td><td align="left" rowspan="1" colspan="1">9.12E-15</td><td align="left" rowspan="1" colspan="1">1.94E-12</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem I subunit IV</td><td align="left" rowspan="1" colspan="1">comp72189_c0</td><td align="left" rowspan="1" colspan="1">1.36</td><td align="left" rowspan="1" colspan="1">1.59E-08</td><td align="left" rowspan="1" colspan="1">1.70E-06</td><td align="left" rowspan="1" colspan="1">1.07</td><td align="left" rowspan="1" colspan="1">1.35E-05</td><td align="left" rowspan="1" colspan="1">1.01E-03</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II oxygen-evolving enhancer protein 1</td><td align="left" rowspan="1" colspan="1">comp73151_c0</td><td align="left" rowspan="1" colspan="1">2.96</td><td align="left" rowspan="1" colspan="1">7.52E-57</td><td align="left" rowspan="1" colspan="1">9.24E-54</td><td align="left" rowspan="1" colspan="1">2.74</td><td align="left" rowspan="1" colspan="1">1.66E-45</td><td align="left" rowspan="1" colspan="1">1.70E-42</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II 22kDa protein</td><td align="left" rowspan="1" colspan="1">comp73229_c0</td><td align="left" rowspan="1" colspan="1">2.28</td><td align="left" rowspan="1" colspan="1">2.04E-14</td><td align="left" rowspan="1" colspan="1">4.05E-12</td><td align="left" rowspan="1" colspan="1">2.14</td><td align="left" rowspan="1" colspan="1">4.45E-12</td><td align="left" rowspan="1" colspan="1">7.67E-10</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II oxygen-evolving enhancer protein 2</td><td align="left" rowspan="1" colspan="1">comp73376_c0</td><td align="left" rowspan="1" colspan="1">1.18</td><td align="left" rowspan="1" colspan="1">2.53E-10</td><td align="left" rowspan="1" colspan="1">3.38E-08</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem I subunit III</td><td align="left" rowspan="1" colspan="1">comp73879_c0</td><td align="left" rowspan="1" colspan="1">1.83</td><td align="left" rowspan="1" colspan="1">2.89E-20</td><td align="left" rowspan="1" colspan="1">8.55E-18</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">1.73E-14</td><td align="left" rowspan="1" colspan="1">3.61E-12</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Photosystem II oxygen-evolving enhancer protein 3</td><td align="left" rowspan="1" colspan="1">comp84774_c2</td><td align="left" rowspan="1" colspan="1">1.83</td><td align="left" rowspan="1" colspan="1">4.26E-10</td><td align="left" rowspan="1" colspan="1">5.51E-08</td><td align="left" rowspan="1" colspan="1">1.55</td><td align="left" rowspan="1" colspan="1">2.13E-07</td><td align="left" rowspan="1" colspan="1">2.10E-05</td></tr><tr><td align="left" rowspan="1" colspan="1">Flavonoid biosynthesis</td><td align="left" rowspan="1" colspan="1">Chalcone synthase</td><td align="left" rowspan="1" colspan="1">comp72843_c0</td><td align="left" rowspan="1" colspan="1">2.23</td><td align="left" rowspan="1" colspan="1">2.61E-41</td><td align="left" rowspan="1" colspan="1">2.07E-38</td><td align="left" rowspan="1" colspan="1">1.63</td><td align="left" rowspan="1" colspan="1">1.96E-29</td><td align="left" rowspan="1" colspan="1">1.03E-26</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp72843_c1</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">6.78E-13</td><td align="left" rowspan="1" colspan="1">1.18E-10</td><td align="left" rowspan="1" colspan="1">1.48</td><td align="left" rowspan="1" colspan="1">4.46E-09</td><td align="left" rowspan="1" colspan="1">5.49E-07</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp84628_c1</td><td align="left" rowspan="1" colspan="1">2.06</td><td align="left" rowspan="1" colspan="1">5.90E-62</td><td align="left" rowspan="1" colspan="1">8.22E-59</td><td align="left" rowspan="1" colspan="1">1.56</td><td align="left" rowspan="1" colspan="1">2.09E-43</td><td align="left" rowspan="1" colspan="1">2.00E-40</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp84628_c0</td><td align="left" rowspan="1" colspan="1">1.88</td><td align="left" rowspan="1" colspan="1">2.55E-78</td><td align="left" rowspan="1" colspan="1">5.51E-75</td><td align="left" rowspan="1" colspan="1">1.56</td><td align="left" rowspan="1" colspan="1">2.23E-68</td><td align="left" rowspan="1" colspan="1">4.13E-65</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Anthocyanidin reductase</td><td align="left" rowspan="1" colspan="1">comp73018_c0</td><td align="left" rowspan="1" colspan="1">1.27</td><td align="left" rowspan="1" colspan="1">1.48E-39</td><td align="left" rowspan="1" colspan="1">1.09E-36</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Leucoanthocyanidin dioxygenase</td><td align="left" rowspan="1" colspan="1">comp75804_c0</td><td align="left" rowspan="1" colspan="1">2.17</td><td align="left" rowspan="1" colspan="1">1.28E-94</td><td align="left" rowspan="1" colspan="1">4.02E-91</td><td align="left" rowspan="1" colspan="1">1.89</td><td align="left" rowspan="1" colspan="1">3.33E-81</td><td align="left" rowspan="1" colspan="1">8.49E-78</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase</td><td align="left" rowspan="1" colspan="1">comp75884_c0</td><td align="left" rowspan="1" colspan="1">1.75</td><td align="left" rowspan="1" colspan="1">2.37E-69</td><td align="left" rowspan="1" colspan="1">3.81E-66</td><td align="left" rowspan="1" colspan="1">1.54</td><td align="left" rowspan="1" colspan="1">8.80E-60</td><td align="left" rowspan="1" colspan="1">1.35E-56</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Trans-cinnamate 4-monooxygenase</td><td align="left" rowspan="1" colspan="1">comp80855_c0</td><td align="left" rowspan="1" colspan="1">1.47</td><td align="left" rowspan="1" colspan="1">1.42E-71</td><td align="left" rowspan="1" colspan="1">2.54E-68</td><td align="left" rowspan="1" colspan="1">1.09</td><td align="left" rowspan="1" colspan="1">7.25E-44</td><td align="left" rowspan="1" colspan="1">7.05E-41</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Flavonoid 3′,5′-hydroxylase</td><td align="left" rowspan="1" colspan="1">comp81439_c0</td><td align="left" rowspan="1" colspan="1">2.50</td><td align="left" rowspan="1" colspan="1">1.74E-93</td><td align="left" rowspan="1" colspan="1">4.97E-90</td><td align="left" rowspan="1" colspan="1">1.86</td><td align="left" rowspan="1" colspan="1">3.08E-57</td><td align="left" rowspan="1" colspan="1">4.39E-54</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Naringenin 3-dioxygenase</td><td align="left" rowspan="1" colspan="1">comp82651_c1</td><td align="left" rowspan="1" colspan="1">1.67</td><td align="left" rowspan="1" colspan="1">5.00E-94</td><td align="left" rowspan="1" colspan="1">1.49E-90</td><td align="left" rowspan="1" colspan="1">1.25</td><td align="left" rowspan="1" colspan="1">6.40E-67</td><td align="left" rowspan="1" colspan="1">1.12E-63</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt103"><label></label><p>Note: Fold = log<sub>2</sub>(fold change). Significant DE genes were determined based on |log<sub>2</sub>(fold change)|>1 and FDR q value<0.005.</p></fn></table-wrap-foot></table-wrap><p>Down-regulated oxidoreductase genes in this study comprise 8 iron/ascorbate family oxidoreductase genes, 1 MDHAR gene, 1 GPX gene and 5 GST genes, which explains the significant enrichment of “drug metabolism” in KEGG analysis of down-regulated genes at the MS stage (<xref ref-type="table" rid="pone-0112320-t004">Table 4</xref>). Within these genes, MDHAR gene had the maximum 17-fold decrease of its transcript abundance level in JA-CMS, which may be a major player of glutathione-ascorbate cycle that causes ROS damages. Thioredoxin 1 and thioredoxin reductase (NADPH) were down-regulated in JA-CMS as well (<xref ref-type="table" rid="pone-0112320-t004">Table 4</xref>). Although not directly involved in glutathione-ascorbate cycle, thioredoxins are involved in oxidative damage avoidance by supplying reducing power to reductases detoxifying lipid hydroperoxides or repairing oxidized proteins <xref rid="pone.0112320-DosSantos1" ref-type="bibr">[37]</xref>. Overall, increased production of ROS through photosynthesis plus decreased enzymatic activity of ROS-scavenging enzymes may dramatically perturb ROS homeostasis in JA-CMS.</p><table-wrap id="pone-0112320-t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0112320.t004</object-id><label>Table 4</label><caption><title>Down-regulated oxidoreductase and alpha-linolenic acid metabolism genes in JA-CMS at the SS and/or MS stages.</title></caption><alternatives><graphic id="pone-0112320-t004-4" xlink:href="pone.0112320.t004"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Function</td><td align="left" rowspan="1" colspan="1">Gene ID</td><td colspan="3" align="left" rowspan="1">SS stage</td><td colspan="3" align="left" rowspan="1">MS stage</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Fold</td><td align="left" rowspan="1" colspan="1">P value</td><td align="left" rowspan="1" colspan="1">Q value</td><td align="left" rowspan="1" colspan="1">Fold</td><td align="left" rowspan="1" colspan="1">P value</td><td align="left" rowspan="1" colspan="1">Q value</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Oxidoreductase</td><td align="left" rowspan="1" colspan="1">Iron/ascorbate family oxidoreductase</td><td align="left" rowspan="1" colspan="1">comp74308_c0</td><td align="left" rowspan="1" colspan="1">−2.17</td><td align="left" rowspan="1" colspan="1">1.52E-19</td><td align="left" rowspan="1" colspan="1">4.38E-17</td><td align="left" rowspan="1" colspan="1">−2.46</td><td align="left" rowspan="1" colspan="1">1.15E-30</td><td align="left" rowspan="1" colspan="1">6.70E-28</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp75973_c0</td><td align="left" rowspan="1" colspan="1">−1.29</td><td align="left" rowspan="1" colspan="1">1.58E-14</td><td align="left" rowspan="1" colspan="1">3.19E-12</td><td align="left" rowspan="1" colspan="1">−1.60</td><td align="left" rowspan="1" colspan="1">6.00E-28</td><td align="left" rowspan="1" colspan="1">2.92E-25</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp76619_c0</td><td align="left" rowspan="1" colspan="1">−1.97</td><td align="left" rowspan="1" colspan="1">2.58E-69</td><td align="left" rowspan="1" colspan="1">4.05E-66</td><td align="left" rowspan="1" colspan="1">−1.58</td><td align="left" rowspan="1" colspan="1">3.31E-40</td><td align="left" rowspan="1" colspan="1">2.66E-37</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp76821_c0</td><td align="left" rowspan="1" colspan="1">−1.55</td><td align="left" rowspan="1" colspan="1">6.12E-11</td><td align="left" rowspan="1" colspan="1">8.77E-09</td><td align="left" rowspan="1" colspan="1">−1.12</td><td align="left" rowspan="1" colspan="1">1.53E-06</td><td align="left" rowspan="1" colspan="1">1.34E-04</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp77552_c0</td><td align="left" rowspan="1" colspan="1">−3.93</td><td align="left" rowspan="1" colspan="1">8.85E-91</td><td align="left" rowspan="1" colspan="1">2.22E-87</td><td align="left" rowspan="1" colspan="1">−2.97</td><td align="left" rowspan="1" colspan="1">1.16E-42</td><td align="left" rowspan="1" colspan="1">1.06E-39</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp78972_c0</td><td align="left" rowspan="1" colspan="1">−4.08</td><td align="left" rowspan="1" colspan="1">8.28E-77</td><td align="left" rowspan="1" colspan="1">1.67E-73</td><td align="left" rowspan="1" colspan="1">−2.89</td><td align="left" rowspan="1" colspan="1">8.24E-43</td><td align="left" rowspan="1" colspan="1">7.64E-40</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp79641_c0</td><td align="left" rowspan="1" colspan="1">−1.96</td><td align="left" rowspan="1" colspan="1">5.14E-05</td><td align="left" rowspan="1" colspan="1">3.22E-03</td><td align="left" rowspan="1" colspan="1">−2.35</td><td align="left" rowspan="1" colspan="1">6.26E-07</td><td align="left" rowspan="1" colspan="1">5.81E-05</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp79893_c0</td><td align="left" rowspan="1" colspan="1">−1.16</td><td align="left" rowspan="1" colspan="1">1.52E-13</td><td align="left" rowspan="1" colspan="1">2.81E-11</td><td align="left" rowspan="1" colspan="1">−1.02</td><td align="left" rowspan="1" colspan="1">9.31E-12</td><td align="left" rowspan="1" colspan="1">1.54E-09</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Monodehydroascorbate reductase (MDHAR)</td><td align="left" rowspan="1" colspan="1">comp73498_c0</td><td align="left" rowspan="1" colspan="1">−4.21</td><td align="left" rowspan="1" colspan="1">1.43E-07</td><td align="left" rowspan="1" colspan="1">1.37E-05</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Glutathione peroxidase (GPX)</td><td align="left" rowspan="1" colspan="1">comp73626_c0</td><td align="left" rowspan="1" colspan="1">−1.45</td><td align="left" rowspan="1" colspan="1">2.69E-80</td><td align="left" rowspan="1" colspan="1">6.02E-77</td><td align="left" rowspan="1" colspan="1">−1.55</td><td align="left" rowspan="1" colspan="1">5.17E-97</td><td align="left" rowspan="1" colspan="1">1.58E-93</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Glutathione S-transferase (GST)</td><td align="left" rowspan="1" colspan="1">comp40493_c0</td><td align="left" rowspan="1" colspan="1">−2.71</td><td align="left" rowspan="1" colspan="1">1.54E-11</td><td align="left" rowspan="1" colspan="1">2.36E-09</td><td align="left" rowspan="1" colspan="1">−2.46</td><td align="left" rowspan="1" colspan="1">3.63E-13</td><td align="left" rowspan="1" colspan="1">6.88E-11</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp74280_c0</td><td align="left" rowspan="1" colspan="1">−1.66</td><td align="left" rowspan="1" colspan="1">3.00E-14</td><td align="left" rowspan="1" colspan="1">5.83E-12</td><td align="left" rowspan="1" colspan="1">−1.32</td><td align="left" rowspan="1" colspan="1">1.37E-10</td><td align="left" rowspan="1" colspan="1">2.00E-08</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp76601_c0</td><td align="left" rowspan="1" colspan="1">−1.85</td><td align="left" rowspan="1" colspan="1">3.65E-15</td><td align="left" rowspan="1" colspan="1">7.60E-13</td><td align="left" rowspan="1" colspan="1">−1.77</td><td align="left" rowspan="1" colspan="1">2.14E-16</td><td align="left" rowspan="1" colspan="1">5.29E-14</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp81581_c0</td><td align="left" rowspan="1" colspan="1">−1.54</td><td align="left" rowspan="1" colspan="1">1.40E-05</td><td align="left" rowspan="1" colspan="1">1.01E-03</td><td align="left" rowspan="1" colspan="1">−1.36</td><td align="left" rowspan="1" colspan="1">4.44E-05</td><td align="left" rowspan="1" colspan="1">2.95E-03</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp92611_c0</td><td align="left" rowspan="1" colspan="1">−2.96</td><td align="left" rowspan="1" colspan="1">1.61E-06</td><td align="left" rowspan="1" colspan="1">1.35E-04</td><td align="left" rowspan="1" colspan="1">−2.32</td><td align="left" rowspan="1" colspan="1">3.14E-06</td><td align="left" rowspan="1" colspan="1">2.60E-04</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Thioredoxin 1</td><td align="left" rowspan="1" colspan="1">comp63917_c0</td><td colspan="3" align="left" rowspan="1">Not significant</td><td align="left" rowspan="1" colspan="1">−1.57</td><td align="left" rowspan="1" colspan="1">4.86E-07</td><td align="left" rowspan="1" colspan="1">4.58E-05</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Thioredoxin reductase (NADPH)</td><td align="left" rowspan="1" colspan="1">comp87684_c1</td><td align="left" rowspan="1" colspan="1">−1.17</td><td align="left" rowspan="1" colspan="1">1.99E-12</td><td align="left" rowspan="1" colspan="1">3.35E-10</td><td align="left" rowspan="1" colspan="1">−1.20</td><td align="left" rowspan="1" colspan="1">1.21E-13</td><td align="left" rowspan="1" colspan="1">2.37E-11</td></tr><tr><td align="left" rowspan="1" colspan="1">Alpha-linolenic acid metabolism</td><td align="left" rowspan="1" colspan="1">Allene oxide synthase (AOS)</td><td align="left" rowspan="1" colspan="1">comp75493_c0</td><td align="left" rowspan="1" colspan="1">−1.45</td><td align="left" rowspan="1" colspan="1">1.04E-127</td><td align="left" rowspan="1" colspan="1">5.93E-124</td><td align="left" rowspan="1" colspan="1">−1.02</td><td align="left" rowspan="1" colspan="1">1.57E-58</td><td align="left" rowspan="1" colspan="1">2.34E-55</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Allene oxide cyclase (AOC)</td><td align="left" rowspan="1" colspan="1">comp70218_c0</td><td align="left" rowspan="1" colspan="1">−2.68</td><td align="left" rowspan="1" colspan="1">1.51E-10</td><td align="left" rowspan="1" colspan="1">2.09E-08</td><td align="left" rowspan="1" colspan="1">−2.03</td><td align="left" rowspan="1" colspan="1">1.63E-08</td><td align="left" rowspan="1" colspan="1">1.83E-06</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">comp80154_c0</td><td align="left" rowspan="1" colspan="1">−1.75</td><td align="left" rowspan="1" colspan="1">4.08E-59</td><td align="left" rowspan="1" colspan="1">5.11E-56</td><td align="left" rowspan="1" colspan="1">−1.61</td><td align="left" rowspan="1" colspan="1">7.73E-53</td><td align="left" rowspan="1" colspan="1">9.09E-50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">OPDA (12-oxophytodienoic acid) reductase</td><td align="left" rowspan="1" colspan="1">comp79354_c1</td><td align="left" rowspan="1" colspan="1">−2.83</td><td align="left" rowspan="1" colspan="1">8.17E-85</td><td align="left" rowspan="1" colspan="1">1.97E-81</td><td align="left" rowspan="1" colspan="1">−2.84</td><td align="left" rowspan="1" colspan="1">1.73E-71</td><td align="left" rowspan="1" colspan="1">3.53E-68</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">OPC-8:0 CoA ligase 1</td><td align="left" rowspan="1" colspan="1">comp82201_c0</td><td align="left" rowspan="1" colspan="1">−1.01</td><td align="left" rowspan="1" colspan="1">8.43E-29</td><td align="left" rowspan="1" colspan="1">3.95E-26</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Acyl-CoA oxidase (AOX)</td><td align="left" rowspan="1" colspan="1">comp86147_c0</td><td align="left" rowspan="1" colspan="1">−1.80</td><td align="left" rowspan="1" colspan="1">4.97E-105</td><td align="left" rowspan="1" colspan="1">1.64E-101</td><td align="left" rowspan="1" colspan="1">−1.63</td><td align="left" rowspan="1" colspan="1">4.80E-79</td><td align="left" rowspan="1" colspan="1">1.18E-75</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Acetyl-CoA acyltransferase 1 (ACAA)</td><td align="left" rowspan="1" colspan="1">comp80628_c2</td><td align="left" rowspan="1" colspan="1">-1.09</td><td align="left" rowspan="1" colspan="1">7.33E-17</td><td align="left" rowspan="1" colspan="1">1.73E-14</td><td colspan="3" align="left" rowspan="1">Not significant</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Alpha-dioxygenase</td><td align="left" rowspan="1" colspan="1">comp66758_c0</td><td align="left" rowspan="1" colspan="1">−5.66</td><td align="left" rowspan="1" colspan="1">2.09E-11</td><td align="left" rowspan="1" colspan="1">3.15E-09</td><td colspan="3" align="left" rowspan="1"></td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt104"><label></label><p>Note: Fold = log<sub>2</sub>(fold change). Significant DE genes were determined based on |log<sub>2</sub>(fold change)|>1 and FDR q value<0.005.</p></fn></table-wrap-foot></table-wrap><p>The “alpha-linolenic acid metabolism” pathway was another one that was significantly over-represented in KEGG analysis of down-regulated DE genes at the SS stage in JA-CMS. Although this pathway was only detected for the SS stage, genes encoding key enzymes in this pathway were significantly down-regulated in JA-CMS at both SS and/or MS stages, which include allene oxide synthase gene (AOS), 2 allene oxide cyclase genes (AOC), OPDA (12-oxo-phytodienoic acid) reductase gene, OPC-8:0 (3-oxo-2((2Z)-pentenyl)-cyclopentane-1-octanoic acid) CoA ligase gene, acyl-CoA oxidase gene (AOX), acetyl-CoA acyltransferase gene (ACAA), and α-dioxygenase gene (<xref ref-type="table" rid="pone-0112320-t004">Table 4</xref>). Among these genes, α-dioxygenase gene is worth noting, because its transcript abundance decreased around 50-fold in JA-CMS comparing with its counterpart at the SS stage; and it was reported to protect against oxidative stress and cell death in <italic>Arabidopsis</italic><xref rid="pone.0112320-DeLeon1" ref-type="bibr">[38]</xref>, lack of which may exacerbate ROS imbalance as aforementioned.</p><p>The “alpha-linolenic acid metabolism” pathway belongs to jasmonic acid biosynthesis module according to KEGG pathway database. Multiple studies have presented the association between jasmonic acid biosynthesis and CMS. McConn and Browse <xref rid="pone.0112320-McConn1" ref-type="bibr">[39]</xref> found that <italic>Arabidopsis</italic> triple mutants that contain negligible levels of trienoic fatty acids were male sterile, and exogenous applications of α-linoleate or jasmonate restored fertility. An <italic>Arabidopsis</italic> knock-out mutant defective in the AOS gene <italic>CYP74A</italic> also showed severe male sterility due to defects in anther and pollen development; this phenotype was completely rescued by exogenous application of methyl jasmonate and by complementation with constitutive expression of the AOS gene <xref rid="pone.0112320-Park1" ref-type="bibr">[40]</xref>. In addition, aberrant jasmonic acid pathway was detected in CMS rice during the development of microspores <xref rid="pone.0112320-Liu2" ref-type="bibr">[41]</xref>.</p><p>Besides photosynthesis related genes, the “flavonoid biosynthesis” pathway were over-represented with up-regulated DE genes in JA-CMS at the SS and MS stages as well. This pathway comprises 10 up-regulated genes in this study, which are 4 chalcone synthase genes, one of the key enzymes for flavonoid biosynthesis, 1 anthocyanidin reductase gene, 1 leucoanthocyanidin dioxygenase gene, 1 bifunctional dihydrofavonol 4-reductase/flavanone 4-reductase gene, 1 trans-cinnamate 4-monooxygenase gene, 1 flavonoid 3′, 5′-hydroxylase gene and 1 naringenin 3-dioxygenase gene (<xref ref-type="table" rid="pone-0112320-t003">Table 3</xref>). However, two flavonoid biosynthesis related genes were found down-regulated in JA-CMS simultaneously; they are 1 chalcone synthase gene with a 4-fold expression decrease (q = 2.21×10<sup>−3</sup>) at the SS stage, and 1 caffeoyl-CoA O-methyltransferase gene down-regulated at both SS (2-fold, q = 1.05×10<sup>−5</sup>) and MS stages (2-fold, q = 3.44×10<sup>−7</sup>).</p><p>In literature, several nuclear genes pertaining to flavonoid biosynthesis were reported to be inhibited in Ogura CMS of <italic>Raphanus sativus</italic>. The expression of chalcone synthase gene was strongly inhibited in the later stages of anther development in sterile cytoplasm <xref rid="pone.0112320-Yang1" ref-type="bibr">[42]</xref>. Wei et al. <xref rid="pone.0112320-Wei1" ref-type="bibr">[43]</xref> implemented digital gene expression analysis to compare gene expressions of wild type and genetic male sterility (GMS) mutant cotton at the meiosis, tetrad, and uninucleate microspore stages, respectively. Genes involved in flavonoid metabolism, such as CHS, flavonoid 3′, 5′-hydroxylase, anthocyanidin reductase, and leucoanthocyanidin reductase were down-regulated at the meiosis and uninucleate microspore stages, but were up-regulated at the tetrad stage, in GMS mutant anthers. Differences between the current study and these findings may lie in diverse plant materials, different sterility mechanisms and distinct development stages. Nevertheless contradictory results, flavonoids have been suggested to constitute a secondary ROS-scavenging system in plants exposed to severe stress conditions <xref rid="pone.0112320-Fini1" ref-type="bibr">[44]</xref> and genes in flavonoid biosynthesis are important in anther and pollen development. Further studies are warranted to disentangle specific functions of these genes in balancing ROS homeostasis at different development stages of sterile plants.</p><p>In addition, in GO analysis of DE genes between the two stages of the same material, ATP biosynthesis related categories were significantly enriched in sterile plants. Three out of the 9 down-regulated genes identified in this comparison were involved in ATP biosynthesis. The three DE genes, only detected in JA-CMS, are proton-transporting ATP synthase complex, coupling factor F<sub>o</sub> (comp73972_c0; 3-fold expression decrease at the MS stage, q = 5.88×10<sup>−85</sup>), ATPase subunit 1 (comp84388_c5; 3-fold, q = 1.58×10<sup>−5</sup>), and ATP synthase CF1 alpha subunit (comp85965_c1; 2-fold, q = 4.19×10<sup>−15</sup>). Same expression patterns of ATP synthase related genes were observed in studies of Dong et al. <xref rid="pone.0112320-Dong1" ref-type="bibr">[45]</xref> and Liu et al. <xref rid="pone.0112320-Liu1" ref-type="bibr">[15]</xref>.</p><p>ATP synthase, located within the plant mitochondria, consists of F<sub>o</sub> and F<sub>1</sub> regions, which plays an important role in providing energy. The F<sub>o</sub> portion is within the membrane and functions as a proton pore; it consists of orf25, orfB, subunit 6 (<italic>atp6</italic>) and subunit 9 (<italic>atp9</italic>) <xref rid="pone.0112320-Heazlewood1" ref-type="bibr">[46]</xref>, among which loci within or up/downstream <italic>atp6</italic> and <italic>atp9</italic> have been associated with CMS in a variety of plant species <xref rid="pone.0112320-Kempken1" ref-type="bibr">[47]</xref>–<xref rid="pone.0112320-Szklarczyk1" ref-type="bibr">[50]</xref>. Comp73972_c0 in this study may be one of these genes. There have been at least 16 CMS-related ORFs containing structural genes, especially the ones encoding ATP-F<sub>o</sub> and ATP-F<sub>1</sub><xref rid="pone.0112320-Chase1" ref-type="bibr">[2]</xref>, <xref rid="pone.0112320-Hanson1" ref-type="bibr">[51]</xref>–<xref rid="pone.0112320-Zhang1" ref-type="bibr">[54]</xref>. Cotton CMS studies have also demonstrated step by step that the <italic>atp1</italic> and <italic>atp6</italic> genes in CMS-D8 could be the candidates of CMS-associated genes in the mitochondrial genome <xref rid="pone.0112320-Wang3" ref-type="bibr">[55]</xref>, the abnormal sequence and expression of <italic>atpA</italic> gene is associated with CMS expression in Upland cotton <xref rid="pone.0112320-Wu1" ref-type="bibr">[56]</xref>, and PCR-based SNP markers in genes encoding ATP synthase subunits are useful in discriminating CMS-D8, CMS-D2 and Upland cotton cytoplasms <xref rid="pone.0112320-Suzuki2" ref-type="bibr">[57]</xref>. Because of the increased energy demand during pollen development, one or a few gene product(s) that interferes with mitochondrial F<sub>o</sub>F<sub>1</sub>-ATP synthase function is likely to induce pollen sterility <xref rid="pone.0112320-Li2" ref-type="bibr">[58]</xref>. But based on the results presented here, dysfunctional ATP synthase may not be the major cause of CMS in this study.</p><p>In summary, as the first study using unbiased transcriptome sequencing technology to investigate CMS-associated genes in cotton, a substantial number of biologically meaningful DE genes, especially the ones related with photosynthesis, redox reactions, alpha-linolenic acid metabolism, and flavonoid biosynthesis, were identified. Profiling expression patterns of these genes across all stages of pollen development may provide additional evidence for their importance. Further experiments are warranted to elucidate molecular mechanisms of these genes that lead to CMS.</p><sec id="s4a"><title>Supplementary Note</title><p>Several CMS systems have been reported in cotton which include <italic>Gossypiumar boreum</italic>, <italic>G. anomalum</italic>, <italic>G. harknessii</italic>, <italic>G. trilobum</italic>, <italic>G. barbadense</italic> and <italic>G. hirsutum</italic>. CMS-D2 and CMS-D8 are the two most commonly studied systems among them.</p><p>The male sterile characteristic of JA-CMS is maintained by <italic>G. hirsutum</italic>, and its <italic>restorer-of-fertility</italic> (<italic>Rf</italic>) genes come from <italic>G. barbadense</italic>. The primitive type and natural outcrossing hybrids of JA-CMS resemble its female parent (<italic>G. hirsutum</italic>) <xref rid="pone.0112320-Yuan1" ref-type="bibr">[59]</xref>. However, the original <italic>Rf</italic> genes of CMS-D2 and 104-7A come from <italic>G. harknessii</italic><xref rid="pone.0112320-Jia1" ref-type="bibr">[17]</xref>, <xref rid="pone.0112320-Wei2" ref-type="bibr">[60]</xref>. And CMS-D8 is known to have its <italic>Rf</italic> genes come from its sister line <xref rid="pone.0112320-Stewart1" ref-type="bibr">[61]</xref>.</p><p>Meanwhile, the mtDNA RFLP and PCR analysis results revealed that JA-CMS possesses a novel mitotype compared with the mitochondrial genome of either 104-7A or <italic>G. harknessii</italic> CMS systems <xref rid="pone.0112320-Huang2" ref-type="bibr">[62]</xref>, <xref rid="pone.0112320-Ma1" ref-type="bibr">[63]</xref>.</p><p>These evidence at genetic, morphological and molecular levels demonstrated that factors and mechanisms associated with CMS induction in JA-CMS may be different from CMS-D2, CMS-D8 and 104-7A.</p></sec></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0112320.s001"><label>Figure S1</label><caption><p>Results of qRT-PCR. Columns and bars represented the means and standard errors (n = 3), respectively. B2 = SS stage of JA-CMS, K2 = SS stage of JB, B3 = MS stage of JA-CMS, K3 = MS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s001.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s002"><label>Figure S2</label><caption><p>Transcript length distribution. Minimum, N90, mean, N50, and maximum lengths of transcripts are marked on the figure. The N90/N50 value is a weighted median and defined as the length of the smallest transcript S in the sorted list of all transcripts where the cumulative length from the largest transcript to S is at least 90%/50% of the total length. There are 206,496 transcripts in total; most of the transcripts (74,998 [36.3%]) have lengths between 200–500 base pairs.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s002.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s003"><label>Figure S3</label><caption><p>Unigene length distribution. Minimum, N90, mean, N50, and maximum lengths of unigenes are marked on the figure. The N90/N50 value is a weighted median and defined as the length of the smallest unigene S in the sorted list of all unigenes where the cumulative length from the largest unigene to S is at least 90%/50% of the total length. A total of 86,093 unigenes were identified; 59.8% (51,491) of the unigenes have lengths between 200–500 base pairs.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s003.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s004"><label>Figure S4</label><caption><p>RPKM density distributions of 854 DE genes. B2 = SS stage of JA-CMS, K2 = SS stage of JB, B3 = MS stage of JA-CMS, K3 = MS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s004.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s005"><label>Figure S5</label><caption><p>GO analysis results of up-regulated DE genes in JA-CMS at the SS stage comparing with JB. BP = Biological process, CC = Cellular component, MF = Molecular function; B2 = SS stage of JA-CMS, K2 = SS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s005.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s006"><label>Figure S6</label><caption><p>GO analysis results of down-regulated DE genes in JA-CMS at the SS stage comparing with JB. BP = Biological process, CC = Cellular component, MF = Molecular function; B2 = SS stage of JA-CMS, K2 = SS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s006.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s007"><label>Figure S7</label><caption><p>GO analysis results of up-regulated DE genes in JA-CMS at the MS stage comparing with JB. BP = Biological process, CC = Cellular component, MF = Molecular function; B3 = MS stage of JA-CMS, K3 = MS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s007.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s008"><label>Figure S8</label><caption><p>GO analysis results of down-regulated DE genes in JA-CMS at the MS stage comparing with JB. BP = Biological process, CC = Cellular component, MF = Molecular function; B3 = MS stage of JA-CMS, K3 = MS stage of JB.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s008.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s009"><label>Figure S9</label><caption><p>GO analysis results of down-regulated DE genes at the MS stage comparing with the SS stage in JA-CMS. BP = Biological process, CC = Cellular component, MF = Molecular function; B3 = MS stage of JA-CMS, B2 = SS stage of JA-CMS.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s009.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s010"><label>Figure S10</label><caption><p>KEGG analysis results of up-regulated DE genes in JA-CMS at the SS stage comparing with JB. Rich factor is the ratio between counts of DE genes and all annotated genes enriched in a certain pathway; qvalue is P value after multiple hypothesis testing correction with a range between 0 and 1. Twenty most significant pathways were plotted, when more than 20 pathways were identified.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s010.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s011"><label>Figure S11</label><caption><p>KEGG analysis results of down-regulated DE genes in JA-CMS at the SS stage comparing with JB. Rich factor is the ratio between counts of DE genes and all annotated genes enriched in a certain pathway; qvalue is P value after multiple hypothesis testing correction with a range between 0 and 1. Twenty most significant pathways were plotted, when more than 20 pathways were identified.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s011.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s012"><label>Figure S12</label><caption><p>KEGG analysis results of up-regulated DE genes in JA-CMS at the MS stage comparing with JB. Rich factor is the ratio between counts of DE genes and all annotated genes enriched in a certain pathway; qvalue is P value after multiple hypothesis testing correction with a range between 0 and 1. Twenty most significant pathways were plotted, when more than 20 pathways were identified.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s012.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s013"><label>Figure S13</label><caption><p>KEGG analysis results of down-regulated DE genes in JA-CMS at the MS stage comparing with JB. Rich factor is the ratio between counts of DE genes and all annotated genes enriched in a certain pathway; qvalue is P value after multiple hypothesis testing correction with a range between 0 and 1. Twenty most significant pathways were plotted, when more than 20 pathways were identified.</p><p>(TIF)</p></caption><media xlink:href="pone.0112320.s013.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s014"><label>Table S1</label><caption><p>Gene-specific primers for qRT-PCR.</p><p>(DOCX)</p></caption><media xlink:href="pone.0112320.s014.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s015"><label>Table S2</label><caption><p>DE genes between JA-CMS and JB plants at the SS stage.</p><p>(XLSX)</p></caption><media xlink:href="pone.0112320.s015.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s016"><label>Table S3</label><caption><p>DE genes between JA-CMS and JB plants at the MS stage.</p><p>(XLSX)</p></caption><media xlink:href="pone.0112320.s016.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s017"><label>Table S4</label><caption><p>DE genes between the SS and MS stages in JA-CMS.</p><p>(XLSX)</p></caption><media xlink:href="pone.0112320.s017.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0112320.s018"><label>Table S5</label><caption><p>DE genes between the SS and MS stages in JB.</p><p>(XLSX)</p></caption><media xlink:href="pone.0112320.s018.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ref-list><title>References</title><ref id="pone.0112320-Young1"><label>1</label><mixed-citation publication-type="journal"><name><surname>Young</surname><given-names>EG</given-names></name>, <name><surname>Hanson</surname><given-names>MR</given-names></name> (<year>1987</year>) <article-title>A fused mitochondrial gene associated with cytoplasmic male sterility is developmentally regulated</article-title>. <source>Cell</source><volume>50</volume>: <fpage>41</fpage>–<lpage>49</lpage>.<pub-id pub-id-type="pmid">2885095</pub-id></mixed-citation></ref><ref id="pone.0112320-Chase1"><label>2</label><mixed-citation publication-type="journal"><name><surname>Chase</surname><given-names>CD</given-names></name> (<year>2007</year>) <article-title>Cytoplasmic male sterility: a window to the world of plant mitochondrial-nuclear interactions</article-title>. <source>Trends Genet</source><volume>23</volume>: <fpage>81</fpage>–<lpage>90</lpage>.<pub-id pub-id-type="pmid">17188396</pub-id></mixed-citation></ref><ref id="pone.0112320-Linke1"><label>3</label><mixed-citation publication-type="journal"><name><surname>Linke</surname><given-names>B</given-names></name>, <name><surname>Borner</surname><given-names>T</given-names></name> (<year>2005</year>) <article-title>Mitochondrial effects on flower and pollen development</article-title>. <source>Mitochondrion</source><volume>5</volume>: <fpage>389</fpage>–<lpage>402</lpage>.<pub-id pub-id-type="pmid">16275170</pub-id></mixed-citation></ref><ref id="pone.0112320-Unseld1"><label>4</label><mixed-citation publication-type="journal"><name><surname>Unseld</surname><given-names>M</given-names></name>, <name><surname>Marienfeld</surname><given-names>JR</given-names></name>, <name><surname>Brandt</surname><given-names>P</given-names></name>, <name><surname>Brennicke</surname><given-names>A</given-names></name> (<year>1997</year>) <article-title>The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides</article-title>. <source>Nat Genet</source><volume>15</volume>: <fpage>57</fpage>–<lpage>61</lpage>.<pub-id pub-id-type="pmid">8988169</pub-id></mixed-citation></ref><ref id="pone.0112320-Kubo1"><label>5</label><mixed-citation publication-type="journal"><name><surname>Kubo</surname><given-names>T</given-names></name>, <name><surname>Nishizawa</surname><given-names>S</given-names></name>, <name><surname>Sugawara</surname><given-names>A</given-names></name>, <name><surname>Itchoda</surname><given-names>N</given-names></name>, <name><surname>Estiati</surname><given-names>A</given-names></name>, <etal>et al</etal> (<year>2000</year>) <article-title>The complete nucleotide sequence of the mitochondrial genome of sugar beet (Beta vulgaris L.) reveals a novel gene for tRNA(Cys)(GCA)</article-title>. <source>Nucleic Acids Res</source><volume>28</volume>: <fpage>2571</fpage>–<lpage>2576</lpage>.<pub-id pub-id-type="pmid">10871408</pub-id></mixed-citation></ref><ref id="pone.0112320-Notsu1"><label>6</label><mixed-citation publication-type="journal"><name><surname>Notsu</surname><given-names>Y</given-names></name>, <name><surname>Masood</surname><given-names>S</given-names></name>, <name><surname>Nishikawa</surname><given-names>T</given-names></name>, <name><surname>Kubo</surname><given-names>N</given-names></name>, <name><surname>Akiduki</surname><given-names>G</given-names></name>, <etal>et al</etal> (<year>2002</year>) <article-title>The complete sequence of the rice (Oryza sativa L.) mitochondrial genome: frequent DNA sequence acquisition and loss during the evolution of flowering plants</article-title>. <source>Mol Genet Genomics</source><volume>268</volume>: <fpage>434</fpage>–<lpage>445</lpage>.<pub-id pub-id-type="pmid">12471441</pub-id></mixed-citation></ref><ref id="pone.0112320-Sugiyama1"><label>7</label><mixed-citation publication-type="journal"><name><surname>Sugiyama</surname><given-names>Y</given-names></name>, <name><surname>Watase</surname><given-names>Y</given-names></name>, <name><surname>Nagase</surname><given-names>M</given-names></name>, <name><surname>Makita</surname><given-names>N</given-names></name>, <name><surname>Yagura</surname><given-names>S</given-names></name>, <etal>et al</etal> (<year>2005</year>) <article-title>The complete nucleotide sequence and multipartite organization of the tobacco mitochondrial genome: comparative analysis of mitochondrial genomes in higher plants</article-title>. <source>Mol Genet Genomics</source><volume>272</volume>: <fpage>603</fpage>–<lpage>615</lpage>.<pub-id pub-id-type="pmid">15583938</pub-id></mixed-citation></ref><ref id="pone.0112320-Rhoads1"><label>8</label><mixed-citation publication-type="journal"><name><surname>Rhoads</surname><given-names>DM</given-names></name>, <name><surname>Subbaiah</surname><given-names>CC</given-names></name> (<year>2007</year>) <article-title>Mitochondrial retrograde regulation in plants</article-title>. <source>Mitochondrion</source><volume>7</volume>: <fpage>177</fpage>–<lpage>194</lpage>.<pub-id pub-id-type="pmid">17320492</pub-id></mixed-citation></ref><ref id="pone.0112320-Suzuki1"><label>9</label><mixed-citation publication-type="journal"><name><surname>Suzuki</surname><given-names>H</given-names></name>, <name><surname>Rodriguez-Uribe</surname><given-names>L</given-names></name>, <name><surname>Xu</surname><given-names>J</given-names></name>, <name><surname>Zhang</surname><given-names>J</given-names></name> (<year>2013</year>) <article-title>Transcriptome analysis of cytoplasmic male sterility and restoration in CMS-D8 cotton</article-title>. <source>Plant Cell Rep</source><volume>32</volume>: <fpage>1531</fpage>–<lpage>1542</lpage>.<pub-id pub-id-type="pmid">23743655</pub-id></mixed-citation></ref><ref id="pone.0112320-Schuster1"><label>10</label><mixed-citation publication-type="journal"><name><surname>Schuster</surname><given-names>SC</given-names></name> (<year>2008</year>) <article-title>Next-generation sequencing transforms today's biology</article-title>. <source>Nat Methods</source><volume>5</volume>: <fpage>16</fpage>–<lpage>18</lpage>.<pub-id pub-id-type="pmid">18165802</pub-id></mixed-citation></ref><ref id="pone.0112320-Ansorge1"><label>11</label><mixed-citation publication-type="journal"><name><surname>Ansorge</surname><given-names>WJ</given-names></name> (<year>2009</year>) <article-title>Next-generation DNA sequencing techniques</article-title>. <source>N Biotechnol</source><volume>25</volume>: <fpage>195</fpage>–<lpage>203</lpage>.<pub-id pub-id-type="pmid">19429539</pub-id></mixed-citation></ref><ref id="pone.0112320-Wang1"><label>12</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>Z</given-names></name>, <name><surname>Gerstein</surname><given-names>M</given-names></name>, <name><surname>Snyder</surname><given-names>M</given-names></name> (<year>2009</year>) <article-title>RNA-Seq: a revolutionary tool for transcriptomics</article-title>. <source>Nat Rev Genet</source><volume>10</volume>: <fpage>57</fpage>–<lpage>63</lpage>.<pub-id pub-id-type="pmid">19015660</pub-id></mixed-citation></ref><ref id="pone.0112320-Zheng1"><label>13</label><mixed-citation publication-type="journal"><name><surname>Zheng</surname><given-names>BB</given-names></name>, <name><surname>Wu</surname><given-names>XM</given-names></name>, <name><surname>Ge</surname><given-names>XX</given-names></name>, <name><surname>Deng</surname><given-names>XX</given-names></name>, <name><surname>Grosser</surname><given-names>JW</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Comparative transcript profiling of a male sterile cybrid pummelo and its fertile type revealed altered gene expression related to flower development</article-title>. <source>PLoS One</source><volume>7</volume>: <fpage>e43758</fpage>.<pub-id pub-id-type="pmid">22952758</pub-id></mixed-citation></ref><ref id="pone.0112320-Yan1"><label>14</label><mixed-citation publication-type="journal"><name><surname>Yan</surname><given-names>X</given-names></name>, <name><surname>Dong</surname><given-names>C</given-names></name>, <name><surname>Yu</surname><given-names>J</given-names></name>, <name><surname>Liu</surname><given-names>W</given-names></name>, <name><surname>Jiang</surname><given-names>C</given-names></name>, <etal>et al</etal> (<year>2013</year>) <article-title>Transcriptome profile analysis of young floral buds of fertile and sterile plants from the self-pollinated offspring of the hybrid between novel restorer line NR1 and Nsa CMS line in Brassica napus</article-title>. <source>BMC Genomics</source><volume>14</volume>: <fpage>26</fpage>.<pub-id pub-id-type="pmid">23324545</pub-id></mixed-citation></ref><ref id="pone.0112320-Liu1"><label>15</label><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>C</given-names></name>, <name><surname>Ma</surname><given-names>N</given-names></name>, <name><surname>Wang</surname><given-names>PY</given-names></name>, <name><surname>Fu</surname><given-names>N</given-names></name>, <name><surname>Shen</surname><given-names>HL</given-names></name> (<year>2013</year>) <article-title>Transcriptome sequencing and De Novo analysis of a cytoplasmic male sterile line and its near-isogenic restorer line in chili pepper (Capsicum annuum L.)</article-title>. <source>PLoS One</source><volume>8</volume>: <fpage>e65209</fpage>.<pub-id pub-id-type="pmid">23750245</pub-id></mixed-citation></ref><ref id="pone.0112320-Huang1"><label>16</label><mixed-citation publication-type="journal"><name><surname>Huang</surname><given-names>J</given-names></name>, <name><surname>Yang</surname><given-names>P</given-names></name>, <name><surname>Li</surname><given-names>B</given-names></name> (<year>2001</year>) <article-title>The Microstructural and Ultrastructural Study on Microsporogenesis of Cytoplasmic Male Sterile Cotton Line Jin A. Cotton Science</article-title>. <volume>13</volume>: <fpage>259</fpage>–<lpage>263</lpage>.</mixed-citation></ref><ref id="pone.0112320-Jia1"><label>17</label><mixed-citation publication-type="journal"><name><surname>Jia</surname><given-names>ZC</given-names></name> (<year>1990</year>) <article-title>Selection and cultivation of CMS 104-7A cotton and its maintainer and restorer line</article-title>. <source>China Cotton</source><volume>17</volume>: <fpage>11</fpage>.</mixed-citation></ref><ref id="pone.0112320-Zhao1"><label>18</label><mixed-citation publication-type="journal"><name><surname>Zhao</surname><given-names>H</given-names></name>, <name><surname>Huang</surname><given-names>J</given-names></name> (<year>2012</year>) <article-title>Study on Microspore Abortion of Male Sterile Cotton Yamian A and Yamian B. Scientia Agricultura Sinica</article-title>. <volume>45</volume>: <fpage>4130</fpage>–<lpage>4140</lpage>.</mixed-citation></ref><ref id="pone.0112320-Hou1"><label>19</label><mixed-citation publication-type="other">Hou L (2001) Detection and isolation of genes differentially expressed during anther development of cotton by cDNA-AFLP.Chongqing: Southwestern Agricultural University.</mixed-citation></ref><ref id="pone.0112320-Grabherr1"><label>20</label><mixed-citation publication-type="journal"><name><surname>Grabherr</surname><given-names>MG</given-names></name>, <name><surname>Haas</surname><given-names>BJ</given-names></name>, <name><surname>Yassour</surname><given-names>M</given-names></name>, <name><surname>Levin</surname><given-names>JZ</given-names></name>, <name><surname>Thompson</surname><given-names>DA</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Full-length transcriptome assembly from RNA-Seq data without a reference genome</article-title>. <source>Nat Biotechnol</source><volume>29</volume>: <fpage>644</fpage>–<lpage>652</lpage>.<pub-id pub-id-type="pmid">21572440</pub-id></mixed-citation></ref><ref id="pone.0112320-Li1"><label>21</label><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>B</given-names></name>, <name><surname>Dewey</surname><given-names>CN</given-names></name> (<year>2011</year>) <article-title>RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome</article-title>. <source>BMC Bioinformatics</source><volume>12</volume>: <fpage>323</fpage>.<pub-id pub-id-type="pmid">21816040</pub-id></mixed-citation></ref><ref id="pone.0112320-Robinson1"><label>22</label><mixed-citation publication-type="journal"><name><surname>Robinson</surname><given-names>MD</given-names></name>, <name><surname>Oshlack</surname><given-names>A</given-names></name> (<year>2010</year>) <article-title>A scaling normalization method for differential expression analysis of RNA-seq data</article-title>. <source>Genome Biol</source><volume>11</volume>: <fpage>R25</fpage>.<pub-id pub-id-type="pmid">20196867</pub-id></mixed-citation></ref><ref id="pone.0112320-Wang2"><label>23</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>L</given-names></name>, <name><surname>Feng</surname><given-names>Z</given-names></name>, <name><surname>Wang</surname><given-names>X</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name> (<year>2010</year>) <article-title>DEGseq: an R package for identifying differentially expressed genes from RNA-seq data</article-title>. <source>Bioinformatics</source><volume>26</volume>: <fpage>136</fpage>–<lpage>138</lpage>.<pub-id pub-id-type="pmid">19855105</pub-id></mixed-citation></ref><ref id="pone.0112320-Young2"><label>24</label><mixed-citation publication-type="journal"><name><surname>Young</surname><given-names>MD</given-names></name>, <name><surname>Wakefield</surname><given-names>MJ</given-names></name>, <name><surname>Smyth</surname><given-names>GK</given-names></name>, <name><surname>Oshlack</surname><given-names>A</given-names></name> (<year>2010</year>) <article-title>Gene ontology analysis for RNA-seq: accounting for selection bias</article-title>. <source>Genome Biol</source><volume>11</volume>: <fpage>R14</fpage>.<pub-id pub-id-type="pmid">20132535</pub-id></mixed-citation></ref><ref id="pone.0112320-Xie1"><label>25</label><mixed-citation publication-type="journal"><name><surname>Xie</surname><given-names>C</given-names></name>, <name><surname>Mao</surname><given-names>X</given-names></name>, <name><surname>Huang</surname><given-names>J</given-names></name>, <name><surname>Ding</surname><given-names>Y</given-names></name>, <name><surname>Wu</surname><given-names>J</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases</article-title>. <source>Nucleic Acids Res</source><volume>39</volume>: <fpage>W316</fpage>–<lpage>322</lpage>.<pub-id pub-id-type="pmid">21715386</pub-id></mixed-citation></ref><ref id="pone.0112320-Livak1"><label>26</label><mixed-citation publication-type="journal"><name><surname>Livak</surname><given-names>KJ</given-names></name>, <name><surname>Schmittgen</surname><given-names>TD</given-names></name> (<year>2001</year>) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. <article-title>Methods</article-title>. <volume>25</volume>: <fpage>402</fpage>–<lpage>408</lpage>.</mixed-citation></ref><ref id="pone.0112320-Subramanian1"><label>27</label><mixed-citation publication-type="journal"><name><surname>Subramanian</surname><given-names>A</given-names></name>, <name><surname>Tamayo</surname><given-names>P</given-names></name>, <name><surname>Mootha</surname><given-names>VK</given-names></name>, <name><surname>Mukherjee</surname><given-names>S</given-names></name>, <name><surname>Ebert</surname><given-names>BL</given-names></name>, <etal>et al</etal> (<year>2005</year>) <article-title>Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles</article-title>. <source>Proc Natl Acad Sci U S A</source><volume>102</volume>: <fpage>15545</fpage>–<lpage>15550</lpage>.<pub-id pub-id-type="pmid">16199517</pub-id></mixed-citation></ref><ref id="pone.0112320-Balk1"><label>28</label><mixed-citation publication-type="journal"><name><surname>Balk</surname><given-names>J</given-names></name>, <name><surname>Leaver</surname><given-names>CJ</given-names></name> (<year>2001</year>) <article-title>The PET1-CMS mitochondrial mutation in sunflower is associated with premature programmed cell death and cytochrome c release</article-title>. <source>Plant Cell</source><volume>13</volume>: <fpage>1803</fpage>–<lpage>1818</lpage>.<pub-id pub-id-type="pmid">11487694</pub-id></mixed-citation></ref><ref id="pone.0112320-Diamond1"><label>29</label><mixed-citation publication-type="other">Diamond M, McCabe P (2011) Mitochondrial Regulation of Plant Programmed Cell Death. In: Kempken F, editor. Plant Mitochondria: Springer New York. 439–465.</mixed-citation></ref><ref id="pone.0112320-Papini1"><label>30</label><mixed-citation publication-type="journal"><name><surname>Papini</surname><given-names>A</given-names></name>, <name><surname>Mosti</surname><given-names>S</given-names></name>, <name><surname>Brighigna</surname><given-names>L</given-names></name> (<year>1999</year>) <article-title>Programmed-cell death events during tapetum development of angiosperms</article-title>. <source>Protoplasma</source><volume>207</volume>: <fpage>213</fpage>–<lpage>221</lpage>.</mixed-citation></ref><ref id="pone.0112320-Doyle1"><label>31</label><mixed-citation publication-type="journal"><name><surname>Doyle</surname><given-names>SM</given-names></name>, <name><surname>Diamond</surname><given-names>M</given-names></name>, <name><surname>McCabe</surname><given-names>PF</given-names></name> (<year>2010</year>) <article-title>Chloroplast and reactive oxygen species involvement in apoptotic-like programmed cell death in Arabidopsis suspension cultures</article-title>. <source>Journal of Experimental Botany</source><volume>61</volume>: <fpage>473</fpage>–<lpage>482</lpage>.<pub-id pub-id-type="pmid">19933317</pub-id></mixed-citation></ref><ref id="pone.0112320-Jiang1"><label>32</label><mixed-citation publication-type="journal"><name><surname>Jiang</surname><given-names>PD</given-names></name>, <name><surname>Zhang</surname><given-names>XQ</given-names></name>, <name><surname>Zhu</surname><given-names>YG</given-names></name>, <name><surname>Zhu</surname><given-names>W</given-names></name>, <name><surname>Xie</surname><given-names>HY</given-names></name>, <etal>et al</etal> (<year>2007</year>) <article-title>Metabolism of reactive oxygen species in cotton cytoplasmic male sterility and its restoration</article-title>. <source>Plant Cell Rep</source><volume>26</volume>: <fpage>1627</fpage>–<lpage>1634</lpage>.<pub-id pub-id-type="pmid">17426978</pub-id></mixed-citation></ref><ref id="pone.0112320-Wan1"><label>33</label><mixed-citation publication-type="journal"><name><surname>Wan</surname><given-names>CX</given-names></name>, <name><surname>Li</surname><given-names>SQ</given-names></name>, <name><surname>Wen</surname><given-names>L</given-names></name>, <name><surname>Kong</surname><given-names>J</given-names></name>, <name><surname>Wang</surname><given-names>K</given-names></name>, <etal>et al</etal> (<year>2007</year>) <article-title>Damage of oxidative stress on mitochondria during microspores development in Honglian CMS line of rice</article-title>. <source>Plant Cell Rep</source><volume>26</volume>: <fpage>373</fpage>–<lpage>382</lpage>.<pub-id pub-id-type="pmid">17053903</pub-id></mixed-citation></ref><ref id="pone.0112320-Gechev1"><label>34</label><mixed-citation publication-type="journal"><name><surname>Gechev</surname><given-names>TS</given-names></name>, <name><surname>Van Breusegem</surname><given-names>F</given-names></name>, <name><surname>Stone</surname><given-names>JM</given-names></name>, <name><surname>Denev</surname><given-names>I</given-names></name>, <name><surname>Laloi</surname><given-names>C</given-names></name> (<year>2006</year>) <article-title>Reactive oxygen species as signals that modulate plant stress responses and programmed cell death</article-title>. <source>Bioessays</source><volume>28</volume>: <fpage>1091</fpage>–<lpage>1101</lpage>.<pub-id pub-id-type="pmid">17041898</pub-id></mixed-citation></ref><ref id="pone.0112320-Moller1"><label>35</label><mixed-citation publication-type="journal"><name><surname>Moller</surname><given-names>IM</given-names></name> (<year>2001</year>) <article-title>Plant mitochondria and oxidative stress: Electron transport, NADPH turnover, and metabolism of reactive oxygen species</article-title>. <source>Annual Review of Plant Physiology and Plant Molecular Biology</source><volume>52</volume>: <fpage>561</fpage>–<lpage>591</lpage>.</mixed-citation></ref><ref id="pone.0112320-Sharma1"><label>36</label><mixed-citation publication-type="journal"><name><surname>Sharma</surname><given-names>P</given-names></name>, <name><surname>Jha</surname><given-names>AB</given-names></name>, <name><surname>Dubey</surname><given-names>RS</given-names></name>, <name><surname>Pessarakli</surname><given-names>M</given-names></name> (<year>2012</year>) <article-title>Reactive Oxygen Species, Oxidative Damage, and Antioxidative Defense Mechanism in Plants under Stressful Conditions</article-title>. <source>Journal of Botany</source><volume>2012</volume>: <fpage>26</fpage>.</mixed-citation></ref><ref id="pone.0112320-DosSantos1"><label>37</label><mixed-citation publication-type="journal"><name><surname>Dos Santos</surname><given-names>CV</given-names></name>, <name><surname>Rey</surname><given-names>P</given-names></name> (<year>2006</year>) <article-title>Plant thioredoxins are key actors in the oxidative stress response</article-title>. <source>Trends in Plant Science</source><volume>11</volume>: <fpage>329</fpage>–<lpage>334</lpage>.<pub-id pub-id-type="pmid">16782394</pub-id></mixed-citation></ref><ref id="pone.0112320-DeLeon1"><label>38</label><mixed-citation publication-type="journal"><name><surname>De Leon</surname><given-names>IP</given-names></name>, <name><surname>Sanz</surname><given-names>A</given-names></name>, <name><surname>Hamberg</surname><given-names>M</given-names></name>, <name><surname>Castresana</surname><given-names>C</given-names></name> (<year>2002</year>) <article-title>Involvement of the Arabidopsis alpha-DOX1 fatty acid dioxygenase in protection against oxidative stress and cell death</article-title>. <source>Plant J</source><volume>29</volume>: <fpage>61</fpage>–<lpage>62</lpage>.<pub-id pub-id-type="pmid">12060227</pub-id></mixed-citation></ref><ref id="pone.0112320-McConn1"><label>39</label><mixed-citation publication-type="journal"><name><surname>McConn</surname><given-names>M</given-names></name>, <name><surname>Browse</surname><given-names>J</given-names></name> (<year>1996</year>) <article-title>The critical requirement for linolenic acid is pollen development, not photosynthesis, in an arabidopsis mutant</article-title>. <source>Plant Cell</source><volume>8</volume>: <fpage>403</fpage>–<lpage>416</lpage>.<pub-id pub-id-type="pmid">12239389</pub-id></mixed-citation></ref><ref id="pone.0112320-Park1"><label>40</label><mixed-citation publication-type="journal"><name><surname>Park</surname><given-names>JH</given-names></name>, <name><surname>Halitschke</surname><given-names>R</given-names></name>, <name><surname>Kim</surname><given-names>HB</given-names></name>, <name><surname>Baldwin</surname><given-names>IT</given-names></name>, <name><surname>Feldmann</surname><given-names>KA</given-names></name>, <etal>et al</etal> (<year>2002</year>) <article-title>A knock-out mutation in allene oxide synthase results in male sterility and defective wound signal transduction in Arabidopsis due to a block in jasmonic acid biosynthesis</article-title>. <source>Plant Journal</source><volume>31</volume>: <fpage>1</fpage>–<lpage>12</lpage>.<pub-id pub-id-type="pmid">12100478</pub-id></mixed-citation></ref><ref id="pone.0112320-Liu2"><label>41</label><mixed-citation publication-type="other">Liu G, Tian H, Huang YQ, Hu J, Ji YX, et al. (2012) Alterations of Mitochondrial Protein Assembly and Jasmonic Acid Biosynthesis Pathway in Honglian (HL)-type Cytoplasmic Male Sterility Rice. Journal of Biological Chemistry 287.</mixed-citation></ref><ref id="pone.0112320-Yang1"><label>42</label><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>S</given-names></name>, <name><surname>Terachi</surname><given-names>T</given-names></name>, <name><surname>Yamagishi</surname><given-names>H</given-names></name> (<year>2008</year>) <article-title>Inhibition of chalcone synthase expression in anthers of Raphanus sativus with Ogura male sterile cytoplasm</article-title>. <source>Annals of Botany</source><volume>102</volume>: <fpage>483</fpage>–<lpage>489</lpage>.<pub-id pub-id-type="pmid">18625698</pub-id></mixed-citation></ref><ref id="pone.0112320-Wei1"><label>43</label><mixed-citation publication-type="journal"><name><surname>Wei</surname><given-names>M</given-names></name>, <name><surname>Song</surname><given-names>M</given-names></name>, <name><surname>Fan</surname><given-names>S</given-names></name>, <name><surname>Yu</surname><given-names>S</given-names></name> (<year>2013</year>) <article-title>Transcriptomic analysis of differentially expressed genes during anther development in genetic male sterile and wild type cotton by digital gene-expression profiling</article-title>. <source>BMC Genomics</source><volume>14</volume>: <fpage>97</fpage>.<pub-id pub-id-type="pmid">23402279</pub-id></mixed-citation></ref><ref id="pone.0112320-Fini1"><label>44</label><mixed-citation publication-type="journal"><name><surname>Fini</surname><given-names>A</given-names></name>, <name><surname>Brunetti</surname><given-names>C</given-names></name>, <name><surname>Di Ferdinando</surname><given-names>M</given-names></name>, <name><surname>Ferrini</surname><given-names>F</given-names></name>, <name><surname>Tattini</surname><given-names>M</given-names></name> (<year>2011</year>) <article-title>Stress-induced flavonoid biosynthesis and the antioxidant machinery of plants</article-title>. <source>Plant Signal Behav</source><volume>6</volume>: <fpage>709</fpage>–<lpage>711</lpage>.<pub-id pub-id-type="pmid">21448007</pub-id></mixed-citation></ref><ref id="pone.0112320-Dong1"><label>45</label><mixed-citation publication-type="journal"><name><surname>Dong</surname><given-names>X</given-names></name>, <name><surname>Feng</surname><given-names>H</given-names></name>, <name><surname>Xu</surname><given-names>M</given-names></name>, <name><surname>Lee</surname><given-names>J</given-names></name>, <name><surname>Kim</surname><given-names>YK</given-names></name>, <etal>et al</etal> (<year>2013</year>) <article-title>Comprehensive analysis of genic male sterility-related genes in Brassica rapa using a newly developed Br300K oligomeric chip</article-title>. <source>PLoS One</source><volume>8</volume>: <fpage>e72178</fpage>.<pub-id pub-id-type="pmid">24039743</pub-id></mixed-citation></ref><ref id="pone.0112320-Heazlewood1"><label>46</label><mixed-citation publication-type="journal"><name><surname>Heazlewood</surname><given-names>JL</given-names></name>, <name><surname>Whelan</surname><given-names>J</given-names></name>, <name><surname>Millar</surname><given-names>AH</given-names></name> (<year>2003</year>) <article-title>The products of the mitochondrial orf25 and orfB genes are FO components in the plant F1FO ATP synthase</article-title>. <source>FEBS Lett</source><volume>540</volume>: <fpage>201</fpage>–<lpage>205</lpage>.<pub-id pub-id-type="pmid">12681508</pub-id></mixed-citation></ref><ref id="pone.0112320-Kempken1"><label>47</label><mixed-citation publication-type="journal"><name><surname>Kempken</surname><given-names>F</given-names></name>, <name><surname>Howard</surname><given-names>W</given-names></name>, <name><surname>Pring</surname><given-names>DR</given-names></name> (<year>1998</year>) <article-title>Mutations at specific atp6 codons which cause human mitochondrial diseases also lead to male sterility in a plant</article-title>. <source>FEBS Lett</source><volume>441</volume>: <fpage>159</fpage>–<lpage>160</lpage>.<pub-id pub-id-type="pmid">9883875</pub-id></mixed-citation></ref><ref id="pone.0112320-delaCanal1"><label>48</label><mixed-citation publication-type="journal"><name><surname>de la Canal</surname><given-names>L</given-names></name>, <name><surname>Crouzillat</surname><given-names>D</given-names></name>, <name><surname>Quetier</surname><given-names>F</given-names></name>, <name><surname>Ledoigt</surname><given-names>G</given-names></name> (<year>2001</year>) <article-title>A transcriptional alteration on the atp9 gene is associated with a sunflower male-sterile cytoplasm</article-title>. <source>Theoretical and Applied Genetics</source><volume>102</volume>: <fpage>1185</fpage>–<lpage>1189</lpage>.</mixed-citation></ref><ref id="pone.0112320-Kurek1"><label>49</label><mixed-citation publication-type="journal"><name><surname>Kurek</surname><given-names>I</given-names></name>, <name><surname>Ezra</surname><given-names>D</given-names></name>, <name><surname>Begu</surname><given-names>D</given-names></name>, <name><surname>Erel</surname><given-names>N</given-names></name>, <name><surname>Litvak</surname><given-names>S</given-names></name>, <etal>et al</etal> (<year>1997</year>) <article-title>Studies on the effects of nuclear background and tissue specificity on RNA editing of the mitochondrial ATP synthase subunits alpha, 6 and 9 in fertile and cytoplasmic male-sterile (CMS) wheat</article-title>. <source>Theoretical and Applied Genetics</source><volume>95</volume>: <fpage>1305</fpage>–<lpage>1311</lpage>.</mixed-citation></ref><ref id="pone.0112320-Szklarczyk1"><label>50</label><mixed-citation publication-type="journal"><name><surname>Szklarczyk</surname><given-names>M</given-names></name>, <name><surname>Oczkowski</surname><given-names>M</given-names></name>, <name><surname>Augustyniak</surname><given-names>H</given-names></name>, <name><surname>Borner</surname><given-names>T</given-names></name>, <name><surname>Linke</surname><given-names>B</given-names></name>, <etal>et al</etal> (<year>2000</year>) <article-title>Organisation and expression of mitochondrial atp9 genes from CMS and fertile carrots</article-title>. <source>Theoretical and Applied Genetics</source><volume>100</volume>: <fpage>263</fpage>–<lpage>270</lpage>.</mixed-citation></ref><ref id="pone.0112320-Hanson1"><label>51</label><mixed-citation publication-type="journal"><name><surname>Hanson</surname><given-names>MR</given-names></name>, <name><surname>Bentolila</surname><given-names>S</given-names></name> (<year>2004</year>) <article-title>Interactions of mitochondrial and nuclear genes that affect male gametophyte development</article-title>. <source>Plant Cell</source><volume>16</volume><comment>Suppl: S154–169</comment></mixed-citation></ref><ref id="pone.0112320-Kim1"><label>52</label><mixed-citation publication-type="journal"><name><surname>Kim</surname><given-names>DH</given-names></name>, <name><surname>Kang</surname><given-names>JG</given-names></name>, <name><surname>Kim</surname><given-names>BD</given-names></name> (<year>2007</year>) <article-title>Isolation and characterization of the cytoplasmic male sterility-associated orf456 gene of chili pepper (Capsicum annuum L.)</article-title>. <source>Plant Mol Biol</source><volume>63</volume>: <fpage>519</fpage>–<lpage>532</lpage>.<pub-id pub-id-type="pmid">17238047</pub-id></mixed-citation></ref><ref id="pone.0112320-Yamamoto1"><label>53</label><mixed-citation publication-type="journal"><name><surname>Yamamoto</surname><given-names>MP</given-names></name>, <name><surname>Kubo</surname><given-names>T</given-names></name>, <name><surname>Mikami</surname><given-names>T</given-names></name> (<year>2005</year>) <article-title>The 5′-leader sequence of sugar beet mitochondrial atp6 encodes a novel polypeptide that is characteristic of Owen cytoplasmic male sterility</article-title>. <source>Mol Genet Genomics</source><volume>273</volume>: <fpage>342</fpage>–<lpage>349</lpage>.<pub-id pub-id-type="pmid">15838639</pub-id></mixed-citation></ref><ref id="pone.0112320-Zhang1"><label>54</label><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>H</given-names></name>, <name><surname>Li</surname><given-names>S</given-names></name>, <name><surname>Yi</surname><given-names>P</given-names></name>, <name><surname>Wan</surname><given-names>C</given-names></name>, <name><surname>Chen</surname><given-names>Z</given-names></name>, <etal>et al</etal> (<year>2007</year>) <article-title>A Honglian CMS line of rice displays aberrant F0 of F0F1-ATPase</article-title>. <source>Plant Cell Rep</source><volume>26</volume>: <fpage>1065</fpage>–<lpage>1071</lpage>.<pub-id pub-id-type="pmid">17226056</pub-id></mixed-citation></ref><ref id="pone.0112320-Wang3"><label>55</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>F</given-names></name>, <name><surname>Feng</surname><given-names>CD</given-names></name>, <name><surname>O'Connell</surname><given-names>MA</given-names></name>, <name><surname>Stewart</surname><given-names>JM</given-names></name>, <name><surname>Zhang</surname><given-names>JF</given-names></name> (<year>2010</year>) <article-title>RFLP analysis of mitochondrial DNA in two cytoplasmic male sterility systems (CMS-D2 and CMS-D8) of cotton</article-title>. <source>Euphytica</source><volume>172</volume>: <fpage>93</fpage>–<lpage>99</lpage>.</mixed-citation></ref><ref id="pone.0112320-Wu1"><label>56</label><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>JY</given-names></name>, <name><surname>Gong</surname><given-names>YC</given-names></name>, <name><surname>Cui</surname><given-names>MH</given-names></name>, <name><surname>Qi</surname><given-names>TX</given-names></name>, <name><surname>Guo</surname><given-names>LP</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Molecular characterization of cytoplasmic male sterility conditioned by Gossypium harknessii cytoplasm (CMS-D2) in upland cotton</article-title>. <source>Euphytica</source><volume>181</volume>: <fpage>17</fpage>–<lpage>29</lpage>.</mixed-citation></ref><ref id="pone.0112320-Suzuki2"><label>57</label><mixed-citation publication-type="journal"><name><surname>Suzuki</surname><given-names>H</given-names></name>, <name><surname>Yu</surname><given-names>JW</given-names></name>, <name><surname>Wang</surname><given-names>F</given-names></name>, <name><surname>Zhang</surname><given-names>JF</given-names></name> (<year>2013</year>) <article-title>Identification of mitochondrial DNA sequence variation and development of single nucleotide polymorphic markers for CMS-D8 in cotton</article-title>. <source>Theoretical and Applied Genetics</source><volume>126</volume>: <fpage>1521</fpage>–<lpage>1529</lpage>.<pub-id pub-id-type="pmid">23435734</pub-id></mixed-citation></ref><ref id="pone.0112320-Li2"><label>58</label><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>J</given-names></name>, <name><surname>Pandeya</surname><given-names>D</given-names></name>, <name><surname>Jo</surname><given-names>YD</given-names></name>, <name><surname>Liu</surname><given-names>WY</given-names></name>, <name><surname>Kang</surname><given-names>BC</given-names></name> (<year>2013</year>) <article-title>Reduced activity of ATP synthase in mitochondria causes cytoplasmic male sterility in chili pepper</article-title>. <source>Planta</source><volume>237</volume>: <fpage>1097</fpage>–<lpage>1109</lpage>.<pub-id pub-id-type="pmid">23274393</pub-id></mixed-citation></ref><ref id="pone.0112320-Yuan1"><label>59</label><mixed-citation publication-type="journal"><name><surname>Yuan</surname><given-names>J</given-names></name>, <name><surname>Zhang</surname><given-names>Y</given-names></name>, <name><surname>Liu</surname><given-names>X</given-names></name>, <name><surname>Hao</surname><given-names>R</given-names></name>, <name><surname>Wang</surname><given-names>H</given-names></name> (<year>1996</year>) <article-title>Discovery and study of JA cytoplasmic male sterility</article-title>. <source>China Cotton</source><volume>23</volume>: <fpage>6</fpage>–<lpage>7</lpage>.</mixed-citation></ref><ref id="pone.0112320-Wei2"><label>60</label><mixed-citation publication-type="other">Wei Z, Luo C (1987) Preliminary studies of the cytoplasm influence on economic traits of upland cotton Cotton Breeding Basic Research Papers: 103–107.</mixed-citation></ref><ref id="pone.0112320-Stewart1"><label>61</label><mixed-citation publication-type="other">Stewart JM (1992) A new cytoplasmic male sterile and restorer. Prod Beltw ide Cotton Prod Res Conf.</mixed-citation></ref><ref id="pone.0112320-Huang2"><label>62</label><mixed-citation publication-type="other">Huang J (2003) Genetic study on JinA cytoplasmic male sterile line in cotton. Taigu: Shanxi Agricultural University.</mixed-citation></ref><ref id="pone.0112320-Ma1"><label>63</label><mixed-citation publication-type="journal"><name><surname>Ma</surname><given-names>W</given-names></name>, <name><surname>Ren</surname><given-names>A</given-names></name>, <name><surname>Cui</surname><given-names>M</given-names></name>, <name><surname>Han</surname><given-names>Q</given-names></name>, <name><surname>Zhang</surname><given-names>Y</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>The difference of mitochondrial DNA of three cotton cytoplasmic male sterility lines</article-title>. <source>Cotton Science</source><volume>34</volume>: <fpage>10</fpage>–<lpage>15</lpage>.</mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Bois/explor/OrangerV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0002161 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0002161 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Bois
|area= OrangerV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.25. |  |