The chloroplast protein HCF164 is predicted to be associated with Coffea SH9 resistance factor against Hemileia vastatrix.

To explore the connection between chloroplast and coffee resistance factors, designated as SH1 to SH9, whole genomic DNA of 42 coffee genotypes was sequenced, and entire chloroplast genomes were de novo assembled. The chloroplast phylogenetic haplotype network clustered individuals per species instead of SH factors. However, for the first time, it allowed the molecular validation of Coffea arabica as the maternal parent of the spontaneous hybrid "Híbrido de Timor". Individual reads were also aligned on the C. arabica reference genome to relate SH factors with chloroplast metabolism, and an in-silico analysis of selected nuclear-encoded chloroplast proteins (132 proteins) was performed. The nuclear-encoded thioredoxin-like membrane protein HCF164 enabled the discrimination of individuals with and without the SH9 factor, due to specific DNA variants linked to chromosome 7c (from C. canephora-derived sub-genome). The absence of both the thioredoxin domain and redox-active disulphide center in the HCF164 protein, observed in SH9 individuals, raises the possibility of potential implications on redox regulation. For the first time, the identification of specific DNA variants of chloroplast proteins allows discriminating individuals according to the SH profile. This study introduces an unexplored strategy for identifying protein/genes associated with SH factors and candidate targets of H. vastatrix effectors, thereby creating new perspectives for coffee breeding programs.

First Report of the Physiological Race (XXIV) of Hemileia vastatrix (Coffee Leaf Rust) in Hawaii.

Hawaii's coffee industry, produced commercially on six islands by over 1,470 growers on ~10,000 acres, is conservatively valued at $100M per year (USDA NASS 2023). Until late October 2020, Hawaii was the only major coffee producing region of the world that was free of Coffee Leaf Rust (CLR). Growers are currently facing their most formidable production challenge with the arrival of Hemileia vastatrix Berk. & Broome, the most economically devastating pathogen of coffee worldwide. Since its introduction (Keith et al. 2022), CLR has rapidly spread throughout the state and can be found on coffee farms and feral coffee throughout the six islands. Implementation of CLR control measures will be difficult in Hawaii, given the extreme environmental heterogeneity, differences in management practices, high production costs, and labor shortages. Compounding these challenges is that all coffee genotypes grown on a large scale in the state are susceptible to CLR. More than 55 different rust races from coffee growing countries worldwide have been identified (Silva et al. 2022). Since key control measures include developing and establishing resistant coffee cultivars, determining the rust race(s) present in Hawaii was imperative. In June 2021, nine spore samples from symptomatic cultivated and feral plants ('Typica') growing on three islands (Hawaii Island: 3, Maui: 5, Molokai: 1) were collected in gelatin capsules using a G-R Electric Manufacturing portable vacuum pump with a mini cyclone spore adapter. The samples were sent to the Coffee Rust Research Center (CIFC) in Portugal. At CIFC, the urediniospores were bulked on susceptible genotype 849/1 Matari and inoculated on a set of coffee differentials following a standard race-typing procedure (Várzea and Marques 2005). The genotype of virulence of rust samples was inferred according to Flor's gene-for-gene theory (Silva et al. 2022). The genes of virulence v2, v4, and v5 (Race XXIV) were identified in all rust samples from all islands in Hawaii, supporting the theory of a single introduction to the state, which subsequently spread (Ramírez-Camejo et al. 2022). Race XXIV was previously characterized at CIFC and is commonly found in the majority of coffee-growing countries in South and Central America, Africa and Asia (CIFC's data base). According to Figueiredo & Arruda (1974), race XXIV is considered highly aggressive with a high spore germination rate, medium germ tube length, and short incubation period required for infection. Race XXIV is pathogenic to all coffee Arabica genotypes with the resistance genes SH5 or SH2,5 like varieties Blue Mountain, Bourbon, Catuaí, Caturra, Kent's, Kona, K7, Mundo Novo, SL 28, SL 39, as well as Accession "Agaro" with resistance genes SH4,5 (CIFC's records). On the other hand, this race is not virulent to some other Arabica genotypes, such as Geisha (SH1,5), S.288 (SH3,5), and Dilla & Alghe (SH1). Race XXIV is unable to infect derivatives of interspecific tetraploid hybrids like the groups Catimor and Sarchimor (Bettencourt and Rodrigues 1988). This is the first report of race XXIV on Coffea arabica in Hawaii. This finding is essential to evaluate the potential resistance of coffee germplasm existing in Hawaii or to be introduced in this region to develop new varieties. Since the emergence of new H. vastatrix races occur preferentially at germplasm collections (Li et al. 2021), proper management is imperative where multiple genotypes/varieties are planted.

The story of coffee: legend and truth.

When we think about coffee, exotic tropical countries such as Colombia, Brazil, and Ethiopia first come to mind. However, the crucial contribution of Portugal and its scientists to each cup of coffee we drink remains either poorly known or overlooked.

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.

Pathological, Morphological, Cytogenomic, Biochemical and Molecular Data Support the Distinction between Colletotrichum cigarro comb. et stat. nov. and Colletotrichum kahawae.

The genus Colletotrichum has witnessed tremendous variations over the years in the number of species recognized, ranging from 11 to several hundreds. Host-specific fungal species, once the rule, are now the exception, with polyphagous behavior regarded as normal in this genus. The species Colletotrichum kahawae was created to accommodate the pathogens that have the unique ability to infect green developing coffee berries causing the devastating Coffee Berry Disease in Africa, but its close phylogenetic relationship to a polyphagous group of fungi in the C. gloeosporioides species complex led some researchers to regard these pathogens as members of a wider species. In this work we combine pathological, morphological, cytogenomic, biochemical, and molecular data of a comprehensive set of phylogenetically-related isolates to show that the Coffee Berry Disease pathogen forms a separate species, C. kahawae, and also to assign the closely related fungi, previously in C. kahawae subsp. cigarro, to a new species, C. cigarro comb. et stat. nov. This taxonomic clarification provides an opportunity to link phylogeny and functional biology, and additionally enables a much-needed tool for plant pathology and agronomy, associating exclusively C. kahawae to the Coffee Berry Disease pathogen.

Genome-Wide Signatures of Selection in Colletotrichum kahawae Reveal Candidate Genes Potentially Involved in Pathogenicity and Aggressiveness.

Plants and their pathogens are engaged in continuous evolutionary battles, with pathogens evolving to circumvent plant defense mechanisms and plants responding through enhanced protection to prevent or mitigate damage induced by pathogen attack. Managed ecosystems are composed of genetically identical populations of crop plants with few changes from year to year. These environments are highly conducive to the emergence and dissemination of pathogens and they exert selective pressure for both qualitative virulence factors responsible for fungal pathogenicity, and quantitative traits linked to pathogen fitness, such as aggressiveness. In this study, we used a comparative genome-wide approach to investigate the genomic basis underlying the pathogenicity and aggressiveness of the fungal coffee pathogen Colletotrichum kahawae infecting green coffee berries. The pathogenicity was investigated by comparing genomic variation between C. kahawae and its non-pathogenic sibling species, while the aggressiveness was studied by a genome-wide association approach with groups of isolates with different phenotypic profiles. High genetic differentiation was observed between C. kahawae and the most closely related species with 5,560 diagnostic SNPs identified, in which a significant enrichment of non-synonymous mutations was detected. Functional annotation of these non-synonymous mutations revealed a significant enrichment mainly in two gene ontology categories, "oxidation-reduction process" and "integral component of membrane." Finally, the annotation of several genes potentially under-selection revealed that C. kahawae's pathogenicity may be a complex biological process, in which important biological functions, such as, detoxification and transport, regulation of host and pathogen gene expression, and signaling are involved. On the other hand, the genome-wide association analyses for aggressiveness were able to identify 10 SNPs and 15 SNPs of small effect in single and multi-association analysis, respectively, from which 7 were common, giving in total 18 SNPs potentially associated. The annotation of these genomic regions allowed the identification of four candidate genes encoding F-box domain-containing, nitrosoguanidine resistance, Fungal specific transcription factor domain-containing and C6 transcription factor that could be associated with aggressiveness. This study shed light, for the first time, on the genetic mechanisms of C. kahawae host specialization.

Population genomic footprints of host adaptation, introgression and recombination in coffee leaf rust.

Coffee leaf rust, caused by Hemileia vastatrix (Hv), represents the biggest threat to coffee production worldwide and ranks amongst the most serious fungal diseases in history. Despite a recent series of outbreaks and emergence of hypervirulent strains, the population evolutionary history and potential of this pathogen remain poorly understood. To address this issue, we used restriction site-associated DNA sequencing (RADseq) to generate ∼19 000 single nucleotide polymorphisms (SNPs) across a worldwide collection of 37 Hv samples. Contrary to the long-standing idea that Hv represents a genetically unstructured and cosmopolitan species, our results reveal the existence of a cryptic species complex with marked host tropism. Using phylogenetic and pathological data, we show that one of these lineages (C3) infects almost exclusively the most economically valuable coffee species (tetraploids that include Coffea arabica and interspecific hybrids), whereas the other lineages (C1 and C2) are severely maladapted to these hosts, but successfully infect diploid coffee species. Population dynamic analyses suggest that the C3 group may be a recent 'domesticated' lineage that emerged via host shift from diploid coffee hosts. We also found evidence of recombination occurring within this group, which could explain the high pace of pathotype emergence despite the low genetic variation. Moreover, genomic footprints of introgression between the C3 and C2 groups were discovered and raise the possibility that virulence factors may be quickly exchanged between groups with different pathogenic abilities. This work advances our understanding of the evolutionary strategies used by plant pathogens in agro-ecosystems with direct and far-reaching implications for disease control.

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.

The coffee leaf rust pathogen Hemileia vastatrix: one and a half centuries around the tropics.

TAXONOMY AND HISTORY: Hemileia vastatrix Berk. and Broome (Basidiomycota, Pucciniales) was described in 1869 in eastern Africa and Ceylon as the agent of coffee leaf rust and has spread to all coffee cultivation areas worldwide. Major disease outbreaks in Asia, Africa and America caused and continue to cause severe yield losses, making this the most important disease of Arabica coffee, a cash crop for many tropical and sub-tropical countries. LIFE CYCLE AND DISEASE SYMPTOMS: Hemileia vastatrix is a hemicyclic fungus with the urediniosporic life cycle as its most important (if not only) source of inoculum. Chlorotic spots are the first macroscopic symptoms, preceding the differentiation of suprastomatal, bouquet-shaped, orange-coloured uredinia. The disease can cause yield losses of up to 35% and have a polyetic epidemiological impact on subsequent years. DISEASE CONTROL: Although the use of fungicides is one of the preferred immediate control measures, the use of resistant cultivars is considered to be the most effective and durable disease control strategy. The discovery of 'Híbrido de Timor' provided sources of resistance that have been used in several breeding programmes and that have been proven to be effective and durable, as some have been in use for more than 30 years. GENETIC DIVERSITY AND MOLECULAR PATHOGENICITY: Although exhibiting limited genetic polymorphism, the very large genome of H. vastatrix (c. 797 Mbp) conceals great pathological diversity, with more than 50 physiological races. Gene expression studies have revealed a very precocious activation of signalling pathways and production of putative effectors, suggesting that the plant-fungus dialogue starts as early as at the germ tube stage, and have provided clues for the identification of avr genes.

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.

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

BACKGROUND: Coffee 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. RESULTS: In 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. CONCLUSION: Taken 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.

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.

Application of the Apn2/MAT locus to improve the systematics of the Colletotrichum gloeosporioides complex: an example from coffee (Coffea spp.) hosts.

To improve phylogenetic resolution of the Colletotrichum gloeosporioides species complex we developed and tested the performance of a new set of primers for the Apn2/MAT locus with a case study of 22 isolates. These were isolated mainly from coffee plants and represent six divergent and well characterized species within the C. gloeosporioides complex. Following previous studies on this locus, we have generated sequence data from an expanded region (> 4600 bp), revealing increased phylogenetic informativeness when compared to other commonly used markers such as ITS, ß-tub2 and GS. Within the Apn2/MAT locus the ApMAT marker alone was almost as informative in terms of phylogenetic resolution as a seven-gene concatenated dataset. Our results further revealed that gene-tree discordance may come to be a common issue in resolving evolutionary relationships in the C. gloeosporioides complex, highlighting the importance of multilocus approaches. The use of state-of-the-art data analysis techniques and a highly informative dataset as employed here may abate this issue and hopefully assist in disentangling the C. gloeosporioides complex.

Genetic diversity in Hemileia vastatrix based on RAPD markers.

Random amplified polymorphic DNA (RAPD) was used to assess the genetic structure of Hemileia vastatrix populations. Forty-five rust isolates with different virulence spectra and from different hosts and geographical regions were analyzed. Out of 45 bands, generated with three RAPD primers, 35 (78%) were polymorphic and scored as molecular markers. Cluster analysis exhibits unstructured variability of this pathogen with regard to physiological race, geographical origin or host. The genotypic diversity (H') inferred from Shannon's index was higher than gene diversity (Ht), suggesting that diversity is distributed among clonal lineages. Estimates of gene diversity in Africa and Asia populations were higher in total (Ht) as compared to within population diversity (Hs). Genetic differentiation was considerable among coffee rust isolates from Africa (Gst = 0.865) and Asia (Gst = 0.768) but not among isolates from South America (Gst = 0.266). We concluded that genetic diversity in H. vastatrix was moderately low and that the genetic differentiation among populations shows that asexual reproduction is likely to play an important role in the population biology of this fungus. This should be taken into account for the development of breeding programs.