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.

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.

Novel insights on colonization routes and evolutionary potential of Colletotrichum kahawae, a severe pathogen of Coffea arabica.

Pathogenic fungi are emerging at an increasing rate on a wide range of host plants, leading to tremendous threats to the global economy and food safety. Several plant pathogens have been considered to be invasive species, rendering large-scale population genomic analyses crucial to better understand their demographic history and evolutionary potential. Colletotrichum kahawae (Ck) is a highly aggressive and specialized pathogen, causing coffee berry disease in Arabica coffee in Africa. This pathogen leads to severe production losses and its dissemination out of Africa is greatly feared. To address this issue, a population genomic approach using thousands of single nucleotide polymorphisms (SNPs) spaced throughout the genome was used to unveil its demographic history and evolutionary potential. The current study confirms that Ck is a true clonal pathogen, perfectly adapted to green coffee berries, with three completely differentiated populations (Angolan, Cameroonian and East African). Two independent clonal lineages were found within the Angolan population as opposed to the remaining single clonal populations. The most probable colonization scenario suggests that this pathogen emerged in Angola and immediately dispersed to East Africa, where these two populations began to differentiate, followed by the introduction in Cameroon from an Angolan population. However, the differentiation between the two Angolan clonal lineages masks the mechanism for the emergence of the Cameroonian population. Our results suggest that Ck is completely differentiated from the ancestral lineage, has a low evolutionary potential and a low dispersion ability, with human transport the most likely scenario for its potential dispersion, which makes the fulfilment of the quarantine measures and management practices implemented crucial.

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.

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.