Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility

Fungi Tryptophan Synthases: What Is the Role of the Linker Connecting the α and β Structural Domains in Hemileia vastatrix TRPS? A Molecular Dynamics Investigation

Tryptophan synthase (TRPS) is a complex enzyme responsible for tryptophan biosynthesis. It occurs in bacteria, plants, and fungi as an αββα heterotetramer. Although encoded by independent genes in bacteria and plants, in fungi, TRPS is generated by a single gene that concurrently expresses the α and β entities, which are linked by an elongated peculiar segment. We conducted 1 µs all-atom molecular dynamics simulations on Hemileia vastatrix TRPS to address two questions: (i) the role of the linker segment and (ii) the comparative mode of action. Since there is not an experimental structure, we started our simulations with homology modeling. Based on the results, it seems that TRPS makes use of an already-existing tunnel that can spontaneously move the indole moiety from the α catalytic pocket to the β one. Such behavior was completely disrupted in the simulation without the linker. In light of these results and the αβ dimer's low stability, the full-working TRPS single genes might be the result of a particular evolution. Considering the significant losses that Hemileia vastatrix causes to coffee plantations, our next course of action will be to use the TRPS to look for substances that can block tryptophan production and therefore control the disease.

Genome sequence of Hemileia vastatrix Berk. and Br. (Race I), the causal agent of coffee leaf rust, isolate from Risaralda, Colombia

Coffee leaf rust, caused by the fungus Hemileia vastatrix (Basidiomycota; Pucciniomycota), is a devastating disease spread worldwide. To improve the available genomes, we use PacBio HiFi sequencing enhanced by Dovetail Omni-C chromatin conformation capture to assemble a highly contiguous 747.98 Mb genome of an isolate collected from Coffea arabica.

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.

Gapless Genome Assembly of Puccinia triticina Provides Insights into Chromosome Evolution in Pucciniales

Chromosome evolution drives species evolution, speciation, and adaptive radiation. Accurate genome assembly is crucial to understanding chromosome evolution of species, such as dikaryotic fungi. Rust fungi (Pucciniales) in dikaryons represent the largest group of plant pathogens, but the evolutionary process of adaptive radiation in Pucciniales remains poorly understood. Here, we report a gapless genome for the wheat leaf rust fungus Puccinia triticina determined using PacBio high-fidelity (HiFi) sequencing. This gapless assembly contains two sets of chromosomes, showing that one contig represents one chromosome. Comparisons of homologous chromosomes between the phased haplotypes revealed that highly frequent small-scale sequence divergence shapes haplotypic variation. Genome analyses of Puccinia triticina along with other rusts revealed that recent transposable element bursts and extensive segmental gene duplications synergistically highlight the evolution of chromosome structures. Comparative analysis of chromosomes indicated that frequent chromosomal rearrangements may act as a major contributor to rapid radiation of Pucciniales. This study presents the first gapless, phased assembly for a dikaryotic rust fungus and provides insights into adaptive evolution and species radiation in Pucciniales. IMPORTANCE Rust fungi (Pucciniales) are the largest group of plant pathogens. Adaptive radiation is a predominant feature in Pucciniales evolution. Chromosome evolution plays an important role in adaptive evolution. Accurate chromosome-scale assembly is required to understand the role of chromosome evolution in Pucciniales. We took advantage of HiFi sequencing to construct a gapless, phased genome for Puccinia triticina. Further analyses revealed that the evolution of chromosome structures in rust lineage is shaped by the combination of transposable element bursts and segmental gene duplications. Chromosome comparisons of Puccinia triticina and other rusts suggested that frequent chromosomal arrangements may make remarkable contributions to high species diversity of rust fungi. Our results present the first gapless genome for Pucciniales and shed light on the feature of chromosome evolution in Pucciniales.

Marker-Assisted Recurrent Selection for Pyramiding Leaf Rust and Coffee Berry Disease Resistance Alleles in Coffea arabica L

In this study, marker-assisted recurrent selection was evaluated for pyramiding resistance gene alleles against coffee leaf rust (CLR) and coffee berry diseases (CBD) in Coffea arabica. A total of 144 genotypes corresponding to 12 hybrid populations from crosses between eight parent plants with desired morphological and agronomic traits were evaluated. Molecular data were used for cross-certification, diversity study and resistance allele marker-assisted selection (MAS) against the causal agent of coffee leaf rust (Hemileia vastatrix) and coffee berry disease (Colletotrichum kahawae). In addition, nine morphological and agronomic traits were evaluated to determine the components of variance, select superior hybrids, and estimate genetic gain. From the genotypes evaluated, 134 were confirmed as hybrids. The genetic diversity between and within populations was 75.5% and 24.5%, respectively, and the cluster analysis revealed three primary groups. Pyramiding of CLR and CBD resistance genes was conducted in 11 genotypes using MAS. A selection intensity of 30% resulted in a gain of over 50% compared to the original population. Selected hybrids with increased gain also showed greater genetic divergence in addition to the pyramided resistance alleles. The strategies used were, therefore, efficient to select superior coffee hybrids for recurrent selection programs and could be used as a source of resistance in various crosses.

Whole genome sequencing of Clarireedia aff. paspali reveals potential pathogenesis factors in Clarireedia species, causal agents of dollar spot in turfgrass

Dollar spot is one of the most damaging diseases in turfgrass, reducing its quality and playability. Two species, Clarireedia monteithiana and C. jacksonii (formerly Sclerotinia homoeocarpa) have been reported so far in the United States To study the Clarireedia genome, two isolates H2 and H3, sampled from seashore paspalum in Hawaii in 2019 were sequenced via Illumina paired-end sequencing by synthesis technology and PacBio SMRT sequencing. Both isolates were identified as C. aff. paspali, a novel species in the United States Using short and long reads, C. aff. paspali H3 contained 193 contigs with 48.6 Mbp and presented the most completed assembly and annotation among Clarireedia species. Out of the 13,428 protein models from AUGUSTUS, 349 cytoplasmic effectors and 13 apoplastic effectors were identified by EffectorP. To further decipher Clarireedia pathogenicity, C. aff. paspali genomes (H2 and H3), as well as available C. jacksonii (LWC-10 and HRI11), C. monteithiana (DRR09 and RB-19) genomes were screened for fifty-four pathogenesis determinants, previously identified in S. sclerotiorum. Seventeen orthologs of pathogenicity genes have been identified in Clarireedia species involved in oxalic acid production (pac1, nox1), mitogen-activated protein kinase cascade (pka1, smk3, ste12), appressorium formation (caf1, pks13, ams2, rgb1, rhs1) and glycolytic pathway (gpd). Within these genes, 366 species-specific SNPs were recorded between Clarireedia species; twenty-eight were non-synonymous and non-conservative. The predicted protein structure of six of these genes showed superimposition of the models among Clarireedia spp. The genomic variations revealed here could potentially lead to differences in pathogenesis and other physiological functions among Clarireedia species.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H. vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H. vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens

Biotrophic plant pathogenic fungi are widely distributed and are among the most damaging pathogenic organisms of agriculturally important crops responsible for significant losses in quality and yield. However, the pathogenesis of obligate parasitic pathogenic microorganisms is still under investigation because they cannot reproduce and complete their life cycle on an artificial medium. The successful lifestyle of biotrophic fungal pathogens depends on their ability to secrete effector proteins to manipulate or evade plant defense response. By integrating genomics, transcriptomics, and effectoromics, insights into how the adaptation of biotrophic plant fungal pathogens adapt to their host populations can be gained. Efficient tools to decipher the precise molecular mechanisms of rust–plant interactions, and standardized routines in genomics and functional pipelines have been established and will pave the way for comparative studies. Deciphering fungal pathogenesis not only allows us to better understand how fungal pathogens infect host plants but also provides valuable information for plant diseases control, including new strategies to prevent, delay, or inhibit fungal development. Our review provides a comprehensive overview of the efforts that have been made to decipher the effector proteins of biotrophic fungal pathogens and demonstrates how rapidly research in the field of obligate biotrophy has progressed.

LysM receptors in Coffea arabica: Identification, characterization, and gene expression in response to Hemileia vastatrix

Pathogen‐associated molecular patterns (PAMPs) are recognized by pattern recognition receptors (PRRs) localized on the host plasma membrane. These receptors activate a broad-spectrum and durable defense, which are desired characteristics for disease resistance in plant breeding programs. In this study, candidate sequences for PRRs with lysin motifs (LysM) were investigated in the Coffea arabica genome. For this, approaches based on the principle of sequence similarity, conservation of motifs and domains, phylogenetic analysis, and modulation of gene expression in response to Hemileia vastatrix were used. The candidate sequences for PRRs in C. arabica (Ca1-LYP, Ca2-LYP, Ca1-CERK1, Ca2-CERK1, Ca-LYK4, Ca1-LYK5 and Ca2-LYK5) showed high similarity with the reference PRRs used: Os-CEBiP, At-CERK1, At-LYK4 and At-LYK5. Moreover, the ectodomains of these sequences showed high identity or similarity with the reference sequences, indicating structural and functional conservation. The studied sequences are also phylogenetically related to the reference PRRs described in Arabidopsis, rice, and other plant species. All candidates for receptors had their expression induced after the inoculation with H. vastatrix, since the first time of sampling at 6 hours post‐inoculation (hpi). At 24 hpi, there was a significant increase in expression, for most of the receptors evaluated, and at 48 hpi, a suppression. The results showed that the candidate sequences for PRRs in the C. arabica genome display high homology with fungal PRRs already described in the literature. Besides, they respond to pathogen inoculation and seem to be involved in the perception or signaling of fungal chitin, acting as receptors or co-receptors of this molecule. These findings represent an advance in the understanding of the basal immunity of this species.

Evaluation of the Physicochemical Properties of Chitosans in Inducing the Defense Response of Coffea arabica against the Fungus Hemileia vastatrix

Chitosan is a natural polymer, and its biological properties depend on factors such as the degree of deacetylation and polymerization, viscosity, molecular mass, and dissociation constant. Chitosan has multiple advantages: it is biodegradable, biocompatible, safe, inexpensive, and non-toxic. Due to these characteristics, it has a wide range of applications. In agriculture, one of the most promising properties of chitosan is as an elicitor in plant defense against pathogenic microorganisms. In this work, four kinds of chitosan (practical grade, low molecular weight, medium molecular weight, and high-density commercial food grade) were used in concentrations of 0.01 and 0.05% to evaluate its protective effect against coffee rust. The best treatment was chosen to evaluate the defense response in coffee plants. The results showed a protective effect using practical-grade and commercial food-grade chitosan. In addition, the activity of enzymes with β-1,3 glucanase and peroxidase was induced, and an increase in the amount of phenolic compounds was observed in plants treated with high-molecular-weight chitosan at 0.05%; therefore, chitosan can be considered an effective molecule for controlling coffee rust.

Spatial and Temporal Genetic Diversity and Population Structure of Hemileia vastatrix from Peruvian Coffee Plantations

Population genetic studies of Hemileia vastatrix have been conducted in order to describe the evolutionary dynamics of the pathogen and the disease epidemiology as consequence of changes in disease management and host distribution occurred in Peru after the 2013 epidemic. These analyses were performed by sequencing the internal transcribed spacers of the nuclear ribosomal DNA (rDNA-ITS) of H. vastatrix collected from two coffee growing areas in 2014 and 2018. H. vastatrix population showed high haplotype diversity (Hd = 0.9373 ± 0.0115) with a low nucleotide diversity (π = 0.00322 ± 0.00018). Likewise, AMOVA indicated that fungus population has behaved as a large population without structuring by geographical origin and sampling years (FST = 0.00180, P = 0.20053 and FST = 0.00241, P = 0.19693, respectively). Additionally, the haplotype network based on intraspecific phylogenetic analysis of H. vastatrix using Peruvian and NCBI sequences revealed that Peruvian ancestral haplotypes, which were maintained in time and space, would correspond to the reported sequences of the races II and XXII. This result suggests that no substantial changes have occurred through time in Peruvian Hemileia vastatrix population.

Recent Advances in Molecular Diagnostics of Fungal Plant Pathogens: A Mini Review

Phytopathogenic fungal species can cause enormous losses in quantity and quality of crop yields and this is a major economic issue in the global agricultural sector. Precise and rapid detection and identification of plant infecting fungi are essential to facilitate effective management of disease. DNA-based methods have become popular methods for accurate plant disease diagnostics. Recent developments in standard and variant polymerase chain reaction (PCR) assays including nested, multiplex, quantitative, bio and magnetic-capture hybridization PCR techniques, post and isothermal amplification methods, DNA and RNA based probe development, and next-generation sequencing provide novel tools in molecular diagnostics in fungal detection and differentiation fields. These molecular based detection techniques are effective in detecting symptomatic and asymptomatic diseases of both culturable and unculturable fungal pathogens in sole and co-infections. Even though the molecular diagnostic approaches have expanded substantially in the recent past, there is a long way to go in the development and application of molecular diagnostics in plant diseases. Molecular techniques used in plant disease diagnostics need to be more reliable, faster, and easier than conventional methods. Now the challenges are with scientists to develop practical techniques to be used for molecular diagnostics of plant diseases. Recent advancement in the improvement and application of molecular methods for diagnosing the widespread and emerging plant pathogenic fungi are discussed in this review.

Effector Biology of Biotrophic Plant Fungal Pathogens: Current Advances and Future Prospects

The interaction of fungal pathogens with their host requires a novel invading mechanism and the presence of various virulence-associated components responsible for promoting the infection. The small secretory proteins, explicitly known as effector proteins, are one of the prime mechanisms of host manipulation utilized by the pathogen to disarm the host. Several effector proteins are known to translocate from fungus to the plant cell for host manipulation. Many fungal effectors have been identified using genomic, transcriptomic, and bioinformatics approaches. Most of the effector proteins are devoid of any conserved signatures, and their prediction based on sequence homology is very challenging, therefore by combining the sequence consensus based upon machine learning features, multiple tools have also been developed for predicting apoplastic and cytoplasmic effectors. Various post-genomics approaches like transcriptomics of virulent isolates have also been utilized for identifying active consortia of effectors. Significant progress has been made in understanding biotrophic effectors; however, most of it is underway due to their complex interaction with host and complicated recognition and signaling networks. This review discusses advances, and challenges in effector identification and highlighted various features of the potential effector proteins and approaches for understanding their genetics and strategies for regulation.

Evolution of virulence in rust fungi - multiple solutions to one problem

Rust fungi are major pathogens that negatively affect crops and ecosystems. Recent rust disease epidemics driven by the emergence of strains with novel virulence profiles demand a better understanding of the evolutionary mechanisms of these organisms. Here, we review research advances in genome-scale analysis coupled with functional validation of effector candidate genes that have been instrumental to elucidate processes that contribute to changes in virulence phenotypes. We highlight how haplotype-phased genome references have paved the road to link these processes to the reproductive phases of rust fungi and have provided evidence for somatic exchange between strains as an important mechanism for generating diversity in asexual populations. With increasing data availability, we envision the future development of molecular virulence diagnostic tools.

The Conifer Root and Stem Rot Pathogen (Heterobasidion parviporum): Effectome Analysis and Roles in Interspecific Fungal Interactions

Heterobasidion parviporum Niemelä & Korhonen is an economically important basidiomycete, causing root and stem rot disease of Norway spruce (Picea abies (L.) Karst) in Northern Europe. The H. parviporum genome encodes numerous small secreted proteins, which might be of importance for interacting with mycorrhiza symbionts, endophytes, and other saprotrophs. We hypothesized that small secreted proteins from H. parviporum (HpSSPs) are involved in interspecific fungal interaction. To identify HpSSP-coding genes potentially involved, we screened the H. parviporum effectome and compared their transcriptomic profiles during fungal development and in planta tree infection. We further conducted phylogenetic analysis, and identified a subset of hypothetical proteins with nonpredicted domain or unknown function as HpSSPs candidates for further characterization. The HpSSPs candidates were selected based on high-quality sequence, cysteine residue frequency, protein size, and in planta expression. We subsequently explored their roles during in vitro interaction in paired cultures of H. parviporum with ectomycorrhizal Cortinarius gentilis, endophytic Phialocephala sphaeroides, saprotrophs (Mycena sp., Phlebiopsis gigantea, and Phanerochaete chrysosporium), respectively. The transcriptomic profile revealed that a large proportion of effector candidates was either barely expressed or highly expressed under all growth conditions. In vitro dual-culture test showed that P. sphaeroides and C. gentilis were overgrown by H. parviporum. The barrage zone formation or no physical contact observed in paired cultures with the saprotrophs suggest they had either combative interaction or antibiosis effect with H. parviporum. Several HpSSPs individuals were up- or downregulated during the nonself interactions. The results of HpSSPs gene expression patterns provide additional insights into the diverse roles of SSPs in tree infection and interspecific fungal interactions.