Discerning the global phylogeographic distribution of Phyllosticta citricarpa by means of whole genome sequencing

Abstract

Phyllosticta citricarpa is a fungal pathogen causing citrus black spot (CBS). As a regulated pest in some countries, the presence of the pathogen limits the export of fruit and is therefore of agricultural and economic importance. In this study, we used high throughput sequencing data to infer the global phylogeographic distribution of this pathogen, including 71 isolates from eight countries, Argentina, Australia, Brazil, China, Cuba, Eswatini, South Africa and the United States of America. We assembled draft genomes and used a pairwise read mapping approach for the detection and enumeration of variants between isolates. We performed SSR marker discovery based on the assembled genome with the best assembly statistics, and generated genotype profiles for all isolates with 1987 SSR markers in silico. Furthermore, we identified 32,560 SNPs relative to a reference sequence followed by population genetic analyses based on the three datasets; pairwise variant counts, SSR genotypes and SNP genotypes. All three analysis approaches gave similar overall results. Possible pathways of dissemination among the populations from China, Australia, southern Africa and the Americas are postulated. The Chinese population is the most diverse, and is genetically the furthest removed from all other populations, and is therefore considered the closest to the origin of the pathogen. Isolates from Australia, Eswatini and the South African province Mpumalanga are closely associated and clustered together with those from Argentina and Brazil. The Eastern Cape, North West, and KwaZulu-Natal populations in South Africa grouped in another cluster, while isolates from Limpopo are distributed between the two aforementioned clusters. Southern African populations showed a close relationship to populations in North America, and could be a possible source of P. citricarpa populations that are now found in North America. This study represents the largest whole genome sequencing survey of P. citricarpa to date and provides a more comprehensive assessment of the population genetic diversity and connectivity of P. citricarpa from different geographic origins. This information could further assist in a better understanding of the epidemiology of the CBS pathogen, its long-distance dispersal and dissemination pathways, and can be used to refine phytosanitary regulations and management programmes for the disease.

Keywords: Citrus black spot; Genomic data analysis; High throughput sequencing; International trade; Phyllostictaceae; Phylogeographic distribution patterns; Population genetic differentiation.