Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation</title>
<author><name sortKey="Li, Feng" sort="Li, Feng" uniqKey="Li F" first="Feng" last="Li">Feng Li</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Upadhyaya, Narayana M" sort="Upadhyaya, Narayana M" uniqKey="Upadhyaya N" first="Narayana M." last="Upadhyaya">Narayana M. Upadhyaya</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sperschneider, Jana" sort="Sperschneider, Jana" uniqKey="Sperschneider J" first="Jana" last="Sperschneider">Jana Sperschneider</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Biological Data Science Institute,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Matny, Oadi" sort="Matny, Oadi" uniqKey="Matny O" first="Oadi" last="Matny">Oadi Matny</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nguyen Phuc, Hoa" sort="Nguyen Phuc, Hoa" uniqKey="Nguyen Phuc H" first="Hoa" last="Nguyen-Phuc">Hoa Nguyen-Phuc</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mago, Rohit" sort="Mago, Rohit" uniqKey="Mago R" first="Rohit" last="Mago">Rohit Mago</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Raley, Castle" sort="Raley, Castle" uniqKey="Raley C" first="Castle" last="Raley">Castle Raley</name>
<affiliation><nlm:aff id="Aff4">Leidos Biomedical Research, Frederick, MD 21702 USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff8"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 9510</institution-id>
<institution-id institution-id-type="GRID">grid.253615.6</institution-id>
<institution>Present Address: The George Washington University,</institution>
</institution-wrap>
Washington, DC 20052 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miller, Marisa E" sort="Miller, Marisa E" uniqKey="Miller M" first="Marisa E." last="Miller">Marisa E. Miller</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff9">Present Address: Pairwise, Durham, NC 27709 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Silverstein, Kevin A T" sort="Silverstein, Kevin A T" uniqKey="Silverstein K" first="Kevin A. T." last="Silverstein">Kevin A. T. Silverstein</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Minnesota Supercomputing Institute,</institution>
</institution-wrap>
Minneapolis, MN 55455 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Henningsen, Eva" sort="Henningsen, Eva" uniqKey="Henningsen E" first="Eva" last="Henningsen">Eva Henningsen</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hirsch, Cory D" sort="Hirsch, Cory D" uniqKey="Hirsch C" first="Cory D." last="Hirsch">Cory D. Hirsch</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Visser, Botma" sort="Visser, Botma" uniqKey="Visser B" first="Botma" last="Visser">Botma Visser</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2284 638X</institution-id>
<institution-id institution-id-type="GRID">grid.412219.d</institution-id>
<institution>University of the Free State,</institution>
</institution-wrap>
Bloemfontein, 9301 Free State South Africa</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Pretorius, Zacharias A" sort="Pretorius, Zacharias A" uniqKey="Pretorius Z" first="Zacharias A." last="Pretorius">Zacharias A. Pretorius</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2284 638X</institution-id>
<institution-id institution-id-type="GRID">grid.412219.d</institution-id>
<institution>University of the Free State,</institution>
</institution-wrap>
Bloemfontein, 9301 Free State South Africa</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Steffenson, Brian J" sort="Steffenson, Brian J" uniqKey="Steffenson B" first="Brian J." last="Steffenson">Brian J. Steffenson</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Schwessinger, Benjamin" sort="Schwessinger, Benjamin" uniqKey="Schwessinger B" first="Benjamin" last="Schwessinger">Benjamin Schwessinger</name>
<affiliation><nlm:aff id="Aff7"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Research School of Biology,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Dodds, Peter N" sort="Dodds, Peter N" uniqKey="Dodds P" first="Peter N." last="Dodds">Peter N. Dodds</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Figueroa, Melania" sort="Figueroa, Melania" uniqKey="Figueroa M" first="Melania" last="Figueroa">Melania Figueroa</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">31699975</idno>
<idno type="pmc">6838127</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC6838127</idno>
<idno type="RBID">PMC:6838127</idno>
<idno type="doi">10.1038/s41467-019-12927-7</idno>
<date when="2019">2019</date>
<idno type="wicri:Area/Pmc/Corpus">000208</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000208</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation</title>
<author><name sortKey="Li, Feng" sort="Li, Feng" uniqKey="Li F" first="Feng" last="Li">Feng Li</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Upadhyaya, Narayana M" sort="Upadhyaya, Narayana M" uniqKey="Upadhyaya N" first="Narayana M." last="Upadhyaya">Narayana M. Upadhyaya</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sperschneider, Jana" sort="Sperschneider, Jana" uniqKey="Sperschneider J" first="Jana" last="Sperschneider">Jana Sperschneider</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Biological Data Science Institute,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Matny, Oadi" sort="Matny, Oadi" uniqKey="Matny O" first="Oadi" last="Matny">Oadi Matny</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Nguyen Phuc, Hoa" sort="Nguyen Phuc, Hoa" uniqKey="Nguyen Phuc H" first="Hoa" last="Nguyen-Phuc">Hoa Nguyen-Phuc</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mago, Rohit" sort="Mago, Rohit" uniqKey="Mago R" first="Rohit" last="Mago">Rohit Mago</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Raley, Castle" sort="Raley, Castle" uniqKey="Raley C" first="Castle" last="Raley">Castle Raley</name>
<affiliation><nlm:aff id="Aff4">Leidos Biomedical Research, Frederick, MD 21702 USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff8"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 9510</institution-id>
<institution-id institution-id-type="GRID">grid.253615.6</institution-id>
<institution>Present Address: The George Washington University,</institution>
</institution-wrap>
Washington, DC 20052 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miller, Marisa E" sort="Miller, Marisa E" uniqKey="Miller M" first="Marisa E." last="Miller">Marisa E. Miller</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff9">Present Address: Pairwise, Durham, NC 27709 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Silverstein, Kevin A T" sort="Silverstein, Kevin A T" uniqKey="Silverstein K" first="Kevin A. T." last="Silverstein">Kevin A. T. Silverstein</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Minnesota Supercomputing Institute,</institution>
</institution-wrap>
Minneapolis, MN 55455 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Henningsen, Eva" sort="Henningsen, Eva" uniqKey="Henningsen E" first="Eva" last="Henningsen">Eva Henningsen</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hirsch, Cory D" sort="Hirsch, Cory D" uniqKey="Hirsch C" first="Cory D." last="Hirsch">Cory D. Hirsch</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Visser, Botma" sort="Visser, Botma" uniqKey="Visser B" first="Botma" last="Visser">Botma Visser</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2284 638X</institution-id>
<institution-id institution-id-type="GRID">grid.412219.d</institution-id>
<institution>University of the Free State,</institution>
</institution-wrap>
Bloemfontein, 9301 Free State South Africa</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Pretorius, Zacharias A" sort="Pretorius, Zacharias A" uniqKey="Pretorius Z" first="Zacharias A." last="Pretorius">Zacharias A. Pretorius</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2284 638X</institution-id>
<institution-id institution-id-type="GRID">grid.412219.d</institution-id>
<institution>University of the Free State,</institution>
</institution-wrap>
Bloemfontein, 9301 Free State South Africa</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Steffenson, Brian J" sort="Steffenson, Brian J" uniqKey="Steffenson B" first="Brian J." last="Steffenson">Brian J. Steffenson</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Schwessinger, Benjamin" sort="Schwessinger, Benjamin" uniqKey="Schwessinger B" first="Benjamin" last="Schwessinger">Benjamin Schwessinger</name>
<affiliation><nlm:aff id="Aff7"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Research School of Biology,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Dodds, Peter N" sort="Dodds, Peter N" uniqKey="Dodds P" first="Peter N." last="Dodds">Peter N. Dodds</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Figueroa, Melania" sort="Figueroa, Melania" uniqKey="Figueroa M" first="Melania" last="Figueroa">Melania Figueroa</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Nature Communications</title>
<idno type="eISSN">2041-1723</idno>
<imprint><date when="2019">2019</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p id="Par1">Parasexuality contributes to diversity and adaptive evolution of haploid (monokaryotic) fungi. However, non-sexual genetic exchange mechanisms are not defined in dikaryotic fungi (containing two distinct haploid nuclei). Newly emerged strains of the wheat stem rust pathogen, <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
(<italic>Pgt</italic>
), such as Ug99, are a major threat to global food security. Here, we provide genomics-based evidence supporting that Ug99 arose by somatic hybridisation and nuclear exchange between dikaryons. Fully haplotype-resolved genome assembly and DNA proximity analysis reveal that Ug99 shares one haploid nucleus genotype with a much older African lineage of <italic>Pgt</italic>
, with no recombination or chromosome reassortment. These findings indicate that nuclear exchange between dikaryotes can generate genetic diversity and facilitate the emergence of new lineages in asexual fungal populations.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Stajich, Je" uniqKey="Stajich J">JE Stajich</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Spatafora, Jw" uniqKey="Spatafora J">JW Spatafora</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aime, Mc" uniqKey="Aime M">MC Aime</name>
</author>
<author><name sortKey="Mctaggart, Ar" uniqKey="Mctaggart A">AR McTaggart</name>
</author>
<author><name sortKey="Mondo, Sj" uniqKey="Mondo S">SJ Mondo</name>
</author>
<author><name sortKey="Duplessis, S" uniqKey="Duplessis S">S Duplessis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watson, I" uniqKey="Watson I">I Watson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ellingboe, A" uniqKey="Ellingboe A">A Ellingboe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Flor, H" uniqKey="Flor H">H Flor</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bartos, P" uniqKey="Bartos P">P Bartos</name>
</author>
<author><name sortKey="Fleischmann, G" uniqKey="Fleischmann G">G Fleischmann</name>
</author>
<author><name sortKey="Samborski, D" uniqKey="Samborski D">D Samborski</name>
</author>
<author><name sortKey="Shipton, W" uniqKey="Shipton W">W Shipton</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burdon, Jj" uniqKey="Burdon J">JJ Burdon</name>
</author>
<author><name sortKey="Marshall, Dr" uniqKey="Marshall D">DR Marshall</name>
</author>
<author><name sortKey="Luig, Nh" uniqKey="Luig N">NH Luig</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, R" uniqKey="Park R">R Park</name>
</author>
<author><name sortKey="Burdon, J" uniqKey="Burdon J">J Burdon</name>
</author>
<author><name sortKey="Jahoor, A" uniqKey="Jahoor A">A Jahoor</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, Rf" uniqKey="Park R">RF Park</name>
</author>
<author><name sortKey="Wellings, Cr" uniqKey="Wellings C">CR Wellings</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, Me" uniqKey="Miller M">ME Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schwessinger, B" uniqKey="Schwessinger B">B Schwessinger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Singh, Rp" uniqKey="Singh R">RP Singh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pretorius, Z" uniqKey="Pretorius Z">Z Pretorius</name>
</author>
<author><name sortKey="Singh, R" uniqKey="Singh R">R Singh</name>
</author>
<author><name sortKey="Wagoire, W" uniqKey="Wagoire W">W Wagoire</name>
</author>
<author><name sortKey="Payne, T" uniqKey="Payne T">T Payne</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Visser, B" uniqKey="Visser B">B Visser</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Olivera, P" uniqKey="Olivera P">P Olivera</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Visser, B" uniqKey="Visser B">B Visser</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pretorius, Za" uniqKey="Pretorius Z">ZA Pretorius</name>
</author>
<author><name sortKey="Pakendorf, Kw" uniqKey="Pakendorf K">KW Pakendorf</name>
</author>
<author><name sortKey="Marais, Gf" uniqKey="Marais G">GF Marais</name>
</author>
<author><name sortKey="Prins, R" uniqKey="Prins R">R Prins</name>
</author>
<author><name sortKey="Komen, Js" uniqKey="Komen J">JS Komen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, Rf" uniqKey="Park R">RF Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Upadhyaya, Nm" uniqKey="Upadhyaya N">NM Upadhyaya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Duplessis, S" uniqKey="Duplessis S">S Duplessis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, J" uniqKey="Chen J">J Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Salcedo, A" uniqKey="Salcedo A">A Salcedo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, C" uniqKey="Anderson C">C Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Belton, Jm" uniqKey="Belton J">JM Belton</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Boehm, Ewa" uniqKey="Boehm E">EWA Boehm</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cuomo, Ca" uniqKey="Cuomo C">CA Cuomo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lewis, Cm" uniqKey="Lewis C">CM Lewis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Figueroa, M" uniqKey="Figueroa M">M Figueroa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lei, Y" uniqKey="Lei Y">Y Lei</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Barr, R" uniqKey="Barr R">R Barr</name>
</author>
<author><name sortKey="Caldwell, R" uniqKey="Caldwell R">R Caldwell</name>
</author>
<author><name sortKey="Amacher, R" uniqKey="Amacher R">R Amacher</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name>
</author>
<author><name sortKey="Mccallum, B" uniqKey="Mccallum B">B McCallum</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kruzel, Ek" uniqKey="Kruzel E">EK Kruzel</name>
</author>
<author><name sortKey="Hull, Cm" uniqKey="Hull C">CM Hull</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Clark, Ta" uniqKey="Clark T">TA Clark</name>
</author>
<author><name sortKey="Anderson, Jb" uniqKey="Anderson J">JB Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Corradi, N" uniqKey="Corradi N">N Corradi</name>
</author>
<author><name sortKey="Brachmann, A" uniqKey="Brachmann A">A Brachmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pretorius, Z" uniqKey="Pretorius Z">Z Pretorius</name>
</author>
<author><name sortKey="Szabo, L" uniqKey="Szabo L">L Szabo</name>
</author>
<author><name sortKey="Boshoff, W" uniqKey="Boshoff W">W Boshoff</name>
</author>
<author><name sortKey="Herselman, L" uniqKey="Herselman L">L Herselman</name>
</author>
<author><name sortKey="Visser, B" uniqKey="Visser B">B Visser</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Koren, S" uniqKey="Koren S">S Koren</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Walker, Bj" uniqKey="Walker B">BJ Walker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
<author><name sortKey="Durbin, R" uniqKey="Durbin R">R Durbin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Quinlan, Ar" uniqKey="Quinlan A">AR Quinlan</name>
</author>
<author><name sortKey="Hall, Im" uniqKey="Hall I">IM Hall</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gurevich, A" uniqKey="Gurevich A">A Gurevich</name>
</author>
<author><name sortKey="Saveliev, V" uniqKey="Saveliev V">V Saveliev</name>
</author>
<author><name sortKey="Vyahhi, N" uniqKey="Vyahhi N">N Vyahhi</name>
</author>
<author><name sortKey="Tesler, G" uniqKey="Tesler G">G Tesler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stanke, M" uniqKey="Stanke M">M Stanke</name>
</author>
<author><name sortKey="Morgenstern, B" uniqKey="Morgenstern B">B Morgenstern</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Simao, Fa" uniqKey="Simao F">FA Simão</name>
</author>
<author><name sortKey="Waterhouse, Rm" uniqKey="Waterhouse R">RM Waterhouse</name>
</author>
<author><name sortKey="Ioannidis, P" uniqKey="Ioannidis P">P Ioannidis</name>
</author>
<author><name sortKey="Kriventseva, Ev" uniqKey="Kriventseva E">EV Kriventseva</name>
</author>
<author><name sortKey="Zdobnov, Em" uniqKey="Zdobnov E">EM Zdobnov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cabanettes, F" uniqKey="Cabanettes F">F Cabanettes</name>
</author>
<author><name sortKey="Klopp, C" uniqKey="Klopp C">C Klopp</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lawrence, M" uniqKey="Lawrence M">M Lawrence</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yin, T" uniqKey="Yin T">T Yin</name>
</author>
<author><name sortKey="Cook, D" uniqKey="Cook D">D Cook</name>
</author>
<author><name sortKey="Lawrence, M" uniqKey="Lawrence M">M Lawrence</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kurtz, S" uniqKey="Kurtz S">S Kurtz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Garrison, E" uniqKey="Garrison E">E Garrison</name>
</author>
<author><name sortKey="Marth, G" uniqKey="Marth G">G Marth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wingett, Steven W" uniqKey="Wingett S">Steven W. Wingett</name>
</author>
<author><name sortKey="Ewels, Philip" uniqKey="Ewels P">Philip Ewels</name>
</author>
<author><name sortKey="Furlan Magaril, Mayra" uniqKey="Furlan Magaril M">Mayra Furlan-Magaril</name>
</author>
<author><name sortKey="Nagano, Takashi" uniqKey="Nagano T">Takashi Nagano</name>
</author>
<author><name sortKey="Schoenfelder, Stefan" uniqKey="Schoenfelder S">Stefan Schoenfelder</name>
</author>
<author><name sortKey="Fraser, Peter" uniqKey="Fraser P">Peter Fraser</name>
</author>
<author><name sortKey="Andrews, Simon" uniqKey="Andrews S">Simon Andrews</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Langmead, B" uniqKey="Langmead B">B Langmead</name>
</author>
<author><name sortKey="Salzberg, Sl" uniqKey="Salzberg S">SL Salzberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghurye, J" uniqKey="Ghurye J">J Ghurye</name>
</author>
<author><name sortKey="Pop, M" uniqKey="Pop M">M Pop</name>
</author>
<author><name sortKey="Koren, S" uniqKey="Koren S">S Koren</name>
</author>
<author><name sortKey="Bickhart, D" uniqKey="Bickhart D">D Bickhart</name>
</author>
<author><name sortKey="Chin, C S" uniqKey="Chin C">C-S Chin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bao, W" uniqKey="Bao W">W Bao</name>
</author>
<author><name sortKey="Kojima, Kk" uniqKey="Kojima K">KK Kojima</name>
</author>
<author><name sortKey="Kohany, O" uniqKey="Kohany O">O Kohany</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, D" uniqKey="Kim D">D Kim</name>
</author>
<author><name sortKey="Langmead, B" uniqKey="Langmead B">B Langmead</name>
</author>
<author><name sortKey="Salzberg, Sl" uniqKey="Salzberg S">SL Salzberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pertea, M" uniqKey="Pertea M">M Pertea</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grabherr, Mg" uniqKey="Grabherr M">MG Grabherr</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Haas, Bj" uniqKey="Haas B">BJ Haas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Testa, Ac" uniqKey="Testa A">AC Testa</name>
</author>
<author><name sortKey="Hane, Jk" uniqKey="Hane J">JK Hane</name>
</author>
<author><name sortKey="Ellwood, Sr" uniqKey="Ellwood S">SR Ellwood</name>
</author>
<author><name sortKey="Oliver, Rp" uniqKey="Oliver R">RP Oliver</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Slater, Gsc" uniqKey="Slater G">GSC Slater</name>
</author>
<author><name sortKey="Birney, E" uniqKey="Birney E">E Birney</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lomsadze, A" uniqKey="Lomsadze A">A Lomsadze</name>
</author>
<author><name sortKey="Ter Hovhannisyan, V" uniqKey="Ter Hovhannisyan V">V Ter-Hovhannisyan</name>
</author>
<author><name sortKey="Chernoff, Yo" uniqKey="Chernoff Y">YO Chernoff</name>
</author>
<author><name sortKey="Borodovsky, M" uniqKey="Borodovsky M">M Borodovsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lowe, Tm" uniqKey="Lowe T">TM Lowe</name>
</author>
<author><name sortKey="Chan, Pp" uniqKey="Chan P">PP Chan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jones, P" uniqKey="Jones P">P Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yin, Y" uniqKey="Yin Y">Y Yin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buchfink, B" uniqKey="Buchfink B">B Buchfink</name>
</author>
<author><name sortKey="Xie, C" uniqKey="Xie C">C Xie</name>
</author>
<author><name sortKey="Huson, Dh" uniqKey="Huson D">DH Huson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Consortium, U" uniqKey="Consortium U">U Consortium</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rawlings, Nd" uniqKey="Rawlings N">ND Rawlings</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huerta Cepas, J" uniqKey="Huerta Cepas J">J Huerta-Cepas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chojnacki, S" uniqKey="Chojnacki S">S Chojnacki</name>
</author>
<author><name sortKey="Cowley, A" uniqKey="Cowley A">A Cowley</name>
</author>
<author><name sortKey="Lee, J" uniqKey="Lee J">J Lee</name>
</author>
<author><name sortKey="Foix, A" uniqKey="Foix A">A Foix</name>
</author>
<author><name sortKey="Lopez, R" uniqKey="Lopez R">R Lopez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gel, B" uniqKey="Gel B">B Gel</name>
</author>
<author><name sortKey="Serra, E" uniqKey="Serra E">E Serra</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bendtsen, Jd" uniqKey="Bendtsen J">JD Bendtsen</name>
</author>
<author><name sortKey="Nielsen, H" uniqKey="Nielsen H">H Nielsen</name>
</author>
<author><name sortKey="Von Heijne, G" uniqKey="Von Heijne G">G von Heijne</name>
</author>
<author><name sortKey="Brunak, S" uniqKey="Brunak S">S Brunak</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Krogh, A" uniqKey="Krogh A">A Krogh</name>
</author>
<author><name sortKey="Larsson, B" uniqKey="Larsson B">B Larsson</name>
</author>
<author><name sortKey="Von Heijne, G" uniqKey="Von Heijne G">G Von Heijne</name>
</author>
<author><name sortKey="Sonnhammer, El" uniqKey="Sonnhammer E">EL Sonnhammer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Emms, Dm" uniqKey="Emms D">DM Emms</name>
</author>
<author><name sortKey="Kelly, S" uniqKey="Kelly S">S Kelly</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bolger, Am" uniqKey="Bolger A">AM Bolger</name>
</author>
<author><name sortKey="Lohse, M" uniqKey="Lohse M">M Lohse</name>
</author>
<author><name sortKey="Usadel, B" uniqKey="Usadel B">B Usadel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Knaus, Bj" uniqKey="Knaus B">BJ Knaus</name>
</author>
<author><name sortKey="Grunwald, Nj" uniqKey="Grunwald N">NJ Grünwald</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A Stamatakis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Paradis, E" uniqKey="Paradis E">E Paradis</name>
</author>
<author><name sortKey="Claude, J" uniqKey="Claude J">J Claude</name>
</author>
<author><name sortKey="Strimmer, K" uniqKey="Strimmer K">K Strimmer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schliep, Kp" uniqKey="Schliep K">KP Schliep</name>
</author>
</analytic>
</biblStruct>
<biblStruct></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">Nat Commun</journal-id>
<journal-id journal-id-type="iso-abbrev">Nat Commun</journal-id>
<journal-title-group><journal-title>Nature Communications</journal-title>
</journal-title-group>
<issn pub-type="epub">2041-1723</issn>
<publisher><publisher-name>Nature Publishing Group UK</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">31699975</article-id>
<article-id pub-id-type="pmc">6838127</article-id>
<article-id pub-id-type="publisher-id">12927</article-id>
<article-id pub-id-type="doi">10.1038/s41467-019-12927-7</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Emergence of the Ug99 lineage of the wheat stem rust pathogen through somatic hybridisation</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" equal-contrib="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-8528-4249</contrib-id>
<name><surname>Li</surname>
<given-names>Feng</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-3052-0416</contrib-id>
<name><surname>Upadhyaya</surname>
<given-names>Narayana M.</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Sperschneider</surname>
<given-names>Jana</given-names>
</name>
<xref ref-type="aff" rid="Aff3">3</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-8447-2886</contrib-id>
<name><surname>Matny</surname>
<given-names>Oadi</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-9329-4287</contrib-id>
<name><surname>Nguyen-Phuc</surname>
<given-names>Hoa</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Mago</surname>
<given-names>Rohit</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Raley</surname>
<given-names>Castle</given-names>
</name>
<xref ref-type="aff" rid="Aff4">4</xref>
<xref ref-type="aff" rid="Aff8">8</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Miller</surname>
<given-names>Marisa E.</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
<xref ref-type="aff" rid="Aff9">9</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-4955-3218</contrib-id>
<name><surname>Silverstein</surname>
<given-names>Kevin A. T.</given-names>
</name>
<xref ref-type="aff" rid="Aff5">5</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-9619-3705</contrib-id>
<name><surname>Henningsen</surname>
<given-names>Eva</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-3409-758X</contrib-id>
<name><surname>Hirsch</surname>
<given-names>Cory D.</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Visser</surname>
<given-names>Botma</given-names>
</name>
<xref ref-type="aff" rid="Aff6">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Pretorius</surname>
<given-names>Zacharias A.</given-names>
</name>
<xref ref-type="aff" rid="Aff6">6</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-7961-5363</contrib-id>
<name><surname>Steffenson</surname>
<given-names>Brian J.</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Schwessinger</surname>
<given-names>Benjamin</given-names>
</name>
<xref ref-type="aff" rid="Aff7">7</xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0003-0620-5923</contrib-id>
<name><surname>Dodds</surname>
<given-names>Peter N.</given-names>
</name>
<address><email>peter.dodds@csiro.au</email>
</address>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0003-2636-661X</contrib-id>
<name><surname>Figueroa</surname>
<given-names>Melania</given-names>
</name>
<address><email>melania.figueroa@csiro.au</email>
</address>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Department of Plant Pathology,</institution>
<institution>University of Minnesota,</institution>
</institution-wrap>
St. Paul, MN 55108 USA</aff>
<aff id="Aff2"><label>2</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.493032.f</institution-id>
<institution>Commonwealth Scientific and Industrial Research Organisation,</institution>
<institution>Agriculture and Food,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</aff>
<aff id="Aff3"><label>3</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Biological Data Science Institute,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</aff>
<aff id="Aff4"><label>4</label>
Leidos Biomedical Research, Frederick, MD 21702 USA</aff>
<aff id="Aff5"><label>5</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000000419368657</institution-id>
<institution-id institution-id-type="GRID">grid.17635.36</institution-id>
<institution>Minnesota Supercomputing Institute,</institution>
</institution-wrap>
Minneapolis, MN 55455 USA</aff>
<aff id="Aff6"><label>6</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2284 638X</institution-id>
<institution-id institution-id-type="GRID">grid.412219.d</institution-id>
<institution>University of the Free State,</institution>
</institution-wrap>
Bloemfontein, 9301 Free State South Africa</aff>
<aff id="Aff7"><label>7</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2180 7477</institution-id>
<institution-id institution-id-type="GRID">grid.1001.0</institution-id>
<institution>Research School of Biology,</institution>
<institution>The Australian National University,</institution>
</institution-wrap>
Canberra, ACT 2601 Australia</aff>
<aff id="Aff8"><label>8</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 9510</institution-id>
<institution-id institution-id-type="GRID">grid.253615.6</institution-id>
<institution>Present Address: The George Washington University,</institution>
</institution-wrap>
Washington, DC 20052 USA</aff>
<aff id="Aff9"><label>9</label>
Present Address: Pairwise, Durham, NC 27709 USA</aff>
</contrib-group>
<pub-date pub-type="epub"><day>7</day>
<month>11</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>7</day>
<month>11</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection"><year>2019</year>
</pub-date>
<volume>10</volume>
<elocation-id>5068</elocation-id>
<history><date date-type="received"><day>27</day>
<month>6</month>
<year>2019</year>
</date>
<date date-type="accepted"><day>3</day>
<month>10</month>
<year>2019</year>
</date>
</history>
<permissions><copyright-statement>© Crown 2019</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
.</license-p>
</license>
</permissions>
<abstract id="Abs1"><p id="Par1">Parasexuality contributes to diversity and adaptive evolution of haploid (monokaryotic) fungi. However, non-sexual genetic exchange mechanisms are not defined in dikaryotic fungi (containing two distinct haploid nuclei). Newly emerged strains of the wheat stem rust pathogen, <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
(<italic>Pgt</italic>
), such as Ug99, are a major threat to global food security. Here, we provide genomics-based evidence supporting that Ug99 arose by somatic hybridisation and nuclear exchange between dikaryons. Fully haplotype-resolved genome assembly and DNA proximity analysis reveal that Ug99 shares one haploid nucleus genotype with a much older African lineage of <italic>Pgt</italic>
, with no recombination or chromosome reassortment. These findings indicate that nuclear exchange between dikaryotes can generate genetic diversity and facilitate the emergence of new lineages in asexual fungal populations.</p>
</abstract>
<abstract id="Abs2" abstract-type="web-summary"><p id="Par2">Strain Ug99 of the wheat stem rust pathogen, <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
, is a major threat to global food security. Here, the authors present genomic analyses supporting that Ug99 arose as a result of non-sexual genetic exchange between dikaryotic ancestors.</p>
</abstract>
<kwd-group kwd-group-type="npg-subject"><title>Subject terms</title>
<kwd>Genome evolution</kwd>
<kwd>Fungal genetics</kwd>
<kwd>Pathogens</kwd>
</kwd-group>
<funding-group><award-group><funding-source><institution>USDA-NIFA (2017-67012-26117)</institution>
</funding-source>
</award-group>
</funding-group>
<funding-group><award-group><funding-source><institution-wrap><institution-id institution-id-type="FundRef">https://doi.org/10.13039/501100000923</institution-id>
<institution>Department of Education and Training | Australian Research Council (ARC)</institution>
</institution-wrap>
</funding-source>
<award-id>FT180100024</award-id>
<principal-award-recipient><name><surname>Schwessinger</surname>
<given-names>Benjamin</given-names>
</name>
</principal-award-recipient>
</award-group>
</funding-group>
<funding-group><award-group><funding-source><institution>TwoBlades Foundation USDA-AFRI (2017-08221 USDA-NIFA Hatch/Figueroa project MIN-22-G19</institution>
</funding-source>
</award-group>
</funding-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2019</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="introduction"><title>Introduction</title>
<p id="Par3">Generation of genetic diversity is crucial for the evolution of new traits, with mutation and sexual recombination as the main drivers of diversity in most eukaryotes. However, many species in the fungal kingdom can propagate asexually for extended periods and therefore understanding alternative mechanisms contributing to genetic diversity in asexual populations has been of great interest<sup><xref ref-type="bibr" rid="CR1">1</xref>
,<xref ref-type="bibr" rid="CR2">2</xref>
</sup>
. Some fungi can use a parasexual mechanism to exchange genetic material independently of meiosis<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
. This process involves anastomosis of haploid hyphae and fusion of two nuclei to generate a single diploid nucleus, which subsequently undergoes progressive chromosome loss to generate recombinant haploid offspring. Parasexuality has been described in members of the ascomycete phylum (64% of described fungal species) in which the dominant asexually propagating form is haploid<sup><xref ref-type="bibr" rid="CR3">3</xref>
</sup>
. However, in basidiomycete fungi (34% of described species), the predominant life stage is generally dikaryotic, with two different haploid nuclei maintained within each individual<sup><xref ref-type="bibr" rid="CR3">3</xref>
</sup>
. The role of non-sexual genetic exchange between such dikaryons in generating genetic diversity is not known.</p>
<p id="Par4">Basidiomycetes include many fungi with critical ecosystem functions, such as wood decay and plant symbiosis, as well as agents of important human and plant diseases<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. Rust fungi (subphylum Pucciniomycotina) comprise over 8000 species including many pathogens of major agricultural and ecological significance<sup><xref ref-type="bibr" rid="CR4">4</xref>
</sup>
. These organisms are obligate parasites with complex life cycles that can include indefinite asexual reproduction through infectious dikaryotic urediniospores. Early researchers speculated that rust fungi can exchange genetic material during the asexual phase based on the isolation of new strains, after co-infection with two potential parental isolates, with novel virulence phenotypes<sup><xref ref-type="bibr" rid="CR5">5</xref>
–<xref ref-type="bibr" rid="CR8">8</xref>
</sup>
. However, these hypotheses could not be tested molecularly at the time. Some naturally occurring rust pathotypes have been suggested to have arisen by somatic hybridisation and genetic exchange based on limited molecular evidence of shared isozyme or random amplified polymorphic DNA (RAPD) markers<sup><xref ref-type="bibr" rid="CR9">9</xref>
,<xref ref-type="bibr" rid="CR10">10</xref>
</sup>
. Mechanisms underlying genetic exchange are unknown, but may involve hyphal anastomosis followed by nuclear exchange and/or nuclear fusion and recombination<sup><xref ref-type="bibr" rid="CR11">11</xref>
</sup>
. Recent advances in assembling complete karyon sequences in rust fungi<sup><xref ref-type="bibr" rid="CR12">12</xref>
,<xref ref-type="bibr" rid="CR13">13</xref>
</sup>
provide the opportunity to definitively detect and discriminate between nuclear exchange and recombination.</p>
<p id="Par5">The Ug99 strain (race TTKSK) of the wheat stem rust pathogen <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
(<italic>Pgt</italic>
) presents a significant threat to global wheat production<sup><xref ref-type="bibr" rid="CR14">14</xref>
</sup>
. It was first detected in Uganda in 1998 and described in 1999<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
, and has since given rise to an asexual lineage that has spread through Africa and the Middle East causing devastating epidemics<sup><xref ref-type="bibr" rid="CR14">14</xref>
</sup>
. The origin of the Ug99 lineage is unknown, although it is genetically distinct from other <italic>Pgt</italic>
races<sup><xref ref-type="bibr" rid="CR16">16</xref>
,<xref ref-type="bibr" rid="CR17">17</xref>
</sup>
. This indicates that Ug99 is not likely derived by mutation of longstanding stem rust asexual lineages such as the race 21 group, which has been predominant in southern Africa at least since the 1920’s and spread to Australia in the 1950’s<sup><xref ref-type="bibr" rid="CR18">18</xref>
–<xref ref-type="bibr" rid="CR20">20</xref>
</sup>
.</p>
<p id="Par6">Here, we generate haplotype-phased genome references for the original Ug99 isolate collected in Uganda<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
and an Australian <italic>Pgt</italic>
isolate of pathotype 21-0<sup><xref ref-type="bibr" rid="CR20">20</xref>
,<xref ref-type="bibr" rid="CR21">21</xref>
</sup>
. We show by genome comparison that Ug99 shares one haploid nucleus genotype with <italic>Pgt</italic>
21-0, with no recombination or chromosome reassortment. This indicates that Ug99 arose by somatic hybridisation and nuclear exchange between an African member of the <italic>Pgt</italic>
21 lineage and an unknown isolate. Thus, nuclear exchanges between dikaryotic fungi can contribute to the emergence of new variants with significant epidemiological impacts.</p>
</sec>
<sec id="Sec2" sec-type="results"><title>Results</title>
<sec id="Sec3"><title>Haplotype-phased genome assembly</title>
<p id="Par7">A single pustule derived from the original Ug99 isolate of <italic>Pgt</italic>
<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
was purified and its pathotype was confirmed using the standard wheat differential set (Supplementary Data <xref rid="MOESM4" ref-type="media">1</xref>
). We generated polished long-read genome assemblies for both Ug99 and the Australian stem rust isolate <italic>Pgt</italic>
21-0<sup><xref ref-type="bibr" rid="CR18">18</xref>
</sup>
using single-molecule real-time (SMRT) and Illumina sequence data (Supplementary Tables <xref rid="MOESM1" ref-type="media">1</xref>
and <xref rid="MOESM1" ref-type="media">2</xref>
). Genome assembly of <italic>Pgt</italic>
21-0 resulted in 410 contigs with a total size of 177 Mbp and a contig N<sub>50</sub>
of 1.26 Mbp (Supplementary Table <xref rid="MOESM1" ref-type="media">3</xref>
). Similarly, the size of the genome assembly of Ug99 was 176 Mbp represented in 514 contigs with contig N<sub>50</sub>
of 0.97 Mbp. Thus, both assemblies were twice the size of a collapsed haploid assembly previously generated for a North American <italic>Pgt</italic>
isolate (88 Mbp)<sup><xref ref-type="bibr" rid="CR22">22</xref>
</sup>
, suggesting that in each assembly the sequences of the two haploid karyons were fully represented independently. Consistent with this, both genome assemblies contained over 96% of conserved fungal BUSCO genes, with the majority present in two copies (Supplementary Table <xref rid="MOESM1" ref-type="media">3</xref>
). Furthermore, the <italic>Pgt</italic>
21-0 assembly contained 69 telomeres, out of a total of 72 expected for a dikaryotic genome given the known haploid chromosome number of eighteen<sup><xref ref-type="bibr" rid="CR23">23</xref>
</sup>
. To identify sequences representing alternate haplotypes within each assembly we developed a gene synteny approach to assign contigs to groups representing paired homologous haplotypes (Fig. <xref rid="Fig1" ref-type="fig">1</xref>
). Using this approach, homologous pairs of sequences from each haplotype were assigned to 44 bins in <italic>Pgt</italic>
21-0 and 62 bins in Ug99, which represented over 94% of each assembly (Supplementary Table <xref rid="MOESM1" ref-type="media">3</xref>
and Supplementary Data <xref rid="MOESM5" ref-type="media">2</xref>
). Three of the 18 chromosomes in <italic>Pgt</italic>
21-0, and two in Ug99 seemed to be fully assembled as these bins contained telomere sequences at each end.<fig id="Fig1"><label>Fig. 1</label>
<caption><p>Strategy to identify homologous contigs in genome assemblies by gene synteny. To detect shared content, <italic>Pgt</italic>
gene models<sup><xref ref-type="bibr" rid="CR21">21</xref>
</sup>
(grey and coloured boxes) were aligned to the genome assemblies and the contig positions of the top two hits of each gene were recorded. Contigs containing at least five shared genes were considered as potential haplotype pairs. Sequence collinearity between contigs was assessed by alignment, and homologous matching contigs were assigned to bins. Examples shown are for Bin04 and Bin12 from <italic>Pgt</italic>
21-0 and Ug99, respectively</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig1_HTML" id="d29e690"></graphic>
</fig>
</p>
<p id="Par8">The <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
genes encode dominant avirulence factors recognised by wheat resistance genes<sup><xref ref-type="bibr" rid="CR23">23</xref>
,<xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. These two genes are located in close proximity (<15 kbp) to each other in the genome assemblies and both haplotypes of this locus were assembled as alternate contigs in <italic>Pgt</italic>
21-0 and Ug99 (Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
). Both isolates were heterozygous for <italic>AvrSr50</italic>
with one allele containing a ~26-kbp-insertion. <italic>Pgt</italic>
21-0 was also heterozygous for <italic>AvrSr35</italic>
, with one allele containing a 400-bp-MITE insertion previously described<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. PCR amplification from the Ug99 strain had previously identified only a single <italic>AvrSr35</italic>
allele, suggesting that it was homozygous in Ug99<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. However, the Ug99 genome assembly contained a second allele of this gene, which contained a 57-kbp-insertion that would have prevented its PCR amplification in the Salcedo et al. study<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. The presence of this insertion in Ug99 was supported by DNA read (PacBio and Illumina) alignments across this genomic region and confirmed by DNA amplification and amplicon sequencing of flanking border regions (Supplementary Fig. <xref rid="MOESM1" ref-type="media">1</xref>
). Thus, Ug99 is also heterozygous for avirulence on <italic>Sr35</italic>
, and may therefore mutate to virulence on this wheat resistance gene more readily than if it were homozygous. This is an important finding as it will inform <italic>Sr35</italic>
deployment strategies against Ug99. Strikingly, the <italic>AvrSr35-virSr50</italic>
haplotype of this locus is very similar in structure in Ug99 and <italic>Pgt</italic>
21-0 (Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
) and shares >99% sequence identity, while the two alternate haplotypes are quite different. We therefore compared the larger genomic regions containing these loci in each isolate: namely bin 06 in <italic>Pgt</italic>
21-0, which is ~3.5 Mbp and includes telomeres at both ends, and bins 15 (1.8 Mbp) and 23 (1.2 Mbp) in Ug99 (Supplementary Fig. <xref rid="MOESM1" ref-type="media">2a</xref>
). One haplotype (designated A) was >99.7% identical in Ug99 and <italic>Pgt</italic>
21-0, while the other two haplotypes (B and C) were highly divergent from each other and from haplotype A (Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
, Supplementary Fig. <xref rid="MOESM1" ref-type="media">2b</xref>
and Supplementary Table <xref rid="MOESM1" ref-type="media">4</xref>
). Only 71–78% of the sequences from each haplotype aligned, with an average identity of ~95% across the aligned regions, yielding total identities of only 68–76% between these three distinct haplotypes. The high similarity between the A haplotypes of this chromosome suggested that Ug99 and <italic>Pgt</italic>
21-0 may share large portions of their genomes, potentially up to an entire haploid genome copy.<fig id="Fig2"><label>Fig. 2</label>
<caption><p>A common haplotype containing <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
is shared between <italic>Pgt</italic>
21-0 and Ug99. <bold>a</bold>
Diagram of genomic regions containing <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
alleles in <italic>Pgt</italic>
21-0 and Ug99. Numbers above tracks correspond to contig coordinates and the sense of the DNA strand is indicated as + or −. Predicted gene models (including introns) are depicted as dark grey boxes and intergenic spaces are shown in light grey. Coloured arrows indicate location and direction of <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
genes, with size and position of insertions shown in yellow. Intergenic distances between <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
are indicated by brackets. <bold>b</bold>
Total sequence identity between contigs representing homologous chromosomes of different haplotypes (coloured bars) containing the <italic>AvrSr50</italic>
/<italic>AvrSr35</italic>
locus (dotted white boxes). Telomere sequences are represented in grey. Chromosome size = ~3.5 Mbp</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig2_HTML" id="d29e831"></graphic>
</fig>
</p>
</sec>
<sec id="Sec4"><title>Whole-genome haplotype assignment and comparison</title>
<p id="Par9">We used a read subtraction and mapping approach (Fig. <xref rid="Fig3" ref-type="fig">3a, b</xref>
) to identify genome regions in the Ug99 and <italic>Pgt</italic>
21-0 assemblies that showed high similarity and may be derived from a shared haplotype. Illumina reads from each isolate were mapped to the genome reference of the other isolate. Reads that failed to map were retained, and this subtracted read set, in which sequences common to both isolates had been removed, was re-mapped to their original genome reference. The subtracted read coverage depth was compared to the coverage depth obtained when using all reads (Supplementary Fig. <xref rid="MOESM1" ref-type="media">3a, b</xref>
). Contigs representing sequences shared by both isolates (designated as haplotype A) displayed a very low subtracted read coverage depth (Supplementary Data <xref rid="MOESM6" ref-type="media">3</xref>
). In contrast, sequences unique to <italic>Pgt</italic>
21-0 or Ug99 (designated as haplotypes B or C, respectively) retained a relatively high subtracted read coverage depth. Some contigs in each assembly appeared to be chimeric with distinct regions assigned to opposite haplotypes, and these contigs were divided into separate fragments (Supplementary Data <xref rid="MOESM7" ref-type="media">4</xref>
) for subsequent haplotype comparisons. Approximately half of each genome assembly was assigned to either the A, B or C haplotypes (Fig. <xref rid="Fig3" ref-type="fig">3c</xref>
) and importantly one set of homologous sequences from each bin was assigned to each haplotype (Supplementary Data <xref rid="MOESM6" ref-type="media">3</xref>
). The A, B and C haplotype sets contained 95–96% of conserved fungal BUSCO genes (Fig. <xref rid="Fig3" ref-type="fig">3d</xref>
), indicating that each represents a full haploid genome equivalent. Consistent with this, the haplotypes were highly contiguous (Fig. <xref rid="Fig4" ref-type="fig">4</xref>
). Overall sequence identity between the A haplotypes of <italic>Pgt</italic>
21-0 and Ug99 was over 99.5%, with structural variation (large insertions/deletions) representing only 0.46% of the haplotype sizes, and with just 0.08% sequence divergence in aligned regions (Fig. <xref rid="Fig4" ref-type="fig">4a</xref>
, Table <xref rid="Tab1" ref-type="table">1</xref>
, and Supplementary Table <xref rid="MOESM1" ref-type="media">5</xref>
). In contrast, only 91–93% of the A, B or C haplotypes could be aligned to each other, with an average sequence identity of ~96% across the aligned regions, giving total identities of 87–91% (Fig. <xref rid="Fig4" ref-type="fig">4b–d</xref>
, Table <xref rid="Tab1" ref-type="table">1</xref>
and Supplementary Table <xref rid="MOESM1" ref-type="media">5</xref>
). Structural variation between these haplotypes accounted for 6.7–8.6% of the haploid genome sizes. There were only ~9000 SNPs (0.1/kbp) between the two A haplotypes, versus 876,000 to 1.4 million SNPs (11–18/kbp) between the A, B and C pairs, which is consistent with estimates of heterozygosity levels in <italic>Pgt</italic>
21-0 (haplotypes A and B) based on variant detection from Illumina read mapping<sup><xref ref-type="bibr" rid="CR21">21</xref>
</sup>
. The high similarity between the A haplotypes, and divergence between A, B and C haplotypes was also supported by Illumina read coverage and SNP calling analysis (Supplementary Fig. <xref rid="MOESM1" ref-type="media">3c–e</xref>
) showing that Ug99 and <italic>Pgt</italic>
21-0 share one nearly identical haploid genome copy.<fig id="Fig3"><label>Fig. 3</label>
<caption><p>Haplotype assignment by read subtraction and mapping process. <bold>a</bold>
Illumina reads from <italic>Pgt</italic>
21-0 were mapped to the Ug99 genome assembly at high stringency. Unmapped reads derived from divergent regions of the B haplotype were retained and then mapped to the <italic>Pgt</italic>
21-0 genome assembly. Read coverage of individual contigs with the original and subtracted reads were compared to designate haplotypes as either A or B. <bold>b</bold>
The same process was followed with reads from Ug99 subtracted against the <italic>Pgt</italic>
21-0 reference to designate the A and C haplotypes. <bold>c</bold>
Pie chart showing proportion and total sizes of contigs assigned to haplotypes A, B or C or unassigned in <italic>Pgt</italic>
21-0 and Ug99 assemblies. <bold>d</bold>
BUSCO analysis to assess completeness of haplotype genome assemblies. Bars represent the percentage of total BUSCOs as depicted by the colour key</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig3_HTML" id="d29e937"></graphic>
</fig>
<fig id="Fig4"><label>Fig. 4</label>
<caption><p><italic>Pgt</italic>
21-0 and Ug99 share one nearly identical haploid genome. <bold>a</bold>
–<bold>d</bold>
Dot plots illustrating sequence alignment of complete haplotypes. <italic>X</italic>
- and <italic>y</italic>
-axes show cumulative size of the haplotype assemblies depicted by coloured bars to the right and top of the graphs. Colour key indicates sequence identity ratios for all dot plots</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig4_HTML" id="d29e960"></graphic>
</fig>
<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Intra- and inter-isolate sequence comparison of entire haplotypes in Ug99 and <italic>Pgt</italic>
21-0</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th></th>
<th>Sequence similarity</th>
<th></th>
<th colspan="3">Structural variation</th>
</tr>
<tr><th></th>
<th></th>
<th></th>
<th></th>
<th colspan="2">Total variant size</th>
</tr>
<tr><th>Isolate comparison</th>
<th>Bases aligned (%)</th>
<th>Sequence divergence (%)</th>
<th>Number of variants</th>
<th>Mbp</th>
<th>% of genome</th>
</tr>
</thead>
<tbody><tr><td>21-0A vs Ug99 A</td>
<td>99.64</td>
<td>0.08</td>
<td>491</td>
<td>0.82</td>
<td>0.46</td>
</tr>
<tr><td>Ug99 A vs Ug99 C</td>
<td>91.52</td>
<td>4.08</td>
<td>2571</td>
<td>13.69</td>
<td>7.88</td>
</tr>
<tr><td>21-0A vs <italic>Pgt</italic>
21-0 B</td>
<td>91.38</td>
<td>4.19</td>
<td>2696</td>
<td>15.01</td>
<td>8.56</td>
</tr>
<tr><td>21-0 B vs Ug99 C</td>
<td>93.44</td>
<td>2.4</td>
<td>1910</td>
<td>11.50</td>
<td>6.69</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
</sec>
<sec id="Sec5"><title>Assessment of inter-nuclear recombination</title>
<p id="Par10">We tested two hypotheses that could explain the shared haplotype between Ug99 and <italic>Pgt</italic>
21-0: (1) Ug99 arose by a somatic hybridisation event in which an isolate of the race 21 lineage donated an intact nucleus of the A haplotype (Fig. <xref rid="Fig5" ref-type="fig">5a</xref>
); and (2) Ug99 arose by a sexual cross in which one haploid pycnial parent was derived from a race 21 lineage isolate after meiosis (Fig. <xref rid="Fig5" ref-type="fig">5b</xref>
). Under both scenarios, the A haplotype of Ug99 represents one entire haploid nucleus that was derived from the race 21 lineage isolate. In the nuclear exchange scenario, the <italic>Pgt</italic>
21-0 A haplotype represents a single nucleus donated intact to generate Ug99. However, under the sexual cross model, this <italic>Pgt</italic>
21-0 haplotype would include segments of both nuclear genomes that were combined by crossing over and chromosome reassortment after karyogamy and meiosis. Although recombination frequency has not been measured in <italic>Pgt</italic>
, an average of 115 recombination events per haploid genome was detected during meiosis in the related flax rust fungus (<italic>Melampsora lini</italic>
), which also has 18 chromosomes<sup><xref ref-type="bibr" rid="CR25">25</xref>
</sup>
. Because the <italic>Pgt</italic>
21-0 and Ug99 genome assemblies represent the phased dikaryotic state of each isolate, all correctly phased contigs in Ug99 should be either A or C haplotype, while those in <italic>Pgt</italic>
21-0 would include mixed haplotype contigs only if the sexual cross hypothesis is correct. In fact, just 19 contigs in the Ug99 assembly contained adjacent regions of either the A or C haplotype. These appeared to result from haplotype phase swap artefacts during genome assembly, since all of the junctions occurred between phase blocks (i.e. opposite gaps between the corresponding alternate contigs). Furthermore, examination of Illumina read mapping to these regions revealed that these sites contained either collapsed haplotype sequences, non-unique sequences or discontinuities in read coverage (Supplementary Fig. <xref rid="MOESM1" ref-type="media">4</xref>
), indicative of assembly errors disrupting phase information across the junction. Likewise, 31 contigs of mixed haplotype in the <italic>Pgt</italic>
21-0 assembly all contained likely phase swap artefacts (Supplementary Fig. <xref rid="MOESM1" ref-type="media">4</xref>
). To experimentally distinguish between phase-swap assembly artefacts and meiotic recombination events, we used Hi-C chromatin cross-linking proximity analysis<sup><xref ref-type="bibr" rid="CR26">26</xref>
</sup>
to assess physical linkage between contigs in the <italic>Pgt</italic>
21-0 assembly. About 90% of all read pair connections were between contigs of the same haplotype. For each of the chimeric contigs, the separated A and B fragments showed significantly more connections to contigs of the same haplotype than to contigs of the other haplotype, including other fragments of the original chimeric contig (Supplementary Data <xref rid="MOESM8" ref-type="media">5</xref>
). These observed physical linkages confirmed that all of the mixed haplotype contigs in <italic>Pgt</italic>
21-0 resulted from phase swap errors during the genome assembly process, and do not correspond to sites of genetic recombination that could have occurred during meiosis under the sexual origin hypothesis. Furthermore, only four structural variants larger than 10 kbp were detected between the A haplotypes of Ug99 and <italic>Pgt</italic>
21-0. Three of these were deletions in Ug99, while the fourth was a tandem repeat duplication, indicating that there was no novel genetic content in Ug99 haplotype A compared to <italic>Pgt</italic>
21-0.<fig id="Fig5"><label>Fig. 5</label>
<caption><p>Models for the emergence of the founder isolate of the <italic>Pgt</italic>
Ug99 lineage. <bold>a</bold>
A somatic hybridisation event and nuclear exchange occurred between an isolate of the <italic>Pgt</italic>
21 lineage and an unknown <italic>Pgt</italic>
isolate. The combination of nuclei A and C yielded the parental isolate of the Ug99 lineage in Africa. Under this scenario, nucleus A of Ug99 is entirely derived from nucleus A in <italic>Pgt</italic>
21-0. <bold>b</bold>
Alternatively, sexual reproduction and mating between these two parental isolates defined the origin of the Ug99 lineage. Under this scenario, meiotic recombination and chromosome reassortment would result in the <italic>Pgt</italic>
21-0-derived A nucleus of Ug99 being composed of a mosaic of the two haploid nuclear genomes of <italic>Pgt</italic>
21-0 (X and Y)</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig5_HTML" id="d29e1195"></graphic>
</fig>
</p>
</sec>
<sec id="Sec6"><title>Chromosome assembly and assessment of reassortment</title>
<p id="Par11">Combining Hi-C scaffolding data with the bin and haplotype assignment information for the <italic>Pgt</italic>
21-0 assembly allowed us to construct 18 chromosome pseudomolecules for each of the A and B haplotypes (Fig. <xref rid="Fig6" ref-type="fig">6a</xref>
, Supplementary Table <xref rid="MOESM1" ref-type="media">6</xref>
and Supplementary Data <xref rid="MOESM9" ref-type="media">6</xref>
). These covered a total of 170 Mbp and ranged from 2.8 to 7.3 Mbp in size, consistent with relative chromosome sizes from karyotype analysis<sup><xref ref-type="bibr" rid="CR27">27</xref>
</sup>
. The A and B chromosomes were collinear except for two translocation events (Fig. <xref rid="Fig6" ref-type="fig">6b</xref>
). In each case these were supported by contigs that spanned the translocation breakpoints. Re-scaffolding the separated fragments of these contigs using Hi-C data supported the original contig assembly, indicating that these are true translocation events within the A or B genomes. The haplotype A chromosomes showed high collinearity with the Ug99 A haplotype contigs (Fig. <xref rid="Fig6" ref-type="fig">6c</xref>
).<fig id="Fig6"><label>Fig. 6</label>
<caption><p>Chromosome sets of haplotype A and B in <italic>Pgt</italic>
21-0. <bold>a</bold>
Schematic representation of assembled chromosomes for <italic>Pgt</italic>
21-0 of each haplotype (scale bar = 1 Mbp). Horizontal bars indicate telomeric repeat sequences. <bold>b</bold>
Dot plot of sequence alignment of <italic>Pgt</italic>
21-0 chromosome pseudomolecules of haplotypes A and B. Two translocation events, one between chromosomes 3 and 5 and one between chromosomes 8 and 16, are evident. <bold>c</bold>
Dot plot of sequence alignment between chromosomes from haplotype A in <italic>Pgt</italic>
21-0 and contigs from haplotype A in Ug99. <bold>d</bold>
Percentage of Hi-C read pairs linking each A haplotype chromosome to other A chromosomes (blue) or to B haplotype chromosomes (orange). <bold>e</bold>
Percentage of Hi-C read pairs linking each B haplotype chromosome to either A (blue) or B (orange) chromosomes</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig6_HTML" id="d29e1260"></graphic>
</fig>
</p>
<p id="Par12">Approximately 65% of the total Hi-C read pairs represented links between physically contiguous sequences on the same chromosome, while the remaining pairs connected sites distributed across the genome. Because Hi-C DNA crosslinking is performed in intact cells, these non-scaffolding linkages should preferentially form between chromosomes that are located in the same nucleus. Indeed, all chromosomes of the A haplotype showed a much higher proportion of Hi-C read pair links to other chromosomes of the A haplotype (~85%) than to chromosomes of the B haplotype (~15%) (Fig. <xref rid="Fig6" ref-type="fig">6d</xref>
), indicating that they are all located in the same nucleus. Thus, the A haplotype of Ug99 derives from a single nucleus of <italic>Pgt</italic>
21-0. Similarly, 17 of the B haplotype chromosomes in <italic>Pgt</italic>
21-0 showed stronger linkage to other B chromosomes (~90%) than to A chromosomes (~10%) (Fig. <xref rid="Fig6" ref-type="fig">6e</xref>
). However, chromosome 11B showed the inverse, suggesting that both homologues of this chromosome are in the same nucleus. This implies that a single chromosome exchange event occurred during asexual propagation of the <italic>Pgt</italic>
21-0 isolate, after its divergence from the race 21 lineage branch giving rise to Ug99.</p>
<p id="Par13">Overall the whole-genome comparison data demonstrate that Ug99 shares one full haploid nuclear genome with the <italic>Pgt</italic>
21-0 isolate with no recombination events within chromosomes and no reassortment of chromosomes from different nuclei. These facts are inconsistent with a sexual origin, and strongly support that the Ug99 lineage arose by a somatic hybridisation event involving one parent derived from the African race 21 lineage and another parent of unknown origin exchanging whole nuclei (Fig. <xref rid="Fig7" ref-type="fig">7</xref>
).<fig id="Fig7"><label>Fig. 7</label>
<caption><p>Model for Ug99 origin by somatic hybridisation and nuclear exchange. The ancestral isolate of the Ug99 lineage acquired the A and C genomes from an isolate of the 21 lineage and an unknown isolate and later gained virulence to wheat cultivars carrying the <italic>Sr31</italic>
resistance gene</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig7_HTML" id="d29e1297"></graphic>
</fig>
</p>
</sec>
<sec id="Sec7"><title>Comparison of gene content between haplotypes</title>
<p id="Par14">Annotation of the <italic>Pgt</italic>
21-0 and Ug99 genome assemblies predicted similar gene numbers (~18,000) in each haplotype (Supplementary Table <xref rid="MOESM1" ref-type="media">7</xref>
). Gene orthology analysis indicated that 65–70% of genes in each of the A, B and C haplotypes were shared and represent a core <italic>Pgt</italic>
gene set, while the remainder were present in only one or two haplotypes (Supplementary Table <xref rid="MOESM1" ref-type="media">8</xref>
). Mapping of orthologous gene pairs supported the overall synteny of the <italic>Pgt</italic>
21-0 A and B chromosome assemblies and confirmed the translocations observed between chromosomes 3 and 5, and between 8 and 16 (Fig. <xref rid="Fig8" ref-type="fig">8a</xref>
). Genes encoding secreted and non-secreted proteins showed a similar distribution across the chromosomes, while repeat sequences displayed an inverse distribution to genes (Fig. <xref rid="Fig8" ref-type="fig">8b</xref>
, Supplementary Fig. <xref rid="MOESM1" ref-type="media">5</xref>
). The location of secreted protein genes in gene-rich rather than repeat-rich regions is consistent with the absence of two-speed genome architecture in the related rust fungal species <italic>P. coronata</italic>
and <italic>P. striiformis</italic>
<sup><xref ref-type="bibr" rid="CR12">12</xref>
,<xref ref-type="bibr" rid="CR13">13</xref>
</sup>
.<fig id="Fig8"><label>Fig. 8</label>
<caption><p>Gene content of <italic>Pgt</italic>
21-0 chromosome pseudomolecules. <bold>a</bold>
Circos plot showing location of orthologous gene pairs in the A and B chromosomes of <italic>Pgt</italic>
21-0. <bold>b</bold>
Gene and repeat density plots for homologous chromosomes 14 A and 14B. Density of genes encoding non-secreted (black) or secreted proteins (red) along the chromosomes are shown, with individual genes indicated by black or red dots. Bottom graph shows density of repeat elements (blue). Positions of <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
genes are indicated</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig8_HTML" id="d29e1368"></graphic>
</fig>
</p>
<p id="Par15">Both Ug99 and <italic>Pgt</italic>
21-0 are heterozygous at the predicted <italic>a</italic>
and <italic>b</italic>
mating type loci (Supplementary Fig. <xref rid="MOESM1" ref-type="media">6</xref>
), despite Ug99 being derived by a non-sexual mechanism. This is consistent with an expectation that formation and maintenance of a stable dikaryon requires two distinct compatible mating types<sup><xref ref-type="bibr" rid="CR11">11</xref>
</sup>
. We observed multiple alleles at the <italic>b</italic>
locus on chromosome 9, which encodes the divergently transcribed transcription factors bE and bW, with variants sharing 70–80% amino acid identity. Cuomo and colleagues<sup><xref ref-type="bibr" rid="CR28">28</xref>
</sup>
previously described two alleles, <italic>b1</italic>
and <italic>b2</italic>
, for this locus in <italic>Pgt</italic>
. <italic>Pgt</italic>
21-0 contained the <italic>b1</italic>
allele and a novel <italic>b3</italic>
allele, while Ug99 contained <italic>b3</italic>
and another novel <italic>b4</italic>
allele. Both isolates were heterozygous for the + and – alleles of the <italic>a</italic>
locus on chromosome 4, which encodes a pheromone (mfa) and pheromone receptor (STE3) pair. However, one of the receptor alleles in Ug99 contained a single-nucleotide deletion that resulted in truncation of the last 48 amino acids of the protein. Thus, the mating type system for <italic>Pgt</italic>
appears to consist of two independent loci, one di-allelic and one multi-allelic.</p>
</sec>
<sec id="Sec8"><title>Phylogenetic analysis of global <italic>Pgt</italic>
isolates</title>
<p id="Par16">We used the haplotype-phased genome references for <italic>Pgt</italic>
21-0 and Ug99 to determine genetic relationships within a set of global <italic>Pgt</italic>
isolates using publicly available sequence data<sup><xref ref-type="bibr" rid="CR21">21</xref>
,<xref ref-type="bibr" rid="CR23">23</xref>
,<xref ref-type="bibr" rid="CR29">29</xref>
</sup>
. Maximum likelihood trees based on whole-genome SNPs (Fig. <xref rid="Fig9" ref-type="fig">9a</xref>
and Supplementary Fig. <xref rid="MOESM1" ref-type="media">7a</xref>
) showed a very similar overall topology to that reported previously for most of these isolates<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
. The 5 isolates of the Ug99 lineage, and the 13 South African and Australian isolates each formed a separate tight clade, consistent with their proposed clonal nature<sup><xref ref-type="bibr" rid="CR18">18</xref>
–<xref ref-type="bibr" rid="CR20">20</xref>
</sup>
. However, tree building using filtered SNPs from just the A haplotype resulted in the formation of a single clade containing the Ug99, South African and Australian isolates, which indicates the clonal derivation of this nucleus among these isolates (Fig. <xref rid="Fig9" ref-type="fig">9b</xref>
, Supplementary Fig. <xref rid="MOESM1" ref-type="media">7b, c</xref>
). The Ug99 group forms a subclade within the race 21 group consistent with a derived origin. In contrast, these groups remained in two distant clades in phylogenies inferred using filtered SNPs from the B genome. Surprisingly, in this case two isolates from the Czech Republic and three isolates from Pakistan were now located in a single clade with the South African and Australian isolates (Fig. <xref rid="Fig9" ref-type="fig">9c</xref>
). This suggests that these isolates contain a haplotype closely related to the B genome of the race 21 lineage and may also have arisen by somatic hybridisation and nuclear exchange. A phylogeny based on the C genome SNPs grouped isolate IR-01 from Iran with the Ug99 lineage (Fig. <xref rid="Fig9" ref-type="fig">9d</xref>
), suggesting that these isolates share the C haplotype. IR-01 could represent a member of the parental lineage that donated the C nucleus to Ug99, or alternatively may have acquired the C nucleus from Ug99. Notably, this was the only isolate that shared the <italic>AvrSr35</italic>
57kbp insertion allele identified in Ug99 (Supplementary Fig. <xref rid="MOESM1" ref-type="media">8a</xref>
). The relationships between these putative hybrid isolates were also supported by the patterns of homozygous and heterozygous SNPs detected in each haplotype (Supplementary Fig. <xref rid="MOESM1" ref-type="media">8b–e</xref>
). The incongruities between phylogenies generated based on different haplotypes highlight the difficulty of inferring relationships between isolates based on whole-genome SNP data without haplotype resolution. Overall, these observations suggest that somatic hybridisation and nuclear exchange may be a common mechanism generating genetic diversity in global populations of <italic>Pgt</italic>
.<fig id="Fig9"><label>Fig. 9</label>
<caption><p>Somatic hybridisation in <italic>Pgt</italic>
evolution. <bold>a</bold>
Phylogenetic analysis of <italic>Pgt</italic>
isolates from diverse countries of origin (colour key) using a RAxML model and biallelic SNPs called against the full dikaryotic genome of <italic>Pgt</italic>
21-0. Scale bar indicates number of nucleotide substitutions per site. Red asterisks indicate <italic>P. graminis</italic>
f. sp. <italic>avenae</italic>
isolates used as an outgroup. <bold>b</bold>
Dendrogram inferred using biallelic SNPs detected against haplotype A of <italic>Pgt</italic>
21-0. <bold>c</bold>
Dendrogram inferred using SNPs detected against haplotype B of <italic>Pgt</italic>
21-0. <bold>d</bold>
Dendrogram inferred from SNPs detected in haplotype C of Ug99</p>
</caption>
<graphic xlink:href="41467_2019_12927_Fig9_HTML" id="d29e1535"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec9" sec-type="discussion"><title>Discussion</title>
<p id="Par17">Although sexual reproduction of <italic>Pgt</italic>
can generate individuals with novel genetic combinations, the completion of the sexual cycle requires infection of an alternate host, <italic>Berberis</italic>
spp. (barberry)<sup><xref ref-type="bibr" rid="CR30">30</xref>
</sup>
. In parts of the world where barberry is scarce or absent, either due to eradication programmes or its natural distribution, <italic>Pgt</italic>
is restricted to asexual propagation with new diversity often arising by mutation or migration<sup><xref ref-type="bibr" rid="CR19">19</xref>
,<xref ref-type="bibr" rid="CR20">20</xref>
</sup>
. Somatic hybridisation provides an alternative explanation for the appearance of new races not derived by stepwise mutation. Hybrids with high adaptive value in agroecosystems may establish new lineages of epidemiological significance, as shown by the emergence of the Ug99 lineage with its substantial impact on East African wheat production and threat to global food security<sup><xref ref-type="bibr" rid="CR14">14</xref>
</sup>
. Results from experiments mostly conducted prior to the molecular biology era suggested the possibility of somatic genetic exchange between rust isolates co-infecting the same host, based on the detection of novel virulence phenotypes<sup><xref ref-type="bibr" rid="CR5">5</xref>
–<xref ref-type="bibr" rid="CR8">8</xref>
,<xref ref-type="bibr" rid="CR31">31</xref>
</sup>
. In some cases only two non-parental phenotypic classes were observed, consistent with simple nuclear exchange<sup><xref ref-type="bibr" rid="CR7">7</xref>
,<xref ref-type="bibr" rid="CR8">8</xref>
</sup>
. The isolation of additional recombinant classes in other experiments was interpreted as recombination between nuclei, although this was also ascribed to isolate contamination<sup><xref ref-type="bibr" rid="CR7">7</xref>
,<xref ref-type="bibr" rid="CR32">32</xref>
</sup>
, which could not be ruled out without molecular data. Our data show that Ug99 arose by a somatic nuclear exchange event with no recombination. Phylogenetic analyses show that at least two different global lineages (from Pakistan and the Czech Republic) share a haplotype similar to the <italic>Pgt</italic>
21-0 B genome, while another (from Iran) shares a haplotype closely related to the Ug99 C genome. This suggests that multiple nuclear exchange events between strains have occurred in the global <italic>Pgt</italic>
population and had significant impacts on genetic diversity. Although we did not observe any recombination between nuclei associated with the Ug99 hybridisation event, we did see evidence for translocation of one complete chromosome between nuclei in <italic>Pgt</italic>
21-0. We also previously found that a <italic>Pgt</italic>
mutant virulent on <italic>Sr50</italic>
arose by exchange of an ~2.5 Mbp region between two haplotypes<sup><xref ref-type="bibr" rid="CR23">23</xref>
</sup>
. Thus, genetic exchange between haploid nuclei may occur as rare events during asexual propagation of a single lineage in <italic>Pgt</italic>
. Whether extensive genetic exchange similar to ascomycete parasexuality<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
can also occur between rust nuclei during hybridisation remains to be determined. This may require controlled infection experiments, as such recombinant hybrids would be difficult to distinguish from the products of sexual recombination in field derived-strains, especially given the potential for long-range spore dispersal.</p>
<p id="Par18">Although there is now clear evidence of nuclear exchange between dikaryons in <italic>Pgt</italic>
, nothing is known of how this process occurs or is regulated. It differs from parasexuality in ascomycetes<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
, as the dikaryotic state is maintained with no nuclear fusion or haploidisation resulting in chromosome reassortment. Wang and McCallum<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
observed the formation of fusion bodies where germ tubes of different <italic>P. triticina</italic>
isolates came into contact, with the potential for nuclear exchange at these junctions. It has been proposed that mating type loci contribute to determining the compatibility between isolates for the formation of hybrids<sup><xref ref-type="bibr" rid="CR6">6</xref>
,<xref ref-type="bibr" rid="CR7">7</xref>
</sup>
, but this also remains to be confirmed experimentally. Our findings provide a new framework to take advantage of haplotype genome resolution to understand the role of somatic exchange in population diversity of rust fungi.</p>
<p id="Par19">Extended dikaryotic developmental stages are common in many other fungi, especially basidiomycetes. Indeed, separation of karyogamy (fusion of haploid nuclei to form a diploid nucleus) from gamete fusion is a feature unique to the fungal kingdom<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. However, it is unclear why fungi maintain an extended dikaryotic stage prior to formation of a diploid nucleus as a precursor to sexual reproduction<sup><xref ref-type="bibr" rid="CR34">34</xref>
</sup>
. One possibility is that the ability to exchange haploid nuclei offers an advantage over the diploid state due to the enhanced genetic variation in long-lived asexual dikaryotes. There is also evidence for somatic exchange of genetic markers in dikaryotes of the mushroom <italic>Schizophyllum commune</italic>
, which belongs to another Basidiomycete subphylum, Agaricomycotina<sup><xref ref-type="bibr" rid="CR35">35</xref>
</sup>
. Arbuscular mycorrhizae (AM) fungi are another ancient fungal lineage whose spores contain hundreds of nuclei and for which no sexual stages have been described, raising questions of how these lineages have survived<sup><xref ref-type="bibr" rid="CR36">36</xref>
</sup>
. Recently some dikaryotic-like AM isolates possessing two divergent classes of nuclei have been observed. Nuclear exchange between dikaryotes could be another driver of genetic variation in these fungi. Evidently, the members of the fungal kingdom display remarkable genetic plasticity and further investigation is required to reveal the mechanism, prevalence and evolutionary importance of nuclear exchange in dikaryotic and multinucleate fungi.</p>
</sec>
<sec id="Sec10"><title>Methods</title>
<sec id="Sec11"><title>Fungal stocks and plant inoculation procedures</title>
<p id="Par20"><italic>Pgt</italic>
isolates Ug99<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
, UVPgt55, UVPgt59, UVPgt60 and UVPgt61 collected in South Africa<sup><xref ref-type="bibr" rid="CR16">16</xref>
,<xref ref-type="bibr" rid="CR37">37</xref>
</sup>
were transferred to the Biosafety Level 3 (BSL-3) containment facility at the University of Minnesota for growth and manipulation. Samples were purified by single pustule isolation and then amplified by 2–3 rounds of inoculation on the susceptible wheat cultivar McNair. Virulence pathotypes and purity of each isolate were confirmed by inoculation onto the standard wheat differential set (Supplementary Data <xref rid="MOESM4" ref-type="media">1</xref>
). Other isolates used in this study were <italic>Pgt</italic>
21-0, which was first isolated in Australia in 1954<sup><xref ref-type="bibr" rid="CR20">20</xref>
,<xref ref-type="bibr" rid="CR21">21</xref>
</sup>
, the North American isolate CRL 75-36-700-3 (pathotype SCCL)<sup><xref ref-type="bibr" rid="CR22">22</xref>
</sup>
and Kenyan isolate 04KEN156/04 (pathotype TTKSK)<sup><xref ref-type="bibr" rid="CR17">17</xref>
</sup>
. For rust inoculations, urediniospores retrieved from −80 °C were heat treated (45 °C for 15 min) and suspended in mineral oil (Soltrol 170, Philips Petroleum, Borger, TX, USA) at 14 mg/ml. Seven day-old seedlings were spray-inoculated at 50 μl/plant and kept in a dark mist chamber at 22–25 °C, 100% humidity for 16 h. Subsequently, plants were exposed to light (150–250 μmol photons s<sup>−1</sup>
m<sup>−2</sup>
) for 5 h in the mist chamber and then transferred to a controlled growth chamber (18 h/6 h of light/dark, 24 °C/18 °C for day/night, 50% relative humidity). Spores were collected 9 and 14 days post inoculation (dpi) and maintained at −80 °C.</p>
</sec>
<sec id="Sec12"><title>DNA extraction and sequencing of rust isolates</title>
<p id="Par21">High molecular weight DNA of Ug99 and <italic>Pgt</italic>
21-0 was extracted from 300 to 350 mg urediniospores as described<sup><xref ref-type="bibr" rid="CR38">38</xref>
</sup>
, with the following modifications: (1) for Phenol:Chloroform:Isoamyl alcohol extractions, samples were centrifuged at 4 °C and 5000 × <italic>g</italic>
for 20 mins; (2) a wide-bore 1-mL pipette tip was used to transfer the DNA pellet; (3) samples were incubated for 1 h at 28 °C with 200–250 rpm shaking to dissolve the final DNA pellet. Double-stranded DNA concentration was quantified using a broad-range assay in a Qubit Fluorometer (Invitrogen, Carlsbad, CA, USA) and a NanoDrop (Thermo Fisher Scientific, Waltham, MA, USA). Approximately 10 µg DNA from Ug99 and <italic>Pgt</italic>
21-0 was sequenced using PacBio SMRT sequencing (Pacific Bioscience, Menlo Park, CA, USA) at either the Frederick National Laboratory for Cancer Research, Leidos Biomedical Research, Inc. (Frederick, MD, USA) or the Ramaciotti Centre (Sydney, Australia), respectively. DNA was concentrated and cleaned using AMPure PB beads for Ug99 or AMPure XP beads for <italic>Pgt</italic>
21-0 (Pacific Biosciences). DNA quantification and size assessment were conducted using a NanoDrop (Thermo Fisher Scientific) and 2200 TapeStation instruments (Agilent Technologies, Santa Clara, CA, USA). DNA was sheared to a targeted average size of 20 kb using G-tubes (Covaris, Woburn, MA, USA). Libraries were constructed following the 20 kb Template Preparation BluePippin Size-Selection System protocol (Pacific Biosciences) using a BluePippin instrument (Sage Science, Beverly, MA, USA) with a 0.75% agarose cassette and a lower cutoff of 15 kbp. For Ug99, 5 SMRT cells were sequenced on a PacBio Sequel platform using P6-C4 chemistry, the Sequel Binding Kit 2.0 (Pacific Biosciences), diffusion loading, 10-h movie lengths and Magbead loading at 2 pM (3 cells) or 4 pM (2 cells). In addition, 4 SMRT cells were run on PacBio RSII sequencer using P6-C4 chemistry, with 0.15 nM MagBead loading and 360-min movie lengths. For <italic>Pgt</italic>
21-0, 17 SMRT cells were run on the RSII platform using P6-C4 chemistry, Magbead loading (0.12–0.18 nM) and 240-min movie lengths.</p>
<p id="Par22">Genomic DNA for Illumina sequencing was extracted from 10 to 20 mg urediniospores of Ug99, UVPgt55, 59, 60 and 61 using the OmniPrep™ kit (G-Biosciences, St. Louis, MO, USA) following the manufacturer’s instructions. TruSeq Nano DNA libraries were prepared from 300 ng of DNA and 150 bp paired-end sequence reads were generated at the University of Minnesota Genomics Center on the Illumina NextSeq 550 platform using Illumina Real-Time Analysis software version 1.18.64 for quality-scored base calling.</p>
</sec>
<sec id="Sec13"><title>De novo long read assembly</title>
<p id="Par23">Genome assemblies of Ug99 and <italic>Pgt</italic>
21-0 were built from PacBio reads using Canu version 1.6<sup><xref ref-type="bibr" rid="CR39">39</xref>
</sup>
with default parameters and an estimated genome size of 170 Mbp. Assemblies were polished with the Arrow algorithm using the raw PacBio reads in the sa3_ds_resequencing pipeline in pbsmrtpipe workflow within SMRTLINK/5.1.0 (Pacific BioSciences). Assemblies were further polished by two rounds of Pilon<sup><xref ref-type="bibr" rid="CR40">40</xref>
</sup>
with the option fix --all using Illumina reads from Ug99 (this work) or <italic>Pgt</italic>
21-0 (NCBI SRA run Accession# SRR6242031). A BLASTN search (version 2.7.1) against the NCBI nr/nt database (downloaded on 4/11/2018) with E-value set as 1e<sup>−10</sup>
identified two contigs in the Ug99 assembly with significant hits to plant rRNA and chloroplast sequences and these were removed.</p>
<p id="Par24">PacBio and Illumina reads were mapped to the assembly using BWA-MEM (version 0.7.17)<sup><xref ref-type="bibr" rid="CR41">41</xref>
</sup>
and BAM files were indexed and sorted using SAMtools (version 1.9)<sup><xref ref-type="bibr" rid="CR42">42</xref>
</sup>
. Read coverage analysis using genomeCoverageBed in BEDtools (version 2.27.1)<sup><xref ref-type="bibr" rid="CR43">43</xref>
</sup>
identified 144 small contigs (<50 kbp) in the Ug99 assembly with low coverage (<2×) for both short and long reads and these contigs were also excluded from the final assembly. Genome assembly metrics were assayed using QUAST (version 4.3)<sup><xref ref-type="bibr" rid="CR44">44</xref>
</sup>
. Genome completeness was assessed via benchmarking universal single-copy orthologs (BUSCOs) of the basidiomycota as fungal lineage and <italic>Ustilago maydis</italic>
as the reference species for AUGUSTUS gene prediction<sup><xref ref-type="bibr" rid="CR45">45</xref>
</sup>
in the software BUSCO v2.0 (genome mode)<sup><xref ref-type="bibr" rid="CR46">46</xref>
</sup>
. Telomeric sequences were identified using either a high stringency BLAST with 32 repeats of TTAGGG as query or a custom python script to detect at least five CCCTAA or TTAGGG repeats in the assemblies. Repeats of at least 60 bp length and occurring within 100 bp of the contig end were defined as telomeric sequences.</p>
</sec>
<sec id="Sec14"><title>Detection of alternate contigs and bin assignment</title>
<p id="Par25">To identify contigs representing corresponding haplotypes (Fig. <xref rid="Fig1" ref-type="fig">1</xref>
) 22,484 predicted <italic>Pgt</italic>
gene coding sequences<sup><xref ref-type="bibr" rid="CR21">21</xref>
</sup>
were screened against the genome assemblies using BLITZ (Blat-like local alignment) in the Biokanga Tool set, (<ext-link ext-link-type="uri" xlink:href="https://github.com/csiro-crop-informatics/biokanga/releases/tag/v4.3.9">https://github.com/csiro-crop-informatics/biokanga/releases/tag/v4.3.9</ext-link>
). For each gene the two best hits (likely alleles) in the assembly were recorded. Contigs sharing best hits for at least five genes were selected as potential haplotype pairs and their sequence collinearity was examined by alignment and similarity plotting using D-genies<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
. Contigs representing contiguous or syntenous haplotypes were grouped together as bins.</p>
</sec>
<sec id="Sec15"><title>Validation of a 57-kbp-insert in <italic>AvrSr35</italic>
</title>
<p id="Par26">Contigs containing the <italic>AvrSr50</italic>
and <italic>AvrSr35</italic>
gene sequences were identified by BLASTN search against customised databases for the Ug99 and <italic>Pgt</italic>
21-0 genome assemblies. Illumina and PacBio reads of Ug99 mapped to the genome assembly were visualised in the Integrative Genomics Viewer (IGV). To validate the presence of the 57-kbp-insert in <italic>AvrSr35</italic>
, flanking and internal sequences were amplified from genomic DNA extracted using the OmniPrep™ kit (G-Biosciences) from urediniospores of Ug99, 04KEN156/04, and CRL 75-36-700-3. PCR was performed using Phusion high-fidelity DNA polymerase according to the manufacturer’s recommendations (New England BioLabs Inc., Ipswich, MA, USA) and primers listed in Supplementary Table <xref rid="MOESM1" ref-type="media">9</xref>
. The amplified PCR products were separated by electrophoresis on a 1% agarose gel and stained using SYBR Safe DNA gel stain (Invitrogen). Specific bands were cleaned using NucleoSpin gel clean-up kit (Takara Bio, Mountain View, CA, USA) for subsequent Sanger sequencing and alignment to <italic>AvrSr35</italic>
alleles. Gene models in the <italic>AvrSr35</italic>
and <italic>AvrSr50</italic>
locus were depicted using GenomicFeatures<sup><xref ref-type="bibr" rid="CR48">48</xref>
</sup>
and ggbio<sup><xref ref-type="bibr" rid="CR49">49</xref>
</sup>
in a custom R script.</p>
</sec>
<sec id="Sec16"><title>Haplotype assignment by read cross-mapping and subtraction</title>
<p id="Par27">Illumina reads from <italic>Pgt</italic>
21-0 (NCBI SRR6242031) were trimmed (“Trim sequences” quality limit = 0.01) and mapped to the Ug99 reference assembly using the “map reads to reference” tool in CLC Genomics Workbench version 10.0.1 or later with high stringency parameters (similarity fraction 0.99, length fraction 0.98, global alignment). Unmapped reads (Ug99-subtracted reads) were retained and then mapped back to the <italic>Pgt</italic>
21-0 assembly contigs using the same parameters. The original <italic>Pgt</italic>
21-0 reads were also mapped to the <italic>Pgt</italic>
21-0 assembly and the read coverage for each contig compared to the Ug99-subtracted reads. Contigs with very low coverage (<2X total and <10% of the original read coverage) with the Ug99-subtracted reads were designated as karyon A (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
, Supplementary Data <xref rid="MOESM6" ref-type="media">3</xref>
). Contigs with substantial coverage of Ug99-subtracted reads (>10% of the original read coverage) were designated as karyon B. Contigs with ambiguous read mapping data, including those with low coverage in the original unsubtracted reads or covered by largely non-uniquely mapping reads were left as unassigned. Read mapping to all contigs was confirmed by visual inspection of coverage graphs and read alignments in the CLC Genomics Workbench browser. Potential chimeric contigs were identified as containing distinct regions with either high or no coverage with the Ug99-subtracted reads (Supplementary Fig. <xref rid="MOESM1" ref-type="media">4</xref>
). For subsequent comparison and analyses, these contigs were manually split into their component fragments which were designated as haplotype A or B accordingly (Supplementary Data <xref rid="MOESM7" ref-type="media">4</xref>
). The same process was followed in reverse for the assignment of the A and C haplotype contigs in Ug99.</p>
</sec>
<sec id="Sec17"><title>Sequence comparisons of genome assemblies</title>
<p id="Par28">Haplotype sequences of the <italic>AvrSr50</italic>
/<italic>AvrSr35</italic>
chromosome as well as the full haploid genomes were aligned using MUMmer4.x<sup>50</sup>
with nucmer -maxmatch and other parameters set as default and the alignment metrics summarised in MUMmer dnadiff. Structural variation between haplotypes was determined using Assemblytics<sup><xref ref-type="bibr" rid="CR50">50</xref>
</sup>
from the MUMmer delta file with a minimum variant size of 50 bp, a maximum variant size of 100 kbp, and a unique sequence length for anchor filtering of 10 kbp. Haplotype dot plot alignments were generated using D-genies<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
.</p>
</sec>
<sec id="Sec18"><title>Read coverage analysis and SNP calling on haplotypes</title>
<p id="Par29">Illumina reads from Ug99 and <italic>Pgt</italic>
21-0 were each mapped against the Ug99 and <italic>Pgt</italic>
21-0 assemblies in CLC Genomics Workbench (similarity fraction 0.98, length fraction 0.95). For each assembly the mean coverage per base was calculated per 1000 bp interval (“window”) using samtools bedcov. Read coverage frequency normalised to the mean coverage of each haplotype was graphed as a violin plot using seaborn 0.9.0 package (<ext-link ext-link-type="uri" xlink:href="https://seaborn.pydata.org/">https://seaborn.pydata.org/</ext-link>
) using a custom python script. To detect SNPs between two haplotypes, Illumina read pairs of <italic>Pgt</italic>
21-0 that mapped uniquely to either the <italic>Pgt</italic>
21-0 A or B haplotype contigs were extracted. Similarly, Ug99-derived read pairs that uniquely mapped to either the A or C haplotype contigs of Ug99 were extracted. These read sets were then separately mapped to the two assemblies in CLC Genomics Workbench (similarity fraction 0.99, length fraction 0.98). Variant calling was performed using FreeBayes v.1.1.0<sup><xref ref-type="bibr" rid="CR51">51</xref>
</sup>
with default parameters in parallel operation and SNPs were filtered using vcffilter of VCFlib (v1.0.0-rc1, <ext-link ext-link-type="uri" xlink:href="https://github.com/vcflib/vcflib">https://github.com/vcflib/vcflib</ext-link>
) with the parameter -f “QUAL > 20 & QUAL / AO > 10 & SAF > 0 & SAR > 0 & RPR > 1 & RPL > 1”. Homozygous and heterozygous SNPs were extracted by vcffilter -f “AC > 0 & AC = 2” and -f “AC > 0 & AC = 1”, respectively. SNP statistics were calculated using vcfstats of VCFlib.</p>
</sec>
<sec id="Sec19"><title>Hi-C data analysis and scaffolding</title>
<p id="Par30">A Hi-C library was constructed with the ProxiMeta Hi-C kit from Phase Genomics v1.0 containing the enzyme Sau3A from ~150 mg of dried urediniospores of <italic>Pgt</italic>
21-0 following the standard protocol with minor modifications. Spores were washed in 1 mL 1X TBS buffer twice before cross-linking. After quenching of the crosslinking reaction, all liquid was removed and spores were frozen in liquid nitrogen. Spores were then lysed using cryogenic bead beating with two 5 mm steel beads shaking twice for 45 sec at 25 Hz using TissueLyser II (Qiagen). Lysis buffer was added and spores vortexed until full suspension. Reverse cross-linking was performed at 65 °C with 700 rpm horizontal shaking for 18 h. The Hi-C library was sequenced (150 bp paired-end reads) on the NextSeq 550 System using the Mid-Output Kit at the Ramaciotti Centre. The raw Hi-C reads were processed with the HiCUP pipeline version 0.7.1<sup><xref ref-type="bibr" rid="CR52">52</xref>
</sup>
(maximum di-tag length 700, minimum di-tag length 100, --re1 ^GATC,Sau3A), using bowtie2 as the aligner<sup><xref ref-type="bibr" rid="CR53">53</xref>
</sup>
and the <italic>Pgt</italic>
21-0 genome assembly as the reference. SAM files of the filtered di-tags were parsed to extract cis-far pairs (pairs located on the same contig and > 10 kbp apart) and trans pairs (located on different contigs). The numbers of trans pairs connecting each pair of contigs was extracted from this data. We compared the number of Hi-C read pair connections between contigs of the same or different haplotypes, either within bins, within chromosomes or at the whole-genome level (all haplotype-assigned contigs). We used a <italic>Χ</italic>
<sup>2</sup>
test to assess the deviation of each read pair distribution from a 1:1 ratio.</p>
<p id="Par31">For scaffolding, the raw Hi-C reads were first mapped to the <italic>Pgt</italic>
21-0 assembly using BWA-MEM<sup><xref ref-type="bibr" rid="CR41">41</xref>
</sup>
and processed using the Arima Genomics pipeline (<ext-link ext-link-type="uri" xlink:href="https://github.com/ArimaGenomics/mapping_pipeline/blob/master/01_mapping_arima.sh">https://github.com/ArimaGenomics/mapping_pipeline/blob/master/01_mapping_arima.sh</ext-link>
). Scaffolding was performed using SALSA 2.2<sup><xref ref-type="bibr" rid="CR54">54</xref>
</sup>
on the full set of contigs, as well as independently on the haplotype A or B sets of contigs (including unassigned contigs). Invalid scaffold linkages between adjacent telomeres, which occur as an artefact of telomere co-location within the nucleus, were discarded. The three sets of scaffolds were compared with the bin and haplotype assignment information to find overlaps and the resulting chromosome pseudomolecules were constructed by concatenating ordered contigs with 100 Ns inserted between contigs. To confirm translocations detected in the A and B chromosome sets, contigs that spanned the translocation site were separated into two fragments at the junction point and the SALSA scaffolding process was repeated on the full genome contig assembly. To detect nucleus-specific cross-links between chromosomes, HiCUP analysis was performed using the chromosome pseudomolecules as the reference assembly and the proportion of trans linkages between chromosomes of the same or different haplotype computed.</p>
</sec>
<sec id="Sec20"><title>Gene prediction and functional annotation</title>
<p id="Par32">The genome assemblies of Ug99 and <italic>Pgt</italic>
21-0 (as chromosome pseudomolecules for <italic>Pgt</italic>
21-0) were annotated using the Funannotate pipeline (<ext-link ext-link-type="uri" xlink:href="https://github.com/nextgenusfs/funannotate">https://github.com/nextgenusfs/funannotate</ext-link>
). Contigs were sorted by length (longest to shortest) and repetitive elements were soft-masked using RepeatModeler (v1.0.11) and RepeatMasker (v4.0.5; <ext-link ext-link-type="uri" xlink:href="http://www.repeatmasker.org/">http://www.repeatmasker.org/</ext-link>
) with RepBase library (v. 23.09)<sup><xref ref-type="bibr" rid="CR55">55</xref>
</sup>
. RNAseq libraries from <italic>Pgt</italic>
21-0 (Supplementary Data <xref rid="MOESM10" ref-type="media">7</xref>
)<sup><xref ref-type="bibr" rid="CR21">21</xref>
,<xref ref-type="bibr" rid="CR23">23</xref>
</sup>
were used for training gene models. In the training step, RNA-seq data were aligned to the genome assembly with HISAT2<sup><xref ref-type="bibr" rid="CR56">56</xref>
</sup>
. Transcripts were reconstructed with Stringtie (v1.3.4d)<sup><xref ref-type="bibr" rid="CR57">57</xref>
</sup>
. Genome-guided Trinity assembly (v2.4.0)<sup><xref ref-type="bibr" rid="CR58">58</xref>
</sup>
and PASA assembly (v2.3.3)<sup><xref ref-type="bibr" rid="CR59">59</xref>
</sup>
were performed. To assist in predicting effector-like genes, stringtie-aligned transcripts were used in CodingQuarry Pathogen Mode (v2.0)<sup><xref ref-type="bibr" rid="CR60">60</xref>
</sup>
. The prediction step of funannotate pipeline (funannotate predict) was run with --ploidy 2, --busco_db basidiomycota and default parameters. Transcript evidence included Trinity transcripts, Pucciniamycotina EST clusters downloaded from the JGI MycoCosm website (<ext-link ext-link-type="uri" xlink:href="http://genome.jgi.doe.gov/pucciniomycotina/pucciniomycotina.info.html">http://genome.jgi.doe.gov/pucciniomycotina/pucciniomycotina.info.html</ext-link>
, April 24, 2017), and predicted transcript sequences of haustorial secreted proteins<sup><xref ref-type="bibr" rid="CR21">21</xref>
</sup>
. Transcript evidence was aligned to the genome using minimap2 v2.1.0<sup><xref ref-type="bibr" rid="CR61">61</xref>
</sup>
and the protein evidence was aligned to genome via Diamond (v0.9.13)/Exonerate (v2.4.0)<sup><xref ref-type="bibr" rid="CR62">62</xref>
</sup>
using the default UniProtKb/SwissProt curated protein database from funannotate. Ab initio gene predictor AUGUSTUS v3.2.3<sup><xref ref-type="bibr" rid="CR45">45</xref>
</sup>
was trained using PASA data and GeneMark-ES v4.32<sup><xref ref-type="bibr" rid="CR63">63</xref>
</sup>
was self-trained using the genome assembly. Evidence Modeler was used to combine all the above evidence using default weight settings except that the weight of PASA and CodingQuarry_PM were both set to 20. tRNA genes were predicted using tRNAscan-SE v1.3.1<sup><xref ref-type="bibr" rid="CR64">64</xref>
</sup>
. Gene models including UTRs and alternative spliced transcripts were updated using RNAseq data based on Annotation Comparisons and Annotation Updates in PASA. Funannotate fix was run to validate gene models and NCBI submission requirements. Genome annotation was assessed using BUSCO v2.0 (transcript and protein modes)<sup><xref ref-type="bibr" rid="CR46">46</xref>
</sup>
. For functional annotation protein coding gene models were firstly parsed using InterProScan5 (v5.23-62.0) to identify InterPro terms, GO ontology and fungal transcription factors<sup><xref ref-type="bibr" rid="CR65">65</xref>
</sup>
. Pfam domains were identified using PFAM v. 32.0, and carbohydrate hydrolysing enzymatic domains (CAZYmes) were annotated using dbCAN v7.0<sup><xref ref-type="bibr" rid="CR66">66</xref>
</sup>
. Diamond blastP<sup><xref ref-type="bibr" rid="CR67">67</xref>
</sup>
was used to search UniProt DB v. 2018_11<sup><xref ref-type="bibr" rid="CR68">68</xref>
</sup>
and MEROPS v. 12.0<sup><xref ref-type="bibr" rid="CR69">69</xref>
</sup>
databases to aid in functional annotation. BUSCO groups were annotated with Basidiomycota models and eggNOG terms were identified using eggNOG-mapper v1.0.3<sup><xref ref-type="bibr" rid="CR70">70</xref>
</sup>
. Putative mating-type loci in <italic>Pgt</italic>
21-0 and Ug99 were identified by BLAST search with the pheromone peptide encoding genes (<italic>mfa2</italic>
or <italic>mfa3</italic>
) and pheromone mating factor receptors (<italic>STE3.2</italic>
and <italic>STE3.3</italic>
) from the <italic>a</italic>
locus, and <italic>bW</italic>
/<italic>bE</italic>
transcription factors from the <italic>b</italic>
locus previously identified in <italic>Pgt</italic>
<sup><xref ref-type="bibr" rid="CR28">28</xref>
</sup>
. Protein sequences were aligned in Clustal Omega<sup><xref ref-type="bibr" rid="CR71">71</xref>
</sup>
.</p>
<p id="Par33">Gene and repeat density plots for chromosomes were generated using karyoploteR<sup><xref ref-type="bibr" rid="CR72">72</xref>
</sup>
. Secreted proteins were predicted using the neural network predictor of SignalP 3.0<sup><xref ref-type="bibr" rid="CR73">73</xref>
</sup>
and retained if they lacked a transmembrane domain outside the first 60 amino acids using TMHMM<sup><xref ref-type="bibr" rid="CR74">74</xref>
</sup>
. RepeatMasker 4.0.6 with the species fungi was used to softmask repeats. Repeats longer than 200 bp were used in the chromosome plotting.</p>
</sec>
<sec id="Sec21"><title>Orthology analysis</title>
<p id="Par34">Gene annotations with multiple isoforms were reduced to a representative isoform by selecting the longest CDS using a custom perl script. Orthologous proteins were identified with Orthofinder<sup><xref ref-type="bibr" rid="CR75">75</xref>
</sup>
using default parameters. Multiple pairwise orthology analyses were run based on within-isolate and cross-isolate comparisons of similar haplotypes (i.e. A versus A or B versus C). Additional comparisons were made between <italic>Pgt</italic>
21-0 A, <italic>Pgt</italic>
21-0 B, and Ug99 C haplotypes, as well as between Ug99 A, Ug99 C and <italic>Pgt</italic>
21-0 B haplotypes.</p>
</sec>
<sec id="Sec22"><title>Phylogenetic analysis of rust isolates</title>
<p id="Par35">For whole-genome SNP calling and phylogenetic analysis we used Illumina DNA sequence data from the five Ug99 lineage isolates described here, seven Australian isolates<sup><xref ref-type="bibr" rid="CR21">21</xref>
,<xref ref-type="bibr" rid="CR23">23</xref>
</sup>
as well as 31 global isolates<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
downloaded from the European Nucleotide Archive (ENA; PRJEB22223) (Supplementary Data <xref rid="MOESM11" ref-type="media">8</xref>
). Read quality was checked using FASTQC (<ext-link ext-link-type="uri" xlink:href="http://www.bioinformatics.babraham.ac.uk/projects/fastqc/">http://www.bioinformatics.babraham.ac.uk/projects/fastqc/</ext-link>
) and reads were trimmed with Trimmomatic v.0.33<sup><xref ref-type="bibr" rid="CR76">76</xref>
</sup>
using default settings and reads <80 bp were discarded. Trimmed reads were aligned to the Ug99 or <italic>Pgt</italic>
21-0 genome assemblies using BWA-MEM v.0.7.17<sup><xref ref-type="bibr" rid="CR41">41</xref>
</sup>
and technical replicates were merged using SAMtools 1.6<sup><xref ref-type="bibr" rid="CR42">42</xref>
</sup>
and the PICARD toolkit (<ext-link ext-link-type="uri" xlink:href="http://broadinstitute.github.io/picard/">http://broadinstitute.github.io/picard/</ext-link>
). Read lengths and coverage were verified by the functions bamtobed and coverage in BEDtools<sup><xref ref-type="bibr" rid="CR43">43</xref>
</sup>
and flagstat in SAMtools. FreeBayes v. 1.1.0<sup><xref ref-type="bibr" rid="CR51">51</xref>
</sup>
was used to call biallelic SNP variants across the 43 samples simultaneously. VCF files were filtered using vcffilter in vcflib (v1.0.0-rc1) with the parameters: QUAL > 20 & QUAL / AO > 10 & SAF > 0 & SAR > 0 & RPR > 1 & RPL > 1 & AC > 0. To verify that each sample consisted of a single genotype free of contamination, read allele frequencies at heterozygous positions were examined using the vcfR package<sup><xref ref-type="bibr" rid="CR77">77</xref>
</sup>
. VCF files were converted to multiple sequence alignment in PHYLIP format using the vcf2phylip script (<ext-link ext-link-type="uri" xlink:href="https://zenodo.org/record/1257058#.XNnE845Kh3g">https://zenodo.org/record/1257058#.XNnE845Kh3g</ext-link>
) and R-package ips/phyloch wrappings (<ext-link ext-link-type="uri" xlink:href="http://www.christophheibl.de/Rpackages.html">http://www.christophheibl.de/Rpackages.html</ext-link>
). Phylogenetic trees were constructed using the maximum likelihood criterion in RAxML v. 8.2.1.pthread<sup><xref ref-type="bibr" rid="CR78">78</xref>
</sup>
, assuming unlinked loci and using 500 bootstrap replicates with a general time reversible model. Convergence and posterior bootstopping (bootstrapping and convergence criterion) were confirmed with the <italic>-I</italic>
parameter in RAxML and also with R-packages ape<sup><xref ref-type="bibr" rid="CR79">79</xref>
</sup>
, ips/phyloch, and phangorn<sup><xref ref-type="bibr" rid="CR80">80</xref>
</sup>
. Dendrograms were drawn using ggplot<italic>2</italic>
<sup><xref ref-type="bibr" rid="CR81">81</xref>
</sup>
and ggbio<sup><xref ref-type="bibr" rid="CR49">49</xref>
</sup>
R-packages.</p>
<p id="Par36">SNPs called against the A, B or C haplotypes were separated from the total SNP sets using the function intersect –header in BEDtools. Homozygous and heterozygous SNPs were extracted by vcffilter -f “TYPE = snp” and -f “AC > 0 & AC = 2” and -f “AC > 0 & AC = 1”, respectively and their frequency was counted using vcfkeepsamples and vcffixup. SNP statistics were calculated using vcfstats of VCFlib (v1.0.0-rc1).</p>
</sec>
<sec id="Sec23"><title>Reporting summary</title>
<p id="Par37">Further information on research design is available in the <xref rid="MOESM12" ref-type="media">Nature Research Reporting Summary</xref>
linked to this article.</p>
</sec>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary information</title>
<sec id="Sec24"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="41467_2019_12927_MOESM1_ESM.pdf"><caption><p>Supplementary Information</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM2"><media xlink:href="41467_2019_12927_MOESM2_ESM.pdf"><caption><p>Peer Review File</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM3"><media xlink:href="41467_2019_12927_MOESM3_ESM.pdf"><caption><p>Description of Additional Supplementary Files</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM4"><media xlink:href="41467_2019_12927_MOESM4_ESM.xlsx"><caption><p>Supplementary Data 1</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM5"><media xlink:href="41467_2019_12927_MOESM5_ESM.xlsx"><caption><p>Supplementary Data 2</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM6"><media xlink:href="41467_2019_12927_MOESM6_ESM.xlsx"><caption><p>Supplementary Data 3</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM7"><media xlink:href="41467_2019_12927_MOESM7_ESM.xlsx"><caption><p>Supplementary Data 4</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM8"><media xlink:href="41467_2019_12927_MOESM8_ESM.xlsx"><caption><p>Supplementary Data 5</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM9"><media xlink:href="41467_2019_12927_MOESM9_ESM.xlsx"><caption><p>Supplementary Data 6</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM10"><media xlink:href="41467_2019_12927_MOESM10_ESM.xlsx"><caption><p>Supplementary Data 7</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM11"><media xlink:href="41467_2019_12927_MOESM11_ESM.xlsx"><caption><p>Supplementary Data 8</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="MOESM12"><media xlink:href="41467_2019_12927_MOESM12_ESM.pdf"><caption><p>Reporting Summary</p>
</caption>
</media>
</supplementary-material>
</p>
</sec>
</sec>
</body>
<back><fn-group><fn><p><bold>Peer review information</bold>
<italic>Nature Communications</italic>
thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.</p>
</fn>
<fn><p><bold>Publisher’s note</bold>
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
</fn>
<fn><p>These authors contributed equally: Feng Li, Narayana M. Upadhyaya</p>
</fn>
</fn-group>
<sec><title>Supplementary information</title>
<p><bold>Supplementary information</bold>
is available for this paper at 10.1038/s41467-019-12927-7.</p>
</sec>
<ack><title>Acknowledgements</title>
<p>We thank J. Ellis, P. van Esse, G. Bakkeren, C. Aime and Y. Jin for valuable discussions, S. Dahl and N. Prenevost for technical support, J. Palmer for gene annotation troubleshooting, and the Minnesota Supercomputing Institute for computational resources. This research was funded by two independent grants from the 2 Blades foundation to P.N.D. and M.F., respectively, by a USDA-Agriculture and Food Research Initiative (AFRI) Competitive Grant (Proposal No. 2017-08221) to M.F, and University of Minnesota Lieberman-Okinow Endowment to B.J.S.; M.F. and M.E.M were supported by the University of Minnesota Experimental Station USDA-NIFA Hatch/Figueroa project MIN-22-G19 and an USDA-NIFA Postdoctoral Fellowship award (2017-67012-26117), respectively. B.S. is supported by an ARC Future Fellowship (FT180100024). J.S. is supported by an ARC DECRA Fellowship (DE190100066).</p>
</ack>
<notes notes-type="author-contribution"><title>Author contributions</title>
<p>M.F. and P.N.D. conceptualised the project, acquired funding and supervised the work. B.V. and Z.A.P. provided study materials. F.L., N.M.U., C.R., O.M., B.S., R.M. and B.J.S. acquired experimental data. F.L., N.M.U., J.S., B.S., B.J.S., H.N.P., P.N.D, K.A.T.S., E.H., M.E.M. and C.D.H. conducted data analysis. M.F., P.N.D. and F.L. drafted the manuscript. All authors contributed to review and editing.</p>
</notes>
<notes notes-type="data-availability"><title>Data availability</title>
<p>All sequence data, assemblies and gene annotation files generated in this study are available in NCBI under BioProject PRJNA516922. Assemblies and annotations are also available at the DOE-JGI Mycocosm Portal (<ext-link ext-link-type="uri" xlink:href="https://mycocosm.jgi.doe.gov/mycocosm/home">https://mycocosm.jgi.doe.gov/mycocosm/home</ext-link>
). Metadata for RNAseq libraries of <italic>Pgt</italic>
21-0 and Illumina DNAseq libraries from all isolates are available in Supplementary Data <xref rid="MOESM10" ref-type="media">7</xref>
and <xref rid="MOESM11" ref-type="media">8</xref>
, respectively. Data underlying Fig. <xref rid="Fig3" ref-type="fig">3c</xref>
and Supplementary Fig. <xref rid="MOESM1" ref-type="media">3a, b</xref>
is provided in Supplementary Data <xref rid="MOESM6" ref-type="media">3</xref>
. Data underlying Fig. <xref rid="Fig6" ref-type="fig">6</xref>
is available in Supplementary Data <xref rid="MOESM9" ref-type="media">6</xref>
. All other relevant data is available upon request from the corresponding authors.</p>
</notes>
<notes notes-type="data-availability"><title>Code availability</title>
<p>Unless specified otherwise, all scripts and files are available at <ext-link ext-link-type="uri" xlink:href="https://github.com/figueroalab/Pgt_genomes">https://github.com/figueroalab/Pgt_genomes</ext-link>
.</p>
</notes>
<notes notes-type="COI-statement"><title>Competing interests</title>
<p id="Par38">The authors declare no competing interests.</p>
</notes>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stajich</surname>
<given-names>JE</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Primer-the fungi</article-title>
<source>Curr. Biol.</source>
<year>2009</year>
<volume>19</volume>
<fpage>R840</fpage>
<pub-id pub-id-type="doi">10.1016/j.cub.2009.07.004</pub-id>
<pub-id pub-id-type="pmid">19788875</pub-id>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<mixed-citation publication-type="other">Fleißner, A., Serrano, A. <italic>Growth, Differentiation and Sexuality</italic>
(eds Wendland, J.) (Springer, 2016).</mixed-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Spatafora</surname>
<given-names>JW</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The fungal tree of life: from molecular systematics to genome-scale phylogenies</article-title>
<source>Microbiol. Spectr.</source>
<year>2017</year>
<volume>5</volume>
<fpage>1</fpage>
<lpage>32</lpage>
<pub-id pub-id-type="doi">10.1128/microbiolspec.FUNK-0053-2016</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aime</surname>
<given-names>MC</given-names>
</name>
<name><surname>McTaggart</surname>
<given-names>AR</given-names>
</name>
<name><surname>Mondo</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Duplessis</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Phylogenetics and phylogenomics of rust fungi</article-title>
<source>Adv. Genet.</source>
<year>2017</year>
<volume>100</volume>
<fpage>267</fpage>
<lpage>307</lpage>
<pub-id pub-id-type="doi">10.1016/bs.adgen.2017.09.011</pub-id>
<pub-id pub-id-type="pmid">29153402</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Watson</surname>
<given-names>I</given-names>
</name>
</person-group>
<article-title>Further studies on the production of new races from mixtures of races of <italic>Puccinia graminis</italic>
var. <italic>tritici</italic>
on wheat seedlings</article-title>
<source>Phytopathology</source>
<year>1957</year>
<volume>47</volume>
<fpage>510</fpage>
<lpage>512</lpage>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ellingboe</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Somatic recombination in <italic>Puccinia graminis</italic>
var <italic>tritici</italic>
</article-title>
<source>Phytopathology</source>
<year>1961</year>
<volume>51</volume>
<fpage>13</fpage>
<lpage>15</lpage>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Flor</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Genetics of somatic variation for pathogenicity in <italic>Melampsora lini</italic>
</article-title>
<source>Phytopathology</source>
<year>1964</year>
<volume>54</volume>
<fpage>823</fpage>
<lpage>826</lpage>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bartos</surname>
<given-names>P</given-names>
</name>
<name><surname>Fleischmann</surname>
<given-names>G</given-names>
</name>
<name><surname>Samborski</surname>
<given-names>D</given-names>
</name>
<name><surname>Shipton</surname>
<given-names>W</given-names>
</name>
</person-group>
<article-title>Studies on asexual variation in the virulence of oat crown rust, <italic>Puccinia coronata</italic>
f. sp. <italic>avenae</italic>
, and wheat leaf rust, <italic>Puccinia recondita</italic>
</article-title>
<source>Can. J. Bot.</source>
<year>1969</year>
<volume>47</volume>
<fpage>1383</fpage>
<lpage>1387</lpage>
<pub-id pub-id-type="doi">10.1139/b69-197</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Burdon</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Marshall</surname>
<given-names>DR</given-names>
</name>
<name><surname>Luig</surname>
<given-names>NH</given-names>
</name>
</person-group>
<article-title>Isozyme analysis indicates that a virulent cereal rust pathogen is a somatic hybrid</article-title>
<source>Nature</source>
<year>1981</year>
<volume>293</volume>
<fpage>565</fpage>
<pub-id pub-id-type="doi">10.1038/293565a0</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname>
<given-names>R</given-names>
</name>
<name><surname>Burdon</surname>
<given-names>J</given-names>
</name>
<name><surname>Jahoor</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Evidence for somatic hybridisation in the leaf rust pathogen of wheat (<italic>Puccinia recondita</italic>
f. sp. <italic>tritici</italic>
)</article-title>
<source>Mycological Res.</source>
<year>1999</year>
<volume>103</volume>
<fpage>715</fpage>
<lpage>723</lpage>
<pub-id pub-id-type="doi">10.1017/S0953756298007631</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname>
<given-names>RF</given-names>
</name>
<name><surname>Wellings</surname>
<given-names>CR</given-names>
</name>
</person-group>
<article-title>Somatic hybridization in the Uredinales</article-title>
<source>Annu. Rev. Phytopathol.</source>
<year>2012</year>
<volume>50</volume>
<fpage>219</fpage>
<lpage>239</lpage>
<pub-id pub-id-type="doi">10.1146/annurev-phyto-072910-095405</pub-id>
<pub-id pub-id-type="pmid">22920559</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname>
<given-names>ME</given-names>
</name>
<etal></etal>
</person-group>
<article-title>De novo assembly and phasing of dikaryotic genomes from two isolates of <italic>Puccinia coronata</italic>
f. sp. <italic>avenae</italic>
, the causal agent of oat crown rust</article-title>
<source>mBio</source>
<year>2018</year>
<volume>9</volume>
<fpage>e01650</fpage>
<lpage>01617</lpage>
<pub-id pub-id-type="pmid">29463655</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schwessinger</surname>
<given-names>B</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A near-complete haplotype-phased genome of the dikaryotic wheat stripe rust fungus <italic>Puccinia striiformis</italic>
f. sp. <italic>tritici</italic>
reveals high interhaplotype diversity</article-title>
<source>mBio</source>
<year>2018</year>
<volume>9</volume>
<fpage>e02275</fpage>
<lpage>02217</lpage>
<pub-id pub-id-type="doi">10.1128/mBio.02275-17</pub-id>
<pub-id pub-id-type="pmid">29463659</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Singh</surname>
<given-names>RP</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Emergence and spread of new races of wheat stem rust fungus: continued threat to food security and prospects of genetic control</article-title>
<source>Phytopathology</source>
<year>2015</year>
<volume>105</volume>
<fpage>872</fpage>
<lpage>884</lpage>
<pub-id pub-id-type="doi">10.1094/PHYTO-01-15-0030-FI</pub-id>
<pub-id pub-id-type="pmid">26120730</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pretorius</surname>
<given-names>Z</given-names>
</name>
<name><surname>Singh</surname>
<given-names>R</given-names>
</name>
<name><surname>Wagoire</surname>
<given-names>W</given-names>
</name>
<name><surname>Payne</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Detection of virulence to wheat stem rust resistance gene <italic>Sr31</italic>
in <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
in Uganda</article-title>
<source>Plant Dis.</source>
<year>2000</year>
<volume>84</volume>
<fpage>203</fpage>
<lpage>203</lpage>
<pub-id pub-id-type="doi">10.1094/PDIS.2000.84.2.203B</pub-id>
<pub-id pub-id-type="pmid">30841334</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Visser</surname>
<given-names>B</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Characterization of two new <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
races within the Ug99 lineage in South Africa</article-title>
<source>Euphytica</source>
<year>2011</year>
<volume>179</volume>
<fpage>E119</fpage>
<lpage>E127</lpage>
<pub-id pub-id-type="doi">10.1007/s10681-010-0269-x</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Olivera</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Phenotypic and genotypic characterization of race TKTTF of <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
that caused a wheat stem rust epidemic in southern Ethiopia in 2013–14</article-title>
<source>Phytopathology</source>
<year>2015</year>
<volume>105</volume>
<fpage>917</fpage>
<lpage>928</lpage>
<pub-id pub-id-type="doi">10.1094/PHYTO-11-14-0302-FI</pub-id>
<pub-id pub-id-type="pmid">25775107</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Visser</surname>
<given-names>B</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Microsatellite analysis and urediniospore dispersal simulations support the movement of <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
from Southern Africa to Australia</article-title>
<source>Phytopathology</source>
<year>2019</year>
<volume>109</volume>
<fpage>133</fpage>
<lpage>144</lpage>
<pub-id pub-id-type="doi">10.1094/PHYTO-04-18-0110-R</pub-id>
<pub-id pub-id-type="pmid">30028232</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pretorius</surname>
<given-names>ZA</given-names>
</name>
<name><surname>Pakendorf</surname>
<given-names>KW</given-names>
</name>
<name><surname>Marais</surname>
<given-names>GF</given-names>
</name>
<name><surname>Prins</surname>
<given-names>R</given-names>
</name>
<name><surname>Komen</surname>
<given-names>JS</given-names>
</name>
</person-group>
<article-title>Challenges for sustainable cereal rust control in South Africa</article-title>
<source>Aust. J. Agric. Res.</source>
<year>2007</year>
<volume>58</volume>
<fpage>593</fpage>
<lpage>601</lpage>
<pub-id pub-id-type="doi">10.1071/AR06144</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname>
<given-names>RF</given-names>
</name>
</person-group>
<article-title>Stem rust of wheat in Australia</article-title>
<source>Aust. J. Agric. Res.</source>
<year>2007</year>
<volume>58</volume>
<fpage>558</fpage>
<lpage>566</lpage>
<pub-id pub-id-type="doi">10.1071/AR07117</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Upadhyaya</surname>
<given-names>NM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Comparative genomics of Australian isolates of the wheat stem rust pathogen <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
reveals extensive polymorphism in candidate effector genes</article-title>
<source>Front. Plant Sci.</source>
<year>2015</year>
<volume>5</volume>
<fpage>759</fpage>
<pub-id pub-id-type="doi">10.3389/fpls.2014.00759</pub-id>
<pub-id pub-id-type="pmid">25620970</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duplessis</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Obligate biotrophy features unraveled by the genomic analysis of rust fungi</article-title>
<source>Proc. Natl Acad. Sci. USA</source>
<year>2011</year>
<volume>108</volume>
<fpage>9166</fpage>
<lpage>9171</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.1019315108</pub-id>
<pub-id pub-id-type="pmid">21536894</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Loss of AvrSr50 by somatic exchange in stem rust leads to virulence for Sr50 resistance in wheat</article-title>
<source>Science</source>
<year>2017</year>
<volume>358</volume>
<fpage>1607</fpage>
<lpage>1610</lpage>
<pub-id pub-id-type="doi">10.1126/science.aao4810</pub-id>
<pub-id pub-id-type="pmid">29269475</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salcedo</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Variation in the <italic>AvrSr35</italic>
gene determines <italic>Sr35</italic>
resistance against wheat stem rust race Ug99</article-title>
<source>Science</source>
<year>2017</year>
<volume>358</volume>
<fpage>1604</fpage>
<lpage>1606</lpage>
<pub-id pub-id-type="doi">10.1126/science.aao7294</pub-id>
<pub-id pub-id-type="pmid">29269474</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname>
<given-names>C</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genome analysis and avirulence gene cloning using a high-density RADseq linkage map of the flax rust fungus, <italic>Melampsora lini</italic>
</article-title>
<source>BMC Genomics</source>
<year>2016</year>
<volume>17</volume>
<fpage>667</fpage>
<pub-id pub-id-type="doi">10.1186/s12864-016-3011-9</pub-id>
<pub-id pub-id-type="pmid">27550217</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Belton</surname>
<given-names>JM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Hi-C: a comprehensive technique to capture the conformation of genomes</article-title>
<source>Methods</source>
<year>2012</year>
<volume>58</volume>
<fpage>268</fpage>
<lpage>276</lpage>
<pub-id pub-id-type="doi">10.1016/j.ymeth.2012.05.001</pub-id>
<pub-id pub-id-type="pmid">22652625</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boehm</surname>
<given-names>EWA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>An ultrastructural pachytene karyotype for <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
</article-title>
<source>Can. J. Bot.</source>
<year>1992</year>
<volume>70</volume>
<fpage>401</fpage>
<lpage>413</lpage>
<pub-id pub-id-type="doi">10.1139/b92-054</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cuomo</surname>
<given-names>CA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Comparative analysis highlights variable genome content of wheat rusts and divergence of the mating loci. <italic>G3: Genes, Genomes</italic>
</article-title>
<source>Genetics</source>
<year>2017</year>
<volume>7</volume>
<fpage>361</fpage>
<lpage>376</lpage>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lewis</surname>
<given-names>CM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Potential for re-emergence of wheat stem rust in the United Kingdom</article-title>
<source>Commun. Biol.</source>
<year>2018</year>
<volume>1</volume>
<fpage>13</fpage>
<pub-id pub-id-type="doi">10.1038/s42003-018-0013-y</pub-id>
<pub-id pub-id-type="pmid">30271900</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Figueroa</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Changing the game: using integrative genomics to probe virulence mechanisms of the stem rust pathogen <italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
</article-title>
<source>Front. Plant Sci.</source>
<year>2016</year>
<volume>7</volume>
<fpage>205</fpage>
<pub-id pub-id-type="doi">10.3389/fpls.2016.00205</pub-id>
<pub-id pub-id-type="pmid">26941766</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lei</surname>
<given-names>Y</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Virulence and molecular characterization of experimental isolates of the stripe rust pathogen (Puccinia striiformis) indicate somatic recombination</article-title>
<source>Phytopathology</source>
<year>2017</year>
<volume>107</volume>
<fpage>329</fpage>
<lpage>344</lpage>
<pub-id pub-id-type="doi">10.1094/PHYTO-07-16-0261-R</pub-id>
<pub-id pub-id-type="pmid">27775498</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barr</surname>
<given-names>R</given-names>
</name>
<name><surname>Caldwell</surname>
<given-names>R</given-names>
</name>
<name><surname>Amacher</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Examination of vegetative recombination of urediospore color and virulence in mixtures of certain races of <italic>Puccinia recondita</italic>
</article-title>
<source>Phytopathology</source>
<year>1964</year>
<volume>54</volume>
<fpage>104</fpage>
<lpage>109</lpage>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>X</given-names>
</name>
<name><surname>McCallum</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Fusion body formation, germ tube anastomosis, and nuclear migration during the germination of urediniospores of the wheat leaf rust fungus, <italic>Puccinia triticina</italic>
</article-title>
<source>Phytopathology</source>
<year>2009</year>
<volume>99</volume>
<fpage>1355</fpage>
<lpage>1364</lpage>
<pub-id pub-id-type="doi">10.1094/PHYTO-99-12-1355</pub-id>
<pub-id pub-id-type="pmid">19900001</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kruzel</surname>
<given-names>EK</given-names>
</name>
<name><surname>Hull</surname>
<given-names>CM</given-names>
</name>
</person-group>
<article-title>Establishing an unusual cell type: how to make a dikaryon</article-title>
<source>Curr. Opin. Microbiol.</source>
<year>2010</year>
<volume>13</volume>
<fpage>706</fpage>
<lpage>711</lpage>
<pub-id pub-id-type="doi">10.1016/j.mib.2010.09.016</pub-id>
<pub-id pub-id-type="pmid">21036099</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Clark</surname>
<given-names>TA</given-names>
</name>
<name><surname>Anderson</surname>
<given-names>JB</given-names>
</name>
</person-group>
<article-title>Dikaryons of the basidiomycete fungus <italic>Schizophyllum commune</italic>
: evolution in long-term culture</article-title>
<source>Genetics</source>
<year>2004</year>
<volume>167</volume>
<fpage>1663</fpage>
<lpage>1675</lpage>
<pub-id pub-id-type="doi">10.1534/genetics.104.027235</pub-id>
<pub-id pub-id-type="pmid">15342506</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Corradi</surname>
<given-names>N</given-names>
</name>
<name><surname>Brachmann</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Fungal mating in the most widespread plant symbionts?</article-title>
<source>Trends Plant Sci.</source>
<year>2017</year>
<volume>22</volume>
<fpage>175</fpage>
<lpage>183</lpage>
<pub-id pub-id-type="doi">10.1016/j.tplants.2016.10.010</pub-id>
<pub-id pub-id-type="pmid">27876487</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pretorius</surname>
<given-names>Z</given-names>
</name>
<name><surname>Szabo</surname>
<given-names>L</given-names>
</name>
<name><surname>Boshoff</surname>
<given-names>W</given-names>
</name>
<name><surname>Herselman</surname>
<given-names>L</given-names>
</name>
<name><surname>Visser</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>First report of a new TTKSF race of wheat stem rust (<italic>Puccinia graminis</italic>
f. sp. <italic>tritici</italic>
) in South Africa and Zimbabwe</article-title>
<source>Plant Dis.</source>
<year>2012</year>
<volume>96</volume>
<fpage>590</fpage>
<lpage>590</lpage>
<pub-id pub-id-type="doi">10.1094/PDIS-12-11-1027-PDN</pub-id>
<pub-id pub-id-type="pmid">30727416</pub-id>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<mixed-citation publication-type="other">Schwessinger, B., Rathjen, J. P. <italic>Wheat Rust Diseases</italic>
(eds Periyannan, S.) (Springer, 2017).</mixed-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koren</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation</article-title>
<source>Genome Res.</source>
<year>2017</year>
<volume>27</volume>
<fpage>722</fpage>
<lpage>736</lpage>
<pub-id pub-id-type="doi">10.1101/gr.215087.116</pub-id>
<pub-id pub-id-type="pmid">28298431</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Walker</surname>
<given-names>BJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement</article-title>
<source>PLoS One</source>
<year>2014</year>
<volume>9</volume>
<fpage>e112963</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0112963</pub-id>
<pub-id pub-id-type="pmid">25409509</pub-id>
</element-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
<name><surname>Durbin</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Fast and accurate short read alignment with Burrows–Wheeler transform</article-title>
<source>Bioinformatics</source>
<year>2009</year>
<volume>25</volume>
<fpage>1754</fpage>
<lpage>1760</lpage>
<pub-id pub-id-type="pmid">19451168</pub-id>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The sequence alignment/map format and SAMtools</article-title>
<source>Bioinformatics</source>
<year>2009</year>
<volume>25</volume>
<fpage>2078</fpage>
<lpage>2079</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btp352</pub-id>
<pub-id pub-id-type="pmid">19505943</pub-id>
</element-citation>
</ref>
<ref id="CR43"><label>43.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Quinlan</surname>
<given-names>AR</given-names>
</name>
<name><surname>Hall</surname>
<given-names>IM</given-names>
</name>
</person-group>
<article-title>BEDTools: a flexible suite of utilities for comparing genomic features</article-title>
<source>Bioinformatics</source>
<year>2010</year>
<volume>26</volume>
<fpage>841</fpage>
<lpage>842</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btq033</pub-id>
<pub-id pub-id-type="pmid">20110278</pub-id>
</element-citation>
</ref>
<ref id="CR44"><label>44.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gurevich</surname>
<given-names>A</given-names>
</name>
<name><surname>Saveliev</surname>
<given-names>V</given-names>
</name>
<name><surname>Vyahhi</surname>
<given-names>N</given-names>
</name>
<name><surname>Tesler</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>QUAST: quality assessment tool for genome assemblies</article-title>
<source>Bioinformatics</source>
<year>2013</year>
<volume>29</volume>
<fpage>1072</fpage>
<lpage>1075</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btt086</pub-id>
<pub-id pub-id-type="pmid">23422339</pub-id>
</element-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stanke</surname>
<given-names>M</given-names>
</name>
<name><surname>Morgenstern</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints</article-title>
<source>Nucleic Acids Res.</source>
<year>2005</year>
<volume>33</volume>
<fpage>W465</fpage>
<lpage>W467</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gki458</pub-id>
<pub-id pub-id-type="pmid">15980513</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simão</surname>
<given-names>FA</given-names>
</name>
<name><surname>Waterhouse</surname>
<given-names>RM</given-names>
</name>
<name><surname>Ioannidis</surname>
<given-names>P</given-names>
</name>
<name><surname>Kriventseva</surname>
<given-names>EV</given-names>
</name>
<name><surname>Zdobnov</surname>
<given-names>EM</given-names>
</name>
</person-group>
<article-title>BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs</article-title>
<source>Bioinformatics</source>
<year>2015</year>
<volume>31</volume>
<fpage>3210</fpage>
<lpage>3212</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btv351</pub-id>
<pub-id pub-id-type="pmid">26059717</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cabanettes</surname>
<given-names>F</given-names>
</name>
<name><surname>Klopp</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>D-GENIES: dot plot large genomes in an interactive, efficient and simple way</article-title>
<source>PeerJ</source>
<year>2018</year>
<volume>6</volume>
<fpage>e4958</fpage>
<pub-id pub-id-type="doi">10.7717/peerj.4958</pub-id>
<pub-id pub-id-type="pmid">29888139</pub-id>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lawrence</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Software for computing and annotating genomic ranges</article-title>
<source>PLoS Comput. Biol.</source>
<year>2013</year>
<volume>9</volume>
<fpage>e1003118</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1003118</pub-id>
<pub-id pub-id-type="pmid">23950696</pub-id>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yin</surname>
<given-names>T</given-names>
</name>
<name><surname>Cook</surname>
<given-names>D</given-names>
</name>
<name><surname>Lawrence</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>ggbio: an R package for extending the grammar of graphics for genomic data</article-title>
<source>Genome Biol.</source>
<year>2012</year>
<volume>13</volume>
<fpage>R77</fpage>
<pub-id pub-id-type="doi">10.1186/gb-2012-13-8-r77</pub-id>
<pub-id pub-id-type="pmid">22937822</pub-id>
</element-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kurtz</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Versatile and open software for comparing large genomes</article-title>
<source>Genome Biol.</source>
<year>2004</year>
<volume>5</volume>
<fpage>R12</fpage>
<pub-id pub-id-type="doi">10.1186/gb-2004-5-2-r12</pub-id>
<pub-id pub-id-type="pmid">14759262</pub-id>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garrison</surname>
<given-names>E</given-names>
</name>
<name><surname>Marth</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Haplotype-based variant detection from short-read sequencing</article-title>
<source>arXiv</source>
<year>2012</year>
<volume>1207</volume>
<fpage>3907</fpage>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wingett</surname>
<given-names>Steven W.</given-names>
</name>
<name><surname>Ewels</surname>
<given-names>Philip</given-names>
</name>
<name><surname>Furlan-Magaril</surname>
<given-names>Mayra</given-names>
</name>
<name><surname>Nagano</surname>
<given-names>Takashi</given-names>
</name>
<name><surname>Schoenfelder</surname>
<given-names>Stefan</given-names>
</name>
<name><surname>Fraser</surname>
<given-names>Peter</given-names>
</name>
<name><surname>Andrews</surname>
<given-names>Simon</given-names>
</name>
</person-group>
<article-title>HiCUP: pipeline for mapping and processing Hi-C data</article-title>
<source>F1000Research</source>
<year>2015</year>
<volume>4</volume>
<fpage>1310</fpage>
<pub-id pub-id-type="doi">10.12688/f1000research.7334.1</pub-id>
<pub-id pub-id-type="pmid">26835000</pub-id>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langmead</surname>
<given-names>B</given-names>
</name>
<name><surname>Salzberg</surname>
<given-names>SL</given-names>
</name>
</person-group>
<article-title>Fast gapped-read alignment with Bowtie 2</article-title>
<source>Nat. Methods</source>
<year>2012</year>
<volume>9</volume>
<fpage>357</fpage>
<pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
<pub-id pub-id-type="pmid">22388286</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghurye</surname>
<given-names>J</given-names>
</name>
<name><surname>Pop</surname>
<given-names>M</given-names>
</name>
<name><surname>Koren</surname>
<given-names>S</given-names>
</name>
<name><surname>Bickhart</surname>
<given-names>D</given-names>
</name>
<name><surname>Chin</surname>
<given-names>C-S</given-names>
</name>
</person-group>
<article-title>Scaffolding of long read assemblies using long range contact information</article-title>
<source>BMC Genomics</source>
<year>2017</year>
<volume>18</volume>
<fpage>527</fpage>
<pub-id pub-id-type="doi">10.1186/s12864-017-3879-z</pub-id>
<pub-id pub-id-type="pmid">28701198</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname>
<given-names>W</given-names>
</name>
<name><surname>Kojima</surname>
<given-names>KK</given-names>
</name>
<name><surname>Kohany</surname>
<given-names>O</given-names>
</name>
</person-group>
<article-title>Repbase Update, a database of repetitive elements in eukaryotic genomes</article-title>
<source>Mob. DNA</source>
<year>2015</year>
<volume>6</volume>
<fpage>11</fpage>
<pub-id pub-id-type="doi">10.1186/s13100-015-0041-9</pub-id>
<pub-id pub-id-type="pmid">26045719</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname>
<given-names>D</given-names>
</name>
<name><surname>Langmead</surname>
<given-names>B</given-names>
</name>
<name><surname>Salzberg</surname>
<given-names>SL</given-names>
</name>
</person-group>
<article-title>HISAT: a fast spliced aligner with low memory requirements</article-title>
<source>Nat. Methods</source>
<year>2015</year>
<volume>12</volume>
<fpage>357</fpage>
<pub-id pub-id-type="doi">10.1038/nmeth.3317</pub-id>
<pub-id pub-id-type="pmid">25751142</pub-id>
</element-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pertea</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>StringTie enables improved reconstruction of a transcriptome from RNA-seq reads</article-title>
<source>Nat. Biotechnol.</source>
<year>2015</year>
<volume>33</volume>
<fpage>290</fpage>
<pub-id pub-id-type="doi">10.1038/nbt.3122</pub-id>
<pub-id pub-id-type="pmid">25690850</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grabherr</surname>
<given-names>MG</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data</article-title>
<source>Nat. Biotechnol.</source>
<year>2011</year>
<volume>29</volume>
<fpage>644</fpage>
<lpage>652</lpage>
<pub-id pub-id-type="doi">10.1038/nbt.1883</pub-id>
<pub-id pub-id-type="pmid">21572440</pub-id>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haas</surname>
<given-names>BJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments</article-title>
<source>Genome Biol.</source>
<year>2008</year>
<volume>9</volume>
<fpage>R7</fpage>
<pub-id pub-id-type="doi">10.1186/gb-2008-9-1-r7</pub-id>
<pub-id pub-id-type="pmid">18190707</pub-id>
</element-citation>
</ref>
<ref id="CR60"><label>60.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Testa</surname>
<given-names>AC</given-names>
</name>
<name><surname>Hane</surname>
<given-names>JK</given-names>
</name>
<name><surname>Ellwood</surname>
<given-names>SR</given-names>
</name>
<name><surname>Oliver</surname>
<given-names>RP</given-names>
</name>
</person-group>
<article-title>CodingQuarry: highly accurate hidden Markov model gene prediction in fungal genomes using RNA-seq transcripts</article-title>
<source>BMC Genomics</source>
<year>2015</year>
<volume>16</volume>
<fpage>170</fpage>
<pub-id pub-id-type="doi">10.1186/s12864-015-1344-4</pub-id>
<pub-id pub-id-type="pmid">25887563</pub-id>
</element-citation>
</ref>
<ref id="CR61"><label>61.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Minimap2: pairwise alignment for nucleotide sequences</article-title>
<source>Bioinformatics</source>
<year>2018</year>
<volume>34</volume>
<fpage>3094</fpage>
<lpage>3100</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/bty191</pub-id>
<pub-id pub-id-type="pmid">29750242</pub-id>
</element-citation>
</ref>
<ref id="CR62"><label>62.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Slater</surname>
<given-names>GSC</given-names>
</name>
<name><surname>Birney</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>Automated generation of heuristics for biological sequence comparison</article-title>
<source>BMC Bioinform.</source>
<year>2005</year>
<volume>6</volume>
<fpage>31</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2105-6-31</pub-id>
</element-citation>
</ref>
<ref id="CR63"><label>63.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lomsadze</surname>
<given-names>A</given-names>
</name>
<name><surname>Ter-Hovhannisyan</surname>
<given-names>V</given-names>
</name>
<name><surname>Chernoff</surname>
<given-names>YO</given-names>
</name>
<name><surname>Borodovsky</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Gene identification in novel eukaryotic genomes by self-training algorithm</article-title>
<source>Nucleic Acids Res.</source>
<year>2005</year>
<volume>33</volume>
<fpage>6494</fpage>
<lpage>6506</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gki937</pub-id>
<pub-id pub-id-type="pmid">16314312</pub-id>
</element-citation>
</ref>
<ref id="CR64"><label>64.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowe</surname>
<given-names>TM</given-names>
</name>
<name><surname>Chan</surname>
<given-names>PP</given-names>
</name>
</person-group>
<article-title>tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes</article-title>
<source>Nucleic Acids Res.</source>
<year>2016</year>
<volume>44</volume>
<fpage>W54</fpage>
<lpage>W57</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkw413</pub-id>
<pub-id pub-id-type="pmid">27174935</pub-id>
</element-citation>
</ref>
<ref id="CR65"><label>65.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jones</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<article-title>InterProScan 5: genome-scale protein function classification</article-title>
<source>Bioinformatics</source>
<year>2014</year>
<volume>30</volume>
<fpage>1236</fpage>
<lpage>1240</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btu031</pub-id>
<pub-id pub-id-type="pmid">24451626</pub-id>
</element-citation>
</ref>
<ref id="CR66"><label>66.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yin</surname>
<given-names>Y</given-names>
</name>
<etal></etal>
</person-group>
<article-title>dbCAN: a web resource for automated carbohydrate-active enzyme annotation</article-title>
<source>Nucleic Acids Res.</source>
<year>2012</year>
<volume>40</volume>
<fpage>W445</fpage>
<lpage>W451</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gks479</pub-id>
<pub-id pub-id-type="pmid">22645317</pub-id>
</element-citation>
</ref>
<ref id="CR67"><label>67.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buchfink</surname>
<given-names>B</given-names>
</name>
<name><surname>Xie</surname>
<given-names>C</given-names>
</name>
<name><surname>Huson</surname>
<given-names>DH</given-names>
</name>
</person-group>
<article-title>Fast and sensitive protein alignment using DIAMOND</article-title>
<source>Nat. Methods</source>
<year>2015</year>
<volume>12</volume>
<fpage>59</fpage>
<pub-id pub-id-type="doi">10.1038/nmeth.3176</pub-id>
<pub-id pub-id-type="pmid">25402007</pub-id>
</element-citation>
</ref>
<ref id="CR68"><label>68.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Consortium</surname>
<given-names>U</given-names>
</name>
</person-group>
<article-title>UniProt: the universal protein knowledgebase</article-title>
<source>Nucleic Acids Res.</source>
<year>2018</year>
<volume>46</volume>
<fpage>2699</fpage>
<pub-id pub-id-type="doi">10.1093/nar/gky092</pub-id>
<pub-id pub-id-type="pmid">29425356</pub-id>
</element-citation>
</ref>
<ref id="CR69"><label>69.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rawlings</surname>
<given-names>ND</given-names>
</name>
</person-group>
<article-title>Peptidase specificity from the substrate cleavage collection in the MEROPS database and a tool to measure cleavage site conservation</article-title>
<source>Biochimie</source>
<year>2016</year>
<volume>122</volume>
<fpage>5</fpage>
<lpage>30</lpage>
<pub-id pub-id-type="doi">10.1016/j.biochi.2015.10.003</pub-id>
<pub-id pub-id-type="pmid">26455268</pub-id>
</element-citation>
</ref>
<ref id="CR70"><label>70.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huerta-Cepas</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<article-title>eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences</article-title>
<source>Nucleic Acids Res.</source>
<year>2015</year>
<volume>44</volume>
<fpage>D286</fpage>
<lpage>D293</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkv1248</pub-id>
<pub-id pub-id-type="pmid">26582926</pub-id>
</element-citation>
</ref>
<ref id="CR71"><label>71.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chojnacki</surname>
<given-names>S</given-names>
</name>
<name><surname>Cowley</surname>
<given-names>A</given-names>
</name>
<name><surname>Lee</surname>
<given-names>J</given-names>
</name>
<name><surname>Foix</surname>
<given-names>A</given-names>
</name>
<name><surname>Lopez</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Programmatic access to bioinformatics tools from EMBL-EBI update: 2017</article-title>
<source>Nucleic Acids Res.</source>
<year>2017</year>
<volume>45</volume>
<fpage>W550</fpage>
<lpage>W553</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkx273</pub-id>
<pub-id pub-id-type="pmid">28431173</pub-id>
</element-citation>
</ref>
<ref id="CR72"><label>72.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gel</surname>
<given-names>B</given-names>
</name>
<name><surname>Serra</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>karyoploteR: an R/Bioconductor package to plot customizable genomes displaying arbitrary data</article-title>
<source>Bioinformatics</source>
<year>2017</year>
<volume>33</volume>
<fpage>3088</fpage>
<lpage>3090</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btx346</pub-id>
<pub-id pub-id-type="pmid">28575171</pub-id>
</element-citation>
</ref>
<ref id="CR73"><label>73.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bendtsen</surname>
<given-names>JD</given-names>
</name>
<name><surname>Nielsen</surname>
<given-names>H</given-names>
</name>
<name><surname>von Heijne</surname>
<given-names>G</given-names>
</name>
<name><surname>Brunak</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Improved prediction of signal peptides: SignalP 3.0</article-title>
<source>J. Mol. Biol.</source>
<year>2004</year>
<volume>340</volume>
<fpage>783</fpage>
<lpage>795</lpage>
<pub-id pub-id-type="doi">10.1016/j.jmb.2004.05.028</pub-id>
<pub-id pub-id-type="pmid">15223320</pub-id>
</element-citation>
</ref>
<ref id="CR74"><label>74.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Krogh</surname>
<given-names>A</given-names>
</name>
<name><surname>Larsson</surname>
<given-names>B</given-names>
</name>
<name><surname>Von Heijne</surname>
<given-names>G</given-names>
</name>
<name><surname>Sonnhammer</surname>
<given-names>EL</given-names>
</name>
</person-group>
<article-title>Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes</article-title>
<source>J. Mol. Biol.</source>
<year>2001</year>
<volume>305</volume>
<fpage>567</fpage>
<lpage>580</lpage>
<pub-id pub-id-type="doi">10.1006/jmbi.2000.4315</pub-id>
<pub-id pub-id-type="pmid">11152613</pub-id>
</element-citation>
</ref>
<ref id="CR75"><label>75.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Emms</surname>
<given-names>DM</given-names>
</name>
<name><surname>Kelly</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy</article-title>
<source>Genome Biol.</source>
<year>2015</year>
<volume>16</volume>
<fpage>157</fpage>
<pub-id pub-id-type="doi">10.1186/s13059-015-0721-2</pub-id>
<pub-id pub-id-type="pmid">26243257</pub-id>
</element-citation>
</ref>
<ref id="CR76"><label>76.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bolger</surname>
<given-names>AM</given-names>
</name>
<name><surname>Lohse</surname>
<given-names>M</given-names>
</name>
<name><surname>Usadel</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Trimmomatic: a flexible trimmer for Illumina sequence data</article-title>
<source>Bioinformatics</source>
<year>2014</year>
<volume>30</volume>
<fpage>2114</fpage>
<lpage>2120</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btu170</pub-id>
<pub-id pub-id-type="pmid">24695404</pub-id>
</element-citation>
</ref>
<ref id="CR77"><label>77.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Knaus</surname>
<given-names>BJ</given-names>
</name>
<name><surname>Grünwald</surname>
<given-names>NJ</given-names>
</name>
</person-group>
<article-title>vcfr: a package to manipulate and visualize variant call format data in R</article-title>
<source>Mol. Ecol. Resour.</source>
<year>2017</year>
<volume>17</volume>
<fpage>44</fpage>
<lpage>53</lpage>
<pub-id pub-id-type="doi">10.1111/1755-0998.12549</pub-id>
<pub-id pub-id-type="pmid">27401132</pub-id>
</element-citation>
</ref>
<ref id="CR78"><label>78.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models</article-title>
<source>Bioinformatics</source>
<year>2006</year>
<volume>22</volume>
<fpage>2688</fpage>
<lpage>2690</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btl446</pub-id>
<pub-id pub-id-type="pmid">16928733</pub-id>
</element-citation>
</ref>
<ref id="CR79"><label>79.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Paradis</surname>
<given-names>E</given-names>
</name>
<name><surname>Claude</surname>
<given-names>J</given-names>
</name>
<name><surname>Strimmer</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>APE: analyses of phylogenetics and evolution in R language</article-title>
<source>Bioinformatics</source>
<year>2004</year>
<volume>20</volume>
<fpage>289</fpage>
<lpage>290</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btg412</pub-id>
<pub-id pub-id-type="pmid">14734327</pub-id>
</element-citation>
</ref>
<ref id="CR80"><label>80.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schliep</surname>
<given-names>KP</given-names>
</name>
</person-group>
<article-title>phangorn: phylogenetic analysis in R</article-title>
<source>Bioinformatics</source>
<year>2010</year>
<volume>27</volume>
<fpage>592</fpage>
<lpage>593</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btq706</pub-id>
<pub-id pub-id-type="pmid">21169378</pub-id>
</element-citation>
</ref>
<ref id="CR81"><label>81.</label>
<mixed-citation publication-type="other">Wickham, H. <italic>ggplot2: Elegant Graphics for Data Analysis</italic>
(Springer-Verlag, 2016).</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Bois/explor/MelampsoraV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0002089 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0002089 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Bois |area= MelampsoraV2 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.38. | ![]() |