Genetic mapping and QTL analysis of horticultural traits in cucumber ( Cucumis sativus L.) using recombinant inbred lines

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.

Targeted Identification of Rice Grain-Associated Gene Allelic Variation Through Mutation Induction, Targeted Sequencing, and Whole Genome Sequencing Combined with a Mixed-Samples Strategy

BACKGROUND

The mining of new allelic variation and the induction of new genetic variability are the basis for improving breeding efficiency.

RESULTS

In this study, in total, 3872 heavy ion-irradiated M₂ generation rice seeds and individual leaves were collected. The grain length was between 8 and 10.22 mm. The grain width was between 1.54 and 2.87 mm. The results showed that there was extensive variation in granulotype. The allelic variation in GS3 and GW5 was detected in 484 mixed samples (8:1) using targeted sequencing technology, and 12 mixed samples containing potential mutations and 15 SNPs were obtained; combined with Sanger sequencing and phenotype data, 13 key mutants and their corresponding SNPs were obtained; protein structural and functional analysis of key mutants screened out 6 allelic variants leading to altered grain shape, as well as the corresponding mutants, including long-grain mutants GS3-2 and GS3-7, short-grain mutants GS3-3 and GS3-5, wide-grain mutant GW5-1 and narrow-grain mutant GW5-4; whole genome sequencing identified new grain length gene allelic variants GS3-G1, GS3-G2 and GS3-G3.

CONCLUSION

Based on the above studies, we found 6 granulotype mutants and 9 granulotype-related allelic variants, which provided new functional gene loci and a material basis for molecular breeding and genotype mutation and phenotype analysis. We propose a method for targeted identification of allelic variation in rice grain type genes by combining targeted sequencing of mixed samples and whole genome sequencing. The method has the characteristics of low detection cost, short detection period, and flexible detection of traits and genes.

The mutation of C-24 reductase, a key enzyme involved in brassinolide biosynthesis, confers a novel compact plant architecture phenotype to cucumber

A novel compact plant architecture mutant, cpa-2, was identified from EMS-induced mutagenesis. Bulked segregant analysis sequencing and map-based cloning revealed CsDWF1 encoding C-24 reductase enzyme as the candidate gene. The compact architecture is a vital and valuable agronomic trait that helps to reduce the labor of plant management, and improve the fruit yield by increasing planting density in cucumbers. However, the molecular basis underlying the regulation of plant architecture in cucumber is complex and largely unknown. In this study, a novel recessive compact allele, designated as cpa-2 (compact plant architecture-2) was fine mapped in a 109 kb region on chromosome 7 by the strategy of bulked segregant analysis sequencing combined with map-based cloning. Gene annotation of the corresponding region revealed that the CsaV3_7G030530 (CsDWF1) gene encoding C-24 reductase, which acts as the key enzyme in brassinosteroids biosynthesis, functions as the candidate gene for cpa-2. Sequence analysis showed that a single-nucleotide mutation (G to A) in the second exon of CsaV3_7G030530 caused an amino acid substitution from E⁵⁰² to K⁵⁰². Compared with wild-type CCMC, CsDWF1 had lower expression levels in the stem, leaf and ovary of cpa-2. In addition, the compact phenotype in cpa-2 could be partially restored by exogenous BR application. Transcriptome analysis revealed that many genes related to plant growth hormones were differentially expressed in cpa-2 plants. This is the first report about the characterization and cloning of the CsDWF1 gene. This work revealed the importance of CsDWF1 in plant development regulation and extended our understanding of the interaction between BRs and other hormones for plant architecture development.

Genetic analysis and QTL mapping of yield and fruit traits in bitter gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an economically important vegetable crop grown in tropical parts of the world. In this study, a high-density linkage map of M. charantia was constructed through genotyping-by-sequencing (GBS) technology using F2:3 mapping population generated from the cross DBGy-201 × Pusa Do Mausami. About 2013 high-quality SNPs were assigned on a total of 20 linkage groups (LGs) spanning over 2329.2 CM with an average genetic distance of 1.16 CM. QTL analysis was performed for six major yield-contributing traits such as fruit length, fruit diameter, fruit weight, fruit flesh thickness, number of fruits per plant and yield per plant. These six quantitative traits were mapped with 19 QTLs (9 QTLs with LOD > 3) using composite interval mapping (CIM). Among 19 QTLs, 12 QTLs derived from 'Pusa Do Mausami' revealed a negative additive effect when its allele increased trait score whereas 7 QTLs derived from 'DBGy-201' revealed a positive additive effect when its allele trait score increased. The phenotypic variation (R2%) elucidated by these QTLs ranged from 0.09% (fruit flesh thickness) on LG 14 to 32.65% (fruit diameter) on LG 16 and a total of six major QTLs detected. Most QTLs detected in the present study were located relatively very close, maybe due to the high correlation among the traits. This information will serve as a significant basis for marker-assisted selection and molecular breeding in bitter gourd crop improvement.

A novel mutation in TFL1 homolog sustaining determinate growth in cucumber (Cucumis sativus L.)

BSA-seq combined with whole-genome resequencing map-based cloning delimited the cucumber det-novel locus into a 44.5 kb region in chromosome 6 harboring a putative candidate gene encoding a phosphatidylethanolamine-binding protein (CsCEN). Determinate and indeterminate growth habits of cucumber can affect plant architecture and crop yield. The TERMINAL FLOWER 1 (TFL1) controls determinate/indeterminate growth in Arabidopsis. In this study, a novel mutation in cucumber TFL1 homolog (CsCEN) has shown to regulate determinate growth and product of terminal flowers in cucumber (Cucumis sativus L.), which is similar to the function of CsTFL1 as previously reported. Genetic analysis in two determinate genotypes (D226 and D082) and indeterminate genotype (CCMC) revealed that a single recessive gene is responsible for this determinate growth trait. With the combination of BSA-seq and whole-genome resequencing, the locus of determinate-novel (det-novel) trait was mapped to a 44.5 kb genomic region in chromosome 6. Sequence alignment identified one non-synonymous SNP mutation (A to T) in the third exon of CsCEN, resulting in an amino acid substitution (Thr to Pro), suggesting that determinate growth might be controlled by a novel gene CsCEN (Csa6G152360) which differed from the reported CsTFL1 gene. The CsCEN expression level in shoot apexes and axillary buds was significantly lower in D226 compared to CCMC, suggesting its essential role in sustaining indeterminate growth habit. Identification and characterization of the CsCEN in the present study provide a new insight into plant architecture modification and development of cucumber cultivars suited to mechanized production system.

Insights on SNP types, detection methods and their utilization in Brassica species: Recent progress and future perspectives

The genus Brassica, family Brassicaceae (Cruciferae), comprises many important species of oil crops, vegetables and medicinal plants including B. rapa, B. oleracea, B. nigra, B. napus, B. juncea, B. carinata. Genomic researches in Brassica species is constrained by polyploidization, mainly due to its complicated genomic structure. However, rapid development of methods for detecting single nucleotide polymorphisms (SNP), such as next generation sequencing and SNP microarray, has accelerated release of reference Brassica species genomes as well as discovery of large numbers and genome-wide SNPs, thus intensifying forward genetics in this genus. In this review, we summarize biological characteristics, classification and various methods for detecting SNPs, focusing on high-throughput techniques. Moreover, we describe the pivotal roles of SNPs in genetic diversity, linkage map construction and QTL mapping, comparative genomics, linkage disequilibrium and genome-wide association studies. These insights are expected to deepen our understanding and guide further advancements in Brassica species research.

Localization of quantitative trait loci for cucumber fruit shape by a population of chromosome segment substitution lines

Cucumber fruit shape, a significant agronomic trait, is controlled by quantitative trait loci (QTLs). Feasibility of chromosome segment substitution lines (CSSLs) is well demonstrated to map QTLs, especially the minor-effect ones. To detect and identify QTLs with CSSLs can provide new insights into the underlying mechanisms regarding cucumber fruit shape. In the present study, 71 CSSLs were built from a population of backcross progeny (BC4F2) by using RNS7 (a round-fruit cucumber) as the recurrent parent and CNS21 (a long-stick-fruit cucumber) as the donor parent in order to globally detect QTLs for cucumber fruit shape. With the aid of 114 InDel markers covering the whole cucumber genome, 21 QTLs were detected for fruit shape-related traits including ovary length, ovary diameter, ovary shape index, immature fruit length, immature fruit diameter, immature fruit shape index, mature fruit length, mature fruit diameter and mature fruit shape index, and 4 QTLs for other traits including fruit ground and flesh color, and seed size were detected as well. Together our results provide important resources for the subsequent theoretical and applied researches on cucumber fruit shape and other traits.

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.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

Genome-wide identification of flowering time genes in cucurbit plants and revealed a gene ClGA2/KS associate with adaption and flowering of watermelon

Watermelon (Citrullus lanatus) is one of the major cucurbit crop that cultivated all over the world. Adaptability and flowering time are important agronomic characteristics that influence the quality and yield of watermelon, however, the molecular basis underlying these traits were still unclear. In this study, we identified 166, 182, 178, and 279 flowering genes in watermelon, melon, cucumber and pumpkin, respectively, and found that a lot of genes in the photoperiodic, autonomous, and vernalization pathways were absence in the four cucurbits. A higher ratio of flowering time genes was identified in the hormone pathway in cucurbits than in Arabidopsis, and a higher average ka/ks value of hormone pathway genes than the photoperiodic and vernalization pathway genes was identified in watermelon. Moreover, a gene ClGA2/KS (Cla005482) were found to associated with ecotype differentiation, flowering time, and whole growth period in watermelon. This study added knowledge to the molecular basis of flowering time regulation in cucurbits, and the molecule marker of ClGA2/KS gene may facilitate the breeding progress for selecting watermelon varieties with superior adaption and flowering time.

CsTFL1 inhibits determinate growth and terminal flower formation through interaction with CsNOT2a in cucumber

Cucumber (Cucumis sativus L.) is an important vegetable crop that carries on vegetative growth and reproductive growth simultaneously. Indeterminate growth is favourable for fresh market under protected environments, whereas determinate growth is preferred for pickling cucumber in the once-over mechanical harvest system. The genetic basis of determinacy is largely unknown in cucumber. In this study, map-based cloning of the de locus showed that the determinate growth habit is caused by a non-synonymous SNP in CsTFL1CsTFL1 is expressed in the subapical regions of the shoot apical meristem, lateral meristem and young stems. Ectopic expression of CsTFL1 rescued the terminal flower phenotype in the Arabidopsis tfl1-11 mutant and delayed flowering in wild-type Arabidopsis Knockdown of CsTFL1 resulted in determinate growth and formation of terminal flowers in cucumber. Biochemical analyses indicated that CsTFL1 interacts with a homolog of the miRNA biogenesis gene CsNOT2a; CsNOT2a interacts with FDP. Cucumber CsFT directly interacts with CsNOT2a and CsFD, and CsFD interacts with two 14-3-3 proteins. These data suggest that CsTFL1 competes with CsFT for interaction with CsNOT2a-CsFDP to inhibit determinate growth and terminal flower formation in cucumber.

Identification of quantitative trait loci governing subgynoecy in cucumber

QTL-seq analysis identified three major QTLs conferring subgynoecy in cucumbers. Furthermore, sequence and expression analyses predicted candidate genes controlling subgynoecy. The cucumber (Cucumis sativus L.) is a typical monoecious having individual male and female flowers, and sex differentiation is an important developmental process that directly affects its fruit yield. Subgynoecy represents a sex form with a high degree of femaleness and would have alternative use as gynoecy under limited resource conditions. Recently, many studies have been reported that QTL-seq, which integrates the advantages of bulked segregant analysis and high-throughput whole-genome resequencing, can be a rapid and cost-effective way of mapping QTLs. Segregation analysis in the F₂ and BC₁ populations derived from a cross between subgynoecious LOSUAS and monoecious BMB suggested the quantitative nature of subgynoecy in cucumbers. Both genome-wide SNP profiling of subgynoecious and monoecious bulks constructed from F₂ and BC₁ plants consistently identified three significant genomic regions, one on chromosome 3 (sg3.1) and another two on short and long arms of chromosome 1 (sg1.1 and sg1.2). Classical QTL analysis using the F₂ confirmed sg3.1 (R² = 42%), sg1.1 (R² = 29%) and sg1.2 (R² = 18%) as major QTLs. These results revealed the unique genetic inheritance of subgynoecious line LOSUAS through two distinct major QTLs, sg3.1 and sg1.1, which mainly increase degree of femaleness, while another QTL, sg1.2, contributes to decrease it. This study demonstrated that QTL-seq allows rapid and powerful detection of QTLs using preliminary generation mapping populations such as F₂ or BC₁ population and further that the identified QTLs could be useful for molecular breeding of cucumber lines with high yield potential.

The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development

Germplasm collections are a crucial resource to conserve natural genetic diversity and provide a source of novel traits essential for sustained crop improvement. Optimal collection, preservation and utilization of these materials depends upon knowledge of the genetic variation present within the collection. Here we use the high-throughput genotyping-by-sequencing (GBS) technology to characterize the United States National Plant Germplasm System (NPGS) collection of cucumber (Cucumis sativus L.). The GBS data, derived from 1234 cucumber accessions, provided more than 23 K high-quality single-nucleotide polymorphisms (SNPs) that are well distributed at high density in the genome (~1 SNP/10.6 kb). The SNP markers were used to characterize genetic diversity, population structure, phylogenetic relationships, linkage disequilibrium, and population differentiation of the NPGS cucumber collection. These results, providing detailed genetic analysis of the U.S. cucumber collection, complement NPGS descriptive information regarding geographic origin and phenotypic characterization. We also identified genome regions significantly associated with 13 horticulturally important traits through genome-wide association studies (GWAS). Finally, we developed a molecularly informed, publicly accessible core collection of 395 accessions that represents at least 96% of the genetic variation present in the NPGS. Collectively, the information obtained from the GBS data enabled deep insight into the diversity present and genetic relationships among accessions within the collection, and will provide a valuable resource for genetic analyses, gene discovery, crop improvement, and germplasm preservation.

Genetic dissection and fine mapping of a novel dt gene associated with determinate growth habit in sesame

BACKGROUND

As an important oil crop, growth habit of sesame (Sesamum indicum L.) is naturally indeterminate, which brings about asynchronous maturity of capsules and causes loss of yield.

RESULTS

The genetic basis of determinate growth habit in sesame was investigated by classical genetic analysis through multiple populations, results revealed that it was controlled by an unique recessive gene. The genotyping by sequencing (GBS) approach was employed for high-throughput SNP identification and genotyping in the F₂ population, then a high density bin map was constructed, the map was 1086.403 cM in length, which consisted of 1184 bins (13,679 SNPs), with an average of 0.918 cM between adjacent bins. Based on bin mapping in conjunction with SSR markers analysis in targeted region, the novel sesame determinacy gene was mapped on LG09 in a genome region of 41 kb.

CONCLUSIONS

This study dissected genetic basis of determinate growth habit in sesame, constructed a new high-density bin map and mapped a novel determinacy gene. Results of this study demonstrate that we employed an optimized approach to get fine-accuracy, high-resolution and high-efficiency mapping result in sesame. The findings provided important foundation for sesame determinacy gene cloning and were expected to be applied in breeding for cultivars suited to mechanized production.

LITTLELEAF (LL) encodes a WD40 repeat domain-containing protein associated with organ size variation in cucumber

Plants employ tight genetic control to integrate intrinsic growth signals and environmental cues to enable organs to grow to a defined size. Many genes contributing to cell proliferation and/or cell expansion, and consequently organ size control, have been identified, but the regulatory pathways are poorly understood. Here we have characterized a cucumber littleleaf (ll) mutant which exhibits smaller organ sizes but more lateral branches than the wild type. The small organ size in ll was due to a reduction of both cell number and cell size. Quantitative trait locus (QTL) analyses revealed co-localization of major-effect QTLs for fruit size, fruit and seed weight, as well as number of lateral branches, with the LL locus indicating pleiotropic effects of the ll mutation. We demonstrate that LL is an ortholog of Arabidopsis STERILE APETALA (SAP) encoding a WD40 repeat domain-containing protein; the mutant protein differed from the wild type by a single amino acid substitution (W264G) in the second WD40 repeat. W264 was conserved in 34 vascular plant genomes examined. Phylogenetic analysis suggested that LL originated before the emergence of flowering plants but was lost in the grass genome lineage. The function of LL in organ size control was confirmed by its overexpression in transgenic cucumbers and ectopic expression in Arabidopsis. Transcriptome profiling in LL and ll bulks revealed a complex regulatory network for LL-mediated organ size variation that involves several known organ size regulators and associated pathways. The data support LL as an important player in organ size control and lateral branch development in cucumber.

QTL-seq for rapid identification of candidate genes for flowering time in broccoli × cabbage

A major QTL controlling early flowering in broccoli × cabbage was identified by marker analysis and next-generation sequencing, corresponding to GRF6 gene conditioning flowering time in Arabidopsis. Flowering is an important agronomic trait for hybrid production in broccoli and cabbage, but the genetic mechanism underlying this process is unknown. In this study, segregation analysis with BC₁P1, BC₁P2, F₂, and F_2:3 populations derived from a cross between two inbred lines "195" (late-flowering) and "93219" (early flowering) suggested that flowering time is a quantitative trait. Next, employing a next-generation sequencing-based whole-genome QTL-seq strategy, we identified a major genomic region harboring a robust flowering time QTL using an F₂ mapping population, designated Ef2.1 on cabbage chromosome 2 for early flowering. Ef2.1 was further validated by indel (insertion or deletion) marker-based classical QTL mapping, explaining 51.5% (LOD = 37.67) and 54.0% (LOD = 40.5) of the phenotypic variation in F₂ and F_2:3 populations, respectively. Combined QTL-seq and classical QTL analysis narrowed down Ef1.1 to a 228-kb genomic region containing 29 genes. A cabbage gene, Bol024659, was identified in this region, which is a homolog of GRF6, a major gene regulating flowering in Arabidopsis, and was designated BolGRF6. qRT-PCR study of the expression level of BolGRF6 revealed significantly higher expression in the early flowering genotypes. Taken together, our results provide support for BolGRF6 as a possible candidate gene for early flowering in the broccoli line 93219. The identified candidate genomic regions and genes may be useful for molecular breeding to improve broccoli and cabbage flowering times.

Comparative analysis of miRNA and mRNA abundance in determinate cucumber by high-throughput sequencing

Determinate cucumber is a type of short vines, fewer nodes, and terminal flowers, it is suitable for high-density planting and available harvesting in field cultivation, whereas the indeterminate cucumber is preferred to cultivate in greenhouses. However, many biotic or abiotic stresses could lead indeterminate cucumber to be determinate in greenhouse cultivation, which may decrease yield and fruit quality. Therefore, it is urgent and essential to investigate the key factors forming determinate and terminal flowering in cucumber. In this study, two close background inbred lines were selected and conducted the miRNA and mRNA high throughput sequencing. Interestingly, ethylene-associated miRNAs and mRNAs were intensively obtained, indicating that the plant hormone ethylene is a key factor impacting determinate and terminal flowering in cucumber. The ethylene metabolites analysis showed that significant higher ethylene was observed in determinate line than that in the indeterminate line. The RT-qPCR validation of ethylene related miRNAs Cas-miR172, Cas-miR396, and Cas-miR414 and their target mRNAs showed a significant negative correlation. These data suggested that ethylene-associated miRNAs might affect determinate and terminal flower phenotypes by regulating their target genes expression. This study not only provides a potential molecular mechanism for determinate formation in cucumber but also establishes a method to demonstrate important physiological processes through the comprehensive association of miRNA and mRNA high-throughput sequencing.

Mapping and QTL Analysis of Gynoecy and Earliness in Bitter Gourd (Momordica charantia L.) Using Genotyping-by-Sequencing (GBS) Technology

A high-density, high-resolution genetic map was constructed for bitter gourd (Momordica charantia L.). A total of 2013 high quality SNP markers binned to 20 linkage groups (LG) spanning a cumulative distance of 2329.2 cM were developed. Each LG ranging from 185.2 cM (LG-12) to 46.2 cM (LG-17) and average LG span of 116.46 cM. The number of SNP markers mapped in each LG varied from 23 markers in LG-20 to 146 markers in LG-1 with an average of 100.65 SNPs per LG. The average distance between markers was 1.16 cM across 20 LGs and average distance between the markers ranged from 0.70 (LG-4) to 2.92 (LG-20). A total of 22 QTLs for four traits (gynoecy, sex ratio, node and days at first female flower appearance) were identified and mapped on 20 LGs. The gynoecious (gy-1) locus is flanked by markers TP_54865 and TP_54890 on LG 12 at a distance of 3.04 cM to TP_54890 and the major QTLs identified for the earliness traits will be extremely useful in marker development and MAS for rapid development of various gynoecious lines with different genetic background of best combiner for development of early and high yielding hybrids in bitter gourd.

A mutant in the CsDET2 gene leads to a systemic brassinosteriod deficiency and super compact phenotype in cucumber (Cucumis sativus L.)

A novel dwarf cucumber mutant, scp-2, displays a typical BR biosynthesis-deficient phenotype, which is due to a mutation in CsDET2 for a steroid 5-alpha-reductase. Brassinosteroids (BRs) are a group of plant hormones that play important roles in the development of plant architecture, and extreme dwarfism is a typical outcome of BR-deficiency. Most cucumber (Cucumis sativus L.) varieties have an indeterminate growth habit, and dwarfism may have its value in manipulation of plant architecture and improve production in certain production systems. In this study, we identified a spontaneous dwarf mutant, super compact-2 (scp-2), that also has dark green, wrinkle leaves. Genetic analyses indicated that scp-2 was different from two previously reported dwarf mutants: compact (cp) and super compact-1 (scp-1). Map-based cloning revealed that the mutant phenotype was due to two single nucleotide polymorphism and a single-base insertion in the CsDET2 gene that resulted in a missense mutation in a conserved amino acid and thus a truncated protein lacking the conserved catalytic domains in the predicted steroid 5α-reductase protein. Measurement of endogenous hormone levels indicated a reduced level of brassinolide (BL, a bioactive BR) in scp-2, and the mutant phenotype could be partially rescued by the application of epibrassinolide (EBR). In addition, scp-2 mutant seedlings exhibited dark-grown de-etiolation, and defects in cell elongation and vascular development. These data support that scp-2 is a BR biosynthesis-deficient mutant, and that the CsDET2 gene plays a key role in BR biosynthesis in cucumber. We also described the systemic BR responses and discussed the specific BR-related phenotypes in cucumber plants.

Transcriptomic analysis of short-fruit 1 (sf1) reveals new insights into the variation of fruit-related traits in Cucumis sativus

Fruit size is an important quality trait in different market classes of Cucumis sativus L., an economically important vegetable cultivated worldwide, but the genetic and molecular mechanisms that control fruit size are largely unknown. In this study, we isolated a natural cucumber mutant, short fruit 1 (sf1), caused by a single recessive Mendelian factor, from the North China-type inbred line CNS2. In addition to significantly decreased fruit length, other fruit-related phenotypic variations were also observed in sf1 compared to the wild-type (WT) phenotype, indicating that sf1 might have pleiotropic effects. Microscopic imaging showed that fruit cell size in sf1 was much larger than that in WT, suggesting that the short fruit phenotype in sf1 is caused by decreased cell number. Fine mapping revealed that sf1 was localized to a 174.3 kb region on chromosome 6. Similarly, SNP association analysis of bulked segregant RNA-Seq data showed increased SNP frequency in the same region of chromosome 6. In addition, transcriptomic analysis revealed that sf1 might control fruit length through the fine-tuning of cytokinin and auxin signalling, gibberellin biosynthesis and signal transduction in cucumber fruits. Overall, our results provide important information for further study of fruit length and other fruit-related features in cucumber.

QTL mapping of parthenocarpic fruit set in North American processing cucumber

Through a novel phenotyping method, four QTLs were consistently associated with increased parthenocarpic fruit set in North American processing cucumber that accounted for over 75 % of observed phenotypic variation. Parthenocarpy is a desirable trait with potential for increasing yield and quality in processing cucumber production. Although many successful parthenocarpic fresh market cucumber varieties have been developed, the genetic and molecular mechanisms behind parthenocarpic expression in cucumber remain largely unknown. Since parthenocarpy is an important yield component, it is difficult to separate the true parthenocarpic character from other yield related traits. In the present study, we developed a novel phenotypic approach for parthenocarpic fruit set focusing on early fruit development. Two hundred and five F₃ families derived from a cross between the highly parthenocarpic line 2A and low parthenocarpic line Gy8 were phenotypically evaluated in three greenhouse experiments. Seven QTLs associated with parthenocarpic fruit set were detected. Among them, one each on chromosomes 5 and 7 (parth5.1 and parth7.1) and two on chromosome 6 (parth6.1 and parth6.2) were consistently identified in all experiments, but their relative contribution to the total phenotypic variation was dependent on plant growth stages. While each of the four QTLs had almost equal contribution to the expression of the trait at commercial harvest stage, parth7.1 played an important role in early parthenocarpic fruit set. The results suggested that parthenocarpic fruit set can be accurately evaluated with as few as 20 nodes of growth. The QTLs identified in this study for parthenocarpic fruit set are a valuable resource for cucumber breeders interested in developing parthenocarpic cultivars and to researchers interested in the genetic and molecular mechanisms of parthenocarpic fruit set.

A High-Density Genetic Linkage Map for Cucumber (Cucumis sativus L.): Based on Specific Length Amplified Fragment (SLAF) Sequencing and QTL Analysis of Fruit Traits in Cucumber

High-density genetic linkage map plays an important role in genome assembly and quantitative trait loci (QTL) fine mapping. Since the coming of next-generation sequencing, makes the structure of high-density linkage maps much more convenient and practical, which simplifies SNP discovery and high-throughput genotyping. In this research, a high-density linkage map of cucumber was structured using specific length amplified fragment sequencing, using 153 F2 populations of S1000 × S1002. The high-density genetic map composed 3,057 SLAFs, including 4,475 SNP markers on seven chromosomes, and spanned 1061.19 cM. The average genetic distance is 0.35 cM. Based on this high-density genome map, QTL analysis was performed on two cucumber fruit traits, fruit length and fruit diameter. There are 15 QTLs for the two fruit traits were detected.

A major quantitative trait locus conferring subgynoecy in cucumber

A major QTL conditioning high degree of femaleness in cucumber was identified by marker analysis and next generation sequencing. Cucumber (Cucumis sativus L.) is a model species for sex determination studies, and its yield is associated with the degree of femaleness. Subgynoecy represents a sex form with a high degree of femaleness for which the genetic basis remains elusive. In this study, genetic analysis in the F2 and BC1 populations developed from a cross between subgynoecious S-2-98 and monoecious M95 suggested a quantitative nature of subgynoecy. Application of simple sequence repeat markers between subgynoecious and monoecious bulks constructed from BC1 plants identified three QTLs: sg3.1, sg6.1, and sg6.2. The major QTL sg3.1 contributed to 54.6% of the phenotypic variation, and its presence was confirmed by genome-wide comparison of SNP profiles between parental lines and a subgynoecious bulk constructed from BC6 plants. Using PCR-based markers developed from the SNP profile, sg3.1 was further delimited to a genomic region of 799 kb. The genetic basis of subgynoecy revealed here shall shed light on the development of elite cultivars with high yield potential.

QTL mapping of cucumber fruit flesh thickness by SLAF-seq

Cucumber is an agriculturally and economically important vegetable crop worldwide. Fruit flesh thickness is an important trait for cucumber and also a central determinant of yield, yet little is known about the underlying mechanism of this trait. In this study, bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) was applied to finely map the gene that underlies fruit flesh thickness in cucumber. A 0.19-Mb-long quantitative trait locus on chromosome 2 controlling fruit flesh thickness (QTL fft2.1) was identified and further confirmed by simple sequence repeat (SSR) marker-based classical QTL mapping in 138 F₂ individuals. Gene prediction of this 0.19-Mb region identified 20 genes. Quantitative RT-PCR revealed higher expression levels of Csa2 M058670.1 (SET domain protein-lysine methyltransferase) in D8 (thick fruit flesh parent) compared with that in XUE1 (thin fruit flesh parent) during fruit development. Sequence alignment analysis of Csa2M058670.1 from thick and thin fruit flesh cucumber lines revealed a 4-bp deletion mutation in the promoter region of this candidate gene, which may result in the loss of Csa2M058670.1 activation in thin fruit flesh lines. The data presented herein suggest that Csa2M058670.1 is a possible candidate gene for controlling flesh thickness in cucumber.

Ultrahigh-density linkage map for cultivated cucumber (Cucumis sativus L.) using a single-nucleotide polymorphism genotyping array

Genotyping arrays are tools for high-throughput genotyping, which is beneficial in constructing saturated genetic maps and therefore high-resolution mapping of complex traits. Since the report of the first cucumber genome draft, genetic maps have been constructed mainly based on simple-sequence repeats (SSRs) or on combinations of SSRs and sequence-related amplified polymorphism (SRAP). In this study, we developed the first cucumber genotyping array consisting of 32,864 single-nucleotide polymorphisms (SNPs). These markers cover the cucumber genome with a median interval of ~2 Kb and have expected genotype calls in parents/F₁ hybridizations as a training set. The training set was validated with Fluidigm technology and showed 96% concordance with the genotype calls in the parents/F₁ hybridizations. Application of the genotyping array was illustrated by constructing a 598.7 cM genetic map based on a ‘9930’ × ‘Gy14’ recombinant inbred line (RIL) population comprised of 11,156 SNPs. Marker collinearity between the genetic map and reference genomes of the two parents was estimated at R² = 0.97. We also used the array-derived genetic map to investigate chromosomal rearrangements, regional recombination rate, and specific regions with segregation distortions. Finally, 82% of the linkage-map bins were polymorphic in other cucumber variants, suggesting that the array can be applied for genotyping in other lines. The genotyping array presented here, together with the genotype calls of the parents/F₁ hybridizations as a training set, should be a powerful tool in future studies with high-throughput cucumber genotyping. An ultrahigh-density linkage map constructed by this genotyping array on RIL population may be invaluable for assembly improvement, and for mapping important cucumber QTLs.

Molecular mapping reveals structural rearrangements and quantitative trait loci underlying traits with local adaptation in semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis Qi et Yuan)

Comparative genetic mapping revealed the origin of Xishuangbanna cucumber through diversification selection after domestication. QTL mapping provided insights into the genetic basis of traits under diversification selection during crop evolution. The Xishuangbanna cucumber, Cucumis sativus L. var. xishuangbannanesis Qi et Yuan (XIS), is a semi-wild landrace from the tropical southwest China with some unique traits that are very useful for cucumber breeding, such as tolerance to low light, large fruit size, heavy fruit weight, and orange flesh color in mature fruits. In this study, using 124 recombinant inbred lines (RILs) derived from the cross of the XIS cucumber with a cultivated cucumber inbred line, we developed a linkage map with 269 microsatellite (or simple sequence repeat) markers which covered 705.9 cM in seven linkage groups. Comparative analysis of orders of common marker loci or marker-anchored draft genome scaffolds among the wild (C. sativus var. hardwickii), semi-wild, and cultivated cucumber genetic maps revealed that the XIS cucumber shares major chromosomal rearrangements in chromosomes 4, 5, and 7 between the wild and cultivated cucumbers suggesting that the XIS cucumber originated through diversifying selection after cucumber domestication. Several XIS-specific minor structural changes were identified in chromosomes 1 and 6. QTL mapping with the 124 RILs in four environments identified 13 QTLs for domestication and diversifying selection-related traits including 2 for first female flowering time (fft1.1, fft6.1), 5 for mature fruit length (fl1.1, fl3.1, fl4.1, fl6.1, and fl7.1), 3 for fruit diameter (fd1.1, fd4.1, and fd6.1), and 3 for fruit weight (fw2.1, fw4.1, and fw6.1). Six of the 12 QTLs were consistently detected in all four environments. Among the 13 QTLs, fft1.1, fl1.1, fl3.1, fl7.1, fd4.1, and fw6.1 were major-effect QTLs for respective traits with each explaining at least 10 % of the observed phenotypic variations. Results from this study provide insights into the cytological and genetic basis of crop evolution leading to the XIS cucumber. The molecular markers associated with the QTLs should be useful in exploring the XIS cucumber genetic resources for cucumber breeding.

An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing

Background

Cucumber, Cucumis sativus L., is an economically important vegetable crop which is processed or consumed fresh worldwide. However, the narrow genetic base in cucumber makes it difficult for constructing high-density genetic maps. The development of massively parallel genotyping methods and next-generation sequencing (NGS) technologies provides an excellent opportunity for developing single nucleotide polymorphisms (SNPs) for linkage map construction and QTL analysis of horticultural traits. Specific-length amplified fragment sequencing (SLAF-seq) is a recent marker development technology that allows large-scale SNP discovery and genotyping at a reasonable cost. In this study, we constructed a high-density SNP map for cucumber using SLAF-seq and detected fruit-related QTLs.

Results

An F₂ population of 148 individuals was developed from an intra-varietal cross between CC3 and NC76. Genomic DNAs extracted from two parents and 148 F₂ individuals were subjected to high-throughput sequencing and SLAF library construction. A total of 10.76 Gb raw data and 75,024,043 pair-end reads were generated to develop 52,684 high-quality SLAFs, out of which 5,044 were polymorphic. 4,817 SLAFs were encoded and grouped into different segregation patterns. A high-resolution genetic map containing 1,800 SNPs was constructed for cucumber spanning 890.79 cM. The average distance between adjacent markers was 0.50 cM. 183 scaffolds were anchored to the SNP-based genetic map covering 46% (168.9 Mb) of the cucumber genome (367 Mb). Nine QTLs for fruit length and weight were detected, a QTL designated fl3.2 explained 44.60% of the phenotypic variance. Alignment of the SNP markers to draft genome scaffolds revealed two mis-assembled scaffolds that were validated by fluorescence in situ hybridization (FISH).

Conclusions

We report herein the development of evenly dispersed SNPs across cucumber genome, and for the first time an SNP-based saturated linkage map. This 1,800-locus map would likely facilitate genetic mapping of complex QTL loci controlling fruit yield, and the orientation of draft genome scaffolds.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1158) contains supplementary material, which is available to authorized users.

Mapping and introgression of QTL involved in fruit shape transgressive segregation into ‘piel de sapo’ melon (cucumis melo l.) [corrected]

A mapping F₂ population from the cross ‘Piel de Sapo’ × PI124112 was selectively genotyped to study the genetic control of morphological fruit traits by QTL (Quantitative Trait Loci) analysis. Ten QTL were identified, five for FL (Fruit Length), two for FD (Fruit Diameter) and three for FS (Fruit Shape). At least one robust QTL per character was found, flqs8.1 (LOD = 16.85, R² = 34%), fdqs12.1 (LOD = 3.47, R² = 11%) and fsqs8.1 (LOD = 14.85, R² = 41%). flqs2.1 and fsqs2.1 cosegregate with gene a (andromonoecious), responsible for flower sex determination and with pleiotropic effects on FS. They display a positive additive effect (a) value, so the PI124112 allele causes an increase in FL and FS, producing more elongated fruits. Conversely, the negative a value for flqs8.1 and fsqs8.1 indicates a decrease in FL and FS, what results in rounder fruits, even if PI124112 produces very elongated melons. This is explained by a significant epistatic interaction between fsqs2.1 and fsqs8.1, where the effects of the alleles at locus a are attenuated by the additive PI124112 allele at fsqs8.1. Roundest fruits are produced by homozygous for PI124112 at fsqs8.1 that do not carry any dominant A allele at locus a (PiPiaa). A significant interaction between fsqs8.1 and fsqs12.1 was also detected, with the alleles at fsqs12.1 producing more elongated fruits. fsqs8.1 seems to be allelic to QTL discovered in other populations where the exotic alleles produce elongated fruits. This model has been validated in assays with backcross lines along 3 years and ultimately obtaining a fsqs8.1-NIL (Near Isogenic Line) in ‘Piel de Sapo’ background which yields round melons.

QTL-seq identifies an early flowering QTL located near Flowering Locus T in cucumber

Next-generation sequencing enabled a fast discovery of a major QTL controlling early flowering in cucumber, corresponding to the FT gene conditioning flowering time in Arabidopsis. Next-generation sequencing technologies are making it faster and more efficient to establish the association of agronomic traits with molecular markers or candidate genes, which is the requirement for marker-assisted selection in molecular breeding. Early flowering is an important agronomic trait in cucumber (Cucumis sativus L.), but the underlying genetic mechanism is unknown. In this study, we identified a candidate gene for early flowering QTL, Ef1.1 through QTL-seq. Segregation analysis in F2 and BC1 populations derived from a cross between two inbred lines "Muromskij" (early flowering) and "9930" (late flowering) suggested quantitative nature of flowering time in cucumber. Genome-wide comparison of SNP profiles between the early and late-flowering bulks constructed from F2 plants identified a major QTL, designated Ef1.1 on cucumber chromosome 1 for early flowering in Muromskij, which was confirmed by microsatellite marker-based classical QTL mapping in the F2 population. Joint QTL-seq and traditional QTL analysis delimited Ef1.1 to an 890 kb genomic region. A cucumber gene, Csa1G651710, was identified in this region, which is a homolog of the FLOWERING LOCUS T (FT), the main flowering switch gene in Arabidopsis. Quantitative RT-PCR study of the expression level of Csa1G651710 revealed significantly higher expression in early flowering genotypes. Data presented here provide support for Csa1G651710 as a possible candidate gene for early flowering in the cucumber line Muromskij.

Genetic linkage mapping and transmission ratio distortion in a three-generation four-founder population of Panicum virgatum (L.)

Switchgrass (Panicum virgatum L.), a warm season, C4, perennial grass, is one of the predominant grass species of the North American tall grass prairies. It is viewed as a high-potential bioenergy feedstock species because it can produce large amounts of lignocellulosic material with relatively few inputs. The objectives of this project were to develop an advanced switchgrass population and use it for the construction of genetic linkage maps and trait characterization. A three-generation, four-founder population was created and a total of 182 progeny of this advanced population were genotyped, including a mixture of self-pollinated and hybrid individuals. The female map integrated both subpopulations and covered 1629 cM of the switchgrass genome, with an average map length of 91 cM per linkage group. The male map of the hybrid progeny covered 1462 cM, with an average map length of 81 cM per linkage group. Average marker density of the female and male maps was 3.9 and 3.5 cM per marker interval, respectively. Based on the parental maps, the genome length of switchgrass was estimated to be 1776 cM and 1596 cM for the female map and male map, respectively. The proportion of the genome within 5 cM of a mapped locus was estimated to be 92% and 93% for the female map and male map, respectively. Thus, the linkage maps have covered most of the switchgrass genome. The assessment of marker transmission ratio distortion found that 26% of the genotyped markers were distorted from either 1:1 or 3:1 ratios expected for segregation of single dose markers in one or both parents, respectively. Several regions affected by transmission ratio distortion were found, with linkage groups Ib-m and VIIIa-f most affected.

A high-density genetic map of cucumber derived from Specific Length Amplified Fragment sequencing (SLAF-seq)

High-density genetic map provides an essential framework for accurate and efficient genome assembly and QTL fine mapping. Construction of high-density genetic maps appears more feasible since the advent of next-generation sequencing (NGS), which eases SNP discovery and high-throughput genotyping of large population. In this research, a high-density genetic map of cucumber (Cucumis sativus L.) was successfully constructed across an F2 population by a recently developed Specific Length Amplified Fragment sequencing (SLAF-seq) method. In total, 18.69 GB of data containing 93,460,000 paired-end reads were obtained after preprocessing. The average sequencing depth was 44.92 in the D8 (female parent), 42.16 in the Jin5-508 (male parent), and 5.01 in each progeny. 79,092 high-quality SLAFs were detected, of which 6784 SLAFs were polymorphic, and 1892 of the polymorphic markers met the requirements for constructing genetic map. The genetic map spanned 845.87 cm with an average genetic distance of 0.45 cm. It is a reliable linkage map for fine mapping and molecular breeding of cucumber for its high marker density and well-ordered markers.

Fine mapping a major QTL for kernel number per row under different phosphorus regimes in maize (Zea mays L.)

Phosphorus (P) is one of the essential macronutrients for plant growth and development. Grain yield is the primary trait of interest in maize breeding programs. Maize grain yield and yield-related traits are seriously affected by P deficiency. Kernel number per row (KN), as one of the major components of grain yield, has attracted the attention of more and more breeders. In our previous study, one major QTL (named qKN), controlling KN under different P regimes was mapped to the interval between molecular markers bnlg1360 and umc1645 on chromosome 10 using a F 2:3 population derived from the cross between maize inbreds 178 and 5,003 (107). In order to understand its genetic basis, we developed a population of near isogenic lines (NILs) and two P regimes were used to fine map and characterize qKN. The QTL qKN was finally localized in a region of ~480 kb. A single qKN allele of inbred 178 increased KN by 6.08-10.76 % in the 5,003 (107) background; qKN acted in a partially dominant manner. Our results will be instrumental for the future identification and isolation of the candidate gene underlying qKN. The tightly linked molecular markers that we developed for qKN will be useful in maize breeding programs for improving KN applying the marker-assisted selection.

Chromosomal Mapping and QTL Analysis of Resistance to Downy Mildew in Cucumis sativus

Downy mildew of cucumber (Cucumis sativus), caused by Pseudoperonospora cubensis, is a major foliar disease worldwide. The cucumber inbred lines K8 (resistant to downy mildew) and K18 (susceptible) were used to study the inheritance of resistance to downy mildew. Chromosomal mapping of the resistance genes was completed to provide a theoretical basis for the resistance mechanisms and for marker assisted selection (MAS). Inoculation was used to test the level of resistance to P. cubensis in the F₂ and F_2:3 families derived from the cross K8 × K18. Simple sequence repeat (SSR) analysis, combined with bulked segregation analysis (BSA), was done with the DNA of F₂ plants using 2,360 pairs of SSR primers. JoinMap Version 3.0 and MapInspect were used to construct SSR linkages and to verify the relationships between these SSR linkages and cucumber chromosomes. Quantitative trait locus (QTL) analysis of downy mildew resistance was done using MapQTL Version 4.0. Inheritance of resistance to downy mildew in K8 was quantitative. Five QTLs for resistance to downy mildew were detected: dm1.1, dm5.1, dm5.2, dm5.3, and dm6.1. The loci of dm1.1 and dm6.1 were on chromosomes 1 and 6, respectively. The loci of dm5.1, dm5.2, and dm5.3 were on chromosome 5, and were linked. Six linked SSR markers for these five QTLs were identified: SSR31116, SSR20705, SSR00772, SSR11012, SSR16882, and SSR16110. Six and four nucleotide binding site (NBS)-type resistance gene analogs (RGAs) were predicted in the region of dm5.2 and dm5.3, respectively. These results will be of benefit for fine-mapping the major QTLs for downy mildew resistance, and for MAS in cucumber.

Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly

Cucumber, Cucumis sativus L. is the only taxon with 2n = 2x = 14 chromosomes in the genus Cucumis. It consists of two cross-compatible botanical varieties: the cultivated C. sativus var. sativus and the wild C. sativus var. hardwickii. There is no consensus on the evolutionary relationship between the two taxa. Whole-genome sequencing of the cucumber genome provides a new opportunity to advance our understanding of chromosome evolution and the domestication history of cucumber. In this study, a high-density genetic map for cultivated cucumber was developed that contained 735 marker loci in seven linkage groups spanning 707.8 cM. Integration of genetic and physical maps resulted in a chromosome-level draft genome assembly comprising 193 Mbp, or 53% of the 367 Mbp cucumber genome. Strategically selected markers from the genetic map and draft genome assembly were employed to screen for fosmid clones for use as probes in comparative fluorescence in situ hybridization analysis of pachytene chromosomes to investigate genetic differentiation between wild and cultivated cucumbers. Significant differences in the amount and distribution of heterochromatins, as well as chromosomal rearrangements, were uncovered between the two taxa. In particular, six inversions, five paracentric and one pericentric, were revealed in chromosomes 4, 5 and 7. Comparison of the order of fosmid loci on chromosome 7 of cultivated and wild cucumbers, and the syntenic melon chromosome I suggested that the paracentric inversion in this chromosome occurred during domestication of cucumber. The results support the sub-species status of these two cucumber taxa, and suggest that C. sativus var. hardwickii is the progenitor of cultivated cucumber.

Construction of a high density integrated genetic map for cucumber (Cucumis sativus L.)

The high-density consensus map was constructed based on the GY14 × PI 183967 map from an inter-subspecific cross and the extended S94 × S06 map from an intra-subspecific cross. The consensus map was composed of 1,369 loci, including 1,152 SSR loci, 192 SRAP loci, 21 SCAR loci and one STS locus as well as three gene loci of fruit external quality traits in seven chromosomes, and spanned 700.5 cM, of which 682.7 cM (97.5%) were covered by SSR markers. The average genetic distance and physical interval between loci were 0.51 cM and ~268 kbp, respectively. Additionally, the physical position of the sequence-associated markers aligned along the assembled cucumber genome sequence established a relationship between genetic maps and cucumber genome sequence and to a great extent validated the order of markers in individual maps and consensus map. This consensus map with a high marker density and well-ordered markers is a saturated and reliable linkage map for genetic analysis of cucumber or the Cucurbitaceae family of plants.

High throughput screening with chlorophyll fluorescence imaging and its use in crop improvement

Marker assisted plant breeding is a powerful technique for targeted crop improvement in horticulture and agriculture. It depends upon the correlation of desirable phenotypic characteristics with specific genetic markers. This can be determined by statistical models that relate the variation in the value of genetic markers to variation in phenotypic traits. It therefore depends upon the convergence of three technologies; the creation of genetically characterised (and thus marked) populations, high throughput screening procedures, and statistical procedures. While a large number of high throughput screening technologies are available, real-time screening techniques are usually based on some kind of imaging technologies, such as chlorophyll fluorescence imaging, that offers physiological data that are eminently suitable as a quantitative trait for genetic mapping.

Syntenic relationships between cucumber (Cucumis sativus L.) and melon (C. melo L.) chromosomes as revealed by comparative genetic mapping

BACKGROUND

Cucumber, Cucumis sativus L. (2n = 2 × = 14) and melon, C. melo L. (2n = 2 × = 24) are two important vegetable species in the genus Cucumis (family Cucurbitaceae). Both species have an Asian origin that diverged approximately nine million years ago. Cucumber is believed to have evolved from melon through chromosome fusion, but the details of this process are largely unknown. In this study, comparative genetic mapping between cucumber and melon was conducted to examine syntenic relationships of their chromosomes.

RESULTS

Using two melon mapping populations, 154 and 127 cucumber SSR markers were added onto previously reported F(2)- and RIL-based genetic maps, respectively. A consensus melon linkage map was developed through map integration, which contained 401 co-dominant markers in 12 linkage groups including 199 markers derived from the cucumber genome. Syntenic relationships between melon and cucumber chromosomes were inferred based on associations between markers on the consensus melon map and cucumber draft genome scaffolds. It was determined that cucumber Chromosome 7 was syntenic to melon Chromosome I. Cucumber Chromosomes 2 and 6 each contained genomic regions that were syntenic with melon chromosomes III+V+XI and III+VIII+XI, respectively. Likewise, cucumber Chromosomes 1, 3, 4, and 5 each was syntenic with genomic regions of two melon chromosomes previously designated as II+XII, IV+VI, VII+VIII, and IX+X, respectively. However, the marker orders in several syntenic blocks on these consensus linkage maps were not co-linear suggesting that more complicated structural changes beyond simple chromosome fusion events have occurred during the evolution of cucumber.

CONCLUSIONS

Comparative mapping conducted herein supported the hypothesis that cucumber chromosomes may be the result of chromosome fusion from a 24-chromosome progenitor species. Except for a possible inversion, cucumber Chromosome 7 has largely remained intact in the past nine million years since its divergence from melon. Meanwhile, many structural changes may have occurred during the evolution of the remaining six cucumber chromosomes. Further characterization of the genomic nature of Cucumis species closely related to cucumber and melon might provide a better understanding of the evolutionary history leading to modern cucumber.

Primary trisomics obtained from autotriploid by diploid reciprocal crosses in cucumber

To promote cytogenetical studies on cucumber (Cucumis sativus L., 2n = 2x = 14), the reciprocal crosses were made between autotriploid and diploid for selecting the primary trisomics. Meanwhile, chromosome behavior during meiosis in autotriploid cucumber was investigated to look for cytological evidences for origin of primary trisomics. Many viable F(1) seeds were obtained from reciprocal crosses between autotriploid and diploid. The number of chromosomes of 56 surviving progenies varied from 14 to 28, with plants having 2n = 15 occurring at the highest frequency (51.8%). Primary trisomics were firstly obtained in this study. Four types of primary trisomics were isolated and they could be distinguished from each other, as well as diploid. Variable chromosome configurations, e.g. univalent, bivalents and trivalents were observed in many pollen mother cells of the autotriploid at metaphase I. Binomial chromosome distribution was observed at anaphase I and frequency of 8/13 was 6.25%. The meiosis of autotriploid, especially the class of gametes with eight chromosomes, gave the cytological evidence of producing 2x + 1 type gamete and could be induced into primary trisomic plants from progeny of autotriploid-diploid crosses. These studies have established a ground work for selecting a series of primary trisomics, and further using them for associating linkage groups with specific chromosomes in cucumber.

Comparative analysis of marker-assisted and phenotypic selection for yield components in cucumber

Theoretical studies suggest that marker-assisted selection (MAS) has case-specific advantages over phenotypic selection (PHE) for selection of quantitative traits. However, few studies have been conducted that empirically compare these selection methods in the context of a plant breeding program. For direct comparison of the effectiveness of MAS and PHE, four cucumber (Cucumis sativus L.; 2n = 2x = 14) inbred lines were intermated and then maternal bulks were used to create four base populations for recurrent mass selection. Each of these populations then underwent three cycles of PHE (open-field evaluations), MAS (genotyping at 18 marker loci), and random mating without selection. Both MAS and PHE were practiced for yield indirectly by selecting for four yield-component traits that are quantitatively inherited with 2-6 quantitative trait loci per trait. These traits were multiple lateral branching, gynoecious sex expression (gynoecy), earliness, and fruit length to diameter ratio. Both MAS and PHE were useful for multi-trait improvement, but their effectiveness depended upon the traits and populations under selection. Both MAS and PHE provided improvements in all traits under selection in at least one population, except for earliness, which did not respond to MAS. The populations with maternal parents that were inferior for a trait responded favorably to both MAS and PHE, while those with maternal parents of superior trait values either did not change or decreased during selection. Generally, PHE was most effective for gynoecy, earliness, and fruit length to diameter ratio, while MAS was most effective for multiple lateral branching and provided the only increase in yield (fruit per plant).

An integrated genetic and cytogenetic map of the cucumber genome

The Cucurbitaceae includes important crops such as cucumber, melon, watermelon, squash and pumpkin. However, few genetic and genomic resources are available for plant improvement. Some cucurbit species such as cucumber have a narrow genetic base, which impedes construction of saturated molecular linkage maps. We report herein the development of highly polymorphic simple sequence repeat (SSR) markers originated from whole genome shotgun sequencing and the subsequent construction of a high-density genetic linkage map. This map includes 995 SSRs in seven linkage groups which spans in total 573 cM, and defines approximately 680 recombination breakpoints with an average of 0.58 cM between two markers. These linkage groups were then assigned to seven corresponding chromosomes using fluorescent in situ hybridization (FISH). FISH assays also revealed a chromosomal inversion between Cucumis subspecies [C. sativus var. sativus L. and var. hardwickii (R.) Alef], which resulted in marker clustering on the genetic map. A quarter of the mapped markers showed relatively high polymorphism levels among 11 inbred lines of cucumber. Among the 995 markers, 49%, 26% and 22% were conserved in melon, watermelon and pumpkin, respectively. This map will facilitate whole genome sequencing, positional cloning, and molecular breeding in cucumber, and enable the integration of knowledge of gene and trait in cucurbits.

Detection of QTL for yield-related traits using recombinant inbred lines derived from exotic and elite US Western Shipping melon germplasm

The inheritance of yield-related traits in melon (Cucumis melo L.; 2n = 2x = 24) is poorly understood, and the mapping of quantitative trait loci (QTL) for such traits has not been reported. Therefore, a set of 81 recombinant inbred lines (RIL) was developed from a cross between the monoecious, highly branched line USDA 846-1 and a standard vining, andromonoecious cultivar, 'Top Mark'. The RIL, parental lines, and three control cultivars ('Esteem', 'Sol Dorado', and 'Hales Best Jumbo') were grown at Hancock, WI and El Centro, CA in 2002, and evaluated for primary branch number (PB), fruit number per plant (FN), fruit weight per plant (FW), average weight per fruit (AWF), and percentage of mature fruit per plot (PMF). A 190-point genetic map was constructed using 114 RAPD, 43 SSR, 32 AFLP markers, and one phenotypic trait. Fifteen linkage groups spanned 1,116 cM with a mean marker interval of 5.9 cM. A total of 37 QTL were detected in both locations (PB = 6, FN = 9, FW = 12, AWF = 5, and PMF = 5). QTL analyses revealed four location-independent factors for PB (pb1.1, pb1.2, pb2.3, and pb10.5), five for FN (fn1.1, fn1.2, fn1.3, fn2.4, and fn8.8), four for FW (fw5.8, fw6.10, fw8.11, and fw8.12), two for AWF (awf1.3 and awf8.5), and one for PMF (pmf10.4). The significant (P </= 0.05) positive phenotypic correlations observed among PB, FN, and FW, and negative phenotypic correlations between PB and AWF and between FN and AWF were consistent with the genomic locations and effects (negative vs. positive) of the QTL detected. Results indicate that genes resident in highly branched melon types have potential for increasing yield in US Western Shipping type germplasm via marker-assisted selection.

Population development by phenotypic selection with subsequent marker-assisted selection for line extraction in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.; 2n = 2x = 14) has a narrow genetic base, and commercial yield of US processing cucumber has plateaued in the last 15 years. Yield may be increased by altering plant architecture to produce unique early flowering (days to flower, DTF), female (gynoecious, GYN), highly branched (multiple lateral branching, MLB), long-fruited (length:diameter ratio, L:D) cultivars with diverse plant statures. The genetic map position of QTL conditioning these quantitatively inherited yield component traits is known, and linked molecular markers may have utility in marker-assisted selection (MAS) programs to increase selection efficiency, and effectiveness. Therefore, a base population (C0), created by intermating four unique but complementary lines, was subjected to three cycles (C1-C3) of phenotypic (PHE) mass selection for DTF, GYN, MLB, and L:D. In tandem, two cycles of marker-assisted backcrossing for these traits began with selected C2 progeny (C2S) to produce families (F1[i.e., C2S x C2S], and BC(1) [i.e., F1 x C2S]) for line extraction, and for comparative analysis of gain from selection by PHE selection, and MAS. Frequencies of marker loci were used to monitor selection-dependent changes during PHE selection, and MAS. Similar gain from selection was detected as a result of PHE selection, and MAS for MLB (approximately 0.3 branches/cycle), and L:D (approximately 0.1 unit increase/cycle) with concomitant changes in frequency at linked marker loci. Although genetic gain was not realized for GYN during PHE selection, the percentage of female flowers of plants subjected to MAS was increased (5.6-9.8% per cycle) depending upon the BC1 population examined. Selection-dependent changes in frequency were also detected at marker loci linked to female sex expression during MAS. MAS operated to fix favorable alleles that were not exploited by PHE selection in this population, indicating that MAS could be applied for altering plant architecture in cucumber to improve its yield potential.

QTL analysis of powdery mildew resistance in cucumber (Cucumis sativus L.)

A population of F7 recombinant inbred lines (RILs) was made from a cross between susceptible ('Santou') and resistant (PI197088-1) lines of cucumber in order to study powdery mildew resistance loci. Susceptibility to powdery mildew in the F7 RIL individuals showed a continuous distribution from susceptible to resistant, suggesting that powdery mildew resistance is controlled by quantitative trait loci (QTLs). A QTL analysis identified two and three loci for powdery mildew resistance under 26 and 20 degrees C conditions, respectively. One QTL was found in the same position under both temperature conditions. Therefore, it is more likely that one major QTL acts under both temperature conditions and that other QTLs are specific to the two temperature conditions. The above results suggest that the four QTLs are controlled in a different temperature manner, and that their combination played an important role in expressing a high level of resistance to powdery mildew in this cucumber population. Sequence-tagged site (STS) markers associated with each QTL were developed and would be useful for breeding a cucumber line with a high level of powdery mildew resistance.

Construction of a cucumber genetic linkage map with SRAP markers and location of the genes for lateral branch traits

Using SRAP (sequence-related amplified polymorphism) markers a genetic linkage map of cucumber was constructed with a population consisting of 138 F2 individuals derived from a cross of the two cucumber lines, S06 and S52. In the survey of parental polymorphisms with 182 primer combinations, 64 polymorphism-revealing primer pairs were screened out, which generated totally 108 polymorphic bands with an average of 1.7 bands per primer pair and at most 6 bands from one primer pair. The constructed molecular linkage map included 92 loci, distributed in seven linkage groups and spanning 1164.2 cM in length with an average genetic distance of 12.6 cM between two neighboring loci. Based on this linkage map, the quantitative trait loci (QTL) for the lateral branch number (Ibn) and the lateral branch average length (Ibl) in cucumber were identified by QTLMapperl.6. A major QTL Ibn1 located between ME11 SA4B and ME5EM5 in LG2 could explain 10.63% of the total variation with its positively effecting allele from S06. A major QTL Ibl1 located between DC1OD3 and DC1EM14 in LG2 could account for 10.38% of the total variation with its positively effecting allele from S06.

Construction of two BAC libraries from cucumber (Cucumis sativus L.) and identification of clones linked to yield component quantitative trait loci

Two bacterial artificial chromosome (BAC) libraries were constructed from an inbred line derived from a cultivar of cucumber (Cucumis sativus L.). Intact nuclei were isolated and embedded in agarose plugs, and high-molecular-weight DNA was subsequently partially digested with BamHI or EcoRI. Ligation of double size-selected DNA fragments with the pECBAC1 vector yielded two libraries containing 23,040 BamHI and 18,432 EcoRI clones. The average BamHI and EcoRI insert sizes were estimated to be 107.0 kb and 100.8 kb, respectively, and BAC clones lacking inserts were 1.3% and 14.5% in the BamHI and EcoRI libraries, respectively. The two libraries together represent approximately 10.8 haploid cucumber genomes. Hybridization with a C(0)t-1 DNA probe revealed that approximately 36% of BAC clones likely carried repetitive sequence-enriched DNA. The frequencies of BAC clones that carry chloroplast or mitochondrial DNA range from 0.20% to 0.47%. Four sequence-characterized amplified region (SCAR), four simple sequence repeat, and an randomly amplified polymorphic DNA marker linked with yield component quantitative trait loci were used either as probes to hybridize high-density colony filters prepared from both libraries or as primers to screen an ordered array of pooled BAC DNA prepared from the BamHI library. Positive BAC clones were identified in predicted numbers, as screening by polymerase chain reaction amplification effectively overcame the problems associated with an overabundance of positives from hybridization with two SCAR markers. The BAC clones identified herein that are linked to the de (determinate habit) and F (gynoecy) locus will be useful for positional cloning of these economically important genes. These BAC libraries will also facilitate physical mapping of the cucumber genome and comparative genome analyses with other plant species.