Comparative transcriptomics and genomic analyses reveal differential gene expression related to Colletotrichum brevisporum resistance in papaya (Carica papaya L.)

Colletotrichum brevisporum is an important causal pathogen of anthracnose that seriously affects the fruit quality and yield of papaya (Carica papaya L.). Although many genes and biological processes involved in anthracnose resistance have been reported in other species, the molecular mechanisms involved in the response or resistance to anthracnose in post-harvest papaya fruits remain unclear. In this study, we compared transcriptome changes in the post-harvest fruits of the anthracnose-susceptible papaya cultivar Y61 and the anthracnose-resistant cultivar G20 following C. brevisporum inoculation. More differentially expressed genes (DEGs) and differentially expressed long non-coding RNAs (DElnRNAs) were identified in G20 than in Y61, especially at 24 h post-inoculation (hpi), suggesting a prompt activation of defense responses in G20 in the first 24 h after C. brevisporum inoculation. These DEGs were mainly enriched in plant-pathogen interaction, phenylpropanoid biosynthesis/metabolism, and peroxisome and flavonoid biosynthesis pathways in both cultivars. However, in the first 24 hpi, the number of DEGs related to anthracnose resistance was greater in G20 than in Y61, and changes in their expression levels were faster in G20 than in Y61. We also identified a candidate anthracnose-resistant gene cluster, which consisted of 12 genes, 11 in G20 and Y61, in response to C. brevisporum inoculation. Moreover, 529 resistance gene analogs were identified in papaya genome, most of which responded to C. brevisporum inoculation and were genetically different between papaya cultivars and wild-type populations. The total expression dose of the resistance gene analogs may help papaya resist C. brevisporum infection. This study revealed the mechanisms underlying different anthracnose resistance between the anthracnose-resistant and anthracnose-susceptible cultivars based on gene expression, and identified some potential anthracnose resistance-related candidate genes/major regulatory factors. Our findings provided potential targets for developing novel genetic strategies to overcome anthracnose in papaya.

[1]  Yuhang Su,et al.  Combined transcriptome and metabolome analysis of Nerium indicum L. elaborates the key pathways that are activated in response to witches’ broom disease , 2022, BMC Plant Biology.

[2]  Q. Ke,et al.  Correction: Comparative transcriptome analysis reveals immunoregulation mechanism of lncRNA-mRNA in gill and skin of large yellow croaker (Larimichthys crocea) in response to Cryptocaryon irritans infection , 2022, BMC genomics.

[3]  Xingtan Zhang,et al.  SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya , 2022, Nature Genetics.

[4]  Madhu,et al.  Long Non-Coding RNAs as Emerging Regulators of Pathogen Response in Plants , 2022, Non-coding RNA.

[5]  Paul Theodor Pyl,et al.  Analysing high-throughput sequencing data in Python with HTSeq 2.0 , 2021, Bioinform..

[6]  F. Ordon,et al.  Quantitative Trait Loci Mapping of Adult Plant and Seedling Resistance to Stripe Rust (Puccinia striiformis Westend.) in a Multiparent Advanced Generation Intercross Wheat Population , 2021, Frontiers in Plant Science.

[7]  M. Le Paslier,et al.  Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits , 2021, Horticulture research.

[8]  Qiang Li,et al.  Transcription factor WRKY22 regulates canker susceptibility in sweet orange (Citrus sinensis Osbeck) by enhancing cell enlargement and CsLOB1 expression , 2021, Horticulture research.

[9]  Chun-Chao Wang,et al.  The landscape of gene-CDS-haplotype diversity in rice (Oryza sativa L.): properties, population organization, footprints of domestication and breeding, and implications in genetic improvement. , 2021, Molecular plant.

[10]  M. Oliva,et al.  An integrative approach reveals five new species of highland papayas (Caricaceae, Vasconcellea) from northern Peru , 2020, PloS one.

[11]  D. Weigel,et al.  Plant NLR diversity: the known unknowns of pan-NLRomes. , 2020, The Plant cell.

[12]  M. C. Gonçalves‐Vidigal,et al.  Integration of anthracnose resistance loci and RLK and NBS‐LRR‐encoding genes in the Phaseolus vulgaris L. genome , 2020 .

[13]  P. Bayer,et al.  Effect of Leptosphaeria maculans Infection on Promoter DNA Methylation of Defence Genes in Brassica napus , 2020, Agronomy.

[14]  R. Kuang,et al.  Genome-wide analysis of basic helix-loop-helix transcription factors in papaya (Carica papaya L.) , 2020, PeerJ.

[15]  Liying Sun,et al.  Regulation of phenylpropanoid biosynthesis by MdMYB88 and MdMYB124 contributes to pathogen and drought resistance in apple , 2020, Horticulture Research.

[16]  Curt R. Fischer,et al.  A Pathogen-Responsive Gene Cluster for Highly Modified Fatty Acids in Tomato , 2020, Cell.

[17]  E. L. de Souza,et al.  Application of coatings formed by chitosan and Mentha essential oils to control anthracnose caused by Colletotrichum gloesporioides and C. brevisporum in papaya (Carica papaya L.) fruit. , 2019, International journal of biological macromolecules.

[18]  Steven L Salzberg,et al.  Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype , 2019, Nature Biotechnology.

[19]  Changle Ma,et al.  Dynamics of Peroxisome Homeostasis and Its Role in Stress Response and Signaling in Plants , 2019, Front. Plant Sci..

[20]  Yuerong Liang,et al.  Involvement of Salicylic Acid in Anthracnose Infection in Tea Plants Revealed by Transcriptome Profiling , 2019, International journal of molecular sciences.

[21]  A. Idnurm,et al.  Identification of a gene cluster for the synthesis of the plant hormone abscisic acid in the plant pathogen Leptosphaeria maculans. , 2019, Fungal genetics and biology : FG & B.

[22]  H. Hirt,et al.  Plant Immunity: From Signaling to Epigenetic Control of Defense. , 2018, Trends in plant science.

[23]  C. Duan,et al.  First Report of Colletotrichum brevisporum Causing Anthracnose on Papaya in Taiwan. , 2018, Plant disease.

[24]  Gerald A Tuskan,et al.  Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus. , 2018, The New phytologist.

[25]  Dong Li,et al.  Full-Length Transcriptome Sequencing and Modular Organization Analysis of the Naringin/Neoeriocitrin-Related Gene Expression Pattern in Drynaria roosii , 2018, Plant & cell physiology.

[26]  O. Franco,et al.  Review: Potential biotechnological assets related to plant immunity modulation applicable in engineering disease-resistant crops. , 2018, Plant science : an international journal of experimental plant biology.

[27]  Nana Li,et al.  Transcriptional analysis and histochemistry reveal that hypersensitive cell death and H2O2 have crucial roles in the resistance of tea plant (Camellia sinensis (L.) O. Kuntze) to anthracnose , 2018, Horticulture Research.

[28]  Jia Gu,et al.  fastp: an ultra-fast all-in-one FASTQ preprocessor , 2018, bioRxiv.

[29]  M. Vega-García,et al.  Effect of the Combination Hot Water - Calcium Chloride on the In Vitro Growth of Colletotrichum gloeosporioides and the Postharvest Quality of Infected Papaya , 2017, The plant pathology journal.

[30]  Bostjan Kobe,et al.  Emerging Insights into the Functions of Pathogenesis-Related Protein 1. , 2017, Trends in plant science.

[31]  Xinnian Dong,et al.  Post-translational regulation of plant immunity. , 2017, Current opinion in plant biology.

[32]  D. Hofius,et al.  Autophagy as an emerging arena for plant-pathogen interactions. , 2017, Current opinion in plant biology.

[33]  Junfan Fu,et al.  Comparative Transcriptomics Reveals Differential Gene Expression Related to Colletotrichum gloeosporioides Resistance in the Octoploid Strawberry , 2017, Front. Plant Sci..

[34]  D. Horvath,et al.  Comprehensive Transcriptome Analyses Reveal Differential Gene Expression Profiles of Camellia sinensis Axillary Buds at Para-, Endo-, Ecodormancy, and Bud Flush Stages , 2017, Front. Plant Sci..

[35]  P. Vijayan,et al.  Transcriptome analysis reveals a complex interplay between resistance and effector genes during the compatible lentil-Colletotrichum lentis interaction , 2017, Scientific Reports.

[36]  J. Kelly,et al.  Transcriptome Profiling of the Phaseolus vulgaris - Colletotrichum lindemuthianum Pathosystem , 2016, PloS one.

[37]  S. Cloutier,et al.  RGAugury: a pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants , 2016, BMC Genomics.

[38]  P. Balint-Kurti,et al.  Maize Homologs of CCoAOMT and HCT, Two Key Enzymes in Lignin Biosynthesis, Form Complexes with the NLR Rp1 Protein to Modulate the Defense Response1 , 2016, Plant Physiology.

[39]  Huigang Hu,et al.  Transcriptome characterization and expression profiles of the related defense genes in postharvest mango fruit against Colletotrichum gloeosporioides. , 2016, Gene.

[40]  Asgar Ali,et al.  Efficacy of ginger oil and extract combined with gum arabic on anthracnose and quality of papaya fruit during cold storage , 2016, Journal of Food Science and Technology.

[41]  N. Hayashi,et al.  CC-NBS-LRR-Type R Proteins for Rice Blast Commonly Interact with Specific WRKY Transcription Factors , 2016, Plant Molecular Biology Reporter.

[42]  N. Hayashi,et al.  CC-NBS-LRR-Type R Proteins for Rice Blast Commonly Interact with Specific WRKY Transcription Factors , 2015, Plant Molecular Biology Reporter.

[43]  Md. Abu Sadat,et al.  Systematic characterization of the peroxidase gene family provides new insights into fungal pathogenicity in Magnaporthe oryzae , 2015, Scientific Reports.

[44]  A. Ali,et al.  Application of lemongrass oil in vapour phase for the effective control of anthracnose of ‘Sekaki’ papaya , 2015, Journal of applied microbiology.

[45]  M. Prasad,et al.  Genetic diversity analysis in a set of Caricaceae accessions using resistance gene analogues , 2014, BMC Genetics.

[46]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[47]  P. Chutichudet,et al.  Effects of Chitosan or Calcium Chloride on External Postharvest Qualities and Shelf-Life of ‘Holland’ Papaya Fruit , 2014 .

[48]  H. Deising,et al.  New gene models and alternative splicing in the maize pathogen Colletotrichum graminicola revealed by RNA-Seq analysis , 2014, BMC Genomics.

[49]  Jie Zhou,et al.  The perplexing role of autophagy in plant innate immune responses. , 2014, Molecular plant pathology.

[50]  T. R. Shojaei,et al.  Effect of pre-harvest calcium chloride applications on fruit calcium level and post-harvest anthracnose disease of papaya , 2014 .

[51]  K. Hyde,et al.  First Report of Papaya Fruit Anthracnose Caused by Colletotrichum brevisporum in Brazil. , 2013, Plant disease.

[52]  J. Kocher,et al.  CPAT: Coding-Potential Assessment Tool using an alignment-free logistic regression model , 2013, Nucleic acids research.

[53]  S. Manickam,et al.  Potential of chitosan‐loaded nanoemulsions to control different Colletotrichum spp. and maintain quality of tropical fruits during cold storage , 2012, Journal of applied microbiology.

[54]  Jian-ye Chen,et al.  Evaluation of New Reference Genes in Papaya for Accurate Transcript Normalization under Different Experimental Conditions , 2012, PloS one.

[55]  Jeremy D. DeBarry,et al.  MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity , 2012, Nucleic acids research.

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

[57]  B. Thomma,et al.  Of PAMPs and Effectors: The Blurred PTI-ETI Dichotomy[OA] , 2011, Plant Cell.

[58]  E. Bosquez-Molina,et al.  Inhibitory effect of essential oils against Colletotrichum gloeosporioides and Rhizopus stolonifer in stored papaya fruit and their possible application in coatings , 2010 .

[59]  Cyril Zipfel,et al.  News from the frontline: recent insights into PAMP-triggered immunity in plants. , 2008, Current opinion in plant biology.

[60]  Yong Zhang,et al.  CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine , 2007, Nucleic Acids Res..

[61]  J. P. Fabi,et al.  Papaya fruit ripening: response to ethylene and 1-methylcyclopropene (1-MCP). , 2007, Journal of agricultural and food chemistry.

[62]  S. Dillon,et al.  Development of a codominant CAPS marker linked to PRSV-P resistance in highland papaya , 2006, Theoretical and Applied Genetics.