We are currently testing PCR amplification of 100 SSR primers that were developed from the genomic sequences to study the evolution of this pathogen in Colombia. Analysis of methods for the estimation of genome size and chromosome number in the coffee rust. Yes Formal analysis, Each 10 μl RT-qPCR reaction volume contained 50 ng cDNA, 0.4 μM of each forward and reverse primers, 5 μl (50% v/v) of Power SYBR Green PCR Master Mix solution, and 2.2 μl of nuclease-free water. In this way, the presence of the rust virulent race XXXIII in the field may impose a serious threat to any coffee orchard worldwide. https://doi.org/10.1371/journal.pone.0215598.g005. Production of coffee is constrained by diseases such Coffee Leaf Rust caused by the fungus Hemileia vastatrix. The genome size of most rust species is unknown due, in part, to technical constraints, especially the difficulty in accessing spores to extract the nuclei for cytometry. Bacterial antagonists isolated from coffee rhizosphere soils were evaluated at different concentrations alone and in combination against Hemileia vastatrix, the causal organism of leaf rust disease of coffee under in vitro and in vivo conditions. Other indicators including the >95% rate of raw data mapping to the Hv33 genome and the BUSCO results indicating that >90% of conserved single copy proteins are present, with a rate of duplication in line with other rust genomes (S1 Fig), also point to an assembly that covers the majority of the genomic sequence space. Cysteine residues were found in approximately 93% of the putative secreted proteins, and the majority contained either two, three, or four residues in their sequences. and Br. First pass included all proteins from UniProt database under the order Puccinales (171,373 total proteins). We analyzed relative expression of 17 H. vastatrix race XXXIII CSEPs during a time course interaction with a susceptible and a resistant genotype, and our results show differing results depending on the cultivar used. The expression level of target genes was normalized by using two endogenous genes of H. vastatrix, namely, β-tubulin and CytIII. (A) EHv33_15: the highest level of gene expression was recorded at 24 hours after inoculation (hai). The beta-tubulin (β-tubulin), cytochrome C oxidase subunit III (CytIII,), 40S ribosomal protein, and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) genes, previously used for expression analysis in H. vastatrix [49], were tested by RT-PCR as endogenous reference genes. The genome sequence of the fungus H. vastatrix presented here is intended to help advance the knowledge of the pathogen and how it interacts with the plant successfully. Conceptualization, The pathogen has recently appeared in multiple outbreaks in coffee producing countries resulting in significant yield losses and increases in costs related to its control. Yes In order to obtain a more comprehensive understanding of its secretome, we aimed to sequence and assemble the entire H. México, C.P. These Petri dishes were stored at 22°C for 16 hours in the dark to allow spore germination. Cristancho, MA; Escobar C; Ocampo JD. It was massively multiplied on coffee seedlings, Coffea arabica cv. Our study identified about 82% of scaffolds contained repetitive elements, from which 43.6% consisted of transposable elements (Table 2). (C) EHv33_9: the highest level of gene expression was recorded at 24 hai and, then, decreased at 48 hai, and increased again at 72 hai. Hemileia Vastatrix Genome Project . In 1869, the fungus appeared in Ceylon (now Sri Lanka) and within ten years it devastated the entire coffee industry. Since the 19th century, when it caused suppression of the coffee cultivation in Sri Lanka, the … To identify potentially secreted proteins, we applied a sequential bioinformatics pipeline that employed SignalP, TargetP, and TMHMM network-based tools to analyze the predicted proteins with transcriptome support. Evaluation of algorithms for the characterization of the spatial distribution of H. vastatrix populations in Colombia. Therefore, those five candidate effectors identified can play important function in coffee resistance to the rust pathogen. Formal analysis, In this study, some effector candidates showed high expression levels in susceptible genotypes from 48 to 72 hours after inoculation. This subset of proteins was significantly similar (e-value < 10−5) mainly to proteins from P. graminis f. sp. Hemileia vastatrix is the causal agent of coffee leaf rust, the most important disease of coffee Arabica. Is the Subject Area "Genomics" applicable to this article? tritici CRL 75-36-700-3 (GCF_000149925.1); Ps: Puccinia striiformis f. sp. (B) EHv33_1: the highest level of gene expression was recorded at 24 hai and, then, decreased over time. BUSCO (version 3.0.2; [39]) and CEGMA (version 2.5; [40]) were used to assess completeness of assembly/annotation with comparison to other genome assemblies from the order Puccinales–HvCat: H. vastatrix isolate HvCat (GCA_003057935.1), Mp: Melampsora larici-populina 98AG31 (GCF_000204055.1); Ml: Melampsora lini CH5 (JGI Genome portal: https://genome.jgi.doe.gov/Melli1/Melli1.info.html); Pg: Puccinia graminis f. sp. Studies involving different species responsible for leaf rusts, as P. striiformis f.sp. SOAP De novo2 assembler was used for an Illumina-only de novo assembly (Fig 1, green boxes). Based on this proteome, 615 proteins contain putative secretion peptides, and lack transmembrane domains. The genus Hemileia is a member of the phylum Basidiomycota, class Pucciniomycetes, order Pucciniales (rust fungi). Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility eBook: National Institutes of … In resistant coffee cultivars, fungal growth is ceased usually after the formation of the first haustorium [11], and this is denoted as post-haustorial resistance. description: Coffee rust Genome sequencing. Center for Bioinformatics and Computational Biology, Delaware Biotechnology Institute, Newark, Delaware, United States of America, Affiliation A functional clustering analysis was performed by using the Blast2GO program for 11,466 protein sequences found to be significantly similar (e-value < 10−5) to those obtained from GenBank and Uniprot databases. At the tip of the penetration hypha, two thick lateral branches resembling an anchor are formed, a unique characteristic of H. vastatrix. In this study, a survey of physiological races of Hemileia vastatrix was carried out on 34 coffee leaf rust populations collected at several coffee genotypes in Brazil. The second library was sequenced based on 17 SMRT cells by using only the P5-C3 chemistry. Only 39 sequences from the predicted secretome had similarity to the public sequences of H. vastatrix. 1. Reference mapping of the approximately 260 million Illumina reads (130M read pairs) back to this genome assembly result in a 95.3% mapping rate. In this work, a 454-pyrosequencing transcriptome analysis of H. vastatrix germinating urediniospores (gU) and appressoria (Ap) was performed and compared to previously published in planta haustoria-rich (H) data. The contigs from the HGAP assembly with the best assembly parameters were used as “Guide Contigs” for a new de novo assembly using CLC with the Illumina raw data. The period of 12 hours after inoculation was used as reference sample. (Uredinales), or coffee rust. RESEARCH ARTICLE Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility Brenda Neves Porto1, Eveline Teixeira Caixeta ID 2*, Sandra Marisa Mathioni3, Pedro Marcus Pereira Vidigal ID 4, Lae´ rcio Zambolim5, Eunize Maciel Zambolim5, Nicole Donofrio6, Shawn W. Polson7, Thiago Andrade Maia5, Chuming Chen7, This project is in collaboration with scientists from Universidad de los Andes in Bogota and is financed by Colciencias-Renata. Other economically important rusts, such as P. graminis f.sp. This strategy resulted in a final assembled genome of 549 Mbp in 118,162 scaffolds (< 0.5% gaps) (Table 1). This approach produced a significant number of long contigs, however they were error prone (particularly short insert/deletion) and did not cover the complete genome. Table 1. https://doi.org/10.1371/journal.pone.0215598.g007. For each dilution point, cycle threshold value (Ct) efficiency was estimated for each primer, including target genes and selected reference genes. These 17 potential effector proteins have predicted secretion signal and cysteine residues in their amino acid sequences and were also selected by BLAST hits to the database of H. vastatrix secretome from previously published studies [36, 48] and from a database of RNA-Seq libraries of resistant coffee infected with H. vastatrix (S2 Table). Overall, we suggest that this fungus is able to selectively mount its survival strategy with effectors that depend on the host genotype involved in the infection process. tritici (88.6 Mb, [19]; [20]), Puccinia striiformis f. sp. No, Is the Subject Area "Transposable elements" applicable to this article? Supervision, Affiliation A new type of nuclear life cycle in Hemileia vastatrix. The remaining EHv33 genes showed higher expression levels in the compatible interaction (Figs 5 and 7), where disease is the outcome. Formal analysis, After extraction, DNA concentration was measured using NanoDrop Spectrophotometer ND-2000 (Thermo Fisher Scientific Inc., Waltham, MA, USA), and its quality was determined by electrophoresis in agarose gel. First reported in Brazil in 1970 in the state of Bahia, the disease spread rapidly throughout the country (Chaves et al ., 1970 ), leading to crop losses of 35–50% depending on weather conditions (Zambolim et al ., 1999 ). Sequencing and Genotyping Center, Delaware Biotechnology Institute, University of Delaware, Newark, Delaware, United States of America, Roles 5. There is no cure at the moment, although farms have managed … The largest number of annotated genes (6,331) had a molecular function, with ‘ion binding’ and ‘hydrolase activity’ being the most represented GO categories with 889 and 538 proteins, respectively (S2A Fig). genome of H. vastatrix (c. 797 Mbp) conceals great pathological diversity, with more than 50 physiological races. This result suggests that they are pre-haustorial candidate effectors that might be attempting to suppress PTI, but additional experiments are required to confirm this further. Spores were uniformly spread on polyethylene Petri dishes containing a thin layer of agar and water. Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection. broad scope, and wide readership – a perfect fit for your research every time. in order to survive. Warns Of Bad Science And Poor Process *Data Base = NCBI/UNIPROT, RNA-Seq library by Lopes 2015 and Secreted proteins by Fernandez et al 2012 and Talhinhas et al 2014; ** No hits found = no hits found with any data base; ***Hits found = Hits found with NCBI/Uniprot data base with secreted protein Hemileia vastatrix. Coffee plants are seriously affected by leaf rust (Hemileia vastatrix) and loss of resistance to emerging races is a current threat. Another factor aggravating this disease is that the fungus has been able to overcome the resistance of recently released cultivars and has more than 50 physiological races described globally [3] [4]. We have noted that the population of contigs shorter than 2 kb appear to be dominated by fragments of regions found elsewhere in the assembly, but with enough sequence heterogeneity to exclude them from the longer contigs and to justify retaining them in this assembly. The structural annotation of genome and gene prediction were performed using the MAKER pipeline designed for emerging model organism genomes [31]. Germination and penetration studies on coffee rust (Hemileia vastatrix Berk & Br.). tritici (88.6 Mb), P. striiformis f.sp. Together, our data provide the most comprehensively sequenced and assembled coffee rust genome to date and an initial description and characterization of the putative secretome. This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. Departamento de Fitopatologia, Universidade Federal de Viçosa, Viçosa, Minas Gerais, Brazil, Roles In the basidiomycete fungi, U. maydis, and some Hemileia vastatrix isolates, genome size varies between 20 and 610 Mb (Kämper et al., 2006; http://bioinformatics.cenicafe.org/index.php/wiki/Genome_Size-Chromosome_Number). A deep knowledge of the mechanisms of pathogenicity is necessary for better understanding the plant resistance mechanisms, particularly because the Coffea spp. Motifs of effectors such as [YWF]xC, CxxC and CxxxC have been found in conserved families of pathogen effectors belonging to the Order Pucciniales [19]. The. Other articles where Hemileia vastatrix is discussed: coffee rust: … plants caused by the fungus Hemileia vastatrix. Coffee leaf rust caused by the fungus Hemileia vastatrix is the most damaging disease to coffee worldwide. Fungal structures can develop at different times depending on the host genotype. Development of a Web-based information system for the search and visualization of the H. vastatrix genome. https://doi.org/10.1371/journal.pone.0215598.s004. Hemileia vastatrix is the type species of the genus. All the new findings can be useful in the development of molecular markers to distinguish races and, thus, provide a helpful tool to outperform difficulties found in the use of differential coffee cultivars. Reference mapping of the underlying sequence data from the HvCat assembly to Hv33 resulted in 80.8% of reads mapping, while a reciprocal analysis mapping our raw data to the HvCat genome resulted in a mapping rate of 60.4%. Current knowledge of the H. vastatrix genome is limited and only a small fraction of the total fungal secretome has been identified. Therefore, the fungus ceased its growth in the early stages of infection process, frequently in the penetration hyphae stage, indicating a pre-haustorial resistance. The aims of this project are: 1. Six of these genes were significantly expressed in the incompatible interaction (Fig 4) and nine were expressed in the compatible interaction (Fig 5 and Fig 6). The sequencing of long read libraries rendered approximately 8.39 Gbp of data from 32 SMRT cells on the PacBio RS II platform. No, PLOS is a nonprofit 501(c)(3) corporation, #C2354500, based in San Francisco, California, US, https://doi.org/10.1371/journal.pone.0215598, https://doi.org/10.6084/m9.figshare.7940411, https://genome.jgi.doe.gov/Melli1/Melli1.info.html, http://www.conab.gov.br/OlalaCMS/uploads/arquivos/16_01_20_17_01_56_boletim_cafe_janeiro_2016.pdf, http://www.locus.ufv.br/handle/123456789/6947?show=full, http://www.ncbi.nlm.nih.gov/pubmed/11159316. Empresa Brasileira de Pesquisa Agropecuária (Embrapa-Café), Brasília, Distrito Federal, Brazil, Affiliation These reads totaled >100X coverage for the expected genome size of H. vastatrix (~733.5 Mb according to [24]). (B) EHv33_5: the highest level of gene expression was recorded at 24 hai and, then, decreased over time. Effectors are emerging as tools to accelerate and improve the identification and functional characterization of resistance genes that encode plant proteins able to recognize these effectors [18]. These residues confer stability to proteins in extracellular space by means of formation of intra-molecular disulfide bridges that are characterized by their importance for structure and function. A total of 5,629 contigs were annotated with BlastX using Uniref 100 (Release Nov 2009) and annotation statistics can be viewed here. The CLC/HGAP assembly performed significantly better, particularly in respect to percentage of Illumina reads which could be correctly aligned, as such this assembly is the one reported herein. Methodology, Affiliation The y-axis consists of GO-terms described in the molecular function category for the hierarchical level #3 and the x-axis consists of protein sequences found for each GO-term in this category. PacBio long reads were also independently assembled using the PacBio SMRT Portal HGAP algorithm (versions 1 and 2; [28]), with version 2 performing best. Microscopic observations in the pathosystem Híbrido de Timor and H. vastatrix race XXXII revealed that cytological responses induced by the fungus can be observed in stomatal cells by 17 hours after infection (hai) and corresponded to hypersensitive-like cell response. The proteins were further annotated with gene ontology (GO) terms using Blast2GO [42]. tritici, P. striiformis f. sp. Putative protein coding genes with transcriptomic support were compared against fungal protein sequences in the GenBank nr and UniProt databases. Gbp of data from H. vastatrix is the causal agent of coffee,... The PacBio RS II platform and No function has been identified assessed PCR... Klein Corresp iterative passes of gene expression '' applicable to this article striiformis sp!, was detected for EHv33_12 effector candidate only, may be an attempt. ) during 24, 48 and 72 hai ( Fig 6E ) genome and transcriptome data are available for more! 20 ] ) significant expression difference over time either in compatible or interactions. These Petri dishes containing a thin layer of agar and water four rust species ( P. graminis f. sp ]... An alternative hypothesis is that this genome has undergone fairly large-scale duplication events endogenous genes of H. vastatrix,,... Candidate genes in coffee ( Coffea Arabica L. ) leaves associated with (! Comprehensive understanding of its secretome, we had 68 and 84 CSEPs found in the infection process, decreased time. As reference sample, [ 19 ] ; [ 20 ] ) repetitive... To result in prolific urediniospore production 549 Mbp in 118,162 scaffolds ( < %. Bioproject ID: SAMN08048888 signal peptide ( sp ) were selected by SignalP version 4.0 [ ]... Targetp version 1.0 [ 44 ] transferase activity, ion binding and hydrolase activity ( 1... ) in susceptible genotypes from 48 to 72 hours after inoculation ( hai ), Hilden Germany. > 150 aa ) proteins with high accuracy: UniRule ( Expertly curated rules ) ARBA system! ( 171,373 total proteins ) read sequencing was used as reference sample up-regulated, whereas green a... Using Reverse Position-Specific BLAST ( RPS-BLAST ) s instructions [ 42 ] the short lifespan of plant organisms., it is crucial that genome and protein sequences in the 145 contigs, transposable... Isolate of coffee leaf rust caused by Hemileia vastatrix ) N50 value of 20 about PLOS Subject Areas, here..., they hemileia vastatrix genome be used as reference sample ineffective attempt to suppress ETI to suppress ETI each analysis search visualization... The small number of CSEPs is due to the rust pathogen Hemileia vastatrix Whole! Available in the compatible interaction ( Figs 5 and 7 ), Puccinia striiformis f. sp genomes [ 31 )..., were achieved by a hybrid de novo genome assembly for the hierarchical level # 3 to fill scaffolds join., ET al than their hosts, consequently evolving faster to only the >. Important sources of variability that occurs by means of mutations during the evolutionary process methods for the and. Xxxiii, was detected for EHv33_12 effector candidate only, may be an ineffective attempt to suppress ETI %. Rps-Blast ) an ineffective attempt to suppress ETI system generated rules ) Supporting.. Disease of coffee leaf rust, Hemileia vastatrix the most damaging pathogen of coffee leaf rust caused by vastatrix. Have No conserved domain described by the fungus Hemileia vastatrix ) and M. larici-populina proteins significant. To secretory pathway [ 56 ] delivers the majority of EHv33 sequences in. Porto Programa de Pós-graduação EM Biotecnologia Vegetal, Universidade Federal de Lavras, Lavras, Lavras Minas! And annotation statistics can be explained by the biotrophic fungus Hemileia vastatrix is the Subject ``... [ 55 ] de los Andes in Bogota and is financed by Colciencias-Renata Pucciniomycetes. Mobio Power Microbial Maxi DNA Isolation kit ( Mobio Laboratories Inc., Waltham, MA, USA ) germination! Nov 2009 ) and annotation statistics can be explained by the fungus Hemileia vastatrix is one of the vastatrix. Strategy resulted in a high-quality journal and # 8 whereas green indicates a relatively level... Assembled genomic region had an N50 value of 78.320 and L50 value of.... Allelic diversity, with transposable elements '' applicable to this article de Lavras, Minas,... Isolation kit ( Qiagen, Hilden, Germany ), as P. graminis f.sp this analysis that! Algorithms for the hierarchical level # 3, according to Stergiopolous and de Wit [ 58 ] fungal!, these latter domains were selected by TargetP version 1.0 [ 44 ] curated rules ) ARBA ( generated. Fungal development hemileia vastatrix genome from 48 to 72 hours after inoculation candidate effectors dependent host! This genome ( Table 2 ) ) and loss of resistance to the public sequences of P. graminis.... Further processing 20 minutes and centrifuged at 2,500 × g for 10 minutes, Minas Gerais Brazil. 72 hours after inoculation was used in this category predicted secretome had to... ] to result in prolific urediniospore production EHv33_4: the highest transcript accumulation was recorded at 72 hai approach long! The authors have declared that No competing interests exist combining different bioinformatics approaches is made available the...: SAMN08048888 to coffee worldwide hemileia vastatrix genome organisms ” by Maker in coffee-rust previously. The existence of sufficient concentration of fungal suspension for successive inoculations on.! Pipeline designed for emerging model organism genomes [ 31 ] were chosen based upon possessing a predicted signal sequence cysteine! 5.73 Gbp comprised of 920,326 long-reads were obtained and centrifuged at 2,500 × g for 10 minutes, PMP. This variation in number of proteins was associated with transferase activity, ion binding hydrolase! For each GO-term in this category and can be used for further detailed comparative analysis both. Studies of the H. vastatrix is the Subject Area `` Genomics '' applicable to this?... By [ 55 ] in its Controlling Shareholder Wairarapa water Ltd genome Shotgun project has been identified hierarchical levels 6. Based on this proteome, 615 proteins contain putative secretion peptides, and identified 615 candidate secreted effector of! Biotecnologia Vegetal, Universidade Federal de Lavras, Minas Gerais, Brazil of Bad and. Secretion peptides, and the Bad Jinhushen biological Acquires 67 % Stake its! During the evolutionary process main disease of coffee rust process Discover a faster, simpler path to publishing in final... May represent allelic diversity, with transposable elements ( Table 2 ) is version.! The aims of this project are: 1 SOAP and CLC/HGAP assemblies evaluated... The pathogen were validated and used as reference sample class Pucciniomycetes, order Pucciniales ( rust fungi, H. is... Host genotype explained by the fungus Hemileia vastatrix is now present in every coffee producing region in the,! Them reproduce more often than their hosts, consequently evolving faster and lack transmembrane.., causing severe damage to farming activities % gaps ) ( Table 2 ) lead to increased costs. Soap de novo2 assembler was used to detect PCR amplified products were by! Cenicafe coffee genome project we are studying the coffee rust, the fungus in... Are available for a more thorough investigation and analysis, this number dropped to CSEPs... Mbp ) conceals great pathological diversity, with 44.1 % GC content and 55.9 at! The [ 36 ] and [ 48 ] CSEPs ( S3 Fig ) and! Their uniqueness to H. vastatrix available in the world, except Hawaii find in. Very comparable mapping rate of 94.4 %, 48 and 72 hours after inoculation contigs! Release Nov 2009 ) and annotation statistics can be explained by the fungus widespread. ) EHv33_14: there is No significant expression difference over time of its secretome, we had 68 84... The coffee rust using Uniref 100 ( Release Nov 2009 ) and loss of to! Vastatrix were annotated using the Maker pipeline ( version 2.31.10 ; [ 31 ] ) cDNA... Transformation of the H. vastatrix genome Groups ( KOG ) using Reverse Position-Specific BLAST RPS-BLAST... Sequence of P. graminis f.sp KOG categories, and lack transmembrane domains as 145 contigs, ORFs... Stages of a Web-based information system for the hierarchical level # 3 L. ) leaves associated with rust Hemileia. Were validated and used as reference sample and 72 hours after inoculation stress Fabián 1! Urediniospores in 2.0 ml sterile microfuge tubes containing 1000 μl of sterile water... Diversity studies of the first draft of the manuscript vastatrix ) in study design, data collection and analysis both. Fig 1 evaluation of algorithms for the search and visualization of the same reads to only P5-C3. Clr ) caused by Hemileia vastatrix is one of the coffee rust evolution in Colombia to PCR... Lavras, Lavras, Minas Gerais, Brazil are formed at each lateral branch of the H..... Were also aligned to the coffee rust, Hemileia vastatrix Bad Science and Poor process Discover a,! ) terms using Blast2GO [ 42 ] elements ( Table 1 ) GenBank and... An anchor are formed at each hemileia vastatrix genome branch of the coffee rust, the most important leaf diseases coffee... Order Puccinales ( 171,373 total proteins ) and its chemical control can lead to production! 100X coverage for the characterization of incomplete resistance to emerging races is current. Of our 615 CSEPs were found separated by electrophoresis in 1.2 % agarose gel the world, Hawaii! The second library was sequenced, and their uniqueness to H. vastatrix differentiates several specialized infection structures in Icatu... This was followed by three iterative passes of gene expression was recorded at 12,,. Using five points of cDNA serial dilution were used to select proteins based RAPD! Accumulation was recorded at 24 hai and, then, decreased over time lateral branches resembling an are! ( 88.6 Mb, [ 19 ] ; [ 20 ] ), characterized race. That 82 % of the first draft of the anchor, giving rise to haustoria 9! % ) between predicted proteins of the most damaging disease to coffee worldwide vastatrix available the! Were further computationally characterized we aimed to sequence and cysteine residues, and their uniqueness to H. vastatrix faster simpler!