Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The Evolution of a High Copy Gene Array in <italic>Arabidopsis</italic></title><author><name sortKey="Kane, Joshua" sort="Kane, Joshua" uniqKey="Kane J" first="Joshua" last="Kane">Joshua Kane</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author><author><name sortKey="Freeling, Michael" sort="Freeling, Michael" uniqKey="Freeling M" first="Michael" last="Freeling">Michael Freeling</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author><author><name sortKey="Lyons, Eric" sort="Lyons, Eric" uniqKey="Lyons E" first="Eric" last="Lyons">Eric Lyons</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">20495794</idno><idno type="pmc">2886086</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2886086</idno><idno type="RBID">PMC:2886086</idno><idno type="doi">10.1007/s00239-010-9350-2</idno><date when="2010">2010</date><idno type="wicri:Area/Pmc/Corpus">000185</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The Evolution of a High Copy Gene Array in <italic>Arabidopsis</italic></title><author><name sortKey="Kane, Joshua" sort="Kane, Joshua" uniqKey="Kane J" first="Joshua" last="Kane">Joshua Kane</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author><author><name sortKey="Freeling, Michael" sort="Freeling, Michael" uniqKey="Freeling M" first="Michael" last="Freeling">Michael Freeling</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author><author><name sortKey="Lyons, Eric" sort="Lyons, Eric" uniqKey="Lyons E" first="Eric" last="Lyons">Eric Lyons</name><affiliation><nlm:aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</nlm:aff></affiliation></author></analytic><series><title level="j">Journal of Molecular Evolution</title><idno type="ISSN">0022-2844</idno><idno type="eISSN">1432-1432</idno><imprint><date when="2010">2010</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Local gene duplication is a prominent mechanism of gene copy number expansion. Elucidating the mechanisms by which local duplicates arise is necessary in understanding the evolution of genomes and their host organisms. Chromosome one of <italic>Arabidopsis thaliana</italic> contains an 81-gene array subdivided into 27 triplet units (t-units), with each t-unit containing three pre-transfer RNA genes. We utilized phylogenetic tree reconstructions and comparative genomics to order the events leading to the array’s formation, and propose a model using unequal crossing-over as the primary mechanism of array formation. The model is supported by additional phylogenetic information from intergenic spacer sequences separating each t-unit, comparative analysis to an orthologous array of 12 t-units in the sister taxa <italic>Arabidopsis lyrata</italic>, and additional modeling using a stochastic simulation of orthologous array divergence. Lastly, comparative phylogenetic analysis demonstrates that the two orthologous t-unit arrays undergo concerted evolution within each taxa and are likely fluctuating in copy number under neutral evolutionary drift. These findings hold larger implications for future research concerning gene and genome evolution.</p><sec><title>Electronic supplementary material</title><p>The online version of this article (doi:10.1007/s00239-010-9350-2) contains supplementary material, which is available to authorized users.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Achaz, G" uniqKey="Achaz G">G Achaz</name></author><author><name sortKey="Coissac, E" uniqKey="Coissac E">E Coissac</name></author><author><name sortKey="Viari, A" uniqKey="Viari A">A Viari</name></author><author><name sortKey="Netter, P" uniqKey="Netter P">P Netter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name></author><author><name sortKey="Gish, W" uniqKey="Gish W">W Gish</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author><author><name sortKey="Myers, Ew" uniqKey="Myers E">EW Myers</name></author><author><name sortKey="Lipman, Dj" uniqKey="Lipman D">DJ Lipman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baumgarten, A" uniqKey="Baumgarten A">A Baumgarten</name></author><author><name sortKey="Cannon, S" uniqKey="Cannon S">S Cannon</name></author><author><name sortKey="Spangler, R" uniqKey="Spangler R">R Spangler</name></author><author><name sortKey="May, G" uniqKey="May G">G May</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Beier, D" uniqKey="Beier D">D Beier</name></author><author><name sortKey="Beier, H" uniqKey="Beier H">H Beier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Beier, D" uniqKey="Beier D">D Beier</name></author><author><name sortKey="Stange, N" uniqKey="Stange N">N Stange</name></author><author><name sortKey="Gross, Hj" uniqKey="Gross H">HJ Gross</name></author><author><name sortKey="Beier, H" uniqKey="Beier H">H Beier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bowers, Je" uniqKey="Bowers J">JE Bowers</name></author><author><name sortKey="Chapman, Ba" uniqKey="Chapman B">BA Chapman</name></author><author><name sortKey="Rong, J" uniqKey="Rong J">J Rong</name></author><author><name sortKey="Paterson, Ah" uniqKey="Paterson A">AH Paterson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brown, Dd" uniqKey="Brown D">DD Brown</name></author><author><name sortKey="Wensink, Pc" uniqKey="Wensink P">PC Wensink</name></author><author><name sortKey="Jordan, E" uniqKey="Jordan E">E Jordan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Eichler, Ee" uniqKey="Eichler E">EE Eichler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Force, A" uniqKey="Force A">A Force</name></author><author><name sortKey="Lynch, M" uniqKey="Lynch M">M Lynch</name></author><author><name sortKey="Pickett, Fb" uniqKey="Pickett F">FB Pickett</name></author><author><name sortKey="Amores, A" uniqKey="Amores A">A Amores</name></author><author><name sortKey="Yan, Yl" uniqKey="Yan Y">YL Yan</name></author><author><name sortKey="Postlethwait, J" uniqKey="Postlethwait J">J Postlethwait</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Freeling, M" uniqKey="Freeling M">M Freeling</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ganley, Ard" uniqKey="Ganley A">ARD Ganley</name></author><author><name sortKey="Kobayashi, T" uniqKey="Kobayashi T">T Kobayashi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gao, L Z" uniqKey="Gao L">L-Z Gao</name></author><author><name sortKey="Innan, H" uniqKey="Innan H">H Innan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gaut, Bs" uniqKey="Gaut B">BS Gaut</name></author><author><name sortKey="Wright, Si" uniqKey="Wright S">SI Wright</name></author><author><name sortKey="Rizzon, C" uniqKey="Rizzon C">C Rizzon</name></author><author><name sortKey="Dvorak, J" uniqKey="Dvorak J">J Dvorak</name></author><author><name sortKey="Anderson, Lk" uniqKey="Anderson L">LK Anderson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goodman, Hm" uniqKey="Goodman H">HM Goodman</name></author><author><name sortKey="Olson, Mv" uniqKey="Olson M">MV Olson</name></author><author><name sortKey="Hall, Bd" uniqKey="Hall B">BD Hall</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hofstetter, H" uniqKey="Hofstetter H">H Hofstetter</name></author><author><name sortKey="Kressman, A" uniqKey="Kressman A">A Kressman</name></author><author><name sortKey="Birnstiel, Ml" uniqKey="Birnstiel M">ML Birnstiel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Horiuchi, T" uniqKey="Horiuchi T">T Horiuchi</name></author><author><name sortKey="Tomizawa, Ji" uniqKey="Tomizawa J">JI Tomizawa</name></author><author><name sortKey="Novick, A" uniqKey="Novick A">A Novick</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Jackson, Ap" uniqKey="Jackson A">AP Jackson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jaillon, O" uniqKey="Jaillon O">O Jaillon</name></author><author><name sortKey="Aury, J M" uniqKey="Aury J">J-M Aury</name></author><author><name sortKey="Noel, B" uniqKey="Noel B">B Noel</name></author><author><name sortKey="Policriti, A" uniqKey="Policriti A">A Policriti</name></author><author><name sortKey="Clepet, C" uniqKey="Clepet C">C Clepet</name></author><author><name sortKey="Casagrande, A" uniqKey="Casagrande A">A Casagrande</name></author><author><name sortKey="Choisne, N" uniqKey="Choisne N">N Choisne</name></author><author><name sortKey="Aubourg, S" uniqKey="Aubourg S">S Aubourg</name></author><author><name sortKey="Vitulo, N" uniqKey="Vitulo N">N Vitulo</name></author><author><name sortKey="Jubin, C" uniqKey="Jubin C">C Jubin</name></author><author><name sortKey="Vezzi, A" uniqKey="Vezzi A">A Vezzi</name></author><author><name sortKey="Legeai, F" uniqKey="Legeai F">F Legeai</name></author><author><name sortKey="Hugueney, P" uniqKey="Hugueney P">P Hugueney</name></author><author><name sortKey="Dasilva, C" uniqKey="Dasilva C">C Dasilva</name></author><author><name sortKey="Horner, D" uniqKey="Horner D">D Horner</name></author><author><name sortKey="Mica, E" uniqKey="Mica E">E Mica</name></author><author><name sortKey="Jublot, D" uniqKey="Jublot D">D Jublot</name></author><author><name sortKey="Poulain, J" uniqKey="Poulain J">J Poulain</name></author><author><name sortKey="Bruyere, C" uniqKey="Bruyere C">C Bruyère</name></author><author><name sortKey="Billault, A" uniqKey="Billault A">A Billault</name></author><author><name sortKey="Segurens, B" uniqKey="Segurens B">B Segurens</name></author><author><name sortKey="Gouyvenoux, M" uniqKey="Gouyvenoux M">M Gouyvenoux</name></author><author><name sortKey="Ugarte, E" uniqKey="Ugarte E">E Ugarte</name></author><author><name sortKey="Cattonaro, F" uniqKey="Cattonaro F">F Cattonaro</name></author><author><name sortKey="Anthouard, V" uniqKey="Anthouard V">V Anthouard</name></author><author><name sortKey="Vico, V" uniqKey="Vico V">V Vico</name></author><author><name sortKey="Fabbro, Cd" uniqKey="Fabbro C">CD Fabbro</name></author><author><name sortKey="Alaux, M" uniqKey="Alaux M">M Alaux</name></author><author><name sortKey="Gaspero, Gd" uniqKey="Gaspero G">GD Gaspero</name></author><author><name sortKey="Dumas, V" uniqKey="Dumas V">V Dumas</name></author><author><name sortKey="Felice, N" uniqKey="Felice N">N Felice</name></author><author><name sortKey="Paillard, S" uniqKey="Paillard S">S Paillard</name></author><author><name sortKey="Juman, I" uniqKey="Juman I">I Juman</name></author><author><name sortKey="Moroldo, M" uniqKey="Moroldo M">M Moroldo</name></author><author><name sortKey="Scalabrin, S" uniqKey="Scalabrin S">S Scalabrin</name></author><author><name sortKey="Canaguier, A" uniqKey="Canaguier A">A Canaguier</name></author><author><name sortKey="Clainche, Il" uniqKey="Clainche I">IL Clainche</name></author><author><name sortKey="Malacrida, G" uniqKey="Malacrida G">G Malacrida</name></author><author><name sortKey="Durand, E" uniqKey="Durand E">E Durand</name></author><author><name sortKey="Pesole, G" uniqKey="Pesole G">G Pesole</name></author><author><name sortKey="Laucou, V" uniqKey="Laucou V">V Laucou</name></author><author><name sortKey="Chatelet, P" uniqKey="Chatelet P">P Chatelet</name></author><author><name sortKey="Merdinoglu, D" uniqKey="Merdinoglu D">D Merdinoglu</name></author><author><name sortKey="Delledonne, M" uniqKey="Delledonne M">M Delledonne</name></author><author><name sortKey="Pezzotti, M" uniqKey="Pezzotti M">M Pezzotti</name></author><author><name sortKey="Lecharny, A" uniqKey="Lecharny A">A Lecharny</name></author><author><name sortKey="Scarpelli, C" uniqKey="Scarpelli C">C Scarpelli</name></author><author><name sortKey="Artiguenave, F" uniqKey="Artiguenave F">F Artiguenave</name></author><author><name sortKey="Pe, Me" uniqKey="Pe M">ME Pè</name></author><author><name sortKey="Valle, G" uniqKey="Valle G">G Valle</name></author><author><name sortKey="Morgante, M" uniqKey="Morgante M">M Morgante</name></author><author><name sortKey="Caboche, M" uniqKey="Caboche M">M Caboche</name></author><author><name sortKey="Adam Blondon, A F" uniqKey="Adam Blondon A">A-F Adam-Blondon</name></author><author><name sortKey="Weissenbach, J" uniqKey="Weissenbach J">J Weissenbach</name></author><author><name sortKey="Quetier, F" uniqKey="Quetier F">F Quétier</name></author><author><name sortKey="Wincker, P" uniqKey="Wincker P">P Wincker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jelesko, Jg" uniqKey="Jelesko J">JG Jelesko</name></author><author><name sortKey="Carter, K" uniqKey="Carter K">K Carter</name></author><author><name sortKey="Thompson, W" uniqKey="Thompson W">W Thompson</name></author><author><name sortKey="Kinoshita, Y" uniqKey="Kinoshita Y">Y Kinoshita</name></author><author><name sortKey="Gruissem, W" uniqKey="Gruissem W">W Gruissem</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koch, Ma" uniqKey="Koch M">MA Koch</name></author><author><name sortKey="Kiefer, M" uniqKey="Kiefer M">M Kiefer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lewis, D" uniqKey="Lewis D">D Lewis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lowe, Tm" uniqKey="Lowe T">TM Lowe</name></author><author><name sortKey="Eddy, Sr" uniqKey="Eddy S">SR Eddy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lynch, M" uniqKey="Lynch M">M Lynch</name></author><author><name sortKey="Conery, Js" uniqKey="Conery J">JS Conery</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lyons, E" uniqKey="Lyons E">E Lyons</name></author><author><name sortKey="Freeling, M" uniqKey="Freeling M">M Freeling</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lyons, E" uniqKey="Lyons E">E Lyons</name></author><author><name sortKey="Pedersen, B" uniqKey="Pedersen B">B Pedersen</name></author><author><name sortKey="Kane, J" uniqKey="Kane J">J Kane</name></author><author><name sortKey="Alam, M" uniqKey="Alam M">M Alam</name></author><author><name sortKey="Ming, R" uniqKey="Ming R">R Ming</name></author><author><name sortKey="Tang, H" uniqKey="Tang H">H Tang</name></author><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name></author><author><name sortKey="Bowers, J" uniqKey="Bowers J">J Bowers</name></author><author><name sortKey="Paterson, A" uniqKey="Paterson A">A Paterson</name></author><author><name sortKey="Lisch, D" uniqKey="Lisch D">D Lisch</name></author><author><name sortKey="Freeling, M" uniqKey="Freeling M">M Freeling</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Nei, M" uniqKey="Nei M">M Nei</name></author><author><name sortKey="Rooney, Ap" uniqKey="Rooney A">AP Rooney</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Ohno, S" uniqKey="Ohno S">S Ohno</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ohta, T" uniqKey="Ohta T">T Ohta</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Posada, D" uniqKey="Posada D">D Posada</name></author><author><name sortKey="Crandall, Ka" uniqKey="Crandall K">KA Crandall</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Richard, G F" uniqKey="Richard G">G-F Richard</name></author><author><name sortKey="Dujon, B" uniqKey="Dujon B">B Dujon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ronquist, F" uniqKey="Ronquist F">F Ronquist</name></author><author><name sortKey="Huelsenbeck, Jp" uniqKey="Huelsenbeck J">JP Huelsenbeck</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schuermann, D" uniqKey="Schuermann D">D Schuermann</name></author><author><name sortKey="Molinier, J" uniqKey="Molinier J">J Molinier</name></author><author><name sortKey="Fritsch, O" uniqKey="Fritsch O">O Fritsch</name></author><author><name sortKey="Hohn, B" uniqKey="Hohn B">B Hohn</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schwartz, S" uniqKey="Schwartz S">S Schwartz</name></author><author><name sortKey="Kent, Wj" uniqKey="Kent W">WJ Kent</name></author><author><name sortKey="Smit, A" uniqKey="Smit A">A Smit</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Baertsch, R" uniqKey="Baertsch R">R Baertsch</name></author><author><name sortKey="Hardison, Rc" uniqKey="Hardison R">RC Hardison</name></author><author><name sortKey="Haussler, D" uniqKey="Haussler D">D Haussler</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Semon, M" uniqKey="Semon M">M Sémon</name></author><author><name sortKey="Wolfe, Kh" uniqKey="Wolfe K">KH Wolfe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, Gp" uniqKey="Smith G">GP Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stange, N" uniqKey="Stange N">N Stange</name></author><author><name sortKey="Beier, D" uniqKey="Beier D">D Beier</name></author><author><name sortKey="Beier, H" uniqKey="Beier H">H Beier</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Tartof, Kd" uniqKey="Tartof K">KD Tartof</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Thompson, Jd" uniqKey="Thompson J">JD Thompson</name></author><author><name sortKey="Higgins, Dg" uniqKey="Higgins D">DG Higgins</name></author><author><name sortKey="Gibson, Tj" uniqKey="Gibson T">TJ Gibson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Yw" uniqKey="Yang Y">YW Yang</name></author><author><name sortKey="Lai, Kn" uniqKey="Lai K">KN Lai</name></author><author><name sortKey="Tai, Py" uniqKey="Tai P">PY Tai</name></author><author><name sortKey="Li, Wh" uniqKey="Li W">WH Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, J" uniqKey="Zhang J">J Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zimmer, Ea" uniqKey="Zimmer E">EA Zimmer</name></author><author><name sortKey="Martin, Sl" uniqKey="Martin S">SL Martin</name></author><author><name sortKey="Beverley, Sm" uniqKey="Beverley S">SM Beverley</name></author><author><name sortKey="Kan, Yw" uniqKey="Kan Y">YW Kan</name></author><author><name sortKey="Wilson, Ac" uniqKey="Wilson A">AC Wilson</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">J Mol Evol</journal-id><journal-title-group><journal-title>Journal of Molecular Evolution</journal-title></journal-title-group><issn pub-type="ppub">0022-2844</issn><issn pub-type="epub">1432-1432</issn><publisher><publisher-name>Springer-Verlag</publisher-name><publisher-loc>New York</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">20495794</article-id><article-id pub-id-type="pmc">2886086</article-id><article-id pub-id-type="publisher-id">9350</article-id><article-id pub-id-type="doi">10.1007/s00239-010-9350-2</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>The Evolution of a High Copy Gene Array in <italic>Arabidopsis</italic></article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Kane</surname><given-names>Joshua</given-names></name><xref ref-type="aff" rid="Aff1"></xref></contrib><contrib contrib-type="author"><name><surname>Freeling</surname><given-names>Michael</given-names></name><xref ref-type="aff" rid="Aff1"></xref></contrib><contrib contrib-type="author" corresp="yes"><name><surname>Lyons</surname><given-names>Eric</given-names></name><address><phone>+1-510-6428058</phone><fax>+1-510-6424995</fax><email>elyons@berkeley.edu</email></address><xref ref-type="aff" rid="Aff1"></xref></contrib><aff id="Aff1">Department of Plant and Microbial Biology, University of California at Berkeley, 311 Koshland Hall, Berkeley, CA 94720 USA</aff></contrib-group><pub-date pub-type="epub"><day>22</day><month>5</month><year>2010</year></pub-date><pub-date pub-type="pmc-release"><day>22</day><month>5</month><year>2010</year></pub-date><pub-date pub-type="ppub"><month>6</month><year>2010</year></pub-date><volume>70</volume><issue>6</issue><fpage>531</fpage><lpage>544</lpage><history><date date-type="received"><day>8</day><month>12</month><year>2009</year></date><date date-type="accepted"><day>3</day><month>5</month><year>2010</year></date></history><permissions><copyright-statement>© The Author(s) 2010</copyright-statement></permissions><abstract><p>Local gene duplication is a prominent mechanism of gene copy number expansion. Elucidating the mechanisms by which local duplicates arise is necessary in understanding the evolution of genomes and their host organisms. Chromosome one of <italic>Arabidopsis thaliana</italic> contains an 81-gene array subdivided into 27 triplet units (t-units), with each t-unit containing three pre-transfer RNA genes. We utilized phylogenetic tree reconstructions and comparative genomics to order the events leading to the array’s formation, and propose a model using unequal crossing-over as the primary mechanism of array formation. The model is supported by additional phylogenetic information from intergenic spacer sequences separating each t-unit, comparative analysis to an orthologous array of 12 t-units in the sister taxa <italic>Arabidopsis lyrata</italic>, and additional modeling using a stochastic simulation of orthologous array divergence. Lastly, comparative phylogenetic analysis demonstrates that the two orthologous t-unit arrays undergo concerted evolution within each taxa and are likely fluctuating in copy number under neutral evolutionary drift. These findings hold larger implications for future research concerning gene and genome evolution.</p><sec><title>Electronic supplementary material</title><p>The online version of this article (doi:10.1007/s00239-010-9350-2) contains supplementary material, which is available to authorized users.</p></sec></abstract><kwd-group><title>Keywords</title><kwd>Comparative genomics</kwd><kwd>Concerted evolution</kwd><kwd>Phylogenetics</kwd><kwd>Tandem duplication</kwd><kwd>Local duplication</kwd><kwd>Gene array</kwd><kwd>tRNA</kwd><kwd>Arabidopsis</kwd><kwd>Synteny</kwd><kwd>Copy number variation</kwd></kwd-group><custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name><meta-value>© Springer Science+Business Media, LLC 2010</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="Sec1"><title>Introduction</title><p>Local gene duplication generates tandemly arrayed genes and is a ubiquitous occurrence during genome evolution. Local duplicates are presumed to contribute to genetic diversity through relaxation of selective pressure on one or both duplicates through redundancy of function. This relaxed selective pressure allows for gene evolution such as subfunctionalization (Force et al. <xref ref-type="bibr" rid="CR9">1999</xref>; Lynch and Conery <xref ref-type="bibr" rid="CR24">2000</xref>), where both the genes specialize in a subset of their original function, and neofunctionalization (Lewis <xref ref-type="bibr" rid="CR22">1953</xref>; Ohno <xref ref-type="bibr" rid="CR32">1970</xref>; Ohta <xref ref-type="bibr" rid="CR33">1988</xref>), where one of the copies evolves novel cellular functions while its counterpart retains ancestral functions. For duplication events that present a dominant negative dosage phenotype, subsequent selection may result in reversion back to the singlet condition, returning gene product dosage to a stoichiometric balance (Freeling <xref ref-type="bibr" rid="CR10">2009</xref>). Unequal crossing-over (UCO) between two duplex loci naturally forms a triplet and a singlet, so that there are no mechanical impediments to returning to the singlet state.</p><p>In contrast, some genomic loci can undergo copy number amplification under positive selection for increase in product dosage. Positive selection for increased dosage has been suggested to explain large tandem arrays of rRNA and histone genes (Hurles <xref ref-type="bibr" rid="CR17">2004</xref>) and is observed experimentally in bacteria to cope with nutrient availability (Horiuchi et al. <xref ref-type="bibr" rid="CR16">1962</xref>). A subset of genes that frequently occur in local tandem arrays is observed to have a high degree of sequence and functional diversity; some classes of plant disease resistance genes (e.g., NB LRR genes) are well-studied examples of this phenomenon. A single cluster of tandemly repeated NB-LRR genes may encode a diverse number of different pathogen specificities (Baumgarten et al. <xref ref-type="bibr" rid="CR3">2003</xref>). These arrays of high sequence diversity stand in contrast to the aforementioned rRNA arrays, for which sequence and functional homogenization (concerted evolution) is observed (Brown et al. <xref ref-type="bibr" rid="CR7">1972</xref>; Nei and Rooney <xref ref-type="bibr" rid="CR30">2005</xref>). Please see Sémon and Wolfe (<xref ref-type="bibr" rid="CR40">2007</xref>) for an excellent review of the evolution and fates of gene duplicates.</p><p>Our study focuses on an unusual pre-tRNA gene array within the <italic>Arabidopsis</italic> lineage. In <italic>A. thaliana</italic> (<italic>At</italic>), this array consists of 81 pre-tRNA genes arranged into 27 repeating triplet gene units, which we call t-units. This pre-tRNA gene array was discovered in <italic>At</italic> by Beier et al. (<xref ref-type="bibr" rid="CR5">1991</xref>) and shown to contain at least 20 copies of t-units. Further study provided data that a number of these pre-tRNA genes produce transcripts (Beier and Beier <xref ref-type="bibr" rid="CR4">1992</xref>; Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>). We show that <italic>Arabidopsis lyrata</italic> (<italic>Al</italic>), a sister taxon of <italic>At</italic>, contains an orthologous array of pre-tRNA genes at the syntenic locus to <italic>At</italic>. Using the complete genome sequences of both the <italic>At</italic> and <italic>Al</italic> genomes, our research has expanded upon the study of Beier et al. using a high-resolution comparative analysis of the t-unit array within two taxa of the <italic>Arabidopsis</italic> lineage to understand the evolution of these t-unit pre-tRNA gene arrays.</p><p>The evolutionary fate of duplicated genes has been the focus of a great deal of research (Zhang <xref ref-type="bibr" rid="CR48">2003</xref>; Sémon and Wolfe <xref ref-type="bibr" rid="CR40">2007</xref>; Freeling <xref ref-type="bibr" rid="CR10">2009</xref>); in contrast, detailed case studies on tandemly arrayed genes containing many repeat units are relatively few in number. While characterizations of the evolution of tandemly arrayed genes using comparative methods have been made (Jackson <xref ref-type="bibr" rid="CR18">2007</xref>), more studies are required in a wider range of taxa to gain an understanding of the general mechanisms that govern the evolution of tandem arrays in genomes. The <italic>Arabidopsis</italic> t-unit pre-tRNA gene array provides an excellent model system for gaining insight into the mechanisms of tandemly arrayed gene evolution due to its recent formation, its rapid expansion, and its conspicuous pattern of multi-gene duplication, the latter hinting at a higher order complexity to the local duplication process. Through a series of phylogenetic and comparative analyses, we demonstrate that complete t-units are almost without exception the basal units of the arrays’ expansion and evolve in concert within each taxon. In addition, we propose a generalized model of t-unit expansion within the array using unequal crossing-over (UCO) as the mechanism of copy number variation and sequence “homogenization,” and designed a custom computer program to stochastically simulate divergence between two orthologous arrays. These simulations demonstrate that relatively few crossover events within an array are necessary to cause lineage-specific divergence of orthologous arrays, as is observed in the <italic>At</italic> and <italic>Al</italic> t-unit arrays, leaving us to conclude that UCO is the primary mechanism of evolution within these tandem gene arrays and sufficient to explain their concerted evolution.</p></sec><sec id="Sec2" sec-type="materials|methods"><title>Materials and Methods</title><p>Many of the sequence comparisons used in this article were performed using the CoGe comparative genomics platform (Lyons and Freeling <xref ref-type="bibr" rid="CR25">2008</xref>; Lyons et al. <xref ref-type="bibr" rid="CR26">2008</xref>).</p><sec id="Sec3"><title>Nomenclature and Triplet Designation</title><p>A nomenclature was devised to denote both element type and location to maintain specificity in language for these highly repetitive genomic structures (Fig. <xref rid="Fig1" ref-type="fig">1</xref>). The pre-tRNA genes are denoted by the amino acid they transfer. A triplet unit, or t-unit, is defined as the genomic sequence from the start of a pre-tRNA<sup>Ser</sup> gene to the end of the second proximal 3′ pre-tRNA<sup>Tyr</sup> gene, with the t-unit oriented in the 5′ to 3′ direction. Each t-unit contains five elements in the following order: a pre-tRNA<sup>Ser</sup> gene, an intergenic region (IGR) B, a pre-tRNA<sup>Tyr</sup>A gene, an IGR C, and a pre-tRNA<sup>Tyr</sup>B gene. Exceptions to this rule are an incomplete t-unit on the 5′ end of the <italic>Al</italic> array, referred to as t-unit 0, and t-unit 1 in <italic>Al</italic> which contains an extra pre-tRNA<sup>Tyr</sup>B gene, referred to as pre-tRNA<sup>Tyr</sup>C, and an extra IGR C sequence, referred to as IGR D1; the <italic>At</italic> array contains no exceptions to the t-unit rule. Preceding each t-unit is an additional intergenic region, IGR A, defined as the sequence from the end of the prior (5′) t-unit to the start of the current t-unit. Suffix letters are used denote relative position within the t-unit—IGR B precedes IGR C, and pre-tRNA<sup>Tyr</sup>A precedes pre-tRNA<sup>Tyr</sup>B. Each element and t-unit is also given a number to indicate its absolute position from the 5′ end of the array. For example, t-unit 23 is the 23rd t-unit in the array, and IGR B12 is the intergenic region between pre-tRNA<sup>Ser</sup> and pre-tRNA<sup>Tyr</sup>A in t-unit 12. The number suffix of an IGR A sequence indicates its position is immediately 5′ of the t-unit of that number. The exception to this rule is IGR A28 in <italic>At</italic>, which is the IGR A sequence immediately 3′ of the last t-unit, t-unit 27. In addition, when t-units and array elements are referred to by the species from which they originate, <italic>At</italic> and <italic>Al</italic> refer to <italic>A. thaliana</italic> and <italic>A. lyrata</italic>, respectively.<fig id="Fig1"><label>Fig. 1</label><caption><p>A visualization of the triplet array in <italic>A. thaliana</italic> (<italic>At</italic>) and <italic>A. lyrata</italic> (<italic>Al</italic>). (<bold>a</bold>) The triplet array in <italic>At</italic>. The pre-tRNA genes are indicated by gray cylinders. The green background indicates G + C content; each G + C basepair is indicated by a single vertical green line. The entire array is 39.7 kb long, from the start of the first pre-tRNA<sup>Ser</sup> on the 5′ end to the final pre-tRNA<sup>Tyr</sup>B on the 3′ end (Chr. 1 21,272,452..21,312,137). (<bold>b</bold>) Two triplet units (t-units) in the array are shown, starting with the twelfth t-unit in the array on the left, and the thirteenth t-unit on the right. Each triplet indicates the two methods of nomenclature used throughout the article. Both letters and numbers are used to distinguish element position in different ways, with all nomenclature increasing from 5′ to 3′. Letters, such as IGR C, indicate relative position of an intergenic region within the t-unit, whereas numbers, such as IGR C13, indicate the element’s absolute position within the array. Thus, IGR C13 indicates it is the third IGR element, C, in the thirteenth triplet. T-unit 12 is indicated using the t-unit nomenclature, simply identifying the entire triplet’s location within the array. The t-unit sequence begins at the start of the pre-tRNA<sup>Ser</sup> gene, and ends at the termination of the pre-tRNA<sup>Tyr</sup>B gene, as indicated by the red, translucent rectangle. T-unit 13, on the right, designates its individual elements by their genetic type—the type of pre-tRNA gene they are (based on their designated amino acid), or the relative location of the intergenic sequence. IGR A lies between the t-units (inter-triplet). IGR A numbering indicates that it precedes the t-unit of the same number, such that IGR A13 is the sequence between t-unit 12 and t-unit 13; note that the exception to this rule is IGR A28 which immediately follows the terminal 3′ t-unit 27. (<bold>c</bold>) The triplet array in <italic>Al</italic> array. As the pre-tRNA genes in the array are not annotated, the genes are represented by colored boxes on top of the genomic G + C background (the orange background indicates an unsequenced region.) Pre-tRNA<sup>Ser</sup> genes are represented by red boxes, and pre-tRNA<sup>Tyr</sup> genes are represented by blue boxes. Note the violation of the three pre-tRNA gene rules in <italic>Al</italic> t-units 0 and 1</p></caption><graphic xlink:href="239_2010_9350_Fig1_HTML" id="MO1"></graphic></fig></p></sec><sec id="Sec4"><title>Gene and Intergenic Sequences</title><p>A complete list of the <italic>At</italic> pre-tRNA genes used in this study can be found at <ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/5bq3o9">http://tinyurl.com/5bq3o9</ext-link> as well as in a spreadsheet in the supplementary file; no annotations existed for the <italic>Al</italic> pre-tRNA genes at the time of this study. All <italic>At</italic> DNA sequences were obtained through The Arabidopsis Information Resource (TAIR), version eight of the <italic>Arabidopsis thaliana</italic> genome. A Perl script was written to parse the array’s DNA sequence, given the annotated gene start/stop positions and knowing the position of the IGRs relative to the pre-tRNA genes. The pre-tRNA<sup>Ser</sup> gene sequences were taken as annotated; however, the pre-tRNA<sup>Tyr</sup> gene models were incomplete. This was apparent through comparisons of the annotated pre-tRNA<sup>Tyr</sup> genes in the array to other, complete gene models both in the <italic>At</italic> genome, and an <italic>At</italic> gene stored in Genbank (M35958.1; <ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/n3f5vh">http://tinyurl.com/n3f5vh</ext-link>). The pre-tRNA<sup>Tyr</sup> gene models were truncated before the gene’s intron sequence (Goodman et al. <xref ref-type="bibr" rid="CR14">1977</xref>). Ten incomplete gene sequences were aligned and compared using blastn (Altschul et al. <xref ref-type="bibr" rid="CR2">1990</xref>; E-value Cutoff: 0.001, Match: 1, Mismatch: −2, Gap Open Penalty: 5, Gap Extension Penalty: 2, Word Size: 7, No Filter on query sequence) to gain a consensus gene sequence and length, and we took the sequence of the most conservative high-scoring segment pair generated by the alignment algorithm. Complete gene models were also predicted using tRNAscan-SE 1.21 (Lowe and Eddy <xref ref-type="bibr" rid="CR23">1997</xref>) and these were concordant with the models obtained using blastn.</p><p>Version one of the <italic>Arabidopsis lyrata</italic> genome was obtained from the Joint Genome Institute (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/l9njsq">http://tinyurl.com/l9njsq</ext-link>). While no pre-tRNA genes were annotated in the region syntenic to the t-unit array in <italic>At</italic>, eleven full length and one partial t-units were found using blastn sequence similarity searches as well as tRNAscan-SE analysis as described above. As this is a draft genome using whole genome shotgun sequencing, we mapped the pair-end sequencing coverage of the region used in our analysis. We found overall coverage to be approximately 8× and had paired-end coverage spanning all parts of this region, but an increase in fold coverage on the 3′ end of the array. (Supplementary Fig. 1.) This may reflect sequencing sampling bias, or additional copies of t-units in unsequenced regions of the genome. All the raw DNA sequences used in this study (genes, intergenic regions, and triplet units) are available in FASTA format in the supplementary file.</p></sec><sec id="Sec5"><title>Syntenic Alignments and Assessment</title><p>We evaluated synteny using the online comparative genomics package CoGe and the method outlined in Lyons et al. (<xref ref-type="bibr" rid="CR26">2008</xref>). First, an approximately 165.4 kb sequence (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/6qghol">http://tinyurl.com/6qghol</ext-link>) encompassing the <italic>At</italic> array (Chr. 1, 21190739-21356171) was retrieved, and annotated coding sequences within the region (5′ and 3′ of the pre-tRNA array) were catalogued (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/5z799u">http://tinyurl.com/5z799u</ext-link>). Using the catalogued CDS sequences as queries, a tblastx search (Altschul et al. <xref ref-type="bibr" rid="CR2">1990</xref>) was then run against the <italic>Carica papaya</italic> (<italic>Cp</italic>) (Ming et al. <xref ref-type="bibr" rid="CR29">2008</xref>), <italic>Vitis vinifera</italic> (<italic>Vv</italic>) (Jaillon et al. <xref ref-type="bibr" rid="CR19">2007</xref>), <italic>Al</italic>, and <italic>At</italic> genomes using CoGeBlast (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/5udsyy">http://tinyurl.com/5udsyy</ext-link>); the genomic regions of the best blast hits were then analyzed for synteny using GEvo (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/d7mzhz">http://tinyurl.com/d7mzhz</ext-link>) and the BlastZ alignment algorithm (Schwartz et al. <xref ref-type="bibr" rid="CR39">2003</xref>). Synteny was determined by identifying a collinear series of putatively homologous genes (Lyons and Freeling <xref ref-type="bibr" rid="CR25">2008</xref>). Four syntenic <italic>At</italic> and <italic>Al</italic> regions (including the region containing the array) were identified in addition to a single syntenic region from both grape and papaya (Supplementary Fig. 2). The ten syntenic genomic regions were as follows: <italic>Cp</italic>—supercontig_48: 1,128,741-1,287,100; <italic>Vv</italic>—Chr. 14: 3,198-528,677 (reverse complement); <italic>Al</italic>—scaffold_1: 3,597,899-3,652,421, scaffold_1: 32,768,291-32,878,607, scaffold_3: 2,260,368-2,319,437, scaffold_6: 12,044,518-12,265,938 (reverse complement); <italic>At</italic>—Chr. 1: 3,039,795-3,100,826, Chr. 1: 21,255,709-21,334,368, Chr. 3: 1,735,349-1,780,617, and Chr. 5: 9,269,098-9,473,551 (reverse complement).</p></sec><sec id="Sec6"><title>Multiple Sequence Alignment</title><p>Sequences were initially aligned using ClustalW 1.83 (Thompson et al. <xref ref-type="bibr" rid="CR46">1994</xref>) using the default settings for a DNA multiple sequence alignment (Open Gap Penalty: 15, Gap Extension Penalty: 6.66, no Weight Transition, Matrix: IUB). The alignment generated by ClustalW was then manually refined in the program Mesquite (Maddison and Maddison <xref ref-type="bibr" rid="CR27">2009</xref>). The final alignment files used in subsequent phylogenetic reconstruction are available in the supplementary materials section in NEXUS format.</p></sec><sec id="Sec7"><title>Phylogenetic Tree Reconstruction</title><p>Maximum likelihood (ML) trees were constructed using Paup* 4.0b10 (Swofford <xref ref-type="bibr" rid="CR43">2003</xref>). Evolutionary models for the likelihood analyses were selected by ModelTest 3.7 (Posada and Crandall <xref ref-type="bibr" rid="CR34">1998</xref>). Trees were visualized using FigTree (Rambaut <xref ref-type="bibr" rid="CR35">2006</xref>–2008). Unrooted ML trees were generated for distinct t-unit elements using a heuristic search with 1,000 repetitions, sequences added randomly, and bisection–reconnection (TBR) branch swapping.</p><p>Bayesian inference trees were also used to reconstruct t-unit phylogeny as well as the <italic>At</italic> IGR A phylogeny. Bayesian analysis was performed by MrBayes v3.1.2 (Ronquist and Huelsenbeck <xref ref-type="bibr" rid="CR37">2003</xref>) using a Markov Chain Monte Carlo (MCMC) run, with t-unit elements partitioned into five unlinked character-sets (genes and IGR sequence), and each partition using a distinct model selected as the best-fit by MrModelTest 2.2 (Nylander <xref ref-type="bibr" rid="CR31">2004</xref>). The simulation was run for ten million generations, distributed over six chains. Burn-in was set at 25% of the trees generated to calculate posterior probabilities. Bayesian inference trees were constructed using the Cyberinfrastructure for Phylogenetic Research (CIPRES) web portal (Miller et al. <xref ref-type="bibr" rid="CR28">2009</xref>).</p><p>Congruence of phylogenetic trees was tested by manually mapping trees for each t-unit element to their genomic locations, and further assessed using an incongruence length difference (ILD) test of homogeneity under a parsimony tree building criterion in Paup*, using 1,000 repetitions for the IGR A to t-unit tests, and 26 repetitions for the intra-t-unit test.</p></sec><sec id="Sec8"><title>Array Pre-tRNA Gene Functionality Data</title><p>Data on functionality of pre-tRNA genes within the array was obtained from previously published data (Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>; Beier and Beier <xref ref-type="bibr" rid="CR4">1992</xref>). Cloned sequences generated from the previous literature (Beier et al. <xref ref-type="bibr" rid="CR5">1991</xref>), found in Genbank (X54368, X54369, X54370, X54371, X543672, X54373, X54374; <ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/kjo7ob">http://tinyurl.com/kjo7ob</ext-link>), were compared to sequences in the array, and functionality data were assigned to identical sequences used in both the cloned gene in vitro analysis and in the array genomic sequence.</p></sec><sec id="Sec9"><title>Tyrosine Usage Measurement</title><p>Data on prevalence of tyrosine codons in genomic coding sequence for both <italic>Arabidopsis</italic> species were obtained using the CoGe application OrganismView and analyzed via a two-sided binomial proportions test using the software package R (Team <xref ref-type="bibr" rid="CR45">2007</xref>).</p></sec><sec id="Sec10"><title>Modeling and Simulating Unequal Crossing-Over</title><p>Unequal crossing-over (UCO) models were generated under the following assumptions:<list list-type="order"><list-item><p>UCO is the sole mechanism responsible for array copy number variation</p></list-item><list-item><p>The array structures were formed using a parsimonious number of expansion events</p></list-item><list-item><p>A given array becomes fixed in the population, such that the crossover events occur between identical arrays (i.e., no distinguishing crossover between homologous chromosomes versus sister chromatids).</p></list-item></list></p><p>Models were devised using these assumptions and the results from the t-unit Bayesian inference tree that t-units group exclusively by species. Only the t-unit phylogeny was taken into account when the models were devised, thus allowing non t-unit data (IGR A sequences) to provide independent validation of the models.</p><p>Using this model, a custom program was written in Perl to stochastically simulate the effects of UCO on orthologous array divergence, named “Array Divergence Simulation.pl” (Supplementary file). The program tracks the divergence of two hypothetical lineages with an initial identical tandem array of gene-unit elements over multiple generations. The stochastic simulation algorithm is as such<list list-type="order"><list-item><p>Initialize two identical arrays of size <italic>n</italic></p></list-item><list-item><p>Choose two crossover sites at random in the set 0 to <italic>n</italic>, inclusive, for each array, giving a total of <italic>n</italic> + 1 possible crossover sites within an array of size <italic>n</italic>.<list list-type="order"><list-item><p>Crossover must have at least one positionally overlapping unit. If the sites chosen do not allow this, sites are reselected by returning to step 2.</p></list-item></list></p></list-item><list-item><p>Generate crossover by swapping array units distal to their respective crossover sites.</p></list-item><list-item><p>Randomly select one of the two resultant arrays for retention into the next cycle.</p></list-item><list-item><p>Return to 2, unless<list list-type="order"><list-item><p>One of the two lineage-specific arrays has been truncated down to a single unit, and thus is unable to perform unequal crossover</p></list-item><list-item><p>Both lineage-specific arrays are homogenous for the same single unit (complete convergence).</p></list-item><list-item><p>Neither lineage-specific array contains any identical units (complete divergence).</p></list-item></list></p></list-item></list></p><p>Each start array consists of gene-units labeled 1 to <italic>n</italic>, where <italic>n</italic> is the initialized size of the array at the onset of the simulation. Each generation permits one crossing-over event (equal or unequal) within the array. The total number of crossover events (CO<sub>T</sub>) possible in any given UCO step is<disp-formula id="Equa"><alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ {\text{CO}}_{\text{T}} = (n + 1) + 2[n + (n - 1) + \cdots + 2] = (n + 1)^{2} - 2$$\end{document}</tex-math><graphic xlink:href="239_2010_9350_Article_Equa.gif" position="anchor"></graphic></alternatives></disp-formula></p><p>We can calculate the expected final array size after a recombination event <italic>E</italic>(<italic>A</italic><sub><italic>n</italic></sub>)<disp-formula id="Equb"><alternatives><tex-math id="M2">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ E(A_{n} ) = {\frac{{\sum {[{\text{CO}}E_{n} \cdot N_{n} \cdot P(I)]} }}{{{\text{CO}}_{\text{T}} }}}$$\end{document}</tex-math><graphic xlink:href="239_2010_9350_Article_Equb.gif" position="anchor"></graphic></alternatives></disp-formula>where COE<sub><italic>n</italic></sub> is the possible number of crossover events given the start size of <italic>n</italic>, <italic>N</italic><sub><italic>n</italic></sub> is the size of the array after the crossover event given start size of <italic>n</italic>, and <italic>P</italic>(<italic>I</italic>) is the probability that a particular product of UCO will be inherited into the next generation. This expands to<disp-formula id="Equc"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ E(A_{n} ) = {\frac{{(n + 1)(1)(n) + 2(n)\left( {\frac{1}{2}} \right)[(n + 1) + (n - 1)] + 2(n - 1)\left( {\frac{1}{2}} \right)[(n + 2) + (n - 2)] + \cdots + 2(2)\left( {\frac{1}{2}} \right)[(1) + (2n - 1)]}}{{(n + 1)^{2} - 2}}} $$\end{document}</tex-math><graphic xlink:href="239_2010_9350_Article_Equc.gif" position="anchor"></graphic></alternatives></disp-formula>which can be further simplified<disp-formula id="Equd"><alternatives><tex-math id="M4">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$ \begin{aligned} E(A_{n} ) & = {\frac{(n + 1)(n) + (n)[2n] + (n - 1)[2n] + \cdots + 2[2n]}{{(n + 1)^{2} - 2}}} \\ E(A_{n} ) & = {\frac{n[(n + 1) + 2[n + (n - 1) + \cdots + 2]]}{{(n + 1)^{2} - 2}}} \\ E(A_{n} ) & = {\frac{{n[(n + 1)^{2} - 2]}}{{(n + 1)^{2} - 2}}} = n \\ \end{aligned} $$\end{document}</tex-math><graphic xlink:href="239_2010_9350_Article_Equd.gif" position="anchor"></graphic></alternatives></disp-formula></p><p>Thus, the expected post-recombination array size <italic>E</italic>(<italic>A</italic><sub><italic>n</italic></sub>) is equal to the initial array, providing our simulation with a parameter for neutral expansion and contraction (no bias toward either), neither mutagenesis of gene-unit elements in the arrays nor any process of gene conversion being permitted, so that simulated array divergence is solely due to gain / loss of gene-unit elements due to UCO events.</p><p>The simulation is repeated for 100,000 replicates for each starting array size of 2–40 elements, and the number of crossover events are computed and statistically analyzed for mean array divergence time. Note that only arrays that completely diverge are used to calculate divergence time, as the simulation is only concerned with the number of crossover events to array divergence and not the likelihood of two orthologous arrays reaching divergence. Another important point is that the two diverging arrays are kept isolated from one another; crossover only occurs between an array within a lineage, not across, as one expects in diverging lineages (i.e., post-speciation events).</p></sec></sec><sec id="Sec11"><title>Results</title><sec id="Sec12"><title>Date Boundaries of T-Unit Array Formation</title><p>A comparative genomics approach was utilized to determine the relative age of the t-unit array (Supplementary Fig. 2). <italic>Carica papaya</italic> (<italic>Cp</italic>) (same order as <italic>At</italic>, Brassicales) and <italic>Vitis vinifera</italic> (<italic>Vv</italic>) (basal eurosid) were chosen as outgroups to the <italic>Arabidopsis</italic> lineage. Subsequent to the divergence of <italic>Cp</italic> and the <italic>Arabidopsis</italic> lineage, the <italic>Arabidopsis</italic> lineage underwent two sequential whole genome duplication (WGD) events (Bowers et al. <xref ref-type="bibr" rid="CR6">2003</xref>; Ming et al. <xref ref-type="bibr" rid="CR29">2008</xref>). These WGD events resulted in ten syntenic genomic regions among these four taxa: four each from <italic>At</italic> and <italic>Al</italic>, and one each from <italic>Cp</italic> and <italic>Vv</italic>. Neither outgroup genome, <italic>Cp</italic> nor <italic>Vv</italic>, contains any pre-tRNA genes within their syntenic regions, nor the three other syntenic regions of <italic>At</italic>. The only other syntenic region containing the array is the orthologous <italic>Al</italic> genomic region. By parsimony, these results indicate that the pre-tRNA genes in the array are not ancestral to the array’s current locus, and likely transposed into the region after the most recent whole genome duplication event in the <italic>Arabidopsis</italic> lineage, and prior to the lineage divergence of <italic>At</italic> and <italic>Al</italic>.</p><p>In order to further refine these dates, a random t-unit was selected to query the NCBI whole-genome shotgun sequence database of <italic>Brassica rapa</italic> (<ext-link ext-link-type="uri" xlink:href="http://tinyurl.com/nh8lgk">http://tinyurl.com/nh8lgk</ext-link>) using the blastn alignment algorithm. No significant hits covering an entire t-unit, nor a series of hits covering the t-units’ tRNA constituents were found (data not shown), likely dating the array’s formation as subsequent to the divergence of <italic>Brassica</italic> and <italic>Arabidopsis</italic> lineages.</p></sec><sec id="Sec13"><title>Origin of tRNA<sup>Tyr</sup> Doublet in T-Unit and Doublet Conservation</title><p>Owing to the tandem appearance of the pre-tRNA<sup>Tyr</sup> doublet, we hypothesized that the doublet in tail-to-head configuration formed from a simple, tandem duplication prior to t-unit expansion (perhaps by inexact reciprocal recombination). A Bayesian inference tree was constructed using all pre-tRNA<sup>Tyr</sup> genes from the <italic>At</italic> and <italic>Al</italic> arrays (Supplementary Fig. 3). The arrays’ pre-tRNA<sup>Tyr</sup> A and B genes cluster into distinct clades. The distinct clustering of the pre-tRNA<sup>Tyr</sup>A and the pre-tRNA<sup>Tyr</sup>B gene clades in both <italic>At</italic> and <italic>Al</italic> indicates that the doublet order is never violated during t-unit duplication; if it were, then intercalation between A and B genes would be observed in the phylogeny due to contradictions in gene locations versus gene homology. Furthermore, while pre-tRNA<sup>Tyr</sup>A genes intercalate between <italic>At</italic> and <italic>Al</italic>, the pre-tRNA<sup>Tyr</sup>B genes do not. The conservation of doublet order evidences that the pre-tRNA<sup>Tyr</sup> subsection of the t-unit is conserved during t-unit duplication in both lineages. The exceptions to this rule are the non-canonical t-units 0 and 1 in <italic>Al</italic>.</p><p>T-unit 0 in <italic>Al</italic> lacks a pre-tRNA<sup>Ser</sup> gene, and thus is in violation of the t-unit conservation rule. However, without further evidence, it is impossible to distinguish between an expansion event that failed to copy the pre-tRNA<sup>Ser</sup> gene and a gene loss event subsequent to duplication. T-unit 1 in <italic>Al</italic> has an additional pre-tRNA<sup>Tyr</sup> gene, and thus is also in violation of the t-unit conservation rule. The last two pre-tRNA<sup>Tyr</sup> genes in t-unit 1 both clade as B genes; however, they do not group together: pre-tRNA<sup>Tyr</sup>B1 clades basally with most of the other pre-tRNA<sup>Tyr</sup>B genes in the <italic>Al</italic> array, and pre-tRNA<sup>Tyr</sup>C1 clades with pre-tRNA<sup>Tyr</sup>B9.</p></sec><sec id="Sec14"><title>Conservation of T-Unit During Expansion</title><p>Manual comparison of ML trees generated from each t-unit element in the <italic>At</italic> array (three pre-tRNA genes and two intergenic regions for a total of five distinct trees) showed that all the trees were congruent with one another, and without conflicts. However, the IGR A tree was not congruent with any of the other trees (Supplementary Fig. 4). This was further assessed using an ILD test of homogeneity. While <italic>At</italic> t-units partitioned into individual elements showed congruence (<italic>P</italic> > 0.25), IGR A sequences were found to be incongruent with t-unit phylogenies either 5′ (<italic>P</italic> = 0.001) or 3′ (<italic>P</italic> = 0.001) of the t-units. The congruence in t-unit element phylogenies, in conjunction with the data confirming conservation of the pre-tRNA<sup>Tyr</sup> doublet order (Supplementary Fig. 3), supports the hypothesis that the t-unit has been the array’s basal unit of duplication and expansion, and that recombination leading to copy number variation occurred in the IGR A sequences. While <italic>Al</italic> t-units exist in fewer copy number and have higher sequence identity to one another, <italic>Al</italic>’s pre-tRNA<sup>Tyr</sup> genes show a similar conservation of gene order to <italic>At</italic>’s. We inferred from this evidence, along with the similar t-unit structure and homology between the two arrays, that t-unit conservation occurs in the <italic>Al</italic> array as well, sans the two aforementioned exceptions in <italic>Al</italic> t-units 0 and 1.</p><p>Using blastn, we searched the At genome for additional t-unit loci; however, none were identified. Of the pre-tRNA genes in this array, seven additional loci in <italic>At</italic> were 100% identical to one or more of the pre-tRNA<sup>Ser</sup> genes. No other loci were 100% identical to any of the pre-tRNA<sup>Tyr</sup> genes. As in <italic>At</italic>, we failed to identify any non-array t-unit loci in the <italic>Al</italic> genome. Of the pre-tRNA genes in the <italic>Al</italic> array, five additional loci in <italic>Al</italic> were 100% identical to one or more of the pre-tRNA<sup>Ser</sup> genes, and none were 100% identical to any of the pre-tRNA<sup>Tyr</sup> genes.</p></sec><sec id="Sec15"><title>Sequence Conservation Within T-Units and Tyrosine tRNA Gene Dosage</title><p>Phylogenetic analysis of all array t-units from both <italic>At</italic> and <italic>Al</italic> shows that t-units clade by paralogs within a taxa rather than orthologs between taxa (Fig. <xref rid="Fig2" ref-type="fig">2</xref>), concurrent with concerted evolution (Zimmer et al. <xref ref-type="bibr" rid="CR49">1980</xref>) occurring within each <italic>Arabidopsis</italic> lineage. Using the previous literature on functionality of t-unit pre-tRNA genes cloned from the <italic>At</italic> t-unit array (Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>; Beier and Beier <xref ref-type="bibr" rid="CR4">1992</xref>), we noted which t-unit genes were capable of producing a functional, mature tRNA product (Supplementary Fig. 5). A number of the genes in the t-unit array are reported not to produce a mature tRNA product in vitro, with non-functional and functional pre-tRNA genes concentrated toward the ends and central region of the <italic>At</italic> array, respectively.<fig id="Fig2"><label>Fig. 2</label><caption><p>Bayesian Inference tree of <italic>At</italic> and <italic>Al</italic> t-units. Bayesian posterior probabilities are indicated by numbers next to nodes. The root is arbitrarily drawn at the tree mid-point and should not be used to infer the actual root of the tree; the “Ancestral Array” and “Speciation” designations are to illustrate common ancestry. Branch lengths are substitutions per site</p></caption><graphic xlink:href="239_2010_9350_Fig2_HTML" id="MO6"></graphic></fig></p><p>In order to test whether positive selection for pre-tRNA<sup>Tyr</sup> gene dosage was the driving force behind array expansion, we analyzed tyrosine amino acid usage across coding regions in <italic>At</italic> and <italic>Al</italic>. A total of 369,223 tyrosine codons were found in the <italic>At</italic> genome (out of 13,133,642 total codons), and 334,108 tyrosine codons were found in the <italic>Al</italic> genome (out of 11,804,715 codons total). A two-sided binomial proportions test showed a very slight lower usage of tyrosine residues in coding regions in <italic>At</italic> compared to <italic>Al</italic> (<italic>P</italic> < 0.005). Using CoGeBlast, we concluded that, aside from the t-unit array, both genomes contain an equal number of pre-tRNA<sup>Tyr</sup> genes (data not shown). As the <italic>At</italic> array contains over two times more pre-tRNA<sup>Tyr</sup> genes compared to the <italic>Al</italic> array (54 versus 24, respectively) without an increase in tyrosine residue usage in coding sequences, it is unlikely that the array’s expansion is due to positive selection for tRNA dosage. However, it should be noted that additional copies of t-units in <italic>Al</italic> could exist as evidenced by the increase in sequencing fold-coverage in the 3′ end of the array (Supplementary Fig. 1).</p></sec><sec id="Sec16"><title>Varied Nature of Array Expansion</title><p>Figure <xref rid="Fig3" ref-type="fig">3</xref> maps <italic>At</italic>’s t-unit phylogenetic tree onto the physical map of the array, revealing the varied nature of the array’s expansion. Multiple putative sites of expansion are observed, involving both single t-units and intercalated sets of two or more t-units. In the array center, t-units 8–14 cluster and, owing to the step-like pattern observed in the t-unit phylogeny, appear to have been duplicated through local tandem duplication events. Interestingly, the array ends show an intercalating pattern between potentially monophyletic clades based on the t-unit phylogeny. T-units 1–6 show two distinct doublets of t-units (t-units 1 and 3; t-units 5 and 7) intercalated by a triplet of t-units (t-units 2, 4, and 6), with the t-unit triplet clading with the t-unit 1 and 3 doublet. T-units 15–23 also show a similar pattern of intercalation, with the even t-units (16, 18, 20, and 22) clustering into a potentially monophyletic group distinct from the odd t-units (15, 17, 19, 21, and 23), which themselves form a potentially monophyletic group. These regions of the array showing intercalation of distinct phylogenetic clusters support expansion events involving multiple t-units as the basal unit of duplication.<fig id="Fig3"><label>Fig. 3</label><caption><p>Mapping of <italic>At</italic> t-unit Bayesian inference tree onto physical map of <italic>At</italic> array. Each triplet is encircled and its position noted numerically within the oval. Imposed on the physical map is the Bayesian inference tree of the triplet units, with each leaf of the tree linked to the appropriate t-unit in the physical array. The array is broken into four regions of duplication by the phylogenetic tree: A, B, C, and D, suggestive of four independent regions of expansion. The branch lengths are arbitrary in this tree</p></caption><graphic xlink:href="239_2010_9350_Fig3_HTML" id="MO7"></graphic></fig></p><p>The <italic>Al</italic> array demonstrates a similar pattern of t-unit expansion (Fig. <xref rid="Fig4" ref-type="fig">4</xref>). T-unit 1 contains two pre-tRNA<sup>Tyr</sup>B genes, each with a different phylogenetic signal—the first clades basally with the other pre-tRNA<sup>Tyr</sup>B genes, and the second clades with pre-tRNA<sup>Tyr</sup>B9. In order to align t-unit 1 with the rest of the t-units in the array, only one of the two pre-tRNA<sup>Tyr</sup>B sequences was included. When the first pre-tRNA<sup>Tyr</sup>B is included in the alignment, t-unit 1 clades with t-unit 0; otherwise, if the second is included, then t-unit 1 clades with t-unit 9. In the Bayesian phylogenetic trees, only the second pre-tRNA<sup>Tyr</sup>B gene, pre-tRNA<sup>Tyr</sup>C1, was included in <italic>Al</italic> t-unit 1. As seen in <italic>At</italic>, <italic>Al</italic>’s array may also involve multiple expansion sites and t-units per expansion step.<fig id="Fig4"><label>Fig. 4</label><caption><p>Mapping of <italic>Al</italic> t-unit Bayesian inference tree onto physical map of <italic>Al</italic> array. Above is the array, with each t-unit encircled, and its position noted numerically below the oval. Imposed on the physical map is the Bayesian inference tree of the <italic>Al</italic> t-units, with each taxa in the tree linked to the appropriate t-unit in the physical array. An intercalated pattern of t-unit duplication, observed in the <italic>At</italic> array (Fig. <xref rid="Fig3" ref-type="fig">3</xref>), is observed here as well. The branch connecting t-units 0, 1, and 9 is drawn as a dotted line due to an inconsistent phylogenetic signal from the pre-tRNA<sup>Tyr</sup> genes in t-unit 1. In order to align the sequences, only two of the three pre-tRNA<sup>Tyr</sup> genes in the t-unit can be used to build the tree. If the last gene is excluded, then t-unit 1 clades with t-unit 0; however, if the middle gene is excluded, then t-unit 1 clades with t-unit 9. The branch lengths are arbitrary in this tree</p></caption><graphic xlink:href="239_2010_9350_Fig4_HTML" id="MO8"></graphic></fig></p></sec><sec id="Sec17"><title>Unequal Crossing-Over as a Model for Array Expansion and Validation by IGR A Phylogeny</title><p>Owing to the step-wise pattern of duplication observed, and prior evidence indicating that unequal crossing-over (UCO) can cause copy number variation in tandem gene arrays (Smith <xref ref-type="bibr" rid="CR41">1976</xref>; Tartof <xref ref-type="bibr" rid="CR44">1988</xref>; Achaz et al. <xref ref-type="bibr" rid="CR1">2000</xref>), we investigated whether a model using UCO as the sole mechanism of t-unit copy number expansion could yield the observed pattern of t-unit formation seen in Fig. <xref rid="Fig3" ref-type="fig">3</xref>. A model of UCO step-wise formation of t-units 1–7 in <italic>At</italic> (Supplementary Fig. 6) was devised under the following three assumptions: (i) UCO is the sole mechanism of t-unit copy number expansion (ii) a parsimonious model of array formation in which the fewest number of UCO events necessary to form the observed expansion patterns occurred, and (iii) crossing-over only occurred in the IGR A region, as the phylogenetic data suggest (namely, the incongruity between the t-unit phylogeny and the IGR A phylogeny). After the parsimonious UCO models were created, the predicted phylogenies of recombinant IGR A sequences created by UCO were compared to the actual IGR A phylogeny (Supplementary Fig. 4). The model-predicted phylogeny was congruent with the phylogenetic reconstruction.</p><p>In order to test whether the observed t-unit sequence divergence between <italic>At</italic> and <italic>Al</italic> could be a result of rapid t-unit turnover from UCO events, a computer program was written to stochastically simulate UCO events between two arrays, each starting with an identical array (common origin) and then evolving independently after a lineage divergence event. Crossover events could only occur between gene units, and equal crossover was permitted. Using array start sizes of 2–40 and with 100,000 repetitions of each simulation for a given start size, the mean number of crossover events to complete array divergence ranged from 7.62 to 16.9, with a local minimum at an array start size of four (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). These simulations show that unequal crossover can cause tandem arrays in two species with a common ancestry to quickly diverge post-speciation, leading to the observed lineage-specific divergence of concerted evolution.<fig id="Fig5"><label>Fig. 5</label><caption><p>Correlation between start array size and number of crossover events until complete divergence of the array following a lineage divergence event, based on simulation data from Array Divergence Simulation.pl. Blue dots denote the mean number of crossover events to obtain complete divergence for 100,000 repetitions. The error bars represent one standard deviation from the mean</p></caption><graphic xlink:href="239_2010_9350_Fig5_HTML" id="MO9"></graphic></fig></p></sec></sec><sec id="Sec18"><title>Discussion</title><p>Syntenic and phylogenetic examination of the <italic>At</italic> t-unit gene array (81 pre-tRNA genes organized into 27 t-units of three pre-tRNA genes each) resulted in an unexpected pattern of gene expansion that was also observed in a smaller, orthologous and syntenic <italic>Al</italic> array. These data suggest a complex series of t-unit duplications, and combined with empirical modeling and simulation data suggest further that unequal crossover events are the likely mechanism of expansion. Interestingly, the array does not appear to have expanded from a single locus within the array. Rather, the array expanded from multiple loci of sub-local duplication consisting of varying units of duplication, nearly all of which involved one or more of the intact t-units (Figs. <xref rid="Fig3" ref-type="fig">3</xref> and <xref rid="Fig4" ref-type="fig">4</xref>).</p><sec id="Sec19"><title>Age of the Array and T-Unit Formation</title><p>Since the radiation of the eurosid super-order, <italic>Vv</italic> (basal rosid) and <italic>Cp</italic> (Brassicales) have experienced no WGD events while the <italic>Arabidopsis</italic> lineage has undergone two sequential WGD events (Ming et al. <xref ref-type="bibr" rid="CR29">2008</xref>). These two WGD events duplicated any ancestral region of the <italic>At</italic> genome twice, thereby generating three syntenic regions for any given genomic region of the <italic>At</italic> genome. By virtue of <italic>At</italic> and <italic>Al</italic> sharing a recent common ancestor, any intra-genomic syntenic regions found in <italic>At</italic> are likely to be shared in <italic>Al</italic>; this is the case for the genomic region in which the <italic>At</italic> t-unit array is found. Comparison of the genomic region in <italic>At</italic> containing the t-unit array to the three <italic>At</italic> intra-genomic syntenic regions, the four syntenic regions of <italic>Al</italic>, and the single syntenic regions of grape and papaya, respectively, revealed the existence of the t-unit array only on chromosome one of <italic>At</italic> and its orthologous syntenic region in <italic>Al</italic> (Supplementary Fig. 2). Neither outgroup syntenic regions, nor the other three syntenic regions in the <italic>At</italic> genome contain any pre-tRNA genes; this is also the case with the <italic>Al</italic> genome. Through parsimony, this indicates that all the pre-tRNA genes in the array did not exist in their current genomic location prior to the <italic>Arabidopsis</italic> lineage’s most recent WGD event. Furthermore, no orthologous t-unit array was found in the <italic>Brassica</italic> lineage, which diverged from the <italic>Arabidopsis</italic> lineage after the aforementioned most recent WGD event. However, until the <italic>B. rapa</italic> genome is completely sequenced and assembled, it remains possible that an orthologous t-unit array does exist in its genome. In addition, more genomic data on lineages diverging prior to and after the <italic>Arabidopsis</italic>–<italic>Brassica</italic> lineage split will aid in narrowing down the actual time of the array’s formation.</p><p>No other intact t-unit was found in either <italic>Arabidopsis</italic> genome indicating that the t-unit formed at this ancestral locus and was likely duplicated prior to the divergence of <italic>At</italic> and <italic>Al</italic>. As the t-unit arrays in <italic>At</italic> and <italic>Al</italic> are derived from the same ancestral array, their divergence time places a lower boundary on the age of the array. Thus, we may infer that the initial formation of the t-unit array likely occurred after the divergence of <italic>Arabidopsis</italic> and <italic>Brassica</italic> [14.5–20.4 million years ago (Yang et al. <xref ref-type="bibr" rid="CR47">1999</xref>; Bowers et al. <xref ref-type="bibr" rid="CR6">2003</xref>)], and necessarily prior to the divergence of <italic>At</italic> and <italic>Al</italic> [5–6 million years ago (Koch and Kiefer <xref ref-type="bibr" rid="CR21">2005</xref>)].</p><p>Unfortunately, a great unanswered question in this study is the origin of the t-unit prior to t-unit copy expansion. Presumably, the ancestral t-unit formed through a series of translocation events of pre-tRNA genes to the array locus, but there is no direct evidence to support this claim. Discovering an outgroup genome with a proto t-unit array may provide the data necessary to elucidate the apparent de novo formation of the ordered and conserved t-unit gene structure.</p></sec><sec id="Sec20"><title>The Triplet Unit as the Basal Unit of Array Expansion</title><p>Apparent conservation of the t-unit during expansion is a fundamental property of the t-unit array. Although the combined <italic>At</italic> and <italic>Al</italic> pre-tRNA<sup>Tyr</sup> gene phylogeny demonstrated conservation of the pre-tRNA<sup>Tyr</sup> gene doublet within the t-units, these data are insufficient to state whether all elements within the t-units are conserved as a single unit during expansion events. By mapping the phylogeny of each sub-element of the t-unit onto their genomic positions within the array, each tree was shown to be congruent; homogeneity tests of congruence provided a statistical basis for using entire t-units in phylogenetic reconstruction. If elements within individual t-units were not conserved during expansion, different elements would show phylogenetic incongruence when mapped to their genomic locations. This conservation demonstrates that the t-unit is the basal unit of expansion in the <italic>At</italic> array, and also, with few exceptions, in the <italic>Al</italic> array. Even within the <italic>Al</italic> array, the gene order of the t-unit containing three pre-tRNA<sup>Tyr</sup> genes (<italic>Al</italic> t-unit 1) is apparently not violated, with the third gene clading with the pre-tRNA<sup>Tyr</sup>B genes. The incongruence of the IGR A phylogeny with that of the t-units is hypothesized to be a result of IGR A being used as the site of recombination during UCO resulting in recombinant inter-t-unit sequences that do not follow the other array element phylogenies.</p><p>Previous research on the <italic>At</italic> array has shown that a number of the pre-tRNA genes fail to produce a mature tRNA product in vitro (Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>; Beier and Beier <xref ref-type="bibr" rid="CR4">1992</xref>). While the nonfunctional pre-tRNAs appear to be a subset of the total pre-tRNA genes in the array, these putative pseudogenes argue that purifying selection is not acting on all the genes in the array, presumably due to redundancy of function. Compounded by the lack of correlation between pre-tRNA<sup>Tyr</sup> gene copy number in the t-unit array and tyrosine residue usage in coding sequences, the numerous deleterious mutations found in t-unit genes argues that the copy number expansions within the arrays is likely non-adaptive and fluctuations in array size are governed by neutral evolution such as drift. If accurate, we predict a range of copy number variation of t-units among individuals in <italic>Arabidopsis</italic> populations.</p><p>The conservation of the t-unit pattern in both <italic>At</italic> and <italic>Al</italic> argues for either selection against breaking up t-units, or an additional mechanism of duplication that preferably copies multi-gene units. Unlike protein coding genes, pre-tRNA genes contain internal promoters (Hofstetter et al. <xref ref-type="bibr" rid="CR15">1981</xref>), thus negating the necessity for gene duplication to carry additional sequence upstream of the gene to allow for proper transcription by RNA polymerase III. Pre-tRNA genes also require short repeating thymine nucleotides as a transcription termination signal (Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>). However, this termination sequence is only a few nucleotides 3′ of the gene, and is unlikely to provide sufficient selection against breaking up a t-unit during duplication.</p><p>An alternative hypothesis to direct selection for t-unit conservation is simply that the sequences flanking the t-units are hot-spots for recombination. In this model, recombination can occur within t-units, but the high affinity for molecular elements responsible for crossover competitively inhibits intra-t-unit recombination. Sequence analysis of the intergenic regions between t-units (IGR A) showed no significant homology to other loci in the <italic>At</italic> genome (data not shown). While this does not negate the possibility of a recombination hot-spot motif indirectly driving t-unit conservation, no current data supports this model. In addition, t-unit conservation may be driven by DNA topology in the form of nucleosome placement, higher order chromatin packaging, or epigenetic marks that would occlude recombination from occurring within the t-units. Future studies of the Arabidopsis arrays assessing the relationship between the chromatin state of the array and recombination rates within the array may aid in understanding the phenomenon of t-unit conservation during array expansion. In any case, we are currently unable to reconcile the neutral evolution of the array as a whole in terms of copy number variation with the apparent selection for conserved t-units.</p><p>In vitro assays of t-unit gene transcription showed that the presence of pre-tRNA<sup>Ser</sup> genes did not significantly affect the transcription or maturation of the pre-tRNA<sup>Tyr</sup> gene product adjacent to them (Stange et al. <xref ref-type="bibr" rid="CR42">1991</xref>). The presence of many putative pseudogenes within the array suggests that the high copy number allows mutated pre-tRNA gene sequences to fly “under the radar” of selection due to redundancy (Ganley and Kobayashi <xref ref-type="bibr" rid="CR11">2007</xref>). Interestingly, most of the noted pseudogenes are located at the ends of the array; given the nature of UCO, it is difficult to remove array ends, so it is no surprise that they contain the most divergent sequences. However, it remains unclear whether one would observe the same pattern of higher sequence divergence toward the array ends given gene conversion as the mechanism of sequence homogenization/concerted evolution. If the array centers are more prone to homogenization due to higher frequencies of UCO, perhaps selection is then acting upon the sites of unequal crossover instead of on the genes as a whole, where selection for expansion or contraction of t-units results in an array containing some minimum number of functional t-units. However, this could also be neutral and the result of UCO occurring with a higher likelihood near the center of the array than at the ends. A similar pattern of multi-gene duplicated units containing deleterious mutations is observed in tandem rRNA gene arrays (Ganley and Kobayashi <xref ref-type="bibr" rid="CR11">2007</xref>), suggesting a similar mechanism of multi-gene duplication.</p></sec><sec id="Sec21"><title>Case Study for Modeling Array Formation Using Unequal Crossing-Over as Mechanism of Expansion</title><p>For a review of homologous recombination in plants and references to more in-depth reviews of recombination in general, please see Schuermann et al. (<xref ref-type="bibr" rid="CR38">2005</xref>). Following an initial tandem gene duplication, UCO is a known mechanism of subsequent tandem gene copy number expansion (Smith <xref ref-type="bibr" rid="CR41">1976</xref>; Tartof <xref ref-type="bibr" rid="CR44">1988</xref>; Achaz et al. <xref ref-type="bibr" rid="CR1">2000</xref>) and, therefore, a plausible mechanism of t-unit array expansion. Experiments studying meiotic UCO in <italic>At</italic> have shown rates of gene copy number variation on the order of ~10<sup>−6</sup> per array per plant per F<sub>1</sub> meiosis (Jelesko et al. <xref ref-type="bibr" rid="CR20">2004</xref>), or at least one copy number variant in one out of approximately 700 seeds due solely to UCO (Gaut et al. <xref ref-type="bibr" rid="CR13">2007</xref>). Assuming a single meiotic event per year per organism that gives rise to a single progeny, a given array can expect a gene copy number variant of 1 array<sup>−1</sup> Mya<sup>−1</sup> per haploid genome. This is two orders of magnitude higher than the proposed rate of an initial tandem gene duplication arising out of 0.01 gene<sup>−1</sup> Mya<sup>−1</sup> (Lynch and Conery <xref ref-type="bibr" rid="CR24">2000</xref>). Note the former value does not correct for the increased frequency of unequal crossover as array size expands, which gives greater opportunity for misalignment (Eichler <xref ref-type="bibr" rid="CR8">2001</xref>), as has been observed in micro-satellite DNA (Richard and Dujon <xref ref-type="bibr" rid="CR36">2006</xref>). Given the ability of UCO to rapidly expand or contract pre-existing tandem gene arrays, UCO is a prime candidate for mechanistically explaining the formation of the t-unit array and its sequence homogenization in <italic>At</italic> and <italic>Al</italic>. In order to test the hypothesis of UCO as the mechanism of t-unit array expansion, a model was devised for the formation of t-units 1–7 (Supplementary Fig. 6). The model assumes that recombination occurs in the IGR A regions, with each duplication from UCO creating recombinant IGR A sequences. The model predicted specific phylogenetic relationships among the derived chimeric IGR A spacer regions, which were subsequently validated using the actual phylogeny of the IGR A sequences (Supplementary Fig. 4; Supplementary Fig. 6). While this parsimonious model does not presume to be an accurate recreation of the true evolutionary history of the array, it does demonstrate that the observed structure of the array can be reconstructed using UCO as the mechanism of copy number expansion.</p><p>If unequal crossing-over causes copy number variation via rapid expansion and contractions, and copy number is selectively neutral, then one can expect expansion events to occur at approximately the same rate as for contraction events. If so, then the model of UCO inherently predicts a continual loss of heterogeneity within the array, as t-units lost via array contraction are replaced by “homogenized” gene duplicates as the array re-expands. This would expand on the Birth-and-Death model of concerted evolution (Nei and Rooney <xref ref-type="bibr" rid="CR30">2005</xref>) by the explicit inclusion of a coupled birth/death mechanism by UCO. Our computer simulation of UCO in orthologous arrays demonstrated this phenomenon of concerted evolution through intra-array sequence homogenization. There appears to be little evidence to show that the copy number size of the array is adaptive given the variation in array size in the two <italic>Arabidopsis</italic>’ genomic regions, the apparent lack of positive selection for pre-tRNA dosage in relation to the codon usage across the genomes of <italic>At</italic> and <italic>Al</italic>, and the presence of many non-functional pre-tRNA genes within the <italic>At</italic> array. Thus, successive rounds of UCO inherently causes a loss of phylogenetic information during contraction events; when an expansion event occurs to increase copy number in a tandem array, the t-units replacing those lost will not replace the lost evolutionary information. Together this leads to the observed sequence homogenization between <italic>At</italic> and <italic>Al</italic> and the phylogenetic grouping of paralogous versus orthologous t-units, and is evidence of concerted evolution (Zimmer et al. <xref ref-type="bibr" rid="CR49">1980</xref>).</p><p>As the original discoverers of the tandem array astutely noted (Beier et al. <xref ref-type="bibr" rid="CR5">1991</xref>), due to the inferred rapid turnover of t-units, the true evolutionary history of the array will almost certainly prove impossible to determine due to erasure of evolutionary history (Gao and Innan <xref ref-type="bibr" rid="CR12">2004</xref>). However, our findings show that the recent evolutionary history of the array can be distilled, and its mechanisms of expansion can be understood. Furthermore, this comparative research between the <italic>At</italic> and <italic>Al</italic> arrays shows the highly dynamic nature of these arrays’ evolution, with new elements likely being added and removed independently in both lineages. Likewise, these arrays may be short-lived genomic phenomena. Without selective pressure to maintain multiple copies, these arrays can theoretically be reduced to a singlet state. More array sequences from closely related taxa will aid in understanding the mechanisms of tandem gene array evolution if a very recent expansion or contraction event is identified. Future genome projects involving both population level sequencing and strategic sequencing of taxa based on their phylogenetic relationships will further elucidate the mechanisms of evolution of the t-unit array, and tandem gene evolution in general. However, one outstanding feature of these arrays, which is not explained by our analyses is why UCO appears to occur only among IGR A regions, and perhaps additional taxa sampling will shed light into this mystery.</p></sec></sec><sec id="Sec22"><title>Conclusion</title><p>The triplet pre-tRNA (t-unit) gene arrays located on chromosome one of <italic>Arabidopsis thaliana</italic> and <italic>Arabidopsis lyrata</italic> have provided a high-resolution case study of the expansion of a unique triplet of tRNA genes. The t-unit formed after the divergence of <italic>Arabidopsis</italic> and <italic>Brassica</italic> lineages, and before the divergence of <italic>A. thaliana</italic> (<italic>At</italic>) and <italic>A. lyrata</italic> (<italic>Al</italic>) lineages. While the mechanism forming the ancestral t-unit almost certainly involved pre-tRNA gene transposition, it is otherwise an unknown saltation. Expansion from the initial t-unit into an array occurred from multiple origins, and unequal crossing-over (UCO) is likely responsible for the generation of at least part, if not the whole, of the array. UCO as the mechanism of t-unit duplication is supported by phylogentic analysis of intergenic sequences between t-units and comparison of t-unit arrays between <italic>At</italic> and <italic>Al</italic>. A consequence of UCO is that such arrays are subject to contraction as well as expansion, which will reduce overall sequence diversity and cause apparent sequence homogenization. This concerted evolution is evidenced by the apparent monophyly of t-units between <italic>At</italic> and <italic>Al</italic>, and is supported by computational simulation.</p></sec><sec sec-type="supplementary-material"><sec id="AppESM1"><title>Electronic supplementary material</title><supplementary-material id="d28e1766" content-type="local-data"><p>Below is the link to the electronic supplementary material.
</p><media xlink:href="239_2010_9350_MOESM1_ESM.doc" id="MOESM1" mimetype="application" mime-subtype="msword"><caption><p>Supplementary material 1 (DOC 32 kb)</p></caption></media><media xlink:href="239_2010_9350_MOESM2_ESM.pdf" id="MOESM2" mimetype="application" mime-subtype="pdf"><caption><p>Supplementary material 2 (PDF 1181 kb)</p></caption></media><media xlink:href="239_2010_9350_MOESM3_ESM.gz" id="MOESM3" mimetype="application" mime-subtype="x-zip-compressed"><caption><p>Supplementary material 3 (gz 71 kb)</p></caption></media></supplementary-material></sec></sec></body><back><ack><p><italic>Arabidopsis lyrata</italic> sequence data were produced by the US Department of Energy Joint Genome Institute (<ext-link ext-link-type="uri" xlink:href="http://www.jgi.doe.gov/">http://www.jgi.doe.gov/</ext-link>) in collaboration with the user community. This study is supported by the UC Berkeley Undergraduate Research Apprentice Program, the UC Berkeley Summer Undergraduate Research Fellowships program, and by the National Science Foundation grant MCB-0820821. We are grateful to Kelly Nissen and Damon Lisch for their contributive discussions and insight, Brent Mishler for his expertise in phylogenetics, Detlef Weigel for his leading the <italic>A. lyrata</italic> genome sequencing project, Jeremy Schmutz for his help with validating the assembly of <italic>A. lyrata</italic>, and the CIPRES portal for their excellent suite of phylogenetic tools.</p><p><bold>Open Access</bold> This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.</p></ack><ref-list id="Bib1"><title>References</title><ref id="CR1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Achaz</surname><given-names>G</given-names></name><name><surname>Coissac</surname><given-names>E</given-names></name><name><surname>Viari</surname><given-names>A</given-names></name><name><surname>Netter</surname><given-names>P</given-names></name></person-group><article-title>Analysis of intrachromosomal duplications in yeast <italic>saccharomyces cerevisiae</italic>: a possible model for their origin</article-title><source>Mol Biol Evol</source><year>2000</year><volume>17</volume><fpage>1268</fpage><lpage>1275</lpage><pub-id pub-id-type="pmid">10908647</pub-id></mixed-citation></ref><ref id="CR2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname><given-names>SF</given-names></name><name><surname>Gish</surname><given-names>W</given-names></name><name><surname>Miller</surname><given-names>W</given-names></name><name><surname>Myers</surname><given-names>EW</given-names></name><name><surname>Lipman</surname><given-names>DJ</given-names></name></person-group><article-title>Basic local alignment search tool</article-title><source>J Mol Biol</source><year>1990</year><volume>215</volume><fpage>403</fpage><lpage>410</lpage><pub-id pub-id-type="pmid">2231712</pub-id></mixed-citation></ref><ref id="CR3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baumgarten</surname><given-names>A</given-names></name><name><surname>Cannon</surname><given-names>S</given-names></name><name><surname>Spangler</surname><given-names>R</given-names></name><name><surname>May</surname><given-names>G</given-names></name></person-group><article-title>Genome-level evolution of resistance genes in <italic>Arabidopsis thaliana</italic></article-title><source>Genetics</source><year>2003</year><volume>165</volume><fpage>309</fpage><lpage>319</lpage><pub-id pub-id-type="pmid">14504238</pub-id></mixed-citation></ref><ref id="CR4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beier</surname><given-names>D</given-names></name><name><surname>Beier</surname><given-names>H</given-names></name></person-group><article-title>Expression of variant nuclear Arabidopsis tRNA(Ser) genes and pre-tRNA maturation differ in HeLa, yeast and wheat germ extracts</article-title><source>Mol General Genetics</source><year>1992</year><volume>233</volume><fpage>201</fpage><lpage>208</lpage><pub-id pub-id-type="doi">10.1007/BF00587580</pub-id></mixed-citation></ref><ref id="CR5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beier</surname><given-names>D</given-names></name><name><surname>Stange</surname><given-names>N</given-names></name><name><surname>Gross</surname><given-names>HJ</given-names></name><name><surname>Beier</surname><given-names>H</given-names></name></person-group><article-title>Nuclear tRNATyr genes are highly amplified at a single chromosomal site in the genome of <italic>Arabidopsis thaliana</italic></article-title><source>Mol General Genetics</source><year>1991</year><volume>225</volume><fpage>72</fpage><lpage>80</lpage></mixed-citation></ref><ref id="CR6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bowers</surname><given-names>JE</given-names></name><name><surname>Chapman</surname><given-names>BA</given-names></name><name><surname>Rong</surname><given-names>J</given-names></name><name><surname>Paterson</surname><given-names>AH</given-names></name></person-group><article-title>Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events</article-title><source>Nature</source><year>2003</year><volume>422</volume><fpage>433</fpage><lpage>438</lpage><pub-id pub-id-type="doi">10.1038/nature01521</pub-id><pub-id pub-id-type="pmid">12660784</pub-id></mixed-citation></ref><ref id="CR7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>DD</given-names></name><name><surname>Wensink</surname><given-names>PC</given-names></name><name><surname>Jordan</surname><given-names>E</given-names></name></person-group><article-title>A comparison of the ribosomal DNA’s of <italic>Xenopus laevis</italic> and <italic>Xenopus mulleri</italic>: the evolution of tandem genes</article-title><source>J Mol Biol</source><year>1972</year><volume>63</volume><fpage>57</fpage><lpage>73</lpage><pub-id pub-id-type="doi">10.1016/0022-2836(72)90521-9</pub-id><pub-id pub-id-type="pmid">5016971</pub-id></mixed-citation></ref><ref id="CR8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eichler</surname><given-names>EE</given-names></name></person-group><article-title>Recent duplication, domain accretion and the dynamic mutation of the human genome</article-title><source>Trends Genet</source><year>2001</year><volume>17</volume><fpage>661</fpage><lpage>669</lpage><pub-id pub-id-type="doi">10.1016/S0168-9525(01)02492-1</pub-id><pub-id pub-id-type="pmid">11672867</pub-id></mixed-citation></ref><ref id="CR9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Force</surname><given-names>A</given-names></name><name><surname>Lynch</surname><given-names>M</given-names></name><name><surname>Pickett</surname><given-names>FB</given-names></name><name><surname>Amores</surname><given-names>A</given-names></name><name><surname>Yan</surname><given-names>YL</given-names></name><name><surname>Postlethwait</surname><given-names>J</given-names></name></person-group><article-title>Preservation of duplicate genes by complementary, degenerative mutations</article-title><source>Genetics</source><year>1999</year><volume>151</volume><fpage>1531</fpage><lpage>1545</lpage><pub-id pub-id-type="pmid">10101175</pub-id></mixed-citation></ref><ref id="CR10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Freeling</surname><given-names>M</given-names></name></person-group><article-title>Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition</article-title><source>Annu Rev Plant Biol</source><year>2009</year><volume>60</volume><fpage>433</fpage><lpage>453</lpage><pub-id pub-id-type="doi">10.1146/annurev.arplant.043008.092122</pub-id><pub-id pub-id-type="pmid">19575588</pub-id></mixed-citation></ref><ref id="CR11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ganley</surname><given-names>ARD</given-names></name><name><surname>Kobayashi</surname><given-names>T</given-names></name></person-group><article-title>Highly efficient concerted evolution in the ribosomal DNA repeats: total rDNA repeat variation revealed by whole-genome shotgun sequence data</article-title><source>Genome Res</source><year>2007</year><volume>17</volume><fpage>184</fpage><lpage>191</lpage><pub-id pub-id-type="doi">10.1101/gr.5457707</pub-id><pub-id pub-id-type="pmid">17200233</pub-id></mixed-citation></ref><ref id="CR12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname><given-names>L-Z</given-names></name><name><surname>Innan</surname><given-names>H</given-names></name></person-group><article-title>Very low gene duplication rate in the yeast genome</article-title><source>Science (New York, NY)</source><year>2004</year><volume>306</volume><fpage>1367</fpage><lpage>1370</lpage></mixed-citation></ref><ref id="CR13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gaut</surname><given-names>BS</given-names></name><name><surname>Wright</surname><given-names>SI</given-names></name><name><surname>Rizzon</surname><given-names>C</given-names></name><name><surname>Dvorak</surname><given-names>J</given-names></name><name><surname>Anderson</surname><given-names>LK</given-names></name></person-group><article-title>Recombination: an underappreciated factor in the evolution of plant genomes</article-title><source>Nat Rev Genet</source><year>2007</year><volume>8</volume><fpage>77</fpage><lpage>84</lpage><pub-id pub-id-type="doi">10.1038/nrg1970</pub-id><pub-id pub-id-type="pmid">17173059</pub-id></mixed-citation></ref><ref id="CR14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodman</surname><given-names>HM</given-names></name><name><surname>Olson</surname><given-names>MV</given-names></name><name><surname>Hall</surname><given-names>BD</given-names></name></person-group><article-title>Nucleotide sequence of a mutant eukaryotic gene: the yeast tyrosine-inserting ochre suppressor SUP4-o</article-title><source>Proc Natl Acad Sci USA</source><year>1977</year><volume>74</volume><fpage>5453</fpage><lpage>5457</lpage><pub-id pub-id-type="doi">10.1073/pnas.74.12.5453</pub-id><pub-id pub-id-type="pmid">341157</pub-id></mixed-citation></ref><ref id="CR15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hofstetter</surname><given-names>H</given-names></name><name><surname>Kressman</surname><given-names>A</given-names></name><name><surname>Birnstiel</surname><given-names>ML</given-names></name></person-group><article-title>A split promoter for a eucaryotic tRNA gene</article-title><source>Cell</source><year>1981</year><volume>24</volume><fpage>573</fpage><lpage>585</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(81)90348-2</pub-id><pub-id pub-id-type="pmid">7237560</pub-id></mixed-citation></ref><ref id="CR16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Horiuchi</surname><given-names>T</given-names></name><name><surname>Tomizawa</surname><given-names>JI</given-names></name><name><surname>Novick</surname><given-names>A</given-names></name></person-group><article-title>Isolation and properties of bacteria capable of high rates of beta-galactosidase synthesis</article-title><source>Biochim Biophys Acta</source><year>1962</year><volume>55</volume><fpage>152</fpage><lpage>163</lpage><pub-id pub-id-type="doi">10.1016/0006-3002(62)90941-1</pub-id><pub-id pub-id-type="pmid">14448925</pub-id></mixed-citation></ref><ref id="CR17"><mixed-citation publication-type="other">Hurles M (2004) Gene duplication: the genomic trade in spare parts. PLoS Biol 2: e206 EP–e206 EP</mixed-citation></ref><ref id="CR18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jackson</surname><given-names>AP</given-names></name></person-group><article-title>Tandem gene arrays in <italic>Trypanosoma brucei</italic>: comparative phylogenomic analysis of duplicate sequence variation</article-title><source>BMC Evol Biol</source><year>2007</year><volume>7</volume><fpage>54</fpage><pub-id pub-id-type="doi">10.1186/1471-2148-7-54</pub-id><pub-id pub-id-type="pmid">17408475</pub-id></mixed-citation></ref><ref id="CR19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jaillon</surname><given-names>O</given-names></name><name><surname>Aury</surname><given-names>J-M</given-names></name><name><surname>Noel</surname><given-names>B</given-names></name><name><surname>Policriti</surname><given-names>A</given-names></name><name><surname>Clepet</surname><given-names>C</given-names></name><name><surname>Casagrande</surname><given-names>A</given-names></name><name><surname>Choisne</surname><given-names>N</given-names></name><name><surname>Aubourg</surname><given-names>S</given-names></name><name><surname>Vitulo</surname><given-names>N</given-names></name><name><surname>Jubin</surname><given-names>C</given-names></name><name><surname>Vezzi</surname><given-names>A</given-names></name><name><surname>Legeai</surname><given-names>F</given-names></name><name><surname>Hugueney</surname><given-names>P</given-names></name><name><surname>Dasilva</surname><given-names>C</given-names></name><name><surname>Horner</surname><given-names>D</given-names></name><name><surname>Mica</surname><given-names>E</given-names></name><name><surname>Jublot</surname><given-names>D</given-names></name><name><surname>Poulain</surname><given-names>J</given-names></name><name><surname>Bruyère</surname><given-names>C</given-names></name><name><surname>Billault</surname><given-names>A</given-names></name><name><surname>Segurens</surname><given-names>B</given-names></name><name><surname>Gouyvenoux</surname><given-names>M</given-names></name><name><surname>Ugarte</surname><given-names>E</given-names></name><name><surname>Cattonaro</surname><given-names>F</given-names></name><name><surname>Anthouard</surname><given-names>V</given-names></name><name><surname>Vico</surname><given-names>V</given-names></name><name><surname>Fabbro</surname><given-names>CD</given-names></name><name><surname>Alaux</surname><given-names>M</given-names></name><name><surname>Gaspero</surname><given-names>GD</given-names></name><name><surname>Dumas</surname><given-names>V</given-names></name><name><surname>Felice</surname><given-names>N</given-names></name><name><surname>Paillard</surname><given-names>S</given-names></name><name><surname>Juman</surname><given-names>I</given-names></name><name><surname>Moroldo</surname><given-names>M</given-names></name><name><surname>Scalabrin</surname><given-names>S</given-names></name><name><surname>Canaguier</surname><given-names>A</given-names></name><name><surname>Clainche</surname><given-names>IL</given-names></name><name><surname>Malacrida</surname><given-names>G</given-names></name><name><surname>Durand</surname><given-names>E</given-names></name><name><surname>Pesole</surname><given-names>G</given-names></name><name><surname>Laucou</surname><given-names>V</given-names></name><name><surname>Chatelet</surname><given-names>P</given-names></name><name><surname>Merdinoglu</surname><given-names>D</given-names></name><name><surname>Delledonne</surname><given-names>M</given-names></name><name><surname>Pezzotti</surname><given-names>M</given-names></name><name><surname>Lecharny</surname><given-names>A</given-names></name><name><surname>Scarpelli</surname><given-names>C</given-names></name><name><surname>Artiguenave</surname><given-names>F</given-names></name><name><surname>Pè</surname><given-names>ME</given-names></name><name><surname>Valle</surname><given-names>G</given-names></name><name><surname>Morgante</surname><given-names>M</given-names></name><name><surname>Caboche</surname><given-names>M</given-names></name><name><surname>Adam-Blondon</surname><given-names>A-F</given-names></name><name><surname>Weissenbach</surname><given-names>J</given-names></name><name><surname>Quétier</surname><given-names>F</given-names></name><name><surname>Wincker</surname><given-names>P</given-names></name></person-group><article-title>The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla</article-title><source>Nature</source><year>2007</year><volume>449</volume><fpage>463</fpage><lpage>467</lpage><pub-id pub-id-type="doi">10.1038/nature06148</pub-id><pub-id pub-id-type="pmid">17721507</pub-id></mixed-citation></ref><ref id="CR20"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jelesko</surname><given-names>JG</given-names></name><name><surname>Carter</surname><given-names>K</given-names></name><name><surname>Thompson</surname><given-names>W</given-names></name><name><surname>Kinoshita</surname><given-names>Y</given-names></name><name><surname>Gruissem</surname><given-names>W</given-names></name></person-group><article-title>Meiotic Recombination Between Paralogous RBCSB Genes on Sister Chromatids of <italic>Arabidopsis thaliana</italic></article-title><source>Genetics</source><year>2004</year><volume>166</volume><fpage>947</fpage><lpage>957</lpage><pub-id pub-id-type="doi">10.1534/genetics.166.2.947</pub-id><pub-id pub-id-type="pmid">15020479</pub-id></mixed-citation></ref><ref id="CR21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koch</surname><given-names>MA</given-names></name><name><surname>Kiefer</surname><given-names>M</given-names></name></person-group><article-title>Genome evolution among cruciferous plants: a lecture from the comparison of the genetic maps of three diploid species–<italic>Capsella rubella</italic>, <italic>Arabidopsis lyrata</italic> subsp. <italic>petraea</italic>, and <italic>A. thaliana</italic></article-title><source>Am J Bot</source><year>2005</year><volume>92</volume><fpage>761</fpage><lpage>767</lpage><pub-id pub-id-type="doi">10.3732/ajb.92.4.761</pub-id></mixed-citation></ref><ref id="CR22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lewis</surname><given-names>D</given-names></name></person-group><article-title>A relationship between dominance, phenotypic stability and variability, and a theory of alternative genetic pathways</article-title><source>Nature</source><year>1953</year><volume>172</volume><fpage>1136</fpage><lpage>1137</lpage><pub-id pub-id-type="doi">10.1038/1721136a0</pub-id><pub-id pub-id-type="pmid">13111267</pub-id></mixed-citation></ref><ref id="CR23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowe</surname><given-names>TM</given-names></name><name><surname>Eddy</surname><given-names>SR</given-names></name></person-group><article-title>tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence</article-title><source>Nucl Acids Res</source><year>1997</year><volume>25</volume><fpage>955</fpage><lpage>964</lpage><pub-id pub-id-type="doi">10.1093/nar/25.5.955</pub-id><pub-id pub-id-type="pmid">9023104</pub-id></mixed-citation></ref><ref id="CR24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lynch</surname><given-names>M</given-names></name><name><surname>Conery</surname><given-names>JS</given-names></name></person-group><article-title>The evolutionary fate and consequences of duplicate genes</article-title><source>Science</source><year>2000</year><volume>290</volume><fpage>1151</fpage><lpage>1155</lpage><pub-id pub-id-type="doi">10.1126/science.290.5494.1151</pub-id><pub-id pub-id-type="pmid">11073452</pub-id></mixed-citation></ref><ref id="CR25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lyons</surname><given-names>E</given-names></name><name><surname>Freeling</surname><given-names>M</given-names></name></person-group><article-title>How to usefully compare homologous plant genes and chromosomes as DNA sequences</article-title><source>Plant J</source><year>2008</year><volume>53</volume><fpage>661</fpage><lpage>673</lpage><pub-id pub-id-type="doi">10.1111/j.1365-313X.2007.03326.x</pub-id><pub-id pub-id-type="pmid">18269575</pub-id></mixed-citation></ref><ref id="CR26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lyons</surname><given-names>E</given-names></name><name><surname>Pedersen</surname><given-names>B</given-names></name><name><surname>Kane</surname><given-names>J</given-names></name><name><surname>Alam</surname><given-names>M</given-names></name><name><surname>Ming</surname><given-names>R</given-names></name><name><surname>Tang</surname><given-names>H</given-names></name><name><surname>Wang</surname><given-names>X</given-names></name><name><surname>Bowers</surname><given-names>J</given-names></name><name><surname>Paterson</surname><given-names>A</given-names></name><name><surname>Lisch</surname><given-names>D</given-names></name><name><surname>Freeling</surname><given-names>M</given-names></name></person-group><article-title>Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids</article-title><source>Plant Physiol</source><year>2008</year><volume>148</volume><fpage>1772</fpage><lpage>1781</lpage><pub-id pub-id-type="doi">10.1104/pp.108.124867</pub-id><pub-id pub-id-type="pmid">18952863</pub-id></mixed-citation></ref><ref id="CR27"><mixed-citation publication-type="other">Maddison W, Maddison D (2009) Mesquite: a modular system for evolutionary analysis. Version 2.6 URL: <ext-link ext-link-type="uri" xlink:href="http://mesquiteproject.org">http://mesquiteproject.org</ext-link></mixed-citation></ref><ref id="CR28"><mixed-citation publication-type="other">Miller MA, Holder MT, Vos R, Midford PE, Liebowitz T, Chan L, Hoover P, Warnow T (2009) The CIPRES Portals URL: <ext-link ext-link-type="uri" xlink:href="http://www.phylo.org/sub_sections/portal">http://www.phylo.org/sub_sections/portal</ext-link></mixed-citation></ref><ref id="CR29"><mixed-citation publication-type="other">Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Saw JH, Senin P, Wang W, Ly BV, Lewis KLT, Salzberg SL, Feng L, Jones MR, Skelton RL, Murray JE, Chen C, Qian W, Shen J, Du P, Eustice M, Tong E, Tang H, Lyons E, Paull RE, Michael TP, Wall K, Rice DW, Albert H, Wang M-L, Zhu YJ, Schatz M, Nagarajan N, Acob RA, Guan P, Blas A, Wai CM, Ackerman CM, Ren Y, Liu C, Wang J, Wang J, Na J-K, Shakirov EV, Haas B, Thimmapuram J, Nelson D, Wang X, Bowers JE, Gschwend AR, Delcher AL, Singh R, Suzuki JY, Tripathi S, Neupane K, Wei H, Irikura B, Paidi M, Jiang N, Zhang W, Presting G, Windsor A, Navajas-Perez R, Torres MJ, Feltus FA, Porter B, Li Y, Burroughs AM, Luo M-C, Liu L, Christopher DA, Mount SM, Moore PH, Sugimura T, Jiang J, Schuler MA, Friedman V, Mitchell-Olds T, Shippen DE, dePamphilis CW, Palmer JD, Freeling M, Paterson AH, Gonsalves D, Wang L, Alam M (2008) The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature 452:991–996</mixed-citation></ref><ref id="CR30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nei</surname><given-names>M</given-names></name><name><surname>Rooney</surname><given-names>AP</given-names></name></person-group><article-title>Concerted and birth-and-death evolution of multigene families</article-title><source>Annu Rev Genet</source><year>2005</year><volume>39</volume><fpage>121</fpage><lpage>152</lpage><pub-id pub-id-type="doi">10.1146/annurev.genet.39.073003.112240</pub-id><pub-id pub-id-type="pmid">16285855</pub-id></mixed-citation></ref><ref id="CR31"><mixed-citation publication-type="other">Nylander JAA (2004) MrModelTest v2. In: Program distributed by the author, Evolutionary Biology Centre, Uppsala University</mixed-citation></ref><ref id="CR32"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Ohno</surname><given-names>S</given-names></name></person-group><source>Evolution by gene duplication</source><year>1970</year><publisher-loc>Heidelberg, Germany</publisher-loc><publisher-name>Springer-Verlag</publisher-name></mixed-citation></ref><ref id="CR33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ohta</surname><given-names>T</given-names></name></person-group><article-title>Time for acquiring a new gene by duplication</article-title><source>Proc Natl Acad Sci USA</source><year>1988</year><volume>85</volume><fpage>3509</fpage><lpage>3512</lpage><pub-id pub-id-type="doi">10.1073/pnas.85.10.3509</pub-id><pub-id pub-id-type="pmid">3368461</pub-id></mixed-citation></ref><ref id="CR34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Posada</surname><given-names>D</given-names></name><name><surname>Crandall</surname><given-names>KA</given-names></name></person-group><article-title>MODELTEST: testing the model of DNA substitution</article-title><source>Bioinformatics</source><year>1998</year><volume>14</volume><fpage>817</fpage><lpage>818</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/14.9.817</pub-id><pub-id pub-id-type="pmid">9918953</pub-id></mixed-citation></ref><ref id="CR35"><mixed-citation publication-type="other">Rambaut A (2006–2008) Tree Figure Drawing Tool Version 1.1.2. In: Institute of evolutionary biology, University of Edinburg</mixed-citation></ref><ref id="CR36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Richard</surname><given-names>G-F</given-names></name><name><surname>Dujon</surname><given-names>B</given-names></name></person-group><article-title>Molecular evolution of minisatellites in hemiascomycetous yeasts</article-title><source>Mol Biol Evol</source><year>2006</year><volume>23</volume><fpage>189</fpage><lpage>202</lpage><pub-id pub-id-type="doi">10.1093/molbev/msj022</pub-id><pub-id pub-id-type="pmid">16177231</pub-id></mixed-citation></ref><ref id="CR37"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ronquist</surname><given-names>F</given-names></name><name><surname>Huelsenbeck</surname><given-names>JP</given-names></name></person-group><article-title>MrBayes 3: Bayesian phylogenetic inference under mixed models</article-title><source>Bioinformatics</source><year>2003</year><volume>19</volume><fpage>1572</fpage><lpage>1574</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btg180</pub-id><pub-id pub-id-type="pmid">12912839</pub-id></mixed-citation></ref><ref id="CR38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schuermann</surname><given-names>D</given-names></name><name><surname>Molinier</surname><given-names>J</given-names></name><name><surname>Fritsch</surname><given-names>O</given-names></name><name><surname>Hohn</surname><given-names>B</given-names></name></person-group><article-title>The dual nature of homologous recombination in plants</article-title><source>Trends Genet</source><year>2005</year><volume>21</volume><fpage>172</fpage><lpage>181</lpage><pub-id pub-id-type="doi">10.1016/j.tig.2005.01.002</pub-id><pub-id pub-id-type="pmid">15734576</pub-id></mixed-citation></ref><ref id="CR39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schwartz</surname><given-names>S</given-names></name><name><surname>Kent</surname><given-names>WJ</given-names></name><name><surname>Smit</surname><given-names>A</given-names></name><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Baertsch</surname><given-names>R</given-names></name><name><surname>Hardison</surname><given-names>RC</given-names></name><name><surname>Haussler</surname><given-names>D</given-names></name><name><surname>Miller</surname><given-names>W</given-names></name></person-group><article-title>Human-mouse alignments with BLASTZ</article-title><source>Genome Res</source><year>2003</year><volume>13</volume><fpage>103</fpage><lpage>107</lpage><pub-id pub-id-type="doi">10.1101/gr.809403</pub-id><pub-id pub-id-type="pmid">12529312</pub-id></mixed-citation></ref><ref id="CR40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sémon</surname><given-names>M</given-names></name><name><surname>Wolfe</surname><given-names>KH</given-names></name></person-group><article-title>Consequences of genome duplication</article-title><source>Curr Opin Genet Dev</source><year>2007</year><volume>17</volume><fpage>505</fpage><lpage>512</lpage><pub-id pub-id-type="doi">10.1016/j.gde.2007.09.007</pub-id><pub-id pub-id-type="pmid">18006297</pub-id></mixed-citation></ref><ref id="CR41"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>GP</given-names></name></person-group><article-title>Evolution of repeated DNA sequences by unequal crossover</article-title><source>Science (New York, NY)</source><year>1976</year><volume>191</volume><fpage>528</fpage><lpage>535</lpage></mixed-citation></ref><ref id="CR42"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stange</surname><given-names>N</given-names></name><name><surname>Beier</surname><given-names>D</given-names></name><name><surname>Beier</surname><given-names>H</given-names></name></person-group><article-title>Expression of nuclear tRNA(Tyr) genes from <italic>Arabidopsis thaliana</italic> in HeLa cell and wheat germ extracts</article-title><source>Plant Mol Biol</source><year>1991</year><volume>16</volume><fpage>865</fpage><lpage>875</lpage><pub-id pub-id-type="doi">10.1007/BF00015078</pub-id><pub-id pub-id-type="pmid">1859869</pub-id></mixed-citation></ref><ref id="CR43"><mixed-citation publication-type="other">Swofford D (2003) PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4b10. In: Sinauer Associates, Sunderland, Massachusetts</mixed-citation></ref><ref id="CR44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tartof</surname><given-names>KD</given-names></name></person-group><article-title>Unequal crossing over then and now</article-title><source>Genetics</source><year>1988</year><volume>120</volume><fpage>1</fpage><lpage>6</lpage><pub-id pub-id-type="pmid">3065137</pub-id></mixed-citation></ref><ref id="CR45"><mixed-citation publication-type="other">Team RDC (2007) R: a language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria</mixed-citation></ref><ref id="CR46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thompson</surname><given-names>JD</given-names></name><name><surname>Higgins</surname><given-names>DG</given-names></name><name><surname>Gibson</surname><given-names>TJ</given-names></name></person-group><article-title>CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice</article-title><source>Nucl Acids Res</source><year>1994</year><volume>22</volume><fpage>4673</fpage><lpage>4680</lpage><pub-id pub-id-type="doi">10.1093/nar/22.22.4673</pub-id><pub-id pub-id-type="pmid">7984417</pub-id></mixed-citation></ref><ref id="CR47"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>YW</given-names></name><name><surname>Lai</surname><given-names>KN</given-names></name><name><surname>Tai</surname><given-names>PY</given-names></name><name><surname>Li</surname><given-names>WH</given-names></name></person-group><article-title>Rates of nucleotide substitution in angiosperm mitochondrial DNA sequences and dates of divergence between Brassica and other angiosperm lineages</article-title><source>J Mol Evol</source><year>1999</year><volume>48</volume><fpage>597</fpage><lpage>604</lpage><pub-id pub-id-type="doi">10.1007/PL00006502</pub-id><pub-id pub-id-type="pmid">10198125</pub-id></mixed-citation></ref><ref id="CR48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>J</given-names></name></person-group><article-title>Evolution by gene duplication: an update</article-title><source>Trends Ecol Evol</source><year>2003</year><volume>18</volume><fpage>292</fpage><lpage>298</lpage><pub-id pub-id-type="doi">10.1016/S0169-5347(03)00033-8</pub-id></mixed-citation></ref><ref id="CR49"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zimmer</surname><given-names>EA</given-names></name><name><surname>Martin</surname><given-names>SL</given-names></name><name><surname>Beverley</surname><given-names>SM</given-names></name><name><surname>Kan</surname><given-names>YW</given-names></name><name><surname>Wilson</surname><given-names>AC</given-names></name></person-group><article-title>Rapid duplication and loss of genes coding for the a chains of hemoglobin</article-title><source>PNAS</source><year>1980</year><volume>77</volume><fpage>2158</fpage><lpage>2162</lpage><pub-id pub-id-type="doi">10.1073/pnas.77.4.2158</pub-id><pub-id pub-id-type="pmid">6929543</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/CyberinfraV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000185 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000185 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Ticri/CIDE
|area= CyberinfraV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.25. |  |