Helena Gil Azinheira

Coffee resistance to the main diseases: leaf rust and coffee berry disease

Considerable success has been obtained in the use of classical breeding to control economically important plant diseases, such as the coffee leaf rust and the coffee berry disease (CBD). There is a strong consensus that growing genetically resistant varieties is the most appropriate cost effective means of managing plant diseases and is one of the key components of crop improvement. It has also been recognized that a better knowledge of both, the pathogens and the plant defence mechanisms will allow the development of novel approaches to enhance the durability of resistance. After a brief description of concepts in the field of plant disease resistance, we attempt to give a view of the research progress on coffee leaf rust and CBD concerned with the pathogens infection and variability, coffee breeding for resistance and coffee resistance mechanisms.

454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant-rust interaction

Coffee (Coffea arabica L.), one of the key export and cash crops in tropical and subtropical countries, suffers severe losses from the rust fungus Hemileia vastatrix. The transcriptome of H. vastatrix was analysed during a compatible interaction with coffee to obtain an exhaustive repertoire of the genes expressed during infection and to identify potential effector genes. Large-scale sequencing (454-GS-FLEX Titanium) of mixed coffee and rust cDNAs obtained from 21-day rust-infected leaves generated 352 146 sequences which assembled into 22 774 contigs. In the absence of any reference genomic sequences for Coffea or Hemileia, specific trinucleotide frequencies within expressed sequence tags (ESTs) and blast homology against a set of dicots and basidiomycete genomes were used to distinguish pathogen from plant sequences. About 30% (6763) of the contigs were assigned to H. vastatrix and 61% (13 951) to C. arabica. The majority (60%) of the rust sequences did not show homology to any genomic database, indicating that they were potential novel fungal genes. In silico analyses of the 6763 H. vastatrix contigs predicted 382 secreted proteins and identified homologues of the flax rust haustorially expressed secreted proteins (HESPs) and bean rust transferred protein 1 (RTP1). These rust candidate effectors showed conserved amino-acid domains and conserved patterns of cysteine positions suggestive of conserved functions during infection of host plants. Quantitative reverse transcription-polymerase chain reaction profiling of selected rust genes revealed dynamic expression patterns during the time course of infection of coffee leaves. This study provides the first valuable genomic resource for the agriculturally important plant pathogen H. vastatrix and the first comprehensive C. arabica EST dataset.

Genome size analyses of Pucciniales reveal the largest fungal genomes.

Rust fungi (Basidiomycota, Pucciniales) are biotrophic plant pathogens which exhibit diverse complexities in their life cycles and host ranges. The completion of genome sequencing of a few rust fungi has revealed the occurrence of large genomes. Sequencing efforts for other rust fungi have been hampered by uncertainty concerning their genome sizes. Flow cytometry was recently applied to estimate the genome size of a few rust fungi, and confirmed the occurrence of large genomes in this order (averaging 225.3 Mbp, while the average for Basidiomycota was 49.9 Mbp and was 37.7 Mbp for all fungi). In this work, we have used an innovative and simple approach to simultaneously isolate nuclei from the rust and its host plant in order to estimate the genome size of 30 rust species by flow cytometry. Genome sizes varied over 10-fold, from 70 to 893 Mbp, with an average genome size value of 380.2 Mbp. Compared to the genome sizes of over 1800 fungi, Gymnosporangium confusum possesses the largest fungal genome ever reported (893.2 Mbp). Moreover, even the smallest rust genome determined in this study is larger than the vast majority of fungal genomes (94%). The average genome size of the Pucciniales is now of 305.5 Mbp, while the average Basidiomycota genome size has shifted to 70.4 Mbp and the average for all fungi reached 44.2 Mbp. Despite the fact that no correlation could be drawn between the genome sizes, the phylogenomics or the life cycle of rust fungi, it is interesting to note that rusts with Fabaceae hosts present genomes clearly larger than those with Poaceae hosts. Although this study comprises only a small fraction of the more than 7000 rust species described, it seems already evident that the Pucciniales represent a group where genome size expansion could be a common characteristic. This is in sharp contrast to sister taxa, placing this order in a relevant position in fungal genomics research.

Cellular and molecular analyses of coffee resistance to Hemileia vastatrix and nonhost resistance to Uromyces vignae in the resistance-donor genotype HDT832/2

In Arabica coffee breeding, some of the most used sources of resistance to leaf rust (Hemileia vastatrix) are natural Coffea arabica x canephora hybrids (“Híbrido de Timor”). To decipher the cellular and molecular nature of that resistance, leaves of genotype HDT832/2, were challenged with H. vastatrix race II, and monitored using light microscopy and RT-qPCR expression analysis of genes involved in plant immunity mechanisms (receptor-like kinase, WRKY transcription factor 1, phenylalanine ammonia-lyase, chalcone synthase, 13-lipoxygenase, glycosyltransferase, pathogenesis related PR1b and PR10). These were compared to the nonhost resistance responses of HDT832/2 to the infection by the cowpea rust fungus (Uromyces vignae). H. vastatrix ceased growth more frequently after stomata penetration, forming few haustoria, inducing a hypersensitive-like response, phenol accumulation and haustorium encasement with callose. U. vignae could enter stomata but failed to form haustoria, while inducing hypersensitive-like responses and phenol accumulation. In host and nonhost interactions, activation of genes involved in signalling coincided with the differentiation of appressoria, and cellular responses (hypersensitive-like responses and accumulation of phenolic compounds) were recorded from the full appressorium or penetration hypha stages onwards. Similarly, a gene related to the JA pathway was first activated at the penetration hypha stage for both interactions, while genes related to the SA pathway were only activated in the host interaction, the latter being the single clear difference between host and nonhost interactions. The cellular and molecular resistance responses of HDT832/2 to these rust fungi suggest that common immunity components are shared between host and nonhost resistance, which may explain the longer durability of this resistance.

Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection

Hemileia vastatrix is the causal agent of coffee leaf rust, the most important disease of coffee Arabica. In this work, a 454-pyrosequencing transcriptome analysis of H. vastatrix germinating urediniospores (gU) and appressoria (Ap) was performed and compared to previously published in planta haustoria-rich (H) data. A total of 9234 transcripts were identified and annotated. Ca. 50% of these transcripts showed no significant homology to international databases. Only 784 sequences were shared by the three conditions, and 75% were exclusive of either gU (2146), Ap (1479) or H (3270). Relative transcript abundance and RT-qPCR analyses for a selection of genes indicated a particularly active metabolism, translational activity and production of new structures in the appressoria and intense signaling, transport, secretory activity and cellular multiplication in the germinating urediniospores, suggesting the onset of a plant-fungus dialogue as early as at the germ tube stage. Gene expression related to the production of carbohydrate-active enzymes and accumulation of glycerol in germinating urediniospores and appressoria suggests that combined lytic and physical mechanisms are involved in appressoria-mediated penetration. Besides contributing to the characterization of molecular processes leading to appressoria-mediated infection by rust fungi, these results point toward the identification of new H. vastatrix candidate virulence factors, with 516 genes predicted to encode secreted proteins.

Expression profiling of genes involved in the biotrophic colonisation of Coffea arabica leaves by Hemileia vastatrix

Coffee Leaf Rust, caused by Hemileia vastatrix, is the most important disease of Arabica coffee (Coffea arabica), which prompts studies aimed at understanding the genetic basis of this pathogen as well as its complex interaction with the host. In this work, 11 genes, putatively involved in signalling, establishment and maintenance of biotrophy (transport and metabolism), were characterised, and their expression profiles during host infection were assessed by RT-qPCR in three compatible coffee-rust interactions comprising two different rust races. The profiles of two chitin deacetylases (CD) and a heterotrimeric G-protein α subunit transcripts suggest that these enzymes are involved in host-pathogen recognition and establishment of biotrophy at early stages of infection, and the late expression of the CD1 gene was also recorded. Different expression profiles were observed for a MAP kinase gene between the two rust races, suggesting that this gene may be involved in the differentiation of infection structures in a race-specific pattern. Two amino acid transporters, an invertase, a hexose transporter and a mannitol dehydrogenase presented expression profiles similar to those reported in other rust fungi, indicating a fairly conserved genetic programme related to host infection in rust fungi. The strong upregulation of a Uromyces fabae rust transferred protein 1 orthologous gene was observed in H. vastatrix in planta structures, suggesting that this gene may also play a role during the establishment and the maintenance of biotrophy in coffee leaves. Overall, our results provide valuable insights to the current understanding of the biotrophic interaction between H. vastatrix – C. arabica at the molecular level and will contribute to a reasoned and sustainable use of resistant genotypes.

Genomic Patterns of Positive Selection at the Origin of Rust Fungi

Understanding the origin and evolution of pathogenicity and biotrophic life-style of rust fungi has remained a conundrum for decades. Research on the molecular mechanisms responsible for rust fungi evolution has been hampered by their biotrophic life-style until the sequencing of some rust fungi genomes. With the availability of multiple whole genomes and EST data for this group, it is now possible to employ genome-wide surveys and investigate how natural selection shaped their evolution. In this work, we employed a phylogenomics approach to search for positive selection and genes undergoing accelerated evolution at the origin of rust fungi on an assembly of single copy genes conserved across a broad range of basidiomycetes. Up to 985 genes were screened for positive selection on the phylogenetic branch leading to rusts, revealing a pervasive signal of positive selection throughout the data set with the proportion of positively selected genes ranging between 19.6–33.3%. Additionally, 30 genes were found to be under accelerated evolution at the origin of rust fungi, probably due to a mixture of positive selection and relaxation of purifying selection. Functional annotation of the positively selected genes revealed an enrichment in genes involved in the biosynthesis of secondary metabolites and several metabolism and transporter classes.

Comparative Validation of Conventional and RNA-Seq Data-Derived Reference Genes for qPCR Expression Studies of Colletotrichum kahawae

Colletotrichum kahawae is an emergent fungal pathogen causing severe epidemics of Coffee Berry Disease on Arabica coffee crops in Africa. Currently, the molecular mechanisms underlying the Coffea arabica—C. kahawae interaction are still poorly understood, as well as the differences in pathogen aggressiveness, which makes the development of functional studies for this pathosystem a crucial step. Quantitative real time PCR (qPCR) has been one of the most promising approaches to perform gene expression analyses. However, proper data normalization with suitable reference genes is an absolute requirement. In this study, a set of 8 candidate reference genes were selected based on two different approaches (literature and Illumina RNA-seq datasets) to assess the best normalization factor for qPCR expression analysis of C. kahawae samples. The gene expression stability of candidate reference genes was evaluated for four isolates of C. kahawae bearing different aggressiveness patterns (Ang29, Ang67, Zim12 and Que2), at different stages of fungal development and key time points of the plant-fungus interaction process. Gene expression stability was assessed using the pairwise method incorporated in geNorm and the model-based method used by NormFinder software. For C. arabica—C. kahawae interaction samples, the best normalization factor included the combination of PP1, Act and ck34620 genes, while for C. kahawae samples the combination of PP1, Act and ck20430 revealed to be the most appropriate choice. These results suggest that RNA-seq analyses can provide alternative sources of reference genes in addition to classical reference genes. The analysis of expression profiles of bifunctional catalase-peroxidase (cat2) and trihydroxynaphthalene reductase (thr1) genes further enabled the validation of the selected reference genes. This study provides, for the first time, the tools required to conduct accurate qPCR studies in C. kahawae considering its aggressiveness pattern, developmental stage and host interaction.

A method for obtaining RNA from Hemileia vastatrix appressoria produced in planta, suitable for transcriptomic analyses.

Appressoria are the first infection structures developed by rust fungi and require specific topographic signals from the host for their differentiation. The ease in obtaining appressoria in vitro for these biotrophic fungi led to studies concerning gene expression and gene discovery at appressorial level, avoiding the need to distinguish plant and fungal transcripts. However, in some pathosystems, it was observed that gene expression in appressoria seems to be influenced by host-derived signals, suggesting that transcriptomic analyses performed from in planta differentiated appressoria would be potentially more informative than those from in vitro differentiated appressoria. Nevertheless analysing appressorial RNA obtained from in planta samples is often hampered by an excessive dilution of fungal RNA within plant RNA, besides uncertainty regarding the fungal or plant origin of RNA from highly conserved genes. To circumvent these difficulties, we have recovered Hemileia vastatrix appressoria from Arabica coffee leaf surface using a film of nitrocellulose dissolved in butyl and ethyl acetates (nail polish), and extracted fungal RNA from the polish peel. RNA thus obtained is of good quality and usable for cDNA synthesis and transcriptomic (quantitative PCR) studies. This method could provide the means to investigate specific host-induced appressoria-related fungal pathogenicity factors.

A first insight into the involvement of phytohormones pathways in coffee resistance and susceptibility to Colletotrichum kahawae

Understanding the molecular mechanisms underlying coffee-pathogen interactions are of key importance to aid disease resistance breeding efforts. In this work the expression of genes involved in salicylic acid (SA), jasmonic acid (JA) and ethylene (ET) pathways were studied in hypocotyls of two coffee varieties challenged with the hemibiotrophic fungus Colletotrichum kahawae, the causal agent of Coffee Berry Disease. Based on a cytological analysis, key time-points of the infection process were selected and qPCR was used to evaluate the expression of phytohormones biosynthesis, reception and responsive-related genes. The resistance to C. kahawae was characterized by restricted fungal growth associated with early accumulation of phenolic compounds in the cell walls and cytoplasmic contents, and deployment of hypersensitive reaction. Similar responses were detected in the susceptible variety, but in a significantly lower percentage of infection sites and with no apparent effect on disease development. Gene expression analysis suggests a more relevant involvement of JA and ET phytohormones than SA in this pathosystem. An earlier and stronger activation of the JA pathway observed in the resistant variety, when compared with the susceptible one, seems to be responsible for the successful activation of defense responses and inhibition of fungal growth. For the ET pathway, the down or non-regulation of ET receptors in the resistant variety, together with a moderate expression of the responsive-related gene ERF1, indicates that this phytohormone may be related with other functions besides the resistance response. However, in the susceptible variety, the stronger activation of ERF1 gene at the beginning of the necrotrophic phase, suggests the involvement of ET in tissue senescence. As far as we know, this is the first attempt to unveil the role of phytohormones in coffee-C. kahawae interactions, thus contributing to deepen our understanding on the complex mechanisms of plant signaling and defense.