Andreia Loureiro

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.

Host‐jump drives rapid and recent ecological speciation of the emergent fungal pathogen Colletotrichum kahawae

Ecological speciation through host‐shift has been proposed as a major route for the appearance of novel fungal pathogens. The growing awareness of their negative impact on global economies and public health created an enormous interest in identifying the factors that are most likely to promote their emergence in nature. In this work, a combination of pathological, molecular and geographical data was used to investigate the recent emergence of the fungus Colletotrichum kahawae. C. kahawae emerged as a specialist pathogen causing coffee berry disease in Coffea arabica, owing to its unparalleled adaptation of infecting green coffee berries. Contrary to current hypotheses, our results suggest that a recent host‐jump underlay the speciation of C. kahawae from a generalist group of fungi seemingly harmless to coffee berries. We posit that immigrant inviability and a predominantly asexual behaviour could have been instrumental in driving speciation by creating pleiotropic interactions between local adaptation and reproductive patterns. Moreover, we estimate that C. kahawae began its diversification at <2200 bp leaving a very short time frame since the divergence from its sibling lineage (c. 5600 bp), during which a severe drop in C. kahawae’s effective population size occurred. This further supports a scenario of recent introduction and subsequent adaptation to C. arabica. Phylogeographical data revealed low levels of genetic polymorphism but provided the first geographically consistent population structure of C. kahawae, inferring the Angolan population as the most ancestral and the East African populations as the most recently derived. Altogether, these results highlight the significant role of host specialization and asexuality in the emergence of fungal pathogens through ecological speciation.

Validation of reference genes for normalization of qPCR gene expression data from Coffea spp. hypocotyls inoculated with Colletotrichum kahawae

BackgroundCoffee production in Africa represents a significant share of the total export revenues and influences the lives of millions of people, yet severe socio-economic repercussions are annually felt in result of the overall losses caused by the coffee berry disease (CBD). This quarantine disease is caused by the fungus Colletotrichum kahawae Waller and Bridge, which remains one of the most devastating threats to Coffea arabica production in Africa at high altitude, and its dispersal to Latin America and Asia represents a serious concern. Understanding the molecular genetic basis of coffee resistance to this disease is of high priority to support breeding strategies. Selection and validation of suitable reference genes presenting stable expression in the system studied is the first step to engage studies of gene expression profiling.ResultsIn this study, a set of ten genes (S24, 14-3-3, RPL7, GAPDH, UBQ9, VATP16, SAND, UQCC, IDE and β-Tub9) was evaluated to identify reference genes during the first hours of interaction (12, 48 and 72 hpi) between resistant and susceptible coffee genotypes and C. kahawae. Three analyses were done for the selection of these genes considering the entire dataset and the two genotypes (resistant and susceptible), separately. The three statistical methods applied GeNorm, NormFinder, and BestKeeper, allowed identifying IDE as one of the most stable genes for all datasets analysed, and in contrast GADPH and UBQ9 as the least stable ones. In addition, the expression of two defense-related transcripts, encoding for a receptor like kinase and a pathogenesis related protein 10, were used to validate the reference genes selected.ConclusionTaken together, our results provide guidelines for reference gene(s) selection towards a more accurate and widespread use of qPCR to study the interaction between Coffea spp. and C. kahawae.

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.

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.

Legitimacy and Implications of Reducing Colletotrichum kahawae to Subspecies in Plant Pathology

Colletotrichum kahawae Waller and Bridge is a highly aggressive and specialized fungal pathogen of coffee, causing the devastating Coffee Berry Disease (CBD), particularly at high altitudes. The disease arises from the unique ability of the pathogen to infect green developing coffee berries. This pathogen is currently confined to the African continent in all countries that grow Arabica coffee (Coffea arabica L.), leading to up to 80% yield losses, if no control measures are applied (Silva et al., 2006; Vossen and Walyaro, 2009; Hindorf and Omondi, 2011). For such huge economic impact, it is ranked as a quarantine pathogen and even as a biological weapon (Australia Group, 2014). Consequently, the pathogens potential dispersal to other Arabica coffee cultivation regions is greatly feared, particularly to those at high altitude also found in Latin America and Asia. Recently, this recognized species was brought down to a subspecific level (C. kahawae subsp. kahawae) based on molecular data (Weir et al., 2012), clustering together with a generalist and cosmopolitan group of Colletotrichum isolates unable to cause CBD (C. kahawae subsp. ciggaro). Since then a growing number of studies have reported the identification of C. kahawae in various hosts and regions of the world (Liu et al., 2013; Afanador-Kafuri et al., 2014; Mosca et al., 2014; Schena et al., 2014; Ismail et al., 2015; Garibaldi et al., 2016a,b; Perrone et al., 2016). Although these reports are referring to C. kahawae subsp. ciggaro, some of them could not distinguish the pathogen at the subspecific level, and this is leading to a wave of confusion of whether the long accepted species C. kahawae, the CBD pathogen, has escaped from Africa and extended its host range. Given the extreme impact that this situation may trigger and the subsequent biosecurity implications, there is a practical need to completely distinguish these pathogens taxonomically as to avoid the risk of misidentification, and caution should be taken on assigning/reassigning taxonomic ranking and nomenclature. Here we consider the evidences sustaining and contradicting the classification proposed by Weir et al. (2012), and discuss the risks and practical implications of changing the CBD pathogens species status in a plant pathology context.

Aggressiveness profiling of the coffee pathogen Colletotrichum kahawae

A. Vieira, I. Diniz* , A. Loureiro, A. P. Pereira, M. C. Silva, V. V arzea and D. Batista CIFC – Centro de Investigac ~ ao das Ferrugens do Cafeeiro, Instituto Superior de Agronomia, Universidade de Lisboa, 2784-505 Oeiras; CoBiG – Computational Biology and Population Genomics Group, cE3c – Centre for Ecology, Evolution and Environmental Changes, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa; and LEAF – Linking Landscape, Environment, Agriculture and Food, Instituto Superior de Agronomia, Universidade de Lisboa, 1349-017 Lisboa, Portugal