Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

Eucalyptus grandis is a commercially important hardwood species and is known to be susceptible to a number of pests and pathogens. Determining mechanisms of defense is therefore a research priority. The published genome for E. grandis has aided the identification of one important class of resistance (R) genes that incorporate nucleotide binding sites and leucine-rich repeat domains (NBS-LRR). Using an iterative search process we identified NBS-LRR gene models within the E. grandis genome. We characterized the gene models and identified their genomic arrangement. The gene expression patterns were examined in E. grandis clones, challenged with a fungal pathogen (Chrysoporthe austroafricana) and insect pest (Leptocybe invasa). 1215 putative NBS-LRR coding sequences were located which aligned into two large classes, Toll or interleukin-1 receptor (TIR) and coiled-coil (CC) based on NB-ARC domains. NBS-LRR gene-rich regions were identified with 76% organized in clusters of three or more genes. A further 272 putative incomplete resistance genes were also identified. We determined that E. grandis has a higher ratio of TIR to CC classed genes compared to other woody plant species as well as a smaller percentage of single NBS-LRR genes. Transcriptome profiles indicated expression hotspots, within physical clusters, including expression of many incomplete genes. The clustering of putative NBS-LRR genes correlates with differential expression responses in resistant and susceptible plants indicating functional relevance for the physical arrangement of this gene family. This analysis of the repertoire and expression of E. grandis putative NBS-LRR genes provides an important resource for the identification of novel and functional R-genes; a key objective for strategies to enhance resilience.

The fungal plant pathogen Austropuccinia psidii is spreading globally and causing myrtle rust disease symptoms on plants in the family Myrtaceae. A. psidii is dikaryotic, with two nuclei that do not exchange genetic material during the dominant phase of its life-cycle. Phased and scaffolded genome resources for rust fungi are important for understanding heterozygosity, mechanisms of pathogenicity, pathogen population structure and for determining the likelihood of disease spread. We have assembled a chromosome-level phased genome for the pandemic biotype of A. psidii and, for the first time, show that each nucleus contains 18 chromosomes, in line with other distantly related rust fungi. We show synteny between the two haplo-phased genomes and provide a new tool, ChromSyn, that enables efficient comparisons between chromosomes based on conserved genes. Our genome resource includes a fully assembled and circularised mitochondrial sequence for the pandemic biotype.

We assembled the genome of Leptospermum scoparium ‘Crimson Glory’ using a combination of Illumina paired‐end sequencing, high‐throughput chromosome conformation capture (Hi‐C) and high density genetic mapping. As ‘Crimson Glory’ is a variety of mānuka, this is the first genome assembly for a plant species culturally recognised as a treasure (taonga) by the indigenous Māori of Aotearoa New Zealand. The mānuka genome spans a total of 297 Mbp organised in 11 pseudo‐chromosomes that are syntenic with the Eucalyptus genome. A large proportion of the genome assembly corresponds to fungal and bacterial sequences, indicating the presence of an associated microbiome. A total of 31,220 protein‐coding gene models were detected throughout the genome, including genes involved in biosynthesis of biologically active phenylpropanoids, triketones and terpenes, as well as genes involved in biotic resistance. The mānuka genome sequence will help shed new light on the genetic control of unique characters such as nectar and foliage biochemical composition, flowering time and disease resistance.

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state. We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains. The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.