Genome Sequencing and Comparative Genomics of the Broad Host-Range Pathogen Rhizoctonia solani AG8

Rhizoctonia solani is a soil-borne basidiomycete fungus with a necrotrophic lifestyle which is classified into fourteen reproductively incompatible anastomosis groups (AGs). One of these, AG8, is a devastating pathogen causing bare patch of cereals, brassicas and legumes. R. solani is a multinucleate heterokaryon containing significant heterozygosity within a single cell. This complexity posed significant challenges for the assembly of its genome. We present a high quality genome assembly of R. solani AG8 and a manually curated set of 13,964 genes supported by RNA-seq. The AG8 genome assembly used novel methods to produce a haploid representation of its heterokaryotic state. The whole-genomes of AG8, the rice pathogen AG1-IA and the potato pathogen AG3 were observed to be syntenic and co-linear. Genes and functions putatively relevant to pathogenicity were highlighted by comparing AG8 to known pathogenicity genes, orthology databases spanning 197 phytopathogenic taxa and AG1-IA. We also observed SNP-level “hypermutation” of CpG dinucleotides to TpG between AG8 nuclei, with similarities to repeat-induced point mutation (RIP). Interestingly, gene-coding regions were widely affected along with repetitive DNA, which has not been previously observed for RIP in mononuclear fungi of the Pezizomycotina. The rate of heterozygous SNP mutations within this single isolate of AG8 was observed to be higher than SNP mutation rates observed across populations of most fungal species compared. Comparative analyses were combined to predict biological processes relevant to AG8 and 308 proteins with effector-like characteristics, forming a valuable resource for further study of this pathosystem. Predicted effector-like proteins had elevated levels of non-synonymous point mutations relative to synonymous mutations (dN/dS), suggesting that they may be under diversifying selection pressures. In addition, the distant relationship to sequenced necrotrophs of the Ascomycota suggests the R. solani genome sequence may prove to be a useful resource in future comparative analysis of plant pathogens.

[1]  Guangrui Huang,et al.  HaploMerger: Reconstructing allelic relationships for polymorphic diploid genome assemblies , 2012, Genome research.

[2]  G. Murray,et al.  Estimating disease losses to the Australian barley industry , 2010, Australasian Plant Pathology.

[3]  Shiv D. Kale.  Oomycete and fungal effector entry, a microbial Trojan horse. , 2012, The New phytologist.

[4]  R. H. Cruickshank,et al.  Pectic zymograms and taxonomy and pathogenicity of the Ceratobasidiaceae , 1986 .

[5]  A. Krogh,et al.  A combined transmembrane topology and signal peptide prediction method. , 2004, Journal of molecular biology.

[6]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[7]  X. Huang,et al.  CAP3: A DNA sequence assembly program. , 1999, Genome research.

[8]  W. Ko,et al.  Factors affecting protoplast formation by Rhizoctonia solani. , 2010, New biotechnology.

[9]  S. Eddy,et al.  tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. , 1997, Nucleic acids research.

[10]  G D Schuler,et al.  Electronic PCR: bridging the gap between genome mapping and genome sequencing. , 1998, Trends in biotechnology.

[11]  A. Goesmann,et al.  Establishment and interpretation of the genome sequence of the phytopathogenic fungus Rhizoctonia solani AG1-IB isolate 7/3/14. , 2013, Journal of biotechnology.

[12]  L. Travassos,et al.  Sialic acids in fungi: a minireview. , 1999, Glycoconjugate journal.

[13]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[14]  E. Selker,et al.  A cytosine methyltransferase homologue is essential for repeat-induced point mutation in Neurospora crassa , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Huanming Yang,et al.  De novo assembly of human genomes with massively parallel short read sequencing. , 2010, Genome research.

[16]  S. Kuninaga,et al.  The advancing identification and classification of Rhizoctonia spp. using molecular and biotechnological methods compared with the classical anastomosis grouping , 2006 .

[17]  H. Robinson,et al.  Protoplast preparation and transient transformation of Rhizoctonia solani , 2001 .

[18]  James K. Hane,et al.  Evolution of Linked Avirulence Effectors in Leptosphaeria maculans Is Affected by Genomic Environment and Exposure to Resistance Genes in Host Plants , 2010, PLoS pathogens.

[19]  G. Hartman,et al.  Response of Ancestral Soybean Lines and Commercial Cultivars to Rhizoctonia Root and Hypocotyl Rot. , 2001, Plant disease.

[20]  S. Brunak,et al.  SignalP 4.0: discriminating signal peptides from transmembrane regions , 2011, Nature Methods.

[21]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[22]  Sean D. Hooper,et al.  Gene Context Analysis in the Integrated Microbial Genomes (IMG) Data Management System , 2009, PloS one.

[23]  Y. Rahbé,et al.  The Phytopathogen Dickeya dadantii (Erwinia chrysanthemi 3937) Is a Pathogen of the Pea Aphid , 2006, Applied and Environmental Microbiology.

[24]  Christopher J. Rawlings,et al.  PHI-base: a new database for pathogen host interactions , 2005, Nucleic Acids Res..

[25]  D. MacLean,et al.  Genome analyses of the wheat yellow (stripe) rust pathogen Puccinia striiformis f. sp. tritici reveal polymorphic and haustorial expressed secreted proteins as candidate effectors , 2013, BMC Genomics.

[26]  Jerzy Jurka,et al.  Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor , 2006, BMC Bioinformatics.

[27]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[28]  M. Hood,et al.  Patterns of Repeat-Induced Point Mutation in Transposable Elements of Basidiomycete Fungi , 2012, Genome biology and evolution.

[29]  Richard A. Moore,et al.  Genome Comparison of Barley and Maize Smut Fungi Reveals Targeted Loss of RNA Silencing Components and Species-Specific Presence of Transposable Elements[W] , 2012, Plant Cell.

[30]  R. Hamelin,et al.  Single‐nucleotide polymorphism discovery in Leptographium longiclavatum, a mountain pine beetle‐associated symbiotic fungus, using whole‐genome resequencing , 2014, Molecular ecology resources.

[31]  U. Gisi,et al.  Variability in the Sensitivity of Rhizoctonia solani Anastomosis Groups to Fungicides , 1991 .

[32]  J. Jurka,et al.  Repbase Update, a database of eukaryotic repetitive elements , 2005, Cytogenetic and Genome Research.

[33]  Jin-Rong Xu,et al.  Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi , 2013, BMC Genomics.

[34]  S. Salzberg,et al.  Versatile and open software for comparing large genomes , 2004, Genome Biology.

[35]  A. Ogoshi,et al.  Identification of Rhizoctonia Species , 1991 .

[36]  Yulin Jia,et al.  Comparative analysis of putative pathogenesis-related gene expression in two Rhizoctonia solani pathosystems , 2011, Current Genetics.

[37]  I. Longden,et al.  EMBOSS: the European Molecular Biology Open Software Suite. , 2000, Trends in genetics : TIG.

[38]  Robert D. Finn,et al.  HMMER web server: interactive sequence similarity searching , 2011, Nucleic Acids Res..

[39]  Sarah Calvo,et al.  Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis , 2006, Nature.

[40]  Mark Borodovsky,et al.  Eukaryotic Gene Prediction Using GeneMark.hmm‐E and GeneMark‐ES , 2011, Current protocols in bioinformatics.

[41]  V. G. García,et al.  Review. Biology and systematics of the form genus Rhizoctonia , 2006 .

[42]  M. Adams,et al.  A tool for analyzing and annotating genomic sequences. , 1997, Genomics.

[43]  Christina A. Cuomo,et al.  Source (or Part of the following Source): Type Article Title Comparative Genomics Reveals Mobile Pathogenicity Chromosomes in Fusarium Author(s) , 2022 .

[44]  Walter Pirovano,et al.  BIOINFORMATICS APPLICATIONS , 2022 .

[45]  Jing Zhang,et al.  The evolution and pathogenic mechanisms of the rice sheath blight pathogen , 2013, Nature Communications.

[46]  M. Cheung,et al.  Rapid genotyping by low-coverage resequencing to construct genetic linkage maps of fungi: a case study in Lentinula edodes , 2013, BMC Research Notes.

[47]  K. Sivasithamparam,et al.  Improved method for protoplast regeneration of rhizoctonia solani , 1993 .

[48]  Rahul M Kohli,et al.  The curious chemical biology of cytosine: deamination, methylation, and oxidation as modulators of genomic potential. , 2012, ACS chemical biology.

[49]  P. Rougé,et al.  The Gal/GalNAc-specific lectin from the plant pathogenic basidiomycete Rhizoctonia solani is a member of the ricin-B family. , 2001, Biochemical and biophysical research communications.

[50]  T. Giraud,et al.  Repeat-Induced Point Mutation and the Population Structure of Transposable Elements in Microbotryum violaceum , 2005, Genetics.

[51]  A. Clutterbuck.  Genomic evidence of repeat-induced point mutation (RIP) in filamentous ascomycetes. , 2011, Fungal genetics and biology : FG & B.

[52]  M. Grynberg,et al.  LTR Retrotransposons in Fungi , 2011, PloS one.

[53]  M. Varrelmann,et al.  Susceptibility of intercrops to infection with Rhizoctonia solani AG 2-2 IIIB and influence on subsequently cultivated sugar beet , 2010 .

[54]  S Rozen,et al.  Primer3 on the WWW for general users and for biologist programmers. , 2000, Methods in molecular biology.

[55]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[56]  Helaine Carrer,et al.  A genome survey of Moniliophthora perniciosa gives new insights into Witches' Broom Disease of cacao , 2008, BMC Genomics.

[57]  James K. Hane,et al.  In silico reversal of repeat-induced point mutation (RIP) identifies the origins of repeat families and uncovers obscured duplicated genes , 2010, BMC Genomics.

[58]  Karam B. Singh,et al.  Interactions of Arabidopsis and M. truncatula with the same pathogens differ in dependence on ethylene and ethylene response factors , 2011, Plant signaling & behavior.

[59]  S. Natarajan,et al.  Gene expression profiling of the plant pathogenic basidiomycetous fungus Rhizoctonia solani AG 4 reveals putative virulence factors , 2012, Mycologia.

[60]  Mikael Bodén,et al.  MEME Suite: tools for motif discovery and searching , 2009, Nucleic Acids Res..

[61]  Matthew R. Pocock,et al.  The Bioperl toolkit: Perl modules for the life sciences. , 2002, Genome research.

[62]  Sonja J. Prohaska,et al.  Proteinortho: Detection of (Co-)orthologs in large-scale analysis , 2011, BMC Bioinformatics.

[63]  Paul Horton,et al.  Nucleic Acids Research Advance Access published May 21, 2007 WoLF PSORT: protein localization predictor , 2007 .

[64]  James K. Hane,et al.  A novel mode of chromosomal evolution peculiar to filamentous Ascomycete fungi , 2011, Genome Biology.

[65]  J. Schwartz,et al.  Current models of the mode of action of Bacillus thuringiensis insecticidal crystal proteins: a critical review. , 2012, Journal of invertebrate pathology.

[66]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

[67]  R. Dean,et al.  Mobile elements and mitochondrial genome expansion in the soil fungus and potato pathogen Rhizoctonia solani AG-3. , 2014, FEMS microbiology letters.

[68]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[69]  S. Schuster,et al.  Integrative analysis of environmental sequences using MEGAN4. , 2011, Genome research.

[70]  A PevznerPavel,et al.  De novo identification of repeat families in large genomes , 2005 .

[71]  E. Uberbacher,et al.  CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. , 2010, Glycobiology.

[72]  D. Ahrén,et al.  Genomic Mechanisms Accounting for the Adaptation to Parasitism in Nematode-Trapping Fungi , 2013, PLoS genetics.

[73]  Ying Li,et al.  Single Nucleus Genome Sequencing Reveals High Similarity among Nuclei of an Endomycorrhizal Fungus , 2014, PLoS genetics.

[74]  Christina A. Cuomo,et al.  Obligate Biotrophy Features Unraveled by the Genomic Analysis of the Rust Fungi, Melampsora larici-populina and Puccinia graminis f. sp. tritici , 2011 .

[75]  Paul R. Ebert,et al.  Antagonistic Interaction between Abscisic Acid and Jasmonate-Ethylene Signaling Pathways Modulates Defense Gene Expression and Disease Resistance in Arabidopsis , 2004, The Plant Cell Online.

[76]  György Abrusán,et al.  TEclass - a tool for automated classification of unknown eukaryotic transposable elements , 2009, Bioinform..

[77]  J. Keijer,et al.  Heterogeneity in electrophoretic karyotype within and between anastomosis groups of Rhizoctonia solani , 1996 .

[78]  T. Paulitz.  Low Input No-till Cereal Production in the Pacific Northwest of the U.S.: The Challenges of Root Diseases , 2006, European Journal of Plant Pathology.

[79]  Andersen,et al.  Biology and Systematics of the form genus Rhizoctonia , 2006 .

[80]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[81]  Christina A. Cuomo,et al.  Obligate biotrophy features unraveled by the genomic analysis of rust fungi , 2011, Proceedings of the National Academy of Sciences.

[82]  Pavel A. Pevzner,et al.  De novo identification of repeat families in large genomes , 2005, ISMB.

[83]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[84]  Karam B. Singh,et al.  Genetic and Genomic Analysis of Rhizoctonia solani Interactions with Arabidopsis; Evidence of Resistance Mediated through NADPH Oxidases , 2013, PloS one.

[85]  Steven Salzberg,et al.  BIOINFORMATICS ORIGINAL PAPER , 2004 .

[86]  Carl Kingsford,et al.  A fast, lock-free approach for efficient parallel counting of occurrences of k-mers , 2011, Bioinform..

[87]  Jonathan E. Allen,et al.  Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments , 2007, Genome Biology.

[88]  S. Schornack,et al.  Recent developments in effector biology of filamentous plant pathogens , 2010, Cellular microbiology.

[89]  E. Birney,et al.  Apollo: a sequence annotation editor , 2002, Genome Biology.

[90]  L. Travassos,et al.  Sialic acids in fungi , 1999, Glycoconjugate Journal.

[91]  P. Schulze-Lefert,et al.  Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts , 2013, Proceedings of the National Academy of Sciences.

[92]  U. Güldener,et al.  Pathogenicity Determinants in Smut Fungi Revealed by Genome Comparison , 2010, Science.

[93]  Karam B. Singh,et al.  The B-3 Ethylene Response Factor MtERF1-1 Mediates Resistance to a Subset of Root Pathogens in Medicago truncatula without Adversely Affecting Symbiosis with Rhizobia1[W][OA] , 2010, Plant Physiology.

[94]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[95]  Daniel Bessette,et al.  Population genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. , 2010, Genome research.

[96]  J. Fraser,et al.  Ploidy variation as an adaptive mechanism in human pathogenic fungi. , 2013, Seminars in cell & developmental biology.

[97]  M. Varrelmann,et al.  Maize genotype susceptibility to Rhizoctonia solani and its effect on sugar beet crop rotations , 2010 .

[98]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[99]  G. Murray,et al.  Estimating disease losses to the Australian wheat industry , 2009, Australasian Plant Pathology.