High-throughput SNP genotyping in Cucurbita pepo for map construction and quantitative trait loci mapping

Genetic mapping reveals a candidate gene CmoFL1 controlling fruit length in pumpkin

Fruit length (FL) is an important economical trait that affects fruit yield and appearance. Pumpkin (Cucurbita moschata Duch) contains a wealth genetic variation in fruit length. However, the natural variation underlying differences in pumpkin fruit length remains unclear. In this study, we constructed a F2 segregate population using KG1 producing long fruit and MBF producing short fruit as parents to identify the candidate gene for fruit length. By bulked segregant analysis (BSA-seq) and Kompetitive Allele-Specific PCR (KASP) approach of fine mapping, we obtained a 50.77 kb candidate region on chromosome 14 associated with the fruit length. Then, based on sequence variation, gene expression and promoter activity analyses, we identified a candidate gene (CmoFL1) encoding E3 ubiquitin ligase in this region may account for the variation of fruit length. One SNP variation in promoter of CmoFL1 changed the GT1CONSENSUS, and DUAL-LUC assay revealed that this variation significantly affected the promoter activity of CmoFL1. RNA-seq analysis indicated that CmoFL1 might associated with the cell division process and negatively regulate fruit length. Collectively, our work identifies an important allelic affecting fruit length, and provides a target gene manipulating fruit length in future pumpkin breeding.

Stress induced production of plant secondary metabolites in vegetables: Functional approach for designing next generation super foods

Plant secondary metabolites are vital for human health leading to the gain the access to natural products. The quality of crops is the result of the interaction of different biotic and abiotic factors. Abiotic stresses during plant growth may reduce the crop performance and quality of the produce. However, abiotic stresses can result in numerous physiological, biochemical, and molecular responses in plants, aiming to deal with these conditions. Abiotic stresses are also elicitors of the biosynthesis of plant secondary metabolites in plants which possess plant defense mechanisms as well as human health benefits such as anti-inflammatory, antioxidative properties etc. Plants either synthesize new compounds or alter the concentration of bioactive compounds. Due to increasing attention towards the production of bioactive compounds, the understanding of crop responses to abiotic stresses in relation to the biosynthesis of bioactive compounds is critical. Plants alter their metabolism at the genetic level in response to different abiotic stresses resulting the changes in secondary metabolite production. Transcriptional factors regulate genes responsible for secondary metabolite biosynthesis in several plants under stress conditions. Understanding the signaling pathways involved in the secondary metabolite biosynthesis has become easy with the use of molecular biology. Therefore, aim of writing the review is to focus on secondary metabolite production in vegetable crops, their health benefits and transcription regulation under various abiotic stresses.

Carotenoid fortification of zucchini fruits using a viral RNA vector

BACKGROUND

Carotenoids are health-promoting metabolites in livestock and human diets. Some important crops have been genetically modified to increase their content. Although the usefulness of transgenic plants to alleviate nutritional deficiencies is obvious, their social acceptance has been controversial.

RESULTS

Here, we demonstrate an alternative biotechnological strategy for carotenoid fortification of edible fruits in which no transgenic DNA is involved. A viral RNA vector derived from Zucchini yellow mosaic virus (ZYMV) was modified to express a bacterial phytoene synthase (crtB), and inoculated to zucchini (Cucurbita pepo L.) leaves nurturing pollinated flowers. After the viral vector moved to the developing fruit and expressed crtB, the rind and flesh of the fruits developed yellow-orange rather than green color. Metabolite analyses showed a substantial enrichment in health-promoting carotenoids, such as α- and β-carotene (provitamin A), lutein and phytoene, in both rind and flesh.

CONCLUSION

Although this strategy is perhaps not free from controversy due to the use of genetically modified viral RNA, our work does demonstrate the possibility of metabolically fortifying edible fruits using an approach in which no transgenes are involved. This article is protected by copyright. All rights reserved.

QTL mapping identifies novel loci and putative candidate genes for rind traits in watermelon

Rind color and pattern are important external fruit quality traits of watermelon influencing the consumer attention and market acceptance. Varying degrees of green and yellow colors in different combinations forming distinct patterns have been observed on watermelon rind. In the current experiment, an attempt has been made to map the QTLs/genes for four important rind traits, namely, stripe color, stripe pattern, interstripe color, and interstripe pattern. The experiment consisted of two mapping populations namely F₂ (Pop I) and BC₁F₂ (Pop II) derived from two parental lines, viz., BIL-53 (characterized by medium green marbled rind pattern and yellowish white blotchy interstripe pattern) and IIHR-140-152 (characterized by dark green solid stripes and greenish wavy interstripe pattern). Linkage mapping identified consistent QTLs across populations on chromosome 9. Comparative genomic analysis of these regions revealed two genes, namely, Cla97C09G175170 and Cla97C09G175150 as potential candidates for stripe and interstripe color. Sequence analysis of Cla97C09G175170 gene in parents along with reference genotypes, viz., 97,103 and Charleston gray suggests a 3 bp deletion on 11th exon to be associated with stripe color polymorphism in watermelon.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03112-7.

A Novel QTL for Resistance to Phytophthora Crown Rot in Squash

Phytophthora capsici Leonian causes significant yield losses in commercial squash (Cucurbita pepo) production worldwide. The deployment of resistant cultivars can complement integrated management practices for P. capsici, but resistant cultivars are currently unavailable for growers. Moderate resistance to Phytophthora crown rot in a selection of accession PI 181761 (C. pepo) (designated line #181761-36P) is controlled by three dominant genes (R4, R5 and R6). Introgression of these loci into elite germplasm through marker-assisted selection (MAS) can accelerate the release of new C. pepo cultivars resistant to crown rot, but these tools are currently unavailable. Here we describe the identification of a quantitative trait locus (QTL), molecular markers and candidate genes associated with crown rot resistance in #181761-36P. Five hundred and twenty-three SNP markers were genotyped in an F₂ (n = 83) population derived from a cross between #181761-36P (R) and Table Queen (S) using targeted genotyping by sequencing. A linkage map (2068.96 cM) consisting of twenty-one linkage groups and an average density of 8.1 markers/cM was developed for the F₂ population. The F₂:₃ families were phenotyped in the greenhouse with a virulent strain of P. capsica, using the spore-spray method. A single QTL (QtlPC-C13) was consistently detected on LG 13 (chromosome 13) across three experiments and explained 17.92-21.47% of phenotypic variation observed in the population. Nine candidate disease resistance gene homologs were found within the confidence interval of QtlPC-C13. Single nucleotide polymorphism (SNP) markers within these genes were converted into Kompetitive Allele Specific PCR (KASP) assays and tested for association with resistance in the F₂ population. One SNP marker (C002686) was significantly associated with resistance to crown rot in the F₂ population (p < 0.05). This marker is a potential target for MAS for crown rot resistance in C. pepo.

Genome Wide Characterization, Comparative and Genetic Diversity Analysis of Simple Sequence Repeats in Cucurbita Species

Simple sequence repeats (SSRs) are widely used in mapping constructions and comparative and genetic diversity analyses. Here, 103,056 SSR loci were found in Cucurbita species by in silico PCR. In general, the frequency of these SSRs decreased with the increase in the motif length, and di-nucleotide motifs were the most common type. For the same repeat types, the SSR frequency decreased sharply with the increase in the repeat number. The majority of the SSR loci were suitable for marker development (84.75% in Cucurbita moschata, 94.53% in Cucurbita maxima, and 95.09% in Cucurbita pepo). Using these markers, the cross-species transferable SSR markers between C. pepo and other Cucurbitaceae species were developed, and the complicated mosaic relationships among them were analyzed. Especially, the main syntenic relationships between C. pepo and C. moschata or C. maxima indicated that the chromosomes in the Cucurbita genomes were highly conserved during evolution. Furthermore, 66 core SSR markers were selected to measure the genetic diversity in 61 C. pepo germplasms, and they were divided into two groups by structure and unweighted pair group method with arithmetic analysis. These results will promote the utilization of SSRs in basic and applied research of Cucurbita species.

A novel introgression line collection to unravel the genetics of climacteric ripening and fruit quality in melon

Introgression lines are valuable germplasm for scientists and breeders, since they ease genetic studies such as QTL interactions and positional cloning as well as the introduction of favorable alleles into elite varieties. We developed a novel introgression line collection in melon using two commercial European varieties with different ripening behavior, the climacteric cantalupensis 'Védrantais' as recurrent parent and the non-climacteric inodorus 'Piel de Sapo' as donor parent. The collection contains 34 introgression lines, covering 99% of the donor genome. The mean introgression size is 18.16 Mb and ~ 3 lines were obtained per chromosome, on average. The high segregation of these lines for multiple fruit quality traits allowed us to identify 27 QTLs that modified sugar content, altered fruit morphology or were involved in climacteric ripening. In addition, we confirmed the genomic location of five major genes previously described, which control mainly fruit appearance, such as mottled rind and external color. Most of the QTLs had been reported before in other populations sharing parental lines, while three QTLs (EAROQP11.3, ECDQP11.2 and FIRQP4.1) were newly detected in our work. These introgression lines would be useful to perform additional genetic studies, as fine mapping and gene pyramiding, especially for important complex traits such as fruit weight and climacteric ripening.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

A combined BSA-Seq and linkage mapping approach identifies genomic regions associated with Phytophthora root and crown rot resistance in squash

Two QTL mapping approaches were used to identify a total of six QTL associated with Phytophthora root and crown rot resistance in a biparental squash population. Phytophthora root and crown rot, caused by the soilborne oomycete pathogen Phytophthora capsici, leads to severe yield losses in squash (Cucurbita pepo). To identify quantitative trait loci (QTL) involved in resistance to this disease, we crossed a partially resistant squash breeding line with a susceptible zucchini cultivar and evaluated over 13,000 F₂ seedlings in a greenhouse screen. Bulked segregant analysis with whole genome resequencing (BSA-Seq) resulted in the identification of five genomic regions-on chromosomes 4, 5, 8, 12, and 16-featuring significant allele frequency differentiation between susceptible and resistant bulks in each of two independent replicates. In addition, we conducted linkage mapping using a population of 176 F₃ families derived from individually genotyped F₂ individuals. Variation in disease severity among these families was best explained by a four-QTL model, comprising the same loci identified via BSA-Seq on chromosomes 4, 5, and 8 as well as an additional locus on chromosome 19, for a combined total of six QTL identified between both methods. Loci, whether those identified by BSA-Seq or linkage mapping, were of small-to-moderate effect, collectively accounting for 28-35% and individually for 2-10% of the phenotypic variance explained. However, a multiple linear regression model using one marker in each BSA-Seq QTL could predict F_2:3 disease severity with only a slight drop in cross-validation accuracy compared to genomic prediction models using genome-wide markers. These results suggest that marker-assisted selection could be a suitable approach for improving Phytophthora crown and root rot resistance in squash.

QTL Mapping for Gummy Stem Blight Resistance in Watermelon (Citrullus spp.)

Watermelon (Citrulluslanatus) is an economically important fruit crop worldwide. Gummy stem blight (GSB) is one of the most damaging diseases encountered during watermelon cultivation. In the present study, we identified quantitative trait loci (QTLs) associated with GSB resistance in an F2 population derived from a cross between maternal-susceptible line '920533' (C. lanatus) and the paternal-resistant line 'PI 189225' (C. amarus). The resistance of 178 F2 plants was assessed by two different evaluation methods, including leaf lesion (LL) and stem blight (SB). To analyze the QTLs associated with GSB resistance, a linkage map was constructed covering a total genetic distance of 1070.2 cM. QTL analysis detected three QTLs associated with GSB resistance on chromosome 8 and 6. Among them, two QTLs, qLL8.1 and qSB8.1 on chromosome 8 identified as major QTLs, explaining 10.5 and 10.0% of the phenotypic variations localizing at same area and sharing the same top markers for both LL and SB traits, respectively. A minor QTL, qSB6.1, explains 9.7% of phenotypic variations detected on chromosome 6 only for the SB trait. High-throughput markers were developed and validated for the selection of resistant QTLs using watermelon accessions, and commercial cultivars. Four potential candidate genes were predicted associated with GSB resistance based on the physical location of flanking markers on chromosome 8. These findings will be helpful for the development of watermelon cultivars resistant to GSB.

Quantitative Trait Locus Analysis in Squash (Cucurbita moschata) Based on Simple Sequence Repeat Markers and Restriction Site-Associated DNA Sequencing Analysis

Squash (Cucurbita moschata) displays wide morphological and genetic variations; however, limited information is available regarding the genetic loci of squash that control its agronomic traits. To obtain basic genetic information for C. moschata, an F2 population was prepared derived from a cross between the Vietnamese cultivar ‘Bí Hồ Lô TN 6 (TN 6)’ and the Japanese cultivar ‘Shishigatani’, and flowering and fruit traits were examined. Overall, the traits showed a continuous distribution in the F2 population, suggesting that they were quantitative traits. A linkage map was constructed based on simple sequence repeat and restriction site-associated DNA (RAD) markers to detect quantitative trait loci (QTLs). Twelve QTLs for flowering and fruit traits, as well as one phenotypic trait locus, were successfully localized on the map. The present QTLs explained the phenotypic variations at a moderate to relatively high level (16.0%–47.3%). RAD markers linked to the QTLs were converted to codominant cleaved amplified polymorphic sequence (CAPS) and derived CAPS markers for the easy detection of alleles. The information reported here provides useful information for understanding the genetics of Cucurbita and other cucurbit species, and for the selection of individuals with ideal traits during the breeding of Cucurbita vegetables.

High-density SNP map facilitates fine mapping of QTLs and candidate genes discovery for Aspergillus flavus resistance in peanut (Arachis hypogaea)

Two novel resistant QTLs mapped and candidate genes identified for Aspergillus flavus resistance in cultivated peanut using SLAF-seq. Aflatoxin contamination in peanuts caused by Aspergillus flavus is a serious food safety issue for human health around the world. Host plant resistance to fungal infection and reduction in aflatoxin are crucial for mitigating this problem. Identification of the resistance-linked markers can be used in marker-assisted breeding for varietal development. Here we report construction of two high-density genetic linkage maps with 1975 SNP loci and 5022 SNP loci, respectively. Two consistent quantitative trait loci (QTL) were identified as qRAF-3-1 and qRAF-14-1, which located on chromosomes A03 and B04, respectively. QTL qRAF-3-1 was mapped within 1.67 cM and had more than 19% phenotypic variance explained (PVE), while qRAF-14-1 was located within 1.34 cM with 5.15% PVE. While comparing with the reference genome, the mapped QTLs, qRAF-3-1 and qRAF-14-1, were located within a physical distance of 1.44 Megabase pair (Mbp) and 2.22 Mbp, harboring 67 and 137 genes, respectively. Among the identified candidate genes, six genes with the same function were found within both QTLs regions. In addition, putative disease resistance RPP13-like protein 1 (RPP13), lipoxygenase (Lox), WRKY transcription factor (WRKY) and cytochrome P450 71B34 genes were also identified. Using microarray analysis, genes responded to A. flavus infection included coding for RPP13, pentatricopeptide repeat-containing-like protein, and Lox which may be possible candidate genes for resistance to A. flavus. The QTLs and candidate genes will further facilitate marker development and validation of genes for deployment in the molecular breeding programs against A. flavus in peanuts.

Identification of Fruit-Associated QTLs in Winter Squash (Cucurbita maxima Duchesne) Using Recombinant Inbred Lines

Cucurbita maxima Duchesne squash and pumpkins are cultivated world-wide. Cucurbita maxima fruits are produced for fresh market and are valuable for food processing. Therefore, fruit characteristics and yield are the traits of high economic importance for breeders. To date, the genetic basis of fruit-associated traits in C. maxima have been poorly understood. In the present study, we evaluated fruit-associated traits and conducted quantitative trait locus (QTL) analysis using recombinant inbred lines (RILs) derived from a cross of two inbred lines with different fruit morphotypes. Phenotypic data for nine fruit traits (earliness, weight, number per plant, yield per plant, length and diameter, shape index, flesh thickness, sucrose content and dry matter content) were collected for RILs in two open-field experiments. Pairwise analysis of the phenotypic data revealed correlations among the fruit and yield-associated traits. Using a previously developed genetic map, we identified 26 QTLs for eight traits. The QTLs were found in 10 locations on eight chromosomes of C. maxima. The QTLs were detected across experiments and explained up to 41.4% of the observed phenotypic variations. Major-effect QTLs for multiple fruit-associated traits were clustered on chromosome 4, suggesting that this genomic region has been under selection during diversification and/or domestication of C. maxima.

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.

Construction, characteristics and high throughput molecular screening methodologies in some special breeding populations: a horticultural perspective

Advanced marker technologies are widely used for evaluation of genetic diversity in cultivated crops, wild ancestors, landraces or any special plant genotypes. Developing agricultural cultivars requires the following steps: (i) determining desired characteristics to be improved, (ii) screening genetic resources to help find a superior cultivar, (iii) intercrossing selected individuals, (iv) generating genetically hybrid populations and screening them for agro-morphological or molecular traits, (v) evaluating the superior cultivar candidates, (vi) testing field performance at different locations, and (vii) certifying. In the cultivar development process valuable genes can be identified by creating special biparental or multiparental populations and analysing their association using suitable markers in given populations. These special populations and advanced marker technologies give us a deeper knowledge about the inherited agronomic characteristics. Unaffected by the changing environmental conditions, these provide a higher understanding of genome dynamics in plants. The last decade witnessed new applications for advanced molecular techniques in the area of breeding,with low costs per sample. These, especially, include next-generation sequencing technologies like reduced representation genome sequencing (genotyping by sequencing, restriction site-associated DNA). These enabled researchers to develop new markers, such as simple sequence repeat and single- nucleotide polymorphism, for expanding the qualitative and quantitative information onpopulation dynamics. Thus, the knowledge acquired from novel technologies is a valuable asset for the breeding process and to better understand the population dynamics, their properties, and analysis methods.

Morphological, Transcriptomic and Hormonal Characterization of Trimonoecious and Subandroecious Pumpkin (Cucurbita maxima) Suggests Important Roles of Ethylene in Sex Expression

Sex expression is a complex process, and in-depth knowledge of its mechanism in pumpkin is important. In this study, young shoot apices at the one-true-leaf stage and 10-leaf stage in Cucurbita maxima trimonoecious line ‘2013–12’ and subandroecious line ‘9–6’ were collected as materials, and transcriptome sequencing was performed using an Illumina HiSeq^TM 2000 System. 496 up-regulated genes and 375 down-regulated genes were identified between shoot apices containing mostly male flower buds and only female flower buds. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the differentially expressed genes were mainly enriched in the ethylene and auxin synthesis and signal transduction pathways. In addition, shoot apices at the 4-leaf stage were treated with the ethylene-releasing agent 2-chloroethylphosphonic acid (Ethrel), aminoethoxyvinyl glycine (AVG), AgNO₃ and indoleacetic acid (IAA). The number of female flowers up to node 20 on the main stem of ‘2013–12’ increased significantly after Ethrel and IAA treatment and decreased significantly after AVG and AgNO₃ treatment. The female flowers in ‘9–6’ showed slight changes after treatment with the exogenous chemicals. The expression of key genes in ethylene synthesis and signal transduction (CmaACS7, CmaACO1, CmaETR1 and CmaEIN3) was determined using quantitative RT-PCR, and the expression of these four genes was positively correlated with the number of female flowers in ‘2013–12’. The variations in gene expression, especially that of CmaACS7, after chemical treatment were small in ‘9–6’. From stage 1 (S1) to stage 7 (S7) of flower development, the expression of CmaACS7 in the stamen was much lower than that in the ovary, stigma and style. These transcriptome data and chemical treatment results indicated that IAA might affect pumpkin sex expression by inducing CmaACS7 expression and indirectly affecting ethylene production, and the ethylene synthesis and signal transduction pathways play crucial roles in pumpkin flower sex expression. A possible reason for the differences in sex expression between pumpkin lines ‘2013–12’ and ‘9–6’ was proposed based on the key gene expression. Overall, these transcriptome data and chemical treatment results suggest important roles for ethylene in pumpkin sex expression.

Application of high-throughput amplicon sequencing-based SSR genotyping in genetic background screening

Background

Host genetic backgrounds affect gene functions. The genetic backgrounds of genetically engineered organisms must be identified to confirm their genetic backgrounds identity with those of recipients. Marker-assisted backcrossing (MAB), transgenesis and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (CRISPR/Cas9) editing are three commonly used genetic engineering techniques. However, methods for genetic background screening between genetically engineered organisms and corresponding recipients suffer from low efficiency, low accuracy or high cost.

Results

Here, we improved our previously reported AmpSeq-SSR method, an amplicon sequencing-based simple sequence repeat (SSR) genotyping method, by selecting SSR loci with high polymorphism among varieties. Ultimately, a set of 396 SSRs was generated and applied to evaluate the genetic backgrounds identity between rice lines developed through MAB, transgenesis, and CRISPR/Cas9 editing and the respective recipient rice. We discovered that the percentage of different SSRs between the MAB-developed rice line and its recipient was as high as 23.5%. In contrast, only 0.8% of SSRs were different between the CRISPR/Cas9-system-mediated rice line and its recipient, while no SSRs showed different genotypes between the transgenic rice line and its recipient. Furthermore, most differential SSRs induced by MAB technology were located in non-coding regions (62.9%), followed by untranslated regions (21.0%) and coding regions (16.1%). Trinucleotide repeats were the most prevalent type of altered SSR. Most importantly, all altered SSRs located in coding regions were trinucleotide repeats.

Conclusions

This method is not only useful for the background evaluation of genetic resources but also expands our understanding of the unintended effects of different genetic engineering techniques. While the work we present focused on rice, this method can be readily extended to other organisms.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5800-4) contains supplementary material, which is available to authorized users.

Quantitative trait loci mapping for the shear force value in breast muscle of F2 chickens

The shear force value is one of the major traits that determine meat quality. In the present study, we performed QTL analysis for chicken breast muscle shear force value at 7 wk of age using 545 single nucleotide polymorphism (SNP) markers developed via restriction-site associated DNA sequencing (RAD-seq). An F2 resource family was generated by mating Oh-Shamo, a native Japanese chicken breed, and the White Plymouth Rock chicken breed. A total of 215 F2 birds were produced. Simple interval mapping revealed one significant main-effect QTL between 6.28 and 8.10 Mb SNPs on the chromosome Z with a logarithm of odds score of 5.53 at the genome-wide 5% level. At this QTL, the confidence interval, phenotypic variance explained, and additive effect were 26 cM, 12.24%, and -0.31 in males and -0.34 in females, respectively. No QTL with epistatic interaction effects were detected. To our knowledge, this is the first report on a QTL affecting the shear force value in the chicken breast muscle, using SNP markers derived from RAD-seq.

Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits

Cucurbita pepo contains two cultivated subspecies, each of which encompasses four fruit-shape morphotypes (cultivar groups). The Pumpkin, Vegetable Marrow, Cocozelle, and Zucchini Groups are of subsp. pepo and the Acorn, Crookneck, Scallop, and Straightneck Groups are of subsp. ovifera. Recently, a de novo assembly of the C. pepo subsp. pepo Zucchini genome was published, providing insights into its evolution. To expand our knowledge of evolutionary processes within C. pepo and to identify variants associated with particular morphotypes, we performed whole-genome resequencing of seven of these eight C. pepo morphotypes. We report for the first time whole-genome resequencing of the four subsp. pepo (Pumpkin, Vegetable Marrow, Cocozelle, green Zucchini, and yellow Zucchini) morphotypes and three of the subsp. ovifera (Acorn, Crookneck, and Scallop) morphotypes. A high-depth resequencing approach was followed, using the BGISEQ-500 platform that enables the identification of rare variants, with an average of 33.5X. Approximately 94.5% of the clean reads were mapped against the reference Zucchini genome. In total, 3,823,977 high confidence single-nucleotide polymorphisms (SNPs) were identified. Within each accession, SNPs varied from 636,918 in green Zucchini to 2,656,513 in Crookneck, and were distributed homogeneously along the chromosomes. Clear differences between subspecies pepo and ovifera in genetic variation and linkage disequilibrium are highlighted. In fact, comparison between subspecies pepo and ovifera indicated 5710 genes (22.5%) with Fst > 0.80 and 1059 genes (4.1%) with Fst = 1.00 as potential candidate genes that were fixed during the independent evolution and domestication of the two subspecies. Linkage disequilibrium was greater in subsp. ovifera than in subsp. pepo, perhaps reflective of the earlier differentiation of morphotypes within subsp. ovifera. Some morphotype-specific genes have been localized. Our results offer new clues that may provide an improved understanding of the underlying genomic regions involved in the independent evolution and domestication of the two subspecies. Comparisons among SNPs unique to particular subspecies or morphotypes may provide candidate genes responsible for traits of high economic importance.

Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita maxima

Seed traits are agronomically important for Cucurbita breeding, but the genes controlling seed size, seed weight and seed number have not been mapped in Cucurbita maxima (C. maxima). In this study, 100 F₂ individual derived from two parental lines, "2013-12" and "9-6", were applied to construct a 3,376.87-cM genetic map containing 20 linkage groups (LGs) with an average genetic distance of 0.47 cM using a total of 8,406 specific length amplified fragment (SLAF) markers in C. maxima. Ten quantitative trait loci (QTLs) of seed width (SW), seed length (SL) and hundred-seed weight (HSW) were identified using the composite interval mapping (CIM) method. The QTLs affecting SW, SL and HSW explained a maximum of 38.6%, 28.9% and 17.2% of the phenotypic variation and were detected in LG6, LG6 and LG17, respectively. To validate these results, an additional 150 F₂ individuals were used for QTL mapping of SW and SL with cleaved amplified polymorphic sequence (CAPS) markers. We found that two major QTLs, SL6-1 and SW6-1, could be detected in both SLAF-seq and CAPS markers in an overlapped region. Based on gene annotation and non-synonymous single-nucleotide polymorphisms (SNPs) in the major SWand SL-associated regions, we found that two genes encoding a VQ motif and an E3 ubiquitin-protein ligase may be candidate genes influencing SL, while an F-box and leucinerich repeat (LRR) domain-containing protein is the potential regulator for SW in C. maxima. This study provides the first high-density linkage map of C. maxima using SNPs developed by SLAF-seq technology, which is a powerful tool for associated mapping of important agronomic traits, map-based gene cloning and marker-assisted selection (MAS)-based breeding in C. maxima.

A High-Density EST-SSR-Based Genetic Map and QTL Analysis of Dwarf Trait in Cucurbita pepo L

As one of the earliest domesticated species, Cucurbita pepo (including squash and pumpkin) is rich in phenotypic polymorphism and has huge economic value. In this research, using 1660 expressed sequence tags-simple sequence repeats (EST-SSRs) and 632 genomic simple sequence repeats (gSSRs), we constructed the highest-density EST-SSR-based genetic map in Cucurbita genus, which spanned 2199.1 cM in total and harbored 623 loci distributed in 20 linkage groups. Using this map as a bridge, the two previous gSSR maps were integrated by common gSSRs and the corresponding relationships around chromosomes in three sets of genomes were also collated. Meanwhile, one large segmental inversion that existed between our map and the C. pepo genome was detected. Furthermore, three Quantitative Trait Loci (QTLs) of the dwarf trait (gibberellin-sensitive dwarf type) in C. pepo were located, and the candidate region that covered the major QTL spanned 1.39 Mb, which harbored a predicted gibberellin 2-β-oxidase gene. Considering the rich phenotypic polymorphism, the important economic value in the Cucurbita genus species and several advantages of the SSR marker were identified; thus, this high-density EST-SSR-based genetic map will be useful in Pumpkin and Squash breeding work in the future.

De Novo Transcriptome Assembly of Cucurbita Pepo L. Leaf Tissue Infested by Aphis Gossypii

Zucchini (Cucurbita pepo L.), extensively cultivated in temperate areas, belongs to the Cucurbitaceae family and it is a species with great economic value. One major threat related to zucchini cultivation is the damage imposed by the cotton/melon aphid Aphis gossypii Glover (Homoptera: Aphididae). We performed RNA-sequencing on cultivar “San Pasquale” leaves, uninfested and infested by A. gossypii, that were collected at three time points (24, 48, and 96 h post infestation). Then, we combined all high-quality reads for de novo assembly of the transcriptome. This resource was primarily established to be used as a reference for gene expression studies in order to investigate the transcriptome reprogramming of zucchini plants following aphid infestation. In addition, raw reads will be valuable for new experiments based on the latest bioinformatic tools and analytical approaches. The assembled transcripts will serve as an important reference for sequence-based studies and for primer design. Both datasets can be used to support/improve the prediction of protein-coding genes in the zucchini genome, which has been recently released into the public domain.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.

High-density genetic mapping of a major QTL for resistance to multiple races of loose smut in a tetraploid wheat cross

Loose smut, caused by Ustilago tritici (Pers.) Rostr., is a systemic disease of tetraploid durum wheat (Triticum turgidum L.). Loose smut can be economically controlled by growing resistant varieties, making it important to find and deploy new sources of resistance. Blackbird, a variety of T. turgidum L. subsp. carthlicum (Nevski) A. Love & D. Love, carries a high level of resistance to loose smut. Blackbird was crossed with the loose smut susceptible durum cultivar Strongfield to produce a doubled haploid (DH) mapping population. The parents and progenies were inoculated with U. tritici races T26, T32 and T33 individually and as a mixture at Swift Current, Canada in 2011 and 2012 and loose smut incidence (LSI) was assessed. Genotyping of the DH population and parents using an Infinium iSelect 90K single nucleotide polymorphism (SNP) array identified 12,952 polymorphic SNPs. The SNPs and 426 SSRs (previously genotyped in the same population) were mapped to 16 linkage groups spanning 3008.4 cM at an average inter-marker space of 0.2 cM in a high-density genetic map. Composite interval mapping analysis revealed three significant quantitative trait loci (QTL) for loose smut resistance on chromosomes 3A, 6B and 7A. The loose smut resistance QTL on 6B (QUt.spa-6B.2) and 7A (QUt.spa-7A.2) were derived from Blackbird. Strongfield contributed the minor QTL on 3A (QUt.spa-3A.2). The resistance on 6B was a stable major QTL effective against all individual races and the mixture of the three races; it explained up to 74% of the phenotypic variation. This study is the first attempt in durum wheat to identify and map loose smut resistance QTL using a high-density genetic map. The QTL QUt.spa-6B.2 would be an effective source for breeding resistance to multiple races of the loose smut pathogen because it provides near-complete broad resistance to the predominant virulence on the Canadian prairies.

Karyotype Stability and Unbiased Fractionation in the Paleo-Allotetraploid Cucurbita Genomes

The Cucurbita genus contains several economically important species in the Cucurbitaceae family. Here, we report high-quality genome sequences of C. maxima and C. moschata and provide evidence supporting an allotetraploidization event in Cucurbita. We are able to partition the genome into two homoeologous subgenomes based on different genetic distances to melon, cucumber, and watermelon in the Benincaseae tribe. We estimate that the two diploid progenitors successively diverged from Benincaseae around 31 and 26 million years ago (Mya), respectively, and the allotetraploidization happened at some point between 26 Mya and 3 Mya, the estimated date when C. maxima and C. moschata diverged. The subgenomes have largely maintained the chromosome structures of their diploid progenitors. Such long-term karyotype stability after polyploidization has not been commonly observed in plant polyploids. The two subgenomes have retained similar numbers of genes, and neither subgenome is globally dominant in gene expression. Allele-specific expression analysis in the C. maxima × C. moschata interspecific F₁ hybrid and their two parents indicates the predominance of trans-regulatory effects underlying expression divergence of the parents, and detects transgressive gene expression changes in the hybrid correlated with heterosis in important agronomic traits. Our study provides insights into polyploid genome evolution and valuable resources for genetic improvement of cucurbit crops.

A high-density linkage map and QTL mapping of fruit-related traits in pumpkin (Cucurbita moschata Duch.)

Pumpkin (Cucurbita moschata) is an economically worldwide crop. Few quantitative trait loci (QTLs) were reported previously due to the lack of genomic and genetic resources. In this study, a high-density linkage map of C. moschata was structured by double-digest restriction site-associated DNA sequencing, using 200 F2 individuals of CMO-1 × CMO-97. By filtering 74,899 SNPs, a total of 3,470 high quality SNP markers were assigned to the map spanning a total genetic distance of 3087.03 cM on 20 linkage groups (LGs) with an average genetic distance of 0.89 cM. Based on this map, both pericarp color and strip were fined mapped to a novel single locus on LG8 in the same region of 0.31 cM with phenotypic variance explained (PVE) of 93.6% and 90.2%, respectively. QTL analysis was also performed on carotenoids, sugars, tuberculate fruit, fruit diameter, thickness and chamber width with a total of 12 traits. 29 QTLs distributed in 9 LGs were detected with PVE from 9.6% to 28.6%. It was the first high-density linkage SNP map for C. moschata which was proved to be a valuable tool for gene or QTL mapping. This information will serve as significant basis for map-based gene cloning, draft genome assembling and molecular breeding.

Resistance to tomato leaf curl New Delhi virus in melon is controlled by a major QTL located in chromosome 11

Identification of three genomic regions and underlying candidate genes controlling the high level of resistance to ToLCNDV derived from a wild melon. SNP markers appropriated for MAS management of resistance. Tomato leaf curl New Delhi virus (ToLCNDV) is a bipartite begomovirus that severely affects melon crop (Cucumis melo) in the main production areas of Spain since 2012. In this work, we evaluated the degree of resistance of four accessions (two belonging to the subsp. agrestis var. momordica and two to the wild agrestis group) and their corresponding hybrids with a susceptible commercial melon belonging to the subsp. melo (Piel de Sapo, PS). The analysis using quantitative PCR (qPCR) allowed us to select one wild agrestis genotype (WM-7) with a high level of resistance and use it to construct segregating populations (F ₂ and backcrosses). These populations were phenotyped for symptom severity and virus content using qPCR, and genotyped with different sets of SNP markers. Phenotyping and genotyping results in the F ₂ and BC1s populations derived from the WM-7 × PS cross were used for QTL analysis. Three genomic regions controlling resistance to ToLCNDV were found, one major locus in chromosome 11 and two additional regions in chromosomes 12 and 2. The highest level of resistance (no or mild symptoms and very low viral titer) was obtained with the homozygous WM-7WM-7 genotype at the major QTL in chromosome 11, even with PSPS genotypes at the other two loci. The resistance derived from WM-7 is useful to develop new melon cultivars and the linked SNPs selected in this paper will be highly useful in marker-assisted breeding for ToLCNDV resistance in melon.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

Background

Genetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.

Results

An F₈ population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F₈ individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F₈ phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.

Conclusions

To our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

An SNP-based saturated genetic map and QTL analysis of fruit-related traits in Zucchini using Genotyping-by-sequencing

Background

Cucurbita pepo is a cucurbit with growing economic importance worldwide. Zucchini morphotype is the most important within this highly variable species. Recently, transcriptome and Simple Sequence Repeat (SSR)- and Single Nucleotide Polymorphism (SNP)-based medium density maps have been reported, however further genomic tools are needed for efficient molecular breeding in the species. Our objective is to combine currently available complete transcriptomes and the Zucchini genome sequence with high throughput genotyping methods, mapping population development and extensive phenotyping to facilitate the advance of genomic research in this species.

Results

We report the Genotyping-by-sequencing analysis of a RIL population developed from the inter subspecific cross Zucchini x Scallop (ssp. pepo x ssp. ovifera). Several thousands of SNP markers were identified and genotyped, followed by the construction of a high-density linkage map based on 7,718 SNPs (average of 386 markers/linkage group) covering 2,817.6 cM of the whole genome, which is a great improvement with respect to previous maps. A QTL analysis was performed using phenotypic data obtained from the RIL population from three environments. In total, 48 consistent QTLs for vine, flowering and fruit quality traits were detected on the basis of a multiple-environment analysis, distributed in 33 independent positions in 15 LGs, and each QTL explained 1.5–62.9% of the phenotypic variance. Eight major QTLs, which could explain greater than 20% of the phenotypic variation were detected and the underlying candidate genes identified.

Conclusions

Here we report the first SNP saturated map in the species, anchored to the physical map. Additionally, several consistent QTLs related to early flowering, fruit shape and length, and rind and flesh color are reported as well as candidate genes for them. This information will enhance molecular breeding in C. pepo and will assist the gene cloning underlying the studied QTLs, helping to reveal the genetic basis of the studied processes in squash.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3439-y) contains supplementary material, which is available to authorized users.

Dicer-Like 4 Is Involved in Restricting the Systemic Movement of Zucchini yellow mosaic virus in Nicotiana benthamiana

Zucchini yellow mosaic virus (ZYMV) induces serious diseases in cucurbits. To create a tool to screen for resistance genes, we cloned a wild ZYMV isolate and inserted the visual marker Rosea1 to obtain recombinant clone ZYMV-Ros1. While in some plant-virus combinations Rosea1 induces accumulation of anthocyanins in infected tissues, ZYMV-Ros1 infection of cucurbits did not lead to detectable anthocyanin accumulation. However, the recombinant virus did induce dark red pigmentation in infected tissues of the model plant Nicotiana benthamiana. In this species, ZYMV-Ros1 multiplied efficiently in local inoculated tissue but only a few progeny particles established infection foci in upper leaves. We used this system to analyze the roles of Dicer-like (DCL) genes, core components of plant antiviral RNA silencing pathways, in ZYMV infection. ZYMV-Ros1 local replication was not significantly affected in single DCL knockdown lines nor in double DCL2/4 and triple DCL2/3/4 knockdown lines. ZYMV-Ros1 systemic accumulation was not affected in knockdown lines DCL1, DCL2, and DCL3. However in DCL4 and also in DCL2/4 and DCL2/3/4 knockdown lines, ZYMV-Ros1 systemic accumulation dramatically increased, which highlights the key role of DCL4 in restricting virus systemic movement. The effect of DCL4 on ZYMV systemic movement was confirmed with a wild-type version of the virus.

Cultivar-Based Introgression Mapping Reveals Wild Species-Derived Pm-0, the Major Powdery Mildew Resistance Locus in Squash

Powdery mildew is a major fungal disease on squash and pumpkin (Cucurbita spp.) in the US and throughout the world. Genetic resistance to the disease is not known to occur naturally within Cucurbita pepo and only infrequently in Cucurbita moschata, but has been achieved in both species through the introgression of a major resistance gene from the wild species Cucurbita okeechobeensis subsp. martinezii. At present, this gene, Pm-0, is used extensively in breeding, and is found in nearly all powdery mildew-resistant C. pepo and C. moschata commercial cultivars. In this study, we mapped C. okeechobeensis subsp. martinezii-derived single nucleotide polymorphism (SNP) alleles in a set of taxonomically and morphologically diverse and resistant C. pepo and C. moschata cultivars bred at Cornell University that, by common possession of Pm-0, form a shared-trait introgression panel. High marker density was achieved using genotyping-by-sequencing, which yielded over 50,000 de novo SNP markers in each of the three Cucurbita species genotyped. A single 516.4 kb wild-derived introgression was present in all of the resistant cultivars and absent in a diverse set of heirlooms that predated the Pm-0 introgression. The contribution of this interval to powdery mildew resistance was confirmed by association mapping in a C. pepo cultivar panel that included the Cornell lines, heirlooms, and 68 additional C. pepo cultivars and with an independent F2 population derived from C. okeechobeensis subsp. martinezii x C. moschata. The interval was refined to a final candidate interval of 76.4 kb and CAPS markers were developed inside this interval to facilitate marker-assisted selection.

Dynamics of sex expression and chromosome diversity in Cucurbitaceae: a story in the making

The family Cucurbitaceae showcases a wide range of sexual phenotypes being variedly regulated by biological and environmental factors. In the present context, we have tried to assemble reports of cytogenetic investigations carried out in cucurbits accompanied by information on sex expression diversities and chromosomal or molecular basis of sex determination in dioecious (or other sexual types, if reported) taxa known so far. Most of the Cucurbitaceae tribes have mixed sexual phenotypes with varying range of chromosome numbers and hence, ancestral conditions become difficult to probe. Occurrence of polyploidy is rare in the family and has no influence on sexual traits. The sex determination mechanisms have been elucidated in some well-studied taxa like Bryonia,Coccinia and Cucumis showing interplay of genic, biochemical, developmental and sometimes chromosomal determinants. Substantial knowledge about genic and molecular sex differentiation has been obtained for genera like Momordica, Cucurbita and Trichosanthes. The detailed information on sex determination schemes, genomic sequences and molecular phylogenetic relationships facilitate further comprehensive investigations in the tribe Bryonieae. The discovery of organ identity genes and sex-specific sequences regulating sexual behaviour in Coccinia,Cucumis and Cucurbita opens up opportunities of relevant investigations to answer yet unaddressed questions pertaining to floral unisexuality, dioecy and chromosome evolution in the family. The present discussion brings the genera in light, previously recognized under subfamily Nhandiroboideae, where the study of chromosome cytology and sex determination mechanisms can simplify our understanding of sex expression pathways and its phylogenetic impacts.

Molecular Response of Crop Plants to Engineered Nanomaterials

Functional toxicology has enabled the identification of genes involved in conferring tolerance and sensitivity to engineered nanomaterial (ENM) exposure in the model plant Arabidopsis thaliana (L.) Heynh. Several genes were found to be involved in metabolic functions, stress response, transport, protein synthesis, and DNA repair. Consequently, analysis of physiological parameters, metal content (through ICP-MS quantification), and gene expression (by RT-qPCR) of A. thaliana orthologue genes were performed across different plant species of agronomic interest to highlight putative biomarkers of exposure and effect related to ENMs. This approach led to the identification of molecular markers in Solanum lycopersicum L. and Cucurbita pepo L. (tomato and zucchini) that might not only indicate exposure to ENMs (CuO, CeO2, and La2O3) but also provide mechanistic insight into response to these materials. Through Gene Ontology (GO) analysis, the target genes were mapped in complex interatomic networks representing molecular pathways, cellular components, and biological processes involved in ENM response. The transcriptional response of 38 (out of 204) candidate genes studied varied according to particle type, size, and plant species. Importantly, some of the genes studied showed potential as biomarkers of ENM exposure and effect and may be useful for risk assessment in foods and in the environment.

Angolan vegetable crops have unique genotypes of potential value for future breeding programmes

Abstract A survey was carried out in Angola with the aim of collecting vegetable crops. Collecting expeditions were conducted in Kwanza-Sul, Benguela, Huíla and Namibe Provinces and a total of 80 accessions belonging to 22 species was collected from farmers and local markets. Species belonging to the Solanaceae (37 accessions) and Cucurbitaceae (36 accessions) families were the most frequently found with pepper and eggplant being the predominant solanaceous crops collected. Peppers were sold in local markets as a mixture of different types, even different species: Capsicum chinense, C. baccatum, C. frutescens and C. pubescens. Most of the eggplant accessions collected belonged to Solanum aethiopicum L. Gilo Group, the so-called ‘scarlet eggplant’. Cucurbita genus was better represented than the other cucurbit crops. A high morphological variation was present in the Cucurbita maxima and C. moschata accessions. A set of 22 Cucurbita accessions from Angola, along with 32 Cucurbita controls from a wide range of origins, was cultivated in Valencia, Spain and characterised based on morphology and molecularity using a set of 15 microsatellite markers. A strong dependence on latitude was found in most of the accessions and as a result, many accessions did not set fruit. The molecular analysis showed high molecular variability and uniqueness in the collected accessions, as shown by their segregation from the set of global controls. In summary, the material collected is quite valuable because of its uniqueness and the potential of the breeding characteristics it possesses.

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 high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.)

Background

Pumpkin (Cucurbita maxima Duch.) is an economically important crop belonging to the Cucurbitaceae family. However, very few genomic and genetic resources are available for this species. As part of our ongoing efforts to sequence the pumpkin genome, high-density genetic map is essential for anchoring and orienting the assembled scaffolds. In addition, a saturated genetic map can facilitate quantitative trait locus (QTL) mapping.

Results

A set of 186 F₂ plants derived from the cross of pumpkin inbred lines Rimu and SQ026 were genotyped using the genotyping-by-sequencing approach. Using the SNPs we identified, a high-density genetic map containing 458 bin-markers was constructed, spanning a total genetic distance of 2,566.8 cM across the 20 linkage groups of C. maxima with a mean marker density of 5.60 cM. Using this map we were able to anchor 58 assembled scaffolds that covered about 194.5 Mb (71.7 %) of the 271.4 Mb assembled pumpkin genome, of which 44 (183.0 Mb; 67.4 %) were oriented. Furthermore, the high-density genetic map was used to identify genomic regions highly associated with an important agronomic trait, dwarf vine. Three QTLs on linkage groups (LGs) 1, 3 and 4, respectively, were recovered. One QTL, qCmB2, which was located in an interval of 0.42 Mb on LG 3, explained 21.4 % phenotypic variations. Within qCmB2, one gene, Cma_004516, encoding the gibberellin (GA) 20-oxidase in the GA biosynthesis pathway, had a 1249-bp deletion in its promoter in bush type lines, and its expression level was significantly increased during the vine growth and higher in vine type lines than bush type lines, supporting Cma_004516 as a possible candidate gene controlling vine growth in pumpkin.

Conclusions

A high-density pumpkin genetic map was constructed, which was used to successfully anchor and orient the assembled genome scaffolds, and to identify QTLs highly associated with pumpkin vine length. The map provided a valuable resource for gene cloning and marker assisted breeding in pumpkin and other related species. The identified vine length QTLs would help to dissect the underlying molecular basis regulating pumpkin vine growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2312-8) contains supplementary material, which is available to authorized users.

RNA-Seq Uncovers SNPs and Alternative Splicing Events in Asian Lotus (Nelumbo nucifera)

RNA-Seq is an efficient way to comprehensively identify single nucleotide polymorphisms (SNPs) and alternative splicing (AS) events from the expressed genes. In this study, we conducted transcriptome sequencing of four Asian lotus (Nelumbo nucifera) cultivars using Illumina HiSeq2000 platform to identify SNPs and AS events in lotus. A total of 505 million pair-end RNA-Seq reads were generated from four cultivars, of which 86% were mapped to the lotus reference genome. Using the four sets of data together, a total of 357,689 putative SNPs were identified with an average density of one SNP per 2.2 kb. These SNPs were located in 1,253 scaffolds and 15,016 expressed genes. A/G and C/T were the two major types of SNPs in the Asian lotus transcriptome. In parallel, a total of 177,540 AS events were detected in the four cultivars and were distributed in 64% of the expressed genes of lotus. The predominant type of AS events was alternative 5’ first exon, which accounted for 41.2% of all the observed AS events, and exon skipping only accounted for 4.3% of all AS. Gene Ontology analysis was conducted to analyze the function of the genes containing SNPs and AS events. Validation of selected SNPs and AS events revealed that 74% of SNPs and 80% of AS events were reliable, which indicates that RNA-Seq is an efficient approach to uncover gene-associated SNPs and AS events. A large number of SNPs and AS events identified in our study will facilitate further genetic and functional genomics research in lotus.

Construction of a high-density DArTseq SNP-based genetic map and identification of genomic regions with segregation distortion in a genetic population derived from a cross between feral and cultivated-type watermelon

Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is an economically important vegetable crop grown extensively worldwide. To facilitate the identification of agronomically important traits and provide new information for genetic and genomic research on this species, a high-density genetic linkage map of watermelon was constructed using an F2 population derived from a cross between elite watermelon cultivar K3 and wild watermelon germplasm PI 189225. Based on a sliding window approach, a total of 1,161 bin markers representing 3,465 SNP markers were mapped onto 11 linkage groups corresponding to the chromosome pair number of watermelon. The total length of the genetic map is 1,099.2 cM, with an average distance between bins of 1.0 cM. The number of markers in each chromosome varies from 62 in chromosome 07 to 160 in chromosome 05. The length of individual chromosomes ranged between 61.8 cM for chromosome 07 and 140.2 cM for chromosome 05. A total of 616 SNP bin markers showed significant (P < 0.05) segregation distortion across all 11 chromosomes, and 513 (83.3 %) of these distorted loci showed distortion in favor of the elite watermelon cultivar K3 allele and 103 were skewed toward PI 189225. The number of SNPs and InDels per Mb varied considerably across the segregation distorted regions (SDRs) on each chromosome, and a mixture of dense and sparse SNPs and InDel SDRs coexisted on some chromosomes suggesting that SDRs were randomly distributed throughout the genome. Recombination rates varied greatly among each chromosome, from 2.0 to 4.2 centimorgans per megabase (cM/Mb). An inconsistency was found between the genetic and physical positions on the map for a segment on chromosome 11. The high-density genetic map described in the present study will facilitate fine mapping of quantitative trait loci, the identification of candidate genes, map-based cloning, as well as marker-assisted selection (MAS) in watermelon breeding programs.

An acorn squash (Cucurbita pepo ssp. ovifera) fruit and seed transcriptome as a resource for the study of fruit traits in Cucurbita

Acorn squash (Cucurbita pepo) is an iconic fall vegetable in the United States, known for its unique fruit shape and also prized for its culinary properties. Little is known about the metabolism that underlies the development of fruit quality attributes such as color, sweetness, texture and nutritional qualities in acorn squash, or any other winter squash grown worldwide. To provide insight into winter squash fruit and seed development and add to the genomic resources in the Cucurbita genus, RNA sequencing was used to generate an acorn squash fruit and seed transcriptome from the cultivar Sweet REBA at critical points throughout fruit development. 141 838 600 high-quality paired-end Illumina reads were assembled into 55 949 unigenes. 85% of unigenes with predicted open reading frames had homology with previously identified genes and over 62% could be functionally annotated. Comparison with the watermelon and cucumber genomes provided confirmation that the unigenes are full-length and comprehensive, covering an average of 90% of the coding sequence of their homologs and 72% of the cucumber and watermelon exomes. Key candidate genes associated with carotenoid and carbohydrate metabolism were identified toward a resource for winter squash fruit quality trait dissection. This transcriptome represents a major advance in C. pepo genomics, providing significant new sequence information and revealing the repertoire of genes expressed throughout winter squash fruit and seed development. Future studies on the genetic basis of fruit quality and future breeding efforts will be enhanced by tools and insights developed from this resource.

First TILLING platform in Cucurbita pepo: a new mutant resource for gene function and crop improvement

Although the availability of genetic and genomic resources for Cucurbita pepo has increased significantly, functional genomic resources are still limited for this crop. In this direction, we have developed a high throughput reverse genetic tool: the first TILLING (Targeting Induced Local Lesions IN Genomes) resource for this species. Additionally, we have used this resource to demonstrate that the previous EMS mutant population we developed has the highest mutation density compared with other cucurbits mutant populations. The overall mutation density in this first C. pepo TILLING platform was estimated to be 1/133 Kb by screening five additional genes. In total, 58 mutations confirmed by sequencing were identified in the five targeted genes, thirteen of which were predicted to have an impact on the function of the protein. The genotype/phenotype correlation was studied in a peroxidase gene, revealing that the phenotype of seedling homozygous for one of the isolated mutant alleles was albino. These results indicate that the TILLING approach in this species was successful at providing new mutations and can address the major challenge of linking sequence information to biological function and also the identification of novel variation for crop breeding.

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.

Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

The genetics and molecular biology of sesame has only recently begun to be studied even though sesame is an important oil seed crop. A high-density genetic map for sesame has not been published yet due to a lack of sufficient molecular markers. Specific length amplified fragment sequencing (SLAF-seq) is a recently developed high-resolution strategy for large-scale de novo SNP discovery and genotyping. SLAF-seq was employed in this study to obtain sufficient markers to construct a high-density genetic map for sesame.

RESULTS

In total, 28.21 Gb of data containing 201,488,285 pair-end reads was obtained after sequencing. The average coverage for each SLAF marker was 23.48-fold in the male parent, 23.38-fold in the female parent, and 14.46-fold average in each F2 individual. In total, 71,793 high-quality SLAFs were detected of which 3,673 SLAFs were polymorphic and 1,272 of the polymorphic markers met the requirements for use in the construction of a genetic map. The final map included 1,233 markers on the 15 linkage groups (LGs) and was 1,474.87 cM in length with an average distance of 1.20 cM between adjacent markers. To our knowledge, this map is the densest genetic linkage map to date for sesame. 'SNP_only' markers accounted for 87.51% of the markers on the map. A total of 205 markers on the map showed significant (P < 0.05) segregation distortion.

CONCLUSIONS

We report here the first high-density genetic map for sesame. The map was constructed using an F2 population and the SLAF-seq approach, which allowed the efficient development of a large number of polymorphic markers in a short time. Results of this study will not only provide a platform for gene/QTL fine mapping, map-based gene isolation, and molecular breeding for sesame, but will also serve as a reference for positioning sequence scaffolds on a physical map, to assist in the process of assembling the sesame genome sequence.

Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.)

BACKGROUND

We have identified a kind of parthenocarpy in zucchini squash which is associated with an incomplete andromonoecy, i.e. a partial conversion of female into bisexual flowers. Given that andromonoecy in this and other cucurbit species is caused by a reduction of ethylene production in the female flower, the associated parthenocarpic development of the fruit suggested the involvement of ethylene in fruit set and early fruit development.

RESULTS

We have compared the production of ethylene as well as the expression of 13 ethylene biosynthesis and signalling genes in pollinated and unpollinated ovaries/fruits of two cultivars, one of which is parthenocarpic (Cavili), while the other is non-parthenocarpic (Tosca). In the latter, unpollinated ovaries show an induction of ethylene biosynthesis and ethylene signal transduction pathway genes three days after anthesis, which is concomitant with the initiation of fruit abortion and senescence. Fruit set and early fruit development in pollinated flowers of both cultivars and unpollinated flowers of Cavili is coupled with low ethylene biosynthesis and signalling, which would also explain the partial andromonoecy in the parthenocarpic genotype. The reduction of ethylene production in the ovary cosegregates with parthenocarpy and partial andromonoecy in the selfing progeny of Cavili. Moreover, the induction of ethylene in anthesis (by ethephon treatments) reduced the percentage of bisexual parthenocarpic flowers in Cavili, while the inhibition of ethylene biosynthesis or response (by AVG and STS treatments) induces not only andromonoecy but also the parthenocarpic development of the fruit in both cultivars.

CONCLUSIONS

Results demonstrate that a reduction of ethylene production or signalling in the zucchini flower is able to induce fruit set and early fruit development, and therefore that ethylene is actively involved in fruit set and early fruit development. Auxin and TIBA treatments, inducing fruit set and early fruit development in this species, also inhibit ethylene production and the expression of ethylene biosynthesis and response genes. A model is presented that discusses the crosstalk between ethylene and auxin in the control of fruit set and early fruit development in zucchini squash.

Carotenogenic gene expression and carotenoid accumulation in three varieties of Cucurbita pepo during fruit development

The control of gene expression is a crucial regulatory mechanism in carotenoid accumulation of fruits and flowers. We investigated the role of transcriptional regulation of nine genes involved in the carotenoid biosynthesis pathway in three varieties of Cucurbita pepo with evident differences in fruit color. The transcriptional levels of the key genes involved in the carotenoid biosynthesis were higher in flower-, leaf-, and fruit skin tissues than flesh tissues. This correlated with higher concentration of carotenoid content in these tissues. The differential expression among the colored and white cultivars detected for some genes, such as LCYe, in combination with other regulatory mechanisms, could explain the large differences found in terms of carotenoid content among the three varieties. These results are a first step to elucidate carotenogenesis in C. pepo and demonstrate that, in general, regulation of the pathway genes is a critical factor that determines the accumulation of these compounds.

SNP genotyping in melons: genetic variation, population structure, and linkage disequilibrium

Novel sequencing technologies were recently used to generate sequences from multiple melon (Cucumis melo L.) genotypes, enabling the in silico identification of large single nucleotide polymorphism (SNP) collections. In order to optimize the use of these markers, SNP validation and large-scale genotyping are necessary. In this paper, we present the first validated design for a genotyping array with 768 SNPs that are evenly distributed throughout the melon genome. This customized Illumina GoldenGate assay was used to genotype a collection of 74 accessions, representing most of the botanical groups of the species. Of the assayed loci, 91 % were successfully genotyped. The array provided a large number of polymorphic SNPs within and across accessions. This set of SNPs detected high levels of variation in accessions from this crop's center of origin as well as from several other areas of melon diversification. Allele distribution throughout the genome revealed regions that distinguished between the two main groups of cultivated accessions (inodorus and cantalupensis). Population structure analysis showed a subdivision into five subpopulations, reflecting the history of the crop. A considerably low level of LD was detected, which decayed rapidly within a few kilobases. Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in melon. Since many of the genotyped accessions are currently being used as the parents of breeding populations in various programs, this set of mapped markers could be used for future mapping and breeding efforts.