Complete resistance to powdery mildew and partial resistance to downy mildew in a Cucumis hystrix introgression line of cucumber were controlled by a co-localized locus

Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.)

KEY MESSAGE

Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

CsPM5.2, a phosphate transporter protein-like gene, promotes powdery mildew resistance in cucumber

Powdery mildew (PM) is one of the most serious fungal diseases affecting cucumbers (Cucumis sativus L.). The mechanism of PM resistance in cucumber is intricate and remains fragmentary as it is controlled by several genes. In this study, we detected the major-effect Quantitative Trait Locus (QTL), PM5.2, involved in PM resistance by QTL mapping. Through fine mapping, the dominant PM resistance gene, CsPM5.2, was cloned and its function was confirmed by transgenic complementation and natural variation identification. In cultivar 9930, a dysfunctional CsPM5.2 mutant resulted from a single nucleotide polymorphism in the coding region and endowed susceptibility to PM. CsPM5.2 encodes a phosphate transporter-like protein PHO1; H3. The expression of CsPM5.2 is ubiquitous and induced by the PM pathogen. In cucumber, both CsPM5.2 and Cspm5.1 (Csmlo1) are required for PM resistance. Transcriptome analysis suggested that the salicylic acid (SA) pathway may play an important role in CsPM5.2-mediated PM resistance. Our findings help parse the mechanisms of PM resistance and provide strategies for breeding PM-resistant cucumber cultivars.

CsMLO8/11 are required for full susceptibility of cucumber stem to powdery mildew and interact with CsCRK2 and CsRbohD

Powdery mildew (PM) is one of the most destructive diseases that threaten cucumber production globally. Efficient breeding of novel PM-resistant cultivars will require a robust understanding of the molecular mechanisms of cucumber resistance against PM. Using a genome-wide association study, we detected a locus significantly correlated with PM resistance in cucumber stem, pm-s5.1. A 1449-bp insertion in the CsMLO8 coding region at the pm-s5.1 locus resulted in enhanced stem PM resistance. Knockout mutants of CsMLO8 and CsMLO11 generated by CRISPR/Cas9 both showed improved PM resistance in the stem, hypocotyl, and leaves, and the double mutant mlo8mlo11 displayed even stronger resistance. We found that reactive oxygen species (ROS) accumulation was higher in the stem of these mutants. Protein interaction assays suggested that CsMLO8 and CsMLO11 could physically interact with CsRbohD and CsCRK2, respectively. Further, we showed that CsMLO8 and CsCRK2 competitively interact with the C-terminus of CsRbohD to affect CsCRK2-CsRbohD module-mediated ROS production during PM defense. These findings provide new insights into the understanding of CsMLO proteins during PM defense responses.

QTL Mapping for Ovary- and Fruit-Related Traits in Cucumis sativus-C. hystrix Introgression Line IL52

IL52 is a valuable introgression line obtained from interspecific hybridization between cultivated cucumber (Cucumis sativus L., 2n = 14) and the wild relative species C. hystrix Chakr. (2n = 24). IL52 exhibits high resistance to a number of diseases, including downy mildew, powdery mildew, and angular leaf spot. However, the ovary- and fruit-related traits of IL52 have not been thoroughly investigated. Here, we conducted quantitative trait loci (QTL) mapping for 11 traits related to ovary size, fruit size, and flowering time using a previously developed 155 F_7:8 RIL population derived from a cross between CCMC and IL52. In total, 27 QTL associated with the 11 traits were detected, distributed on seven chromosomes. These QTL explained 3.61% to 43.98% of the phenotypic variance. Notably, we identified a major-effect QTL (qOHN4.1) on chromosome 4 associated with the ovary hypanthium neck width and further delimited it into a 114-kb candidate region harboring 13 candidate genes. Furthermore, the QTL qOHN4.1 is co-localized with the QTL detected for ovary length, mature fruit length, and fruit neck length, all residing within the consensus QTL FS4.1, suggesting a plausible pleiotropic effect.

Pre- and Post-Zygotic Barriers Contribute to Reproductive Isolation and Correlate with Genetic Distance in Cucumis

Hybridization between Cucumis species, including cultivated melon (C. melo), is hampered by Interspecific Reproductive Barriers (IRBs). However, the nature of IRBs in Cucumis is largely unknown. This study explores locations, timing, and contribution to reproductive isolation (RI) of pre- and post-zygotic IRBs in Cucumis. To do this, we assessed crossability among Cucumis African wild species and C. melo at the pre-zygotic level by visualizing pollen tubes under fluorescence microscopy and, post-zygotically, by evaluating fruit/seed set and F₁ hybrid fertility. Genetic distances among Cucumis species were inferred from Genotyping-by-Sequencing, and its correlation with RI stages was analyzed. Observed pre- and post-zygotic IRBs included pollen tube arrest, fruit set failure, and hybrid male sterility. Unilateral cross-incongruity/incompatibility (UCI) was detected in some hybridizations, and dominant gene action is suggested for pistil-side UCI in interspecific F₁ hybrids. Notably, the allotetraploid C. ficifolius was very fertile as a seed parent but infertile in all reciprocal crosses. Contribution to RI was found significant for both pre- and post-zygotic IRBs. Additionally, a significant positive correlation was detected between genetic distance and pre- and post-zygotic RI stages. Interestingly, UCI offers an accessible system to dissect the genetics of IRBs in Cucumis, which may facilitate the use of wild relatives in breeding.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

A chitinase CsChi23 promoter polymorphism underlies cucumber resistance against Fusarium oxysporum f. sp. cucumerinum

Fusarium wilt disease, caused by Fusarium oxysporum f. sp. cucumerinum (Foc), leads to widespread yield loss and quality decline in cucumber. However, the molecular mechanisms underlying Foc resistance remain poorly understood. We report the mapping and functional characterisation of CsChi23, encoding a cucumber class I chitinase with antifungal properties. We assessed sequence variations at CsChi23 and the associated defence response against Foc. We functionally characterised CsChi23 using transgenic assay and expression analysis. The mechanism regulating CsChi23 expression was assessed by genetic and molecular approaches. CsChi23 was induced by Foc infection, which led to rapid upregulation in resistant cucumber lines. Overexpressing CsChi23 enhanced fusarium wilt resistance and reduced fungal biomass accumulation, whereas silencing CsChi23 causes loss of resistance. CsHB15, a homeodomain leucine zipper (HD-Zip) III transcription factor, was found to bind to the CsChi23 promoter region and activate its expression. Furthermore, silencing of CsHB15 reduces CsChi23 expression. A single-nucleotide polymorphism variation -400 bp upstream of CsChi23 abolished the HD-Zip III binding site in a susceptible cucumber line. Collectively, our study indicates that CsChi23 is sufficient to enhance fusarium wilt resistance and reveals a novel function of an HD-Zip III transcription factor in regulating chitinase expression in cucumber defence against fusarium wilt.

Characterization of the mitochondrial genome of Cucumis hystrix and comparison with other cucurbit crops

The mitochondria ofCucumis genus contain several intriguing features such as paternal inheritance and three-ring genome structure. However, the evolutionary relationships of mitochondria inCucumisremain elusive. Here, we assembled the mitochondrial genome ofC. hystrixand performed a comparative genomic analysis with other crops inthe Cucurbitaceae. The mitochondrial genome ofC. hystrixhas three circular-mapping chromosomes of lengths 1,113,461 bp, 110,683 bp, and 92,288 bp, which contain 73 genes including 38 protein-coding genes, 31tRNAgenes, and 4rRNAgenes. Repeat sequences, RNA editing, and horizontal gene transfer events were identified. The results of phylogenetic analyses, collinearity and gene clusters revealed thatC. hystrixis closer toC. sativus than to C. melo. Meanwhile, wedemonstrated mitochondrial paternal inheritance inC. hystrixbymolecular markers. In comparison with other cucurbitcrops, wefound amarker foridentification of germplasm resources ofCucumis. Collectively, our findings provide a tool to help clarify the paternal lineage within that genus in the evolution of Cucumis.

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.

Genetic Analysis of the Grapevine GATA Gene Family and Their Expression Profiles in Response to Hormone and Downy Mildew Infection

Grapevine (Vitis. vinifera L.) is one of the most economically important fruit crops throughout the world. However, grape production is increasingly impacted by numerous diseases, including downy mildew, caused by the oomycete Plasmopara viticola. In grapevine and other plants, members of the GATA family of transcription factors play key roles in light and phytohormone signaling. However, little is known about their potential roles in biotic defense responses. As a first step, we identified 27 GATA transcription factors in grapevine and defined their transcriptional responses to three biotic stress-related phytohormones (SA, MeJA, and BR) in callus cells, and challenge with P. viticola in a downy mildew-sensitive cultivar, V. vinifera ‘Pinot noir’, and a resistant cultivar, V. piasezkii ‘Liuba-8′. Many of the VvGATA genes had higher expression at 0.5 h after hormones treatments. Moreover, a group of VvGATAs was dramatically induced in ‘Liuba-8′ at 24 post infection by P. viticola. However, the same genes were significantly repressed and showed low expression levels in ‘Pinot noir’. Additionally, VvGATA27 was located in the nucleus and had transcriptional activity. Taken together, the study identified the GATA full gene families in grapes on phylogenetic analysis and protein structure. Moreover, this study provided a basis for discussing the roles of VvGATAs in response to hormones and P. viticola infection. Our results provided evidence for the selection of candidate genes against downy mildew and lay the foundation for further investigation of VvGATA transcription factors.

Characterization of lncRNAs and mRNAs Involved in Powdery Mildew Resistance in Cucumber

Powdery mildew (PM) is a severe fungal disease of cucumber worldwide. Identification of genetic factors resistant to PM is of great importance for marker-assisted breeding to ensure cucumber production. Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) have been shown to play important roles in plant development and immunity; however, whether they have a role in PM response in cucurbit crops remains unknown. We performed strand-specific RNA sequencing and miRNA sequencing using RNA from cucumber leaves of two near-isogenic lines (NILs), S1003 and NIL (Pm5.1) infected with PM, and systematically characterized the profiles of cucumber lncRNAs and messenger RNA (mRNAs) responsive to PM. In total, we identified 12,903 lncRNAs and 25,598 mRNAs responsive to PM. Differential expression (DE) analysis showed that 119 lncRNAs and 136 mRNAs correlated with PM resistance. Functional analysis of these DE lncRNAs and DE mRNAs revealed that they are significantly associated with phenylpropanoid biosynthesis, phenylalanine metabolism, ubiquinone and other terpenoid-quinone biosynthesis, and endocytosis. Particularly, two lncRNAs, LNC_006805 and LNC_012667, might play important roles in PM resistance. In addition, we also predicted mature miRNAs and competing endogenous RNA (ceRNA) networks of lncRNA-miRNA-mRNA involved in PM resistance. A total of 49 DE lncRNAs could potentially act as target mimics for 106 miRNAs. Taken together, our results provide an abundant resource for further exploration of cucumber lncRNAs, mRNAs, miRNAs, and ceRNAs in PM resistance, and will facilitate the molecular breeding for PM-resistant varieties to control this severe disease in cucumber.

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.

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 quantitative trait loci for powdery mildew resistance in highly resistant cucumber (Cucumis sativus L.) using ddRAD-seq analysis

Powdery mildew, caused by Podosphaera xanthii (syn. Sphaerotheca fuliginea ex Fr. Poll.), is one of the most economically important foliar diseases in cucumber (Cucumis sativus L.). Cucumber parental line 'Kyuri Chukanbohon Nou 5 Go', developed from weedy cucumber line CS-PMR1, is highly resistant to powdery mildew and is promising breeding material. We performed quantitative trait locus (QTL) analysis using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in a population from a cross between 'Kyuri Chukanbohon Nou 5 Go' and the Japanese native cultivar 'Kaga-aonaga-fushinari', which is susceptible to powdery mildew. The resistance of the population and its parents was evaluated using leaf disc assays and image analysis. We detected one major QTL on Chr. 5 that was effective at both 20°C and 25°C and one minor QTL on Chr. 1 effective at 20°C. We detected two additional QTLs in subpopulation: one on Chr. 3 effective at 20°C and one on Chr. 5 effective at both 20°C and 25°C in a position different from the major QTL. The resistance alleles at all four QTLs were contributed by 'Kyuri Chukanbohon Nou 5 Go'. The results of this study can be used to develop practical DNA markers tightly linked to genes for powdery mildew resistance.

Whole-Genome Sequence of Synthesized Allopolyploids in Cucumis Reveals Insights into the Genome Evolution of Allopolyploidization

The importance of allopolyploidy in plant evolution has been widely recognized. The genetic changes triggered by allopolyploidy, however, are not yet fully understood due to inconsistent phenomena reported across diverse species. The construction of synthetic polyploids offers a controlled approach to systematically reveal genomic changes that occur during the process of polyploidy. This study reports the first fully sequenced synthetic allopolyploid constructed from a cross between Cucumis sativus and C. hystrix, with high-quality assembly. The two subgenomes are confidently partitioned and the C. sativus-originated subgenome predominates over the C. hystrix-originated subgenome, retaining more sequences and showing higher homeologous gene expression. Most of the genomic changes emerge immediately after interspecific hybridization. Analysis of a series of genome sequences from several generations (S0, S4-S13) of C. ×hytivus confirms that genomic changes occurred in the very first generations, subsequently slowing down as the process of diploidization is initiated. The duplicated genome of the allopolyploid with double genes from both parents broadens the genetic base of C. ×hytivus, resulting in enhanced phenotypic plasticity. This study provides novel insights into plant polyploid genome evolution and demonstrates a promising strategy for the development of a wide array of novel plant species and varieties through artificial polyploidization.

Identification of Novel Loci and Candidate Genes for Resistance to Powdery Mildew in a Resequenced Cucumber Germplasm

Powdery mildew (PM) is one of the most serious diseases in cucumber and causes huge yield loss. Multiple quantitative trait loci (QTLs) for PM resistance have been reported in previous studies using a limited number of cucumber accessions. In this study, a cucumber core germplasm (CG) consisting of 94 resequenced lines was evaluated for PM resistance in four trials across three years (2013, 2014, and 2016). These trials were performed on adult plants in the field with natural infection. Using genome-wide association study (GWAS), 13 loci (pmG1.1, pmG1.2, pmG2.1, pmG2.2, pmG3.1, pmG4.1, pmG4.2, pmG5.1, pmG5.2, pmG5.3, pmG5.4, pmG6.1, and pmG6.2) associated with PM resistance were detected on all chromosomes except for Chr.7. Among these loci, ten were mapped to chromosomal intervals where QTLs had been reported in previous studies, while, three (pmG2.1, pmG3.1, and pmG4.1) were novel. The loci of pmG2.1, pmG5.2, pmG5.3 showed stronger signal in four trials. Based on the annotation of homologous genes in Arabidopsis and pairwise LD correlation analysis, candidate genes located in the QTL intervals were predicted. SNPs in these candidate genes were analyzed between haplotypes of highly resistant (HR) and susceptible (HS) CG lines, which were defined based on combing disease index data of all trials. Furthermore, candidate genes (Csa5G622830 and CsGy5G015660) reported in previous studies for PM resistance and cucumber orthologues of several PM susceptibility (S) genes (PMR5, PMR-6, and MLO) that are colocalized with certain QTLs, were analyzed for their potential contribution to the QTL effect on both PM and DM in the CG population. This study shows that the CG germplasm is a very valuable resource carrying known and novel QTLs for both PM and DM resistance, which can be exploited in cucumber breeding.

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.

Global Profiling of lncRNAs Expression Responsive to Allopolyploidization in Cucumis

Long non-coding RNAs (lncRNAs) play critical regulatory roles in various biological processes. However, the presence of lncRNAs and how they function in plant polyploidy are still largely unknown. Hence, we examined the profile of lncRNAs in a nascent allotetraploid Cucumis hytivus (S₁₄), its diploid parents, and the F₁ hybrid, to reveal the function of lncRNAs in plant-interspecific hybridization and whole genome duplication. Results showed that 2206 lncRNAs evenly transcribed from all 19 chromosomes were identified in C. hytivus, 44.6% of which were from intergenic regions. Based on the expression trend in allopolyploidization, we found that a high proportion of lncRNAs (94.6%) showed up-regulated expression to varying degrees following hybridization. However, few lncRNAs (33, 2.1%) were non-additively expressed after genome duplication, suggesting the significant effect of hybridization on lncRNAs, rather than genome duplication. Furthermore, 253 cis-regulated target genes were predicted for these differentially expressed lncRNAs in S₁₄, which mainly participated in chloroplast biological regulation (e.g., chlorophyll synthesis and light harvesting system). Overall, this study provides new insight into the function of lncRNAs during the processes of hybridization and polyploidization in plant evolution.

Identification of Novel Loci and Candidate Genes for Cucumber Downy Mildew Resistance Using GWAS

Downy mildew (DM) is one of the most serious diseases in cucumber. Multiple quantitative trait loci (QTLs) for DM resistance have been detected in a limited number of cucumber accessions. In this study we applied genome-wide association analysis (GWAS) to detected genetic loci for DM resistance in a core germplasm (CG) of cucumber lines that represent diverse origins and ecotypes. Phenotypic data on responses to DM infection were collected in four field trials across three years, 2014, 2015, and 2016. With the resequencing data of these CG lines, GWAS for DM resistance was performed and detected 18 loci that were distributed on all the seven cucumber chromosomes. Of these 18 loci, only six (dmG1.4, dmG4.1, dmG4.3, dmG5.2, dmG7.1, and dmG7.2) were detected in two experiments, and were considered as loci with a stable effect on DM resistance. Further, 16 out of the 18 loci colocalized with the QTLs that were reported in previous studies and two loci, dmG2.1 and dmG7.1, were novel ones identified only in this study. Based on the annotation of homologous genes in Arabidopsis and pairwise LD correlation analysis, several candidate genes were identified as potential causal genes underlying the stable and novel loci, including Csa1G575030 for dmG1.4, Csa2G060360 for dmG2.1, Csa4G064680 for dmG4.1, Csa5G606470 for dmG5.2, and Csa7G004020 for dmG7.1. This study shows that the CG germplasm is a very valuable resource carrying known and novel QTLs for DM resistance. The potential of using these CG lines for future allele-mining of candidate genes was discussed in the context of breeding cucumber with resistance to DM.

Fine Mapping and Transcriptome Analysis of Virescent Leaf Gene v-2 in Cucumber (Cucumis sativus L.)

Leaf color mutants are the ideal materials to explore the pathways of chlorophyll metabolism, chloroplast development and photosynthesis system. In this study, a new virescent leaf mutant 104Y was identified by spontaneous mutation, whose cotyledon and upper five true leaves were yellow color. The yellow true leaves gradually turned green from top to bottom with increased chlorophyll contents. Genetic analysis indicated that the virescent leaf was controlled by one single recessive gene v-2, which was accurately mapped into 36.0-39.7 Mb interval on chromosome 3 by using BSA-seq and linkage analysis. Fine mapping analysis further narrowed v-2 into 73-kb genomic region including eight genes with BC₁ and F₂ populations. Through BSA-seq and cDNA sequencing analysis, only one nonsynonymous mutation existed in the Csa3G890020 gene encoding auxin F-box protein was identified, which was predicted as the candidate gene controlling virescent leaf. Comparative transcriptome analysis and quantitative real-time PCR analysis revealed that the expression level of Csa3G890020 was not changed between EC1 and 104Y. However, RNA-seq analysis identified that the key genes involved in chlorophyll biosynthesis and auxin signaling transduction network were mainly down-regulated in 104Y compared with EC1, which indicated that the regulatory functions of Csa3G890020 could be performed at post-transcriptional level rather than transcriptional level. This is the first report to map-based clone an auxin F-box protein gene related to virescent leaf in cucumber. The results will exhibit a new insight into the chlorophyll biosynthesis regulated by auxin signaling transduction network.

Development of alien addition lines from Cucumis hystrix in Cucumis sativus: cytological and molecular marker analyses

Transferring desired genes from wild species to cultivars through alien addition lines (AALs) has been shown to be an effective method for genetic improvement. Cucumis hystrix Chakr. (HH, 2n = 24) is a wild species of Cucumis that possesses many resistant genes. A synthetic allotetraploid species, C. hytivus (HHCC, 2n = 38), was obtained from the cross between cultivated cucumber, C. sativus (CC, 2n = 14), and C. hystrix followed by chromosome doubling. Cucumis sativus - C. hystrix AALs were developed by continuous backcrossing to the cultivated cucumbers. In this study, 10 different types of AALs (CC-H01, CC-H06, CC-H08, CC-H10, CC-H12, CC-H06+H09, CC-H06+H10, CC-H06+H12, CC-H08+H10, CC-H01+H06+H10) were identified based on the analysis of fluorescence in situ hybridization (FISH) and molecular markers specific to C. hystrix chromosomes. And the behavior of the alien chromosomes in three AALs (CC-H01, CC-H06+H10, CC-H01+H06+H10) at meiosis was investigated. The results showed that alien chromosomes paired with C. sativus chromosome in few pollen mother cells (PMCs). Further, disomic alien addition lines (DAALs) carrying a pair of C. hystrix chromosome H10 were screened from the selfed progenies of CC-H10. Chromosome pairing between genomes provides cytological evidence for the possible introgression of alien chromosome segments. The development of AALs could serve as a key step for exploiting and utilizing valuable genes from C. hystrix.

PyBSASeq: a simple and effective algorithm for bulked segregant analysis with whole-genome sequencing data

BACKGROUND

Bulked segregant analysis (BSA), coupled with next-generation sequencing, allows the rapid identification of both qualitative and quantitative trait loci (QTL), and this technique is referred to as BSA-Seq here. The current SNP index method and G-statistic method for BSA-Seq data analysis require relatively high sequencing coverage to detect significant single nucleotide polymorphism (SNP)-trait associations, which leads to high sequencing cost.

RESULTS

We developed a simple and effective algorithm for BSA-Seq data analysis and implemented it in Python; the program was named PyBSASeq. Using PyBSASeq, the significant SNPs (sSNPs), SNPs likely associated with the trait, were identified via Fisher's exact test, and then the ratio of the sSNPs to total SNPs in a chromosomal interval was used to detect the genomic regions that condition the trait of interest. The results obtained this way are similar to those generated via the current methods, but with more than five times higher sensitivity. This approach was termed the significant SNP method here.

CONCLUSIONS

The significant SNP method allows the detection of SNP-trait associations at much lower sequencing coverage than the current methods, leading to ~ 80% lower sequencing cost and making BSA-Seq more accessible to the research community and more applicable to the species with a large genome.

An irregularly striped rind mutant reveals new insight into the function of PG1β in cucumber (Cucumis sativus L.)

Via bulked segregant analysis sequencing combined with linkage mapping, the ist gene responsible for the irregularly striped rind mutation was delimited to a 144-kb region in cucumber. Sequencing and expression analysis identified Csa1G005490 as the candidate gene. The rind appearance of cucumber is one of the most important commercial quality traits. Usually, an immature cucumber fruit has a uniform rind that varies from green to yellow to white among different cultivated varieties. In the present paper, we isolated a novel fruit appearance cucumber mutant, ist, that has an irregularly striped rind pattern. The mutant displayed green irregular stripes on a yellow-green background at the immature fruit stage. Genetic analysis revealed that a single recessive gene, ist, is responsible for this mutation. A BSA (bulked segregant analysis) sequencing approach combined with genetic mapping delimited the ist locus to an interval with a length of 144 kb, and 21 predicted genes were annotated in the region. Based on mutation site screening and expression analysis, two single-nucleotide polymorphisms within the candidate gene, Csa1G005490, were identified as constituting the mutation. Csa1G005490 encodes a polygalacturonase-1 noncatalytic subunit beta protein (PG1β) known to be involved in fruit softening. The expression of Csa1G005490 was significantly lower in the ist mutant than in the wild type. Transcriptome analysis identified 1796 differentially expressed genes (DEGs) between the ist mutant and wild type. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that these DEGs were enriched mostly in photosynthesis and chlorophyll metabolism pathways. Decreased expression patterns of several chlorophyll synthesis genes in the mutant suggest that ist plays a key role in chlorophyll biosynthesis. These results will provide new insight into the molecular mechanism underlying rind appearance polymorphisms in cucumber.

Molecularly tagged genes and quantitative trait loci in cucumber with recommendations for QTL nomenclature

Cucumber, Cucumis sativus L. (2n = 2x = 14), is an important vegetable crop worldwide. It was the first specialty crop with a publicly available draft genome. Its relatively small, diploid genome, short life cycle, and self-compatible mating system offers advantages for genetic studies. In recent years, significant progress has been made in molecular mapping, and identification of genes and QTL responsible for key phenotypic traits, but a systematic review of the work is lacking. Here, we conducted an extensive literature review on mutants, genes and QTL that have been molecularly mapped or characterized in cucumber. We documented 81 simply inherited trait genes or major-effect QTL that have been cloned or fine mapped. For each gene, detailed information was compiled including chromosome locations, allelic variants and associated polymorphisms, predicted functions, and diagnostic markers that could be used for marker-assisted selection in cucumber breeding. We also documented 322 QTL for 42 quantitative traits, including 109 for disease resistances against seven pathogens. By alignment of these QTL on the latest version of cucumber draft genomes, consensus QTL across multiple studies were inferred, which provided insights into heritable correlations among different traits. Through collaborative efforts among public and private cucumber researchers, we identified 130 quantitative traits and developed a set of recommendations for QTL nomenclature in cucumber. This is the first attempt to systematically summarize, analyze and inventory cucumber mutants, cloned or mapped genes and QTL, which should be a useful resource for the cucurbit research community.

Inheritance and Quantitative Trait Locus Mapping of Fusarium Wilt Resistance in Cucumber

Fusarium wilt (FW) is a very serious soil-borne disease worldwide, which usually results in huge yield losses in cucumber production. However, the inheritance and molecular mechanism of the response to FW are still unknown in cucumber (Cucumis sativus L.). In this study, two inbred cucumber lines Superina (P₁) and Rijiecheng (P₂) were used as the sensitive and resistant lines, respectively. A mixed major gene plus polygene inheritance model was used to analyze the resistance to FW in different generations of cucumber, namely, P₁, P₂, F₁ (P₁×P₂), B₁, and B₂, obtained by backcrossing F₁ plants with Superina (B₁) or Rijiecheng (B₂), and F₂, obtained by self-crossing the F₁ plants. After screening 18 genetic models, we chose the E-1 model, which included two pairs of additive-dominance-epistatic major genes and additive-dominance polygenes, as the optimal model for resistance to FW on the basis of fitness tests. The major effect quantitative trait locus (QTL) fw2.1 was detected in a 1.91-Mb-long region of chromosome 2 by bulked-segregant analysis. We used five insertion/deletion markers to fine-map the fw2.1 to a 0.60 Mb interval from 1,248,093 to 1,817,308 bp on chromosome 2 that contained 80 candidate genes. We also used the transcriptome data of Rijiecheng inoculated with Fusarium oxysporum f. sp. cucumerinum (Foc) to screen the candidate genes. Twelve differentially expressed genes were detected in fw2.1, and five of them were significantly induced by FW. The expression levels of the five genes were higher in FW-resistant Rijiecheng inoculated with Foc than in the control inoculated with water. Our results will contribute to a better understanding of the genetic basis of FW resistance in cucumber, which may help in breeding FW-resistant cucumber lines in the future.