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

GWAS reveals novel loci and identifies a pentatricopeptide repeat-containing protein (CsPPR) that improves low temperature germination in cucumber

Low temperatures (LTs) negatively affect the percentage and rate of cucumber (Cucumis sativus L.) seed germination, which has deleterious effects on yield. Here, a genome-wide association study (GWAS) was used to identify the genetic loci underlying low temperature germination (LTG) in 151 cucumber accessions that represented seven diverse ecotypes. Over two years, phenotypic data for LTG i.e., relative germination rate (RGR), relative germination energy (RGE), relative germination index (RGI) and relative radical length (RRL), were collected in two environments, and 17 of the 151 accessions were found to be highly cold tolerant using cluster analysis. A total of 1,522,847 significantly associated single-nucleotide polymorphism (SNP) were identified, and seven loci associated with LTG, on four chromosomes, were detected: gLTG1.1, gLTG1.2, gLTG1.3, gLTG4.1, gLTG5.1, gLTG5.2, and gLTG6.1 after resequencing of the accessions. Of the seven loci, three, i.e., gLTG1.2, gLTG4.1, and gLTG5.2, showed strong signals that were consistent over two years using the four germination indices, and are thus strong and stable for LTG. Eight candidate genes associated with abiotic stress were identified, and three of them were potentially causal to LTG: CsaV3_1G044080 (a pentatricopeptide repeat-containing protein) for gLTG1.2, CsaV3_4G013480 (a RING-type E3 ubiquitin transferase) for gLTG4.1, and CsaV3_5G029350 (a serine/threonine-protein kinase) for gLTG5.2. The function for CsPPR (CsaV3_1G044080) in regulating LTG was confirmed, as Arabidopsis lines ectopically expressing CsPPR showed higher germination and survival rates at 4°C compared to the wild-type, which preliminarily illustrates that CsPPR positively regulates cucumber cold tolerance at the germination stage. This study will provide insights into cucumber LT-tolerance mechanisms and further promote cucumber breeding development.

Sigma factor binding protein 1 (CsSIB1) is a putative candidate of the major-effect QTL dm5.3 for downy mildew resistance in cucumber (Cucumis sativus)

The dm5.3 major-effect QTL in cucumber encodes a homolog of Arabidopsis sigma factor binding protein 1 (CsSIB1). CsSIB1 positively regulates defense responses against downy mildew in cucumber through the salicylic acid (SA) biosynthesis/signaling pathway. Downy mildew (DM) caused by the oomycete pathogen Pseudoperonospora cubensis is an important disease of cucumber and other cucurbits. Our knowledge on molecular mechanisms of DM resistance is still limited. In this study, we reported identification and functional characterization of the candidate gene for the major-effect QTL, dm5.3 for DM resistance originated from PI 197088. The dm5.3 QTL was Modelized through marker-assisted development of near isogenic lines (NILs). NIL-derived segregating populations were used for fine mapping which narrowed the dm5.3 locus down to a 144 kb region. Based on multiple lines of evidence, we show that CsSIB1 (CsGy5G027140) that encodes the VQ motif-containing sigma factor binding protein 1 as the most likely candidate for dm5.3. Local association analysis identified a haplotype consisting of 7 SNPs inside the coding and promoter region of CsSIB1 that was associated with DM resistance. Expression of CsSIB1 was up-regulated with P. cubensis infection. Transcriptome profiling of NILs in response to P. cubensis inoculation revealed key players and associated gene networks in which increased expression of CsSIB1 antagonistically promoted salicylic acid (SA) but suppressed jasmonic acid (JA) biosynthesis/signaling pathways. Our work provides novel insights into the function of CsSIB1/dm5.3 as a disease resistance (R) gene. The roles of sigma factor binding protein genes in pathogen defense in cucumber were also discussed.

Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods

Downy mildew, commercially the most important disease of spinach, is caused by the obligate oomycete Peronospora effusa. In the past two decades, new pathogen races have repeatedly overcome the resistance used in newly released cultivars, urging the need for more durable resistance. Commercial spinach cultivars are bred with major R genes to impart resistance to downy mildew pathogens and are effective against some pathogen races/isolates. This work aimed to evaluate the worldwide USDA spinach germplasm collections and commercial cultivars for resistance to downy mildew pathogen in the field condition under natural inoculum pressure and conduct genome wide association analysis (GWAS) to identify resistance-associated genomic regions (alleles). Another objective was to evaluate the prediction accuracy (PA) using several genomic prediction (GP) methods to assess the potential implementation of genomic selection (GS) to improve spinach breeding for resistance to downy mildew pathogen. More than four hundred diverse spinach genotypes comprising USDA germplasm accessions and commercial cultivars were evaluated for resistance to downy mildew pathogen between 2017-2019 in Salinas Valley, California and Yuma, Arizona. GWAS was performed using single nucleotide polymorphism (SNP) markers identified via whole genome resequencing (WGR) in GAPIT and TASSEL programs; detected 14, 12, 5, and 10 significantly associated SNP markers with the resistance from four tested environments, respectively; and the QTL alleles were detected at the previously reported region of chromosome 3 in three of the four experiments. In parallel, PA was assessed using six GP models and seven unique marker datasets for field resistance to downy mildew pathogen across four tested environments. The results suggest the suitability of GS to improve field resistance to downy mildew pathogen. The QTL, SNP markers, and PA estimates provide new information in spinach breeding to select resistant plants and breeding lines through marker-assisted selection (MAS) and GS, eventually helping to accumulate beneficial alleles for durable disease resistance.

Identification and Expression Profile Analysis of the OSCA Gene Family Related to Abiotic and Biotic Stress Response in Cucumber

Simple Summary

Hyperosmolality-gated calcium-permeable channels (OSCAs) are calcium nonselective cation channel proteins involved in multiple biological processes. In this work, the members of the OSCA family in cucumber were systematically analyzed, including their sequence characteristics, phylogenetic relationships, conserved motifs, gene structures, promoter regions, and tissue expression patterns. In addition, the effects of different osmotic-related abiotic stresses [salt (NaCl), drought (PEG), and abscisic acid (ABA)] and three biotic stresses [powdery mildew (PM), downy mildew (DM), and root-knot nematode (RKN)] on OSCA family genes were also determined. The results indicated that cucumber OSCA genes play important roles in response to osmotic-related abiotic stresses and pathogen invasion. Overall, this study lays a foundation for research on the biological function and evolutionary process of OSCA family genes in cucumber.

Abstract

Calcium ions are important second messengers, playing an important role in the signal transduction pathways. Hyperosmolality gated calcium-permeable channels (OSCA) gene family members play critical modulating roles in response to osmotic-related abiotic stress as well as other abiotic and biotic stresses, which has been reported in many plant species such as Arabidopsis, rice, maize, and wheat. However, there has been no report about the identification and expression profile of the OSCA genes in cucumber. In this study, a total of nine OSCA genes were identified, which are unevenly distributed on the six chromosomes of cucumber. Phylogenetic analysis revealed that the OSCAs of cucumber, Arabidopsis, rice and maize were clustered into four clades. The motif arrangement of CsOSCAs was strongly conserved, and the CsOSCA genes in each group had similar genetic structure. A total of 11 and 10 types of cis-elements related to hormone and stress, respectively, were identified in the promoter regions of CsOSCA genes. Gene expression analysis indicated that the CsOSCA genes have different expression patterns in various tissues, and some of them were regulated by three osmotic-related abiotic stresses (salt, drought and ABA) and three biotic stresses (powdery mildew, downy mildew, and root-knot nematode). As the first genome-wide identification and characterization of the OSCA gene family in cucumber, this study lays a foundation for research on the biological function and evolutionary process of this gene family, which is of great significance for exploiting stress resistant cucumber varieties.

A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings

Salt stress seriously restricts plant growth and development, affects yield and quality, and thus becomes an urgent problem to be solved in cucumber stress resistance breeding. Mining salt tolerance genes and exploring the molecular mechanism of salt tolerance could accelerate the breeding of cucumber germplasm with excellent salt stress tolerance. In this study, 220 cucumber core accessions were used for Genome-Wide Association Studies (GWAS) and the identification of salt tolerance genes. The salinity injury index that was collected in two years showed significant differences among the core germplasm. A total of seven loci that were associated with salt tolerance in cucumber seedlings were repeatedly detected, which were located on Chr.2 (gST2.1), Chr.3 (gST3.1 and gST3.2), Chr.4 (gST4.1 and gST4.2), Chr.5 (gST5.1), and Chr.6 (gST6.1). Within these loci, 62 genes were analyzed, and 5 candidate genes (CsaV3_2G035120, CsaV3_3G023710, CsaV3_4G033150, CsaV3_5G023530, and CsaV3_6G009810) were predicted via the functional annotation of Arabidopsis homologous genes, haplotype of extreme salt-tolerant accessions, and qRT-PCR. These results provide a guide for further research on salt tolerance genes and molecular mechanisms of cucumber seedlings.

Genome-Wide Association Studies Reveal Candidate Genes Related to Stem Diameter in Cucumber (Cucumis sativus L.)

The stem diameter, an important agronomic trait, affects cucumber growth and yield. However, no genes responsible for cucumber stem diameter have been identified yet. In this study, the stem diameter of 88 cucumber core germplasms were measured in spring 2020, autumn 2020 and autumn 2021, and a genome-wide association study (GWAS) was carried out based on the gene sequence and stem diameter of core germplasms. A total of eight loci (gSD1.1, gSD2.1, gSD3.1, gSD3.2, gSD4.1, gSD5.1, gSD5.2, and gSD6.1) significantly associated with cucumber stem diameter were detected. Of these, five loci (gSD1.1, gSD2.1, gSD3.1, gSD5.2, and gSD6.1) were repeatedly detected in two or more seasons and were considered as robust and reliable loci. Based on the linkage disequilibrium sequences of the associated SNP loci, 37 genes were selected. By further investigating the five loci via analyzing Arabidopsis homologous genes and gene haplotypes, five genes (CsaV3_1G028310, CsaV3_2G006960, CsaV3_3G009560, CsaV3_5G031320, and CsaV3_6G031260) showed variations in amino acid sequence between thick stem lines and thin stem lines. Expression pattern analyses of these genes also showed a significant difference between thick stem and thin stem lines. This study laid the foundation for gene cloning and molecular mechanism study of cucumber stem development.

The amino acid permease (AAP) genes CsAAP2A and SlAAP5A/B are required for oomycete susceptibility in cucumber and tomato

Cucurbit downy mildew (DM), caused by the obligate biotroph Pseudoperonospora cubensis, is a destructive disease in cucumber. A valuable source of DM resistance is the Indian cucumber accession PI 197088, which harbours several quantitative trait loci (QTLs) contributing to quantitatively inherited DM resistance. With a combination of fine-mapping and transcriptomics, we identified Amino Acid Permease 2A (CsAAP2A) as a candidate gene for QTL DM4.1.3. Whole-genome and Sanger sequencing revealed the insertion of a Cucumis Mu-like element (CUMULE) transposon in the allele of the resistant near-isogenic line DM4.1.3. To confirm whether loss of CsAAP2A contributes to partial DM resistance, we performed targeting induced local lesions in genomes on a DM-susceptible cucumber genotype to identify an additional csaap2a mutant, which indeed was partially DM resistant. In view of the loss of the putative function as amino acid transporter, we measured amino acids in leaves. We found that DM-inoculated leaves of line DM4.1.3 (with the csaap2a mutation) contained significantly fewer amino acids than wild-type cucumber. The decreased flow of amino acids towards infected leaves in csaap2a plants compared to the wild type might explain the resistant phenotype of the mutant, as this would limit the available nutrients for the pathogen and thereby its fitness. To examine whether AAP genes play a conserved role as susceptibility factors in plant-oomycete interactions, we made targeted mutations in two AAP genes from tomato and studied the effect on susceptibility to Phytophthora infestans. We conclude that not only CsAAP2A but also SlAAP5A/SlAAP5B are susceptibility genes for oomycete pathogens.

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.