# Genome-Wide Analysis of Corynespora cassiicola Leaf Fall Disease Putative Effectors

David Lopez, Sébastien Ribeiro, Philippe Label, Boris Fumanal, Jean-Stéphane Venisse, Annegret Kohler, Ricardo R. de Oliveira, Kurt Labutti, Anna Lipzen, Kathleen Lail, et al.

► To cite this version:

David Lopez, Sébastien Ribeiro, Philippe Label, Boris Fumanal, Jean-Stéphane Venisse, et al. Genome-Wide Analysis of Corynespora cassiicola Leaf Fall Disease Putative Effectors. Frontiers in Microbiology, 2018, 9, pp.1-21. 10.3389/fmicb.2018.00276. hal-01753929

HAL Id: hal-01753929

https://hal.science/hal-01753929v1

Submitted on 29 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY-SA 4.0 - Attribution - ShareAlike - International
License
<figure><img src="image_2.png" /></figure>

# Genome-Wide Analysis of *Corynespora cassiicola* Leaf Fall Disease Putative Effectors

David Lopez<fnref n="1" />, Sébastien Ribeiro<fnref n="1" />,<fnref n="2" />,³, Philippe Label<fnref n="1" />, Boris Fumanal<fnref n="1" />, Jean-Stéphane Venisse<fnref n="1" />, Annegret Kohler⁴, Ricardo R. de Oliveira⁵, Kurt Labutti⁶, Anna Lipzen⁶, Kathleen Lail⁶, Diane Bauer⁶, Robin A. Ohm⁶,⁷, Kerrie W. Barry⁶, Joseph Spatafora⁸, Igor V. Grigoriev⁶,⁹, Francis M. Martin⁴ and Valérie Pujade-Renaud<fnref n="1" />,<fnref n="2" />,³*

**OPEN ACCESS**

**Edited by:**
Gail Preston,
University of Oxford, United Kingdom

**Reviewed by:**
Guus Bakkeren,
Agriculture and Agri-Food Canada
(AAFC), Canada
Bryan Bailey,
Agricultural Research Service (USDA),
United States

**\*Correspondence:**
Valérie Pujade-Renaud
valerie.pujade-renaud@cirad.fr

**Specialty section:**
This article was submitted to Plant Microbe Interactions, a section of the journal Frontiers in Microbiology

**Received:** 28 July 2017
**Accepted:** 07 February 2018
**Published:** 02 March 2018

**Citation:**
Lopez D, Ribeiro S, Label P, Fumanal B, Venisse J-S, Kohler A, de Oliveira RR, Labutti K, Lipzen A, Lail K, Bauer D, Ohm RA, Barry KW, Spatafora J, Grigoriev IV, Martin FM and Pujade-Renaud V (2018) Genome-Wide Analysis of Corynespora cassiicola Leaf Fall Disease Putative Effectors. *Front. Microbiol.* 9:276.
doi: 10.3389/fmicb.2018.00276

*Corynespora cassiicola* is an Ascomycetes fungus with a broad host range and diverse life styles. Mostly known as a necrotrophic plant pathogen, it has also been associated with rare cases of human infection. In the rubber tree, this fungus causes the *Corynespora* leaf fall (CLF) disease, which increasingly affects natural rubber production in Asia and Africa. It has also been found as an endophyte in South American rubber plantations where no CLF outbreak has yet occurred. The *C. cassiicola* species is genetically highly diverse, but no clear relationship has been evidenced between phylogenetic lineage and pathogenicity. Cassiicolin, a small glycosylated secreted protein effector, is thought to be involved in the necrotrophic interaction with the rubber tree but some virulent *C. cassiicola* isolates do not have a cassiicolin gene. This study set out to identify other putative effectors involved in CLF. The genome of a highly virulent *C. cassiicola* isolate from the rubber tree (CCP) was sequenced and assembled. *In silico* prediction revealed 2870 putative effectors, comprising CAZymes, lipases, peptidases, secreted proteins and enzymes associated with secondary metabolism. Comparison with the genomes of 44 other fungal species, focusing on effector content, revealed a striking proximity with phylogenetically unrelated species (*Colletotrichum acutatum*, *Colletotrichum gloesporioides*, *Fusarium oxysporum*, *nectria hematococca*, and *Botrosphaeria dothidea*) sharing life style plasticity and broad host range. Candidate effectors involved in the compatible interaction with the rubber tree were identified by transcriptomic analysis. Differentially expressed genes included 92 putative effectors, among which cassiicolin and two other secreted singleton proteins. Finally, the genomes of 35 *C. cassiicola* isolates representing the genetic diversity of the species were sequenced and assembled, and putative effectors identified. At the intraspecific level, effector-based classification was found to be highly consistent with the phylogenomic trees. Identification of lineage-specific effectors is a key step toward understanding *C. cassiicola* virulence and host specialization mechanisms.

**Keywords:** *Corynespora cassiicola*, *Hevea brasiliensis*, plant-pathogens interaction, effectors, cassiicolin, genomics, gene expression

<footnote n="1">¹ Université Clermont Auvergne, Institut National de la Recherche Agronomique, UMR PIAF, Clermont-Ferrand, France,</footnote>

<footnote n="2">² CIRAD, UMR AGAP, Clermont-Ferrand, France, ³ AGAP, Université Montpellier, CIRAD, Institut National de la Recherche Agronomique, Montpellier SupAgro, Montpellier, France, ⁴ Institut National de la Recherche Agronomique, UMR INRA-Université de Lorraine "Interaction Arbres/Microorganismes," Champenoux, France, ⁵ Departamento de Agronomía, Universidade Estadual de Maringá, Maringá, Brazil, ⁶ United States Department of Energy Joint Genome Institute, Walnut Creek, CA, United States, ⁷ Department of Microbiology, Utrecht University, Utrecht, Netherlands, ⁸ Department of Botany and Plant Pathology, Oregon State University, Corvallis, OR, United States, ⁹ Department of Plant and Microbial Biology, University of California, Berkeley, Berkeley, CA, United States</footnote>
# INTRODUCTION

*Corynespora cassiicola* (Berk. & M. A. Curtis) is an Ascomycetes fungus responsible for diseases in a wide range of plants (Farr and Rossman, 2016), mainly in tropical and subtropical areas or greenhouses. It has also been found in nematodes (Carris et al., 1986), sponges (Zhao et al., 2015) and in rare cases of human infections (Mahgoub, 1969; Huang et al., 2010; Yamada et al., 2013; Yan et al., 2016). In the rubber tree, *C. cassiicola* causes the *Corynespora* leaf fall (CLF) disease, characterized by necrotic lesions on the leaves and massive defoliation in susceptible cultivars. CLF was first observed in Sierra Leone (Deighton, 1936). Initially limited to nurseries, CLF was recognized as a threat to rubber farming after a severe outbreak in Sri Lanka in the late 1980's, which led to the uprooting of more than 4,000 ha of the highly susceptible rubber cultivar RRIC103 (Liyanage et al., 1986). Since then, CLF has gradually spread to most rubber producing areas in Asia and Africa, impairing natural rubber production yields (Kuruvilla Jacob, 2006).

*Corynespora* is a polyphyletic genus in the large Pleosporales order, mostly composed of plant pathogens. *C. cassiicola* forms with *Corynespora smithi* a specific family within Pleosporales (Schoch et al., 2009). Like most species of the *Corynespora* genus, *C. cassiicola* is anamorphic.

*Corynespora cassiicola* isolates display various life styles, from endophyte (Collado et al., 1999; Gond et al., 2007; Promputtha et al., 2007; Suryanarayanan et al., 2011; Déon et al., 2012b) to saprophyte (Kingsland, 1985; Lee et al., 2004; Cai et al., 2006), but are mostly reported as necrotrophic pathogens.

The high genetic diversity of the species is largely documented (Silva et al., 1998, 2003; Atan and Hamid, 2003; Nghia et al., 2008; Dixon et al., 2009; Qi et al., 2011; Déon et al., 2014; Hieu et al., 2014; Shuib et al., 2015). We have previously analyzed the diversity of 129 *C. cassiicola* isolates from various hosts (including 71 from the rubber tree) and of various geographical origins, based on the polymorphism of four combined *loci*, as proposed by Dixon et al. (Dixon et al., 2009). We identified eight major clades (A-H), with no clear geographic or host speciation (Déon et al., 2014).

*Corynespora cassiicola* diversity is not structured by strict host specificities, but rather host specialization (Dixon et al., 2009; Hieu et al., 2014), which suggests the involvement of specialized effectors recognized by a limited range of compatible hosts. In addition, in the rubber tree, *C. cassiicola* virulence profiles vary significantly depending on the cultivars (Breton et al., 2000; Atan and Hamid, 2003; Nghia et al., 2008; Qi et al., 2009; Déon et al., 2014).

The only *C. cassiicola* effector characterized so far is cassiicolin. This necrotrophic toxin was purified from the filtrate of the highly virulent isolate CCP, collected from diseased rubber trees in the Philippines (Breton et al., 2000; de Lamotte et al., 2007). Cassiicolin is a 27 amino acid glycoprotein with six cysteines engaged in three disulphide bonds (Barthe et al., 2007). It is matured from a pre-protein encoded by the three-exon *Cas* gene (GenBank EF667973) (Déon et al., 2012a). *Cas* is transiently upregulated 24-48 h after inoculation of rubber tree leaves with isolate CCP, just before the first symptoms occur in

a susceptible cultivar, suggesting a role during the first phase of infection. Various cassiicolin isoforms (Cas1 to Cas6) were later identified by PCR screening of 70 isolates from various hosts and geographical origins. *Cas1* refers to the initially described *Cas* gene, from CCP. We thus adopted a typology based on (i) the genetic group (A-H) and (ii) the toxin class determined according to the *Cas* gene polymorphism (Déon et al., 2014). Isolates of type C/Cas1 were found to be generally more aggressive than those of other types, but further study will be necessary to describe the possible associations between genetic group, toxin class and virulence profiles. Notably, some *C. cassiicola* isolates without *Cas* gene (isolates of toxin class Cas0) were nevertheless found to be moderately virulent in the rubber tree, evidencing effectors other than cassiicolin. Phenotypic diversity among *C. cassiicola* isolates was also evidenced by measuring the leaf damage (electrolyte leakages) caused by the application of fungal culture filtrates (Tran et al., 2016). Both qualitative and quantitative variations among the secreted effectors, together with multiple plant targets or sensitivity factors, may explain the observed phenotypic diversity.

This study set out to identify the putative effectors involved in CLF in the rubber tree, from the genomes of *C. corynespora* isolates representative of the intraspecific diversity. The genome of the virulent reference isolate CCP was first sequenced, assembled, and annotated by the US Department of Energy Joint Genome Institute (JGI) as part of the 1,000 Fungal Genomes Project (Grigoriev et al., 2014). Putative effectors were predicted and compared with the effector repertoires of 44 other fungal species. Transcript profiling was conducted to identify the functional effectors differentially expressed during a compatible interaction with the rubber tree. Lastly, we sequenced and assembled 35 genomes of new *C. cassiicola* isolates to compare their respective effector repertoires. Our findings are discussed in relation to the phylogeny and known physiological specificities of the selected isolates, with a special focus on the cassiicolin effector.

# MATERIALS AND METHODS

## Biological Material

*Corynespora cassiicola* (Berk. & M.A. Curtis) isolate CCP, highly virulent in the rubber tree (Breton et al., 2000), was chosen as phytopathogenic reference for the species. Genomic DNA and RNA samples were prepared from CCP mycelium grown for 3 weeks in liquid Czapeck medium, at 26/25°C (day/night), 12/12h light cycle. This culture was set up from 7-day-old cultures grown on PDA medium (potato dextrose agar, DIFCO), at 26/25°C (day/night), in the dark. For long-term conservation, mycelium plugs were kept at -80°C in 20% glycerol.

Interspecies comparative genomic analysis involved the CCP genome (http://genome.jgi.doe.gov/Corca1/Corca1.home.html) along with 44 other fungal genomes available from the JGI MycoCosm Portal (Table S1). Genome assembly from UM591, *C. cassiicola* isolate of human origin, was obtained from the NCBI Whole Genome Shotgun project under the accession http://www.ncbi.nlm.nih.gov/nuccore/JAQF00000000.1.
Intraspecies comparative genomic analysis was carried out on CCP along with 35 other *C. cassiicola* strains previously isolated from various host plants (Déon et al., 2014; **Table 1**) with CCP as reference isolate.

Transcript profiling was conducted in the context of a compatible interaction between the susceptible *H. brasiliensis*

<table><caption>**TABLE 1 | (A)** Metrics of the CCP genome assembly (B) Repetitive DNA content and transposable elements.</caption><thead><tr><th colspan="2">A</th></tr></thead><tbody><tr><td colspan="2"><strong>Genome assembly</strong></td></tr><tr><td>Genome Assembly size (Mbp)</td><td>44.85</td></tr><tr><td>Sequencing read coverage depth</td><td>58.8x</td></tr><tr><td>No. of contigs</td><td>500</td></tr><tr><td>No. of scaffolds</td><td>244</td></tr><tr><td>No. of scaffolds &gt;= 2Kbp</td><td>173</td></tr><tr><td>Scaffold N50</td><td>8</td></tr><tr><td>Scaffold L50 (Mbp)</td><td>2.51</td></tr><tr><td>No. of gaps</td><td>256</td></tr><tr><td>% of scaffold length in gaps</td><td>0.3%</td></tr><tr><td>Three largest Scaffolds (Mbp)</td><td>4.76, 2.92, 2.67</td></tr><tr><td>Assembled_RNAseq_reads</td><td>98.5%</td></tr><tr><td>No. of gene models</td><td>17,167</td></tr><tr><td colspan="2"><strong>Gene models</strong></td></tr><tr><td>Gene length (bp)</td><td>1,610</td></tr><tr><td>Transcript length (bp)</td><td>1,480</td></tr><tr><td>Exon length (bp)</td><td>570</td></tr><tr><td>Intron length (bp)</td><td>84</td></tr><tr><td>Protein length (aa)</td><td>421</td></tr><tr><td>Exons per gene</td><td>2.60</td></tr></tbody></table>

<table><thead><tr><th colspan="3">B</th></tr><tr><th>Repeats and transposable elements</th><th>Number of elements</th><th>Length occupied</th></tr></thead><tbody><tr><td>SINEs</td><td>25</td><td>1,679 bp</td></tr><tr><td>MIRs</td><td>6</td><td>433 bp</td></tr><tr><td>LINEs</td><td>192</td><td>14,456 bp</td></tr><tr><td>LINE1</td><td>11</td><td>882 bp</td></tr><tr><td>LINE2</td><td>28</td><td>2,238 bp</td></tr><tr><td>L3/CR1</td><td>68</td><td>4,700 bp</td></tr><tr><td>LTR elements</td><td>2</td><td>106 bp</td></tr><tr><td>ERVL-MaLRs</td><td>1</td><td>64 bp</td></tr><tr><td>ERV_classII</td><td>1</td><td>42 bp</td></tr><tr><td>DNA elements</td><td>45</td><td>3,689 bp</td></tr><tr><td>hAT-Charlie</td><td>8</td><td>560 bp</td></tr><tr><td>Total interspersed repeats</td><td></td><td>19,930 bp</td></tr><tr><td>Small RNA</td><td>88</td><td>22,102 bp</td></tr><tr><td>Simple repeats</td><td>14,030</td><td>590,502 bp (1.32%)</td></tr><tr><td>Low complexity</td><td>1,929</td><td>98,961 bp (0.22%)</td></tr></tbody></table>

Scaffold N50, scaffold length at which 50% of total bases in the assembly are in scaffolds of that length or greater. Scaffold L50, number of scaffolds whose sum length is N50.

clone PB260 and *C. cassiicola* reference isolate CCP. Rubber trees were grown in a greenhouse at 28°C during the day (16 h) and 26°C at night (8 h), with 80% relative humidity. Leaves were sampled at physiological stage C according to Hallé and Martin (Hallé and Martin, 1968), and placed in large Petri dishes, on moist paper. Spore suspensions were obtained from 7-day-old mycelium cultures on PDA medium, washed with sterile water (4-5 ml per Petri dish), and filtered through a 100 µm cell strainer (Biologix). Inoculations were as follows: six 20 µl droplets of spore suspension at a concentration of 5,000 spores/ml were spotted on the abaxial face of the detached leaflets and incubated at 26°C, 60% relative humidity (Déon et al., 2012a). Two leaf discs (2.2 cm²) per leaflet were collected at the inoculation spots using a cork borer, 24 and 48 h after inoculation (from distinct leaflets), with three independent biological repeats for each time point. In addition, we used two independent suspensions of germinating spores (kept in water for 4 h after collection) as reference for differential gene expression analysis.

## Nucleic Acid Preparation and Sequencing

Mycelium was flash frozen and ground in liquid nitrogen. For the reference genome, total DNA was extracted from 500 mg of ground material according to the 1,000 Fungal Genomes Project recommended protocol (http://1000.fungalgenomes.org/home/protocols/high-quality-genomic-dna-extraction/). In brief, the protocol consisted in a crude extraction using CTAB buffer followed by a purification step using Qiagen 5,000/G Genomic Tips (Qiagen, Courtaboeuf, France). Quality control was ensured by agarose gel migration and by NanoDrop1000 spectrophotometer quality assay (Thermo Fisher Scientific, Wilmington, U.S.A.).

The *C. cassiicola* CCP genome and transcriptome were sequenced using Illumina, the former in a combination of fragment and long mate pair (LMP) libraries.

For the fragment library, 100 ng of DNA was sheared to 300 bp using the Covaris LE220 and size-selected using SPRI beads (Beckman Coulter). The fragments were treated with end repair, A-tailing, and ligation of Illumina-compatible adapters (IDT, Inc.,) using the KAPA-Illumina library creation kit (KAPA biosystems).

LMP was produced from 6 µg of DNA sheared using the Covaris g-TUBE™ (Covaris) and gel size selected for 4 kb. The sheared DNA was treated with end repair and ligated with biotinylated adapters containing loxP. The adapter ligated DNA fragments were circularized via recombination by a Cre excision reaction (NEB). The circularized DNA templates were then randomly sheared using the Covaris LE220 (Covaris). The sheared fragments were treated with end repair and A-tailing using the KAPA-Illumina library creation kit (KAPA biosystems) followed by immobilization of mate pair fragments on strepavidin beads (Invitrogen). Illumina-compatible adapters (IDT, Inc) were ligated to the mate pair fragments and 10 cycles of PCR were used to enrich for the final library (KAPA Biosystems).

Stranded cDNA libraries were generated using the Illumina Truseq Stranded RNA LT kit. mRNA was purified from 1 µg of total RNA using magnetic beads containing poly-T oligos. mRNA was fragmented and reverse-transcribed using random hexamers
and SSII (Invitrogen) followed by second strand synthesis. The fragmented cDNA was treated with end-pair, A-tailing, and adapter ligation, and 10 cycles of PCR were used to enrich for the final library.

The prepared library was quantified using KAPA Biosystem's next-generation sequencing library qPCR kit and run on a Roche LightCycler 480 real-time PCR instrument. The quantified library was then multiplexed with other libraries, and the pool of libraries was then prepared for sequencing on the Illumina HiSeq sequencing platform using a TruSeq paired-end cluster kit, v3, and Illumina's cBot instrument to generate a clustered flowcell for sequencing. Sequencing of the flowcell was performed on the Illumina HiSeq2000 sequencer using a TruSeq SBS sequencing kit, v3, following a 2 × 100 bp (LPM) or 2 × 150 bp (fragments and transcriptome) indexed run recipe.

For the other isolates, DNAseq library preparation (Illumina TruSeq v3, 2 × 100 bp) and sequencing was carried out at the GenoToul GeT Platform (INRA Auzeville, France) on the Illumina Hiseq2000 platform.

For the reference transcriptome (CCP mycelium grown in liquid culture), RNA was extracted from the flash-frozen mycelium as previously described (Chang et al., 1993) and sequenced by JGI using Illumina.

For RNA-Seq kinetic analysis, the collected material (leaf discs from inoculated leaves and germinating spores used as reference) was immediately flash-frozen. For the germinating spores, 60 Petri dishes were required to yield sufficient material. Spores were collected in water as described above and the suspension was maintained for 4 h at room temperature and concentrated by centrifugation (20 min at 6,000 g). The supernatant was discarded and the pellet was ground in liquid nitrogen. In both cases, total RNA was extracted as previously described (Chang et al., 1993). RNA-Seq library preparation (Illumina TruSeq v3, 2 × 150 bp) and sequencing was carried out at the GenoToul GeT Platform (INRA Auzeville, France) starting from 2 µg of total RNA, on the Illumina Hiseq2000 platform.

# In Silico Analyses

## Genome de Novo Assembly

For the reference genome, Illumina genomic reads from two libraries were filtered and assembled with AllPathsLG (Gnerre et al., 2011). RNA-seq data for each genome were de novo assembled into consensus sequences using Rnnotator (Martin et al., 2010). The genome assembly was annotated using the JGI Annotation Pipeline and made available via the JGI fungal genome portal MycoCosm (jgi.doe.gov/fungi; Grigoriev et al., 2014). The Whole Genome Shotgun project of the reference isolate CCP has been deposited at DDBJ/ENA/GenBank under the accession NSJI00000000. The version described in this paper is version NSJI010000000.

For additional strains, after quality control, short reads were assembled in contigs using Velvet v.1.2.10 (Zerbino and Birney, 2008) with k-mers hash tables of 57 and automatic coverage cutoff. Mapping back the short reads to the generated contigs validated assembly parameters, yielding >97% mapped and >90% properly paired reads with the BWA SAMPE algorithm (Li and Durbin, 2009).

Post-assembly/annotation analysis relied on a parallel implementation (Tange, 2011) hosted by an Ubuntu 12.10 x64 server running on a 32-thread Intel Xeon Dell PowerEdge R720 server.

The Whole Genome Shotgun project of the additional *C. cassiicola* strains (BioProject PRJNA428435) has been deposited at DDBJ/ENA/GenBank under the accessions POQW00000000 to POSE00000000. The versions described in this paper are versions POQW01000000 to POSE01000000.

## Functional Annotation

Transposable elements (TEs) and repeats were searched for in the *C. cassiicola* reference genome assembly by first using RepeatScout v1.05 for de novo identification of simple repeats and then RepeatMasker open-4.0.5 (Smit et al., 2013–2015) with RepBase Update 20140131, RMLib 20140131 and Dfam v1.3 (Wheeler et al., 2013) and NCBI rmblastn v.2.2.27+ (Altschul et al., 1990; Camacho et al., 2009).

Carbohydrate-active enzymes (CAZymes) were predicted in fungal genomes using predicted protein sequences as queries and CAZyme Analysis Toolkit (Park et al., 2010) and the CAZy database (Lombard et al., 2014) (http://www.cazy.org). Proteins were considered as putative CAZymes if detected by both the Blast (E = 1e-5) and Pfam methods.

Lipases were identified in fungal genomes using predicted protein sequences as queries and hmm profiles from the Lipase Engineering Database v3 (http://www.led.uni-stuttgart.de) and hmmsearch (HMMER 3.1b1; http://hmmer.org/) with an E-value inclusion threshold set at 0.01.

Proteases were identified in fungal genomes using predicted protein sequences as queries and the MEROPS v9.12 blast database (Rawlings et al., 2014) (http://merops.sanger.ac.uk/) with Blastp (E = 1e-4).

Genes putatively involved in the secondary metabolism were searched for: non-ribosomal peptide synthases (NRPSs), polyketide synthases (PKSs) and terpene synthases (TSs) were identified using predicted protein sequences as queries and a local distribution of antiSMASH v2.0 (Blin et al., 2013).

Proteins were predicted as potentially secreted when the following two conditions were satisfied: (i) absence of a transmembrane domain identified by our TMHMM v2.0c local server (Sonnhammer et al., 1998), and (ii) presence of a secretion peptide signal as predicted by SignalP v4.1 with the Eukaryote database (Petersen et al., 2011). Proteins potentially secreted through ER/Golgi-independent pathways were not taken into account in this study.

## Protein Clustering

Clusters of homologous proteins were computed using the latest available standalone version of OrthoMCL 2.0.9 (Li et al., 2003; Fischer et al., 2011). Protein sequences were first filtered using the default complexity setting (orthomclFilterFasta <PATH> 10 20) yielding 578,863 proteins concatenated into a single blast database. All-versus-all blastp was then carried out with NCBI Blastp v2.2.30+ with the E-value filter set at 1e-5. MCL (Enright et al., 2002) was then used with the inflation parameter set at 1.5. We determined this value experimentally using a manually
curated cluster of orthologous major intrinsic proteins (MIPs), X-intrinsic proteins (XIPs). Using inflation set at 1.5, we obtained the least number of protein clusters while keeping XIPs distinct from other MIP families. Cluster dynamics was assessed using CAFE v3.1 (Han et al., 2013).

## Phylogenies

To build an interspecific phylogeny of the 45 fungal species of this study, we used predicted protein sequences found in exactly one copy in all fungal representatives, on the basis of cluster analysis. These 651 protein sequences of each fungal representative were individually aligned using MAFFT v7.215 (Katoh and Standley, 2013) with the default algorithm. Informative residues were extracted using Gblocks 0.91b (Castresana, 2000) using default settings considered by the author as stringent (minimum number of sequences for a conserved position 9; minimum number of sequences for a flank position 14; maximum number of contiguous non-conserved positions 8; minimum length of a block 10; gaps not allowed). After trimming, individual protein alignments were tested to obtain the best evolutionary model according to AIC scores using ProtTest v3.4 (Guindon and Gascuel, 2003; Darriba et al., 2011). We found that the LG model was the best substitution matrix with our data, and so it was selected for maximum likelihood tree reconstruction using RAxML v8.1.16 (Stamatakis, 2014) with 1,000 bootstraps on the concatenated chimeric protein sequence. Interspecific phylogeny was computed on 120,063 residues for each of the 45 species.

For intraspecific phylogeny of the 37 *C. cassiicola* isolates, the core set of 651 protein sequences previously used for the interspecific phylogeny was insufficient to discriminate closely related isolates (ie. the chimeric sequences were identical between CCAM1 and CCAM4, between IA and JQ, between CSRI1, CSRI2, and TSB1, and between CTHA1, CTHA3, and CTHA4). We thus used a larger set of 12,112 protein sequences found in exactly one copy in each of the 37 *C. cassiicola* isolates studied. After Gblocks selection, only the polymorphic residues were conserved in the alignment using SNP-sites (Page et al., 2016) yielding 37 chimeric concatenated sequences of 300,213 residues for final maximum likelihood tree calculation using RAxML with LG substitution matrix and 1,000 bootstraps. Tree rendering was carried out using FigTree v1.42 (http://tree.bio.ed.ac.uk/software/figtree/).

## Statistical Analysis

Data management and statistical analysis were carried out using R version 3.2.1 (R Core Team, 2015) (http://www.R-project.org/). PCA analyses used FactoMineR package (Lê et al., 2008). RNA-Seq expression was rendered using the heatmap.2 function of the plots package (Warnes et al., 2016). Hierarchical clustering was computed using the R base package using Euclidean distances and the Ward method. Gene ontology enrichments were computed by a chi-squared test between observed and expected abundances and *p*-values were adjusted with a 0.01 FDR under a dependent multiple hypothesis (GO abundances are not independent since they are graph-related) (Benjamini and Yekutieli, 2001).

## RNA-Seq Analysis

RNA-Seq analysis was carried out using our custom pipeline (Garavillon-Tournayre et al., 2015) running on the GenoToul computing cluster. Briefly, 125-bp pair-end shorts reads were first decontaminated by mapping against a database of potential contaminants, then finally mapped on annotated CCP transcripts using the Burrows-Wheeler alignment algorithm BWA-MEM (Li and Durbin, 2009). Owing to the very small proportion of short reads that targeted fungal RNA (0.05%), counts for most genes were very low. Transcripts with fewer than three counts in at least two of the eight biological replicates (25%) were filtered out, allowing 1,052 transcripts for differential analysis. Differential expression analysis was carried out using the DEseq2 package (Love et al., 2014) with 0.01 FDR yielding 353 genes.

## RESULTS AND DISCUSSION

CCP was selected as the reference among a collection of *C. cassiicola* isolates from various hosts and of diverse geographical origins, for its high virulence in the rubber tree (Breton et al., 2000; Déon et al., 2012a,b, 2014). CCP mycelium grown on PDA was fluffy and whitish gray when young, turning dark gray when older (Figures 1A,B). The conidia varied in size (from 10 to 145 µm) and shape (from aseptate to six septates) (Figure 1C). Such variability is a known feature of the species (Qi et al., 2011). Conidia inoculation of detached leaves from sensitive rubber clones induces the so-called fishbone necrosis with typical darkening of the veins within 48 h (Figure 1D).

Sequencing of the *C. cassiicola* (alias Corc) CCP genome was performed by JGI as part of the 1000 Fungal Genomes Project using Illumina sequencing. The assembled genome covers 44.85 Mbp in 500 contigs, with an average depth of 58x (Table 1). Among the 17,167 predicted gene models, 12,261 (71%) have at least one type of annotation, whether Interpro domain (56.8% of the genes), Gene Ontology (46%), Kyoto Encyclopedia of Genes and Genomes (18.1%), or EuKaryotic Orthologous Groups (13.6%) annotations.

The relatively compact genome has small amounts of repetitive DNA and transposable elements. The 44.85 Mbp CCP genome is larger than the average of 36.91 Mb reported for Ascomycota but close to that of the Dothideomycetes (44.59 Mb) (Mohanta and Bae, 2015). Furthermore, the CCP genome has 17,167 gene models, which is very high compared with the ~11,000 found on average in Ascomycota and the ~13,000 in Dothideomycetes (Mohanta and Bae, 2015). Recently, the genome of *C. cassiicola* strain UM591, isolated from the contact lens of a patient diagnosed with keratomycosis in Malaysia, was published (Looi et al., 2017). The authors used both long and short insert sequencing libraries, allowing the assembly of 41.88 Mb with 1,941 contigs. The authors identified 13,531 protein coding genes in UM591, significantly fewer than the 17,167 identified by the DOE-JGI annotation pipeline. Taken together, both analyses suggest that the *C. cassiicola* genome has undergone expansion.
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | *Corynespora cassiicola* strain CCP and the *Corynespora* Leaf Fall (CLF) disease symptoms. (A) CCP mycelium colony on PDA medium 7 days after subculture. (B) Optical microscopy view of CCP conidia in water. (C) *Hevea brasiliensis* clone PB217 leaflet inoculated with CCP spores (7 days post-inoculation, 500 spores per spot). (D) CLF symptoms on a susceptible rubber tree clone (left, RRIC103) next to a tolerant clone (right, GT1), in Nigeria.</figcaption></figure>

# PART 1. Interspecific Comparative Genomics

In the first part of this study, we present an interspecific comparative analysis of the CCP predicted proteome along with 44 species sequenced by the DOE-JGI (Table S1), representatives of diverse fungal lifestyles (i.e., saprotrophic, necrotrophic, hemibiotrophic, biotrophic, ectomycorrhizal). The selected species included six Basidiomycota and 38 Ascomycota. Special emphasis was placed on the Dothideomycetes class (23 species), which comprises many plant pathogens, including *C. cassiicola*. Altogether, 31 genera from 10 classes were represented in this study, thus maximizing the diversity of gene molecular functions.

## Interspecific Genome-Wide Phylogeny

To evaluate the dynamics of the *C. cassiicola* gene families, and more especially of those involved in virulence (i.e. effectors), we clustered the 578,823 predicted protein sequences identified in the 45 fungal genomes studied here. All-versus-all blastp analysis gave 68 million hits with *E* > 1e-5. We used this input to

carry out a Markov clustering (MCL algorithm). Since there is a trade-off between sensitivity and specificity, we experimentally tuned clustering parameters to yield the most meaningful clusters using a manually expertized gene family, the X-intrinsic proteins (XIPs), a subfamily of the large major intrinsic proteins (MIP) family (Lopez et al., 2012). Clustering parameters were tuned to obtain the fewest protein clusters while keeping XIPs distinct from other MIP families. As a result, 480,841 proteins predicted from the 45 genomes fell into 40,701 clusters, composed of a maximum of 508 genes with a mean of 11.81 genes and a median of 4 genes per cluster. The CCP genome has a high percentage of singletons (23.2%) compared with the other species (16.5% on average across all 45 species, and 10.5% across the 13 Pleosporales), the highest (37.9%) being found in *Botrytis cinerea* (Figure 2). As expected, few CCP singletons (461 only) have predicted GO functional annotations (Tables S2, S3).

We built a phylogenomic tree of the 45 fungal species based on proteins sequences having exactly one representative in each species. With our parameters, 651 clusters met with
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Interspecific phylogenomic tree, protein clusters and putative effectors of 45 fungal species. Details of the species and their genome access links are listed in Table S1. Life styles are indicated by a letter: S, saprotrophic; N, necrotrophic; H, hemibiotrophic; B, biotrophic; E, ectomycorrhizal. The maximum likelihood phylogenetic tree was based on 651 concatenated core protein sequences. Branch lengths are indicated by the bar (substitutions/site); 1,000 bootstrap values are shown as percentages. Clustered and non-clustered (species-specific) protein numbers were predicted based on OrthoMCL clustering. Putative effectors: number of gene models encoding CAZymes, peptidases, lipases and other secreted proteins (left diagram), or involved in the secondary-metabolism (NRPS, PKS, terpenes synthases, right diagram).

this criterion. These core sequences were concatenated into a
chimeric sequence used for phylogenetic distance calculation
and tree construction (Figure 2). The computed output is in
agreement with the Fungal Tree of Life (Celio et al., 2006; Schoch
et al., 2009) and places CCP in a distinct clade among the
Pleosporales. Our genome-wide phylogeny confirms with very
high support that C. cassiicola speciation has occurred before the
speciation of the Pleosporinae suborder, which includes several
families of plant pathogens such as the Phaeosphaeriaceae, the
Laeplosphaeriaceae and the Pleosporaceae.

CCP Putative Effectors

In this study, effectors are broadly defined as fungal proteins
that can modulate plant immunity. Two main classes of effectors
were distinguished (Figure 2). One includes enzymes interacting
with plant substrates (CAZymes, peptidases and lipases),
together with secreted protein of still unknown function, with

a distinction between SSPs (“small secreted proteins” of 300
amino acids or fewer) and LSPs (“large secreted proteins”).
Only classically-predicted secreted proteins (i.e., proteins with a
signal peptide but no transmembrane domain) were considered
in this study, but we cannot rule out the possibility that non-
canonical proteins may also act as effectors. The second class of
effectors corresponds to enzymes associated with the secondary
metabolites. To allow comparative genomics analysis, we
predicted the putative effectors simultaneously for C. cassiicola
and the other 44 species. In terms of quantity, C. cassiicola
appeared to be particularly rich in putative effectors compared
with the other species. A detailed list of CCP putative effectors is
given (Table S4) and discussed below.

CAZymes

CAZymes play a substantial role in carbon acquisition and
metabolism, and thus in the heterotrophic lifestyle of fungi (Zhao
et al., 2014). They have been grouped into six classes (Lombard et al., 2014): glycoside hydrolases (GHs), glycosyl transferases (GTs), polysaccharide lyases (PLs), carbohydrate esterases (CEs), auxiliary activities (AAs), and carbohydrate-binding modules (CBM). In pathogens, GTs are said to mainly ensure carbohydrate assembly, while CEs, GHs, and PLs are cell wall degrading enzymes (CWDEs).

In CCP, we found a total of 1,152 proteins sharing homology with CAZymes (Figure 2, Table S4). Among these, 326 (29%) have a secretion signal peptide, including 47 SSPs and 279 LSPs. Among the 45 species of our study, CCP ranked second richest in CAZymes after the hemibiotroph *Colletotrichum gloesporioides* (Glci). CCP predicted CAZyme comprised 417 putative GHs, 305 GTs, 209 AAs, 148 CEs, 42 CBMs, and 31 PLs. Of the predominant family, the GHs, 40% belonged to five classes: GH92 (52), GH76 (40), GH3 (28), GH43 (27), GH18 (23). GH92s are exo-acting α-mannosidases and GH76s are endo-acting α-mannanases. Mannans being among the most abundant components of hemicellulose, the mannan endo- and exohydrolysing enzymes could together play a role in softening host tissues (Moreira and Filho, 2008; Rodríguez-Gacio Mdel et al., 2012).

The most represented CAZyme modules were GT34 (121), CE10 (91), AA7 (64), AA3 (56), GH92 (52), GT2 (44), and AA9 (37). GT34s ensure hemicellulose degradation. CE10s may be involved in fungal pathogenesis (Islam et al., 2012; Aragona et al., 2014; Zhao et al., 2014; Kuan et al., 2015), although most of them have been described as acting on non-carbohydrate substrates (Cantarel et al., 2009). AA7s and AA3s have catalytic modules involved in enhancing plant cell wall degradation (Levasseur et al., 2013); GH92s are mannosidases; GT2s ensure chitin synthesis; AA9s are believed to act directly on cellulose, working in concert with AA3s for oxidative cleavage, but some of them also degrade chitin and starch (Quinlan et al., 2011; Levasseur et al., 2013). In addition, Corc has 58 pectin degrading enzymes (14 PL1s, 14 PL3s, 17 GH28s, and 9 GH78s) that could also participate in host cell wall degradation. Finally, in an extensive survey of CAZyme contents conducted on 103 diverse fungi (Zhao et al., 2014), the authors noted the expansion of PL1 and PL3 in necrotrophs.

The high CAZyme content of CCP is consistent with its life
style as a necrotrophic pathogen.

## Lipases

Lipases form an emerging field of research into virulence
mechanisms in fungi. Lipases (triacylglycerol acylhydrolases E.C.
3.1.1.3) hydrolyse carboxyl ester bonds on triacylglycerols, and
belong to the class of α/β hydrolase fold enzymes. Notably,
secreted lipases may act against host morpho-anatomical barriers
(wax and cuticle). We identified 104 putative lipases in CCP,
including 40 with secretion signals (27 SSPs and 13 LSPs)
(Figure 2, Table S4). Compared with the other species, with a
mean value of 127 lipases, CCP had one of the lowest numbers
of lipase representatives, the maximum being observed for Glci
with 271 lipases.

Cutinases (E.C. 3.1.1.74) are esterases that hydrolyse cutin
polymers into cutin monomers. They can be annotated as
members of the CE5 class of CAZymes or as lipases, since they

also belong to the α/β hydrolase fold class of enzymes, and are
reported as the abH36 superfamily in the Lipase Engineering
Database (http://www.led.uni-stuttgart.de). As leaf cuticles are
partly composed of cutin, cutinases are thought to participate
in pathogenesis through host recognition, spore adhesion,
enhancing mycelium, growth, appressorium differentiation,
pathogen penetration, and carbon acquisition (van Kan et al.,
1997; Comménil et al., 1998; Davies et al., 2000; Reis et al., 2005;
Skamnioti and Gurr, 2009; Liu et al., 2016). Globally, necrotrophs
have more cutinases than hemibiotrophs while saprotrophs and
ectomycorrhizal fungi tend to have few to none. Our analysis
again ranks CCP among the best equipped species, with 11
putative cutinases (Table S4).

## Peptidases

Peptidases, also known as proteinases, proteases and proteolytic
enzymes, are responsible for protein degradation (proteolysis)
and also belong to the α/β hydrolase superfamily. In fungal
pathogens, secreted peptidases could serve for both nutrition and
the degradation of host defense proteins, and could thus interfere
with plant innate immunity (Carlile et al., 2000; Olivieri et al.,
2002; Rodríguez-Herva et al., 2012). The role of peptidases in
pathogen virulence is supported by the upregulation of peptidase
inhibitors in *Fusarium* head blight resistant cultivars (Gottwald
et al., 2012). Peptidases are considered potential targets for
disease control using peptidase inhibitors (Lowe et al., 2015).
We identified 630 peptidases in the CCP genome, of which 124
harbored a secretion signal (Figure 2, Table S4). Here again,
CCP ranked as one of the richest for this class of effectors
after *F. oxysporum* (Fuox), *N. haematococca* (Neha) and the two
*Colletotrichum* species (Glci and Glac). Overall, serine-, metallo-
and cysteine-peptidase families dominate CCP peptidase content,
with respectively 328, 173, and 88 representatives. For Ohm et al.
(2012), Dothideomycetes have more serine-peptidases S10 and
metallo-peptidases M14, but fewer aspartic peptidases (A01) than
other plant pathogenic necrotrophs. In the case of CCP, these
families were close to the average over all necrotrophs (12 vs. 11
for S10, 19 vs. 21 for A01, and 7 vs. 6 for M14).

## Other Secreted Proteins

The CCP genome encodes 1,411 proteins with a predicted
secretion signal, which represents 8% of the total proteome. Of
these secreted proteins, only 35% fall into one of the effector
categories described above (326 CAZymes, 124 peptidases and
40 lipases, Figure 2). The others comprise 411 LSPs and 510
SSPs: 50% of all LSPs and 85% of all SSPs (Table S4). We carried
out a GO enrichment analysis on the "other secreted proteins"
subcategory (Table S5) and found an over representation
of hydrolase activities (GO:0016787, GO:0016788), related
to carbohydrates (GO:0005975), and also lipid metabolism
(GO:0006629). Six genes annotated with serine-type peptidase
activity (GO:0008236) contribute to the overrepresentation of
the term (five of the subfamily S41 and one S26). However,
these genes were not included in the peptidase category, since
they did not meet the threshold E-value in our analysis. Redox
processes largely contribute to the overrepresented GO terms
of the "other secreted proteins" category. Among them, seven
oxidases with FAD binding (GO:0050660) annotation were
found including two members with berberine/berberine-like domains thought to play a role in ROS production. The electron transport biological process (GO:0005506) was also found significantly overrepresented, notably due to the contribution of seven heme-binding proteins (GO:0020037) among which four cytochrome P450, one cytochrome B5, one peroxidase and one catalase.

As expected, very few “other SSPs” had annotations compared with “other LSPs.” In *C. cassiicola*, the only SSP characterized to date as an effector is the cassiicolin toxin (Barthe et al., 2007; de Lamotte et al., 2007) encoded by *Cas1* in CCP (Déon et al., 2012a). This special case will be addressed in Part4.

## Secondary Metabolism Genes

Genes involved in the secondary metabolism play an important role in fungal pathogenic lifestyle (Ohm et al., 2012). Polyketide synthases (PKSs), non-ribosomal peptide synthetases (NRPSs) and terpene synthases are involved in the production of polyketide, peptides and terpenes respectively. Our predictions revealed that the CCP genome encodes 20 NRPS, 32 type-1 polyketide synthases (T1PKS), 1 type-3 PKS (T3PKS), 2 hybrid NRPS-PKS proteins and 8 terpenes synthases (Figure 2, Table S4). With a total of 63 genes related to the secondary metabolism, against 34 on average, CCP ranked fourth among all 45 species, second among Dothideomycetes, and first among Pleosporales. Several of these genes were found physically clustered in the genome (Table S4).

## Comparison of CCP and UM591 Putative Effectors

We compared our results with those recently published by Looi et al describing the *in silico* identification of putative effectors in the human-hosted *C. cassiicola* isolate UM591 (Looi et al., 2017). Both studies (UM591 vs. CCP) yielded similar quantities of predicted lipases (105 vs. 104), cutinases (9 vs. 11) and secondary metabolism genes (50 vs. 63), but differed in peptidases (216 vs. 630) and CAZymes (973 vs. 1152). Notably, the most strongly represented CAZyme module in CCP (GT34, 121 members) was poorly represented in UM591 with seven members only. Conversely, the xylan esterases CE1s were among the most represented in UM591 (49) while CCP had only five. For the other CAZymes, counts were in similar ranges in both isolates. We are nonetheless confident in our counts because our annotation pipeline yielded results similar to those described for the other plant-hosted fungi studied here, as described in Part 1. Furthermore, numerous effector subfamily counts were similar between UM591 and CCP. Although annotation differences between the two studies cannot be ruled out, the hypothesis that UM591 gene loss for some but not all CAZymes families and peptidases could reflect a specialization toward non-plant hosts deserves further examination.

## Interspecific Comparison Based on All Putative Effectors

We carried out a principal component analysis (PCA), based on the whole repertoire of putative effectors, for the 45 plant-associated species of our study (Figure 3). This analysis identified

four main clusters, with no clear consistency for phylogeny or lifestyle (Figure 2 and Table S1).

CCP is the only Pleosporales member (1/24) in a cluster grouping the only Botryosphaeriales representative, *Botryosphaeria dothidea* (Bodo), and four Sordarimycetes (4/9), namely the two *Colletotrichum* species *C. gloesporioides* (Glci) and *C. acutatum* (Glac), *Nectria haematococca* (Nema), and *Fusarium oxysporum* (Fuox) (Figure 3, blue cluster). These six species are described as necrotrophic or hemibiotrophic. They are characterized as having a high number of total genes and the highest number of putative effector genes in the CAZyme, peptidase, lipase and other secreted protein categories (Figure 2).

Most (10/13) of the other Pleosporales members form a cluster (green) with seven other Ascomycetes belonging to three classes (Sordariomycetes, Eurotiomycete and Dothideomycetes). Most lifestyles are represented in this cluster: *Trichoderma asperellum* (Tras) and *T. harzianum* (Trha) are mycoparasitic; *Aspergilus nidulans* (Aspn), *Hysterium pulicare* (Hypu) and *Rhytidhysteron rufulum* (Rhru) are saprophytic, and all the others are pathogenic (biotrophic, necrotrophic or hemibiotrophic).

The last two Pleosporales species, *Leptosphaeria.maculans* (Lemu) and *Alternaria brassicicola* (Altb), stand in a phylogenetically highly diverse cluster (red), including both Basidiomycetes (5/6) and Ascomycetes belonging to five classes. Here again, most lifestyles are represented, with both pathogenic and non-pathogenic species.

The last cluster (black) is composed of the four species with the lowest numbers of total and effector genes in this selection (Figure 2), i.e., the Ascomycetes *Saccharomyces pombe* (Scpo), *Saccharomyces cerevisiae* (Sace), and *Tuber melanosporum* (Tume), together with the Basidiomycetes *Ustilago maydis* (Usma). While the three Ascomycetes are non-pathogenic (saprotrophic or ectomycorrhizal), *U. maydis* is a plant pathogen of maize and teosinte, emphasizing that even species with a small genome and a small effector arsenal can be pathogenic, though with a narrow host range (Kämper et al., 2006).

Of the 15 components of this PCA analysis, the first two dimensions covered over 62% of the total variance of the dataset. The first dimension, which accounts for 52% of the total variance, shows a large contribution from CAZymes and a smaller one from lipases and peptidases, while the second dimension shows the contribution of secondary metabolism genes. When PCA was based on CAZymes only (Figure S1), CCP was grouped with the six pathogenic species of the Sordariomycetes order, but remained separated from the other Pleosporales; *B. dothidea*, previously grouped with CPP, formed a single-species cluster. When PCA was based on the content of genes associated with the secondary metabolism (Figure S2), CCP was found closely associated with *P. teres* f. sp. *teres* (Pytt), another necrotrophic species belonging to the Pleosporales, although the overall clustering could be related neither to phylogeny nor life style.

In conclusion, *C. cassiicola* isolate CCP was consistently associated with *C. acutatum* (Glac) and *F. oxysporum* (Fuox) whatever the effector categories, and to a lesser extent with *C. gloesporioides* (Glci), *N. haematococca* (Neha) and *B. dothidea* (Bodo). All these genomes have in common expansions in most categories of effectors. This clustering agrees weakly with
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Interspecific principal component analysis (PCA) of 45 fungal species based on their putative effector composition. Species are represented by their alias name (Table S1). PCA was computed from the counts in each effector category as represented in **Figure 2** (14 dimensions). Different colors represent different clusters.</figcaption></figure>

phylogeny but seems more relevant to lifestyle considerations such as pathogenicity and nutrient uptake mechanisms: all six species are known to have broad host ranges and potentially multiple trophic modes.

*Colletotrichum* is a widespread genus of some 600 species causing necrotic spots and rot in a wide range of plants. It is a hemibiotrophic fungus able to operate lifestyle transitions during its interaction with the host plant, switching from epiphytic growth to biotrophy and finally to necrotrophy. The duration of the symptomless phase may vary significantly between species. Some *Colletotrichum* species have been clearly described as non-pathogenic plant endophytes conferring fitness benefits to their host, while others have been described as saprophytes ensuring survival on soil debris (Cannon et al., 2012; O'Connell et al., 2012; Hiruma et al., 2016). *C. gloesporioides* and *C. acutatum* species have been described as pathogens of the rubber tree (Brown and Soepena, 1994), although species of the *C. gloesporioides* complex have also been isolated from asymptomatic rubber tree (Gazis et al., 2011). Interestingly, *C. cassiicola* and *Colletotrichum* spp. are frequently co-isolated from diseased rubber leaves (unpublished result) indicating that both are adapted to the same host and cellular environment. Their proximity in terms of putative effector contents is consistent with this observation.

*F. oxysporum* and *N. haematococca* (teleomorph of *F. solani*) are also cosmopolite species able to colonize a wide variety of environments, as saprobes or necrotrophic pathogens (Fravel et al., 2003; Coleman et al., 2009). *B. dothidea* is also a widespread latent pathogen, virtually present as an endophyte in all examined woody plants, and switching to pathogenicity under stress conditions (Slippers and Wingfield, 2007). However, none of these species has yet been found in the rubber tree, although species of the same genus have occasionally been described (Gazis and Chaverri, 2010). Most fungi species so far associated with rubber were isolated after surface-sterilization of the explant. A number of species with an epiphytic phase may thus have been dismissed. A recent study emphasized the role of epiphytic growth prior to invasion in the pathogenic fungal species *Zymoseptoria tritici* (Fones et al., 2017), and the bias it may create for interpreting the life style and infection strategies of the pathogen.

*Corynespora cassiicola* fits with this model of polymorph species with a broad host range. The CCP isolate is known to be necrotrophic and highly virulent in the rubber tree, but the conditions required for disease onset remain unclear and symptom intensities vary according to the cultivar. The fungus may survive as an endophyte or latent pathogen in living tissues,
or as a saprobe on decaying material, and may potentially switch to necrotrophy under favorable conditions. Very numerous and diverse hosts of *C. cassiicola* have been described, including some outside the plant kingdom, with animal- and human-associated isolates. Nevertheless, individual *C. cassiicola* isolates display some host specialization, being restricted to a limited host range (Barthe et al., 2007; Dixon et al., 2009). Large effector repertoires may confer a greater chance of manipulating the defenses and escaping immunity in a wider range of hosts. By contrast, pathogens adapted to a narrow range of hosts maintain only a reduced set of effectors (O'Connell et al., 2012).

## PART 2. Transcriptome Analysis of CCP Compatible Interaction With Rubber Tree

As presented above, CCP contains 2,870 genes encoding putative effectors. To identify the candidates most likely to be involved in virulence in the rubber tree, we monitored the CCP transcriptome using RNA-Seq during the compatible interaction between CCP and the susceptible rubber tree clone PB260. Leaflets were inoculated with droplets of spore suspension. Germinating spores were used as reference for the differential analysis. Leaf disk samples were harvested 24 and 48 h post-inoculation. This early time frame was chosen to cover the initial stages of successful colonization: spore germination, mycelium penetration and development in the plant tissues, until the onset of the first necrosis symptoms. These symptoms are preceded by transient activation of the cassiicolin-encoding gene (Déon et al., 2012a,b).

In the germinating spores, we found that most of the 17,167 genes were expressed (Table 2, Table S6); only 2,179 transcripts had no sequencing reads coverage at all. However, the transcriptome mean coverage per gene was 457.7 while the median was only 49.5. This indicates that most CCP genes, though expressed, had low transcriptional activity in the context of spore germination outside a host.

In the inoculated leaf samples, only a tiny portion of the sequencing reads (0.05%) were mapped to the CCP transcriptome (Table 2), owing to the low fungal/plant cell ratio at this early stage of infection. Consequently, only the most

differentially expressed genes were detectable. As a result, we found that 353 genes had their expression significantly modified during the early stages of infection (Table 2 and Table S6), among which 92 (26%) were categorized as putative effectors (Figure 4). More precisely, 52 were annotated with CAZyme modules and 40 categorized as lipases/phospholipases, peptidases or other secreted proteins. No gene of the secondary metabolism was found differentially expressed in this context.

### Differentially Expressed Cazymes (Figure 4A)

The 52 differentially expressed genes with CAZyme modules included 20 GHs, 11 AAs, 10 CEs, 6 GTs, 3 CBMs, and 2 PLs (respectively 4.7, 5.2, 6.7, 1.9, 7.1, and 6.5% of the total count in each category). Interestingly, fewer GTs (only 1.9%), responsible for carbohydrate assembly, were differentially expressed compared with CEs, GHs and PLs, described as cell wall degrading enzymes (CWDEs). Most of the differentially expressed CAZyme genes (84%) were upregulated. Twenty differentially expressed CAZyme genes (40%) had a predicted secretion signal (they are marked with an asterisk). All of them were up-regulated, and notably five with a very high fold change (between 7.1 and 9.6): the two putative pectin lyases (PL1 and PL3), one glucanase, one FAD-linked oxidase and one lectin beta domain-containing protein. They can be considered as good effector candidates potentially involved in cell wall degradation (Kubicek et al., 2014; Ismail and Able, 2016). Other up-regulated candidate effectors were potential actors of cell wall degradation: a secreted cellulase, two GH16s, four non-secreted CE10 carbohydrate esterases, two AA9s cellulose-binding glycoside hydrolases, and two AA3s glucose-methanol-choline (GMC) oxidoreductases. Two upregulated CE4s with secretion signals annotated as "chitin-binding, glycoside hydrolase-deacetylase" could play a role in the camouflage of chitin oligomers that would normally be acting as elicitors in the host plant, as previously hypothesized (Tsigos and Bouriotis, 1995). A secreted alpha-amylase (185,576), along with 417,524, a secreted GH31 with predicted glycoside hydrolase activity, could ensure carbon acquisition. Finally, other CAZyme activities such as glycosylation (GTs) could be associated with post-translational modification of proteins delivered to the host that could modulate pathogen virulence (Kubicek et al., 2014). Notably, in CCP, the effector cassiicolin was shown to be glycosylated (Barthe et al., 2007).

None of the eight repressed CAZYme candidates had a predicted secretion signal. We can however underline the repression of a chitin synthase gene, potentially as part of a strategy to avoid chitin-induced elicitation of the plant defenses.

### Differentially Expressed Lipases, Peptidases and Other Secreted Proteins (Figure 4B)

Besides CAZymes, we found that 11 peptidases, seven lipases and 22 other secreted proteins were differentially expressed, representing respectively 1.7, 6.7, and 2.4% of the total count in each effector category. They were mostly up-regulated.

<table><caption>TABLE 2 | Metrics of the CCP gene transcripts.</caption><tbody><tr><td>Number of transcripts expressed in germinating spores</td><td>14,988</td></tr><tr><td>Differentially expressed in inoculated rubber tree leaves</td><td>353</td></tr><tr><td>Putative effectors</td><td>92</td></tr><tr><td>Others</td><td>261</td></tr><tr><td>Non differentially expressed in inoculated rubber tree leaves</td><td>14,635</td></tr><tr><td>Putative effectors</td><td>2,397</td></tr><tr><td>Others</td><td>12,238</td></tr><tr><td>Number of transcripts non expressed in germinating spores</td><td>2,179</td></tr><tr><td>Differentially expressed in inoculated rubber tree leaves</td><td>0</td></tr><tr><td>Non differentially expressed in inoculated rubber tree leaves</td><td>2,179</td></tr><tr><td>Putative effectors</td><td>381</td></tr><tr><td>Others</td><td>1,798</td></tr><tr><td>Total</td><td>17,169</td></tr></tbody></table>
A CAZymes (52)

<figure><img src="image_6.png" /></figure>

* 185576 (GH13/CBM20)

* 584290 (AA5)

* 523708 (GH76)

* 597988 (GH16)

* 577365 (AA6)

* 568081 (CBM50)

* 627480 (GH76)

* 383757 (CE4)

* 577841 (CE10)

* 491715 (CBM13)

* 478760 (GH92)

* 620773 (CE10)

* 568400 (GH92)

* 257707 (AA3)

* 417524 (GH31)

* 636352 (GH109)

* 604042 (GH65)

* 629841 (GH5)

* 574023 (AA4/CE10)

* 576856 (GH3)

* 574440 (CE4)

* 577285 (GH17)

* 456712 (CE4)

* 631521 (GT35)

* 497436 (AA9)

* 567012 (GH16)

* 554214 (CE10)

* 81004 (AA3)

* 23468 (GH72)

* 679644 (GH18)

* 557752 (GT41)

* 376070 (GT2)

* 308190 (GT50)

* 585576 (AA7)

* 563211 (GH47)

* 569483 (AA4/AA7)

* 517619 (CE10)

* 568685 (PL3)

* 577973 (PL1)

* 553845 (GH131)

* 473551 (AA9)

* 680806 (AA7)

* 599951 (AA1)

* 629997 (GH18)

* 569089 (GT2)

* 600578 (CE11)

* 70313 (GT4)

* 577402 (CE11)

* 569191 (GH18)

* 503688 (GH76)

* 618747 (CBM13)

* 637055 (AA10)

B Lipases - Proteases and other secreted proteins (40)

<figure><img src="image_7.png" /></figure>

* 19745 Cupredoxin (SSP)

* 573335 Peptidase

* 243454 Lactohydrolase (LSP)

* 601858 Necrosis inducing protein NPP1 (SSP)

* 513290 Alkaline phosphatase (LSP)

* 578510 Alkaline phosphatase (SSP)

* 508575 SSP

* 492140 Phospholipase A2 (SSP)

* 673479 LSP

* 571426 Lipase

* 212468 Peptidase

* 485487 Cutinase

* 577878 Cell surface antigen (LSP)

* 605944 (LSP)

* 574375 Tyrosinase/Di-copper centre-containing (LSP)

* 573880 LSP

* 491861 Lipase

* 518433 LSP

* 514165 LSP

* 94041 Cupredoxin (SSP)

* 84402 Mannose-binding lectin (MBL) / DNase I-like (LSP)

* 149614 Peptidase

* 24345 SSP

* 576232 Peptidase

* 646750 SSP

* 308561 Tyrosinase/Di-copper centre-containing (LSP)

* 555383 Lipase (LSP)

* 888888 Cassiicolin (Cystein rich SSP)

* 495460 Peptidase

* 498213 Catalase (LSP)

* 9411 Cupin_5 (SSP)

* 632054 Lipase

* 551674 Peptidase

* 545393 Peptidase

* 366051 Peptidase

* 640080 Peptidase

* 568026 Lipase

* 576717 SSP

* 514786 Peptidase

* 248763 Peptidase (LSP)

Log2 fold change

<figure><img src="image_8.png" /></figure>

**FIGURE 4** | RNA-Seq heatmaps of CCP putative effectors differentially expressed 24 and 48 h after spore inoculation on detached rubber tree leaves (susceptible clone PB260). Differential expression was calculated as the Log2 fold change at each time point using the spore suspension as reference. The code numbers refer to JGI transcript IDs (http://genome.jgi.doe.gov/Corca1) except for cassiicolin Cas1 (888888) which was manually annotated. Genes with predicted secretion signals are noted with an asterisk. Colored bars on the left side define the three clusters obtained by hierarchical classification of the differential expression values. (A) CAZymes; (B) other effectors.

Of the 11 peptidases, nine were found to be upregulated, but none had a secretion signal, in contrast with previous studies on *F. graminearum* showing an enrichment in secreted peptidases among the upregulated gene clusters (Lowe et al., 2015). Two peptidases were down-regulated, including one potentially secreted.

Among the lipases, one was annotated as cutinase and another as phospholipase. Both genes were strongly upregulated, but only phospholipase was predicted to be secreted. The other lipases had no secretion signal, except for one (555,383).

The other secreted proteins comprised 14 LSPs and 11 SSPs. All were upregulated except for one down-regulated SSP. One of the upregulated SSPs is a necrosis-inducing peptide 1 homolog (NPP1) reported to be an elicitor of plant defenses (Fellbrich et al., 2002). This peptide is required for virulence in barley (Fiegen and Knogge, 2002). We also found a cupin 5 SSP, reportedly involved in oxidative stress reactions (Dunwell, 1998; Cechin et al., 2008). Two upregulated cupredoxins (SSPs, 19,745, and 94,041) could participate in copper stabilization in concert with tyrosinase

(574,375 and 308,561) for pigment production and/or with the multi copper-dependent laccase (599,951, CAZyme AA1) for lignin degradation (Benoit et al., 2015). A lactonehydrolase LSP (243,454) was found highly upregulated. Lactonehydrolase was shown to detoxify the zearalenone mycotoxin (Takahashi-Ando et al., 2002). Brown et al (Brown et al., 2012) found two detoxifying lactonehydrolases in *F. graminearum* and hypothesized a self-defense mechanism. Finally, the cysteine-rich SSP cassiicolin, only effector characterized to date in *C. cassiicola* (Déon et al., 2012a,b) was also found among the upregulated SSPs, as expected.

Other Differentially Expressed Genes

Beside the 92 differentially expressed candidate effectors (Figure 4), 261 other genes were significantly modulated in our RNA-Seq experiment (Table 2, Table S6), evenly distributed between upregulated and downregulated genes (respectively 129 and 132 genes), in contrast with candidate effectors that were mostly upregulated during the interaction with the host.
We carried out GO term enrichment analysis on this set of genes (Table S7). In the downregulated gene cluster, terms associated with transcriptional and translational activities were overrepresented. In the upregulated gene cluster, terms associated with membrane transport (notably of carbohydrates) and oxidoreductase activity were over-represented. We can speculate that the decrease in both transcription and translation of these fungal genes is part of the response to plant defense mechanisms, while the upregulation of transporter activities is in line with increased nutrient acquisition by the pathogen.

Among the genes identified as singletons according to our interspecific protein clustering (Figure 2, Table S2) eleven were differentially expressed (i.e., 93,300, 567,860, 646,750, 119,471, 579,535, 594,582, 641,745, 673,479, 372,273, 349,923, and the cassiicolin gene *Cas1*). However, no annotation was currently available for these, besides cassiicolin and two other genes annotated as encoding secreted proteins (the LSP 673479 and the SSP 646750). Genes of still unknown function, specific to our pathogen and differentially expressed during its compatible interaction with the rubber tree, are clearly of interest and warrant further study.

This transcriptomic analysis has allowed the identification of differentially expressed genes, among which 45 putative secreted candidate effectors may play a key role in this compatible interaction with the rubber tree. This now needs validation by functional analysis.

## PART 3. Intraspecific Comparative Genomics

In the first and second part of this study we identified all putative effectors in the CCP genome and experimentally narrowed the list to those modulated during the compatible interaction with the rubber tree. The third objective of our study was to identify the effectors that could explain specific features in terms of pathogenicity and virulence profiles among the highly diverse *C. cassiicola* isolates. We thus sequenced the genomes of 35 plant-associated *C. cassiicola* isolates (Table 3), using Illumina technology with a theoretical average depth of ~45x, to conduct an intra-species comparison of all putative effectors. The 35 isolates were mostly collected from *Hevea brasiliensis* (74%), but also from other host plants: *Cucumis sativus* (ATI11, CBS129.5, EDIG, IA, JQ, PB), *Glycine max* (777AA, RUD), *Piper hispidinervum* (LP07) and *Vernonia cinerea* (GSO2). All were collected from lesioned leaves, except for three (E55, E79, and E139) which originate from asymptomatic rubber tree leaves from Brazil, where CLF is absent or very limited. We selected members belonging to various phylogenetic clades according to a prior phylogeny based on four concatenated *loci* sequences (Déon et al., 2014). Among each clade, whenever available, we selected isolates with and without cassiicolin genes, with diversified geographical origins. Isolate CCP was used as the reference.

Our initial strategy was to map the 767 million WGS short reads from the 35 isolates onto the CCP genome assembly. This mapping strategy can reveal structural variations of candidate effectors, such as sequence polymorphism, presence/absence

variants (PAVs) and copy number variants (CNVs). However, the mapping rates were found to be too low for most isolates (Table 3), ranging from 65.4 to 93.7% of the reads (73.3% on average). In addition, this strategy failed to reveal the cassiicolin locus in 15 isolates of our selection from which this gene was previously amplified by PCR (Déon et al., 2014). We thus moved on to *de novo* assembly. The proportion of sequence reads used in the assemblies of the 35 genomes (min 94.2%, max 96.2%, average 95.7%) was significantly higher than the proportion of reads mapped onto CCP. The assembled genome sizes are very similar to CCP, ranging from 40 to 44.7 Mb. No correlation was observed between assembled genome size and depth, and so we considered that the sequencing effort was sufficient. We then searched for orthologs of all predicted CCP genes in the different assemblies. We found values ranging from 15,160 to 23,231 genes where CCP had 17,167 genes. No clear relation could be found between genome size or gene number and features such as host or geographical origin.

### Intraspecific Genome-Wide Phylogeny

We built an intraspecific genome-wide phylogeny of 37 *C. cassiicola* isolates (i.e. CCP, the 35 resequenced plant-associated isolates and the human-hosted isolate UM591). The 651-sequence core set previously used for the intraspecific phylogeny was found insufficient to discriminate closely related *C. cassiicola* isolates. We thus used a new core set of 12,112 protein sequences (70% of CCP predicted proteome) found in exactly one copy in each isolate. Overall, this genome-wide intraspecific phylogeny (Figure 5) was in close agreement with the one previously built using four *loci* only (Déon et al., 2014). All previously identified clades and sub-clades were confirmed. Moreover, eight isolates from Clade A (ATI11, EDIG, RUD, CBS129.5, IA, JQ, PB, and 777AA) now formed a new highly supported subclade that was named A5. These were Brazilian isolates from cucumber and soybean, while all other isolates in Clade A originated from the rubber tree (except for GSO2, isolated from *Vernonia*). Several previously unclassified isolates could be positioned: CCAM1 and CCAM4 were placed in Clade D with E139 and E79; LP07 joined E55 in Clade G and GSO2, while slightly apart, was included in Subclade A4. This robust genome-wide phylogeny shows that an early differentiation event occurred between Clade B and the other clades.

UM591 was found to be distantly related to plant isolates of Clade B4. Whether UM591 belongs to a distinct Subclade in B or a new clade needs further investigation. Although it shares a common ancestor with plant isolates, there is no current evidence that UM591 is a plant pathogen itself, or conversely that plant isolates could also infect humans.

### C. cassiicola Core and Accessory Effectors

We searched for homologs of the 2,870 CCP putative effectors in the 35 re-sequenced genomes of the plant-associated *C. cassiicola* isolates. Among them, 2,258 genes (78.6%) were found in exactly one copy in each isolate and were considered as putative "core effectors" for this set of isolates. The remaining 612, either absent or present in multiple copies in at least one isolate (e.g. *Cas*), were considered as putative "accessory effectors" (Table 4, Table
<table><caption>TABLE 3 | Resequenced corynespora cassiicola isolates and Whole Genome Shotgun (WGS) sequencing metrics.</caption><thead><tr><th>Isolate name</th><th>Geographical origin</th><th>Host plant</th><th>Sanitary status of host tissue</th><th>aToxin class</th><th>WGS mapping on CCP</th><th>Property mapped reads (%)</th><th>Assembled genome size</th><th>bContigs</th><th>cN50</th><th>Used reads</th><th>Used reads (%)</th><th>Coverage depth</th><th>Number of predicted genes</th></tr></thead><tbody><tr><td>777AA</td><td>Brazil (BR)</td><td>Glycine max (Gm)</td><td>Lesioned</td><td>Cas2</td><td>25,946</td><td>69.6</td><td>43,337,818</td><td>54,421</td><td>15,652,205</td><td>96.2</td><td>36</td><td>16,745</td></tr><tr><td>AT111</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2</td><td>18,580</td><td>69.7</td><td>42,232,695</td><td>55,794</td><td>18,402,398</td><td>96.0</td><td>34</td><td>17,583</td></tr><tr><td>OBS129,25</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2</td><td>18,580</td><td>73.4</td><td>41,411,750</td><td>56,420</td><td>21,960,627</td><td>96.0</td><td>43</td><td>16,848</td></tr><tr><td>CCAM1</td><td>Cameroon (CA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>25,726</td><td>8.4</td><td>43,770,965</td><td>44,219</td><td>15,395,570</td><td>96.2</td><td>28</td><td>20,398</td></tr><tr><td>CCAM2</td><td>Cameroon (CA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas1</td><td>21,575</td><td>93.7</td><td>43,143,070</td><td>45,751</td><td>19,761,908</td><td>95.3</td><td>36</td><td>20,926</td></tr><tr><td>CCAM3</td><td>Cameroon (CA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas1</td><td>34,046</td><td>91.5</td><td>44,706,950</td><td>44,020</td><td>25,509,053</td><td>95.2</td><td>43</td><td>23,231</td></tr><tr><td>CCAM4</td><td>Cameroon (CA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>24,864</td><td>82.1</td><td>43,735,702</td><td>50,281</td><td>23,499,794</td><td>96.0</td><td>43</td><td>20,437</td></tr><tr><td>CC113</td><td>Côte d'Ivoire (CI)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>26,299</td><td>68</td><td>43,462,677</td><td>57,107</td><td>20,289,471</td><td>95.8</td><td>35</td><td>17,139</td></tr><tr><td>CC16</td><td>Côte d'Ivoire (CI)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>25,732</td><td>68.8</td><td>43,521,463</td><td>58,564</td><td>21,920,208</td><td>95.6</td><td>38</td><td>17,106</td></tr><tr><td>CCP</td><td>Philippines (PH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas 1</td><td>244</td><td></td><td>44,850,000</td><td></td><td></td><td></td><td></td><td>17,167</td></tr><tr><td>CGAB1</td><td>Gabon (GA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>23,001</td><td>70.9</td><td>42,622,471</td><td>54,421</td><td>15,652,205</td><td>96.2</td><td>29</td><td>17,143</td></tr><tr><td>CGAB2</td><td>Gabon (GA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>25,233</td><td>70</td><td>42,753,115</td><td>56,175</td><td>15,803,351</td><td>96.1</td><td>29</td><td>17,320</td></tr><tr><td>CIND3</td><td>India (IN)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>25,446</td><td>68.4</td><td>43,142,275</td><td>55,647</td><td>18,552,642</td><td>95.9</td><td>33</td><td>17,065</td></tr><tr><td>CLN16</td><td>Malaysia (MA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>26,480</td><td>69.6</td><td>43,236,492</td><td>49,745</td><td>14,478,070</td><td>95.8</td><td>26</td><td>17,203</td></tr><tr><td>CSB16</td><td>Malaysia (MA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas5</td><td>8,172</td><td>66.9</td><td>41,090,012</td><td>58,209</td><td>17,911,142</td><td>95.2</td><td>34</td><td>15,160</td></tr><tr><td>CSR11</td><td>Sri Lanka (SL)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>16,169</td><td>66.2</td><td>41,512,866</td><td>57,706</td><td>20,366,586</td><td>95.8</td><td>39</td><td>16,143</td></tr><tr><td>CSR12</td><td>Sri Lanka (SL)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas5</td><td>22,858</td><td>65.4</td><td>42,036,764</td><td>57,382</td><td>19,938,115</td><td>95.7</td><td>37</td><td>16,247</td></tr><tr><td>CSR15</td><td>Sri Lanka (SL)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>24,753</td><td>69</td><td>43,093,641</td><td>56,383</td><td>19,074,836</td><td>95.5</td><td>34</td><td>17,128</td></tr><tr><td>CTHA1</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>16,615</td><td>84.4</td><td>42,179,887</td><td>44,014</td><td>19,483,933</td><td>95.7</td><td>38</td><td>19,069</td></tr><tr><td>CTHA2</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>18,777</td><td>67</td><td>41,852,029</td><td>54,006</td><td>19,105,681</td><td>96.0</td><td>37</td><td>16,189</td></tr><tr><td>CTHA3</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>6,035</td><td>84.2</td><td>40,558,447</td><td>33,197</td><td>12,008,130</td><td>95.5</td><td>24</td><td>15,705</td></tr><tr><td>CTHA4</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>6,794</td><td>86.3</td><td>40,998,259</td><td>43,846</td><td>25,526,040</td><td>95.2</td><td>51</td><td>16,564</td></tr><tr><td>CTHA5</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>19,267</td><td>71.8</td><td>41,979,272</td><td>65,277</td><td>15,226,899</td><td>95.7</td><td>29</td><td>17,056</td></tr><tr><td>CTHA6</td><td>Thailand (TH)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas0</td><td>24,161</td><td>70.6</td><td>42,283,119</td><td>72,894</td><td>14,104,323</td><td>94.8</td><td>27</td><td>17,145</td></tr><tr><td>E139</td><td>Brazil (BR)</td><td>Hevea brasiliensis (Hb)</td><td>Healthy</td><td>Cas4</td><td>27,903</td><td>80.2</td><td>44,431,849</td><td>45,812</td><td>23,838,352</td><td>96.2</td><td>43</td><td>21,058</td></tr><tr><td>E55</td><td>Brazil (BR)</td><td>Hevea brasiliensis (Hb)</td><td>Healthy</td><td>Cas0</td><td>17,402</td><td>74.3</td><td>41,269,705</td><td>49,888</td><td>25,968,811</td><td>94.2</td><td>48</td><td>15,772</td></tr><tr><td>E79</td><td>Brazil (BR)</td><td>Hevea brasiliensis (Hb)</td><td>Healthy</td><td>Cas4</td><td>17,402</td><td>81.2</td><td>44,028,261</td><td>47,985</td><td>22,067,764</td><td>96.2</td><td>40</td><td>21,866</td></tr><tr><td>EDIG</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2</td><td>17,402</td><td>70.5</td><td>42,249,625</td><td>37,754</td><td>12,130,622</td><td>96.2</td><td>22</td><td>17,357</td></tr><tr><td>GSO2</td><td>Brazil (BR)</td><td>Vernonia cinerea (Vo)</td><td>Lesioned</td><td>Cas2</td><td>22,836</td><td>67.4</td><td>43,198,292</td><td>30,770</td><td>12,431,866</td><td>95.7</td><td>22</td><td>17,060</td></tr><tr><td>IA</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2 + 7</td><td>23,053</td><td>68.7</td><td>42,992,747</td><td>58,869</td><td>29,543,154</td><td>96.0</td><td>53</td><td>17,059</td></tr><tr><td>JQ</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2 + 7</td><td>20,611</td><td>70.5</td><td>42,862,184</td><td>55,346</td><td>17,233,550</td><td>96.1</td><td>32</td><td>17,001</td></tr><tr><td>LPO7</td><td>Brazil (BR)</td><td>Piper hispidinervum (PH)</td><td>Lesioned</td><td>Cas0</td><td>30,661</td><td>71.8</td><td>43,904,214</td><td>49,813</td><td>25,807,090</td><td>95.3</td><td>47</td><td>19,014</td></tr><tr><td>PB</td><td>Brazil (BR)</td><td>Cucumis sativus (Cs)</td><td>Lesioned</td><td>Cas2 + 7</td><td>21,830</td><td>69.4</td><td>42,924,591</td><td>55,818</td><td>22,793,115</td><td>96.1</td><td>41</td><td>16,977</td></tr><tr><td>RUD</td><td>Brazil (BR)</td><td>Glycine max (Gm)</td><td>Lesioned</td><td>Cas2 + 6</td><td>18,250</td><td>70.8</td><td>42,698,174</td><td>49,215</td><td>15,754,738</td><td>95.8</td><td>29</td><td>16,905</td></tr><tr><td>SS1</td><td>Malaysia (MA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas5</td><td>21,219</td><td>66.9</td><td>41,929,664</td><td>52,434</td><td>17,449,048</td><td>95.7</td><td>34</td><td>16,152</td></tr><tr><td>TSB1</td><td>Malaysia (MA)</td><td>Hevea brasiliensis (Hb)</td><td>Lesioned</td><td>Cas5</td><td>21,882</td><td>66.3</td><td>42,091,408</td><td>55,474</td><td>20,226,399</td><td>95.7</td><td>38</td><td>16,238</td></tr></tbody></table>

In bold characters: reference CCP isolate.

aToxin class according to this analysis.

bNumber of contigs in each genome (scaffolds in the case of CCP).

cContig length at which 50% of total bases in the assembly are in contigs of that length or greater.
<figure><img src="image_9.png" /><figcaption>**FIGURE 5** | Intraspecific phylogenomic tree and putative accessory effectors of *C. cassiicola* isolates. The maximum likelihood phylogenomic tree was based on 12,420 conserved protein sequences. Branch lengths are indicated by the bar (substitutions/site); 1,000 bootstrap values are shown as percentages. The isolates names are followed by their country and host codes, and by the toxin class, as indicated in **Table 1**. The diagram on the right represents the composition in putative accessory effectors (i.e. either absent or varying in copy number in at least one isolate).</figcaption></figure>

S8). Accessory effectors were found in all functional categories, in variable ratios: 24 were associated with the secondary metabolism (38.1% of all putative effectors in this category), 195 were CAZymes (16.9%), 24 lipases (3.9%), 135 peptidases (22.1%) and 234 “other secreted proteins” divided into 108 LSPs (26.3%) and 126 SSPs (24.7%). The distribution was evenly balanced between secreted (307) and non-secreted (305) accessory effectors. The distribution of the putative accessory effectors in the different isolates is depicted in **Figure 5** as the sum of orthologs found in each functional category. CCP isolate, being the reference, has all 612 putative accessory effectors. Of these, only 20 were specific to CCP (1 secondary metabolism gene, 6 CAZymes, 4 peptidases, 1 lipase, 4SSPs, 3 LSPs, Table S8). Twenty-two others, shared by CCP, CCAM2 and CCAM3, were specific to Clade C isolates (2 secondary metabolism genes, 8 CAZymes, 7 peptidases, 1 lipase, 4 LSPs) (Table S8). However, no effector specific to an isolate other than CCP (i.e., not sharing significant homology with a CCP protein) could be identified, since CCP was used as the reference.

We carried out an intraspecific PCA based on the composition in accessory effectors found among the 36 *C. cassiicola* isolates (**Figure 6**). Despite the high number of dimensions, the first two

explained >32% of the total variance. The classification obtained was fully consistent with the genome-wide phylogeny depicted in **Figure 5**. We note that the phylogram and the accessory effector PCA cladogram were obtained by totally independent datasets: in the first case, the ~12,000 core genes (found in all isolates in one copy), and in the second case, the accessory effector gene subset. This result suggests that new accessory effectors have been acquired mostly through evolution. However, isolated acquisitions by horizontal gene transfer cannot be excluded, though not detected here. Whatever the classification approach, no clear structure by geographic origin or host plant could be established.

Further work could seek to relate specific sets of predicted effectors with biological features such as pathogenicity specificities on various hosts, or virulence profiles on different rubber clones. However, such biological information is still scarce. CCP and its closest relatives, CCAM2 and CCAM3 (Clade C), show the highest effectors diversity. They were previously found to be more aggressive than isolates of other types when spore-inoculated on the susceptible rubber clone PB217 (Déon et al., 2014). Their filtrates were also found to
<table><caption>TABLE 4 | Putative core and accessory effectors identified among 36 *C. cassiicola* isolates.</caption><thead><tr><th rowspan="3">Effector categories</th><th colspan="3">CCP total effectors</th><th colspan="3">Accessory effectors</th></tr><tr><th rowspan="2">Total</th><th colspan="2">Secreted</th><th rowspan="2">Total<br>Nbr (% of CCP)</th><th colspan="2">Secreted<br>Nbr (% of CCP)</th></tr><tr><th>Secreted</th><th>Not secreted</th><th>Not secreted</th><th>Nbr (% of CCP)</th></tr></thead><tbody><tr><td>Secondary metabolism</td><td>63</td><td></td><td>63</td><td>24 (38.1)</td><td></td><td>24 (38.1)</td></tr><tr><td>Cazymes</td><td>1,152</td><td>326</td><td>826</td><td>195 (16.9)</td><td>45 (13.8)</td><td>150 (18.2)</td></tr><tr><td>Lipases</td><td>104</td><td>40</td><td>64</td><td>24 (3.9)</td><td>5 (12.5)</td><td>19 (29.7)</td></tr><tr><td>Peptidases</td><td>630</td><td>124</td><td>506</td><td>135 (22.1)</td><td>23 (18.5)</td><td>112 (22.1)</td></tr><tr><td>Other LSPs</td><td>411</td><td>411</td><td></td><td>108 (26.3)</td><td>108 (26.3)</td><td></td></tr><tr><td>Other SSPs</td><td>510</td><td>510</td><td></td><td>126 (24.7)</td><td>126 (24.7)</td><td></td></tr><tr><td>Total</td><td>2,870</td><td>1,411</td><td>1,459</td><td>612 (21.3)</td><td>307 (21.8)</td><td>305 (20.9)</td></tr></tbody></table>

Accessory effectors are effectors either absent or varying in copy number in at least one isolate, as opposed to the core effectors, found in exactly one copy in each of the 36 isolates. Only conventional ER/Golgi-dependent secretion was considered; LSP (large secreted proteins, >300 amino acids), SSP (small secreted proteins, <300 amino acids).

be highly toxic when tested on detached leaves over a range of rubber clones, using a method based on electrolyte-leakage measurements (Tran et al., 2016). We advance the hypothesis that the effector expansion in Clade C isolates may favor higher virulence and/or a broader host range. Some isolates in Clade A4 also displayed high filtrate toxicity on various clones in Africa (Tran et al., 2016). In our analysis, four putative accessory effectors (626988, 580294, 580296 and 580661) were found in all A4 isolates, but not in other types besides two Clade C and two Clade A5 isolates. None was differentially expressed in CCP, but their expression in A4 isolates remains to be investigated. Isolates in Clade B are also frequently associated with CLF in Asia (Shuib et al., 2015). In our study, nine candidate effectors, all potentially secreted, were shared by Clade B and C isolates only, among which one LSP (577878) differentially expressed in CCP. Finally, three isolates (E55, E79 and E139) categorized as rubber tree endophytes (Déon et al., 2012b) had effector repertoires equivalent to or even larger than that of pathogenic isolates of the same clades. These endophytic isolates may be adapted to a broader host range thanks to their large effector repertoires. Whether they may be necrotrophs on some hosts and which effector(s) might be involved is unknown.

Within the 353 differentially transcribed CCP genes revealed by the RNA-Seq analysis (Table S6), 19 were categorized as potential accessory effectors (Table S8). Fourteen were upregulated and five downregulated. Nine are CAZymes (597,988; 577,841; 574,440; 577,285; 568,400; 569,191; 577,402; 618,747; 600,578), four are peptidases (149614; 495460; 366051; 514786), three are “other LSPs” (513,290; 673,479; 577,878) and three are up-regulated SSPs, among which cassiicolin, shared by 18 isolates, and two uncharacterized SSPs (24,345 and 646,750), specific to CCP.

The interspecific clustering revealed 4,047 singletons in the CCP genome (i.e. without homology in the other species) (Table S2), of which 158 were putative effectors. Our intraspecific analysis revealed that 29 of them were accessory, including two associated with the secondary metabolism, five CAZymes, five lipases, one peptidase, nine SSPs and seven LSPs (Table S8). Interestingly, the SSP 646750 discussed above is the only

putative accessory effector differentially expressed (upregulated) in the compatible interaction that has no ortholog at either the interspecific or intraspecific scale. However, its specificity to the CCP isolate could be artefactual, due to the lack of assembly of the genome region carrying this gene.

## PART 4: The Case of Cassiicolin

The most emblematic and best characterized *C. cassiicola* effector to date is the cassiicolin toxin. The *Cas1* gene (GenBank EF667973) was not predicted by the JGI annotation pipeline despite high RNA-Seq reading coverage (>3,900), possibly due to its short size and its intron-exon organization. We could nevertheless locate it on scaffold 130:2337-2848, minus strand.

<figure><img src="image_10.png" /><figcaption>FIGURE 6 | Intraspecific principal component analysis (PCA) of 36 plant-associated *C. cassiicola* isolates based on their composition in putative accessory effectors (i.e. either absent or varying in copy number in at least one isolate). PCA was computed from the counts in each effector subcategory (611 dimensions). Different colors represent different clusters.</figcaption></figure>
According to the JGI VISTA plot of scaffold 130 (4,131 bp), compared with the other genomes of our interspecific study, no nucleotide sequence homology could be found by blast against the entire NCBI nr and MycoCosm databases. This suggests that the genomic region bearing *Cas1* was acquired recently, either through evolution or horizontal gene transfer from still unknown organisms. Interestingly, a fragment homologous to a MOLLY transposon from *S. nodurum* (GeneBank AJ488502) was located at the end of scaffold 130, only ~2 kb downstream from *Ca1*. Another gene (580614) encoding a predicted protein with a SAM-methyltransferase conserved domain was identified on the same scaffold, between *Cas1* and the MOLLY transposon sequence. Proximity with the transposon suggests that *Cas1* and neighboring genes may have been multiplied through transposition events, although draft genome assembly alone could not support this hypothesis. The very high RNA-Seq reads coverage of *Cas1* (>3,900) supports this hypothesis. The neighboring SAM-methyltransferase gene may be involved in the regulation of cassiicolin activity through methylation: the mature cassiicolin was shown to carry a methylated mannose on the second residue (Barthe et al., 2007), but whether this post-transcriptional modification is required for virulence is unknown.

In this study, we could confirm that the cassiicolin gene is transcriptionally up-regulated during the early phase of the compatible interaction with the rubber tree (Figure 4). PCR screening had previously evidenced that it was an accessory effector, since it was detected in only 47% of the tested *C. cassiicola* isolates, with various isoforms (*Cas1* to *Cas6*) used to delineate toxin classes (Déon et al., 2014). Our intraspecific comparative genomic analysis confirmed the previously defined toxin classes for most isolates except IA, JQ and PB, three phylogenetically close relatives isolated from cucumber (Figure 5). A *Cas2* gene was detected in isolate JQ, previously categorized as *Cas0*. Additionally, we were able to identify a new cassiicolin isoform, shared by IA, JQ and PB, that was named *Cas7*. To rule out any artifact due to the assembly process, we confirmed the sequence by PCR amplification and Sanger sequencing of that region. The *Cas7* sequence was registered in GenBank (BankIt2035475 Seq1 MF564202). It is quite divergent from the previously identified cassicolin genes, with only 72.3–76.7% nucleotide sequence identity and 71.9–77.8% deduced amino acid sequence identity. However, we were not able to amplify *Cas7* transcripts from germinating spores of the *Cas2*+7 isolates IA, JQ and PB. RT-qPCR analysis of *Cas2*+7 isolates during interactions with plant hosts may demonstrate the functionality of *Cas7* and confirm its protein sequence.

Cassicolin is the perfect illustration of a species-specific effector responsible for intraspecific diversity due to sequence polymorphism, presence/absence and copy number variations. However, the phenotypic traits associated with these variations are still unclear. Our transcriptomic analysis revealed the existence of other candidate effectors with similar features, i.e. secreted proteins differentially expressed during the compatible interaction between *C. cassicola* isolate CCP and the rubber tree, with similar or higher orders of fold change magnitude. Further

experimentation will thus be needed to clarify the role of each candidate in CLF.

## GENERAL CONCLUSION

This study provides the genomic description of a plant pathogenic *C. cassiicola* isolate and the first transcriptomic analysis for this species. By combining *in silico* mining of all putative effectors and experimental identification of the fungal genes differentially expressed during the compatible interaction with the rubber tree, we were able to identify pertinent candidate effectors potentially involved in CLF disease, in addition to the already characterized effector cassiicolin. Resequencing and *de novo* assembly of a large set of *C. cassiicola* isolates allowed very accurate refining of the intraspecific phylogeny and the genome-wide description of their respective effector repertoires. In addition to core effectors shared by all isolates, accessory effector repertoires were generally consistent with the phylogeny, besides a few isolate-specific variations. At this stage, it is difficult to offer a cogent interpretation of such specificities owing to the lack of phenotypic information in terms of virulence profiles or host ranges. Further transcriptomic, proteomic or metabolomic analyses may help identify the most pertinent effectors. These could then be used to develop phenotyping tests or genotyping tools, after identification of their plant molecular targets, for the selection of tolerant cultivars.

## AUTHOR CONTRIBUTIONS

VP-R conceived and managed the research, with scientific support from FM. DL performed data mining of the reference genome, assemblies of the newly sequenced genomes, inter- and intraspecific comparative genome analyses, and RNAseq analyses. SR managed the biological material and performed the plant-pathogen interaction experiment. DL and SR performed the nucleic acid extractions. FM and AK mediated the collaboration with DOE-JGI and provided expertise for the interspecific analysis. VP-R and RdO provided isolates from rubber tree and other hosts respectively. PL built and managed the RNAseq pipeline used in Clermont-Ferrand. BF provided expertise in phylogeny. J-SV provided expertise for the manual tuning of clustering parameters. KuL, AL, KaL, DB, RO, and KB, performed sequencing, assembly and annotation of the reference genome at DOE-JGI. JS and IG managed the 1,000 Fungal Genomes Project. DL and VP-R wrote the draft article. All authors read, corrected and approved the article.

## FUNDING

DL's salary was funded by IFC (Institut Français du Caoutchouc, Paris, France) and Michelin, SIPH and Socfin Companies. *C. cassiicola* isolate CCP sequencing, assembly, and annotation were performed by U.S. Department of Energy Joint Genome Institute. The work conducted by the U.S. Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, is supported by the Office of Science of the
U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Research in the laboratory of FM is funded by the Laboratory of Excellence Advanced Research on the Biology of Tree and Forest Ecosystems (ARBRE grant ANR 11 LABX 0002 01).

## ACKNOWLEDGMENTS

We thank the Metatranscriptomics of Forest Soil Ecosystems consortium for access to unpublished genomes of *Colletotrichum*

*gloesporioides* and *C. acutatum*, Irina S. Druzhinina, Tu Wien and Christian P. Kubicek for access to the unpublished genomes of *Trichoderma asperellum* and *T. harzianum*, and Michel Dron for scientific advices and support.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00276/full#supplementary-material

## REFERENCES

Altschul, S. F., Gish, W., Miller, W., Meyers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. *J. Mol. Biol.* 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

Aragona, M., Minio, A., Ferrarini, A., Valente, M. T., Bagnaresi, P., Orrù, L., et al. (2014). *De novo* genome assembly of the soil-borne fungus and tomato pathogen *Pyrenochaeta lycopersici*. *BMC Genomics* 15:313. doi: 10.1186/1471-2164-15-313

Atan, S., and Hamid, N. H. (2003). Differentiating races of *Corynespora cassiicola* using RAPD and internal transcribed spacer markers. *J. Rubber Res.* 6, 58–64.

Barthe, P., Pujade-Renaud, V., Breton, F., Gargani, D., Thai, R., Roumestand, C., et al. (2007). Structural analysis of cassiicolin, a host-selective protein toxin from *Corynespora cassiicola*. *J. Mol. Biol.* 367, 89–101. doi: 10.1016/j.jmb.2006.11.086

Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. *Ann. Stat.* 29, 1165–1188. doi: 10.1214/aos/1013699998

Benoit, I., Cullenon, H., Zhou, M., DiFalco, M., Aguilar-Osorio, G., Battaglia, E., et al. (2015). Closely related fungi employ diverse enzymatic strategies to degrade plant biomass. *Biotechnol. Biofuels* 8:107. doi: 10.1186/s13068-015-0285-0

Blin, K., Medema, M. H., Kazempour, D., Fischbach, M. A., Breitling, R., Takano, E., et al. (2013). antiSMASH 2.0-a versatile platform for genome mining of secondary metabolite producers. *Nucleic Acids Res.* 41, W204–W212. doi: 10.1093/nar/gkt449

Breton, F., Sanier, C., and d'Auzac, J. (2000). Role of cassiicolin, a host-selective toxin, in pathogenicity of *Corynespora cassiicola*, causal agent of a leaf fall disease of Hevea. *J. Rubber Res.* 3, 115–128.

Brown, A. E., and Soepena, H. (1994). Pathogenicity of *Colletotrichum acutatum* and *C. gloesporioides* on leaves of hevea spp. *Mycol. Res.* 98, 264–266. doi: 10.1016/S0953-7562(09)80453-X

Brown, K. J., Formolo, C. A., Seol, H., Marathi, R. L., Duguez, S., An, E., et al. (2012). Advances in the proteomic investigation of the cell secretome. *Expert Rev. Proteomics* 9, 337–345. doi: 10.1586/epr.12.21

Cai, L., Ji, K. F., and Hyde, K. D. (2006). Variation between freshwater and terrestrial fungal communities on decaying bamboo culms. *Antonie Van Leeuwenhoek* 89, 293–301. doi: 10.1007/s10482-005-9030-1

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. *BMC Bioinformatics* 10:421. doi: 10.1186/1471-2105-10-421

Cannon, P. F., Damm, U., Johnston, P. R., and Weir, B. S. (2012). *Colletotrichum* - current status and future directions. *Stud. Mycol.* 73, 181–213. doi: 10.3114/sim0014

Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. *Nucleic Acids Res.* 37, D233–D238. doi: 10.1093/nar/gkn663

Carlile, A. J., Bindschedler, L. V., Bailey, A. M., Bowyer, P., Clarkson, J. M., and Cooper, R. M. (2000). Characterization of SNPI, a cell wall-degrading trypsin, produced during infection by *Stagonospora nodorum*. *Mol. Plant Microbe Interact.* 13, 538–550. doi: 10.1094/MPMI.2000.13.5.538

Carris, L. M., Glawe, D. A., and Gray, L. E. (1986). Isolation of the soybean pathogens *Corynespora cassiicola* and *Phialophora gregata* from cysts of *Heterodera glycines* in Illinois. *Mycologia* 78, 503–506. doi: 10.2307/3793061

Castresana, J. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. *Mol. Biol. Evol.* 17, 540–552. doi: 10.1093/oxfordjournals.molbev.a026334

Cechin, A. L., Sinigaglia, M., Lemke, N., Echeverrigaray, S., Cabrera, O. G., Pereira, G. A., et al. (2008). Cupin: a candidate molecular structure for the Nep1-like protein family. *BMC Plant Biol.* 8:50. doi: 10.1186/1471-2229-8-50

Celio, G. J., Padamsee, M., Dentinger, B. T., Bauer, R., and McLaughlin, D. J. (2006). Assembling the fungal tree of life: constructing the structural and biochemical database. *Mycologia* 98, 850–859. doi: 10.1080/15572536.2006.11832615

Chang, S., Puryear, J., and Cairney, J. (1993). A simple and efficient method for isolating RNA from pine trees. *Plant Mol. Biol. Rep.* 11, 113–116. doi: 10.1007/BF02670468

Coleman, C., Copetti, D., Cipriani, G., Hoffmann, S., Kozma, P., Kovács, L., et al. (2009). The powdery mildew resistance gene REN1 co-segregates with an NBS-LRR gene cluster in two central Asian grapevines. *BMC Genet.* 10:89. doi: 10.1186/1471-2156-10-89

Collado, J., Platas, G., Gonzalez, I., and Pelaez, F. (1999). Geographical and seasonal influences on the distribution of fungal endophytes in *Quercus ilex*. *New Phytol.* 144, 525–532. doi: 10.1046/j.1469-8137.1999.00533.x

Comménil, P., Belingheri, L., and Dehorter, B. (1998). Antilipase antibodies prevent infection of tomato leaves by *Botrytis cinerea*. *Physiol. Mol. Plant Pathol.* 52, 1–14. doi: 10.1006/pmpp.1997.0130

Darriba, D., Taboada, G. L., Doallo, R., and Posada, D. (2011). ProtTest 3: fast selection of best-fit models of protein evolution. *Bioinformatics* 27, 1164–1165. doi: 10.1093/bioinformatics/btr088

Davies, K. A., De Lorono, I., Foster, S. J., Li, D., Johnstone, K., and Ashby, A. M. (2000). Evidence for a role of cutinase in pathogenicity of *Pyrenopeziza brassicae* on brassicas. *Physiol. Mol. Plant Pathol.* 57, 63–75. doi: 10.1006/pmpp.2000.0282

Deighton, F. (1936). A preliminary list of fungi and diseases of plants from Sierra Leone. *Kew Bull.* 7, 397–424.

de Lamotte, F., Duviau, M. P., Sanier, C., Thai, R., Poncet, J., Bieysse, D., et al. (2007). Purification and characterization of cassiicolin, the toxin produced by *Corynespora cassiicola*, causal agent of the leaf fall disease of rubber tree. *J. Chromatogr. B Analyt. Technol. Biomed. Life Sci.* 849, 357–362. doi: 10.1016/j.jchromb.2006.10.051

Déon, M., Bourré, Y., Gimenez, S., Berger, A., Bieysse, D., de Lamotte, F., et al. (2012a). Characterization of a cassiicolin-encoding gene from *Corynespora cassiicola*, pathogen of rubber tree (*Hevea brasiliensis*). *Plant Sci.* 185–186, 227–237. doi: 10.1016/j.plantsci.2011.10.017

Déon, M., Fumanal, B., Gimenez, S., Bieysse, D., Oliveira, R. R., Shuib, S. S., et al. (2014). Diversity of the cassiicolin gene in *Corynespora cassiicola* and relation with the pathogenicity in *Hevea brasiliensis*. *Fungal Biol.* 118, 32–47. doi: 10.1016/j.funbio.2013.10.011

Déon, M., Scomparin, A., Tixier, A., Mattos, C., Leroy, T., Seguin, M., et al. (2012b). First characterization of endophytic *Corynespora cassiicola* isolates with variant cassiicolin genes recovered from rubber trees in Brazil. *Fungal Divers.* 54, 87–99. doi: 10.1007/s13225-012-0169-6

Dixon, L. J., Schlub, R. L., Pernezny, K., and Datnoff, L. E. (2009). Host specialization and phylogenetic diversity of *Corynespora*
*casiicola*. *Phytopathology* 99, 1015-1027. doi: 10.1094/PHYTO-99-9-1015

Dunwell, J. M. (1998). Cupins: a new superfamily of functionally diverse proteins that include germins and plant storage proteins. *Biotechnol. Genet. Eng. Rev.* 15, 1-32. doi: 10.1080/02648725.1998.10647950

Enright, A. J., Van Dongen, S., and Ouzounis, C. A. (2002). An efficient algorithm for large-scale detection of protein families. *Nucleic Acids Res.* 30, 1575-1584. doi: 10.1093/nar/30.7.1575

Farr, D. F., and Rossman, A. Y. (2016). *Fungal Databases, U.S. National Fungus Collections, ARS, USDA*. Available online at: https://nt.ars-grin.gov/fungaldatabases/

Fellbrich, G., Romanski, A., Varet, A., Blume, B., Brunner, F., Engelhardt, S., et al. (2002). NPP1, a Phytophthora-associated trigger of plant defense in parsley and Arabidopsis. *Plant J.* 32, 375-390. doi: 10.1046/j.1365-313X.2002.01454.x

Fiegen, M., and Knogge, W. (2002). Amino acid alterations in isoforms of the effector protein NIP1 from *Rhynchosporium secalis* have similar effects on its avirulence- and virulence-associated activities on barley. *Physiol. Mol. Plant Pathol.* 61, 299-302. doi: 10.1006/pmpm.2002.0442

Fischer, S., Brunk, B. P., Chen, F., Gao, X., Harb, O. S., Iodice, J. B., et al. (2011). Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups. *Curr. Protoc. Bioinformatics* Chapter 6, Unit 6.12.1-19. doi: 10.1002/0471250935 bi0612s35

Fones, H. N., Eyles, C. J., Kay, W., Cowper, J., and Gurr, S. J. (2017). A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by *Zymoseptoria tritici*. *Fungal Genet. Biol.* 106, 51-60. doi: 10.1016/j.fgb.2017.07.002

Fravel, D., Olivain, C., and Alabouvette, C. (2003). *Fusarium oxysporum* and its biocontrol. *New Phytologist* 157, 493-502. doi: 10.1046/j.1469-8137.2003.00700.x

Garavillon-Tournayre, M., Lopez, D., Landi, M., Cochard, H., Venisse, J. S., Fumanal, B., et al. (2015). "High throuput workflow for RNAseq data treatment linking laboratory data server and remote parallel calculation platform," in *JOBIM*, eds P. Leroy and P. Peyret (Clermont-Ferrand). Available on line at: https://www6.inra.fr/jobim2015/Post-JOBIM/Affiches-Posters

Gazis, R., and Chaverri, P. (2010). Diversity of fungal endophytes in leaves and stems of wild rubber trees (*Hevea brasiliensis*) in Peru. *Fungal Ecol.* 3, 240-254. doi: 10.1016/j.funeco.2009.12.001

Gazis, R., Rehner, S., and Chaverri, P. (2011). Species delimitation in fungal endophyte diversity studies and its implications in ecological and biogeographic inferences. *Mol. Ecol.* 20, 3001-3013. doi: 10.1111/j.1365-294X.2011.05110.x

Gnerre, S., MacCallum, I., Przybylski, D., Ribeiro, F. J., Burton, J. N., Walker, B. J., et al. (2011). High-quality draft assemblies of mammalian genomes from massively parallel sequence data. *Proc. Natl. Acad. Sci.* 108, 1513-1518. doi: 10.1073/pnas.1017351108

Gond, S., Verma, V., Kumar, A., Kumar, V., and Kharwar, R. (2007). Study of endophytic fungal community from different parts of *Aegle marmelos* Correae (Rutaceae) from Varanasi (India). *World J. Microbiol. Biotechnol.* 23, 1371-1375. doi: 10.1007/s11274-007-9375-x

Gottwald, S., Samans, B., Lück, S., and Friedt, W. (2012). Jasmonate and ethylene dependent defence gene expression and suppression of fungal virulence factors: two essential mechanisms of Fusarium head blight resistance in wheat? *BMC Genomics* 13:369. doi: 10.1186/1471-2164-13-369

Grigoriev, I. V., Nikitin, R., Haridas, S., Kuo, A., Ohm, R., Otillar, R., et al. (2014). MycoCosm portal: gearing up for 1000 fungal genomes. *Nucleic Acids Res.* 42, D699-D704. doi: 10.1093/nar/gkt1183

Guindon, S., and Gascuel, O. (2003). A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. *Syst. Biol.* 52, 696-704. doi: 10.1080/10635150390235520

Hallé, F., and Martin, R. (1968). Etude de la croissance rythmique chez l'hévéa (*Hevea brasiliensis*) Müll. Arg., (Euphorbiacées, crotonoïdées). *Adansonia* 8, 475-503.

Han, M. V., Thomas, G. W., Lugo-Martinez, J., and Hahn, M. W. (2013). Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. *Mol. Biol. Evol.* 30, 1987-1997. doi: 10.1093/molbev/mst100

Hieu, N. D., Nghia, N. A., Chi, V. T. Q., and Dung, P. (2014). Genetic diversity and pathogenicity of *Corynespora cassiicola* isolates from rubber trees and other hosts in Vietnam. *J. Rubber Res.* 17, 187-203.

Hiruma, K., Gerlach, N., Sacristán, S., Nakano, R. T., Hacquard, S., Kracher, B., et al. (2016). Root endophyte *Colletotrichum tofieldiae* confers plant fitness benefits that are phosphate status dependent. *Cell* 165, 464-474. doi: 10.1016/j.cell.2016.02.028

Huang, H. K., Liu, C. E., Liou, J. H., Hsiue, H. C., Hsiao, C. H., and Hsueh, P. R. (2010). Subcutaneous infection caused by *Corynespora cassiicola*, a plant pathogen. *J. Infect.* 60, 188-190. doi: 10.1016/j.jinf.2009.11.002

Islam, M. S., Haque, M. S., Islam, M. M., Emdad, E. M., Halim, A., Hossen, Q. M., et al. (2012). Tools to kill: genome of one of the most destructive plant pathogenic fungi *Macrophomina phaseolina*. *BMC Genomics* 13:493. doi: 10.1186/1471-2164-13-493

Ismail, I. A., and Able, A. J. (2016). Secretome analysis of virulent *Pyrenophora teres* f. teres isolates. *Proteomics* 16, 2625-2636. doi: 10.1002/pmic.201500498

Kämper, J., Kahmann, R., Bölker, M., Ma, L. J., Brefort, T., Saville, B. J., et al. (2006). Insights from the genome of the biophytic fungal plant pathogen *Ustilago maydis*. *Nature* 444, 97-101. doi: 10.1038/nature05248

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* 30, 772-780. doi: 10.1093/molbev/mst010

Kingsland, G. C. (1985). Pathogenicity and epidemiology of *Corynespora cassiicola* in the Republic of the Seychelles. *Acta Horticult.* 153, 229-230. doi: 10.17660/ActaHortic.1985.153.31

Kuan, C. S., Yew, S. M., Toh, Y. F., Chan, C. L. Ngeow, Y. F., Lee, K. W., et al. (2015). Dissecting the fungal biology of *Bipolaris papendorfii*: from phylogenetic to comparative genomic analysis. *DNA Res.* 22, 219-232. doi: 10.1093/dnares/dsv007

Kubicek, C. P., Starr, T. L., and Glass, N. L. (2014). Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. *Annu. Rev. Phytopathol.* 52, 427-451. doi: 10.1146/annurev-phyto-102313-045831

Kuruvilla Jacob, C. (2006). *Corynespora Leaf Disease of Hevea Brasiliensis, Strategies for Management*. Kottayam: Rubber Research Institute of India.

Lee, S., Melnik, V., Taylor, J., and Crous, P. (2004). Diversity of saprobic hyphomycetes on Proteaceae and Restionaceae from South Africa. *Fungal Divers.* 17, 91-114.

Lé, S., Josse, J., and Husson, F. (2008). FactoMineR: An R package for multivariate analysis. *J. Stat. Softw.* 25, 1-18. doi: 10.18637/jss.v025.i01

Levasseur, A., Drula, E., Lombard, V., Coutinho, P. M., and Henrissat, B. (2013). Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. *Biotechnol. Biofuels* 6:41. doi: 10.1186/1754-6834-6-41

Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* 25, 1754-1760. doi: 10.1093/bioinformatics/btp324

Li, L., Stoeckert, C. J. Jr., and Roos, D. S. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes. *Genome Res.* 13, 2178-2189. doi: 10.1101/gr.1224503

Liu, T., Hou, J., Wang, Y., Jin, Y., Borth, W., Zhao, F., et al. (2016). Genome-wide identification, classification and expression analysis in fungal plant interactions of cutinase gene family and functional analysis of a putative ClCUT7 in *Curvularia lunata*. *Mol. Genet. Genomics* 291, 1105-1115. doi: 10.1007/s00438-016-1168-1

Liyanage, A. D. S., Jayasinghe, C. K., Liyanage, N. I. S., and Jayaratne, A. H. R. (1986). Corynespora leaf spot disease of rubber (*Hevea brasiliensis*) - a new record. *J. Rubber Res. Inst. Sri Lanka* 65, 47-50.

Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Res.* 42, D490-D495. doi: 10.1093/nar/gkt1178

Looi, H. K., Toh, Y. F., Yew, S. M., Na, S. L., Tan, Y. C., Chong, P. S., et al. (2017). Genomic insight into pathogenicity of dematiaceous fungus *Corynespora cassiicola*. *Peer J.* 5:e2841. doi: 10.7717/peerj.2841

Lopez, D., Bronner, G., Brunel, N., Auguin, D., Bourgerie, S., Brignolas, F., et al. (2012). Insights into Populus XIP aquaporins: evolutionary expansion, protein functionality, and environmental regulation. *J. Exp. Bot.* 63, 2217-2230. doi: 10.1093/jxb/err404
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* 15:550. doi: 10.1186/s13059-014-0550-8

Lowe, R. G., McCorkelle, O., Bleackley, M., Collins, C., Faou, P., Mathivanan, S., et al. (2015). Extracellular peptidases of the cereal pathogen *Fusarium graminearum*. *Front. Plant Sci.* 6:962. doi: 10.3389/fpls.2015.00962

Mahgoub, E. (1969). *Corynespora cassiicola*, a new agent of maduromycetoma. *J. Trop. Med. Hyg.* 72, 218-221.

Martin, J., Bruno, V. M., Fang, Z., Meng, X., Blow, M., Zhang, T., et al. (2010). Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads. *BMC Genomics* 11:663. doi: 10.1186/1471-2164-11-663

Mohanta, T. K., and Bae, H. (2015). The diversity of fungal genome. *Biol. Proced. Online* 17:8. doi: 10.1186/s12575-015-0020-z

Moreira, L. R., and Filho, E. X. (2008). An overview of mannan structure and mannan-degrading enzyme systems. *Appl. Microbiol. Biotechnol.* 79, 165-178. doi: 10.1007/s00253-008-1423-4

Nghia, N. A., Kadir, J., Sunderasan, E., Puad Abdullah, M., Malik, A., and Napis, S. (2008). Morphological and inter simple sequence repeat (ISSR) markers analyses of *Corynespora cassiicola* isolates from rubber plantations in Malaysia. *Mycopathologia* 166, 189-201. doi: 10.1007/s11046-008-9138-8

O'Connell, R. J., Thon, M. R., Hacquard, S., Amyotte, S. G., Kleemann, J., Torres, M. F., et al. (2012). Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses. *Nat. Genet.* 44, 1060-1065. doi: 10.1038/ng.2372

Ohm, R. A., Feau, N., Henrissat, B., Schoch, C. L., Horwitz, B. A., Barry, K. W., et al. (2012). Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. *PLoS Pathog.* 8:e1003037. doi: 10.1371/journal.ppat.1003037

Olivieri, F., Eugenia Zanetti, M., Oliva, C. R., Covarrubias, A. A., and Casalongué, C. A. (2002). Characterization of an extracellular serine protease of *Fusarium eumartii* and its action on pathogenesis related proteins. *Eur. J. Plant Pathol.* 108, 63-72. doi: 10.1023/A:1013920929965

Page, A. J., Taylor, B., Delaney, A. J., Soares, J., Seemann, T., Keane, J. A., et al. (2016). SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. *bioRxiv.* 2:e000056. doi: 10.1099/mgen.0.000056

Park, B. H., Karpinski, T. V., Syed, M. H., Leuze, M. R., and Uberbacher, E. C. (2010). CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. *Glycobiology* 20, 1574-1584. doi: 10.1093/glycob/cwq106

Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. *Nat. Methods* 8, 785-786. doi: 10.1038/nmeth.1701

Promputtha, I., Lumyong, S., Dhanasekaran, V., McKenzie, E. H., Hyde, K. D., and Jeewon, R. (2007). A phylogenetic evaluation of whether endophytes become saprotrophs at host senescence. *Microb. Ecol.* 53, 579-590. doi: 10.1007/s00248-006-9117-x

Qi, Y., Xie, Y., Zhang, X., Pu, J., Zhang, H., Huang, S., et al. (2009). Molecular and pathogenic variation identified among isolates of *Corynespora cassiicola*. *Mol. Biotechnol.* 41, 145-151. doi: 10.1007/s12033-008-9109-9

Qi, Y.-X., Zhang, X., Pu, J.-J., Liu, X.-M., Lu, Y., Zhang, H., et al. (2011). Morphological and molecular analysis of genetic variability within isolates of *Corynespora cassiicola* from different hosts. *Eur. J. Plant Pathol.* 130, 83-95. doi: 10.1007/s10658-010-9734-6

Quinlan, R. J., Sweeney, M. D., Lo Leggio, L., Otten, H., Poulsen, J. C., Johansen, K. S., et al. (2011). Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. *Proc. Natl. Acad. Sci. U.S.A.* 108, 15079-15084. doi: 10.1073/pnas.1105776108

Rawlings, N. D., Waller, M., Barrett, A. J., and Bateman, A. (2014). MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. *Nucleic Acids Res.* 42, D503-D509. doi: 10.1093/nar/gkt953

R Core Team. (2015). *R: A Language and Environment for Statistical Computing, in: R Foundation for Statistical Computing*. (Vienna, Austria).

Reis, H., Pfi ff i, S., and Hahn, M. (2005). Molecular and functional characterization of a secreted lipase from *Botrytis cinerea*. *Mol. Plant Pathol.* 6, 257-267. doi: 10.1111/j.1364-3703.2005.00280.x

Rodríguez-Gacio Mdel, C., Iglesias-Fernández, R., Carbonero, P., and Matilla, A. J. (2012). Softening-up mannan-rich cell walls. *J. Exp. Bot.* 63, 3976-3988. doi: 10.1093/jxb/ers096

Rodríguez-Herva, J. J., González-Melendi, P., Cuartas-Lanza, R., Antúnez-Lamas, M., Río-Alvarez, I., Li, Z., et al. (2012). A bacterial cysteine protease effector protein interferes with photosynthesis to suppress plant innate immune responses. *Cell. Microbiol.* 14, 669-681. doi: 10.1111/j.1462-5822.2012.01749.x

Schoch, C. L., Crous, P. W., Groenewald, J. Z., Boehm, E. W., Burgess, T. I., de Gruyter, J., et al. (2009). A class-wide phylogenetic assessment of Dothideomycetes. *Stud. Mycol.* 64, 1-15-S10. doi: 10.3114/sim.2009.64.01

Shuib, S. S., Déon, M., Mahyuddin, M. M., Izhar, A., Fumanal, B., Sunderasan, E., et al. (2015). Cassiicolin genes among Corynespora isolates from rubber plantations in Malaysia. *J. Rubber Res.* 18, 109-126.

Silva, W. P., Karunanayake, E. H., Wijesundera, R. L., and Priyanka, U. M. (2003). Genetic variation in *Corynespora cassiicola*: a possible relationship between host origin and virulence. *Mycol. Res.* 107, 567-571. doi: 10.1017/S0953756203007755

Silva, W. P. K., Deverall, B. J., and Lyon, B. R. (1998). Molecular, physiological and pathological characterization of Corynespora leaf spot fungi from rubber plantations in Sri Lanka. *Plant Pathol.* 47, 267-277. doi: 10.1046/j.1365-3059.1998.00245.x

Skamnioti, P., and Gurr, S. J. (2009). Against the grain: safeguarding rice from rice blast disease. *Trends Biotechnol.* 27, 141-150. doi: 10.1016/j.tibtech.2008.12.002

Slippers, B., and Wingfield, M. J. (2007). Botryosphaeriaceae as endophytes and latent pathogens of woody plants: diversity, ecology and impact. *Fungal Biol. Rev.* 21, 90-106. doi: 10.1016/j.fbr.2007.06.002

Smit, A. F. A., Hubley, R., and Green, P. (2013-2015). "RepeatMasker Open-4.0". Available online at http://repeatmasker.org

Sonnhammer, E. L., von Heijne, G., and Krogh, A. (1998). A hidden Markov model for predicting transmembrane helices in protein sequences. *Proc. Int. Conf. Intell. Syst. Mol. Biol.* 6, 175-182.

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. *Bioinformatics* 30, 1312-1313. doi: 10.1093/bioinformatics/btu033

Suryanarayanan, T., Murali, T., Thirunavukkarasu, N., Govinda Rajulu, M., Venkatesan, G., and Sukumar, R. (2011). Endophytic fungal communities in woody perennials of three tropical forest types of the Western Ghats, southern India. *Biodivers. Conserv.* 20, 913-928. doi: 10.1007/s10531-011-0004-5

Takahashi-Ando, N., Kimura, M., Kakeya, H., Osada, H., and Yamaguchi, I. (2002). A novel lactonohydrolase responsible for the detoxification of zearalenone: enzyme purification and gene cloning. *Biochem. J.* 365, 1-6. doi: 10.1042/bj20020450

Tange, O. (2011). GNU parallel: the command-line power tool. *login: The USENIX Magazine* 36, 42-47.

Tran, D. M., Clément-Demange, A., Déon, M., Garcia, D., Le Guen, V., Clément-Vidal, A., et al. (2016). Genetic determinism of sensitivity to *Corynespora cassiicola* exudates in rubber tree (*Hevea brasiliensis*). *PLoS ONE* 11:e0162807. doi: 10.1371/journal.pone.0162807

Tsigos, I., and Bouriotis, V. (1995). Purification and characterization of chitin deacetylase from Colletotrichum lindemuthianum. *J. Biol. Chem.* 270, 26286-26291. doi: 10.1074/jbc.270.44.26286

van Kan, J. A. L., van 't Klooster, J. W., Wagemakers, C. A. M., Dees, D. C. T., and van der Vlugt-Bergmans, C. J. B. (1997). Cutinase A of *Botrytis cinerea* is expressed, but not essential, during penetration of gerbera and tomato. *Mol. Plant-Microbe Interact.* 10, 30-38. doi: 10.1094/MPMI.1997.10.1.30

Warnes, M. G. R., Bolker, B., and Bonebakker, L. (2016). *Package 'gplots'*. Various R Programming Tools for Plotting Data.

Wheeler, T. J., Clements, J., Eddy, S. R., Hubley, R., Jones, T. A., Jurka, J., et al. (2013). Dfam: a database of repetitive DNA based on profile hidden Markov models. *Nucleic Acids Res.* 41, D70-D82. doi: 10.1093/nar/gks1265

Yamada, H., Takahashi, N., Hori, N., Asano, Y., Mochizuki, K., Ohkusu, K., et al. (2013). Rare case of fungal keratitis caused by *Corynespora cassiicola*. *J. Infect. Chemother.* 19, 1167-1169. doi: 10.1007/s10156-013-0579-8
Yan, X. X., Yu, C. P., Fu, X. A., Bao, F. F., Du, D. H., Wang, C., et al. (2016). CARD9 mutation linked to *Corynespora cassiicola* infection in a Chinese patient. *Br. J. Dermatol.* 174, 176–179. doi: 10.1111/bjd.14082

Zerbino, D. R., and Birney, E. (2008). Velvet: algorithms for de novo short read assembly using de Bruijn graphs. *Genome Res.* 18, 821–829. doi: 10.1101/gr.074492.107

Zhao, D. L., Shao, C. L., Gan, L. S., Wang, M., and Wang, C. Y. (2015). Chromone derivatives from a sponge-derived strain of the fungus *Corynespora cassiicola*. *J. Nat. Prod.* 78, 286–293. doi: 10.1021/np5009152

Zhao, Z., Liu, H., Wang, C., and Xu, J. R. (2014). Correction: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. *BMC Genomics* 15:6. doi: 10.1186/1471-2164-15-6

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Lopez, Ribeiro, Label, Fumanal. Venisse, Kohler, de Oliveira, Labutti, Lipzen, Lail, Bauer, Ohm, Barry, Spatafora, Grigoriev, Martin and Pujade-Renaud. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.