Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Genome-Wide Comparative Analysis Reveals Similar Types of NBS Genes in Hybrid <italic>Citrus sinensis</italic> Genome and Original <italic>Citrus clementine</italic> Genome and Provides New Insights into Non-TIR NBS Genes</title><author><name sortKey="Wang, Yunsheng" sort="Wang, Yunsheng" uniqKey="Wang Y" first="Yunsheng" last="Wang">Yunsheng Wang</name><affiliation><nlm:aff id="aff001"><addr-line>College of Plant Protection, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Zhou, Lijuan" sort="Zhou, Lijuan" uniqKey="Zhou L" first="Lijuan" last="Zhou">Lijuan Zhou</name><affiliation><nlm:aff id="aff003"><addr-line>Agricultural Research Service, U.S. Horticultural Research Laboratory, Fort Pierce, Florida, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Li, Dazhi" sort="Li, Dazhi" uniqKey="Li D" first="Dazhi" last="Li">Dazhi Li</name><affiliation><nlm:aff id="aff004"><addr-line>National Center for Citrus Improvement, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Dai, Liangying" sort="Dai, Liangying" uniqKey="Dai L" first="Liangying" last="Dai">Liangying Dai</name><affiliation><nlm:aff id="aff001"><addr-line>College of Plant Protection, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lawton Rauh, Amy" sort="Lawton Rauh, Amy" uniqKey="Lawton Rauh A" first="Amy" last="Lawton-Rauh">Amy Lawton-Rauh</name><affiliation><nlm:aff id="aff005"><addr-line>Department of Genetics and Biochemistry, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Srimani, Pradip K" sort="Srimani, Pradip K" uniqKey="Srimani P" first="Pradip K." last="Srimani">Pradip K. Srimani</name><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Duan, Yongping" sort="Duan, Yongping" uniqKey="Duan Y" first="Yongping" last="Duan">Yongping Duan</name><affiliation><nlm:aff id="aff003"><addr-line>Agricultural Research Service, U.S. Horticultural Research Laboratory, Fort Pierce, Florida, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Luo, Feng" sort="Luo, Feng" uniqKey="Luo F" first="Feng" last="Luo">Feng Luo</name><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">25811466</idno><idno type="pmc">4374887</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4374887</idno><idno type="RBID">PMC:4374887</idno><idno type="doi">10.1371/journal.pone.0121893</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000136</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Genome-Wide Comparative Analysis Reveals Similar Types of NBS Genes in Hybrid <italic>Citrus sinensis</italic> Genome and Original <italic>Citrus clementine</italic> Genome and Provides New Insights into Non-TIR NBS Genes</title><author><name sortKey="Wang, Yunsheng" sort="Wang, Yunsheng" uniqKey="Wang Y" first="Yunsheng" last="Wang">Yunsheng Wang</name><affiliation><nlm:aff id="aff001"><addr-line>College of Plant Protection, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Zhou, Lijuan" sort="Zhou, Lijuan" uniqKey="Zhou L" first="Lijuan" last="Zhou">Lijuan Zhou</name><affiliation><nlm:aff id="aff003"><addr-line>Agricultural Research Service, U.S. Horticultural Research Laboratory, Fort Pierce, Florida, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Li, Dazhi" sort="Li, Dazhi" uniqKey="Li D" first="Dazhi" last="Li">Dazhi Li</name><affiliation><nlm:aff id="aff004"><addr-line>National Center for Citrus Improvement, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Dai, Liangying" sort="Dai, Liangying" uniqKey="Dai L" first="Liangying" last="Dai">Liangying Dai</name><affiliation><nlm:aff id="aff001"><addr-line>College of Plant Protection, Hunan Agricultural University, Changsha, China</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lawton Rauh, Amy" sort="Lawton Rauh, Amy" uniqKey="Lawton Rauh A" first="Amy" last="Lawton-Rauh">Amy Lawton-Rauh</name><affiliation><nlm:aff id="aff005"><addr-line>Department of Genetics and Biochemistry, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Srimani, Pradip K" sort="Srimani, Pradip K" uniqKey="Srimani P" first="Pradip K." last="Srimani">Pradip K. Srimani</name><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Duan, Yongping" sort="Duan, Yongping" uniqKey="Duan Y" first="Yongping" last="Duan">Yongping Duan</name><affiliation><nlm:aff id="aff003"><addr-line>Agricultural Research Service, U.S. Horticultural Research Laboratory, Fort Pierce, Florida, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Luo, Feng" sort="Luo, Feng" uniqKey="Luo F" first="Feng" last="Luo">Feng Luo</name><affiliation><nlm:aff id="aff002"><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>In this study, we identified and compared nucleotide-binding site (NBS) domain-containing genes from three <italic>Citrus</italic> genomes (<italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> from USA and <italic>C</italic>. <italic>sinensis</italic> from China). Phylogenetic analysis of all <italic>Citrus</italic> NBS genes across these three genomes revealed that there are three approximately evenly numbered groups: one group contains the Toll-Interleukin receptor (TIR) domain and two different Non-TIR groups in which most of proteins contain the Coiled Coil (CC) domain. Motif analysis confirmed that the two groups of CC-containing NBS genes are from different evolutionary origins. We partitioned NBS genes into clades using NBS domain sequence distances and found most clades include NBS genes from all three <italic>Citrus</italic> genomes. This suggests that three <italic>Citrus</italic> genomes have similar numbers and types of NBS genes. We also mapped the re-sequenced reads of three pomelo and three mandarin genomes onto the <italic>C</italic>. <italic>sinensis</italic> genome. We found that most NBS genes of the hybrid <italic>C</italic>. <italic>sinensis</italic> genome have corresponding homologous genes in both pomelo and mandarin genomes. The homologous NBS genes in pomelo and mandarin suggest that the parental species of <italic>C</italic>. <italic>sinensis</italic> may contain similar types of NBS genes. This explains why the hybrid <italic>C</italic>. <italic>sinensis</italic> and original <italic>C</italic>. <italic>clementina</italic> have similar types of NBS genes in this study. Furthermore, we found that sequence variation amongst <italic>Citrus</italic> NBS genes were shaped by multiple independent and shared accelerated mutation accumulation events among different groups of NBS genes and in different <italic>Citrus</italic> genomes. Our comparative analyses yield valuable insight into the structure, organization and evolution of NBS genes in <italic>Citrus</italic> genomes. Furthermore, our comprehensive analysis showed that the non-TIR NBS genes can be divided into two groups that come from different evolutionary origins. This provides new insights into non-TIR genes, which have not received much attention.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Dangl, Jl" uniqKey="Dangl J">JL Dangl</name></author><author><name sortKey="Jones, Jd" uniqKey="Jones J">JD Jones</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jia, Y" uniqKey="Jia Y">Y Jia</name></author><author><name sortKey="Mcadams, Sa" uniqKey="Mcadams S">SA McAdams</name></author><author><name sortKey="Bryan, Gt" uniqKey="Bryan G">GT Bryan</name></author><author><name sortKey="Hershey, Hp" uniqKey="Hershey H">HP Hershey</name></author><author><name sortKey="Valent, B" uniqKey="Valent B">B Valent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ellis, Jg" uniqKey="Ellis J">JG Ellis</name></author><author><name sortKey="Lawrence, Gj" uniqKey="Lawrence G">GJ Lawrence</name></author><author><name sortKey="Luck, Je" uniqKey="Luck J">JE Luck</name></author><author><name sortKey="Dodds, Pn" uniqKey="Dodds P">PN Dodds</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ellis, J" uniqKey="Ellis J">J Ellis</name></author><author><name sortKey="Dodds, P" uniqKey="Dodds P">P Dodds</name></author><author><name sortKey="Pryor, T" uniqKey="Pryor T">T Pryor</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mun, J H" uniqKey="Mun J">J-H Mun</name></author><author><name sortKey="Yu, H J" uniqKey="Yu H">H-J Yu</name></author><author><name sortKey="Park, S" uniqKey="Park S">S Park</name></author><author><name sortKey="Park, B S" uniqKey="Park B">B-S Park</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meyers, Bc" uniqKey="Meyers B">BC Meyers</name></author><author><name sortKey="Kozik, A" uniqKey="Kozik A">A Kozik</name></author><author><name sortKey="Griego, A" uniqKey="Griego A">A Griego</name></author><author><name sortKey="Kuang, H" uniqKey="Kuang H">H Kuang</name></author><author><name sortKey="Michelmore, Rw" uniqKey="Michelmore R">RW Michelmore</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wan, H" uniqKey="Wan H">H Wan</name></author><author><name sortKey="Yuan, W" uniqKey="Yuan W">W Yuan</name></author><author><name sortKey="Bo, K" uniqKey="Bo K">K Bo</name></author><author><name sortKey="Shen, J" uniqKey="Shen J">J Shen</name></author><author><name sortKey="Pang, X" uniqKey="Pang X">X Pang</name></author><author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sun, X" uniqKey="Sun X">X Sun</name></author><author><name sortKey="Wang, G L" uniqKey="Wang G">G-L Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, S" uniqKey="Yang S">S Yang</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author><author><name sortKey="Yue, J X" uniqKey="Yue J">J-X Yue</name></author><author><name sortKey="Tian, D" uniqKey="Tian D">D Tian</name></author><author><name sortKey="Chen, J Q" uniqKey="Chen J">J-Q Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ameline Torregrosa, C" uniqKey="Ameline Torregrosa C">C Ameline-Torregrosa</name></author><author><name sortKey="Wang, B B" uniqKey="Wang B">B-B Wang</name></author><author><name sortKey="O Bleness, Ms" uniqKey="O Bleness M">MS O'Bleness</name></author><author><name sortKey="Deshpande, S" uniqKey="Deshpande S">S Deshpande</name></author><author><name sortKey="Zhu, H" uniqKey="Zhu H">H Zhu</name></author><author><name sortKey="Roe, B" uniqKey="Roe B">B Roe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pan, Q" uniqKey="Pan Q">Q Pan</name></author><author><name sortKey="Wendel, J" uniqKey="Wendel J">J Wendel</name></author><author><name sortKey="Fluhr, R" uniqKey="Fluhr R">R Fluhr</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Michelmore, Rw" uniqKey="Michelmore R">RW Michelmore</name></author><author><name sortKey="Meyers, Bc" uniqKey="Meyers B">BC Meyers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guo, Y L" uniqKey="Guo Y">Y-L Guo</name></author><author><name sortKey="Fitz, J" uniqKey="Fitz J">J Fitz</name></author><author><name sortKey="Schneeberger, K" uniqKey="Schneeberger K">K Schneeberger</name></author><author><name sortKey="Ossowski, S" uniqKey="Ossowski S">S Ossowski</name></author><author><name sortKey="Cao, J" uniqKey="Cao J">J Cao</name></author><author><name sortKey="Weigel, D" uniqKey="Weigel D">D Weigel</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Nicolosi, E" uniqKey="Nicolosi E">E Nicolosi</name></author><author><name sortKey="Deng, Z" uniqKey="Deng Z">Z Deng</name></author><author><name sortKey="Gentile, A" uniqKey="Gentile A">A Gentile</name></author><author><name sortKey="La Malfa, S" uniqKey="La Malfa S">S La Malfa</name></author><author><name sortKey="Continella, G" uniqKey="Continella G">G Continella</name></author><author><name sortKey="Tribulato, E" uniqKey="Tribulato E">E Tribulato</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, X" uniqKey="Li X">X Li</name></author><author><name sortKey="Xie, R" uniqKey="Xie R">R Xie</name></author><author><name sortKey="Lu, Z" uniqKey="Lu Z">Z Lu</name></author><author><name sortKey="Zhou, Z" uniqKey="Zhou Z">Z Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, Q" uniqKey="Xu Q">Q Xu</name></author><author><name sortKey="Chen, Ll" uniqKey="Chen L">LL Chen</name></author><author><name sortKey="Ruan, X" uniqKey="Ruan X">X Ruan</name></author><author><name sortKey="Chen, D" uniqKey="Chen D">D Chen</name></author><author><name sortKey="Zhu, A" uniqKey="Zhu A">A Zhu</name></author><author><name sortKey="Chen, C" uniqKey="Chen C">C Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, Q" uniqKey="Xu Q">Q Xu</name></author><author><name sortKey="Biswas, Mk" uniqKey="Biswas M">MK Biswas</name></author><author><name sortKey="Lan, H" uniqKey="Lan H">H Lan</name></author><author><name sortKey="Zeng, W" uniqKey="Zeng W">W Zeng</name></author><author><name sortKey="Liu, C" uniqKey="Liu C">C Liu</name></author><author><name sortKey="Xu, J" uniqKey="Xu J">J Xu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Ga" uniqKey="Wu G">GA Wu</name></author><author><name sortKey="Prochnik, S" uniqKey="Prochnik S">S Prochnik</name></author><author><name sortKey="Jenkins, J" uniqKey="Jenkins J">J Jenkins</name></author><author><name sortKey="Salse, J" uniqKey="Salse J">J Salse</name></author><author><name sortKey="Hellsten, U" uniqKey="Hellsten U">U Hellsten</name></author><author><name sortKey="Murat, F" uniqKey="Murat F">F Murat</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Finn, Rd" uniqKey="Finn R">RD Finn</name></author><author><name sortKey="Clements, J" uniqKey="Clements J">J Clements</name></author><author><name sortKey="Eddy, Sr" uniqKey="Eddy S">SR Eddy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Punta, M" uniqKey="Punta M">M Punta</name></author><author><name sortKey="Coggill, Pc" uniqKey="Coggill P">PC Coggill</name></author><author><name sortKey="Eberhardt, Ry" uniqKey="Eberhardt R">RY Eberhardt</name></author><author><name sortKey="Mistry, J" uniqKey="Mistry J">J Mistry</name></author><author><name sortKey="Tate, J" uniqKey="Tate J">J Tate</name></author><author><name sortKey="Boursnell, C" uniqKey="Boursnell C">C Boursnell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boeckmann, B" uniqKey="Boeckmann B">B Boeckmann</name></author><author><name sortKey="Bairoch, A" uniqKey="Bairoch A">A Bairoch</name></author><author><name sortKey="Apweiler, R" uniqKey="Apweiler R">R Apweiler</name></author><author><name sortKey="Blatter, Mc" uniqKey="Blatter M">MC Blatter</name></author><author><name sortKey="Estreicher, A" uniqKey="Estreicher A">A Estreicher</name></author><author><name sortKey="Gasteiger, E" uniqKey="Gasteiger E">E Gasteiger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name></author><author><name sortKey="Madden, Tl" uniqKey="Madden T">TL Madden</name></author><author><name sortKey="Schaffer, Aa" uniqKey="Schaffer A">AA Schaffer</name></author><author><name sortKey="Zhang, J" uniqKey="Zhang J">J Zhang</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Birney, E" uniqKey="Birney E">E Birney</name></author><author><name sortKey="Clamp, M" uniqKey="Clamp M">M Clamp</name></author><author><name sortKey="Durbin, R" uniqKey="Durbin R">R Durbin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Moreno Hagelsieb, G" uniqKey="Moreno Hagelsieb G">G Moreno-Hagelsieb</name></author><author><name sortKey="Latimer, K" uniqKey="Latimer K">K Latimer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Suyama, M" uniqKey="Suyama M">M Suyama</name></author><author><name sortKey="Torrents, D" uniqKey="Torrents D">D Torrents</name></author><author><name sortKey="Bork, P" uniqKey="Bork P">P Bork</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Katoh, K" uniqKey="Katoh K">K Katoh</name></author><author><name sortKey="Standley, Dm" uniqKey="Standley D">DM Standley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Price, Mn" uniqKey="Price M">MN Price</name></author><author><name sortKey="Dehal, Ps" uniqKey="Dehal P">PS Dehal</name></author><author><name sortKey="Arkin, Ap" uniqKey="Arkin A">AP Arkin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Retief, Jd" uniqKey="Retief J">JD Retief</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Prosperi, Mc" uniqKey="Prosperi M">MC Prosperi</name></author><author><name sortKey="Ciccozzi, M" uniqKey="Ciccozzi M">M Ciccozzi</name></author><author><name sortKey="Fanti, I" uniqKey="Fanti I">I Fanti</name></author><author><name sortKey="Saladini, F" uniqKey="Saladini F">F Saladini</name></author><author><name sortKey="Pecorari, M" uniqKey="Pecorari M">M Pecorari</name></author><author><name sortKey="Borghi, V" uniqKey="Borghi V">V Borghi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Delorenzi, M" uniqKey="Delorenzi M">M Delorenzi</name></author><author><name sortKey="Speed, T" uniqKey="Speed T">T Speed</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lupas, A" uniqKey="Lupas A">A Lupas</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bailey, Tl" uniqKey="Bailey T">TL Bailey</name></author><author><name sortKey="Boden, M" uniqKey="Boden M">M Boden</name></author><author><name sortKey="Buske, Fa" uniqKey="Buske F">FA Buske</name></author><author><name sortKey="Frith, M" uniqKey="Frith M">M Frith</name></author><author><name sortKey="Grant, Ce" uniqKey="Grant C">CE Grant</name></author><author><name sortKey="Clementi, L" uniqKey="Clementi L">L Clementi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Carriero, N" uniqKey="Carriero N">N Carriero</name></author><author><name sortKey="Zheng, D" uniqKey="Zheng D">D Zheng</name></author><author><name sortKey="Karro, J" uniqKey="Karro J">J Karro</name></author><author><name sortKey="Harrison, Pm" uniqKey="Harrison P">PM Harrison</name></author><author><name sortKey="Gerstein, M" uniqKey="Gerstein M">M Gerstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, Z" uniqKey="Xu Z">Z Xu</name></author><author><name sortKey="Wang, H" uniqKey="Wang H">H Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Y" uniqKey="Wang Y">Y Wang</name></author><author><name sortKey="Tang, H" uniqKey="Tang H">H Tang</name></author><author><name sortKey="Debarry, Jd" uniqKey="Debarry J">JD Debarry</name></author><author><name sortKey="Tan, X" uniqKey="Tan X">X Tan</name></author><author><name sortKey="Li, J" uniqKey="Li J">J Li</name></author><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sawyer, S" uniqKey="Sawyer S">S Sawyer</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Quinlan, Ar" uniqKey="Quinlan A">AR Quinlan</name></author><author><name sortKey="Hall, Im" uniqKey="Hall I">IM Hall</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, L" uniqKey="Zhou L">L Zhou</name></author><author><name sortKey="Powell, Ca" uniqKey="Powell C">CA Powell</name></author><author><name sortKey="Hoffman, Mt" uniqKey="Hoffman M">MT Hoffman</name></author><author><name sortKey="Li, W" uniqKey="Li W">W Li</name></author><author><name sortKey="Fan, G" uniqKey="Fan G">G Fan</name></author><author><name sortKey="Liu, B" uniqKey="Liu B">B Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kohler, A" uniqKey="Kohler A">A Kohler</name></author><author><name sortKey="Rinaldi, C" uniqKey="Rinaldi C">C Rinaldi</name></author><author><name sortKey="Duplessis, S" uniqKey="Duplessis S">S Duplessis</name></author><author><name sortKey="Baucher, M" uniqKey="Baucher M">M Baucher</name></author><author><name sortKey="Geelen, D" uniqKey="Geelen D">D Geelen</name></author><author><name sortKey="Duchaussoy, F" uniqKey="Duchaussoy F">F Duchaussoy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, T" uniqKey="Zhou T">T Zhou</name></author><author><name sortKey="Wang, Y" uniqKey="Wang Y">Y Wang</name></author><author><name sortKey="Chen, Jq" uniqKey="Chen J">JQ Chen</name></author><author><name sortKey="Araki, H" uniqKey="Araki H">H Araki</name></author><author><name sortKey="Jing, Z" uniqKey="Jing Z">Z Jing</name></author><author><name sortKey="Jiang, K" uniqKey="Jiang K">K Jiang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tan, S" uniqKey="Tan S">S Tan</name></author><author><name sortKey="Wu, S" uniqKey="Wu S">S Wu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Punta, M" uniqKey="Punta M">M Punta</name></author><author><name sortKey="Coggill, Pc" uniqKey="Coggill P">PC Coggill</name></author><author><name sortKey="Eberhardt, Ry" uniqKey="Eberhardt R">RY Eberhardt</name></author><author><name sortKey="Mistry, J" uniqKey="Mistry J">J Mistry</name></author><author><name sortKey="Tate, J" uniqKey="Tate J">J Tate</name></author><author><name sortKey="Boursnell, C" uniqKey="Boursnell C">C Boursnell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, S" uniqKey="Yang S">S Yang</name></author><author><name sortKey="Feng, Z" uniqKey="Feng Z">Z Feng</name></author><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name></author><author><name sortKey="Jiang, K" uniqKey="Jiang K">K Jiang</name></author><author><name sortKey="Jin, X" uniqKey="Jin X">X Jin</name></author><author><name sortKey="Hang, Y" uniqKey="Hang Y">Y Hang</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, CA USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">25811466</article-id><article-id pub-id-type="pmc">4374887</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0121893</article-id><article-id pub-id-type="publisher-id">PONE-D-14-48519</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Genome-Wide Comparative Analysis Reveals Similar Types of NBS Genes in Hybrid <italic>Citrus sinensis</italic> Genome and Original <italic>Citrus clementine</italic> Genome and Provides New Insights into Non-TIR NBS Genes</article-title><alt-title alt-title-type="running-head">Comparative Analysis of Citrus NBS Genes</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Wang</surname><given-names>Yunsheng</given-names></name><xref ref-type="aff" rid="aff001"><sup>1</sup></xref><xref ref-type="aff" rid="aff002"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Zhou</surname><given-names>Lijuan</given-names></name><xref ref-type="aff" rid="aff003"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Li</surname><given-names>Dazhi</given-names></name><xref ref-type="aff" rid="aff004"><sup>4</sup></xref></contrib><contrib contrib-type="author"><name><surname>Dai</surname><given-names>Liangying</given-names></name><xref ref-type="aff" rid="aff001"><sup>1</sup></xref><xref rid="cor001" ref-type="corresp">*</xref></contrib><contrib contrib-type="author"><name><surname>Lawton-Rauh</surname><given-names>Amy</given-names></name><xref ref-type="aff" rid="aff005"><sup>5</sup></xref></contrib><contrib contrib-type="author"><name><surname>Srimani</surname><given-names>Pradip K.</given-names></name><xref ref-type="aff" rid="aff002"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Duan</surname><given-names>Yongping</given-names></name><xref ref-type="aff" rid="aff003"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Luo</surname><given-names>Feng</given-names></name><xref ref-type="aff" rid="aff002"><sup>2</sup></xref><xref rid="cor001" ref-type="corresp">*</xref></contrib></contrib-group><aff id="aff001"><label>1</label><addr-line>College of Plant Protection, Hunan Agricultural University, Changsha, China</addr-line></aff><aff id="aff002"><label>2</label><addr-line>School of Computing, Clemson University, Clemson, South Carolina, United States of America</addr-line></aff><aff id="aff003"><label>3</label><addr-line>Agricultural Research Service, U.S. Horticultural Research Laboratory, Fort Pierce, Florida, United States of America</addr-line></aff><aff id="aff004"><label>4</label><addr-line>National Center for Citrus Improvement, Hunan Agricultural University, Changsha, China</addr-line></aff><aff id="aff005"><label>5</label><addr-line>Department of Genetics and Biochemistry, Clemson University, Clemson, South Carolina, United States of America</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Liu</surname><given-names>Ji-Hong</given-names></name><role>Academic Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>Key Laboratory of Horticultural Plant Biology (MOE), CHINA</addr-line></aff><author-notes><fn fn-type="conflict" id="coi001"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con" id="contrib001"><p>Conceived and designed the experiments: FL YD PS LD. Performed the experiments: LZ. Analyzed the data: YW DL ALR FL. Wrote the paper: ALR PS LD YD FL.</p></fn><corresp id="cor001">* E-mail: <email>daily@hunau.net</email> (LD); <email>luofeng@clemson.edu</email> (FL)</corresp></author-notes><pub-date pub-type="epub"><day>26</day><month>3</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>10</volume><issue>3</issue><elocation-id>e0121893</elocation-id><history><date date-type="received"><day>28</day><month>10</month><year>2014</year></date><date date-type="accepted"><day>4</day><month>2</month><year>2015</year></date></history><permissions><license xlink:href="https://creativecommons.org/publicdomain/zero/1.0/"><license-p>This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/publicdomain/zero/1.0/">Creative Commons CC0</ext-link> public domain dedication</license-p></license></permissions><self-uri content-type="pdf" xlink:type="simple" xlink:href="pone.0121893.pdf"></self-uri><abstract><p>In this study, we identified and compared nucleotide-binding site (NBS) domain-containing genes from three <italic>Citrus</italic> genomes (<italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> from USA and <italic>C</italic>. <italic>sinensis</italic> from China). Phylogenetic analysis of all <italic>Citrus</italic> NBS genes across these three genomes revealed that there are three approximately evenly numbered groups: one group contains the Toll-Interleukin receptor (TIR) domain and two different Non-TIR groups in which most of proteins contain the Coiled Coil (CC) domain. Motif analysis confirmed that the two groups of CC-containing NBS genes are from different evolutionary origins. We partitioned NBS genes into clades using NBS domain sequence distances and found most clades include NBS genes from all three <italic>Citrus</italic> genomes. This suggests that three <italic>Citrus</italic> genomes have similar numbers and types of NBS genes. We also mapped the re-sequenced reads of three pomelo and three mandarin genomes onto the <italic>C</italic>. <italic>sinensis</italic> genome. We found that most NBS genes of the hybrid <italic>C</italic>. <italic>sinensis</italic> genome have corresponding homologous genes in both pomelo and mandarin genomes. The homologous NBS genes in pomelo and mandarin suggest that the parental species of <italic>C</italic>. <italic>sinensis</italic> may contain similar types of NBS genes. This explains why the hybrid <italic>C</italic>. <italic>sinensis</italic> and original <italic>C</italic>. <italic>clementina</italic> have similar types of NBS genes in this study. Furthermore, we found that sequence variation amongst <italic>Citrus</italic> NBS genes were shaped by multiple independent and shared accelerated mutation accumulation events among different groups of NBS genes and in different <italic>Citrus</italic> genomes. Our comparative analyses yield valuable insight into the structure, organization and evolution of NBS genes in <italic>Citrus</italic> genomes. Furthermore, our comprehensive analysis showed that the non-TIR NBS genes can be divided into two groups that come from different evolutionary origins. This provides new insights into non-TIR genes, which have not received much attention.</p></abstract><funding-group><funding-statement>This work was supported in part by the Citrus Research and Development Foundation, Inc. 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><fig-count count="6"></fig-count><table-count count="4"></table-count><page-count count="24"></page-count></counts><custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name><meta-value>All relevant data are within the paper and its Supporting Information files.</meta-value></custom-meta></custom-meta-group></article-meta><notes><title>Data Availability</title><p>All relevant data are within the paper and its Supporting Information files.</p></notes></front><body><sec sec-type="intro" id="sec001"><title>Introduction</title><p>Disease resistance genes (<italic>R</italic> genes) are essential components of plant immune systems. Amongst five diverse classes of disease resistance genes [<xref rid="pone.0121893.ref001" ref-type="bibr">1</xref>], the largest class of <italic>R</italic> genes includes genes that encode proteins with a nucleotide-binding site (NBS) domain and a Leucine-Rich Repeat (LRR) domain. The NBS domain is highly sequence constrained and is typically used to identify and characterize plant <italic>R</italic> genes. The LRR domains in <italic>R</italic> genes that mediate direct or indirect interactions [<xref rid="pone.0121893.ref002" ref-type="bibr">2</xref>, <xref rid="pone.0121893.ref003" ref-type="bibr">3</xref>] with pathogen molecules are usually rapidly evolving to adapt to the change of pathogen ligands [<xref rid="pone.0121893.ref004" ref-type="bibr">4</xref>]. NBS domain-containing <italic>R</italic> genes were classified into two major types based on their domain structures in N terminus: proteins with a Toll-Interleukin receptor (TIR) domain and proteins without TIR domains, which usually contain a Coiled Coil (CC) domain. While the TIR domain is usually well defined, the CC domain has higher sequence variation and is less well characterized.</p><p>Genome-wide analyses of NBS genes in many genomes have shown that NBS <italic>R</italic> genes are diverse in number, structure and organization [<xref rid="pone.0121893.ref005" ref-type="bibr">5</xref>–<xref rid="pone.0121893.ref010" ref-type="bibr">10</xref>]. For example, grass genomes studied to date do not have TIR NBS genes; on the other hand, dicot genomes usually contain more TIR NBS genes than non-TIR NBS genes [<xref rid="pone.0121893.ref011" ref-type="bibr">11</xref>]. Comparison of NBS <italic>R</italic> genes in multiple genomes indicated that <italic>R</italic> genes are shaped by dynamic birth-and-death processes [<xref rid="pone.0121893.ref012" ref-type="bibr">12</xref>]. Comparative analyses of NBS genes in two <italic>Arabidopsis</italic> genomes (<italic>A</italic>. <italic>thaliana</italic> and <italic>A</italic>. <italic>lyrata</italic>) indicate that mating system shift from outcrossing to inbreeding has had a limited impact on the numbers of NBS genes, at least in the time since <italic>A</italic>. <italic>lyrata</italic> and <italic>A</italic>. <italic>thaliana</italic> diverged from their common ancestor approximately 5 million years ago [<xref rid="pone.0121893.ref013" ref-type="bibr">13</xref>]. Comparative analyses of NBS genes of diploid <italic>Phaseolus</italic> and tetraploid <italic>Glycine</italic> species concluded that whole genome duplication did not result in NBS <italic>R</italic> gene number increasing. Recently, comparison of <italic>R</italic> genes in diverse grass genomes by Yang <italic>et al</italic>. showed that rapid evolution in <italic>R</italic> genes in <italic>Zea</italic>, <italic>Sorghum</italic>, and <italic>Brachypodium</italic> is associated with rice blast disease resistance [<xref rid="pone.0121893.ref014" ref-type="bibr">14</xref>].</p><p><italic>Citrus</italic> species are amongst the most important fruit trees and have been cultivated for more than 4000 years [<xref rid="pone.0121893.ref015" ref-type="bibr">15</xref>, <xref rid="pone.0121893.ref016" ref-type="bibr">16</xref>]. Phylogenetic analyses using molecular markers showed that cultivated <italic>Citrus</italic> species (sweet orange, grapefruit, and lemon) are derived from three original cultivated <italic>Citrus</italic> species: <italic>C</italic>. <italic>medica</italic> (citron), <italic>C</italic>. <italic>reticulata</italic> (mandarin) and <italic>C</italic>. <italic>maxima</italic> (pomelo) [<xref rid="pone.0121893.ref017" ref-type="bibr">17</xref>, <xref rid="pone.0121893.ref018" ref-type="bibr">18</xref>]. For example, the <italic>C</italic>. <italic>sinensis</italic> (sweet orange) is suggested to be the backcross hybrid of <italic>C</italic>. <italic>maxima</italic> (pomelo) and <italic>C</italic>. <italic>reticulata</italic> (mandarin) [<xref rid="pone.0121893.ref017" ref-type="bibr">17</xref>–<xref rid="pone.0121893.ref019" ref-type="bibr">19</xref>]. Phylogenetic analyses of partial NBS genes of <italic>Poncirus trifoliata</italic> (trifoliate orange), <italic>C</italic>. <italic>reticulata</italic> (tangerine) and their F1 progeny showed that NBS genes of <italic>Poncirus trifoliata</italic> (trifoliate orange), <italic>C</italic>. <italic>reticulata</italic> formed genus-specific clades. Additionally, NBS genes of their F1 progeny had sister relationships to only one of the parents [<xref rid="pone.0121893.ref020" ref-type="bibr">20</xref>]. This suggests that NBS genes in crossing hybrid <italic>Citrus</italic> species are also different from those in original <italic>Citrus</italic> species.</p><p>Recently available draft whole genome sequences of three <italic>Citrus</italic> species: <italic>C</italic>. <italic>clementina</italic> [<xref rid="pone.0121893.ref021" ref-type="bibr">21</xref>], <italic>C</italic>. <italic>sinensis</italic> from USA and <italic>C</italic>. <italic>sinensis</italic> from China [<xref rid="pone.0121893.ref019" ref-type="bibr">19</xref>], have made it possible to scan and identify all NBS genes in those genomes. In this study, we performed a genome-wide comparative analysis of NBS genes in three <italic>Citrus</italic> genomes (<italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> from China and <italic>C</italic>. <italic>sinensis</italic> from USA) to address the following questions: (1) What are the features of NBS genes in <italic>Citrus</italic>, such as numbers, physical locations, within-gene domain structures, and evolutionary dynamics of NBS genes?; (2) Do NBS genes differ between hybrid and non-hybrid <italic>Citrus</italic> species?; and (3) Do <italic>Citrus</italic> NBS genes differ from NBS genes in other plant genomes?</p></sec><sec sec-type="materials|methods" id="sec002"><title>Materials and Methods</title><sec id="sec003"><title>Sequences Used</title><p>We downloaded the draft genome sequences and the original gene annotations of <italic>C</italic>. <italic>clementina</italic> (clementine) [<xref rid="pone.0121893.ref021" ref-type="bibr">21</xref>] and <italic>C</italic>. <italic>sinensis</italic> USA (sweet orange from USA) from the <italic>Citrus</italic> Genome Database (<ext-link ext-link-type="uri" xlink:href="http://www.citrusgenomedb.org/">http://www.citrusgenomedb.org/</ext-link>) and those of <italic>C</italic>. <italic>sinensis</italic> China (Chinese sweet orange [<xref rid="pone.0121893.ref019" ref-type="bibr">19</xref>]) from the <italic>C</italic>. <italic>sinensis</italic> annotation project (<ext-link ext-link-type="uri" xlink:href="http://citrus.hzau.edu.cn/orange/">http://citrus.hzau.edu.cn/orange/</ext-link>). The sizes of assembled genomes of <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA are 301, 328 and 319 million base pairs, respectively. The <italic>C</italic>. <italic>clementina</italic> genomes assembled into nine major scaffolds and 95.8% of the sequences were assigned to those nine scaffolds. About 72.9% of the <italic>C</italic>. <italic>sinensis</italic> China genome assigned to nine chromosomes. The genome of <italic>C</italic>. <italic>sinensis</italic> USA only assembled to 12,574 scaffolds and the N50 of scaffolds is 250 kb. The original gene annotations of <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA have 24,533, 29,385 and 25,397 genes respectively. We also downloaded the re-sequence data of three mandarin and three pomelo genomes from the <italic>C</italic>. <italic>sinensis</italic> annotation project (<ext-link ext-link-type="uri" xlink:href="http://citrus.hzau.edu.cn/orange/download/data.php">http://citrus.hzau.edu.cn/orange/download/data.php</ext-link>). The re-sequence data is single end Illumina reads with length of 101 bp.</p></sec><sec id="sec004"><title>Identification of NBS Genes in <italic>Citrus</italic> Genomes</title><p>We first screened the original predicted <italic>Citrus</italic> open reading frames (ORFs) using hmmsearch [<xref rid="pone.0121893.ref022" ref-type="bibr">22</xref>] with the hidden Markov models (HMM) of Pfam [<xref rid="pone.0121893.ref023" ref-type="bibr">23</xref>] for NBS domain presence (NB-ARC, PF00931) using an e-value cut-off of 0.1 for the <bold>hmmsearch</bold>. Then, the proteins selected by the HMM were searched against the Swiss protein database [<xref rid="pone.0121893.ref024" ref-type="bibr">24</xref>] using BLASTP [<xref rid="pone.0121893.ref025" ref-type="bibr">25</xref>]. Only the proteins that have a significant match (e-value < 1E-5 in <bold>BLASTP</bold> search) with the NBS proteins or resistant proteins in the Swiss protein database were identified as potential NBS-containing proteins. To recover possible NBS genes that may be missed in the original gene annotations, we mapped the identified NBS genes to the draft genome using <bold>TBLASTN</bold>. The matched sequences with e-value < 1E-5 were then predicted using <bold>Genewise</bold> [<xref rid="pone.0121893.ref026" ref-type="bibr">26</xref>]. The new potential NBS genes predicted by Genewise were also confirmed by <bold>BLASTP</bold> searching through the Swiss protein database. Finally, we scanned all potential NBS genes for NBS domain using <bold>hmmsearch</bold> [<xref rid="pone.0121893.ref022" ref-type="bibr">22</xref>] and only keep the NBS genes with e value lower than more stringent cutoff of 1E-5.</p></sec><sec id="sec005"><title>Identification of Orthologous NBS genes in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic></title><p>Orthologous NBS genes of <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> were identified using the reciprocal best blast method [<xref rid="pone.0121893.ref027" ref-type="bibr">27</xref>]. We used NBS genes of <italic>C</italic>. <italic>clementina</italic> as query sequences to search against the NBS genes of <italic>C</italic>. <italic>sinensis</italic> and <italic>vice versa</italic>. Protein pairs with reciprocal best hits of e-value < 1E-20 were defined as orthologs.</p><p>To calculate the rates of nonsynonymous, synonymous and their rate ratio (dN, dS and dN/dS) of orthologous pairs, we first aligned orthologous protein sequence pairs using <bold>mafft</bold>, and then converted the protein alignments to codon-based alignments using <bold>PAL2NAL</bold> [<xref rid="pone.0121893.ref028" ref-type="bibr">28</xref>]. We calculated the dN, dS and dN/dS rate ratios using the <bold>codeml</bold> program in <bold>PAML</bold> version 4.7 [<xref rid="pone.0121893.ref029" ref-type="bibr">29</xref>].</p></sec><sec id="sec006"><title>Phylogenetic Analysis</title><p>We constructed the best fit phylogenetic tree of NBS genes of the three <italic>Citrus</italic> genomes using only the conserved NBS domain. Only 1,099 NBS domain sequences that have both the P-loop and the MHDV motifs and were longer than 70% (200 amino acids) of the full-length NBS domain were included. Next, NBS domain sequences were aligned in <bold>mafft</bold> [<xref rid="pone.0121893.ref030" ref-type="bibr">30</xref>] using an auto alignment model and a best fit maximum likelihood phylogenetic tree of NBS genes was constructed using <bold>FastTree</bold> [<xref rid="pone.0121893.ref031" ref-type="bibr">31</xref>] with default parameters (JTT+CAT). The resultant best fit phylogenetic tree was divided into 3 main groups based upon clade support. For each group, we constructed a new ML tree with <italic>Streptomyces coelicolor</italic> protein P25941 as the outgroup using <bold>FastTree</bold> [<xref rid="pone.0121893.ref031" ref-type="bibr">31</xref>]. The average identity of each group was calculated using alistat implemented in <bold>SQUID</bold> (<ext-link ext-link-type="uri" xlink:href="http://selab.janelia.org/software.html">http://selab.janelia.org/software.html</ext-link>).</p><p>To evaluate the confidence level of the phylogenetic tree, we resampled the alignments 1000 times using <bold>SEQBOOT</bold> [<xref rid="pone.0121893.ref032" ref-type="bibr">32</xref>]. Then, we constructed the phylogenetic tree from each resampled alignments using <bold>FastTree</bold> [<xref rid="pone.0121893.ref031" ref-type="bibr">31</xref>]. The bootstrap values were calculated by comparing resampled phylogenetic trees to the original tree using <bold>CompareToBootstrap</bold> program in the <bold>FastTree</bold> [<xref rid="pone.0121893.ref031" ref-type="bibr">31</xref>].</p><p>We further partitioned the NBS gene tree into clades using the depth-first phylogeny partition method in <bold>PhyloPart</bold> [<xref rid="pone.0121893.ref033" ref-type="bibr">33</xref>] with distance threshold 0.025. This clustered the NBS genes into 114 clades, which contain more than 2 genes, and 8 orphan genes.</p></sec><sec id="sec007"><title>Domain and Motif Annotation</title><p>The Toll-Interleukin receptor (TIR) domains in <italic>Citrus</italic> NBS-containing proteins were identified using <bold>hmmsearch</bold> [<xref rid="pone.0121893.ref022" ref-type="bibr">22</xref>] with the HMM model of Pfam domain PF01582 and a 0.1 e-value cut-off. The Leucine-Rich Repeat (LRR) domains in <italic>Citrus</italic> NBS containing proteins were identified using the HMM models of Pfam LRR domains with e-value cut-off of 0.1. As long as there is a significant hit to one of LRR domain models (e-value < 0.1), it was define as an LRR-containing protein. We used <bold>MARCOIL</bold> [<xref rid="pone.0121893.ref034" ref-type="bibr">34</xref>] with a threshold probability of 90 and <bold>COILS</bold> [<xref rid="pone.0121893.ref035" ref-type="bibr">35</xref>] with a threshold of 0.9 to search for Coiled Coil (CC) domains in the N-terminal region of <italic>Citrus</italic> NBS-containing protein. We considered a protein as a CC-containing protein if either <bold>MARCOI</bold>L or <bold>COILS</bold> reported a CC domain in it.</p><p>We identified 20 motifs amongst the NBS genes in each of the three main phylogenetic groups separately using <bold>MEME SUITE</bold> [<xref rid="pone.0121893.ref036" ref-type="bibr">36</xref>]. The motif width was set to between 6 and 50 for <bold>MEME</bold>. Then, we searched the motif structure of all genes in each group using <bold>MAST</bold> with default parameters (-ev 10—mt 0.0001).</p></sec><sec id="sec008"><title>Pseudogene Identification</title><p>We identified possible pseudogenes using <bold>PseudoPipe</bold> [<xref rid="pone.0121893.ref037" ref-type="bibr">37</xref>] with default parameters (-e 0.1). The <bold>PseudoPipe</bold> algorithm identifies pseudo genes by integrating sequence similarity, intron-exon structure, plus presence of stop codons and frame-shifts. We used all <italic>Citrus</italic> NBS genes to search <italic>Citrus</italic> genomes for potential NGS pseudo genes.</p></sec><sec id="sec009"><title>Transposon Identification</title><p>All long terminal repeat (LTR) retrotransponsons in each <italic>Citrus</italic> genome were identified using LTR finder [<xref rid="pone.0121893.ref038" ref-type="bibr">38</xref>] with default parameters (-o 3-t 1—e 1-m 2-u -2). Then, we used a script program to match the location of the LTR transposons to the NBS-LRR genes in each <italic>Citrus</italic> genome.</p></sec><sec id="sec010"><title>NBS Gene Synteny Identification</title><p>Gene synteny was identified and defined using <bold>MCScanX</bold> [<xref rid="pone.0121893.ref039" ref-type="bibr">39</xref>] with default parameters (-A-u 5000). First, NBS genes between two genomes were aligned using <bold>BLASTP</bold> and matches with E-value < 1e-5 were sorted according to their chromosome positions. Synteny scores were then calculated for each block based upon gene position. Two genes were considered in the same block if there were fewer than 25 genes separating them. <bold>MCScanX</bold> reported blocks with at least 5 collinear gene pairs.</p></sec><sec id="sec011"><title>Gene Cluster Analysis</title><p>We grouped the NBS genes in each <italic>Citrus</italic> genome into the same cluster if the genome location between two genes was within 200 kb. We also identified the conserved gene clusters between the <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic>. If all genes in a cluster of <italic>C</italic>. <italic>clementina</italic> have orthologs in the corresponding cluster of <italic>C</italic>. <italic>sinensis</italic> and vice versa, then these two clusters were called completely conserved clusters. If only part of genes in the clusters have orthologs, we called these two clusters partially conserved.</p></sec><sec id="sec012"><title>Gene Conversion Detection</title><p>We first aligned the sequences of NBS genes in the same cluster using <bold>mafft</bold> [<xref rid="pone.0121893.ref030" ref-type="bibr">30</xref>]. Then, we used <bold>GENECONV</bold> [<xref rid="pone.0121893.ref040" ref-type="bibr">40</xref>] version 1.81a with default settings (N = 10,000) to detect gene conversions. <bold>GENECONV</bold> identifies gene conversions by finding identical fragments between pairs of sequences in a nucleotide alignment. A global <italic>P</italic> value ≤ 0.05 was used to assess the statistical significance of the observed conversions. <bold>GENECONV</bold> requires at least three sequences for analyses in order to account for shared ancestral states. Thus, we only detected conversions in clusters containing three or more genes.</p></sec><sec id="sec013"><title>Tests for Sites under Positive Selection</title><p>The amino acid sequences of NBS genes from the same clade were aligned with <bold>mafft</bold> [<xref rid="pone.0121893.ref030" ref-type="bibr">30</xref>]. Then we converted the protein alignments to codon-based alignment using <bold>PAL2NAL</bold> [<xref rid="pone.0121893.ref028" ref-type="bibr">28</xref>]. The positively selected sites were statistically identified using the Bayesian approach implemented in <bold>codeml</bold> within <bold>PAML</bold> [<xref rid="pone.0121893.ref029" ref-type="bibr">29</xref>]. We also further examined sites in the ω > 1 class with >90% posterior probability.</p></sec><sec id="sec014"><title>Mapping of Re-sequencing Data of Mandarin and Pomelo</title><p>We mapped the raw reads of each re-sequenced sample to the draft genome of <italic>C</italic>. <italic>sinensis</italic> China using <bold>BWA</bold> [<xref rid="pone.0121893.ref041" ref-type="bibr">41</xref>] with default parameters (-k 19-d 100-A 1-B 4-O 6). Then, the mapped reads of NBS genes regions were extracted using <bold>BEDtools</bold> [<xref rid="pone.0121893.ref042" ref-type="bibr">42</xref>]. Two types of coverage of each NBS gene in <italic>C</italic>. <italic>sinensis</italic> were calculated. One coverage type divides the length of mapped sequences by the whole gene sequence length and the other type divides the length of mapped sequences by the length of exon sequence only.</p></sec><sec id="sec015"><title><italic>Citrus</italic> DNA Extraction and PCR Amplication</title><p><italic>Citrus</italic> leaf samples were collected from the six <italic>Citrus</italic> plants in USHRL's (USDA Horticultural Research Laboratory, Fort Pierce, Florida): <italic>C</italic>. <italic>sinensis</italic> (sweet, navel orange), <italic>C</italic>. <italic>aurantium</italic> (karum jamir, sour orange), <italic>C</italic>. <italic>reticulata</italic> (mandarin orange), <italic>C</italic>. <italic>clementina</italic> (clementina), <italic>C</italic>. <italic>aurantiifolia</italic> (sweet lime), <italic>C</italic>. <italic>japonica</italic> (Yuzu, kumquat), and <italic>C</italic>. <italic>maxima</italic> (pomelo). Total DNA was extracted from leaf midribs following the Plant Mini Kit standard protocol from Qiagen Inc. (Valencia, CA), followed by DNA quantity and quality evaluation with Nanodrop. We chose the NBS gene, Cs1g09350, which was conserved in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> for validating the conservation of NBS gene among different <italic>Citrus</italic> genomes. Primers used in this study were designed using Oligo 7.23 (Molecular Biology Insights, Inc., Cascade, CO, USA). DNA Polymerase (Invitrogen, Carlsbad, CA, USA) was used to amplify the NBS-LRR genes from <italic>Citrus</italic> DNA. For PCR, 20 μL reactions using standard conditions provided by the manufacturer for DNA Polymerase. PCR was performed using an initial denaturation at 95°C for 3 minutes, 35 cycles of 94°C for 20 seconds, 50–52°C for 20 seconds (specified by different primer sets) and 68°C for 3 minutes, follow by final extension at 68°C for 10 minutes. The cloning and sequencing analysis of amplified PCR products were conducted as previously described [<xref rid="pone.0121893.ref043" ref-type="bibr">43</xref>]</p></sec></sec><sec sec-type="results" id="sec016"><title>Results</title><sec id="sec017"><title>Identification and classification of <italic>Citrus</italic> NBS Genes</title><p>We searched the <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA genomes for genes containing the NBS domain using <bold>hmmsearch</bold>. Then, the NBS-containing genes were confirmed through homology searches against the Swiss protein database (see “<xref rid="sec002" ref-type="sec">Materials and Methods</xref>”). We identified similar numbers of NBS domain-containing genes amongst these three genomes. We found 618, 650 and 508 NBS genes from <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA genomes respectively (<xref rid="pone.0121893.t001" ref-type="table">Table 1</xref>). Among them, 413, 499 and 484 NBS genes were predicted in the original gene annotations and 205, 151 and 24 NBS genes in <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA respectively were newly predicted in this project (<xref rid="pone.0121893.s012" ref-type="supplementary-material">S1 Table</xref>).</p><table-wrap id="pone.0121893.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.t001</object-id><label>Table 1</label><caption><title>Classification of <italic>Citrus</italic> NBS genes.</title></caption><alternatives><graphic id="pone.0121893.t001g" xlink:href="pone.0121893.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Class</th><th align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic></th><th align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> China</th><th align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> USA</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">CNL</td><td align="left" rowspan="1" colspan="1">117</td><td align="left" rowspan="1" colspan="1">113</td><td align="left" rowspan="1" colspan="1">60</td></tr><tr><td align="left" rowspan="1" colspan="1">TNL</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">30</td></tr><tr><td align="left" rowspan="1" colspan="1">NL</td><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">85</td></tr><tr><td align="left" rowspan="1" colspan="1">Others</td><td align="left" rowspan="1" colspan="1">351</td><td align="left" rowspan="1" colspan="1">398</td><td align="left" rowspan="1" colspan="1">333</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">618</td><td align="left" rowspan="1" colspan="1">650</td><td align="left" rowspan="1" colspan="1">508</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="t001fn001"><p>The CNL class of NBS genes contains three domains: Coiled Coil (CC), nucleotide-binding sites (NBS) and Leucine-Rich Repeat (LRR). The TNL class of NBS genes contains three domains: Toll-Interleukin receptor (TIR), nucleotide-binding sites (NBS) and Leucine-Rich Repeat (LRR). The NL class of NBS genes contains two domains: nucleotide-binding sites (NBS) and Leucine-Rich Repeat (LRR). The other NBS genes contain no LRR domain.</p></fn></table-wrap-foot></table-wrap><p>NBS genes could be classified into different classes based on their domain structures [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>]. We searched the Toll-Interleukin receptor (TIR) and Coiled Coil (CC) domains in the N-terminal region and the Leucine-Rich Repeat (LRR) domains in the C-terminal of NBS genes. We identified 117 NBS genes with CC-NBS-LRR (CNL) domains, 82 NBS genes with TIR-NBS-LRR (TNL) domains, 68 NBS genes with NBS-LRR (NL) and 351 others NBS genes without LRR domains (N, CN and TN) from <italic>C</italic>. <italic>clementina</italic> (<xref rid="pone.0121893.t001" ref-type="table">Table 1</xref>). We also identified 113 CNL, 77 TNL, 62 NL and 398 others NBS genes from <italic>C</italic>. <italic>sinensis</italic> China and 60 CNL, 30 TNL, 85 NL and 333 others NBS genes from <italic>C</italic>. <italic>sinensis</italic> USA (<xref rid="pone.0121893.t001" ref-type="table">Table 1</xref>). In comparison to other genomes, there are many more <italic>Citrus</italic> NBS genes without the LRR domain. There are only 43.2%, 38.8% and 34.4% NBS genes with existing LRR domains in <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA, respectively, while 72% NBS genes in <italic>Arabidopsis thaliana</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and 78% NBS genes in <italic>Populus trichocarpa</italic> [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>] have LRR domain.</p><p>The structures of TNL and CNL NBS genes are significantly different. The TNL NBS genes tend to have more introns than that of CNL NBS genes as previously found in <italic>Arabidopsis</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and <italic>Populus</italic> [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>]. The average numbers of introns in TNL NBS genes are 4.39, 4.70 and 4.47, while the numbers of introns in CNL NBS genes are 0.89, 1.15 and 1.38 in <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA respectively (<xref rid="pone.0121893.s001" ref-type="supplementary-material">S1 Fig</xref>.). The median numbers of introns of different types of NBS genes in the three <italic>Citrus</italic> genomes are similar: four in TNL NBS genes and one in CNL NBS genes.</p><p>The NBS genes were unevenly distributed in the <italic>Citrus</italic> scaffolds/chromosomes. There were 158, 217 and 110 NBS genes, totaling 78.3% of 618 NGS genes, distributed in scaffolds_3, _5 and _7 of <italic>C</italic>. <italic>clementina</italic>, respectively. There were 125, 121 and 75 NBS genes distributed in chromosomes 1, 3 and 5 of <italic>C</italic>. <italic>sinensis</italic> China. The majority (95 out of 107) of NBS genes with TIR domain in <italic>C</italic>. <italic>clementina</italic> was distributed in scaffold 3. There were 33 and 51 (out of 102) NBS genes with TIR domain in <italic>C</italic>. <italic>sinensis</italic> China distributed on chromosome 5 and chromosome unknown, respectively.</p></sec><sec id="sec018"><title>Phylogenetic and Clade Analysis of the <italic>Citrus</italic> NBS Genes</title><p>To avoid mutation saturation effects in nucleotide sequences over time which can lead to an underestimation of the number of mutation events due to higher substitution rates versus proteins, we used the protein sequences of NBS domains, which are the most conserved part of NBS genes, to construct a phylogenetic tree of NBS genes. We only selected NBS domain sequences longer than 200 amino acid residues and contain both P-loop and MHDV motifs. Finally, 442 <italic>C</italic>. <italic>clementina</italic>, 393 <italic>C</italic>. <italic>sinensis</italic> China, and 264 <italic>C</italic>. <italic>sinensis</italic> USA NBS domain sequences were used for phylogenetic analysis (<xref rid="pone.0121893.s013" ref-type="supplementary-material">S2 Table</xref>). A Maximum Likelihood (ML) phylogenetic tree of these 1,099 NBS genes was then constructed using <bold>FastTree</bold> [<xref rid="pone.0121893.ref031" ref-type="bibr">31</xref>]. As shown in <xref rid="pone.0121893.g001" ref-type="fig">Fig. 1</xref>, the un-rooted phylogenetic tree can be divided into three main groups. Most NBS genes containing the TIR domain are in one branch and the non-TIR NBS genes comprise the other two branches (<xref rid="pone.0121893.g001" ref-type="fig">Fig. 1</xref>), in which most of them contain a CC domain. Therefore, we denoted the three branches as TIR, CC1 and CC2.</p><fig id="pone.0121893.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g001</object-id><label>Fig 1</label><caption><title>Maximum likelihood phylogenetic tree of <italic>Citrus</italic> NBS-LRR genes constructed from multiple sequences alignment of NBS domain.</title><p>There are three gene groups in the phylogenetic tree: CC1 (pink fan-shape), CC2 (orange fan-shape) and TIR (light grey fan-shape). Clades were classified using PhyloPart as shown in alternating color and represented by alternative light yellow and grey colors in the middle circle. The outer circle shows species with <italic>C</italic>. <italic>clementina</italic> in red, <italic>C</italic>. <italic>sinensis</italic> China in green and <italic>C</italic>. <italic>sinensis</italic> USA in blue. Nodes with bootstrap support of 100 are indicated with grey circles on the tree.</p></caption><graphic xlink:href="pone.0121893.g001"></graphic></fig><p>There were 452, 382 and 265 <italic>Citrus</italic> NBS genes in CC1, CC2 and TIR groups, respectively. <xref rid="pone.0121893.s013" ref-type="supplementary-material">S2 Table</xref> lists the group classification of 1,099 <italic>Citrus</italic> NBS genes. For each sub-group, we constructed a new ML phylogenetic tree rooted by <italic>Streptomyces coelicolor</italic> protein P25941 (<xref rid="pone.0121893.s002" ref-type="supplementary-material">S2 Fig</xref>.). The NBS domains in CC groups, especially in the CC1 group, were relatively more diverged in sequence versus those in TIR group. The average Poisson corrected distances between sequence pairs within each group were 0.947, 0.827 and 0.663 for CC1, CC2 and TIR groups, respectively. We aligned the NBS domains sequences in each group using mafft with default parameters (—auto). The average percentages of identities for NBS domain sequences were 40%, 44% and 52% for CC1, CC2 and TIR groups, respectively.</p><p>The number of LRR domains varied amongst the three classes of <italic>Citrus</italic> NBS genes. The number of LRR domains in the TIR group is significantly higher than the numbers of LRR domains in the CC1 and CC2 groups (<xref rid="pone.0121893.s003" ref-type="supplementary-material">S3 Fig</xref>.). The majority of CC1 NBS genes have only one LRR domain. On other hand, the TIR NBS genes have 2.4 LRR domains on average. The most frequent type of LRR domain in the CC1 group is LRR_8 (PF13855), but LRR_1 (PF00560) is the most frequent domain in both the CC2 group and the TIR group (<xref rid="pone.0121893.s014" ref-type="supplementary-material">S3 Table</xref>).</p><p>NBS genes in the three <italic>Citrus</italic> genomes are polyphyletic amongst branches of phylogenetic tree as shown in the <xref rid="pone.0121893.g001" ref-type="fig">Fig. 1</xref>. We partitioned the NBS gene tree into clades based on NBS domain sequence distance using PhyloPart with distance threshold of 0.025. This resulted in 114 clades (with at least two genes) and eight orphan genes (<xref rid="pone.0121893.s013" ref-type="supplementary-material">S2 Table</xref>). The numbers of clades in each NBS gene group are similar. There are 39, 42, 32 clades in CC1, CC2 and TIR groups, respectively. We calculated average sequence similarities between NBS domains in each clade using <bold>alistat</bold> in <bold>SQUID</bold> (<ext-link ext-link-type="uri" xlink:href="http://selab.janelia.org/software.html">http://selab.janelia.org/software.html</ext-link>). The minimum average sequence similarity of NBS domains in each clade was 70%. The largest clade included 90 NBS genes. On average, there were 9.56 genes per clade. There were seven clades with more than 40 genes. Among these clades, four clades were in CC1; two clades were in CC2 and one clade was in the TIR group. Most of the clades contain NBS genes from the three <italic>Citrus</italic> genomes (<xref rid="pone.0121893.g001" ref-type="fig">Fig. 1</xref>, <xref rid="pone.0121893.t002" ref-type="table">Table 2</xref>). There are 89 clades with at least three NBS genes and 87 of these clades contain members from each of the three <italic>Citrus</italic> genomes. Together, these results imply that three <italic>Citrus</italic> genomes may have similar types of NBS genes.</p><table-wrap id="pone.0121893.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.t002</object-id><label>Table 2</label><caption><title>Types of <italic>Citrus</italic> NBS gene clades.</title></caption><alternatives><graphic id="pone.0121893.t002g" xlink:href="pone.0121893.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">NBS genes from Citrus genomes</th><th align="left" rowspan="1" colspan="1">Number of Clade</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China, <italic>C</italic>. <italic>sinensis</italic> USA</td><td align="left" rowspan="1" colspan="1">87</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China</td><td align="left" rowspan="1" colspan="1">13</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> USA</td><td align="left" rowspan="1" colspan="1">7</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> China, <italic>C</italic>. <italic>sinensis</italic> USA</td><td align="left" rowspan="1" colspan="1">5</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic>,</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> China,</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> USA</td><td align="left" rowspan="1" colspan="1">1</td></tr></tbody></table></alternatives></table-wrap><p>To understand the evolutionary dynamics maintaining functional constraint while maintaining so many gene family members, we tested for evidence of accelerated sequence evolution (positive selection) by searching for positively-selected sites within each clade. We detected positively selected sites in 53.5% (61) of clades. There were 18, 24 and 19 clades with sites under positive selection in CC1, CC2 and TIR groups, respectively. In total, we detected 541 positively selected sites. Consistent with previous reports [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>], there are more positively selected sites in the C-terminal region of NBS genes (LRR domains) than those in the N-terminal and NBS domains (<xref rid="pone.0121893.s015" ref-type="supplementary-material">S4 Table</xref>). Approximately 56.4% (305 out of 541) of the positively selected sites were located in the C-terminal region (LRR domain). While most clades have positively selected sites in the LRR domain, there were clades with a greater number of positively selected sites within the NBS domain. For example, there were 46 positively selected sites in Clade_1260 that are located in NBS domain while there are only two positively selected sites in LRR domain (<xref rid="pone.0121893.s004" ref-type="supplementary-material">S4 Fig</xref>.).</p></sec><sec id="sec019"><title>Mapping Re-sequencing Reads Showed the Conservation of NBS Genes in <italic>Citrus</italic> Genomes</title><p>To verify that <italic>C</italic>. <italic>sinensis</italic> (sweet orange) is a backcross hybrid of <italic>C</italic>. <italic>maxima</italic> (pomelo) and <italic>C</italic>. <italic>reticulata</italic> (mandarin) [<xref rid="pone.0121893.ref019" ref-type="bibr">19</xref>], Xu <italic>et al</italic>. previously re-sequenced three pomelo cultivars and three mandarin cultivars and showed that SSR and SNP markers in the <italic>C</italic>. <italic>sinensis</italic> genome derive from the pomelo and mandarin genomes at an approximately ratio of 1:3 [<xref rid="pone.0121893.ref019" ref-type="bibr">19</xref>]. Since our comparison of <italic>C</italic>. <italic>sinensis</italic> and <italic>C</italic>. <italic>clementina</italic> genomes showed that most of their NBS genes are very similar in sequence, we investigated the relationships amongst the NBS genes of <italic>C</italic>. <italic>sinensis</italic> (sweet orange), <italic>C</italic>. <italic>maxima</italic> (pomelo), and <italic>C</italic>. <italic>reticulata</italic> (mandarin) to compare and contrast how mutations accumulated and determine if there are differences in sequence, and hence potentially functional, constraints amongst NBS gene types in <italic>Citrus</italic> genomes.</p><p>We mapped the re-sequenced reads of three <italic>C</italic>. <italic>maxima</italic> (pomelo) genomes and three <italic>C</italic>. <italic>reticulata</italic> (mandarin) genomes to hybridized genome of <italic>C</italic>. <italic>sinensis</italic> (sweet orange) China using BWA with mem method [<xref rid="pone.0121893.ref041" ref-type="bibr">41</xref>] and calculated the coverage of each NBS gene in <italic>C</italic>. <italic>sinensis</italic> China. We detected 524 of 650 <italic>C</italic>. <italic>sinensis</italic> China NBS genes with more than 50% coverage in all of the six re-sequenced genomes. Additionally, 249 of 650 <italic>C</italic>. <italic>sinensis</italic> China NBS genes with more than 85% coverage in all six resequenced genomes and 592 (91%) of <italic>C</italic>. <italic>sinensis</italic> China NBS genes with coverage over 85% in at least one of the re-sequenced <italic>Citrus</italic> genomes. Using a cutoff of 15% coverage as missing calls over the entire length of the gene, only 25 (<4%) of <italic>C</italic>. <italic>sinensis</italic> China NBS genes seem to be deleted in at least one of the six re-sequenced genomes. No NBS genes were lost amongst all of the six re-sequenced genomes.</p><p>Re-sequenced reads tended to map to exons rather than to introns (<xref rid="pone.0121893.s016" ref-type="supplementary-material">S5 Table</xref>). For example, the entire coding region of NBS gene Cs9g13310 was 3,693 bp in length and covered by all reads from the re-sequenced samples. However, the 16 kb introns were not consistently present in all of the re-sequenced reads. When we used exon presence instead of whole gene coverage, 403 out of 650 (62%) <italic>C</italic>. <italic>sinensis</italic> China NBS genes have more than 85% coverage in all six re-sequenced <italic>Citrus</italic> genomes and 645 (99%) <italic>C</italic>. <italic>sinensis</italic> China NBS genes have more than 85% coverage in at least one of the resequenced <italic>Citrus</italic> genomes. Most <italic>C</italic>. <italic>sinensis</italic> China NBS genes have high exon coverage on both the <italic>C</italic>. <italic>maxima</italic> (pomelo) and <italic>C</italic>. <italic>reticulata</italic> (mandarin) genomes (<xref rid="pone.0121893.s005" ref-type="supplementary-material">S5 Fig</xref>.). There were only three <italic>C</italic>. <italic>sinensis</italic> China NBS genes, Cs1g02140.1, Cs6g02120.1 and Cs7g02220.1, with low (<40%) exon coverage in all three <italic>C</italic>. <italic>maxima</italic> genomes, but high coverage was achieved in at least two of the three <italic>C</italic>. <italic>reticulata</italic> genomes (<xref rid="pone.0121893.s005" ref-type="supplementary-material">S5 Fig</xref>.). Our mapping results showed that most NBS genes of <italic>C</italic>. <italic>sinensis</italic> (sweet orange) have a corresponding copy in the re-sequenced <italic>C</italic>. <italic>maxima</italic> (pomelo) and <italic>C</italic>. <italic>reticulata</italic> (mandarin) genomes.</p></sec><sec id="sec020"><title>Motif Patterns in <italic>Citrus</italic> NBS Genes</title><p>To further examine the gene structures of three <italic>Citrus</italic> NBS gene groups, we searched the motifs in the sequences of each group using <bold>MEME</bold>. <xref rid="pone.0121893.g002" ref-type="fig">Fig. 2</xref> lists consensus sequences of the top 20 motifs identified by MEME from the CC1, CC2 and TIR <italic>Citrus</italic> NBS gene groups.</p><fig id="pone.0121893.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g002</object-id><label>Fig 2</label><caption><title>Architecture of NBS gene motifs in <italic>Citrus</italic>.</title><p>The consensus sequences of the top 20 motifs identified by MEME for each group of <italic>Citrus</italic> NBS genes are listed (A: CC1 group, B: CC2 group, C: TIR group). The motifs in blue are identified from N-terminal (CC or TIR domain). The motifs in red are identified from NBS domain. The motifs in green are identified from C-terminal (LRR domain). The motifs in black are identified from other region of NBS genes. The bold letters highlight the conserved motifs in <italic>Arabidopsis</italic>. The motifs underlined in green were repeated in C-terminal (LRR domain).</p></caption><graphic xlink:href="pone.0121893.g002"></graphic></fig><p>The motifs in the N-terminal of CC1 and CC2 groups showed little similarity to each other. The MEME found two motifs from CC1 group. One of them, EEQQQMRRLNQVQGWLSRVEA, was present in 356 of 452 CC1 NBS genes. The other motif, GSQEIDKLCLGGYCSKNCKSSYKFGKKVA, was present in 189 of 452 CC1 NBS genes. Both motifs have high level of sequence similarities with the CC NBS genes of <italic>Arabidopsis</italic> (<xref rid="pone.0121893.s006" ref-type="supplementary-material">S6 Fig</xref>. A). The MEME found three motifs from CC2 group. Two of the motifs, KKLTNMLEMIKAVLDDAEEKQ and AVKLWLGKLKDAAYDVEDVLDEFQTEALR were identified in 325 and 340 of 382 CC2 NBS genes, respectively. The motif AVKLWLGKLKDAAYDVEDVLDEFQTEALR has high levels of similarity with the motif identified from the CC NBS genes of <italic>Oryza sativa japonica</italic> (Japonica rice) [<xref rid="pone.0121893.ref045" ref-type="bibr">45</xref>] (<xref rid="pone.0121893.s006" ref-type="supplementary-material">S6 Fig</xref>. B). The different motifs in CC1 and CC2 NBS genes implied that they may be from different evolutionary origins.</p><p>For the TIR group, the MEME identified five motifs in the N-terminal. They can be found in 76% to 90% of 265 TIR <italic>Citrus</italic> NBS genes. The first four motifs YDVFLSFRGEDTRDNFTSHLY, AIEASAISVIIFSEGYASSRWCLDELVKI, GQIVIPVFYRVDPSDVRKQTG, and ENPEKVQKWRDALKEA were very similar to the TIR1-4 motifs identified from TIR NBS genes in <italic>Arabidopsis</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and <italic>Populus trichocarpa</italic> [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>] (<xref rid="pone.0121893.s006" ref-type="supplementary-material">S6 Fig</xref>. C).</p><p>The MEME algorithm identified nine motifs from NBS domains of CC2 and TIR groups, which are similar to the motif structures of NBS domains in <italic>Arabidopsis thaliana</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and <italic>Populus trichocarpa</italic> NBS genes [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>]. MEME results also showed that the motif structure of NBS domain of CC1 groups is slightly different from those of CC2 and TIR groups. Eleven motifs were identified from NBS domains of CC1 group with two extra RNBS motifs (RNBS-F, RNBS-G) between GLPL and MHDV motifs. Five motifs: P-loop, Kinase-2, RNBS-B, GLPL and MHDV, showed high levels of sequence similarity amongst NBS genes of three groups, which were also similar to those motifs from <italic>Arabidopsis</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and <italic>Populus trichocarpa</italic> NBS genes [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>]. The MHDV motif in <italic>Citrus</italic> was often slightly modified to MHDL, as found in <italic>Arabidopsis</italic> and <italic>Populus</italic> previously. Meanwhile, the motifs RNBS-A, RNBS-C, and RNBS-D were quite dissimilar to each other amongst three groups.</p><p>The LRR domains in the C-terminal of NBS genes usually have high sequence diversity as they play a role in recognizing pathogen virulence proteins. The MEME algorithm identified five, six, and four motifs from LRR domains of CC1, CC2 and TIR NBS genes, respectively. Most of the motifs contained LxxL repeats. The LRR motifs from different groups were highly variable in sequence, which implied that the three groups of NBS genes play different roles in the <italic>Citrus</italic> immune system. Some LRR motifs had repeated several times in the same gene. Examples include: motif APNLKSLEVSSCxxMEEIISV found in 367 of 452 CC1 group with an average 4.8 motifs per gene; motif IKTLPESVCELYNLQTLDLEGCRRLKKLP found in 332 of 382 CC2 group with an average 2.8 times per gene and motifs CKRLKSLPSSLCKLKSLGxLNLSGCSNLE and GNISELFLDGTAIEELPSSIE found in 208 and 209 of 265 TIR group with 2.5 and 2.1 times per gene, respectively.</p></sec><sec id="sec021"><title>Analysis of NBS Gene Clusters in <italic>Citrus</italic> Genomes</title><p>The majority of <italic>Citrus</italic> NBS genes were physically clustered in genome (<xref rid="pone.0121893.t003" ref-type="table">Table 3</xref>). 525 of 618 NBS genes in <italic>C</italic>. <italic>clementina</italic> were found in 108 clusters and 500 of 650 NBS genes in <italic>C</italic>. <italic>sinensis</italic> China were found in 126 clusters (<xref rid="pone.0121893.s013" ref-type="supplementary-material">S2 Table</xref>). Although the assembly of <italic>C</italic>. <italic>sinensis</italic> USA is more fragmented, there still are 207 of 508 NBS genes present in 72 clusters. The largest number of gene clusters in C. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA contain 55, 37 and 13 NBS genes, respectively.</p><table-wrap id="pone.0121893.t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.t003</object-id><label>Table 3</label><caption><title><italic>Citrus</italic> NBS genes in clusters.</title></caption><alternatives><graphic id="pone.0121893.t003g" xlink:href="pone.0121893.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Citrus species</th><th align="left" rowspan="1" colspan="1">Total NBS genes</th><th align="left" rowspan="1" colspan="1">Number of NBS genes in cluster</th><th align="left" rowspan="1" colspan="1">Number of clusters</th><th align="left" rowspan="1" colspan="1">Percent of NBS genes in cluster (%)</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>clementina</italic></td><td align="left" rowspan="1" colspan="1">618</td><td align="left" rowspan="1" colspan="1">525</td><td align="left" rowspan="1" colspan="1">108</td><td align="char" char="." rowspan="1" colspan="1">84.9</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> China</td><td align="left" rowspan="1" colspan="1">650</td><td align="left" rowspan="1" colspan="1">500</td><td align="left" rowspan="1" colspan="1">126</td><td align="char" char="." rowspan="1" colspan="1">76.9</td></tr><tr><td align="left" rowspan="1" colspan="1"><italic>C</italic>. <italic>sinensis</italic> USA</td><td align="left" rowspan="1" colspan="1">508</td><td align="left" rowspan="1" colspan="1">207</td><td align="left" rowspan="1" colspan="1">72</td><td align="char" char="." rowspan="1" colspan="1">40.7</td></tr></tbody></table></alternatives></table-wrap><p>Most clusters contain NBS genes from the same group. In <italic>C</italic>. <italic>clementina</italic>, there were 38 clusters with NBS genes of CC1 group, 40 clusters with NBS genes of CC2 group and 23 clusters with NBS genes of TIR group. Only seven out of 108 clusters contain the NBS genes from two or three groups. In <italic>C</italic>. <italic>sinensis</italic> China, there were 41 clusters with NBS genes of CC1 group, 44 clusters with NBS genes of CC2 group and 26 clusters with NBS genes of TIR group. There were 15 clusters containing NBS genes from two or three groups in <italic>C</italic>. <italic>sinensis</italic> China. In <italic>C</italic>. <italic>sinensis</italic> USA, there were 27 clusters with NBS genes of CC1 group, 25 clusters with NBS genes of CC2 group and 15 clusters with NBS genes of TIR group. Only five clusters contain NBS genes from different groups in <italic>C</italic>. <italic>sinensis</italic> USA. The lower number of clusters in <italic>C</italic>. <italic>sinensis</italic> USA may due to its fragmented assembly.</p><p>The sequences of NBS genes within clusters are much more similar to each other than those between clusters (T-test p value < 2.2e-16). The mean identities between genes within and between clusters were 0.628 and 0.393 respectively (<xref rid="pone.0121893.g003" ref-type="fig">Fig. 3</xref>). Furthermore, NBS genes in the same cluster tend to be in the same strand, which indicates that the NBS genes in the clusters are due to tandem duplication. We detected 204 pairs of tandem duplications in <italic>C</italic>. <italic>clementina</italic> and 217 pairs of tandem duplications in <italic>C</italic>. <italic>sinensis</italic> China using MCScanX [<xref rid="pone.0121893.ref039" ref-type="bibr">39</xref>]. The numbers of tandem duplications within CC1 and CC2 groups are much greater than that within the TIR group. Among 204 pairs in <italic>C</italic>. <italic>clementina</italic>, 90 and 85 pairs were present in the CC1 and CC2 groups and only 29 pairs from TIR group. Among 217 tandem gene pairs in <italic>C</italic>. <italic>sinensis</italic> China, 88 and 87 pairs were present in the CC1 and CC2 groups and only 42 pairs in the TIR group. The fewer tandem duplications of NBS genes in the TIR group may be the reason that there are fewer clusters of TIR groups in <italic>Citrus</italic> genomes.</p><fig id="pone.0121893.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g003</object-id><label>Fig 3</label><caption><title>Percentage identity distribution of NBS-encoding genes in <italic>C</italic>. <italic>clementina</italic> (Cc) and <italic>C</italic>. <italic>sinensis</italic> China (Csc).</title><p>Green lines indicate pairwise identity distribution of inter-cluster NBS-encoding genes. Orange lines indicate intro-cluster NBS-encoding genes.</p></caption><graphic xlink:href="pone.0121893.g003"></graphic></fig><p>We identified 254 and 246 gene conversion events from 116 NBS gene clusters in <italic>C</italic>. <italic>clementina</italic> and 144 NBS gene clusters in <italic>C</italic>. <italic>sinensis</italic> China, respectively (<xref rid="pone.0121893.t004" ref-type="table">Table 4</xref>). The gene conversion events in <italic>C</italic>. <italic>sinensis</italic> USA is much less due to fragmented assembly. It is interesting that most of the conversion events (483 out of 520) were identified from the relatively small clusters with less than 10 NBS genes. Among these conversions, 119 events located in N-terminal, 178 in NBS domains and 223 in C-terminal, which indicating that there was no significantly bias in the location of conversion. Most gene conversions were between genes from the same group. We identified 101, 184 and 229 conversion events from NBS genes in TIR, CC1 and CC2 groups respectively and only four conversion events were identified between NBS genes of TIR and CC2 groups. While we could identify similar total numbers of conversion events in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China, we found almost double such events in <italic>C</italic>. <italic>sinensis</italic> China TIR clusters versus in <italic>C</italic>. <italic>clementina</italic> TIR clusters. Meanwhile, we also found many more conversion events in <italic>C</italic>. <italic>clementina</italic> CC2 clusters than that in <italic>C</italic>. <italic>sinensis</italic> China CC2 clusters (<xref rid="pone.0121893.t004" ref-type="table">Table 4</xref>).</p><table-wrap id="pone.0121893.t004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.t004</object-id><label>Table 4</label><caption><title>Gene conversion events found in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China.</title></caption><alternatives><graphic id="pone.0121893.t004g" xlink:href="pone.0121893.t004"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Citrus species</th><th align="left" rowspan="1" colspan="1">Cluster type</th><th align="left" rowspan="1" colspan="1">Number of clusters (with conversion events)</th><th align="left" rowspan="1" colspan="1">Number of conversion events</th><th align="left" rowspan="1" colspan="1">Affected genes</th><th align="left" rowspan="1" colspan="1">Mean size of conversion tracts</th></tr></thead><tbody><tr><td rowspan="3" align="left" colspan="1"><italic>C</italic>. <italic>clementina</italic></td><td align="left" rowspan="1" colspan="1">TIR</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">111</td></tr><tr><td align="left" rowspan="1" colspan="1">CC1</td><td align="left" rowspan="1" colspan="1">43</td><td align="left" rowspan="1" colspan="1">86</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">137</td></tr><tr><td align="left" rowspan="1" colspan="1">CC2</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">133</td><td align="left" rowspan="1" colspan="1">83</td><td align="left" rowspan="1" colspan="1">146</td></tr><tr><td rowspan="3" align="left" colspan="1"><italic>C</italic>. <italic>sinensis</italic> China</td><td align="left" rowspan="1" colspan="1">TIR</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">69</td><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">137</td></tr><tr><td align="left" rowspan="1" colspan="1">CC1</td><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">89</td><td align="left" rowspan="1" colspan="1">48</td><td align="left" rowspan="1" colspan="1">162</td></tr><tr><td align="left" rowspan="1" colspan="1">CC2</td><td align="left" rowspan="1" colspan="1">57</td><td align="left" rowspan="1" colspan="1">92</td><td align="left" rowspan="1" colspan="1">75</td><td align="left" rowspan="1" colspan="1">116</td></tr></tbody></table></alternatives></table-wrap></sec><sec id="sec022"><title>Analysis of NBS Orthologs</title><p>We identified 719 <italic>Citrus</italic> NBS gene pairs of orthologs amongst <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA. 270 orthologs were shared between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China (<xref rid="pone.0121893.s007" ref-type="supplementary-material">S7 Fig</xref>.); 227 orthologous gene pairs were shared between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> USA and 222 orthologous gene pairs were shared between <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA. The percentages of identities between orthologous genes range from 53.8% to 100% (mean of 93.6%, median of 96.8%). The percentages of identities between orthologous genes from TIR group (TNL and TN, 88.42 ± 10.04) were significantly (T-test: P<2.6e-6) lower than those from CC groups (CNL and CN, 92.92 ± 7.03). The nonsynonymous divergence (dN) values of the 719 orthologs ranged from 1.46e-6 to 0.31 (mean of 0.03, median of 0.08) and the synonymous divergence (dS) values ranged from 1.14e-5 to 0.75 (mean of 0.05, median of 0.019) (<xref rid="pone.0121893.g004" ref-type="fig">Fig. 4A, B</xref>). The dN/dS rate ratios ranged from 0.001 to 50 with a median of 0.62. Similar to <italic>Arabidopsis</italic>, the dN and dS of orthologs from TIR group (TNL and TN) were higher than those from CC groups (CNL and CN). The dN values of orthologs from TIR group ranged from 0.001 to 0.24 (mean of 0.038, median of 0.017) while those of orthologs from CC group ranged from 1.4e-6 to 0.23 (mean of 0.025, median of 0.007). The dS rate of orthologs from TIR group ranged from 3.27e-5 to 0.76 (mean of 0.078, median of 0.029) while those of orthologs from CC group ranged from 1e-5 to 0.67 (mean of 0.04, median of 0.016). However, the dN/dS rate ratios of orthologs from TIR group were relatively lower than that of orthologs from CC group (<xref rid="pone.0121893.g004" ref-type="fig">Fig. 4C</xref>).</p><fig id="pone.0121893.g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g004</object-id><label>Fig 4</label><caption><title>Divergence of orthologous NBS-encoding genes amongst <italic>Citrus</italic> species.</title><p><bold>A-C.</bold> Divergence of different classes. <bold>D-F.</bold> Divergence of singletons and cluster members. dN indicates the nonsynonymous mutation rate (rate of nonsynonymous mutations per nonsynonymous site). dS indicates the synonymous mutation rate. CNL, CN, TNL, TN, NL and N indicate orthologs between two proteins belonging to same cluster with CC-NBS-LRR domains (CNL), CC-NBS domains (CN), TIR-NBS-LRR domains (TNL), TIR-NBS domains (TN), NBS-LRR domains (NL), NBS domain (N), respectively. The others indicate orthologs between different clusters.</p></caption><graphic xlink:href="pone.0121893.g004"></graphic></fig><p>The dN and dS rates as well as dN/dS rate ratios of ortholog residues in clusters are generally greater than those of orthologous singletons (<xref rid="pone.0121893.g004" ref-type="fig">Fig. 4D-F</xref>). The median dN/dS rate ratio of 148 orthologous singletons was 0.55 while that of 361 orthologs in clusters was 0.63. There were 124 orthologs with dN/dS rate ratios above 1. Variation in these orthologs may be shaped by positive selection pressure. Most of the positive selected orthologs (87.9%, 109 out of 124) belonged to CC groups (49 orthologs from CC1 group and 50 from CC2 group).</p><p>We also detected 38 NBS gene syntenic blocks between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China using MCScanX [<xref rid="pone.0121893.ref039" ref-type="bibr">39</xref>]. Totally, there are 416 syntenic orthologs of NBS genes between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China (<xref rid="pone.0121893.s008" ref-type="supplementary-material">S8 Fig</xref>.). On average, there are 11 genes per block. The biggest block contains 45 NBS genes and was located in scaffold_5 of <italic>C</italic>. <italic>clementina</italic> and chromosome 3 of <italic>C</italic>. <italic>sinensis</italic> China. There were 80, 77 and 71 syntenic NBS genes in scaffold_5 of <italic>C</italic>. <italic>clementina</italic>/chromosome 5 of <italic>C</italic>. <italic>sinensis</italic> China, scaffold_3 of <italic>C</italic>. <italic>clementina</italic>/chromosome 5 of <italic>C</italic>. <italic>sinensis</italic> China and scaffold_7 of <italic>C</italic>. <italic>clementina</italic>/chromosome 1 of <italic>C</italic>. <italic>sinensis</italic> China. Furthermore, there were 68, 66 and 8 syntenic NBS genes in the unknown chromosome of <italic>C</italic>. <italic>sinensis</italic> China corresponding to NBS genes in scaffold_3, _5 and _7 of <italic>C</italic>. <italic>clementina</italic>, respectively.</p></sec><sec id="sec023"><title>Analysis of Conserved NBS Gene Clusters</title><p>We identified 118 pairs of conserved NBS gene clusters between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China (<xref rid="pone.0121893.s009" ref-type="supplementary-material">S9 Fig</xref>.). There were 19 NBS gene clusters completely conserved in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China. For example, all 14 NBS genes in cluster CL225 of <italic>C</italic>. <italic>sinensis</italic> China have orthologs in the cluster CL168 <italic>C</italic>. <italic>clementina</italic> which contains nine NBS genes and vice versa. Some clusters (in both <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China) have several corresponding conserved clusters. This may due to either the genome arrangement or the incomplete assembly of genomes. Furthermore, there were 19 clusters of <italic>C</italic>. <italic>clementina</italic> with no conserved clusters in <italic>C</italic>. <italic>sinensis</italic> China and 18 clusters of <italic>C</italic>. <italic>sinensis</italic> China with no conserved clusters in <italic>C</italic>. <italic>clementina</italic>.</p><p>The conserved clusters provided additional insights into NBS gene evolution within and between <italic>Citrus</italic> genomes. For example, cluster CL142 (9 NBS genes) in <italic>C</italic>. <italic>clementina</italic> and cluster CL282 (7 NBS genes) in <italic>C</italic>. <italic>sinensis</italic> China were highly conserved. The phylogenetic tree of these 16 genes suggested division into two subgroups (depicted in blue and orange color in <xref rid="pone.0121893.g005" ref-type="fig">Fig. 5A</xref>). This division suggests two ancestral genes for this conserved cluster. Using the phylogenetic tree as a framework, we reconstructed the evolutionary history of these two clusters. There were several tandem duplication events in the evolutionary history of the conserved clusters, and extra tandem duplication was observed in <italic>C</italic>. <italic>clementina</italic> after it separated from <italic>C</italic>. <italic>sinensis</italic> (<xref rid="pone.0121893.g005" ref-type="fig">Fig. 5B</xref>). Two NBS genes were lost in <italic>C</italic>. <italic>sinensis</italic> and one NBS gene was lost in <italic>C</italic>. <italic>clementina</italic>. Furthermore, there was a recombination event within <italic>C</italic>. <italic>clementina</italic>.</p><fig id="pone.0121893.g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g005</object-id><label>Fig 5</label><caption><title>Duplication histories of one of the conserved NBS-encoding genes cluster in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic>.</title><p><bold>A.</bold> Phylogenetic tree of the NBS-encoding genes; <bold>B.</bold> Orthologs in the conserved cluster of <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic>. Black connectors indicate tandem duplications and red dotted connectors indicate gene loss. A gene recombination event occurred in the ancestor of <italic>C</italic>. <italic>clementina</italic>, as indicated with an orange connector.</p></caption><graphic xlink:href="pone.0121893.g005"></graphic></fig></sec><sec id="sec024"><title>Mutations and Transposons in <italic>Citrus</italic> NBS Genes</title><p>Plant NBS genes are continuously evolving. Sequence variation and structural constraints are shaped by gene birth-and-death processes [<xref rid="pone.0121893.ref046" ref-type="bibr">46</xref>]. Besides possible interallelic recombination and gene conversion, gene mutations and transposable elements appear to play important roles in NBS gene evolution.</p><p>We compared the DNA sequences of NBS genes from different <italic>Citrus</italic> genomes to identify mutations. Approximately half of <italic>Citrus</italic> NBS genes have mutations maintained amongst corresponding orthologs. When a mutation took place in an exon and resulted in stop-codon gaining or frame-shift, the target gene often became a pseudogene. For example, Cs1g18610.1_cc_32, Cs1g18610.1 and orange1.1g003367m are orthologs of <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA. There were 2 stop-codon gaining mutations in Cs1g18610.1 (<xref rid="pone.0121893.s010" ref-type="supplementary-material">S10 Fig</xref>. A), which means that Cs1g18610.1 might become a pseudo (nonfunctional) gene after gaining these mutations.</p><p>LTR retrotransposons are widespread in eukaryotic genomes, especially plant genomes. We predicted 19014, 6296 and 1479 LTR retrotransponsons with typical LTR characters in the draft genomes of <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA, respectively. Then, we filtered out LTR retrotransponsons with low similarity with known TE proteins using BLASTX with e-value greater than 1e-5. Finally, 4920, 3726 and 1240 LTR retrotransponsons remained in <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA, respectively. We identified 33, 32 and 4 NBS genes that were inserted with LTR retrotransposons in <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA respectively (<xref rid="pone.0121893.s017" ref-type="supplementary-material">S6 Table</xref>). Most of these genes will likely become pseudogenes due to these insertions. For an instance, orange1.1g043039m, orange1.1g043039m_cc_116 and orange1.1g043039m_csc_123 were orthologs of <italic>C</italic>. <italic>sinensis</italic> USA, <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China respectively (<xref rid="pone.0121893.s010" ref-type="supplementary-material">S10 Fig</xref>. B). The results of gene structure analysis showed that structure in orange1.1g043039m in <italic>C</italic>. <italic>sinensis</italic> USA seems relatively well maintained. But there is an about 10 kb fragment of LTR-retrotransposons in the corresponding homologous gene in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China (<xref rid="pone.0121893.s010" ref-type="supplementary-material">S10 Fig</xref>. B). We identified both orange1.1g043039m_cc_116 and orange1.1g043039m_csc_123 as pseudogenes using PseudoPipe [<xref rid="pone.0121893.ref037" ref-type="bibr">37</xref>].</p></sec><sec id="sec025"><title>Experimental Validation of One NBS Gene in <italic>Citrus</italic> Genomes</title><p>We validated the orthologs of a conserved NB gene, Cs1g09350.1, in a wide range of <italic>Citrus</italic> species. Cs1g09350.1 conserved in the 3 sequenced <italic>Citrus</italic> genomes. It is a CNL NBS gene and has 5 exons with about 4 kb in length. We designed the primers targeting about 3.5 kb fragment and could amplify a 3.5 kb fragment using PCR in different <italic>Citrus</italic> species, including <italic>C</italic>. <italic>sinensis</italic> (sweet orange), <italic>C</italic>. <italic>clementina</italic> (clementine), <italic>C</italic>. <italic>japonica</italic> (kumquat), <italic>C</italic>. <italic>sinensis</italic> Navelina (navel orange), <italic>C</italic>. <italic>maxima</italic> (pomelo), <italic>C</italic>. <italic>aurantiifolia</italic> (lime), and <italic>C</italic>. <italic>aurantium</italic> (sour orange). The orthologs from <italic>C</italic>. <italic>sinensis</italic> and <italic>C</italic>. <italic>clementina</italic> were identical as expected. There are a few mutations in the orthologs of <italic>C</italic>. <italic>japonica</italic>, <italic>C</italic>. <italic>sinensis</italic> Navelina, and <italic>C</italic>. <italic>maxima</italic>, but the orthologs of <italic>C</italic>. <italic>aurantiifolia</italic> and <italic>C</italic>. <italic>aurantium</italic> should be undergoing a pseudogenization process. There was a deletion plus several mutations in the NBS domain in <italic>C</italic>. <italic>aurantiifolia</italic>. There was an eight-base deletion in the LRR domain in <italic>C</italic>. <italic>aurantiifolia</italic> and <italic>C</italic>. <italic>aurantium</italic>. Most <italic>C</italic>. <italic>aurantiifolia</italic> mutations were shared with <italic>C</italic>. <italic>aurantium</italic> (<xref rid="pone.0121893.s011" ref-type="supplementary-material">S11 Fig</xref>.), and these mutations were not shared with other species. The results suggest that in <italic>C</italic>. <italic>aurantiifolia</italic> and <italic>C</italic>. <italic>aurantium</italic>, this gene is likely derived from the same common ancestor and was inherited as a pseudogene.</p></sec></sec><sec sec-type="conclusions" id="sec026"><title>Discussion</title><sec id="sec027"><title>Similar NBS Genes in Hybrid <italic>Citrus sinensis</italic> and Original <italic>Citrus clementina</italic> Genome</title><p>After carefully reannotating the <italic>Citrus</italic> genome sequences, we found similar numbers of NBS genes in <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China. There are slightly fewer NBS genes in <italic>C</italic>. <italic>sinensis</italic> USA, possibly due to the more fragmented assembly of this genome. In phylogenetic tree using NBS domains, the NBS genes from three different <italic>Citrus</italic> genomes are spread across clades. After partitioning NBS genes on the tree into clades, we found that 97.7% of the clades containing at least three genes had members from each of three <italic>Citrus</italic> genomes. This pattern suggests that these three <italic>Citrus</italic> genomes have similar types of NBS genes derived from common ancestors. Because <italic>C</italic>. <italic>sinensis</italic> is the hybrid of <italic>C</italic>. <italic>reticulata</italic> (mandarin) and <italic>C</italic>. <italic>maxima</italic> (pomelo), it should be heterozygous and some of NBS genes in <italic>C</italic>. <italic>sinensis</italic> are expected have different genetic distances from NBS genes in <italic>C</italic>. <italic>clementina</italic>. This would support previous observations of NBS genes in F1 progeny of <italic>Poncirus trifoliata</italic> (trifoliate orange) and <italic>Citrus reticulata</italic> (tangerine). However, this is clearly not the case.</p><p>Furthermore, we mapped the re-sequenced reads of three <italic>C</italic>. <italic>maxima</italic> (pomelo) genomes and three <italic>C</italic>. <italic>reticulata</italic> (Mandarin orange) genomes onto the genome of <italic>C</italic>. <italic>sinensis</italic> China. In this case, 62% of <italic>C</italic>. <italic>sinensis</italic> China NBS genes have a copy present in all six resequenced genomes and 99% of <italic>C</italic>. <italic>sinensis</italic> China NBS genes have a copy in at least one of the re-sequenced genomes. The mapping results confirmed that a significant percentage of NBS genes of hybrid <italic>C</italic>. <italic>sinensis</italic> genomes have corresponding homologous genes in both the <italic>C</italic>. <italic>maxima</italic> and the <italic>C</italic>. <italic>reticulata</italic> genomes. Because the reference genome sequence of <italic>C</italic>. <italic>maxima</italic> is not yet available, the total number of NBS genes in <italic>C</italic>. <italic>maxima</italic> genomes is still unknown. However, we can at least conclude from the mapping of re-sequenced genomes that the <italic>C</italic>. <italic>maxima</italic> genome has homologous copies of NBS genes in <italic>C</italic>. <italic>reticulata</italic> and <italic>C</italic>. <italic>sinensis</italic> genomes. The homologous NBS genes in <italic>C</italic>. <italic>maxima</italic> and <italic>C</italic>. <italic>reticulata</italic> may be the reason that NBS genes in their hybrid <italic>C</italic>. <italic>sinensis</italic> are similar to those in <italic>C</italic>. <italic>clementina</italic> in this study.</p></sec><sec id="sec028"><title>Three Groups of <italic>Citrus</italic> NBS Genes</title><p>We identified 442, 393 and 264 genes with full length NBS domains from <italic>C</italic>. <italic>clementina</italic>, <italic>C</italic>. <italic>sinensis</italic> China and <italic>C</italic>. <italic>sinensis</italic> USA reference genomes, respectively. There are also many genes with short NBS domains in three <italic>Citrus</italic> genomes. The <italic>Citrus</italic> NBS genes can be divided into three groups according to the phylogenetic tree of NBS domains: two of them without TIR domain (CC1 and CC2 groups) and the other group with TIR domain (TIR group). The number of non-TIR NBS genes is three times of the number of TIR NBS genes. In most of the TIR NBS genes, we can find the LRR domains defined in Pfam database [<xref rid="pone.0121893.ref047" ref-type="bibr">47</xref>]. We only can identify LRR domains in small part of non-TIR NBS genes using Pfam LRR domain definition. However, we can find the LxxLs repeats in most of the non-TIR NBS genes as shown in motifs from MEME. This implied that there may be other types of LRR domain in <italic>Citrus</italic> non-TIR NBS genes.</p><p>Our motif analyses also showed that motifs of TIR domains in <italic>Citrus</italic> TIR NBS genes are similar to those of TIR domains in <italic>Arabidopsis</italic> [<xref rid="pone.0121893.ref006" ref-type="bibr">6</xref>] and <italic>Populus trichocarpa</italic> [<xref rid="pone.0121893.ref044" ref-type="bibr">44</xref>] TIR NBS genes. Furthermore, the motifs of <italic>Citrus</italic> CC domains in CC1 NBS genes are similar to those motifs of CC domains in <italic>Arabidopsis</italic> CC NBS genes and the motifs of <italic>Citrus</italic> CC domains in <italic>Citrus</italic> CC2 NBS genes are similar to those motifs of CC domains in japonica rice CC NBS genes [<xref rid="pone.0121893.ref045" ref-type="bibr">45</xref>]. The different structure of motifs in NBS domain and CC domain between <italic>Citrus</italic> CC1 NBS genes and <italic>Citrus</italic> CC2 NBS genes suggested that they are likely from separate evolutionary origins.</p><p>To further confirm the three groups of NBS genes, we identified the NBS genes from <italic>Arabidopsis</italic>, <italic>Populus</italic>, <italic>Oryza sativa</italic> and grape. We used the same criteria to select the NBS domains from NBS genes of these genomes. Finally, we selected 152, 216, 209, 126 NBS domains from <italic>Arabidopsis</italic>, <italic>Populus</italic>, <italic>Oryza sativa</italic> and grape, respectively. Then, together with 442 NBS domains selected from <italic>C</italic>. <italic>clementina</italic>, we constructed a phylogenetic tree of those 1145 NBS domain sequences from five genomes. As shown in <xref rid="pone.0121893.g006" ref-type="fig">Fig. 6</xref>, the un-rooted phylogenetic tree was divided into three main branches. The <italic>Populus</italic>, grape and <italic>Citrus</italic> genomes have significant amount of NBS genes in all three branches. However, the NBS genes of <italic>Oryza sativa</italic> dominated in CC2 branch and most of <italic>Arabidopsis</italic> NBS genes located in CC1 and TIR branches. Our study showed that the NBS genes can be divided into three major groups as the non-TIR NBS genes are separated into two groups. The three groups of NBS genes underwent divergent evolution in different genomes. Further comparison of NBS genes of more genomes may help to understand the evolution of NBS genes and will help elucidate how plants maintain and adapt their defense system against pathogens.</p><fig id="pone.0121893.g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0121893.g006</object-id><label>Fig 6</label><caption><title>The maximum likelihood phylogenetic tree of NBS genes of <italic>C</italic>. <italic>clementina</italic>, <italic>Arabidopsis thaliana</italic>, <italic>Populus trichocarpa</italic>, <italic>Oryza sativa</italic> and <italic>V</italic>. <italic>vinifera</italic> (grape).</title><p>Species are indicated by the color of the outer circle and branches. Nodes with bootstrap support of 100 are indicated with grey circles on the tree.</p></caption><graphic xlink:href="pone.0121893.g006"></graphic></fig></sec><sec id="sec029"><title>Highly Clustering of <italic>Citrus</italic> NBS Genes</title><p>The <italic>Citrus</italic> NBS genes are highly clustered in the genome. 84.9% of NBS genes in <italic>C</italic>. <italic>clementina</italic> and 76.9% of NBS genes in <italic>C</italic>. <italic>sinensis</italic> China were found in clusters. Previous studies showed 76% of rice [<xref rid="pone.0121893.ref048" ref-type="bibr">48</xref>], 64% of <italic>A</italic>. <italic>lyrata</italic> and 71% <italic>A</italic>. <italic>thaliana</italic> [<xref rid="pone.0121893.ref013" ref-type="bibr">13</xref>], 83.2% of grapevine and 67.5% of poplar [<xref rid="pone.0121893.ref009" ref-type="bibr">9</xref>] NBS genes are found in clusters. The percentage of <italic>Citrus</italic> NBS genes in clusters is in the high level comparing to other genomes. The average number of genes per clusters is 4.86 in <italic>C</italic>. <italic>clementina</italic> and 3.97 in <italic>C</italic>. <italic>sinensis</italic> China. These numbers are similar to those in other genomes [<xref rid="pone.0121893.ref009" ref-type="bibr">9</xref>]. Furthermore, most <italic>Citrus</italic> NBS genes in the same cluster belong to the same phylogenetic group, suggesting that tandem duplication is the primary mechanism for the expansion of NBS genes in the <italic>Citrus</italic> genus.</p></sec><sec id="sec030"><title>Molecular Evolution of <italic>Citrus</italic> NBS Genes</title><p>Similar to NBS genes in other genomes, <italic>Citrus</italic> NBS genes are highly dynamic and are shaped by several evolutionary processes leading to several differences amongst NBS genes, including domain presence and mutation constraints and genome organization. Our results revealed multiple molecular evolution events amongst <italic>Citrus</italic> NBS genes including gene duplications, gene conversions, mutation constraint changes, recombination and transposable element insertions. Likely these events support a birth-and-death process leading to the origins of new NBS genes as well as nonfunctionalization and total loss of other NBS genes [<xref rid="pone.0121893.ref012" ref-type="bibr">12</xref>]. We found more than 200 tandem duplications in both <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> China genomes alone. We also found that NBS genes become pseudogenes following original frame-shift mutations leading to mutation accumulations plus disruption of gene constraint leading to loss of function through transposable element insertions.</p><p>Most molecular evolutionary events occurred within NBS gene groups, suggesting that gene birth-death processes have been occurring since divergence from a common ancestral gene copy. Interestingly, the molecular evolution processes occurred differently amongst NBS gene groups and differently within each <italic>Citrus</italic> genome. For example, the number of gene conversion events in <italic>C</italic>. <italic>sinensis</italic> China TIR clusters is almost double than that in <italic>C</italic>. <italic>clementina</italic> TIR clusters, while there are many more conversion events in <italic>C</italic>. <italic>clementina</italic> CC2 clusters than in <italic>C</italic>. <italic>sinensis</italic> China CC2 clusters (<xref rid="pone.0121893.t004" ref-type="table">Table 4</xref>). Numbers of tandem duplications in CC1 and CC2 groups are much greater than that from the TIR group.</p></sec></sec><sec sec-type="conclusions" id="sec031"><title>Conclusions</title><p>Our comparative analyses yield valuable insight into the understanding of the structure, evolution and organization of NBS genes in <italic>Citrus</italic> genomes. <italic>Citrus</italic> NBS genes are structurally highly clustered. The hybrid <italic>C</italic>. <italic>sinensis</italic> genomes have similar types of NBS genes as those progenitor <italic>C</italic>. <italic>clementina</italic> genomes have. Furthermore, our comprehensive analysis also showed that there are three groups of plant NBS genes while non-TIR NBS genes can be divided into two groups. The distributions of three groups of NBS genes among genomes are different implied that the three groups of NBS genes underwent divergent evolution in different genomes.</p></sec><sec sec-type="supplementary-material" id="sec032"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0121893.s001"><label>S1 Fig</label><caption><title>Average intron number of CNL and TNL.</title><p>Cc: <italic>C</italic>. <italic>clementina</italic>, CsCN: <italic>C</italic>. <italic>sinensis</italic> China and CsUSA: <italic>C</italic>. <italic>sinensis</italic> USA.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s002"><label>S2 Fig</label><caption><title>Phylogenetic trees of three groups of <italic>Citrus</italic> NBS gene with P25941 as outgroup.</title><p><bold>A.</bold> CC1 group; <bold>B.</bold> CC2 group; and <bold>C.</bold> TIR group. Nodes with bootstrap support of 100 are indicated with grey circles on the tree.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s003"><label>S3 Fig</label><caption><title>Structure of <italic>Citrus</italic> NBS gene domains.</title><p>The red rectangles indicate NBS domains, the green ellipses indicate CC domains, the left pointing pentagons indicate TIR domains and the left pointing triangles indicate LRR domains.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s004"><label>S4 Fig</label><caption><title>Positively selection sites of Clade_1260.</title><p>The grey box indicate the NBS domain.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s005"><label>S5 Fig</label><caption><title>The exon coverage of NBS-encoding genes in the resequences of <italic>Citrus</italic> genome (3 mandarin and 3 pomelo).</title><p>(PDF)</p></caption><media xlink:href="pone.0121893.s005.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s006"><label>S6 Fig</label><caption><title>Multiple alignments of motifs from CC1 domain (A), CC2 domain (B) and TIR domain (C).</title><p>(PDF)</p></caption><media xlink:href="pone.0121893.s006.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s007"><label>S7 Fig</label><caption><title>Orthologs of <italic>Citrus</italic> NBS-encoding genes of <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> and their genome locations.</title><p>The outer blue circles indicate the nine largest scaffolds of Cc and orange circles indicate nine chromosomes plus a non-localized pseudo-chromosome. The NBS-encoding genes were arranged by their chromosomal positions. The NBS genes in the TIR group are indicated with orange links, CC1 in blue links and CC2 in light blue links. The NBS-encoding orthologs of <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> are indicated by link lines.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s007.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s008"><label>S8 Fig</label><caption><title>Dot plots of syntenic orthologs of the NBS genes between <italic>C</italic>. <italic>clementina</italic> (Cc) and <italic>C</italic>. <italic>sinensis</italic> China (CsCN) produced by MCScanX.</title><p>The numbers in the grids indicate the number of syntenic ortholog pairs in the corresponding scaffolds. The scaffolds of <italic>C</italic>. <italic>clementina</italic> (Cc) are arranged along the x-axis and that of <italic>C</italic>. <italic>sinensis</italic> China (CsCN) are arranged on the y-axis.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s008.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s009"><label>S9 Fig</label><caption><title>Diagram of conserved cluster between <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic>.</title><p>The nodes in orange color indicate clusters of <italic>C</italic>. <italic>clementina</italic>, the blue nodes indicate clusters of <italic>C</italic>. <italic>sinensis</italic>, and the sizes of the nodes are shown relative to cluster size. The numbers on the edges indicate the number of orthologs identified in other species.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s009.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s010"><label>S10 Fig</label><caption><title>Pseudogenization of NBS-encoding genes of <italic>Citrus</italic> due to mutation (A) and retrotransposon insertion (B).</title><p>A: Sequences alignment of Cs1g18610.1, Cs1g18610.1_cc_32 and orange1.1g003367m. The bases indicated the variation bases in each <italic>Citrus</italic> species. The star in red color represented stop-codon gaining mutation and may result in pseudogene of Cs1g18610.1. B: LTR retrotransposon insertion in NBS-encoding genes orange1.1g043039m_cc_116 and orange1.1g043039m_csc_123 from <italic>C</italic>. <italic>clementina</italic> and <italic>C</italic>. <italic>sinensis</italic> respectively, and the multiple sequences alignment of the orthologs. LTR, Long terminal repeats; TSR, Target site repeats.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s010.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s011"><label>S11 Fig</label><caption><title>Validation and phylogenetic analysis of NBS-LRR gene, Cs1g09350.1, in different <italic>Citrus</italic> species.</title><p>A) multiple sequence alignment, with blue bars along exons and thin lines along introns. The vertical lines in each alignment blocks indicate mutations that differ from the reference sequences of <italic>C</italic>. <italic>sinensis</italic>. <bold>B)</bold> neighbor-joining phylogenetic tree of the orthologs of Cs1g09350.1 in different <italic>Citrus</italic> species and bootstrap support displayed on the branches.</p><p>(PDF)</p></caption><media xlink:href="pone.0121893.s011.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s012"><label>S1 Table</label><caption><title>The NBS genes in three <italic>Citrus</italic> genomes.</title><p>(DOCX)</p></caption><media xlink:href="pone.0121893.s012.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s013"><label>S2 Table</label><caption><title>Classification, cluster and clade of <italic>Citrus</italic> NBS genes.</title><p>(XLSX)</p></caption><media xlink:href="pone.0121893.s013.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s014"><label>S3 Table</label><caption><title>LRR domain distribution of <italic>Citrus</italic> NBS gene.</title><p>(XLSX)</p></caption><media xlink:href="pone.0121893.s014.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s015"><label>S4 Table</label><caption><title>Positively selection sites distribution.</title><p>(XLSX)</p></caption><media xlink:href="pone.0121893.s015.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s016"><label>S5 Table</label><caption><title>The exome coverage of NBS genes from the 3 mandarin and 3 pomelo resequenced genomes.</title><p>(XLSX)</p></caption><media xlink:href="pone.0121893.s016.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0121893.s017"><label>S6 Table</label><caption><title><italic>Citrus</italic> NBS gene interrupted with LTR transposons.</title><p>(XLSX)</p></caption><media xlink:href="pone.0121893.s017.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.0121893.ref001"><label>1</label><mixed-citation publication-type="journal"><name><surname>Dangl</surname><given-names>JL</given-names></name>, <name><surname>Jones</surname><given-names>JD</given-names></name>. <article-title>Plant pathogens and integrated defence responses to infection</article-title>. <source>Nature</source>. <year>2001</year>;<volume>411</volume>(<issue>6839</issue>):<fpage>826</fpage>–<lpage>33</lpage>. <pub-id pub-id-type="pmid">11459065</pub-id></mixed-citation></ref><ref id="pone.0121893.ref002"><label>2</label><mixed-citation publication-type="journal"><name><surname>Jia</surname><given-names>Y</given-names></name>, <name><surname>McAdams</surname><given-names>SA</given-names></name>, <name><surname>Bryan</surname><given-names>GT</given-names></name>, <name><surname>Hershey</surname><given-names>HP</given-names></name>, <name><surname>Valent</surname><given-names>B</given-names></name>. <article-title>Direct interaction of resistance gene and avirulence gene products confers rice blast resistance</article-title>. <source>The EMBO Journal</source>. <year>2000</year>;<volume>19</volume>(<issue>15</issue>):<fpage>4004</fpage>–<lpage>14</lpage>. <pub-id pub-id-type="pmid">10921881</pub-id></mixed-citation></ref><ref id="pone.0121893.ref003"><label>3</label><mixed-citation publication-type="journal"><name><surname>Ellis</surname><given-names>JG</given-names></name>, <name><surname>Lawrence</surname><given-names>GJ</given-names></name>, <name><surname>Luck</surname><given-names>JE</given-names></name>, <name><surname>Dodds</surname><given-names>PN</given-names></name>. <article-title>Identification of regions in alleles of the flax rust resistance gene L that determine differences in gene-for-gene specificity</article-title>. <source>The Plant Cell Online</source>. <year>1999</year>;<volume>11</volume>(<issue>3</issue>):<fpage>495</fpage>–<lpage>506</lpage>. <pub-id pub-id-type="pmid">10072407</pub-id></mixed-citation></ref><ref id="pone.0121893.ref004"><label>4</label><mixed-citation publication-type="journal"><name><surname>Ellis</surname><given-names>J</given-names></name>, <name><surname>Dodds</surname><given-names>P</given-names></name>, <name><surname>Pryor</surname><given-names>T</given-names></name>. <article-title>The generation of plant disease resistance gene specificities</article-title>. <source>Trends in plant science</source>. <year>2000</year>;<volume>5</volume>(<issue>9</issue>):<fpage>373</fpage>–<lpage>9</lpage>. <pub-id pub-id-type="pmid">10973092</pub-id></mixed-citation></ref><ref id="pone.0121893.ref005"><label>5</label><mixed-citation publication-type="journal"><name><surname>Mun</surname><given-names>J-H</given-names></name>, <name><surname>Yu</surname><given-names>H-J</given-names></name>, <name><surname>Park</surname><given-names>S</given-names></name>, <name><surname>Park</surname><given-names>B-S</given-names></name>. <article-title>Genome-wide identification of NBS-encoding resistance genes in Brassica rapa</article-title>. <source>Molecular Genetics and Genomics</source>. <year>2009</year>;<volume>282</volume>(<issue>6</issue>):<fpage>617</fpage>–<lpage>31</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s00438-009-0492-0">10.1007/s00438-009-0492-0</ext-link></comment><pub-id pub-id-type="pmid">19838736</pub-id></mixed-citation></ref><ref id="pone.0121893.ref006"><label>6</label><mixed-citation publication-type="journal"><name><surname>Meyers</surname><given-names>BC</given-names></name>, <name><surname>Kozik</surname><given-names>A</given-names></name>, <name><surname>Griego</surname><given-names>A</given-names></name>, <name><surname>Kuang</surname><given-names>H</given-names></name>, <name><surname>Michelmore</surname><given-names>RW</given-names></name>. <article-title>Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis</article-title>. <source>The Plant cell</source>. <year>2003</year>;<volume>15</volume>(<issue>4</issue>):<fpage>809</fpage>–<lpage>34</lpage>. <pub-id pub-id-type="pmid">12671079</pub-id></mixed-citation></ref><ref id="pone.0121893.ref007"><label>7</label><mixed-citation publication-type="journal"><name><surname>Wan</surname><given-names>H</given-names></name>, <name><surname>Yuan</surname><given-names>W</given-names></name>, <name><surname>Bo</surname><given-names>K</given-names></name>, <name><surname>Shen</surname><given-names>J</given-names></name>, <name><surname>Pang</surname><given-names>X</given-names></name>, <name><surname>Chen</surname><given-names>J</given-names></name>. <article-title>Genome-wide analysis of NBS-encoding disease resistance genes in Cucumis sativus and phylogenetic study of NBS-encoding genes in Cucurbitaceae crops</article-title>. <source>BMC genomics</source>. <year>2013</year>;<volume>14</volume>(<issue>1</issue>):<fpage>109</fpage>.<pub-id pub-id-type="pmid">23418910</pub-id></mixed-citation></ref><ref id="pone.0121893.ref008"><label>8</label><mixed-citation publication-type="journal"><name><surname>Sun</surname><given-names>X</given-names></name>, <name><surname>Wang</surname><given-names>G-L</given-names></name>. <article-title>Genome-wide identification, characterization and phylogenetic analysis of the rice LRR-kinases</article-title>. <source>PloS one</source>. <year>2011</year>;<volume>6</volume>(<issue>3</issue>):<fpage>e16079</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0016079">10.1371/journal.pone.0016079</ext-link></comment><pub-id pub-id-type="pmid">21408199</pub-id></mixed-citation></ref><ref id="pone.0121893.ref009"><label>9</label><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>S</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name>, <name><surname>Yue</surname><given-names>J-X</given-names></name>, <name><surname>Tian</surname><given-names>D</given-names></name>, <name><surname>Chen</surname><given-names>J-Q</given-names></name>. <article-title>Recent duplications dominate NBS-encoding gene expansion in two woody species</article-title>. <source>Molecular Genetics and Genomics</source>. <year>2008</year>;<volume>280</volume>(<issue>3</issue>):<fpage>187</fpage>–<lpage>98</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s00438-008-0355-0">10.1007/s00438-008-0355-0</ext-link></comment><pub-id pub-id-type="pmid">18563445</pub-id></mixed-citation></ref><ref id="pone.0121893.ref010"><label>10</label><mixed-citation publication-type="journal"><name><surname>Ameline-Torregrosa</surname><given-names>C</given-names></name>, <name><surname>Wang</surname><given-names>B-B</given-names></name>, <name><surname>O'Bleness</surname><given-names>MS</given-names></name>, <name><surname>Deshpande</surname><given-names>S</given-names></name>, <name><surname>Zhu</surname><given-names>H</given-names></name>, <name><surname>Roe</surname><given-names>B</given-names></name>, <etal>et al</etal><article-title>Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula</article-title>. <source>Plant physiology</source>. <year>2008</year>;<volume>146</volume>(<issue>1</issue>):<fpage>5</fpage>–<lpage>21</lpage>. <pub-id pub-id-type="pmid">17981990</pub-id></mixed-citation></ref><ref id="pone.0121893.ref011"><label>11</label><mixed-citation publication-type="journal"><name><surname>Pan</surname><given-names>Q</given-names></name>, <name><surname>Wendel</surname><given-names>J</given-names></name>, <name><surname>Fluhr</surname><given-names>R</given-names></name>. <article-title>Divergent evolution of plant NBS-LRR resistance gene homologues in dicot and cereal genomes</article-title>. <source>Journal of Molecular Evolution</source>. <year>2000</year>;<volume>50</volume>(<issue>3</issue>):<fpage>203</fpage>–<lpage>13</lpage>. <pub-id pub-id-type="pmid">10754062</pub-id></mixed-citation></ref><ref id="pone.0121893.ref012"><label>12</label><mixed-citation publication-type="journal"><name><surname>Michelmore</surname><given-names>RW</given-names></name>, <name><surname>Meyers</surname><given-names>BC</given-names></name>. <article-title>Clusters of resistance genes in plants evolve by divergent selection and a birth-and-death process</article-title>. <source>Genome Research</source>. <year>1998</year>;<volume>8</volume>(<issue>11</issue>):<fpage>1113</fpage>–<lpage>30</lpage>. <pub-id pub-id-type="pmid">9847076</pub-id></mixed-citation></ref><ref id="pone.0121893.ref013"><label>13</label><mixed-citation publication-type="journal"><name><surname>Guo</surname><given-names>Y-L</given-names></name>, <name><surname>Fitz</surname><given-names>J</given-names></name>, <name><surname>Schneeberger</surname><given-names>K</given-names></name>, <name><surname>Ossowski</surname><given-names>S</given-names></name>, <name><surname>Cao</surname><given-names>J</given-names></name>, <name><surname>Weigel</surname><given-names>D</given-names></name>. <article-title>Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis</article-title>. <source>Plant physiology</source>. <year>2011</year>;<volume>157</volume>(<issue>2</issue>):<fpage>757</fpage>–<lpage>69</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1104/pp.111.181990">10.1104/pp.111.181990</ext-link></comment><pub-id pub-id-type="pmid">21810963</pub-id></mixed-citation></ref><ref id="pone.0121893.ref014"><label>14</label><mixed-citation publication-type="other">Yang S, Li J, Zhang X, Zhang Q, Huang J, Chen J-Q, et al. Rapidly evolving R genes in diverse grass species confer resistance to rice blast disease. Proceedings of the National Academy of Sciences. 2013. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1073/pnas.1318211110">10.1073/pnas.1318211110</ext-link></comment>.</mixed-citation></ref><ref id="pone.0121893.ref015"><label>15</label><mixed-citation publication-type="other">Barrett H, Rhodes A. A numerical taxonomic study of affinity relationships in cultivated Citrus and its close relatives. Systematic Botany. 1976:105–36.</mixed-citation></ref><ref id="pone.0121893.ref016"><label>16</label><mixed-citation publication-type="other">Scora RW. On the history and origin of Citrus. Bulletin of the Torrey Botanical Club. 1975:369–75.</mixed-citation></ref><ref id="pone.0121893.ref017"><label>17</label><mixed-citation publication-type="journal"><name><surname>Nicolosi</surname><given-names>E</given-names></name>, <name><surname>Deng</surname><given-names>Z</given-names></name>, <name><surname>Gentile</surname><given-names>A</given-names></name>, <name><surname>La Malfa</surname><given-names>S</given-names></name>, <name><surname>Continella</surname><given-names>G</given-names></name>, <name><surname>Tribulato</surname><given-names>E</given-names></name>. <article-title>Citrus phylogeny and genetic origin of important species as investigated by molecular markers</article-title>. <source>Theoretical and Applied Genetics</source>. <year>2000</year>;<volume>100</volume>(<issue>8</issue>):<fpage>1155</fpage>–<lpage>66</lpage>.</mixed-citation></ref><ref id="pone.0121893.ref018"><label>18</label><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>X</given-names></name>, <name><surname>Xie</surname><given-names>R</given-names></name>, <name><surname>Lu</surname><given-names>Z</given-names></name>, <name><surname>Zhou</surname><given-names>Z</given-names></name>. <article-title>The origin of cultivated citrus as inferred from internal transcribed spacer and chloroplast DNA sequence and amplified fragment length polymorphism fingerprints</article-title>. <source>Journal of the American Society for Horticultural Science</source>. <year>2010</year>;<volume>135</volume>(<issue>4</issue>):<fpage>341</fpage>–<lpage>50</lpage>.</mixed-citation></ref><ref id="pone.0121893.ref019"><label>19</label><mixed-citation publication-type="journal"><name><surname>Xu</surname><given-names>Q</given-names></name>, <name><surname>Chen</surname><given-names>LL</given-names></name>, <name><surname>Ruan</surname><given-names>X</given-names></name>, <name><surname>Chen</surname><given-names>D</given-names></name>, <name><surname>Zhu</surname><given-names>A</given-names></name>, <name><surname>Chen</surname><given-names>C</given-names></name>, <etal>et al</etal><article-title>The draft genome of sweet orange (Citrus sinensis)</article-title>. <source>Nature genetics</source>. <year>2013</year>;<volume>45</volume>(<issue>1</issue>):<fpage>59</fpage>–<lpage>66</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/ng.2472">10.1038/ng.2472</ext-link></comment><pub-id pub-id-type="pmid">23179022</pub-id></mixed-citation></ref><ref id="pone.0121893.ref020"><label>20</label><mixed-citation publication-type="journal"><name><surname>Xu</surname><given-names>Q</given-names></name>, <name><surname>Biswas</surname><given-names>MK</given-names></name>, <name><surname>Lan</surname><given-names>H</given-names></name>, <name><surname>Zeng</surname><given-names>W</given-names></name>, <name><surname>Liu</surname><given-names>C</given-names></name>, <name><surname>Xu</surname><given-names>J</given-names></name>, <etal>et al</etal><article-title>Phylogenetic and evolutionary analysis of NBS-encoding genes in Rutaceae fruit crops</article-title>. <source>Molecular Genetics and Genomics</source>. <year>2011</year>;<volume>285</volume>(<issue>2</issue>):<fpage>151</fpage>–<lpage>61</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s00438-010-0593-9">10.1007/s00438-010-0593-9</ext-link></comment><pub-id pub-id-type="pmid">21153735</pub-id></mixed-citation></ref><ref id="pone.0121893.ref021"><label>21</label><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>GA</given-names></name>, <name><surname>Prochnik</surname><given-names>S</given-names></name>, <name><surname>Jenkins</surname><given-names>J</given-names></name>, <name><surname>Salse</surname><given-names>J</given-names></name>, <name><surname>Hellsten</surname><given-names>U</given-names></name>, <name><surname>Murat</surname><given-names>F</given-names></name>, <etal>et al</etal><article-title>Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication</article-title>. <source>Nat Biotech</source>. <year>2014</year>;<volume>32</volume>(<issue>7</issue>):<fpage>656</fpage>–<lpage>62</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/nbt.2906">10.1038/nbt.2906</ext-link></comment><pub-id pub-id-type="pmid">24908277</pub-id></mixed-citation></ref><ref id="pone.0121893.ref022"><label>22</label><mixed-citation publication-type="journal"><name><surname>Finn</surname><given-names>RD</given-names></name>, <name><surname>Clements</surname><given-names>J</given-names></name>, <name><surname>Eddy</surname><given-names>SR</given-names></name>. <article-title>HMMER web server: interactive sequence similarity searching</article-title>. <source>Nucleic acids research</source>. <year>2011</year>;<volume>39</volume>(Web Server issue):<fpage>W29</fpage>–<lpage>37</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/nar/gkr367">10.1093/nar/gkr367</ext-link></comment><pub-id pub-id-type="pmid">21593126</pub-id></mixed-citation></ref><ref id="pone.0121893.ref023"><label>23</label><mixed-citation publication-type="journal"><name><surname>Punta</surname><given-names>M</given-names></name>, <name><surname>Coggill</surname><given-names>PC</given-names></name>, <name><surname>Eberhardt</surname><given-names>RY</given-names></name>, <name><surname>Mistry</surname><given-names>J</given-names></name>, <name><surname>Tate</surname><given-names>J</given-names></name>, <name><surname>Boursnell</surname><given-names>C</given-names></name>, <etal>et al</etal><article-title>The Pfam protein families database</article-title>. <source>Nucleic acids research</source>. <year>2012</year>;<volume>40</volume>(Database issue):<fpage>D290</fpage>–<lpage>301</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/nar/gkr1065">10.1093/nar/gkr1065</ext-link></comment><pub-id pub-id-type="pmid">22127870</pub-id></mixed-citation></ref><ref id="pone.0121893.ref024"><label>24</label><mixed-citation publication-type="journal"><name><surname>Boeckmann</surname><given-names>B</given-names></name>, <name><surname>Bairoch</surname><given-names>A</given-names></name>, <name><surname>Apweiler</surname><given-names>R</given-names></name>, <name><surname>Blatter</surname><given-names>MC</given-names></name>, <name><surname>Estreicher</surname><given-names>A</given-names></name>, <name><surname>Gasteiger</surname><given-names>E</given-names></name>, <etal>et al</etal><article-title>The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003</article-title>. <source>Nucleic acids research</source>. <year>2003</year>;<volume>31</volume>(<issue>1</issue>):<fpage>365</fpage>–<lpage>70</lpage>. <pub-id pub-id-type="pmid">12520024</pub-id></mixed-citation></ref><ref id="pone.0121893.ref025"><label>25</label><mixed-citation publication-type="journal"><name><surname>Altschul</surname><given-names>SF</given-names></name>, <name><surname>Madden</surname><given-names>TL</given-names></name>, <name><surname>Schaffer</surname><given-names>AA</given-names></name>, <name><surname>Zhang</surname><given-names>J</given-names></name>, <name><surname>Zhang</surname><given-names>Z</given-names></name>, <name><surname>Miller</surname><given-names>W</given-names></name>, <etal>et al</etal><article-title>Gapped BLAST and PSI-BLAST: a new generation of protein database search programs</article-title>. <source>Nucleic acids research</source>. <year>1997</year>;<volume>25</volume>(<issue>17</issue>):<fpage>3389</fpage>–<lpage>402</lpage>. <pub-id pub-id-type="pmid">9254694</pub-id></mixed-citation></ref><ref id="pone.0121893.ref026"><label>26</label><mixed-citation publication-type="journal"><name><surname>Birney</surname><given-names>E</given-names></name>, <name><surname>Clamp</surname><given-names>M</given-names></name>, <name><surname>Durbin</surname><given-names>R</given-names></name>. <article-title>GeneWise and Genomewise</article-title>. <source>Genome research</source>. <year>2004</year>;<volume>14</volume>(<issue>5</issue>):<fpage>988</fpage>–<lpage>95</lpage>. <pub-id pub-id-type="pmid">15123596</pub-id></mixed-citation></ref><ref id="pone.0121893.ref027"><label>27</label><mixed-citation publication-type="journal"><name><surname>Moreno-Hagelsieb</surname><given-names>G</given-names></name>, <name><surname>Latimer</surname><given-names>K</given-names></name>. <article-title>Choosing BLAST options for better detection of orthologs as reciprocal best hits</article-title>. <source>Bioinformatics</source>. <year>2008</year>;<volume>24</volume>(<issue>3</issue>):<fpage>319</fpage>–<lpage>24</lpage>. <pub-id pub-id-type="pmid">18042555</pub-id></mixed-citation></ref><ref id="pone.0121893.ref028"><label>28</label><mixed-citation publication-type="journal"><name><surname>Suyama</surname><given-names>M</given-names></name>, <name><surname>Torrents</surname><given-names>D</given-names></name>, <name><surname>Bork</surname><given-names>P</given-names></name>. <article-title>PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments</article-title>. <source>Nucleic acids research</source>. <year>2006</year>;<volume>34</volume>(Web Server issue):<fpage>W609</fpage>–<lpage>12</lpage>. <pub-id pub-id-type="pmid">16845082</pub-id></mixed-citation></ref><ref id="pone.0121893.ref029"><label>29</label><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>Z</given-names></name>. <article-title>PAML 4: phylogenetic analysis by maximum likelihood</article-title>. <source>Molecular biology and evolution</source>. <year>2007</year>;<volume>24</volume>(<issue>8</issue>):<fpage>1586</fpage>–<lpage>91</lpage>. <pub-id pub-id-type="pmid">17483113</pub-id></mixed-citation></ref><ref id="pone.0121893.ref030"><label>30</label><mixed-citation publication-type="journal"><name><surname>Katoh</surname><given-names>K</given-names></name>, <name><surname>Standley</surname><given-names>DM</given-names></name>. <article-title>MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability</article-title>. <source>Molecular biology and evolution</source>. <year>2013</year>;<volume>30</volume>(<issue>4</issue>):<fpage>772</fpage>–<lpage>80</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/molbev/mst010">10.1093/molbev/mst010</ext-link></comment><pub-id pub-id-type="pmid">23329690</pub-id></mixed-citation></ref><ref id="pone.0121893.ref031"><label>31</label><mixed-citation publication-type="journal"><name><surname>Price</surname><given-names>MN</given-names></name>, <name><surname>Dehal</surname><given-names>PS</given-names></name>, <name><surname>Arkin</surname><given-names>AP</given-names></name>. <article-title>FastTree: computing large minimum evolution trees with profiles instead of a distance matrix</article-title>. <source>Molecular biology and evolution</source>. <year>2009</year>;<volume>26</volume>(<issue>7</issue>):<fpage>1641</fpage>–<lpage>50</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/molbev/msp077">10.1093/molbev/msp077</ext-link></comment><pub-id pub-id-type="pmid">19377059</pub-id></mixed-citation></ref><ref id="pone.0121893.ref032"><label>32</label><mixed-citation publication-type="book"><name><surname>Retief</surname><given-names>JD</given-names></name>. <chapter-title>Phylogenetic analysis using PHYLIP</chapter-title><source>Bioinformatics methods and protocols</source>: <publisher-name>Springer</publisher-name>; <year>1999</year> p. <fpage>243</fpage>–<lpage>58</lpage>.</mixed-citation></ref><ref id="pone.0121893.ref033"><label>33</label><mixed-citation publication-type="journal"><name><surname>Prosperi</surname><given-names>MC</given-names></name>, <name><surname>Ciccozzi</surname><given-names>M</given-names></name>, <name><surname>Fanti</surname><given-names>I</given-names></name>, <name><surname>Saladini</surname><given-names>F</given-names></name>, <name><surname>Pecorari</surname><given-names>M</given-names></name>, <name><surname>Borghi</surname><given-names>V</given-names></name>, <etal>et al</etal><article-title>A novel methodology for large-scale phylogeny partition</article-title>. <source>Nature communications</source>. <year>2011</year>;<volume>2</volume>:<fpage>321</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/ncomms1325">10.1038/ncomms1325</ext-link></comment><pub-id pub-id-type="pmid">21610724</pub-id></mixed-citation></ref><ref id="pone.0121893.ref034"><label>34</label><mixed-citation publication-type="journal"><name><surname>Delorenzi</surname><given-names>M</given-names></name>, <name><surname>Speed</surname><given-names>T</given-names></name>. <article-title>An HMM model for coiled-coil domains and a comparison with PSSM-based predictions</article-title>. <source>Bioinformatics</source>. <year>2002</year>;<volume>18</volume>(<issue>4</issue>):<fpage>617</fpage>–<lpage>25</lpage>. <pub-id pub-id-type="pmid">12016059</pub-id></mixed-citation></ref><ref id="pone.0121893.ref035"><label>35</label><mixed-citation publication-type="journal"><name><surname>Lupas</surname><given-names>A</given-names></name>. <article-title>Prediction and analysis of coiled-coil structures</article-title>. <source>Methods in enzymology</source>. <year>1996</year>;<volume>266</volume>:<fpage>513</fpage>–<lpage>25</lpage>. <pub-id pub-id-type="pmid">8743703</pub-id></mixed-citation></ref><ref id="pone.0121893.ref036"><label>36</label><mixed-citation publication-type="journal"><name><surname>Bailey</surname><given-names>TL</given-names></name>, <name><surname>Boden</surname><given-names>M</given-names></name>, <name><surname>Buske</surname><given-names>FA</given-names></name>, <name><surname>Frith</surname><given-names>M</given-names></name>, <name><surname>Grant</surname><given-names>CE</given-names></name>, <name><surname>Clementi</surname><given-names>L</given-names></name>, <etal>et al</etal><article-title>MEME SUITE: tools for motif discovery and searching</article-title>. <source>Nucleic acids research</source>. <year>2009</year>;<volume>37</volume>(Web Server issue):<fpage>W202</fpage>–<lpage>8</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/nar/gkp335">10.1093/nar/gkp335</ext-link></comment><pub-id pub-id-type="pmid">19458158</pub-id></mixed-citation></ref><ref id="pone.0121893.ref037"><label>37</label><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>Z</given-names></name>, <name><surname>Carriero</surname><given-names>N</given-names></name>, <name><surname>Zheng</surname><given-names>D</given-names></name>, <name><surname>Karro</surname><given-names>J</given-names></name>, <name><surname>Harrison</surname><given-names>PM</given-names></name>, <name><surname>Gerstein</surname><given-names>M</given-names></name>. <article-title>PseudoPipe: an automated pseudogene identification pipeline</article-title>. <source>Bioinformatics</source>. <year>2006</year>;<volume>22</volume>(<issue>12</issue>):<fpage>1437</fpage>–<lpage>9</lpage>. <pub-id pub-id-type="pmid">16574694</pub-id></mixed-citation></ref><ref id="pone.0121893.ref038"><label>38</label><mixed-citation publication-type="journal"><name><surname>Xu</surname><given-names>Z</given-names></name>, <name><surname>Wang</surname><given-names>H</given-names></name>. <article-title>LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons</article-title>. <source>Nucleic acids research</source>. <year>2007</year>;<volume>35</volume>(Web Server issue):<fpage>W265</fpage>–<lpage>8</lpage>. <pub-id pub-id-type="pmid">17485477</pub-id></mixed-citation></ref><ref id="pone.0121893.ref039"><label>39</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>Y</given-names></name>, <name><surname>Tang</surname><given-names>H</given-names></name>, <name><surname>Debarry</surname><given-names>JD</given-names></name>, <name><surname>Tan</surname><given-names>X</given-names></name>, <name><surname>Li</surname><given-names>J</given-names></name>, <name><surname>Wang</surname><given-names>X</given-names></name>, <etal>et al</etal><article-title>MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity</article-title>. <source>Nucleic acids research</source>. <year>2012</year>;<volume>40</volume>(<issue>7</issue>):<fpage>e49</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/nar/gkr1293">10.1093/nar/gkr1293</ext-link></comment><pub-id pub-id-type="pmid">22217600</pub-id></mixed-citation></ref><ref id="pone.0121893.ref040"><label>40</label><mixed-citation publication-type="journal"><name><surname>Sawyer</surname><given-names>S</given-names></name>. <article-title>Statistical tests for detecting gene conversion</article-title>. <source>Molecular biology and evolution</source>. <year>1989</year>;<volume>6</volume>(<issue>5</issue>):<fpage>526</fpage>–<lpage>38</lpage>. <pub-id pub-id-type="pmid">2677599</pub-id></mixed-citation></ref><ref id="pone.0121893.ref041"><label>41</label><mixed-citation publication-type="other">Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:13033997v2. 2013.</mixed-citation></ref><ref id="pone.0121893.ref042"><label>42</label><mixed-citation publication-type="journal"><name><surname>Quinlan</surname><given-names>AR</given-names></name>, <name><surname>Hall</surname><given-names>IM</given-names></name>. <article-title>BEDTools: a flexible suite of utilities for comparing genomic features</article-title>. <source>Bioinformatics</source>. <year>2010</year>;<volume>26</volume>(<issue>6</issue>):<fpage>841</fpage>–<lpage>2</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bioinformatics/btq033">10.1093/bioinformatics/btq033</ext-link></comment><pub-id pub-id-type="pmid">20110278</pub-id></mixed-citation></ref><ref id="pone.0121893.ref043"><label>43</label><mixed-citation publication-type="journal"><name><surname>Zhou</surname><given-names>L</given-names></name>, <name><surname>Powell</surname><given-names>CA</given-names></name>, <name><surname>Hoffman</surname><given-names>MT</given-names></name>, <name><surname>Li</surname><given-names>W</given-names></name>, <name><surname>Fan</surname><given-names>G</given-names></name>, <name><surname>Liu</surname><given-names>B</given-names></name>, <etal>et al</etal><article-title>Diversity and plasticity of the intracellular plant pathogen and insect symbiont "Candidatus Liberibacter asiaticus" as revealed by hypervariable prophage genes with intragenic tandem repeats</article-title>. <source>Applied and environmental microbiology</source>. <year>2011</year>;<volume>77</volume>(<issue>18</issue>):<fpage>6663</fpage>–<lpage>73</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/AEM.05111-11">10.1128/AEM.05111-11</ext-link></comment><pub-id pub-id-type="pmid">21784907</pub-id></mixed-citation></ref><ref id="pone.0121893.ref044"><label>44</label><mixed-citation publication-type="journal"><name><surname>Kohler</surname><given-names>A</given-names></name>, <name><surname>Rinaldi</surname><given-names>C</given-names></name>, <name><surname>Duplessis</surname><given-names>S</given-names></name>, <name><surname>Baucher</surname><given-names>M</given-names></name>, <name><surname>Geelen</surname><given-names>D</given-names></name>, <name><surname>Duchaussoy</surname><given-names>F</given-names></name>, <etal>et al</etal><article-title>Genome-wide identification of NBS resistance genes in Populus trichocarpa</article-title>. <source>Plant molecular biology</source>. <year>2008</year>;<volume>66</volume>(<issue>6</issue>):<fpage>619</fpage>–<lpage>36</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s11103-008-9293-9">10.1007/s11103-008-9293-9</ext-link></comment><pub-id pub-id-type="pmid">18247136</pub-id></mixed-citation></ref><ref id="pone.0121893.ref045"><label>45</label><mixed-citation publication-type="journal"><name><surname>Zhou</surname><given-names>T</given-names></name>, <name><surname>Wang</surname><given-names>Y</given-names></name>, <name><surname>Chen</surname><given-names>JQ</given-names></name>, <name><surname>Araki</surname><given-names>H</given-names></name>, <name><surname>Jing</surname><given-names>Z</given-names></name>, <name><surname>Jiang</surname><given-names>K</given-names></name>, <etal>et al</etal><article-title>Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes</article-title>. <source>Molecular genetics and genomics: MGG</source>. <year>2004</year>;<volume>271</volume>(<issue>4</issue>):<fpage>402</fpage>–<lpage>15</lpage>. <pub-id pub-id-type="pmid">15014983</pub-id></mixed-citation></ref><ref id="pone.0121893.ref046"><label>46</label><mixed-citation publication-type="journal"><name><surname>Tan</surname><given-names>S</given-names></name>, <name><surname>Wu</surname><given-names>S</given-names></name>. <article-title>Genome Wide Analysis of Nucleotide-Binding Site Disease Resistance Genes in Brachypodium distachyon</article-title>. <source>Comparative and functional genomics</source>. <year>2012</year>;<volume>2012</volume>:<fpage>418208</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1155/2012/418208">10.1155/2012/418208</ext-link></comment><pub-id pub-id-type="pmid">22693425</pub-id></mixed-citation></ref><ref id="pone.0121893.ref047"><label>47</label><mixed-citation publication-type="journal"><name><surname>Punta</surname><given-names>M</given-names></name>, <name><surname>Coggill</surname><given-names>PC</given-names></name>, <name><surname>Eberhardt</surname><given-names>RY</given-names></name>, <name><surname>Mistry</surname><given-names>J</given-names></name>, <name><surname>Tate</surname><given-names>J</given-names></name>, <name><surname>Boursnell</surname><given-names>C</given-names></name>, <etal>et al</etal><article-title>The Pfam protein families database</article-title>. <source>Nucleic Acids Research</source>. <year>2012</year>;<volume>40</volume>(<issue>D1</issue>):<fpage>D290</fpage>–<lpage>D301</lpage>.<pub-id pub-id-type="pmid">22127870</pub-id></mixed-citation></ref><ref id="pone.0121893.ref048"><label>48</label><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>S</given-names></name>, <name><surname>Feng</surname><given-names>Z</given-names></name>, <name><surname>Zhang</surname><given-names>X</given-names></name>, <name><surname>Jiang</surname><given-names>K</given-names></name>, <name><surname>Jin</surname><given-names>X</given-names></name>, <name><surname>Hang</surname><given-names>Y</given-names></name>, <etal>et al</etal><article-title>Genome-wide investigation on the genetic variations of rice disease resistance genes</article-title>. <source>Plant molecular biology</source>. <year>2006</year>;<volume>62</volume>(<issue>1–2</issue>):<fpage>181</fpage>–<lpage>93</lpage>.<pub-id pub-id-type="pmid">16915523</pub-id></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 000136 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000136 | 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. |  |