Elucidation of the molecular responses of a cucumber segment substitution line carrying Pm5.1 and its recurrent parent triggered by powdery mildew by comparative transcriptome profiling

Genome-Wide Reidentification and Expression Analysis of MADS-Box Gene Family in Cucumber

MADS-box transcription factors play a crucial role in plant growth and development. Although previous genome-wide analyses have investigated the MADS-box family in cucumber, this study provides the first comprehensive reannotation of the MADS-box gene family in Cucumis sativus using updated Cucurbitaceae genome data, offering novel insights into the gene family's evolution and functional diversity. The results show that a total of 48 CsMADS-box genes were identified in the V3 version of cucumber, while 3 of the 43 genes identified in the V1 version were duplicated. The V1 version actually has only 40 genes. Additionally, we analyzed the variability in protein sequences and found that the amino acid sequences of 14 genes showed no differences between the two versions of the database, while the amino acid sequences of 29 genes exhibited significant differences. The further analysis of conserved motifs revealed that although the amino acid lengths of 15 genes had changed, their conserved motifs remained unchanged; however, the conserved motifs of 12 genes had altered. Furthermore we found that motif1 and motif2 were present in most proteins, indicating that they are highly conserved. Gene structure analysis revealed that most type I (Mα, Mβ) MADS-box genes lack introns, whereas type II (MIKC) genes exhibit a similar structure with a higher number of introns. Chromosomal localization analysis indicated that CsMADS-box genes are unevenly distributed across the seven chromosomes of cucumber. Promoter region analysis showed that the promoter regions of CsMADS-box genes contain response elements related to plant growth and development, suggesting that CsMADS-box genes may be extensively involved in plant growth and development. Different CsMADS-box genes exhibit specific high expression in roots, stems, leaves, tendrils, male flowers, female flowers, and ovaries, suggesting that these genes play crucial roles in the growth, development, reproduction and morphogenesis of cucumber. Moreover, 26, 18, 8, and 10 CsMADS-box genes were differentially expressed under high temperature, NaCl and/or silicon, downy mildew, and powdery mildew treatments, respectively. Interestingly, CsMADS07 and CsMADS16 responded to all tested stress conditions. These findings provide a reference and basis for further investigation into the function and mechanisms of the MADS-box genes for resistance breeding in cucumber.

In-silico identification of putatively functional intergenic small open reading frames in the cucumber genome and their predicted response to biotic and abiotic stresses

The availability of high-throughput sequencing technologies increased our understanding of different genomes. However, the genomes of all living organisms still have many unidentified coding sequences. The increased number of missing small open reading frames (sORFs) is due to the length threshold used in most gene identification tools, which is true in the genic and, more importantly and surprisingly, in the intergenic regions. Scanning the cucumber genome intergenic regions revealed 420 723 sORF. We excluded 3850 sORF with similarities to annotated cucumber proteins. To propose the functionality of the remaining 416 873 sORF, we calculated their codon adaptation index (CAI). We found 398 937 novel sORF (nsORF) with CAI ≥ 0.7 that were further used for downstream analysis. Searching against the Rfam database revealed 109 nsORFs similar to multiple RNA families. Using SignalP-5.0 and NLS, identified 11 592 signal peptides. Five predicted proteins interacting with Meloidogyne incognita and Powdery mildew proteins were selected using published transcriptome data of host-pathogen interactions. Gene ontology enrichment interpreted the function of those proteins, illustrating that nsORFs' expression could contribute to the cucumber's response to biotic and abiotic stresses. This research highlights the importance of previously overlooked nsORFs in the cucumber genome and provides novel insights into their potential functions.

Identification of MYC genes in four Cucurbitaceae species and their roles in the response to temperature stress

BACKGROUND

Myelocytomatosis (MYC) transcription factors are crucial mediators of the response of plants to environmental stresses through via binding to DNA regulatory regions. However, few systematic characterizations of MYC genes are available in Cucurbitaceae species.

RESULTS

In this study, we identified 10, 8, 12, and 10 MYC genes in Cucumis sativus, Cucumis melo, Citrullus lanatus, and Benincasa hispida, respectively. Characterization revealed that all of the MYC proteins contain a highly conserved H4-V5-E6-E8-R9-R11-R12 sequence, which is essential for the binding of DNA regulatory regions. Evolutionary analysis enabled us to categorize 40 predicted MYC proteins from seven species into five distinct groups and revealed that the expansion of the MYC genes occurred before the divergence of monocots and dicots. The upstream promoter regions of the MYC genes contain a variety of developmental, stress, and hormone-responsive regulatory elements. The expression of cucumber MYC genes varies significantly across organs, with particularly high expression of CsaV3_3G001710 observed across all organs. Transcriptomic analysis revealed that certain cucumber MYC genes undergo specific upregulation or downregulation in response to both biotic and abiotic stressors. In particular, under temperature stress, the cucumber genes CsaV3_3G007980 and CsaV3_3G001710 were significantly upregulated. Interestingly, the homologs of these two genes in C. lanatus presented a similar expression pattern to that in C. sativus, whereas in B. hispida, they presented the opposite pattern, i.e., significant downregulation. These findings indicated that these two genes indeed respond to temperature stress but with different expression patterns, highlighting the divergent functions of homologous genes across different species.

CONCLUSIONS

This study analyzed the size and composition of the MYC gene family in four Cucurbitaceae species and investigated stress-responsive expression profiles, especially under temperature stress. All the results showed that MYC genes play important roles in development and stress responses, laying a theoretical foundation for further investigations of these response mechanisms.

Genome-Wide Identification and Expression Analysis of the Cucumber FKBP Gene Family in Response to Abiotic and Biotic Stresses

The FKBP (FK506-binding protein) gene family is an important member of the PPlase protease family and plays a vital role during the processes of plant growth and development. However, no studies of the FKBP gene family have been reported in cucumber. In this study, 19 FKBP genes were identified in cucumber, which were located on chromosomes 1, 3, 4, 6, and 7. Phylogenetic analysis divided the cucumber FKBP genes into three subgroups. The FKBP genes in the same subgroup exhibited similar structures and conserved motifs. The cis-acting elements analysis revealed that the promoters of cucumber FKBP genes contained hormone-, stress-, and development-related cis-acting elements. Synteny analysis of the FKBP genes among cucumber, Arabidopsis, and rice showed that 12 kinds of syntenic relationships were detected between cucumber and Arabidopsis FKBP genes, and 3 kinds of syntenic relationships were observed between cucumber and rice FKBP genes. The tissue-specific expression analysis showed that some FKBP genes were expressed in all tissues, while others were only highly expressed in part of the 10 types of tissues. The expression profile analysis of cucumber FKBP genes under 13 types of stresses showed that the CsaV3_1G007080 gene was differentially expressed under abiotic stresses (high temperature, NaCl, silicon, and photoperiod) and biotic stresses (downy mildew, green mottle mosaic virus, Fusarium wilt, phytophthora capsica, angular leaf spot, and root-knot nematode), which indicated that the CsaV3_1G007080 gene plays an important role in the growth and development of cucumber. The interaction protein analysis showed that most of the proteins in the FKBP gene family interacted with each other. The results of this study will lay the foundation for further research on the molecular biological functions of the cucumber FKBP gene family.

Genome-wide characterization, phylogenetic and expression analysis of ABCG gene subfamily in cucumber (Cucumis sativus L.)

The ABCG is the largest subfamily of the ABC family with extensive functions, and only a few members have been identified in detail. However, more and more studies have shown that the members of this family are very important and are involved in many life processes such as plant development and response to various stresses. Cucumber is an important vegetable crops around the world. The cucumber development is essential for its production and quality. Meanwhile, various stresses have caused serious losses of cucumber. However, the ABCG genes were not well characterized and functioned in cucumber. In this study, the cucumber CsABCG gene family were identified and characterized, and their evolutionary relationship and functions were analyzed. The cis-acting elements and expression analysis showed that they played important role in development and responding to various biotic and abiotic stresses in cucumber. Phylogenetic analysis, sequence alignment and MEME (Multiple Em for Motif Elicitation) analysis indicated that the functions of ABCG proteins in different plants are evolutionarily conserved. Collinear analysis revealed that the ABCG gene family was highly conserved during the evolution. In addition, the potential binding sites of the CsABCG genes targeted by miRNA were predicted. These results will lay a foundation for further research on the function of the CsABCG genes in cucumber.

Genome-wide characterization, phylogenetic and expression analysis of Histone gene family in cucumber (Cucumis sativus L.)

Histones are essential components of chromatin and play an important role in regulating gene transcription and participating in DNA replication. Here, we performed a comprehensive analysis of this gene family. In this study, we identified 37 CsHistones that were classified into five groups (H1, H2A, H2B, H3 and H4). The closely linked subfamilies exhibited more similarity in terms of motifs and intron/exon numbers. Segmental duplication (SD) is the main driving force of cucumber CsHistones expansion. Analysis of cis-regulatory elements in the promoter region of CsHistones showed that CsHistones can respond to a variety of stresses. RNA-Seq analysis indicated that the expression of most CsHistones was associated with different stresses, including downy mildew, powdery mildew, wilt, heat, cold, salt stress, and waterlogging. Expression analysis showed that several genes of H3 group were highly expressed in different reproductive organs. Notably, CsCENH3 (CsHistone30) has the characteristics of a variant histone, and we demonstrated that CsCENH3 was localized on the nucleus and its proteins were expressed in centromere region. These findings provide valuable information for the identification and potential functions of Histone genes and ideas for the cultivation of CENH3-mediated haploid induction lines in cucumber.

CRISPR/Cas9 based mlo-mediated resistance against Podosphaera xanthii in cucumber (Cucumis sativus L.)

Powdery mildews (PM) are common and severe pathogen groups that threaten plants, and PM resistance is complex and polygenic in cucumbers. Previously mlo-based resistance was reported in various plants, including cucumber, with generated loss-of CsaMLO function mutants. However, mlo-based resistance in cucumber is also complex and involves additional mechanisms such as hypersensitive response (HR) and papillae formation. For this reason, we focused on determining the mlo-based powdery mildew resistance mechanism in cucumber. CRISPR/Cas9 was used in the present study to generate loss-of-function mutants for CsaMLO1, CsaMLO8, and CsaMLO11 of PM susceptible ADR27 cucumber inbred lines and CsaMLO mutants were obtained and validated. Trypan Blue and DAB staining were performed to detect Podosphaera xanthii germination/penetration rates and accumulation of Reactive Oxygen Species (ROS). Our results indicate that PM-susceptibility associated CsaMLOs in cucumber are negative regulators in different defense mechanisms against powdery mildew at early and late stages of infection. Further, the experiment results indicated that CsaMLO8 mutation-based resistance was associated with the pre-invasive response, while CsaMLO1 and CsaMLO11 could be negative regulators in the post-invasive defense response in cucumber against P. xanthii. Although the loss-of CsaMLO8 function confers the highest penetration resistance, CsaMLO1 and CsaMLO11 double mutations could be potential candidates for HR-based resistance against PM pathogen in cucumber. These results highlighted the crucial role of CRISPR/Cas9 to develop PM resistant cucumber cultivars, possessing strong pre-invasive defense with CsaMLO8 or post-invasive with CsaMLO1/CsaMLO11 mutations.

Identification of a major QTL, Parth6.1 associated with parthenocarpic fruit development in slicing cucumber genotype, Pusa Parthenocarpic Cucumber-6

Parthenocarpy is an extremely important trait that revolutionized the worldwide cultivation of cucumber under protected conditions. Pusa Parthenocarpic Cucumber-6 (PPC-6) is one of the important commercially cultivated varieties under protected conditions in India. Understanding the genetics of parthenocarpy, molecular mapping and the development of molecular markers closely associated with the trait will facilitate the introgression of parthenocarpic traits into non-conventional germplasm and elite varieties. The F1, F2 and back-crosses progenies with a non-parthenocarpic genotype, Pusa Uday indicated a single incomplete dominant gene controlling parthenocarpy in PPC-6. QTL-seq comprising of the early parthenocarpy and non-parthenocarpic bulks along with the parental lines identified two major genomic regions, one each in chromosome 3 and chromosome 6 spanning over a region of 2.7 Mb and 7.8 Mb, respectively. Conventional mapping using F2:3 population also identified two QTLs, Parth6.1 and Parth6.2 in chromosome 6 which indicated the presence of a major effect QTL in chromosome 6 determining parthenocarpy in PPC-6. The flanking markers, SSR01148 and SSR 01012 for Parth6.1 locus and SSR10476 and SSR 19174 for Parth6.2 locus were identified and can be used for introgression of parthenocarpy through the marker-assisted back-crossing programme. Functional annotation of the QTL-region identified two major genes, Csa_6G396640 and Csa_6G405890 designated as probable indole-3-pyruvate monooxygenase YUCCA11 and Auxin response factor 16, respectively associated with auxin biosynthesis as potential candidate genes. Csa_6G396640 showed only one insertion at position 2179 in the non-parthenocarpic parent. In the case of Csa_6G405890, more variations were observed between the two parents in the form of SNPs and InDels. The study provides insight about genomic regions, closely associated markers and possible candidate genes associated with parthenocarpy in PPC-6 which will be instrumental for functional genomics study and better understanding of parthenocarpy in cucumber.

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.

Research Advances in Genetic Mechanisms of Major Cucumber Diseases Resistance

Cucumber (Cucumis sativus L.) is an important economic vegetable crop worldwide that is susceptible to various common pathogens, including powdery mildew (PM), downy mildew (DM), and Fusarium wilt (FM). In cucumber breeding programs, identifying disease resistance and related molecular markers is generally a top priority. PM, DM, and FW are the major diseases of cucumber in China that cause severe yield losses and the genetic-based cucumber resistance against these diseases has been developed over the last decade. Still, the molecular mechanisms of cucumber disease resistance remain unclear. In this review, we summarize recent findings on the inheritance, molecular markers, and quantitative trait locus mapping of cucumber PM, DM, and FM resistance. In addition, several candidate genes, such as PM, DM, and FM resistance genes, with or without functional verification are reviewed. The data help to reveal the molecular mechanisms of cucumber disease resistance and provide exciting new opportunities for further resistance breeding.

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.

Construction of a collection of introgression lines of "Texas" almond DNA fragments in the "Earlygold" peach genetic background

Peach [Prunus persica L. Batsch] is one of the major temperate fruit tree species, the commercial materials of which have a low level of genetic variability. Almond [P. dulcis (Mill) DA Webb], a close relative of peach cultivated for its kernels, has a much higher level of diversity. The species are inter-compatible and often produce fertile hybrids, almond being a possible source of new genes for peach that could provide biotic and abiotic stress tolerance traits. In this paper we describe the development of a collection of peach-almond introgression lines (ILs) having a single fragment of almond (cv. Texas) in the peach background (cv. Earlygold). Lines with few introgressions were selected with markers from successive generations from a "Texas" × "Earlygold" F1 hybrid, initially using a set of SSRs and later with the 18 k peach SNP chip, allowing for the final extraction of 67 lines, 39 with almond heterozygous introgressions covering 99% of the genome, and 28 with homozygous introgressions covering 83% of the genome. As a proof of concept, four major genes and four quantitative characters were examined in the selected ILs giving results generally consistent with previous information on the genetics of these characters. This collection is the first of its kind produced in a woody perennial species and promises to be a valuable tool for genetic analyses, including dissection of quantitative traits, positional cloning, epistasis and as prebreeding material to introgress almond genes of interest into the peach commercial gene pool.

Genome-Wide Identification, Phylogenetic and Expression Pattern Analysis of GATA Family Genes in Cucumber (Cucumis sativus L.)

GATA transcription factors are a class of transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain, which play important roles in plant growth and development. The GATA gene family has been characterized in various plant species. However, GATA family genes have not been identified in cucumber. In this study, 26 GATA family genes were identified in cucumber genome, whose physicochemical characteristics, chromosomal distributions, phylogenetic tree, gene structures conserved motifs, cis-regulatory elements in promoters, homologous gene pairs, downstream target genes were analyzed. Tissue expression profiles of cucumber GATA family genes exhibited that 17 GATA genes showed constitutive expression, and some GATA genes showed tissue-specific expression patterns. RNA-seq analysis of green and virescent leaves revealed that seven GATA genes might be involved in the chloroplast development and chlorophyll biosynthesis. Importantly, expression patterns analysis of GATA genes in response to abiotic and biotic stresses indicated that some GATA genes respond to either abiotic stress or biotic stress, some GATA genes such as Csa2G162660, Csa3G017200, Csa3G165640, Csa4G646060, Csa5G622830 and Csa6G312540 were simultaneously functional in resistance to abiotic and biotic stresses. Overall, this study will provide useful information for further analysis of the biological functions of GATA factors in cucumber.

Transgressive resistance to Heterodera glycines in chromosome segment substitution lines derived from susceptible soybean parents

Chromosome segment substitution lines (CSSLs) are valuable genetic resources for quantitative trait loci (QTL) mapping of complex agronomic traits especially suitable for minor effect QTL. Here, 162 BC3 F7 -BC7 F3 CSSLs derived from crossing two susceptible parent lines, soybean [Glycine max (L.) Merr.] 'Suinong14' (recurrent parent) × wild soybean (G. soja Siebold & Zucc.) ZYD00006, were used for QTL mapping of soybean cyst nematode (SCN, Heterodera glycine Ichinohe) resistance based on female index (FI) and cysts per gram root (CGR) through phenotypic screening and whole-genome resequencing of CSSLs. Phenotypic results displayed a wide range of distribution and transgressive lines in both HG Type 2.5.7 FI and CGR and demonstrated a higher correlation between CGR and root weight (R2 = .5424) compared with than between FI and CGR (R2 = .0018). Using the single-marker analysis nonparametric mapping test, 33 significant QTL were detected on 18 chromosomes contributing resistance to FI and CGR. Fourteen QTL contributing 5.6-15.5% phenotypic variance (PVE) to FI were revealed on 11 chromosomes, and 16 QTL accounting for 6.1-36.2% PVE in CGR were detected on 14 chromosomes with strong additive effect by multiple-QTL model (MQM) mapping. Twenty-five and 13 out of all 38 QTL identified for FI and CGR on 20 chromosomes were from ZYD00006 and Suinong14, respectively. The CSSLs with the combination of positive alleles for FI, CGR, and root weight exhibited low nematode reproduction. For the first time, QTL associated with CGR have been detected, and both FI and CGR should be considered for breeding purposes in the absence of strong resistance genes such as rhg1 and Rhg4.

Transcriptome profiling reveals response genes for downy mildew resistance in cucumber

We discovered a potential defense pathway of cucumber to downy mildew. The signaling that activates the pathways of ROS and lignin accumulation may play an important role in the defense response. Many resistance genes were identified by transcriptome analysis. Downy mildew (DM), caused by Pseudoperonospora cubensis, is one of the most destructive diseases and causes severe yield losses of cucumber. However, the genes and pathways involved in regulating DM resistance were still poorly understood. In our study, we observed that the highly sensitive inbred line 53 (IL53) exhibited more severe disease symptoms than the highly resistant inbred line 51 (IL51) under P. cubensis infection. Furthermore, lignin, limiting the germination and extension of P. cubensis, and H₂O₂, as a signaling molecule during the resistant process, were both shown to increase, indicating that the signaling that activates these pathways might be responsible for the resistance divergence between IL51 and IL53. Transcriptome analysis, using the resistant and susceptible pools in F₂ populations with IL51 and IL53 as parents, showed that a series of differentially expressed genes was involved in multiple functions of defense response: pathogen-associated molecular pattern recognition, signal transduction, reactive oxygen species and lignin accumulation, and transcription regulators. Combining physiological data with transcriptomes, we predicted a potential molecular mechanism of cucumber resistance to DM. Our research provided a foundation for further studies on the mechanism of cucumber resistance to DM.

Connecting the dots: from nanodomains to physiological functions of REMORINs

REMORINs (REMs) are a plant-specific protein family, proposed regulators of membrane-associated molecular assemblies and well-established markers of plasma membrane nanodomains. REMs play a diverse set of functions in plant interactions with pathogens and symbionts, responses to abiotic stresses, hormone signaling and cell-to-cell communication. In this review, we highlight the established and more putative roles of REMs throughout the literature. We discuss the physiological functions of REMs, the mechanisms underlying their nanodomain-organization and their putative role as regulators of nanodomain-associated molecular assemblies. Furthermore, we discuss how REM phosphorylation may regulate their functional versatility. Overall, through data-mining and comparative analysis of the literature, we suggest how to further study the molecular mechanisms underpinning the functions of REMs.

Transcriptome profiling analysis reveals distinct resistance response of cucumber leaves infected with powdery mildew

Powdery mildew is the main disease affecting cucumber cultivation and causes severe economic loss. So far, research on cucumber resistance to powdery mildew has not yielded feasible solutions. This study selected two inbred cucumber lines, XY09-118 (resistant) and Q10 (susceptible) and investigated their responses to powdery mildew infection (harvested 24 and 48 h after inoculation) using RNA sequencing. More than 20,000 genes were detected in cucumber leaves both with and without powdery mildew infection at the above two time points. Among these, 5478 genes were identified as differently expressed genes (DEGs) between XY09-118 and Q10. Based on the databases GO and KEGG, the functions of DEGs were analysed. Moreover, the complex regulatory network for powdery mildew resistance was assessed, which involves plant hormone signal transduction, phenylpropanoid biosynthesis, plant-pathogen interaction and the MAPK signalling pathway. In particular, genes encoding WRKY, NAC and TCP were highlighted. In addition, genes involved in plant hormone biosynthesis, metabolism and signal transduction, pathogen resistance and abiotic stress response were analysed. Co-expression analysis indicated that the transcription factors correlated with plant hormone signal pathway and metabolism, defence and abiotic response. The expression of several genes was validated by qRT-PCR. The pathogen resistance regulatory network was identified by comparing resistant and susceptible inbred lines infected with powdery mildew. The transcriptome data provide novel insights into cucumber response to powdery mildew infection and the identified pathogen resistance genes will be highly useful for breeding efforts to enhance the resistance of cucumber to powdery mildew.

Genome-Wide Characterization of HSP90 Gene Family in Cucumber and Their Potential Roles in Response to Abiotic and Biotic Stresses

Heat shock protein 90 (HSP90) possesses critical functions in plant developmental control and defense reactions. The HSP90 gene family has been studied in various plant species. However, the HSP90 gene family in cucumber has not been characterized in detail. In this study, a total of six HSP90 genes were identified from the cucumber genome, which were distributed to five chromosomes. Phylogenetic analysis divided the cucumber HSP90 genes into two groups. The structural characteristics of cucumber HSP90 members in the same group were similar but varied among different groups. Synteny analysis showed that only one cucumber HSP90 gene, Csa1G569290, was conservative, which was not collinear with any HSP90 gene in Arabidopsis and rice. The other five cucumber HSP90 genes were collinear with five Arabidopsis HSP90 genes and six rice HSP90 genes. Only one pair of paralogous genes in the cucumber HSP90 gene family, namely one pair of tandem duplication genes (Csa1G569270/Csa1G569290), was detected. The promoter analysis showed that the promoters of cucumber HSP90 genes contained hormone, stress, and development-related cis-elements. Tissue-specific expression analysis revealed that only one cucumber HSP90 gene Csa3G183950 was highly expressed in tendril but low or not expressed in other tissues, while the other five HSP90 genes were expressed in all tissues. Furthermore, the expression levels of cucumber HSP90 genes were differentially induced by temperature and photoperiod, gibberellin (GA), downy mildew, and powdery mildew stimuli. Two cucumber HSP90 genes, Csa1G569270 and Csa1G569290, were both differentially expressed in response to abiotic and biotic stresses, which means that these two HSP90 genes play important roles in the process of cucumber growth and development. These findings improve our understanding of cucumber HSP90 family genes and provide preliminary information for further studies of cucumber HSP90 gene functions in plant growth and development.

QTL-seq analysis of powdery mildew resistance in a Korean cucumber inbred line

QTL mapping and RT-PCR analyses identified the CsGy5G015660 as a strong powdery mildew resistance candidate gene and natural variation of CsGy5G015660 allele was observed using 115 core germplasm. Powdery mildew (PM) is among the most serious fungal diseases encountered in the cultivation of cucurbits. The development of PM-resistant inbred lines is thus of considerable significance for cucumber breeding programs. In this study, we applied bulked segregant analysis combined with QTL-seq to identify PM resistance loci using F₂ population derived from a cross between two Korean cucumber inbred lines, PM-R (resistant) and PM-S (susceptible). Genome-wide SNP profiling using bulks of the two extreme phenotypes identified two QTLs on chromosomes 5 and 6, designated pm5.2 and pm6.1, respectively. The two PM resistance loci were validated using molecular marker-based classical QTL analysis: pm5.2 (30% R² at LOD 11) and pm6.1 (11% R² at LOD 3.2). Furthermore, reverse transcriptase-PCR analyses, using genes found to be polymorphic between PM-R and PM-S, were conducted to identify the candidate gene(s) responsible for PM resistance. We found that transcripts of the gene CsGy5G015660, encoding a putative leucine-rich repeat receptor-like serine/threonine-protein kinase (RPK2), showed specific accumulation in PM-R prior to the appearance of disease symptoms, and was accordingly considered a strong candidate gene for PM resistance. In addition, cleaved amplified polymorphic sequence markers from CsGy5G015660 were developed and used to screen 35 inbred lines. Natural variation in the CsGy5G015660 allele was also observed based on analysis of a core collection of 115 cucumber accessions. Our results provide new genetic insights for gaining a better understanding of the genetic basis of PM resistance in cucumber, and pave the way for further utilization in cucumber PM resistance breeding programs.

Comparative analysis of powdery mildew resistant and susceptible cultivated cucumber (Cucumis sativus L.) varieties to reveal the metabolic responses to Sphaerotheca fuliginea infection

BACKGROUND

Cucumber (Cucumis sativus L.) is a widely planted vegetable crop that suffers from various pathogen infections. Powdery mildew (PM) is typical disease caused by Sphaerotheca fuliginea infection and destroys the production of cucumber. However, the metabolic responses to S. fuliginea infection are largely unknown.

RESULTS

In our study, a PM resistant variety 'BK2' and a susceptible variety 'H136' were used to screen differentially accumulated metabolites (DAMs) and differentially expressed genes (DEGs) under S. fuliginea infection. Most of DEGs and DAMs were enriched in several primary and secondary metabolic pathways, including flavonoid, hormone, fatty acid and diterpenoid metabolisms. Our data showed that many flavonoid-related metabolites were significantly accumulated in BK2 rather than H136, suggesting an essential role of flavonoids in formation of resistant quality. Changes in expression of CYP73A, CYP81E1, CHS, F3H, HCT and F3'M genes provided a probable explanation for the differential accumulation of flavonoid-related metabolites. Interestingly, more hormone-related DEGs were detected in BK2 compared to H136, suggesting a violent response of hormone signaling pathways in the PM-resistant variety. The number of fatty acid metabolism-related DAMs in H136 was larger than that in BK2, indicating an active fatty acid metabolism in the PM-susceptible variety.

CONCLUSIONS

Many differentially expressed transcription factor genes were identified under S. fuliginea infection, providing some potential regulators for the improvement of PM resistance. PM resistance of cucumber was controlled by a complex network consisting of various hormonal and metabolic pathways.

Optimized combination methods for exploring and verifying disease-resistant transcription factors in melon

A large amount of omics data and number of bioinformatics tools has been produced. However, the methods for further exploring omics data are simple, in particular, to mine key regulatory genes, which are a priority concern in biological systems, and most of the specific functions are still unknown. First, raw data of two genotypes of melon (susceptible and resistant) were obtained by transcriptome analysis. Second, 391 transcription factors (TFs) were identified from the plant transcription factor database and cucurbit genomics database. Then, functional enrichment analysis indicated that these genes were mainly annotated in the process of transcription regulation. Third, 243 and 230 module-specific TFs were screened by weighted gene coexpression network analysis and short time series expression miner, respectively. Several TF genes, such as WRKYs and bHLHs, were regarded as key regulatory genes according to the values of significantly different modules. The coexpression network showed that these TF genes were significant correlated with resistance (R) genes, such as DRP2, RGA3, DRP1 and NB-ARC. Fourth, cis-acting element analysis illustrated that these R genes may bind to WRKY and bHLH. Finally, the expression of WRKY genes was verified by quantitative reverse transcription PCR (RT-qPCR). Phylogenetic analysis was carried out to further confirm that these TFs may play a critical role in Curcurbitaceae disease resistance. This study provides a new optimized combination strategy to explore the functions of TFs in a wide spectrum of biological processes. This strategy may also effectively predict potential relationships in the interactions of essential genes.

Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses

BACKGROUND

Cucumber (Cucumis sativus L.) is an economically important vegetable crop species. However, it is susceptible to various abiotic and biotic stresses. WRKY transcription factors play important roles in plant growth and development, particularly in the plant response to biotic and abiotic stresses. However, little is known about the expression pattern of WRKY genes under different stresses in cucumber.

RESULTS

In the present study, an analysis of the new assembly of the cucumber genome (v3.0) allowed the identification of 61 cucumber WRKY genes. Phylogenetic and synteny analyses were performed using related species to investigate the evolution of the cucumber WRKY genes. The 61 CsWRKYs were classified into three main groups, within which the gene structure and motif compositions were conserved. Tissue expression profiles of the WRKY genes demonstrated that 24 CsWRKY genes showed constitutive expression (FPKM > 1 in all samples), and some WRKY genes showed organ-specific expression, suggesting that these WRKYs might be important for plant growth and organ development in cucumber. Importantly, analysis of the CsWRKY gene expression patterns revealed that five CsWRKY genes strongly responded to both salt and heat stresses, 12 genes were observed to be expressed in response to infection from downy mildew and powdery mildew, and three CsWRKY genes simultaneously responded to all treatments analysed. Some CsWRKY genes were observed to be induced/repressed at different times after abiotic or biotic stress treatment, demonstrating that cucumber WRKY genes might play different roles during different stress responses and that their expression patterns vary in response to stresses.

CONCLUSIONS

Sixty-one WRKY genes were identified in cucumber, and insight into their classification, evolution, and expression patterns was gained in this study. Responses to different abiotic and biotic stresses in cucumber were also investigated. Our results provide a better understanding of the function of CsWRKY genes in improving abiotic and biotic stress resistance in cucumber.

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.

Identification of lncRNAs and their regulatory relationships with target genes and corresponding miRNAs in melon response to powdery mildew fungi

Melon (Cucumis melo L.), an economically beneficial crop widely cultivated around the world, is vulnerable to powdery mildew (PM). However, the studies on molecular mechanism of melon response to PM fungi is still limited. Long non coding RNAs (lncRNAs) have emerged as new regulators in plants response to biotic stresses. We predicted and identified the intricate regulatory roles of lncRNAs in melon response to PM fungi. A total of 539 lncRNAs were identified from PM-resistant (MR-1) and susceptible melon (Top Mark), in which 254 were significantly altered after PM fungi infection. Multiple target genes of lncRNAs were found to be involved in the hydrolysis of chitin, callose deposition and cell wall thickening, plant-pathogen interaction and plant hormone signal transduction pathway. Additionally, a total of 42 lncRNAs possess the various functions with microRNAs (miRNAs), including lncRNAs that are targeted by miRNAs and function as miRNA precursors or miRNA sponges. These findings provide a comprehensive view of potentially functional lncRNAs, corresponding target genes and related lncRNA-miRNA pairs, which will greatly increase our knowledge of the mechanism underlying susceptibility and resistance to PM in melon.

Comparative proteomic analysis of cucumber powdery mildew resistance between a single-segment substitution line and its recurrent parent

Powdery mildew (PM) is considered a major cause of yield losses and reduced quality in cucumber worldwide, but the molecular basis of PM resistance remains poorly understood. A segment substitution line, namely, SSL508-28, was developed with dominant PM resistance in the genetic background of PM-susceptible cucumber inbred line D8. The substituted segment contains 860 genes. An iTRAQ-based comparative proteomic technology was used to map the proteomes of PM-inoculated and untreated (control) D8 and SSL508-28. The number of differentially regulated proteins (DRPs) in SSL508-28 was almost three times higher than that in D8. Fourteen DRPs were located in the substituted segment interval. Comparative gene expression analysis revealed that nodulin-related protein 1 (NRP1) may be a good candidate for PM resistance. Gene Ontology enrichment analysis showed that DRPs functioning in tetrapyrrole biosynthetic process, sulfur metabolic process and cell redox homeostasis were specifically enriched in the resistant line SSL508-28. DRPs categorized in the KEGG term photosynthesis increased in both lines upon PM infection, suggesting that the strategies used by cucumber may be different from those used by other crops to react to PM attacks at the initial stage. The measurement of hydrogen peroxide and superoxide anion production and net photosynthetic rate were consistent with the changes in protein abundance, suggesting that the proteomic results were reliable. There was a poor correlation between DRPs measured by iTRAQ and the corresponding gene expression changes measured by RNA-seq with the same experimental design. Taken together, these findings improve the understanding of the molecular mechanisms underlying the response of cucumber to PM infection.

Development and use of chromosome segment substitution lines as a genetic resource for crop improvement

CSSLs are a complete library of introgression lines with chromosomal segments of usually a distant genotype in an adapted background and are valuable genetic resources for basic and applied research on improvement of complex traits. Chromosome segment substitution lines (CSSLs) are genetic stocks representing the complete genome of any genotype in the background of a cultivar as overlapping segments. Ideally, each CSSL has a single chromosome segment from the donor with a maximum recurrent parent genome recovered in the background. CSSL development program requires population-wide backcross breeding and genome-wide marker-assisted selection followed by selfing. Each line in a CSSL library has a specific marker-defined large donor segment. CSSLs are evaluated for any target phenotype to identify lines significantly different from the parental line. These CSSLs are then used to map quantitative trait loci (QTLs) or causal genes. CSSLs are valuable prebreeding tools for broadening the genetic base of existing cultivars and harnessing the genetic diversity from the wild- and distant-related species. These are resources for genetic map construction, mapping QTLs, genes or gene interactions and their functional analysis for crop improvement. In the last two decades, the utility of CSSLs in identification of novel genomic regions and QTL hot spots influencing a wide range of traits has been well demonstrated in food and commercial crops. This review presents an overview of how CSSLs are developed, their status in major crops and their use in genomic studies and gene discovery.

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

Cucurbit powdery mildew (PM) is one of the most severe fungal diseases, but the molecular mechanisms underlying PM resistance remain largely unknown, especially in pumpkin (Cucurbita moschata Duch.). The goal of this study was to identify gene expression differences in PM-treated plants (harvested at 24 h and 48 h after inoculation) and untreated (control) plants of inbred line “112–2” using RNA sequencing (RNA-Seq). The inbred line “112–2” has been purified over 8 consecutive generations of self-pollination and shows high resistance to PM. More than 7600 transcripts were examined in pumpkin leaves, and 3129 and 3080 differentially expressed genes (DEGs) were identified in inbred line “112–2” at 24 and 48 hours post inoculation (hpi), respectively. Based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database and GO (Gene Ontology) database, a complex regulatory network for PM resistance that may involve hormone signal transduction pathways, transcription factors and defense responses was revealed at the transcription level. In addition, the expression profiles of 16 selected genes were analyzed using quantitative RT-PCR. Among these genes, the transcript levels of 6 DEGs, including bHLH87 (Basic Helix-loop-helix transcription factor), ERF014 (Ethylene response factor), WRKY21 (WRKY domain), HSF (heat stress transcription factor A), MLO3 (Mildew Locus O), and SGT1 (Suppressor of G-Two Allele of Skp1), in PM-resistant “112–2” were found to be significantly up- or down-regulated both before 9 hpi and at 24 hpi or 48 hpi; this behavior differed from that observed in the PM-susceptible material (cultivar “Jiujiangjiaoding”). The transcriptome data provide novel insights into the response of Cucurbita moschata to PM stress and are expected to be highly useful for dissecting PM defense mechanisms in this major vegetable and for improving pumpkin breeding with enhanced resistance to PM.