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, Puccinia graminis f. sp. tritici (Pgt), 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 Pgt, 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.

The genetic distance between hybridizing parents affects heterosis; however, the mechanisms for this remain unclear. Here we report that this genetic distance correlates with natural variation and epigenetic regulation of circadian clock-mediated stress responses. In intraspecific hybrids of Arabidopsis thaliana, genome-wide expression of many biotic and abiotic stress-responsive genes is diurnally repressed and this correlates with biomass heterosis and biomass quantitative trait loci. Expression differences of selected stress-responsive genes among diverse ecotypes are predictive of heterosis in their hybrids. Stress-responsive genes are repressed in the hybrids under normal conditions but are induced to mid-parent or higher levels under stress at certain times of the day, potentially balancing the tradeoff between stress responses and growth. Consistent with this hypothesis, repression of two candidate stress-responsive genes increases growth vigour. Our findings may therefore provide new criteria for effectively selecting parents to produce high-or lowyield hybrids.

A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales. In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies.

Both ploidy and hybridity affect cell size and growth vigor in plants and animals, but the relative effects of genome dosage and hybridization on biomass, fitness, and gene expression changes have not been systematically examined. Here we performed the first comparative analysis of seed, cell, and flower sizes, starch and chlorophyll content, biomass, and gene expression changes in diploid, triploid, and tetraploid hybrids and their respective parents in three Arabidopsis thaliana ecotypes: Columbia, C24, and Landsberg erecta (Ler). Ploidy affects many morphological and fitness traits, including stomatal size, flower size, and seed weight, whereas hybridization between the ecotypes leads to altered expression of central circadian clock genes and increased starch and chlorophyll content, biomass, and seed weight. However, varying ploidy levels has subtle effects on biomass, circadian clock gene expression, and chlorophyll and starch content. Interestingly, biomass, starch content, and seed weight are significantly different between the reciprocal hybrids at all ploidy levels tested, with the lowest and highest levels found in the reciprocal triploid hybrids, suggesting parent-of-origin effects on biomass, starch content, and seed weight. These findings provide new insights into molecular events of polyploidy and heterosis, as well as complex agronomic traits that are important to biomass and seed production in hybrid and polyploid crops.

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenae.