Genome sequencing and comparative genomics reveal a repertoire of putative pathogenicity genes in chilli anthracnose fungus Colletotrichum truncatum

Unveiling the Arsenal of Apple Bitter Rot Fungi: Comparative Genomics Identifies Candidate Effectors, CAZymes, and Biosynthetic Gene Clusters in Colletotrichum Species

The bitter rot of apple is caused by Colletotrichum spp. and is a serious pre-harvest disease that can manifest in postharvest losses on harvested fruit. In this study, we obtained genome sequences from four different species, C. chrysophilum, C. noveboracense, C. nupharicola, and C. fioriniae, that infect apple and cause diseases on other fruits, vegetables, and flowers. Our genomic data were obtained from isolates/species that have not yet been sequenced and represent geographic-specific regions. Genome sequencing allowed for the construction of phylogenetic trees, which corroborated the overall concordance observed in prior MLST studies. Bioinformatic pipelines were used to discover CAZyme, effector, and secondary metabolic (SM) gene clusters in all nine Colletotrichum isolates. We found redundancy and a high level of similarity across species regarding CAZyme classes and predicted cytoplastic and apoplastic effectors. SM gene clusters displayed the most diversity in type and the most common cluster was one that encodes genes involved in the production of alternapyrone. Our study provides a solid platform to identify targets for functional studies that underpin pathogenicity, virulence, and/or quiescence that can be targeted for the development of new control strategies. With these new genomics resources, exploration via omics-based technologies using these isolates will help ascertain the biological underpinnings of their widespread success and observed geographic dominance in specific areas throughout the country.

Detection and Quantification of Anthracnose Pathogen Colletotrichum fructicola in Cultivated Tea-Oil Camellia Species from Southern China Using a DNA-Based qPCR Assay

Tea-oil Camellia species as edible-oil producing trees are widely cultivated in southern China. Camellia anthracnose that is mainly caused by Colletotrichum fructicola is a major disease of tea-oil trees. However, rapid detection and precise quantification of C. fructicola in different Camellia species that are crucial for the fundamental study of this pathosystem and effective disease management remain largely unexplored. Here, we developed a sensitive, rapid, and accurate method for quantifying C. fructicola growth in different Camellia species using a quantitative PCR assay. Amplified C. fructicola DNA using ITS-specific primers is relatively compared with the amplification of Camellia oleifera using the TUB gene. We determined that the fungal growth is tightly associated with the disease development in Ca. oleifera following C. fructicola infection in a time-course manner. This assay is highly sensitive, as fungal growth was detected in six different inoculated tea-oil Camellia species without visible disease lesion symptoms. Additionally, this method was validated by quantifying the Camellia anthracnose in orchards that did not show any disease symptoms. This assay enables the rapid, highly sensitive, and precise detection and quantification of C. fructicola growth in different tea-oil Camellia species, which will have a practical application for early diagnosis of anthracnose disease under asymptomatic conditions in Camellia breeding and field and will facilitate the development of tea-oil trees and C. fructicola interaction as a mold system to study woody plant and fungal pathogens interaction.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Comparative transcriptomic provides novel insights into the soybean response to Colletotrichum truncatum infection

Introduction

Soybean (Glycine max) is among the most important crops in the world, and its production can be threatened by biotic diseases, such as anthracnose. Soybean anthracnose is a seed-borne disease mainly caused by the hemibiotrophic fungus Colletotrichum truncatum. Typical symptoms are pre- and post-emergence damping off and necrotic lesions on cotyledons, petioles, leaves, and pods. Anthracnose symptoms can appear early in the field, causing major losses to soybean production.

Material and Methods

In preliminary experiments, we observed that the same soybean cultivar can have a range of susceptibility towards different strains of C. truncatum, while the same C. truncatum strain can cause varying levels of disease severity in different soybean cultivars. To gain a better understanding of the molecular mechanisms regulating the early response of different soybean cultivars to different C. truncatum strains, we performed pathogenicity assays to select two soybean cultivars with significantly different susceptibility to two different C. truncatum strains and analyzed their transcriptome profiles at different time points of interaction (0, 12, 48, and 120 h post-inoculation, hpi).

Results and Discussion

The pathogenicity assays showed that the soybean cultivar Gm1 is more resistant to C. truncatum strain 1080, and it is highly susceptible to strain 1059, while cultivar Gm2 shows the opposite behavior. However, if only trivial anthracnose symptoms appeared in the more resistant phenotype (MRP; Gm1-1080; Gm2-1059) upon 120 hpi, in the more susceptible phenotype (MSP; Gm-1059; Gm2- 1080) plants show mild symptoms already at 72 hpi, after which the disease evolved rapidly to severe necrosis and plant death. Interestingly, several genes related to different cellular responses of the plant immune system (pathogen recognition, signaling events, transcriptional reprogramming, and defense-related genes) were commonly modulated at the same time points only in both MRP. The list of differentially expressed genes (DEGs) specific to the more resistant combinations and related to different cellular responses of the plant immune system may shed light on the important host defense pathways against soybean anthracnose.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Genomic Characteristics and Comparative Genomics Analysis of Parafenestella ontariensis sp. nov

A new ascomycetous species of Parafenestella was isolated from Acer negundo during the survey of diseased trees in Southern Ontario, Canada. The species is morphologically similar to other taxa of Cucurbitariacea (Pleosporales). The new species is different from the extant species in the morphology of ascospores, culture characteristics and molecular data. The novel species is described as Parafenestella ontariensis sp. nov. based on morphological and multi-gene phylogenetic analyses using a combined set of ITS, LSU, tef1 and tub2 loci. Additionally, the genome of P. ontariensis was sequenced and analyzed. The phylogenomic analysis confirmed the close relationship of the species to the fenestelloid clades of Cucurbitariaceae. The comparative genomics analysis revealed that the species lifestyle appears to be multitrophic (necrotrophic or hemi-biotrophic) with a capability to turn pathogenic on a corresponding plant host.

Genome-wide comparison deciphers lifestyle adaptation and glass biodeterioration property of Curvularia eragrostidis C52

Glass biodeterioration by fungi has caused irreversible damage to valuable glass materials such as cultural heritages and optical devices. To date, knowledge about metabolic potential and genomic profile of biodeteriorative fungi is still scarce. Here, we report for the first time the whole genome sequence of Curvularia eragrostidis C52 that strongly degraded silica-based glasses coated with fluorine and hafnium, as expressed by the hyphal surface coverage of 46.16 ± 3.3% and reduced light transmission of 50.93 ± 1.45%. The genome of C. eragrostidis C52 is 36.9 Mb long with a GC content of 52.1% and contains 14,913 protein-coding genes, which is the largest genome ever recorded in the genus Curvularia. Phylogenomic analysis revealed C. eragrostidis C52 formed a distinct cluster with Curvularia sp. IFB-Z10 and was not evolved from compared genomes. Genome-wide comparison showed that strain C52 harbored significantly higher proportion of proteins involved in carbohydrate-active enzymes, peptidases, secreted proteins, and transcriptional factors, which may be potentially attributed to a lifestyle adaptation. Furthermore, 72 genes involved in the biosynthesis of 6 different organic acids were identified and expected to be crucial for the fungal survival in the glass environment. To form biofilm against stress, the fungal strain utilized 32 genes responsible for exopolysaccharide production. These findings will foster a better understanding of the biology of C. eragrostidis and the mechanisms behind fungal biodeterioration in the future.

Structured Framework and Genome Analysis of Magnaporthe grisea Inciting Pearl Millet Blast Disease Reveals Versatile Metabolic Pathways, Protein Families, and Virulence Factors

Magnaporthe grisea (T.T. Herbert) M.E. Barr is a major fungal phytopathogen that causes blast disease in cereals, resulting in economic losses worldwide. An in-depth understanding of the basis of virulence and ecological adaptation of M. grisea is vital for devising effective disease management strategies. Here, we aimed to determine the genomic basis of the pathogenicity and underlying biochemical pathways in Magnaporthe using the genome sequence of a pearl millet-infecting M. grisea PMg_Dl generated by dual NGS techniques, Illumina NextSeq 500 and PacBio RS II. The short and long nucleotide reads could be draft assembled in 341 contigs and showed a genome size of 47.89 Mb with the N50 value of 765.4 Kb. Magnaporthe grisea PMg_Dl showed an average nucleotide identity (ANI) of 86% and 98% with M. oryzae and Pyricularia pennisetigena, respectively. The gene-calling method revealed a total of 10,218 genes and 10,184 protein-coding sequences in the genome of PMg_Dl. InterProScan of predicted protein showed a distinct 3637 protein families and 695 superfamilies in the PMg_Dl genome. In silico virulence analysis revealed the presence of 51VFs and 539 CAZymes in the genome. The genomic regions for the biosynthesis of cellulolytic endo-glucanase and beta-glucosidase, as well as pectinolytic endo-polygalacturonase, pectin-esterase, and pectate-lyases (pectinolytic) were detected. Signaling pathways modulated by MAPK, PI3K-Akt, AMPK, and mTOR were also deciphered. Multicopy sequences suggestive of transposable elements such as Type LTR, LTR/Copia, LTR/Gypsy, DNA/TcMar-Fot1, and Type LINE were recorded. The genomic resource presented here will be of use in the development of molecular marker and diagnosis, population genetics, disease management, and molecular taxonomy, and also provide a genomic reference for ascomycetous genome investigations in the future.

Proteinaceous Effector Discovery and Characterization in Plant Pathogenic Colletotrichum Fungi

Anthracnose caused by plant pathogenic Colletotrichum fungi results in large economic losses in field crop production worldwide. To aid the establishment of plant host infection, Colletotrichum pathogens secrete numerous effector proteins either in apoplastic space or inside of host cells for effective colonization. Understanding these effector repertoires is critical for developing new strategies for resistance breeding and disease management. With the advance of genomics and bioinformatics tools, a large repertoire of putative effectors has been identified in Colletotrichum genomes, and the biological functions and molecular mechanisms of some studied effectors have been summarized. Here, we review recent advances in genomic identification, understanding of evolutional characteristics, transcriptional profiling, and functional characterization of Colletotrichum effectors. We also offer a perspective on future research.

Fruit Fly Larval Survival in Picked and Unpicked Tomato Fruit of Differing Ripeness and Associated Gene Expression Patterns

The larvae of frugivorous tephritid fruit flies feed within fruit and are global pests of horticulture. With the reduced use of pesticides, alternative control methods are needed, of which fruit resistance is one. In the current study, we explicitly tested for phenotypic evidence of induced fruit defences by running concurrent larval survival experiments with fruit on or off the plant, assuming that defence induction would be stopped or reduced by fruit picking. This was accompanied by RT-qPCR analysis of fruit defence and insect detoxification gene expression. Our fruit treatments were picking status (unpicked vs. picked) and ripening stage (colour break vs. fully ripe), our fruit fly was the polyphagous Bactrocera tryoni, and larval survival was assessed through destructive fruit sampling at 48 and 120 h, respectively. The gene expression study targeted larval and fruit tissue samples collected at 48 h and 120 h from picked and unpicked colour-break fruit. At 120 h in colour-break fruit, larval survival was significantly higher in the picked versus unpicked fruit. The gene expression patterns in larval and plant tissue were not affected by picking status, but many putative plant defence and insect detoxification genes were upregulated across the treatments. The larval survival results strongly infer an induced defence mechanism in colour-break tomato fruit that is stronger/faster in unpicked fruits; however, the gene expression patterns failed to provide the same clear-cut treatment effect. The lack of conformity between these results could be related to expression changes in unsampled candidate genes, or due to critical changes in gene expression that occurred during the unsampled periods.

Carbohydrate active enzymes (CAZy) regulate cellulolytic and pectinolytic enzymes in Colletotrichum falcatum causing red rot in sugarcane

Colletotrichum falcatum, an ascomycete pathogen causes red rot of sugarcane which is specialized to infect cane stalks. Cellulolytic and pectinolytic enzymes are necessary for degradation of plant cell wall which stands as barrier for successful fungal pathogenesis. In the study, we have confined to the CAZy genes that regulate cellulolytic and pectinolytic enzymes in two distinctive pathotypes of C. falcatum. Comparative transcriptome analysis revealed that a number of CAZy genes producing cellulolytic and pectinolytic enzyme were present in the virulent (Cf671) and least virulent (RoC) pathotypes. Two consecutive transcriptome analyses (in vitro) were performed using Illumina Hi Seq 2500, further analysis was done with various bioinformatic tools. In vitro expression analysis of cutinase, glycoside hydrolyase and pectin-related genes revealed number of genes that attributes virulence. Numerous pectin-related genes involved in degradation of plant cell wall, pectinase and pectin lyase are considered to be key precursor in degradation of pectin in sugarcane. These results suggest that cellulolytic enzymes, cutinase and pectin-related genes are essential for degradation of sugarcane cell wall and considered to be an important pathogenic factor in C. falcatum. This is the first detailed report on sugarcane cell wall-degrading enzymes during its interaction with C. falcatum and also this comparative transcriptome analysis provided more insights into pathogen mechanism on C. falcatum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03113-6.

Gene selection for studying frugivore-plant interactions: a review and an example using Queensland fruit fly in tomato

Fruit production is negatively affected by a wide range of frugivorous insects, among them tephritid fruit flies are one of the most important. As a replacement for pesticide-based controls, enhancing natural fruit resistance through biotechnology approaches is a poorly researched but promising alternative. The use of quantitative reverse transcription PCR (RT-qPCR) is an approach to studying gene expression which has been widely used in studying plant resistance to pathogens and non-frugivorous insect herbivores, and offers a starting point for fruit fly studies. In this paper, we develop a gene selection pipe-line for known induced-defense genes in tomato fruit, Solanum lycopersicum, and putative detoxification genes in Queensland fruit fly, Bactrocera tryoni, as a basis for future RT-qPCR research. The pipeline started with a literature review on plant/herbivore and plant/pathogen molecular interactions. With respect to the fly, this was then followed by the identification of gene families known to be associated with insect resistance to toxins, and then individual genes through reference to annotated B. tryoni transcriptomes and gene identity matching with related species. In contrast for tomato, a much better studied species, individual defense genes could be identified directly through literature research. For B. tryoni, gene selection was then further refined through gene expression studies. Ultimately 28 putative detoxification genes from cytochrome P450 (P450), carboxylesterase (CarE), glutathione S-transferases (GST), and ATP binding cassette transporters (ABC) gene families were identified for B. tryoni, and 15 induced defense genes from receptor-like kinase (RLK), D-mannose/L-galactose, mitogen-activated protein kinase (MAPK), lipoxygenase (LOX), gamma-aminobutyric acid (GABA) pathways and polyphenol oxidase (PPO), proteinase inhibitors (PI) and resistance (R) gene families were identified from tomato fruit. The developed gene selection process for B. tryoni can be applied to other herbivorous and frugivorous insect pests so long as the minimum necessary genomic information, an annotated transcriptome, is available.

Genomic Characteristics and Comparative Genomics Analysis of Two Chinese Corynespora cassiicola Strains Causing Corynespora Leaf Fall (CLF) Disease

Rubber tree Corynespora leaf fall (CLF) disease, caused by the fungus Corynespora cassiicola, is one of the most damaging diseases in rubber tree plantations in Asia and Africa, and this disease also threatens rubber nurseries and young rubber plantations in China. C. cassiicola isolates display high genetic diversity, and virulence profiles vary significantly depending on cultivar. Although one phytotoxin (cassicolin) has been identified, it cannot fully explain the diversity in pathogenicity between C. cassiicola species, and some virulent C. cassiicola strains do not contain the cassiicolin gene. In the present study, we report high-quality gapless genome sequences, obtained using short-read sequencing and single-molecule long-read sequencing, of two Chinese C. cassiicola virulent strains. Comparative genomics of gene families in these two stains and a virulent CPP strain from the Philippines showed that all three strains experienced different selective pressures, and metabolism-related gene families vary between the strains. Secreted protein analysis indicated that the quantities of secreted cell wall-degrading enzymes were correlated with pathogenesis, and the most aggressive CCP strain (cassiicolin toxin type 1) encoded 27.34% and 39.74% more secreted carbohydrate-active enzymes (CAZymes) than Chinese strains YN49 and CC01, respectively, both of which can only infect rubber tree saplings. The results of antiSMASH analysis showed that all three strains encode ~60 secondary metabolite biosynthesis gene clusters (SM BGCs). Phylogenomic and domain structure analyses of core synthesis genes, together with synteny analysis of polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) gene clusters, revealed diversity in the distribution of SM BGCs between strains, as well as SM polymorphisms, which may play an important role in pathogenic progress. The results expand our understanding of the C. cassiicola genome. Further comparative genomic analysis indicates that secreted CAZymes and SMs may influence pathogenicity in rubber tree plantations. The findings facilitate future exploration of the molecular pathogenic mechanism of C. cassiicola.

Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi

Fungi produce a wealth of pharmacologically bioactive secondary metabolites (SMs) from biosynthetic gene clusters (BGCs). It is common practice for drug discovery efforts to treat species' secondary metabolomes as being well represented by a single or a small number of representative genomes. However, this approach misses the possibility that intraspecific population dynamics, such as adaptation to environmental conditions or local microbiomes, may harbor novel BGCs that contribute to the overall niche breadth of species. Using 94 isolates of Aspergillus flavus, a cosmopolitan model fungus, sampled from seven states in the United States, we dereplicate 7,821 BGCs into 92 unique BGCs. We find that more than 25% of pangenomic BGCs show population-specific patterns of presence/absence or protein divergence. Population-specific BGCs make up most of the accessory-genome BGCs, suggesting that different ecological forces that maintain accessory genomes may be partially mediated by population-specific differences in secondary metabolism. We use ultra-high-performance high-resolution mass spectrometry to confirm that these genetic differences in BGCs also result in chemotypic differences in SM production in different populations, which could mediate ecological interactions and be acted on by selection. Thus, our results suggest a paradigm shift that previously unrealized population-level reservoirs of SM diversity may be of significant evolutionary, ecological, and pharmacological importance. Last, we find that several population-specific BGCs from A. flavus are present in Aspergillus parasiticus and Aspergillus minisclerotigenes and discuss how the microevolutionary patterns we uncover inform macroevolutionary inferences and help to align fungal secondary metabolism with existing evolutionary theory.

Soybean anthracnose caused by Colletotrichum species: Current status and future prospects

Soybean (Glycine max) is one of the most important cultivated plants worldwide as a source of protein-rich foods and animal feeds. Anthracnose, caused by different lineages of the hemibiotrophic fungus Colletotrichum, is one of the main limiting factors to soybean production. Losses due to anthracnose have been neglected, but their impact may threaten up to 50% of the grain production.

TAXONOMY

While C. truncatum is considered the main species associated with soybean anthracnose, recently other species have been reported as pathogenic on this host. Until now, it has not been clear whether the association of new Colletotrichum species with the disease is related to emerging species or whether it is due to the undergoing changes in the taxonomy of the genus.

DISEASE SYMPTOMS

Typical anthracnose symptoms are pre- and postemergence damping-off; dark, depressed, and irregular spots on cotyledons, stems, petioles, and pods; and necrotic laminar veins on leaves that can result in premature defoliation. Symptoms may evolve to pod rot, immature opening of pods, and premature germination of grains.

CHALLENGES

As accurate species identification of the causal agent is decisive for disease control and prevention, in this work we review the taxonomic designation of Colletotrichum isolated from soybean to understand which lineages are pathogenic on this host. We also present a comprehensive literature review of soybean anthracnose, focusing on distribution, symptomatology, epidemiology, disease management, identification, and diagnosis. We consider the knowledge emerging from population studies and comparative genomics of Colletotrichum spp. associated with soybean providing future perspectives in the identification of molecular factors involved in the pathogenicity process.

USEFUL WEBSITE

Updates on Colletotrichum can be found at http://www.colletotrichum.org/. All available Colletotrichum genomes on GenBank can be viewed at http://www.colletotrichum.org/genomics/.

Genome Sequence Resources of Colletotrichum truncatum, C. plurivorum, C. musicola, and C. sojae: Four Species Pathogenic to Soybean (Glycine max)

Colletotrichum is a large genus of plant pathogenic fungi comprising more than 200 species. In this work, we present the genome sequences of four Colletotrichum species pathogenic to soybean: C. truncatum, C. plurivorum, C. musicola, and C. sojae. While C. truncatum is globally considered the most important pathogen, the other three species have been described and associated with soybean only recently. The genome sequences will provide insights into factors that contribute to pathogenicity toward soybean and will be useful for further research into the evolution of Colletotrichum.

Pathogenic Adaptations Revealed by Comparative Genome Analyses of Two Colletotrichum spp., the Causal Agent of Anthracnose in Rubber Tree

Colletotrichum siamense and Colletotrichum australisinense cause Colletotrichum leaf disease that differ in their symptoms in rubber tree (Hevea brasiliensis), and pathogenicity of these two fungal species is also not identical on different cultivars of rubber tree. This divergence is often attributed to pathogen virulence factors, namely carbohydrate-active enzymes (CAZymes), secondary metabolites (SM), and small-secreted protein (SSP) effectors. The draft genome assembly and functional annotation of potential pathogenicity genes of both species obtained here provide an important and timely genomic resource for better understanding the biology and lifestyle of Colletotrichum spp. This should pave the way for designing more efficient disease control strategies in plantations of rubber tree. In this study, the genes associated with these categories were manually annotated in the genomes of C. australisinense GX1655 and C. siamense HBCG01. Comparative genomic analyses were performed to address the evolutionary relationships among these gene families in the two species. First, the size of genome assembly, number of predicted genes, and some of the functional categories differed significantly between the two congeners. Second, from the comparative genomic analyses, we identified some specific genes, certain higher abundance of gene families associated with CAZymes, CYP450, and SM in the genome of C. siamense, and Nep1-like proteins (NLP) in the genome of C. australisinense.

Integrated transcriptomic and metabolomic analyses reveal the effects of callose deposition and multihormone signal transduction pathways on the tea plant-Colletotrichum camelliae interaction

Colletotrichum infects diverse hosts, including tea plants, and can lead to crop failure. Numerous studies have reported that biological processes are involved in the resistance of tea plants to Colletotrichum spp. However, the molecular and biochemical responses in the host during this interaction are unclear. Cuttings of the tea cultivar Longjing 43 (LJ43) were inoculated with a conidial suspension of Colletotrichum camelliae, and water-sprayed cuttings were used as controls. In total, 10,592 differentially expressed genes (DEGs) were identified from the transcriptomic data of the tea plants and were significantly enriched in callose deposition and the biosynthesis of various phytohormones. Subsequently, 3,555 mass spectra peaks were obtained by LC-MS detection in the negative ion mode, and 27, 18 and 81 differentially expressed metabolites (DEMs) were identified in the tea leaves at 12 hpi, 24 hpi and 72 hpi, respectively. The metabolomic analysis also revealed that the levels of the precursors and intermediate products of jasmonic acid (JA) and indole-3-acetate (IAA) biosynthesis were significantly increased during the interaction, especially when the symptoms became apparent. In conclusion, we suggest that callose deposition and various phytohormone signaling systems play important roles in the tea plant-C. camelliae interaction.

Colletotrichum: species complexes, lifestyle, and peculiarities of some sources of genetic variability

The genus Colletotrichum comprises species with different lifestyles but is mainly known for phytopathogenic species that infect crops of agronomic relevance causing considerable losses. The fungi of the genus Colletotrichum are distributed in species complexes and within each complex some species have particularities regarding their lifestyle. The most commonly found and described lifestyles in Colletotrichum are endophytic and hemibiotrophic phytopathogenic. Several of these phytopathogenic species show wide genetic variability, which makes long-term maintenance of resistance in plants difficult. Different mechanisms may play an important role in the emergence of genetic variants but are not yet fully understood in this genus. These mechanisms include heterokaryosis, a parasexual cycle, sexual cycle, transposable element activity, and repeat-induced point mutations. This review provides an overview of the genus Colletotrichum, the species complexes described so far and the most common lifestyles in the genus, with a special emphasis on the mechanisms that may be responsible, at least in part, for the emergence of new genotypes under field conditions.

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.

Comparative genome analysis indicates high evolutionary potential of pathogenicity genes in Colletotrichum tanaceti

Colletotrichum tanaceti is an emerging foliar fungal pathogen of commercially grown pyrethrum (Tanacetum cinerariifolium). Despite being reported consistently from field surveys in Australia, the molecular basis of pathogenicity of C. tanaceti on pyrethrum is unknown. Herein, the genome of C. tanaceti (isolate BRIP57314) was assembled de novo and annotated using transcriptomic evidence. The inferred putative pathogenicity gene suite of C. tanaceti comprised a large array of genes encoding secreted effectors, proteases, CAZymes and secondary metabolites. Comparative analysis of its putative pathogenicity gene profiles with those of closely related species suggested that C. tanaceti likely has additional hosts to pyrethrum. The genome of C. tanaceti had a high repeat content and repetitive elements were located significantly closer to genes inferred to influence pathogenicity than other genes. These repeats are likely to have accelerated mutational and transposition rates in the genome, resulting in a rapid evolution of certain CAZyme families in this species. The C. tanaceti genome showed strong signals of Repeat Induced Point (RIP) mutation which likely caused its bipartite nature consisting of distinct gene-sparse, repeat and A-T rich regions. Pathogenicity genes within these RIP affected regions were likely to have a higher evolutionary rate than the rest of the genome. This "two-speed" genome phenomenon in certain Colletotrichum spp. was hypothesized to have caused the clustering of species based on the pathogenicity genes, to deviate from taxonomic relationships. The large repertoire of pathogenicity factors that potentially evolve rapidly due to the plasticity of the genome, indicated that C. tanaceti has a high evolutionary potential. Therefore, C. tanaceti poses a high-risk to the pyrethrum industry. Knowledge of the evolution and diversity of the putative pathogenicity genes will facilitate future research in disease management of C. tanaceti and other Colletotrichum spp.

Differences in the Characteristics and Pathogenicity of Colletotrichum camelliae and C. fructicola Isolated From the Tea Plant [Camellia sinensis (L.) O. Kuntze]

Colletotrichum, the causative agent of anthracnose, is an important pathogen that invades the tea plant (Camellia sinensis). In this study, 38 isolates were obtained from the diseased leaves of tea plants collected in different areas of Zhejiang Province, China. A combination of multigene (ITS, ACT, GAPDH, TUB2, CAL, and GS) and morphology analyses showed that the 38 strains belonged to two different species, namely, C. camelliae (CC), and C. fructicola (CF). Pathogenicity tests revealed that CC was more invasive than CF. In vitro inoculation experiments demonstrated that CC formed acervuli at 72 hpi and developed appressoria on wound edges, but CF did not develop these structures. Under treatment with catechins and caffeine, the growth inhibition rates of CF were remarkably higher than those of CC, indicating that the nonpathogenic species CF was more vulnerable to catechins and caffeine. Growth condition testing indicated that CF grew at a wide temperature range of 15-35°C and that the optimum temperature for CC growth was 25°C. Growth of both CC and CF did not differ between acidic and weakly alkaline environments (pH 5-8), but the growth of CC was significantly reduced at pH values of 9 and 10. Furthermore, the PacC/RIM101 gene, which associated with pathogenicity, was identified from CC and CF genomes, and its expression was suppressed in the hyphae of both species under pH value of 5 and 10, and much lower expression level was detected in CC than that in CF at pH 6. These results indicated that temperature has more important effect than pH for the growth of two Colletotrichum species. In conclusion, the inhibition by secondary metabolite is an important reason why the pathogenicity by CC and CF are different to tea plant, although the environmental factors including pH and temperature effect the growth of two Colletotrichum species.

Inherent Resistance to 14α-Demethylation Inhibitor Fungicides in Colletotrichum truncatum Is Likely Linked to CYP51A and/or CYP51B Gene Variants

Anthracnose disease, caused by Colletotrichum truncatum, affects marketable yield during preharvest production and postharvest storage of fruits and vegetables worldwide. Demethylation inhibitor (DMI) fungicides are among the very few chemical classes of single-site mode of action fungicides that are effective in controlling anthracnose disease. However, some species are inherently resistant to DMIs and more information is needed to understand this phenomenon. Isolates of C. truncatum were collected from the United States and China from peach, soybean, citrus, and begonia and sensitivity to six DMIs (difenoconazole, propiconazole, metconazole, tebuconazole, flutriafol, and fenbuconazole) was determined. Compared with DMI sensitive isolates of C. fructicola, C. siamense, and C. fioriniae (EC50 value ranging from 0.03 to 16.2 µg/ml to six DMIs), C. truncatum and C. nymphaeae were resistant to flutriafol and fenbuconazole (with EC50 values of more 50 µg/ml). Moreover, C. truncatum was resistant to tebuconazole and metconazole (with resistance factors of 27.4 and 96.0) and displayed reduced sensitivity to difenoconazole and propiconazole (with resistance factors of 5.1 and 5.2). Analysis of the Colletotrichum spp. genome revealed two potential DMI targets, CYP51A and CYP51B, that putatively encode P450 sterol 14α-demethylases. Both genes were identified and sequenced from C. truncatum and other species and no correlation between CYP51 gene expression levels and fungicide sensitivity was found. Four amino acid variations L208Y, H238R, S302A, and I366L in CYP51A, and three variations H373 N, M376L, and S511T in CYP51B correlated with the DMI resistance phenotype. CYP51A structure model analysis suggested the four alterations may reduce azole affinity. Likewise, CYP51B structure analysis suggested the H373 N and M376L variants may change the conformation of the DMI binding pocket, thereby causing differential sensitivity to DMI fungicides in C. truncatum.

The Landscape of Repetitive Elements in the Refined Genome of Chilli Anthracnose Fungus Colletotrichum truncatum

The ascomycete fungus Colletotrichum truncatum is a major phytopathogen with a broad host range which causes anthracnose disease of chilli. The genome sequencing of this fungus led to the discovery of functional categories of genes that may play important roles in fungal pathogenicity. However, the presence of gaps in C. truncatum draft assembly prevented the accurate prediction of repetitive elements, which are the key players to determine the genome architecture and drive evolution and host adaptation. We re-sequenced its genome using single-molecule real-time (SMRT) sequencing technology to obtain a refined assembly with lesser and smaller gaps and ambiguities. This enabled us to study its genome architecture by characterising the repetitive sequences like transposable elements (TEs) and simple sequence repeats (SSRs), which constituted 4.9 and 0.38% of the assembled genome, respectively. The comparative analysis among different Colletotrichum species revealed the extensive repeat rich regions, dominated by Gypsy superfamily of long terminal repeats (LTRs), and the differential composition of SSRs in their genomes. Our study revealed a recent burst of LTR amplification in C. truncatum, C. higginsianum, and C. scovillei. TEs in C. truncatum were significantly associated with secretome, effectors and genes in secondary metabolism clusters. Some of the TE families in C. truncatum showed cytosine to thymine transitions indicative of repeat-induced point mutation (RIP). C. orbiculare and C. graminicola showed strong signatures of RIP across their genomes and "two-speed" genomes with extensive AT-rich and gene-sparse regions. Comparative genomic analyses of Colletotrichum species provided an insight into the species-specific SSR profiles. The SSRs in the coding and non-coding regions of the genome revealed the composition of trinucleotide repeat motifs in exons with potential to alter the translated protein structure through amino acid repeats. This is the first genome-wide study of TEs and SSRs in C. truncatum and their comparative analysis with six other Colletotrichum species, which would serve as a useful resource for future research to get insights into the potential role of TEs in genome expansion and evolution of Colletotrichum fungi and for development of SSR-based molecular markers for population genomic studies.