Transcriptome profiling of pumpkin (Cucurbita moschata Duch.) leaves infected with powdery mildew

Variability of plant transcriptomic responses under stress acclimation: a review from high throughput studies

Plant transcriptomes are complex entities shaped spatially and temporally by a multitude of stressors. The aim of this review was to summarize the most relevant transcriptomic responses to selected abiotic (UV radiation, chemical compounds, drought, suboptimal temperature) and biotic (bacteria, fungi, viruses, viroids) stress conditions in a variety of plant species, including model species, crops, and medicinal plants. Selected basic and applicative studies employing RNA-seq from various sequencing platforms and single-cell RNA-seq were involved. The transcriptomic responsiveness of various plant species and the diversity of affected gene families were discussed. Under stress acclimation, plant transcriptomes respond particularly dynamically. Stress response involved both distinct, but also similar gene families, depending on the species, tissue, and the quality and dosage of the stressor. We also noted the over-representation of transcriptomic data for some plant organs. Studies on plant transcriptomes allow for a better understanding of response strategies to environmental conditions. Functional analyses reveal the multitude of stress-affected genes as well as acclimatory mechanisms and suggest metabolome diversity, particularly among medicinal species. Extensive characterization of transcriptomic responses to stress would result in the development of new cultivars that would cope with stress more efficiently. These actions would include modern methodological tools, including advanced genetic engineering, as well as gene editing, especially for the expression of selected stress proteins in planta and for metabolic modifications that allow more efficient synthesis of secondary metabolites.

Transcriptome-based analysis of candidate gene markers associated with resistance mechanism to Phytophthora melonis that causes root and crown rot in pumpkin

Root and crown rot incited by an oomycete, Phytophthora melonis , causes significant yield losses in commercial pumpkin (Cucurbita pepo ) production worldwide. Currently, resistant cultivars and knowledge of molecular mechanism of C. pepo against P. melonis are scarce. Here, we analysed the quantitative gene expression changes of 10 candidate gene markers (bHLH87, ERF014, HSF, MYB, PR-1, WRKY21, CPI, POD, PSK, SGT ) in pumpkin roots and leaves at three time points (h post-inoculation, hpi) following inoculation with P. melonis in two resistant (Ghelyani and Tanbal), and two susceptible (Marmari and Khoreshti) varieties of pumpkin. Gene expression using quantitative real time PCR along a time course revealed the strongest transcriptomic response at 48 and 72hpi in resistant genotypes, 1.1-2.7-fold in roots and leaves, respectively, with a high significant correlation (r =0.857**-0.974**). We also found that CPI , PSK, SGT1 and POD act as a dual regulator that similarly modulate immunity not only against P. melonis , but also against other diseases such as early blight (Alternaria cucumerina) , powdery mildew (Podosphaera xanthii ), downy mildews (Pseudoperonospora cubensis ), and pathogenic plant nematodes (Meloidogyne javanica ). Furthermore, significantly higher activities of the ROS scavenging defence enzymes, catalase (1.6-fold increase) and peroxidase (6-fold increase) were observed in the roots of resistant cultivars at different hpi compared with non-inoculated controls. In addition, the biomass growth parameters including leaf and root length, stem and root diameter, root fresh weight and volume were significantly different among studied genotypes. Cumulatively, the transcriptome data provide novel insights into the response of pumpkins for improving pumpkin breeding to P. melonis .

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

A pathogenesis-related protein 1 of Cucurbita moschata responds to powdery mildew infection

Pumpkin (Cucurbita moschata Duch.) productivity is severely hindered by powdery mildew (PM) worldwide. The causative agent of pumpkin PM is Podosphaera xanthii, a biotrophic fungus. Pathogenesis-related protein 1 (PR1) homolog was previously identified from transcriptomic analysis of a PM-resistant pumpkin. Here, we investigated the effects of CmPR1 gene from pumpkin for resistance to PM. Subcellular localization assay revealed that CmPR1 is a cytoplasmic protein in plants. The expression of CmPR1 gene was strongly induced by P. xanthii inoculation at 48 h and exogenous ethylene (ET), jasmonic acid (JA) and NaCl treatments, but repressed by H2O2 and salicylic acid (SA) treatments. Visual disease symptoms, histological observations of fungal growth and host cell death, and accumulation of H2O2 in transgenic tobacco plants indicated that CmPR1 overexpression significantly enhanced the resistance to Golovinomyces cichoracearum compared to wild type plants during PM pathogens infection, possibly due to inducing cell death and H2O2 accumulation near infected sites. The expression of PR1a was significantly induced in transgenic tobacco plants in response to G. cichoracearum, suggesting that CmPR1 overexpression positively modulates the resistance to PM via the SA signaling pathway. These findings indicate that CmPR1 is a defense response gene in C. moschata and can be exploited to develop disease-resistant crop varieties.

Suppression of Chitin-Triggered Immunity by Plant Fungal Pathogens: A Case Study of the Cucurbit Powdery Mildew Fungus Podosphaera xanthii

Fungal pathogens are significant plant-destroying microorganisms that present an increasing threat to the world's crop production. Chitin is a crucial component of fungal cell walls and a conserved MAMP (microbe-associated molecular pattern) that can be recognized by specific plant receptors, activating chitin-triggered immunity. The molecular mechanisms underlying the perception of chitin by specific receptors are well known in plants such as rice and Arabidopsis thaliana and are believed to function similarly in many other plants. To become a plant pathogen, fungi have to suppress the activation of chitin-triggered immunity. Therefore, fungal pathogens have evolved various strategies, such as prevention of chitin digestion or interference with plant chitin receptors or chitin signaling, which involve the secretion of fungal proteins in most cases. Since chitin immunity is a very effective defensive response, these fungal mechanisms are believed to work in close coordination. In this review, we first provide an overview of the current understanding of chitin-triggered immune signaling and the fungal proteins developed for its suppression. Second, as an example, we discuss the mechanisms operating in fungal biotrophs such as powdery mildew fungi, particularly in the model species Podosphaera xanthii, the main causal agent of powdery mildew in cucurbits. The key role of fungal effector proteins involved in the modification, degradation, or sequestration of immunogenic chitin oligomers is discussed in the context of fungal pathogenesis and the promotion of powdery mildew disease. Finally, the use of this fundamental knowledge for the development of intervention strategies against powdery mildew fungi is also discussed.

Arbuscular mycorrhizal fungus changes alfalfa (Medicago sativa) metabolites in response to leaf spot (Phoma medicaginis) infection, with subsequent effects on pea aphid (Acyrthosiphon pisum) behavior

Plant disease occurs simultaneously with insect attack. Arbuscular mycorrhizal fungi (AMF) modify plant biotic stress response. Arbuscular mycorrhizal fungi and pathogens may modify plant volatile organic compound (VOC) production and insect behavior. Nevertheless, such effects are rarely studied, particularly for mesocosms where component organisms interact with each other. Plant-mediated effects of leaf pathogen (Phoma medicaginis) infection on aphid (Acyrthosiphon pisum) infestation, and role of AMF (Rhizophagus intraradices) in modifying these interactions were elucidated in a glasshouse experiment. We evaluated alfalfa disease occurrence, photosynthesis, phytohormones, trypsin inhibitor (TI) and total phenol response to pathogen and aphid attack, with or without AMF, and aphid behavior towards VOCs from AMF inoculated and non-mycorrhizal alfalfa, with or without pathogen infection. AM fungus enhanced alfalfa resistance to pathogen and aphid infestation. Plant biomass, root : shoot ratio, net photosynthetic rate, transpiration rate, stomatal conductance, salicylic acid, and TI were significantly increased in AM-inoculated alfalfa. Arbuscular mycorrhizal fungi and pathogen significantly changed alfalfa VOCs. Aphids preferred VOCs of AM-inoculated and nonpathogen-infected to nonmycorrhizal and pathogen-infected alfalfa. We propose that AMF alter plant response to multiple biotic stresses in ways both beneficial and harmful to the plant host, providing a basis for strategies to manage pathogens and herbivore pests.

Comparative RNA-seq analysis reveals a critical role for ethylene in rose (Rosa hybrida) susceptible response to Podosphera pannosa

Rose is one of the most important ornamental flowers, accounting for approximately one-third of the world's cut flower market. Powdery mildew caused by Podosphera pannosa is a devastating fungal disease in rose, mainly infecting the young leaves and causing serious economic losses. Therefore, a study on the mechanism of the fungus infecting the rose leaves and the possibility to improve resistance hereby is interesting and meaningful. Accordingly, we conducted transcriptome sequencing of rose leaves infected by P. pannosa at different time points to reveal the molecular mechanism of resistance to powdery mildew. The high-quality reads were aligned to the reference genome of Rosa chinensis, yielding 51,230 transcripts. A total of 1,181 differentially expressed genes (DEGs) were identified in leaves during P. pannosa infection at 12, 24, and 48 hpi. The transcription factors of ERF, MYB, bHLH, WRKY, etc., family were identified among DEGs, and most of them were downregulated during P. pannosa infection. The Kyoto Encyclopedia of Genes and Genomes analysis showed that the hormone signal transduction pathway, especially ethylene signal-related genes, was consistently showing a downregulated expression during powdery mildew infection. More importantly, exogenous 1-MCP (inhibitor of ethylene) treatment could improve the rose leaves' resistance to P. pannosa. In summary, our transcriptome of rose leaf infected by powdery mildew gives universal insights into the complex gene regulatory networks mediating the rose leaf response to P. pannosa, further demonstrating the positive role of 1-MCP in resistance to biotrophic pathogens.

Transcriptomic Insight into Viviparous Growth in Water Lily

Water lily is an important ornamental flower plant which is capable of viviparous plantlet development. But no study has been reported on the molecular basis of viviparity in water lily. Hence, we performed a comparative transcriptome study between viviparous water lily Nymphaea micrantha and a nonviviparous species Nymphaea colorata at four developmental stages. The higher expression of highly conserved AUX/IAA, ARF, GH3, and SAUR gene families in N. micrantha compared to N. colorata is predicted to have a major impact on the development and evolution of viviparity in water lily. Likewise, differential regulation of hormone signaling, brassinosteroid, photosynthesis, and energy-related pathways in the two species provide clues of their involvement in viviparity phenomenon. This study revealed the complex mechanism of viviparity trait in water lily. The transcriptomic signatures identified are important basis for future breeding and research of viviparity in water lily and other plant species.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

GWAS Reveals a Novel Candidate Gene CmoAP2/ERF in Pumpkin (Cucurbita moschata) Involved in Resistance to Powdery Mildew

Pumpkin (Cucurbita moschata Duchesne ex Poir.) is a multipurpose cash crop rich in antioxidants, minerals, and vitamins; the seeds are also a good source of quality oils. However, pumpkin is susceptible to the fungus Podosphaera xanthii, an obligate biotrophic pathogen, which usually causes powdery mildew (PM) on both sides of the leaves and reduces photosynthesis. The fruits of infected plants are often smaller than usual and unpalatable. This study identified a novel gene that involves PM resistance in pumpkins through a genome-wide association study (GWAS). The allelic variation identified in the CmoCh3G009850 gene encoding for AP2-like ethylene-responsive transcription factor (CmoAP2/ERF) was proven to be involved in PM resistance. Validation of the GWAS data revealed six single nucleotide polymorphism (SNP) variations in the CmoAP2/ERF coding sequence between the resistant (IT 274039 [PMR]) and the susceptible (IT 278592 [PMS]). A polymorphic marker (dCAPS) was developed based on the allelic diversity to differentiate these two haplotypes. Genetic analysis in the segregating population derived from PMS and PMR parents provided evidence for an incomplete dominant gene-mediated PM resistance. Further, the qRT-PCR assay validated the elevated expression of CmoAP2/ERF during PM infection in the PMR compared with PMS. These results highlighted the pivotal role of CmoAP2/ERF in conferring resistance to PM and identifies it as a valuable molecular entity for breeding resistant pumpkin cultivars.

Transcriptome analysis reveals the molecular mechanisms of response to an emergent yellow-flower disease in green Chinese prickly ash (Zanthoxylum schinifolium)

Chinese prickly ash (Zanthoxylum) is extensively used as spice and traditional medicine in eastern Asian countries. Recently, an emergent yellow-flower disease (YFD) break out in green Chinese prickly ash (Zanthoxylum schinifolium, Qinghuajiao in Chinese) at Chongqing municipality, and then leads to a sharp reduction in the yield of Qinghuajiao, and thus results in great economic losses for farmers. To address the molecular response for the emergent YFD of Qinghuajiao, we analyzed the transcriptome of 12 samples including the leaves and inflorescences of asymptomatic and symptomatic plants from three different towns at Chongqing by high-throughput RNA-Seq technique. A total of 126,550 genes and 229,643 transcripts were obtained, and 21,054 unigenes were expressed in all 12 samples. There were 56 and 164 different expressed genes (DEGs) for the AL_vs_SL (asymptomatic leaf vs symptomatic leaf) and AF_vs_SF (asymptomatic flower vs symptomatic flower) groups, respectively. The results of KEGG analysis showed that the “phenylpropanoid biosynthesis” pathway that related to plant–pathogen interaction were found in AL_vs_SL and AF_vs_SF groups, and the “Plant–pathogen interaction” found in AF_vs_SF group, implying that this Qinghuajiao YFD might cause by plant pathogen. Interestingly, we detected 33 common unigenes for the 2 groups, and almost these unigenes were up-regulated in the symptomatic plants. Moreover, most of which were homologs to virus RNA, the components of viruses, implying that this YFD was related to virus. Our results provided a primary molecular basis for the prevention and treatment of YFD of Qinghuajiao trees.

Gene Mining for Conserved, Non-Annotated Proteins of Podosphaera xanthii Identifies Novel Target Candidates for Controlling Powdery Mildews by Spray-Induced Gene Silencing

The powdery mildew fungus Podosphaera xanthii is one of the most important limiting factors for cucurbit production worldwide. Despite the significant efforts made by breeding and chemical companies, effective control of this pathogen remains elusive to growers. In this work, we examined the suitability of RNAi technology called spray-induced gene silencing (SIGS) for controlling cucurbit powdery mildew. Using leaf disc and cotyledon infiltration assays, we tested the efficacy of dsRNA applications to induce gene silencing in P. xanthii. Furthermore, to identify new target candidate genes, we analyzed sixty conserved and non-annotated proteins (CNAPs) deduced from the P. xanthii transcriptome in silico. Six proteins presumably involved in essential functions, specifically respiration (CNAP8878, CNAP9066, CNAP10905 and CNAP30520), glycosylation (CNAP1048) and efflux transport (CNAP948), were identified. Functional analysis of these CNAP coding genes by dsRNA-induced gene silencing resulted in strong silencing phenotypes with large reductions in fungal growth and disease symptoms. Due to their important contributions to fungal development, the CNAP1048, CNAP10905 and CNAP30520 genes were selected as targets to conduct SIGS assays under plant growth chamber conditions. The spray application of these dsRNAs induced high levels of disease control, supporting that SIGS could be a sustainable approach to combat powdery mildew diseases.

Identification and mapping of CpPM10.1, a major gene involved in powdery mildew (race 2 France of Podosphaera xanthii) resistance in zucchini (Cucurbita pepo L.)

Powdery mildew resistance in zucchini is controlled by one major dominant locus, CpPM10.1. CpPM10.1 was fine mapped. The expression of candidate gene Cp4.1LG10g02780 in resistant individuals was significantly upregulated after inoculation with the powdery mildew. Powdery mildew (PM) is one of the most destructive fungal diseases, reducing the productivity of Cucurbita crops globally. PM influences the photosynthesis, growth and development of infected zucchini and seriously reduces fruit yield and quality. In the present study, the zucchini inbred line 'X10' had highly stable PM resistance, and the inbred line 'Jin234' was highly susceptible to PM in the seedling stage and adult stages. Genetic analysis revealed that PM resistance in 'X10' is controlled by one major dominant locus. Based on the strategy of QTL-seq combined with linkage analysis and developed molecular markers, the major locus was found to be located in a 382.9-kb candidate region on chromosome 10; therefore, the major locus was named CpPM10.1. Using 1,400 F₂ individuals derived from a cross between 'X10' and 'JIN234' and F_2:3 offspring of the recombinants, the CpPM10.1 locus was defined in a region of approximately 20.9 kb that contained 5 coding genes. Among them, Cp4.1LG10g02780 contained a conserved domain (RPW8), which controls resistance to a broad range of PM pathogens. Cp4.1LG10g02780 also had nonsynonymous SNPs between the resistant 'X10' and susceptible 'Jin234.' Furthermore, the expression of Cp4.1LG10g02780 was strongly positively involved in PM resistance in the key period of inoculation. Further allelic diversity analysis in zucchini germplasm resources indicated that PM resistance was associated with two SNPs in the Cp4.1LG10g02780 RPW8 domain. This study not only provides highly stable PM resistance gene resources for cucurbit crops but also lays the foundation for the functional analysis of PM resistance and resistance breeding in zucchini.

Identification of microRNA-like RNAs from Trichoderma asperellum DQ-1 during its interaction with tomato roots using bioinformatic analysis and high-throughput sequencing

MicroRNA-like small RNAs (milRNAs) and their regulatory roles in the interaction between plant and fungus have recently aroused keen interest of plant pathologists. Trichoderma spp., one of the widespread biocontrol fungi, can promote plant growth and induce plant disease resistance. To investigate milRNAs potentially involved in the interaction between Trichoderma and tomato roots, a small RNA (sRNA) library expressed during the interaction of T. asperellum DQ-1 and tomato roots was constructed and sequenced using the Illumina HiSeqTM 2500 sequencing platform. From 13,464,142 sRNA reads, we identified 21 milRNA candidates that were similar to other known microRNAs in the miRBase database and 22 novel milRNA candidates that possessed a stable microRNA precursor hairpin structure. Among them, three milRNA candidates showed different expression level in the interaction according to the result of stem-loop RT-PCR indicating that these milRNAs may play a distinct regulatory role in the interaction between Trichoderma and tomato roots. The potential transboundary milRNAs from T. asperellum and their target genes in tomato were predicted by bioinformatics analysis. The results revealed that several interesting proteins involved in plant growth and development, disease resistance, seed maturation, and osmotic stress signal transduction might be regulated by the transboundary milRNAs. To our knowledge, this is the first report of milRNAs taking part in the process of interaction of T. asperellum and tomato roots and associated with plant promotion and disease resistance. The results might be useful to unravel the mechanism of interaction between Trichoderma and tomato.

Advances in Transcriptomics in the Response to Stress in Plants

Adverse stress influences the normal growth and development of plants. With the development of molecular biology technology, understanding the molecular mechanism of plants in response to adverse stress has gradually become an important topic for academic exploration. The expression of the transcriptome is dynamic, which reflects the level of expression of all genes in a particular cell, tissue, or organ of an individual organism at a particular stage of growth and development. Transcriptomics can disclose the expression at the whole genome level under stress from the whole transcriptional level, which can be useful in understanding the complex regulatory network associated with the adaptability and tolerance of plants to stress. In this article, we review the application of transcriptomics in understanding the response of plants to biotic stresses such as diseases and insect infestation and abiotic stresses such as water, temperature, salt, and heavy metals to provide a guideline for related research.

Expression of Pumpkin CmbHLH87 Gene Improves Powdery Mildew Resistance in Tobacco

Powdery mildew (PM), caused by Podosphaera xanthii, is a major threat to the global cucurbit yield. The molecular mechanisms underlying the PM resistance of pumpkin (Cucurbita moschata Duch.) are largely unknown. A homolog of the basic helix-loop-helix (bHLH) transcription factor was previously identified through a transcriptomic analysis of a PM-resistant pumpkin. In this study, this bHLH homolog in pumpkin has been functionally characterized. CmbHLH87 is present in the nucleus. CmbHLH87 expression in the PM-resistant material was considerably downregulated by PM; and abscisic acid, methyl jasmonate, ethephon, and NaCl treatments induced CmbHLH87 expression. Ectopic expression of CmbHLH87 in tobacco plants alleviated the PM symptoms on the leaves, accelerated cell necrosis, and enhanced H2O2 accumulation. The expression levels of PR1a, PR5, and NPR1 were higher in the PM-infected transgenic plants than in PM-infected wild-type plants. Additionally, the chlorosis and yellowing of plant materials were less extensive and the concentration of bacteria at infection sites was lower in the transgenic tobacco plants than in the wild-type plants in response to bacterial wilt and scab pathogens. CmbHLH87 may be useful for genetic engineering of novel pumpkin cultivars in the future.

Identification of MicroRNAs and Their Targets That Respond to Powdery Mildew Infection in Cucumber by Small RNA and Degradome Sequencing

Powdery mildew (PM) is a prevalent disease known to limit cucumber production worldwide. MicroRNAs (miRNAs) are single-stranded molecules that regulate host defense responses through posttranscriptional gene regulation. However, which specific miRNAs are involved and how they regulate cucumber PM resistance remain elusive. A PM-resistant single-segment substitution line, SSSL508-28, was developed previously using marker-assisted backcrossing of the PM-susceptible cucumber inbred D8 line. In this study, we applied small RNA and degradome sequencing to identify PM-responsive miRNAs and their target genes in the D8 and SSSL508-28 lines. The deep sequencing resulted in the identification of 156 known and 147 novel miRNAs. Among them, 32 and six differentially expressed miRNAs (DEMs) were detected in D8 and SSSL508-28, respectively. The positive correlation between DEMs measured by small RNA sequencing and stem-loop quantitative real-time reverse transcription-polymerase chain reaction confirmed the accuracy of the observed miRNA abundances. The 32 DEMs identified in the PM-susceptible D8 were all upregulated, whereas four of the six DEMs identified in the PM-resistant SSSL508-28 were downregulated. Using in silico and degradome sequencing approaches, 517 and 20 target genes were predicted for the D8 and SSSL508-28 DEMs, respectively. Comparison of the DEM expression profiles with the corresponding mRNA expression profiles obtained in a previous study with the same experimental design identified 60 and three target genes in D8 and SSSL508-28, respectively, which exhibited inverse expression patterns with their respective miRNAs. In particular, five DEMs were located in the substituted segment that contained two upregulated DEMs, Csa-miR172c-3p and Csa-miR395a-3p, in D8 and two downregulated DEMs, Csa-miR395d-3p and Csa-miR398b-3p, in SSSL508-28. One gene encoding L-aspartate oxidase, which was targeted by Csa-miR162a, was also located on the same segment and was specifically downregulated in PM-inoculated D8 leaves. Our results will facilitate the future use of miRNAs in breeding cucumber varieties with enhanced resistance to PM.

Transcriptomic Profiling of Acute Cold Stress-Induced Disease Resistance (SIDR) Genes and Pathways in the Grapevine Powdery Mildew Pathosystem

Temperatures from 2 to 8°C transiently induce quantitative resistance to powdery mildew in several host species (cold stress-induced disease resistance [SIDR]). Although cold SIDR events occur in vineyards worldwide an average of 14 to 21 times after budbreak of grapevine and can significantly delay grapevine powdery mildew (Erysiphe necator) epidemics, its molecular basis was poorly understood. We characterized the biology underlying the Vitis vinifera cold SIDR phenotype-which peaks at 24 h post-cold (hpc) treatment and results in a 22 to 28% reduction in spore penetration success-through highly replicated (n = 8 to 10) RNA sequencing experiments. This phenotype was accompanied by a sweeping transcriptional downregulation of photosynthesis-associated pathways whereas starch and sugar metabolism pathways remained largely unaffected, suggesting a transient imbalance in host metabolism and a suboptimal target for pathogen establishment. Twenty-six cold-responsive genes peaked in their differential expression at the 24-hpc time point. Finally, a subset of genes associated with nutrient and amino acid transport accounted for four of the eight most downregulated transcripts, including two nodulin 1A gene precursors, a nodulin MtN21 precursor, and a Dynein light chain 1 motor protein precursor. Reduced transport could exacerbate localized nutrient sinks that would again be transiently suboptimal for pathogen growth. This study links the transient cold SIDR phenotype to underlying transcriptional changes and provides an experimental framework and library of candidate genes to further explore cold SIDR in several systems, with an ultimate goal of identifying novel breeding or management targets for reduced disease.

Arbuscular mycorrhizal fungus alleviates alfalfa leaf spots caused by Phoma medicaginis revealed by RNA-seq analysis

AIMS

One of the major limitations to the production of alfalfa (Medicago sativa) is the fungus Phoma medicaginis, which infects alfalfa and causes leaf spots. This study aims to understand alfalfa's response to P. medicaginis infection, the colonization of arbuscular mycorrhizal fungus (AMF) and the effect of AMF on plant-pathogen interactions.

METHODS AND RESULTS

Transcriptome analysis (RNA-seq) was used to identify differentially expressed genes (DEGs) in alfalfa infected by P. medicaginis and colonized by AMF Rhizophagus intraradices. AMF ameliorated the effects of P. medicaginis infection on alfalfa by reducing leaf spot incidence and disease index by 39·48 and 56·18% respectively. Inoculation with pathogen and AMF induced the activity of defence pathways, including peroxidase (POD), polyphenol oxidase activities and jasmonic acid (JA), salicylic acid concentration. Plants showed differential expression of P. medicaginis resistance-related genes, including genes belonging to pathogenesis-related (PR) proteins, chitinase activity, flavonoid biosynthesis, phenylpropanoid biosynthesis, glutathione metabolism, phenylalanine metabolism and photosynthesis. Inoculation with AMF led to changes in the expression of genes involved in PR proteins, chitinase activity, phenylalanine metabolism and photosynthesis.

CONCLUSION

The physiological and transcriptional changes caused by P. medicaginis infection in non-mycorrhizal and mycorrhizal alfalfa provides crucial information for understanding AMF's association with pathogenic systems.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study showed that AMF alleviated alfalfa leaf spots demonstrating that AMF can serve as a biocontrol strategy for alfalfa disease management.

Improved Powdery Mildew Resistance of Transgenic Nicotiana benthamiana Overexpressing the Cucurbita moschata CmSGT1 Gene

Powdery mildew (PM), which is mainly caused by Podosphaera xanthii, is a serious biotrophic pathogen disease affecting field-grown and greenhouse-grown cucurbit crops worldwide. Because fungicides poorly control PM, the development and cultivation of PM-resistant varieties is critical. A homolog of SGT1 (suppressor of the G2 allele of skp1), which encodes a key component of the plant disease-associated signal transduction pathway, was previously identified through a transcriptomic analysis of a PM-resistant pumpkin (Cucurbita moschata) inbred line infected with PM. In this study, we have characterized this SGT1 homolog in C. moschata, and investigated its effects on biotic stress resistance. Subcellular localization results revealed that CmSGT1 is present in the nucleus. Additionally, CmSGT1 expression levels in the PM-resistant material was strongly induced by PM, salicylic acid (SA) and hydrogen peroxide (H2O2). In contrast, SA and H2O2 downregulated CmSGT1 expression in the PM-susceptible material. The ethephon (Eth) and methyl jasmonate (MeJA) treatments upregulated CmSGT1 expression in both plant materials. The constitutive overexpression of CmSGT1 in Nicotiana benthamiana (N. benthamiana) minimized the PM symptoms on the leaves of PM-infected seedlings, accelerated the onset of cell necrosis, and enhanced the accumulation of H2O2. Furthermore, the expression levels of PR1a and PR5, which are SA signaling transduction markers, were higher in the transgenic plants than in wild-type plants. Thus, the transgenic N. benthamiana plants were significantly more resistant to Erysiphe cichoracearum than the wild-type plants. This increased resistance was correlated with cell death, H2O2 accumulation, and upregulated expression of SA-dependent defense genes. However, the chlorosis and yellowing of plant materials and the concentration of bacteria at infection sites were greater in the transgenic N. benthamiana plants than in the wild-type plants in response to infections by the pathogens responsible for bacterial wilt and scab. Therefore, CmSGT1-overexpressing N. benthamiana plants were hypersensitive to these two diseases. The results of this study may represent valuable genetic information for the breeding of disease-resistant pumpkin varieties, and may also help to reveal the molecular mechanism underlying CmSGT1 functions.

Morphological, Transcriptomic and Hormonal Characterization of Trimonoecious and Subandroecious Pumpkin (Cucurbita maxima) Suggests Important Roles of Ethylene in Sex Expression

Sex expression is a complex process, and in-depth knowledge of its mechanism in pumpkin is important. In this study, young shoot apices at the one-true-leaf stage and 10-leaf stage in Cucurbita maxima trimonoecious line ‘2013–12’ and subandroecious line ‘9–6’ were collected as materials, and transcriptome sequencing was performed using an Illumina HiSeq^TM 2000 System. 496 up-regulated genes and 375 down-regulated genes were identified between shoot apices containing mostly male flower buds and only female flower buds. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the differentially expressed genes were mainly enriched in the ethylene and auxin synthesis and signal transduction pathways. In addition, shoot apices at the 4-leaf stage were treated with the ethylene-releasing agent 2-chloroethylphosphonic acid (Ethrel), aminoethoxyvinyl glycine (AVG), AgNO₃ and indoleacetic acid (IAA). The number of female flowers up to node 20 on the main stem of ‘2013–12’ increased significantly after Ethrel and IAA treatment and decreased significantly after AVG and AgNO₃ treatment. The female flowers in ‘9–6’ showed slight changes after treatment with the exogenous chemicals. The expression of key genes in ethylene synthesis and signal transduction (CmaACS7, CmaACO1, CmaETR1 and CmaEIN3) was determined using quantitative RT-PCR, and the expression of these four genes was positively correlated with the number of female flowers in ‘2013–12’. The variations in gene expression, especially that of CmaACS7, after chemical treatment were small in ‘9–6’. From stage 1 (S1) to stage 7 (S7) of flower development, the expression of CmaACS7 in the stamen was much lower than that in the ovary, stigma and style. These transcriptome data and chemical treatment results indicated that IAA might affect pumpkin sex expression by inducing CmaACS7 expression and indirectly affecting ethylene production, and the ethylene synthesis and signal transduction pathways play crucial roles in pumpkin flower sex expression. A possible reason for the differences in sex expression between pumpkin lines ‘2013–12’ and ‘9–6’ was proposed based on the key gene expression. Overall, these transcriptome data and chemical treatment results suggest important roles for ethylene in pumpkin sex expression.

RNA-seq analysis and fluorescence imaging of melon powdery mildew disease reveal an orchestrated reprogramming of host physiology

The cucurbit powdery mildew elicited by Podosphaera xanthii is one of the most important limiting factors in cucurbit production. Our knowledge of the genetic and molecular bases underlying the physiological processes governing this disease is very limited. We used RNA-sequencing to identify differentially expressed genes in leaves of Cucumis melo upon inoculation with P. xanthii, using RNA samples obtained at different time points during the early stages of infection and their corresponding uninfected controls. In parallel, melon plants were phenotypically characterized using imaging techniques. We found a high number of differentially expressed genes (DEGs) in infected plants, which allowed for the identification of many plant processes that were dysregulated by the infection. Among those, genes involved in photosynthesis and related processes were found to be upregulated, whereas genes involved in secondary metabolism pathways, such as phenylpropanoid biosynthesis, were downregulated. These changes in gene expression could be functionally validated by chlorophyll fluorescence imaging and blue-green fluorescence imaging analyses, which corroborated the alterations in photosynthetic activity and the suppression of phenolic compound biosynthesis. The powdery mildew disease in melon is a consequence of a complex and multifaceted process that involves the dysregulation of many plant pathways such as primary and secondary metabolism.

Transcriptome Analyses Provide Novel Insights into Heat Stress Responses in Chieh-Qua (Benincasa hispida Cogn. var. Chieh-Qua How)

Temperature rising caused by global warming has imposed significant negative effects on crop qualities and yields. To get the well-known molecular mechanism upon the higher temperature, we carefully analyzed the RNA sequencing-based transcriptomic responses of two contrasting chieh-qua genotypes: A39 (heat-tolerant) and H5 (heat-sensitive). In this study, twelve cDNA libraries generated from A39 and H5 were performed with a transcriptome assay under normal and heat stress conditions, respectively. A total of 8705 differentially expressed genes (DEGs) were detected under normal conditions (3676 up-regulated and 5029 down-regulated) and 1505 genes under heat stress (914 up-regulated and 591 down-regulated), respectively. A significant positive correlation between RNA-Seq data and qRT-PCR results was identified. DEGs related to heat shock proteins (HSPs), ubiquitin-protein ligase, transcriptional factors, and pentatricopeptide repeat-containing proteins were significantly changed after heat stress. Several genes, which encoded HSPs (CL2311.Contig3 and CL6612.Contig2), cytochrome P450 (CL4517.Contig4 and CL683.Contig7), and bHLH TFs (CL914.Contig2 and CL8321.Contig1) were specifically induced after four days of heat stress. DEGs detected in our study between these two contrasting cultivars would provide a novel basis for isolating useful candidate genes of heat stress responses in chieh-qua.