Folding Features and Dynamics of 3D Genome Architecture in Plant Fungal Pathogens

Topologically associating domains and the evolution of three-dimensional genome architecture in rice

We examined the nature and evolution of three-dimensional (3D) genome conformation, including topologically associating domains (TADs), in five genomes within the genus Oryza. These included three varieties from subspecies within domesticated Asian rice O. sativa as well as their closely related wild relatives O. rufipogon and O. meridionalis. We used the high-resolution chromosome conformation capture technique Micro-C, which we modified for use in rice. Our analysis of rice TADs shows that TAD boundaries have high transcriptional activity, low methylation levels, low transposable element (TE) content, and increased gene density. We also find a significant correlation of expression levels for genes within TADs, suggesting that they do function as genomic domains with shared regulatory features. Our findings indicate that animal and plant TADs may share more commonalities than were initially thought, as evidenced by similar genetic and epigenetic signatures associated with TADs and boundaries. To examine 3D genome divergence, we employed a computer vision-based algorithm for the comparison of chromatin contact maps and complemented this analysis by assessing the evolutionary conservation of individual TADs and their boundaries. We conclude that overall chromatin organization is conserved in rice, and 3D structural divergence correlates with evolutionary distance between genomes. We also note that individual TADs are not well conserved, even at short evolutionary timescales.

Chromosome-level genome and annotation of the tetraploid Rhodiola kirilowii

Rhodiola kirilowii, a perennial medical herb native to China, is highly valued for its detoxification and anti-swelling properties, as well as its role as an adaptogen, making it an intriguing subject for understanding its medicinal potential and molecular biochemistry. In this study, we generated a high-quality chromosome-level reference genome of R. kirilowii achieved through a combination of Illumina short-read sequencing, PacBio long-read sequencing, and Hi-C sequencing techniques. The final assembly spans 1.92 Gb, including 40 homoeologous chromosomes and one sex chromosome, with a scaffold NG50 of 46.03 Mb, and a BUSCO completeness of 98.9%. Additionally, we annotated a total of 1.23 Gb of repetitive sequences, encompassing 63.88% of the entire genome, and identified 122,035 protein-coding genes. Each sub-genome achieved similar completeness and continuity. This high-quality reference genome provides critical insights into the genetic underpinnings of R. kirilowii's pharmacological properties, facilitating comparative genomics and the enhancement of its medicinal applications.

Genomics Research on the Road of Studying Biology and Virulence of Cereal Rust Fungi

Rust fungi are highly destructive pathogens that pose a significant threat to crop production worldwide, especially cereals. Obligate biotrophy and, in many cases, complex life cycles make rust fungi particularly challenging to study. However, recent rapid advances in sequencing technologies and genomic analysis tools have revolutionised rust fungal research. It is anticipated that the increasing availability and ongoing substantial improvements in genome assemblies will propel the field of rust biology into the post-genomic era, instigating a cascade of research endeavours encompassing multi-omics and gene discoveries. This is especially the case for many cereal rust pathogens, for which continental-scale studies of virulence have been conducted over many years and historical collections of viable isolates have been sequenced and assembled. Genomic analysis plays a crucial role in uncovering the underlying causes of the high variability of virulence and the complexity of population dynamics in rust fungi. Here, we provide an overview of progress in rust genomics, discuss the strategies employed in genomic analysis, and elucidate the strides that will drive cereal rust biology into the post-genomic era.

The Three-Dimensional Structure of the Genome of the Dark Septate Endophyte Exophiala tremulae and Its Symbiosis Effect on Alpine Meadow Plant Growth

The establishment of artificial grassland is a good pathway for resolving serious social and economic problems in the Qinghai-Tibet Plateau. Some beneficial indigenous microbes may be used to improve productivity in artificial grassland. The genome of the indigenous dark septate fungus, Exophiala tremulae CICC2537, was sequenced and assembled at the chromosome level using the PacBio sequencing platform, with the assistance of the Hi-C technique for scaffolding, and its 3D genome structures were investigated. The genome size of E. tremulae is 51.903848 Mb, and it contains eight chromosomes. A total of 12,277 protein-coding genes were predicted, and 11,932 genes (97.19%) were annotated. As for the distribution of exon and intron number and the distribution of gene GC and CDS GC, E. tremulae showed similar distribution patterns to the other investigated members of the genus Exophiala. The analysis of carbohydrate-active enzymes showed that E. tremulae possesses the greatest number of enzymes with auxiliary activities and the lowest number of enzymes with carbohydrate-binding modules among the investigated fungi. The total number of candidate effector proteins was 3337, out of which cytoplasmic and apoplastic effector proteins made up 3100 and 163, respectively. The whole genome of E. tremulae contained 40 compartment As and 76 compartment Bs, and there was no significant difference in GC content in its compartment As and Bs. The whole genome of E. tremulae was predicted to contain 155 topologically associating domains (TADs), and their average length was 250,000 bp, but there were no significant differences in the numbers of genes and the GC content per bin localized within the boundaries and interiors of TADs. Comparative genome analysis showed that E. tremulae diverged from Exophiala mesophila about 34.1 (30.0-39.1) Myr ago, and from Exophiala calicioides about 85.6 (76.1-90.6) Myr ago. Compared with all the investigated fungi, the numbers of contraction and expansion gene families in the E. tremulae genome were 13 and 89, respectively, and the numbers of contraction and expansion genes were 14 and 670, respectively. Our work provides a basis for the use of the dark septate fungus in alpine artificial grassland and further research into its symbiosis mechanisms, which may improve the growth of plant species used in the Qinghai-Tibet Plateau.

Complete telomere-to-telomere genomes uncover virulence evolution conferred by chromosome fusion in oomycete plant pathogens

Variations in chromosome number are occasionally observed among oomycetes, a group that includes many plant pathogens, but the emergence of such variations and their effects on genome and virulence evolution remain ambiguous. We generated complete telomere-to-telomere genome assemblies for Phytophthora sojae, Globisporangium ultimum, Pythium oligandrum, and G. spinosum. Reconstructing the karyotype of the most recent common ancestor in Peronosporales revealed that frequent chromosome fusion and fission drove changes in chromosome number. Centromeres enriched with Copia-like transposons may contribute to chromosome fusion and fission events. Chromosome fusion facilitated the emergence of pathogenicity genes and their adaptive evolution. Effectors tended to duplicate in the sub-telomere regions of fused chromosomes, which exhibited evolutionary features distinct to the non-fused chromosomes. By integrating ancestral genomic dynamics and structural predictions, we have identified secreted Ankyrin repeat-containing proteins (ANKs) as a novel class of effectors in P. sojae. Phylogenetic analysis and experiments further revealed that ANK is a specifically expanded effector family in oomycetes. These results revealed chromosome dynamics in oomycete plant pathogens, and provided novel insights into karyotype and effector evolution.

The jet-like chromatin structure defines active secondary metabolism in fungi

Eukaryotic genomes are spatially organized within the nucleus in a nonrandom manner. However, fungal genome arrangement and its function in development and adaptation remain largely unexplored. Here, we show that the high-order chromosome structure of Fusarium graminearum is sculpted by both H3K27me3 modification and ancient genome rearrangements. Active secondary metabolic gene clusters form a structure resembling chromatin jets. We demonstrate that these jet-like domains, which can propagate symmetrically for 54 kb, are prevalent in the genome and correlate with active gene transcription and histone acetylation. Deletion of GCN5, which encodes a core and functionally conserved histone acetyltransferase, blocks the formation of the domains. Insertion of an exogenous gene within the jet-like domain significantly augments its transcription. These findings uncover an interesting link between alterations in chromatin structure and the activation of fungal secondary metabolism, which could be a general mechanism for fungi to rapidly respond to environmental cues, and highlight the utility of leveraging three-dimensional genome organization in improving gene transcription in eukaryotes.

A fully haplotype-resolved and nearly gap-free genome assembly of wheat stripe rust fungus

Stripe rust fungus Puccinia striiformis f. sp. tritici (Pst) is a destructive pathogen of wheat worldwide. Pst has a macrocyclic-heteroecious lifecycle, in which one-celled urediniospores are dikaryotic, each nucleus containing one haploid genome. We successfully generated the first fully haplotype-resolved and nearly gap-free chromosome-scale genome assembly of Pst by combining PacBio HiFi sequencing and trio-binning strategy. The genome size of the two haploid assemblies was 75.59 Mb and 75.91 Mb with contig N50 of 4.17 Mb and 4.60 Mb, and both had 18 pseudochromosomes. The high consensus quality values of 55.57 and 59.02 for both haplotypes confirmed the correctness of the assembly. Of the total 18 chromosomes, 15 and 16 were gapless while there were only five and two gaps for the remaining chromosomes of the two haplotypes, respectively. In total, 15,046 and 15,050 protein-coding genes were predicted for the two haplotypes, and the complete BUSCO scores achieved 97.7% and 97.9%, respectively. The genome will lay the foundation for further research on genetic variations and the evolution of rust fungi.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.