Single-nucleotide polymorphisms detected in expressed sequence tags of melon (Cucumis melo L.).

A search was performed for single-nucleotide polymorphisms (SNP) and short insertions-deletions (indels) in 34 melon (Cucumis melo L.) expressed sequence tag (EST) fragments between two distantly related melon genotypes, a group Inodorus 'Piel de sapo' market class breeding line T111 and the Korean accession PI 161375. In total, we studied 15 kb of melon sequence. The average frequency of SNPs between the two genotypes was one every 441 bp. One indel was also found every 1666 bp. Seventy-five percent of the polymorphisms were located in introns and the 3'untranslated regions. On average, there were 1.26 SNPs plus indels per amplicon. We explored three different SNP detection systems to position five of the SNPs in a melon genetic map. Three of the SNPs were mapped using cleaved amplified polymorphic sequence (CAPS) markers, one SNP was mapped using the single primer extension reaction with fluorescent-labelled dideoxynucleotides, and one indel was mapped using polyacrilamide gel electrophoresis separation. The discovery of SNPs based on ESTs and a suitable system for SNP detection has broad potential utility in melon genome mapping.     

Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.)

A set of 118 simple sequence repeat (SSR) markers has been developed in melon from two different sources: genomic libraries (gSSR) and expressed sequence-tag (EST) databases (EST-SSR). Forty-nine percent of the markers showed polymorphism between the Piel de Sapo (PS) and PI161375 melon genotypes used as parents for the mapping populations. Similar polymorphism levels were found in gSSR (51.2%) and EST-SSR (45.5%). Two populations, F₂ and a set of double haploid lines (DHLs), developed from the same parent genotypes were used for map construction. Twenty-three SSRs and 79 restriction fragment length polymorphisms (RFLPs), evenly distributed through the melon genome, were used to anchor the maps of both populations. Ten cucumber SSRs, 41 gSSRs, 16 EST-SSR, three single nucleotide polymorphism (SNP) markers, and the Nsv locus were added in the DHL population. The maps developed in the F₂ and DHL populations were co-linear, with similar lengths, except in linkage groups G1, G9, and G10. There was segregation distortion in a higher proportion of markers in the DHL population compared with the F₂, probably caused by selection during the construction of DHLs through in vitro culture. After map merging, a composite genetic map was obtained including 327 transferable markers: 226 RFLPs, 97 SSRs, three SNPs, and the Nsv locus. The map length is 1,021 cM, distributed in 12 linkage groups, and map density is 3.11 cM/marker. SSR markers alone cover nearly 80% of the map length. This map is proposed as a basis for a framework melon map to be merged with other maps and as an anchor point for map comparison between species of the Cucurbitaceae family.Electronic Supplementary Material Supplementary material is available for this article at     

Development of EST-based SNP and InDel markers and their utilization in tetraploid cotton genetic mapping

Background

Availability of molecular markers has proven to be an efficient tool in facilitating progress in plant breeding, which is particularly important in the case of less researched crops such as cotton. Considering the obvious advantages of single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), expressed sequence tags (ESTs) were analyzed in silico to identify SNPs and InDels in this study, aiming to develop more molecular markers in cotton.

Results

A total of 1,349 EST-based SNP and InDel markers were developed by comparing ESTs between Gossypium hirsutum and G. barbadense, mining G. hirsutum unigenes, and analyzing 3′ untranslated region (3′UTR) sequences. The marker polymorphisms were investigated using the two parents of the mapping population based on the single-strand conformation polymorphism (SSCP) analysis. Of all the markers, 137 (10.16%) were polymorphic, and revealed 142 loci. Linkage analysis using a BC₁ population mapped 133 loci on the 26 chromosomes. Statistical analysis of base variations in SNPs showed that base transitions accounted for 55.78% of the total base variations and gene ontology indicated that cotton genes varied greatly in harboring SNPs ranging from 1.00 to 24.00 SNPs per gene. Sanger sequencing of three randomly selected SNP markers revealed discrepancy between the in silico predicted sequences and the actual sequencing results.

Conclusions

In silico analysis is a double-edged blade to develop EST-SNP/InDel markers. On the one hand, the designed markers can be well used in tetraploid cotton genetic mapping. And it plays a certain role in revealing transition preference and SNP frequency of cotton genes. On the other hand, the developmental efficiency of markers and polymorphism of designed primers are comparatively low.

Electronic supplementary material

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

Use of VeraCode 384-plex assays for watermelon diversity analysis and integrated genetic map of watermelon with single nucleotide polymorphisms and simple sequence repeats

Watermelon (Citrullus lanatus var. lanatus) is one of the most important vegetable crops in the world. Molecular markers have become the tools of choice for resolving watermelon taxonomic relationships and evolution. Increased numbers of single nucleotide polymorphism (SNP) markers together with simple sequence repeat (SSR) markers would be useful for phylogenetic analyses of germplasm accessions and for linkage mapping for marker-assisted breeding with quantitative trait loci and single genes. We aimed to construct a genetic map based on SNPs (generated by Illumina Veracode multiplex assays for genotyping) and SSR markers and evaluate relationships inferred from SNP genotypes between 130 watermelon accessions collected throughout the world. We incorporated 282 markers (232 SNPs and 50 SSRs) into the linkage map. The genetic map consisted of 11 linkage groups spanning 924.72 cM with an average distance of 3.28 cM between markers. Because all of the SNP-containing sequences were assembled with the whole-genome sequence draft for watermelon, chromosome numbers could be readily assigned for all the linkage groups. We found that 134 SNPs were polymorphic in 130 watermelon accessions chosen for diversity studies. The current 384-plex SNP set is a powerful tool for characterizing genetic relatedness and for developing medium-resolution genetic maps.     

Highly efficient genotyping of rice biparental populations by GoldenGate assays based on parental resequencing

Key message

A new time- and cost-effective strategy was developed for medium-density SNP genotyping of rice biparental populations, using GoldenGate assays based on parental resequencing.

Abstract

Since the advent of molecular markers, crop researchers and breeders have dedicated huge amounts of effort to detecting quantitative trait loci (QTL) in biparental populations for genetic analysis and marker-assisted selection (MAS). In this study, we developed a new time- and cost-effective strategy for genotyping a population of progeny from a rice cross using medium-density single nucleotide polymorphisms (SNPs). Using this strategy, 728,362 “high quality” SNPs were identified by resequencing Teqing and Lemont, the parents of the population. We selected 384 informative SNPs that were evenly distributed across the genome for genotyping the biparental population using the Illumina GoldenGate assay. 335 (87.2 %) validated SNPs were used for further genetic analyses. After removing segregation distortion markers, 321 SNPs were used for linkage map construction and QTL mapping. This strategy generated SNP markers distributed more evenly across the genome than previous SSR assays. Taking the GW5 gene that controls grain shape as an example, our strategy provided higher accuracy (0.8 Mb) and significance (LOD 5.5 and 10.1) in QTL mapping than SSR analysis. Our study thus provides a rapid and efficient strategy for genetic studies and QTL mapping using SNP genotyping assays.     

Use of targeted SNP selection for an improved anchoring of the melon (Cucumis melo L.) scaffold genome assembly

Background

The genome of the melon (Cucumis melo L.) double-haploid line DHL92 was recently sequenced, with 87.5 and 80.8% of the scaffold assembly anchored and oriented to the 12 linkage groups, respectively. However, insufficient marker coverage and a lack of recombination left several large, gene rich scaffolds unanchored, and some anchored scaffolds unoriented. To improve the anchoring and orientation of the melon genome assembly, we used resequencing data between the parental lines of DHL92 to develop a new set of SNP markers from unanchored scaffolds.

Results

A high-resolution genetic map composed of 580 SNPs was used to anchor 354.8 Mb of sequence, contained in 141 scaffolds (average size 2.5 Mb) and corresponding to 98.2% of the scaffold assembly, to the 12 melon chromosomes. Over 325.4 Mb (90%) of the assembly was oriented. The genetic map revealed regions of segregation distortion favoring SC alleles as well as recombination suppression regions coinciding with putative centromere, 45S, and 5S rDNA sites. New chromosome-scale pseudomolecules were created by incorporating to the previous v3.5 version an additional 38.3 Mb of anchored sequence representing 1,837 predicted genes contained in 55 scaffolds. Using fluorescent in situ hybridization (FISH) with BACs that produced chromosome-specific signals, melon chromosomes that correspond to the twelve linkage groups were identified, and a standardized karyotype of melon inbred line T111 was developed.

Conclusions

By utilizing resequencing data and targeted SNP selection combined with a large F2 mapping population, we significantly improved the quantity of anchored and oriented melon scaffold genome assembly. Using genome information combined with FISH mapping provided the first cytogenetic map of an inodorus melon type. With these results it was possible to make inferences on melon chromosome structure by relating zones of recombination suppression to centromeres and 45S and 5S heterochromatic regions. This study represents the first steps towards the integration of the high-resolution genetic and cytogenetic maps with the genomic sequence in melon that will provide more information on genome organization and allow for the improvement of the melon genome draft sequence.

Electronic supplementary material

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

In Vitro vs In Silico Detected SNPs for the Development of a Genotyping Array: What Can We Learn from a Non-Model Species?

Background

There is considerable interest in the high-throughput discovery and genotyping of single nucleotide polymorphisms (SNPs) to accelerate genetic mapping and enable association studies. This study provides an assessment of EST-derived and resequencing-derived SNP quality in maritime pine (Pinus pinaster Ait.), a conifer characterized by a huge genome size (∼23.8 Gb/C).

Methodology/Principal Findings

A 384-SNPs GoldenGate genotyping array was built from i/ 184 SNPs originally detected in a set of 40 re-sequenced candidate genes (in vitro SNPs), chosen on the basis of functionality scores, presence of neighboring polymorphisms, minor allele frequencies and linkage disequilibrium and ii/ 200 SNPs screened from ESTs (in silico SNPs) selected based on the number of ESTs used for SNP detection, the SNP minor allele frequency and the quality of SNP flanking sequences. The global success rate of the assay was 66.9%, and a conversion rate (considering only polymorphic SNPs) of 51% was achieved. In vitro SNPs showed significantly higher genotyping-success and conversion rates than in silico SNPs (+11.5% and +18.5%, respectively). The reproducibility was 100%, and the genotyping error rate very low (0.54%, dropping down to 0.06% when removing four SNPs showing elevated error rates).

Conclusions/Significance

This study demonstrates that ESTs provide a resource for SNP identification in non-model species, which do not require any additional bench work and little bio-informatics analysis. However, the time and cost benefits of in silico SNPs are counterbalanced by a lower conversion rate than in vitro SNPs. This drawback is acceptable for population-based experiments, but could be dramatic in experiments involving samples from narrow genetic backgrounds. In addition, we showed that both the visual inspection of genotyping clusters and the estimation of a per SNP error rate should help identify markers that are not suitable to the GoldenGate technology in species characterized by a large and complex genome.     

Genetic diversity among Chinese Hami melon and its relationship with melon germplasm of diverse origins revealed by microsatellite markers

Thick-skinned melon called Hami melon is the most widely cultivated and exported type of melon in China, and mainly grown in Xinjiang province. Here the genetic variation of 64 melon genotypes including 43 Xinjiang Hami melon accessions was analyzed using 36 simple sequence repeat (SSR) markers yielding 145 alleles. The polymorphic information content of SSR markers ranged from 0.09 to 0.83 (average 0.45). Based on the SSR markers, the melon accessions were clustered into 2 major groups (thick and thin-skinned melons). In addition, the sub-cluster analysis based on SSR markers partitioned different botanical groups, even separating similar agronomic trait groups (Xinjiang landraces var. ameri and var. inodorus). SSR analysis showed that 4 SSR markers (CMBR150, CMCTT144, CMBR84 and CMBR12) produced polymorphic bands of different sizes between these two botanical groups. Those four molecular markers might be related to melon fruit maturing time. A considerably low level of genetic diversity was detected in Xinjiang melon accessions. Genetic distances indicated the relatively narrower genetic base but specific taxonomic status of Xinjiang landraces compared with foreign reference accessions.     

Molecular authentication of three Italian melon accessions by ARMS-PCR and ITS1 (internal transcribed spacer 1) secondary structure prediction

Genetic assessment was carried out on three Italian melon accessions by sequence and structural analysis of the internal transcribed spacer 1 (ITS1) from three populations belonging to two Cucumis melo L. varieties (madras and tendral). Alignment of the 18S-5.8S-26S sequences from three melon accessions showed that there were three single-nucleotide polymorphisms (SNPs) and one short insertion-deletion (indel) at the 5'end ITS1. An amplification refractory mutation system (ARMS)-PCR-based analysis was successfully applied to the SNP markers of the ITS1 sequences for the fingerprinting analysis of three melon populations. Secondary structure models for each ITS1 were derived. The prediction of ITS1 RNA secondary structure from each accession was improved by detecting key functional elements shared by all sequences in the alignments. Our results demonstrated that the ITS1secondary structure models can be used to improve the preliminary genetic assessment of the three melon accessions, suggesting a new tool in plant fingerprinting analysis.     

Whole genome resequencing in tomato reveals variation associated with introgression and breeding events

Background

One of the goals of genomics is to identify the genetic loci responsible for variation in phenotypic traits. The completion of the tomato genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of genetic variation present in the tomato genome. Like many self-pollinated crops, cultivated tomato accessions show a low molecular but high phenotypic diversity. Here we describe the whole-genome resequencing of eight accessions (four cherry-type and four large fruited lines) chosen to represent a large range of intra-specific variability and the identification and annotation of novel polymorphisms.

Results

The eight genomes were sequenced using the GAII Illumina platform. Comparison of the sequences with the reference genome yielded more than 4 million single nucleotide polymorphisms (SNPs). This number varied from 80,000 to 1.5 million according to the accessions. Almost 128,000 InDels were detected. The distribution of SNPs and InDels across and within chromosomes was highly heterogeneous revealing introgressions from wild species and the mosaic structure of the genomes of the cherry tomato accessions. In-depth annotation of the polymorphisms identified more than 16,000 unique non-synonymous SNPs. In addition 1,686 putative copy-number variations (CNVs) were identified.

Conclusions

This study represents the first whole genome resequencing experiment in cultivated tomato. Substantial genetic differences exist between the sequenced tomato accessions and the reference sequence. The heterogeneous distribution of the polymorphisms may be related to introgressions that occurred during domestication or breeding. The annotated SNPs, InDels and CNVs identified in this resequencing study will serve as useful genetic tools, and as candidate polymorphisms in the search for phenotype-altering DNA variations.

Electronic supplementary material

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

Construction of a high-density SNP genetic map in flue-cured tobacco based on SLAF-seq

Tobacco (Nicotiana tabacum L., 2n = 48) is an important agronomic crop and model plant. Flue-cured tobacco is the most important type and accounts for approximately 80 % of tobacco production worldwide. The low genetic diversity of flue-cured tobacco impedes the construction of a high-density genetic linkage map using simple sequence repeat (SSR) markers and warrants the exploitation of single nucleotide polymorphic (SNP) markers from genomic regions. In this article, initially using specific locus-amplified fragment sequencing, we discovered 10,891 SNPs that were subsequently used as molecular markers for genetic map construction. Combined with SSR markers, a final high-density genetic map was generated containing 4215 SNPs and 194 SSRs distributed on 24 linkage groups (LGs). The genetic map was 2662.43 cM in length, with an average distance of 0.60 cM between adjacent markers. Furthermore, by mapping the SNP markers to the ancestral genomes of Nicotiana tomentosiformis and Nicotiana sylvestris, a large number of genome rearrangements were identified as occurring after the polyploidization event. Finally, using this novel integrated map and mapping population, two major quantitative trait loci (QTLs) were identified for flue-curing and mapped to the LG6 of tobacco. This is the first report of SNP markers and a SNP-based linkage map being developed in tobacco. The high-density genetic map and QTLs related to tobacco curing will support gene/QTL fine mapping, genome sequence assembly and molecular breeding in tobacco.     

High-throughput single nucleotide polymorphism (SNP) identification and mapping in the sesame (<Emphasis Type="Italic">Sesamum indicum</Emphasis> L.) genome with genotyping by sequencing (GBS) analysis

Sesame (Sesamum indicum L. syn. Sesamum orientale L.) is considered to be the first oil seed crop known to man. Despite its versatile use as an oil seed and a leafy vegetable, sesame is a neglected crop and has not been a subject of molecular genetic research until the last decade. There is thus limited knowledge regarding genome-specific molecular markers that are indispensible for germplasm enhancement, gene identification, and marker-assisted breeding in sesame. In this study, we employed a genotyping by sequencing (GBS) approach to a sesame recombinant inbred line (RIL) population for high-throughput single nucleotide polymorphism (SNP) identification and genotyping. A total of 15,521 SNPs were identified with 14,786 SNPs (95.26 %) located along sesame genome assembly pseudomolecules. By incorporating sesame-specific simple sequence repeat (SSR) markers developed in our previous work, 230.73 megabases (99 %) of sequence from the genome assembly were saturated with markers. This large number of markers will be available for sesame geneticists as a resource for candidate polymorphisms located along the physical chromosomes of sesame. Defining SNP loci in genome assembly sequences provides the flexibility to utilize any genotyping strategy to survey any sesame population. SNPs selected through a high stringency filtering protocol (770 SNPs) for improved map accuracy were used in conjunction with SSR markers (50 SSRs) in linkage analysis, resulting in 13 linkage groups that encompass a total genetic distance of 914 cM with 432 markers (420 SNPs, 12 SSRs). The genetic linkage map constitutes the basis for future work that will involve quantitative trait locus (QTL) analyses of metabolic and agronomic traits in the segregating RIL population.     

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.     

Simple sequence repeats in Cucumis mapping and map merging.

Thirty-four polymorphic simple-sequence repeats (SSRs) were evaluated for length polymorphism in melon (Cucumis melo L.) and cucumber (Cucumis sativus L.). SSR markers were located on three melon maps (18 on the map of 'Vedrantais' and PI 161375, 23 on the map of 'Piel de Sapo' and PI 161375, and 16 on the map of PI 414723 and 'Dulce'). In addition, 14 of the markers were located on the cucumber map of GY14 and PI 183967. SSRs proved to be randomly distributed throughout the melon and cucumber genomes. Mapping of the SSRs in the different maps led to the cross-identification of seven linkage groups in all melon maps. In addition, nine SSRs were common to both melon and cucumber maps. The potential of SSR markers as anchor points for melon-map merging and for comparative mapping with cucumber was demonstrated.     

Insertion-deletion polymorphisms in 3′ regions of maize genes occur frequently and can be used as highly informative genetic markers

Single-nucleotide polymorphisms (SNPs) are the most frequent variations in the genome of any organism. SNP discovery approaches such as resequencing or data mining enable the identification of insertion deletion (indel) polymorphisms. These indels can be treated as biallelic markers and can be utilized for genetic mapping and diagnostics. In this study 655 indels have been identified by resequencing 502 maize (Zea mays) loci across 8 maize inbreds (selected for their high allelic variation). Of these 502 loci, 433 were polymorphic, with indels identified in 215 loci. Of the 655 indels identified, single-nucleotide indels accounted for more than half (54.8%) followed by two- and three-nucleotide indels. A high frequency of 6-base (3.4%) and 8-base (2.3%) indels were also observed. When analysis is restricted to the B73 and Mo17 genotypes, 53% of the loci analyzed contained indels, with 42% having an amplicon size difference. Three novel miniature inverted-repeat transposable element (MITE)-like sequences were identified as insertions near genes. The utility of indels as genetic markers was demonstrated by using indel polymorphisms to map 22 loci in a B73 × Mo17 recombinant inbred population. This paper clearly demonstrates that the resequencing of 3 EST sequence and the discovery and mapping of indel markers will position corresponding expressed genes on the genetic map.     

Construction of a SNP and SSR linkage map in autotetraploid blueberry using genotyping by sequencing

The construction of the first genetic map in autotetraploid blueberry has been made possible by the development of new SNP markers developed using genotyping by sequencing in a mapping population created from a cross between two key highbush blueberry cultivars, Draper × Jewel (Vaccinium corymbosum). The novel SNP markers were supplemented with existing SSR markers to enable the alignment of parental maps.  In total, 1794 single nucleotide polymorphic (SNP) markers and 233 simple sequence repeat (SSR) markers exhibited segregation patterns consistent with a random chromosomal segregation model for meiosis in an autotetraploid. Of these, 700 SNPs and 85 SSRs were utilized for construction of the ‘Draper’ genetic map, and 450 SNPs and 86 SSRs for the ‘Jewel’ map.  The ‘Draper’ map comprises 12  linkage groups (LG), associated with the haploid chromosome number for blueberry, and totals 1621 cM while the ‘Jewel’ map comprises 20 linkage groups totalling 1610 cM. Tentative alignments of the two parental maps have been made on the basis of shared SSR alleles and linkages to double-simplex markers segregating in both parents. Tentative alignments of the two parental maps have been made on the basis of shared SSR alleles and linkages to double-simplex markers segregating in both parents.