Genetic characterization of cucumber genetic resources in the NARO Genebank indicates their multiple dispersal trajectories to the East.

KEY MESSAGE: Genotyping-by-sequencing of 723 worldwide cucumber genetic resources revealed that cucumbers were dispersed eastward via at least three distinct routes, one to Southeast Asia and two from different directions to East Asia. The cucumber (Cucumis sativus) is an economically important vegetable crop cultivated and consumed worldwide. Despite its popularity, the manner in which cucumbers were dispersed from their origin in South Asia to the rest of the world, particularly to the east, remains a mystery due to the lack of written records. In this study, we performed genotyping-by-sequencing (GBS) on 723 worldwide cucumber accessions, mainly deposited in the Japanese National Agriculture and Food Research Organization (NARO) Genebank, to characterize their genetic diversity, relationships, and population structure. Analyses based on over 60,000 genome-wide single-nucleotide polymorphisms identified by GBS revealed clear genetic differentiation between Southeast and East Asian populations, suggesting that they reached their respective region independently, not progressively. A deeper investigation of the East Asian population identified two subpopulations with different fruit characteristics, supporting the traditional classification of East Asian cucumbers into two types thought to have been introduced by independent routes. Finally, we developed a core collection of 100 accessions representing at least 93.2% of the genetic diversity present in the entire collection. The genetic relationships and population structure, their associations with geographic distribution and phenotypic traits, and the core collection presented in this study are valuable resources for elucidating the dispersal history and promoting the efficient use and management of genetic resources for research and breeding in cucumber.

Melon diversity on the Silk Road by molecular phylogenetic analysis in Kazakhstan melons.

To uncover population structure, phylogenetic relationship, and diversity in melons along the famous Silk Road, a seed size measurement and a phylogenetic analysis using five chloroplast genome markers, 17 RAPD markers and 11 SSR markers were conducted for 87 Kazakh melon accessions with reference accessions. Kazakh melon accessions had large seed with exception of two accessions of weedy melon, Group Agrestis, and consisted of three cytoplasm types, of which Ib-1/-2 and Ib-3 were dominant in Kazakhstan and nearby areas such as northwestern China, Central Asia and Russia. Molecular phylogeny showed that two unique genetic groups, STIa-2 with Ib-1/-2 cytoplasm and STIa-1 with Ib-3 cytoplasm, and one admixed group, STIAD combined with STIa and STIb, were prevalent across all Kazakh melon groups. STIAD melons that phylogenetically overlapped with STIa-1 and STIa-2 melons were frequent in the eastern Silk Road region, including Kazakhstan. Evidently, a small population contributed to melon development and variation in the eastern Silk Road. Conscious preservation of fruit traits specific to Kazakh melon groups is thought to play a role in the conservation of Kazakh melon genetic variation during melon production, where hybrid progenies were generated through open pollination.

Elucidation of genetic variation and population structure of melon genetic resources in the NARO Genebank, and construction of the World Melon Core Collection.

Numerous genetic resources of major crops have been introduced from around the world and deposited in Japanese National Agriculture and Food Research Organization (NARO) Genebank. Understanding their genetic variation and selecting a representative subset ("core collection") are essential for optimal management and efficient use of genetic resources. In this study, we conducted genotyping-by-sequencing (GBS) to characterize the genetic relationships and population structure in 755 accessions of melon genetic resources. The GBS identified 39,324 single-nucleotide polymorphisms (SNPs) that are distributed throughout the melon genome with high density (one SNP/10.6 kb). The phylogenetic relationships and population structure inferred using this SNP dataset are highly associated with the cytoplasm type and geographical origin. Our results strongly support the recent hypothesis that cultivated melon was established in Africa and India through multiple independent domestication events. Finally, we constructed a World Melon Core Collection that covers at least 82% of the genetic diversity and has a wide range of geographical origins and fruit morphology. The genome-wide SNP dataset, phylogenetic relationships, population structure, and the core collection provided in this study should largely contribute to genetic research, breeding, and genetic resource preservation in melon.

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.

Host range and complete genome sequence of Cucurbit chlorotic yellows virus, a new member of the genus Crinivirus.

Cucurbit chlorotic yellows virus (CCYV) causes chlorotic yellows on cucumber (Cucumis sativus) and melon (Cucumis melo) and is transmitted by Bemisia tabaci biotype B and Q whiteflies. To characterize the host range of CCYV, 21 cucurbitaceous and 12 other plant species were inoculated using whitefly vectors. All tested Cucumis spp. except Cucumis anguria and Cucumis zeyheri were systemically infected with CCYV, although infection rates varied among species. Citrullus lanatus, Cucurbita pepo, and Luffa cylindrica were susceptible to CCYV; however, the infection rates were low and symptoms were unclear. In addition to the cucurbitaceous plants, Beta vulgaris, Chenopodium amaranticolor, Chenopodium quinoa, Spinacia oleracea, Lactuca sativa, Datura stramonium, and Nicotiana benthamiana were also systemically infected by CCYV. Complete RNA1 and RNA2 were reverse-transcribed, cloned, and sequenced. CCYV RNA1 was found to be 8,607 nucleotides (nt) long and contained four open reading frames (ORFs). The first ORF spanned methyltransferase and RNA helicase domains followed by an RNA-dependent RNA polymerase domain, presumably translated by a +1 ribosomal frameshift. CCYV RNA2 was found to be 8,041 nt long and contained eight ORFs, including the hallmark gene array of the family Closteroviridae. Phylogenetic analysis demonstrated that CCYV was genetically close to Lettuce chlorosis virus, Bean yellow disorder virus, and Cucurbit yellow stunting disorder virus. Amino acid sequence similarities of representative proteins with these viruses indicated that CCYV should be classified as a distinct crinivirus species.

Identification of QTLs for resistance to powdery mildew and SSR markers diagnostic for powdery mildew resistance genes in melon (Cucumis melo L.).

Powdery mildew caused by Podosphaera xanthii is an important foliar disease in melon. To find molecular markers for marker-assisted selection, we constructed a genetic linkage map of melon based on a population of 93 recombinant inbred lines derived from crosses between highly resistant AR 5 and susceptible 'Earl's Favourite (Harukei 3)'. The map spans 877 cM and consists of 167 markers, comprising 157 simple sequence repeats (SSRs), 7 sequence characterized amplified region/cleavage amplified polymorphic sequence markers and 3 phenotypic markers segregating into 20 linkage groups. Among them, 37 SSRs and 6 other markers were common to previous maps. Quantitative trait locus (QTL) analysis identified two loci for resistance to powdery mildew. The effects of these QTLs varied depending on strain and plant stage. The percentage of phenotypic variance explained for resistance to the pxA strain was similar between QTLs (R (2) = 22-28%). For resistance to pxB strain, the QTL on linkage group (LG) XII was responsible for much more of the variance (41-46%) than that on LG IIA (12-13%). The QTL on LG IIA was located between two SSR markers. Using an independent population, we demonstrated the effectiveness of these markers. This is the first report of universal and effective markers linked to a gene for powdery mildew resistance in melon.