De novo transcriptome assembly, functional annotation, and expression profiling of rye (Secale cereale L.) hybrids inoculated with ergot (Claviceps purpurea)

Rye is used as food, feed, and for bioenergy production and remain an essential grain crop for cool temperate zones in marginal soils. Ergot is known to cause severe problems in cross-pollinated rye by contamination of harvested grains. The molecular response of the underlying mechanisms of this disease is still poorly understood due to the complex infection pattern. RNA sequencing can provide astonishing details about the transcriptional landscape, hence we employed a transcriptomic approach to identify genes in the underlying mechanism of ergot infection in rye. In this study, we generated de novo assemblies from twelve biological samples of two rye hybrids with identified contrasting phenotypic responses to ergot infection. The final transcriptome of ergot susceptible (DH372) and moderately ergot resistant (Helltop) hybrids contain 208,690 and 192,116 contigs, respectively. By applying the BUSCO pipeline, we confirmed that these transcriptome assemblies contain more than 90% of gene representation of the available orthologue groups at Virdiplantae odb10. We employed a de novo assembled and the draft reference genome of rye to count the differentially expressed genes (DEGs) between the two hybrids with and without inoculation. The gene expression comparisons revealed that 228 genes were linked to ergot infection in both hybrids. The genome ontology enrichment analysis of DEGs associated them with metabolic processes, hydrolase activity, pectinesterase activity, cell wall modification, pollen development and pollen wall assembly. In addition, gene set enrichment analysis of DEGs linked them to cell wall modification and pectinesterase activity. These results suggest that a combination of different pathways, particularly cell wall modification and pectinesterase activity contribute to the underlying mechanism that might lead to resistance against ergot in rye. Our results may pave the way to select genetic material to improve resistance against ergot through better understanding of the mechanism of ergot infection at molecular level. Furthermore, the sequence data and de novo assemblies are valuable as scientific resources for future studies in rye.

[1]  Juan Miguel García-Gómez,et al.  BIOINFORMATICS APPLICATIONS NOTE Sequence analysis Manipulation of FASTQ data with Galaxy , 2005 .

[2]  Longhua Zhou,et al.  Transcriptome Analysis Identifies Candidate Genes and Functional Pathways Controlling the Response of Two Contrasting Barley Varieties to Powdery Mildew Infection , 2019, International journal of molecular sciences.

[3]  Anthony M. Bolger,et al.  De Novo Transcriptome Analysis of Durum Wheat Flag Leaves Provides New Insights Into the Regulatory Response to Elevated CO2 and High Temperature , 2019, Front. Plant Sci..

[4]  L. Roden,et al.  A small-scale RNA isolation protocol useful for high-throughput extractions from recalcitrant plants. , 2010 .

[5]  Michael Gribskov,et al.  Comprehensive evaluation of de novo transcriptome assembly programs and their effects on differential gene expression analysis , 2016, Bioinform..

[6]  Vincenzo Lionetti,et al.  Overexpression of Pectin Methylesterase Inhibitors in Arabidopsis Restricts Fungal Infection by Botrytis cinerea1[C][W] , 2007, Plant Physiology.

[7]  Isabelle P. Oswald,et al.  Current Situation of Mycotoxin Contamination and Co-occurrence in Animal Feed—Focus on Europe , 2012, Toxins.

[8]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[9]  M. Geisler,et al.  Plant ABC Transporters , 2014, Signaling and Communication in Plants.

[10]  Carroll P. Vance,et al.  An RNA-Seq Transcriptome Analysis of Orthophosphate-Deficient White Lupin Reveals Novel Insights into Phosphorus Acclimation in Plants1[W][OA] , 2012, Plant Physiology.

[11]  P. Tudzynski,et al.  Polygalacturonase is a pathogenicity factor in the Claviceps purpurea/rye interaction. , 2002, Fungal genetics and biology : FG & B.

[12]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[13]  Marshall Nichols,et al.  Comparing reference-based RNA-Seq mapping methods for non-human primate data , 2014, BMC Genomics.

[14]  I. Nookaew,et al.  A comprehensive comparison of RNA-Seq-based transcriptome analysis from reads to differential gene expression and cross-comparison with microarrays: a case study in Saccharomyces cerevisiae , 2012, Nucleic acids research.

[15]  T. Miedaner,et al.  Biology, Genetics, and Management of Ergot (Claviceps spp.) in Rye, Sorghum, and Pearl Millet , 2015, Toxins.

[16]  Kenneth K. Lopiano,et al.  RNA-seq: technical variability and sampling , 2011, BMC Genomics.

[17]  Monika S. Doblin,et al.  De novo transcriptome assembly and analysis of differentially expressed genes of two barley genotypes reveal root-zone-specific responses to salt exposure , 2016, Scientific Reports.

[18]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[19]  V. Lionetti,et al.  Methyl esterification of pectin plays a role during plant-pathogen interactions and affects plant resistance to diseases. , 2012, Journal of plant physiology.

[20]  Transcriptome analysis of wheat seedling and spike tissues in the hybrid Jingmai 8 uncovered genes involved in heterosis , 2018, Planta.

[21]  Uwe Scholz,et al.  Towards a whole‐genome sequence for rye (Secale cereale L.) , 2017, The Plant journal : for cell and molecular biology.

[22]  T. Ouellet,et al.  Comparative Transcriptome Profiles of Near-Isogenic Hexaploid Wheat Lines Differing for Effective Alleles at the 2DL FHB Resistance QTL , 2018, Front. Plant Sci..

[23]  T. Miedaner,et al.  Ergot infection in winter rye hybrids shows differential contribution of male and female genotypes and environment , 2020, Euphytica.

[24]  K. S. Kim,et al.  A cytochrome P450 gene is differentially expressed in compatible and incompatible interactions between pepper (Capsicum annuum) and the anthracnose fungus, Colletotrichum gloeosporioides. , 1999, Molecular plant-microbe interactions : MPMI.

[25]  A. Ciccodicola,et al.  Uncovering the Complexity of Transcriptomes with RNA-Seq , 2010, Journal of biomedicine & biotechnology.

[26]  P. Tudzynski,et al.  Claviceps purpurea: molecular aspects of a unique pathogenic lifestyle. , 2004, Molecular plant pathology.

[27]  X. Zhang,et al.  Arabidopsis COBRA-LIKE 10, a GPI-anchored protein, mediates directional growth of pollen tubes. , 2013, The Plant journal : for cell and molecular biology.

[28]  P. Pevzner,et al.  An Eulerian path approach to DNA fragment assembly , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[29]  Hyong Woo Choi,et al.  Pepper pectin methylesterase inhibitor protein CaPMEI1 is required for antifungal activity, basal disease resistance and abiotic stress tolerance , 2008, Planta.

[30]  R. D'Ovidio,et al.  Pyramiding PvPGIP2 and TAXI-III But Not PvPGIP2 and PMEI Enhances Resistance Against Fusarium graminearum. , 2016, Molecular plant-microbe interactions : MPMI.

[31]  Hongwei Wang,et al.  Cloning and characterization of a novel UDP-glycosyltransferase gene induced by DON from wheat , 2015 .

[32]  S. Kelly,et al.  TransRate: reference-free quality assessment of de novo transcriptome assemblies , 2015, bioRxiv.

[33]  Justin Preece,et al.  De Novo Transcriptome Assembly and Analyses of Gene Expression during Photomorphogenesis in Diploid Wheat Triticum monococcum , 2014, PloS one.

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

[35]  T. Miedaner,et al.  Rye ( Secale cereale L.) , 2009 .

[36]  Michael Ott,et al.  De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity , 2013 .

[37]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[38]  T. Miedaner,et al.  Resistance to Ergot in Self‐incompatible Germplasm Resources of Winter Rye , 2009 .

[39]  S. Park,et al.  De novo transcriptome assembly of the Chinese pearl barley, adlay, by full-length isoform and short-read RNA sequencing , 2018, PloS one.

[40]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[41]  J. Fangel,et al.  The role of the cell wall in plant immunity , 2014, Front. Plant Sci..

[42]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[43]  P. Schröder,et al.  Glutathione S-Transferase Enzymes in Plant-Pathogen Interactions , 2018, Front. Plant Sci..

[44]  Runxuan Zhang,et al.  High-quality reference transcript datasets hold the key to transcript-specific RNA-sequencing analysis in plants. , 2017, The New phytologist.

[45]  R. D'Ovidio,et al.  Claviceps purpurea expressing polygalacturonases escaping PGIP inhibition fully infects PvPGIP2 wheat transgenic plants but its infection is delayed in wheat transgenic plants with increased level of pectin methyl esterification. , 2013, Plant physiology and biochemistry : PPB.

[46]  Tim H. Brom,et al.  A Reference-Free Algorithm for Computational Normalization of Shotgun Sequencing Data , 2012, 1203.4802.

[47]  A. Molina,et al.  Plant cell wall-mediated immunity: cell wall changes trigger disease resistance responses. , 2018, The Plant journal : for cell and molecular biology.

[48]  A. Dandekar,et al.  Deep RNA-Seq profile reveals biodiversity, plant–microbe interactions and a large family of NBS-LRR resistance genes in walnut (Juglans regia) tissues , 2016, AMB Express.

[49]  W. Bushuk.  Rye: Production, Chemistry, and Technology , 1976 .

[50]  Xuan Li,et al.  Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study , 2011, BMC Bioinformatics.

[51]  Richard J. Edwards,et al.  Plant adaptation or acclimation to rising CO2? Insight from first multigenerational RNA‐Seq transcriptome , 2016, Global change biology.

[52]  T. Miedaner,et al.  Genetic variation for resistance to ergot (Claviceps purpurea [Fr.] Tul.) among full-sib families of five populations of winter rye (Secale cereale L.) , 2008, Theoretical and Applied Genetics.

[53]  Caroline C. Friedel,et al.  A Comprehensive Evaluation of Alignment Algorithms in the Context of RNA-Seq , 2012, PloS one.

[54]  T. Miedaner,et al.  Genetic variation of winter rye cultivars for their ergot (Claviceps purpurea) reaction tested in a field design with minimized interplot interference , 2010 .

[55]  Jialei Duan,et al.  Optimizing de novo common wheat transcriptome assembly using short-read RNA-Seq data , 2012, BMC Genomics.

[56]  B. Moerschbacher,et al.  Immuno-affinity purification of PglPGIP1, a polygalacturonase-inhibitor protein from pearl millet: studies on its inhibition of fungal polygalacturonases and role in resistance against the downy mildew pathogen , 2015, Molecular Biology Reports.

[57]  Uwe Scholz,et al.  From RNA-seq to large-scale genotyping - genomics resources for rye (Secale cereale L.) , 2011, BMC Plant Biology.

[58]  R. Hückelhoven.  Cell wall-associated mechanisms of disease resistance and susceptibility. , 2007, Annual review of phytopathology.

[59]  V. Lionetti,et al.  The ectopic expression of a pectin methyl esterase inhibitor increases pectin methyl esterification and limits fungal diseases in wheat. , 2011, Molecular plant-microbe interactions : MPMI.

[60]  P. Tudzynski,et al.  Biotechnology and genetics of ergot alkaloids , 2001, Applied Microbiology and Biotechnology.

[61]  G. Piro,et al.  Three Pectin Methylesterase Inhibitors Protect Cell Wall Integrity for Arabidopsis Immunity to Botrytis1 , 2017, Plant Physiology.

[62]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[63]  Rodrigo L. Silveira,et al.  Plant biomass recalcitrance: effect of hemicellulose composition on nanoscale forces that control cell wall strength. , 2013, Journal of the American Chemical Society.