Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing
© Wang et al.; licensee BioMed Central Ltd. 2012
Received: 23 February 2012
Accepted: 18 July 2012
Published: 21 August 2012
Genetic mapping and QTL detection are powerful methodologies in plant improvement and breeding. Construction of a high-density and high-quality genetic map would be of great benefit in the production of superior grapes to meet human demand. High throughput and low cost of the recently developed next generation sequencing (NGS) technology have resulted in its wide application in genome research. Sequencing restriction-site associated DNA (RAD) might be an efficient strategy to simplify genotyping. Combining NGS with RAD has proven to be powerful for single nucleotide polymorphism (SNP) marker development.
An F1 population of 100 individual plants was developed. In-silico digestion-site prediction was used to select an appropriate restriction enzyme for construction of a RAD sequencing library. Next generation RAD sequencing was applied to genotype the F1 population and its parents. Applying a cluster strategy for SNP modulation, a total of 1,814 high-quality SNP markers were developed: 1,121 of these were mapped to the female genetic map, 759 to the male map, and 1,646 to the integrated map. A comparison of the genetic maps to the published Vitis vinifera genome revealed both conservation and variations.
The applicability of next generation RAD sequencing for genotyping a grape F1 population was demonstrated, leading to the successful development of a genetic map with high density and quality using our designed SNP markers. Detailed analysis revealed that this newly developed genetic map can be used for a variety of genome investigations, such as QTL detection, sequence assembly and genome comparison.
KeywordsGrape Genetic map Next generation sequencing (NGS) Restriction-site associated DNA (RAD)
Grape (2n = 38) is one of the most important fruits worldwide, with a production of ~68 million tons over a harvested area of 7.2 million ha in 2010 (FAOSTAT, 2010). Grapes can be classified into either table or wine varieties, based on their intended mode of consumption, i.e., eaten raw or used to make wine. Consumption of grapes and wine has proven to be greatly beneficial for human health [1–4], and there has been a recent rise in the demand for high-quality grapes for human consumption. There is therefore a need to focus on grape improvement to optimize their attractive characteristics, such as contents of secondary metabolites, sugars and organic acids, resistance and yield. This can be achieved by using different germplasms from domesticated or wild-type grapes and then selecting for the genetic components that control the superior traits. However, it takes decades to produce advanced high-performing grape cultivars with the required traits, and there is still not enough resource which can produce grapes with high quality and quantity. High-density genetic map, one of the most valuable genomic resources, can largely reveal genome compositions and meet the requirement of high throughput superior traits selection among a lot of germplasms in most species, including plant and animal. Thus, construction of a high-quality genetic map for grape is necessary for its further studies and production.
In the past two decades, there have been a number of reports on the construction of grape genetic maps. Lodhi et al.  developed a genetic map for Vitis with 422 random amplified polymorphic DNA (RAPD) and 16 restriction fragment length polymorphism (RFLP) molecular markers, as well as a number of isozyme markers , possibly the first report of a complete genetic map for grape. From that study, a number of new genetic maps were developed, several of them based on the framework of that map. The latter studies generally made use of an F1 population as the plant material, with amplified fragment length polymorphisms (AFLP), simple sequence repeats (SSR), and single nucleotide polymorphisms (SNP) being the three major molecular marker types for map construction [6–16]. Although some genetic maps for grapes already exist, the total marker number on the linkage groups (LGs) of these existing maps is generally < 1,000 and some of these mapped markers have no sequence information. Thus a high-density genetic map for grape is still lacking, and one that covers a large number of molecular markers with sufficient sequence information is needed to meet the demand for improvement.
A key step in genetic map construction is the development of a set of testable molecular markers. In the last decade, a number of molecular marker technologies have been developed, including RAPD, AFLP, SSR and SNP. RAPD and AFLP have proven to be unstable due to many uncontrollable experimental conditions . SSRs are considered to be one of the most stable and reliable markers for genetic map construction, but the experiments are time- and cost-consuming . Thus, these markers are not suitable for high-density genetic map construction with high throughput. SNPs are single nucleotide polymorphisms or small InDels in the genome. They can be more numerous than other types of markers, but this is difficult to test. Before next generation sequencing (NGS) technique was developed, a number of other platforms were available for their identification, such as SNP Gene-Chip , high-resolution melt (HRM) analysis , TILLING and EcoTILLING [21, 22]. With the improved sequencing technology, the last two years have seen the development of NGS combining restriction-site associated DNA (RAD) for SNP testing . Pfender et al.  successfully used RAD markers to construct a high-density genetic map, which was subsequently employed to detect the QTL for resistance to stem rust in Lolium perenne. Using 2,383 RAD prior markers, an ultra-high-density genetic map was also developed for barley by Chutimanitsakun et al. , who showed that next generation RAD sequencing is a powerful high-throughput technique. Next generation RAD sequencing has also been successfully applied in other plants, including globe artichoke  and eggplant .
In this study, an F1 population of grape was constructed by crossing two interspecies hybridization progeny, Z180 and Beihong. Analyses of resveratrol content in the fruit skin, sugar and acid contents in the berry, berry size and cold resistance over several years revealed stable segregation of these traits in this F1 population. To take advantage of this F1 population, a high-density genetic map was constructed using next generation RAD sequencing for genotyping. The > 1,500 SNP markers contained in this map were analyzed, and aligned with the reference grapevine genome. Consequently, additional information on the genomic structures of different Vitis species was obtained, and the map can also be used to identify marker-linked loci that potentially control the superior traits of the two parents.
Materials and methods
Mapping population and DNA extraction
The F1 mapping population consisted of 100 progeny from a cross of Z180 (V. monticola × V. riparia) and Beihong (V. vinifera × V. amurensis) in 2003. Since pollen abortion occurred in Z180, Beihong was employed as the male parent. The seedlings of the two parents and their progeny were planted in the vineyard of the Germplasm Repository at the Institute of Botany of the Chinese Academy of Sciences in Beijing.
Young leaf samples (second and third leaves from the apex) were harvested from each individual F1 plant and the two parents at the beginning of the vegetative period (late spring). The samples were immediately stored in liquid nitrogen and transferred to a −70°C freezer. Young leaves (0.5 g) from each plant were ground in liquid nitrogen and their DNA extracted using the DNeasy plant mini prep kit (Qiagen). DNA concentration was measured and adjusted to the same level.
In-silico analysis of restriction enzyme-recognition sites on the reference grape genome
The sequence of the Vitis vinifera Pinot noir PN40024 12x genome assembly was downloaded from the international Grape Genome Browser (http://www.genoscope.cns.fr/externe/GenomeBrowser). Recognition sequences of 30 common restriction enzymes (data not shown) were chosen to investigate their digestion sites in the reference genome using Perl script. Total number of digestion sites, length of the resultant fragments, and their distribution were calculated from the results of the in-silico analysis.
Sample preparation and data analysis
Sample preparation for sequencing followed that in a number of published papers for NGS combined with RAD [23–25, 28], with a few modifications. Illumina Solexa adapters (2006 Illumina, Inc., all right reserved.), largely unmodified, were used for library construction. In brief, 2 μg genomic DNA from each sample (100 F1 progeny and both parents) was treated with 20 units (U) MseI (New England Biolabs [NEB]) for 60 min at 37°C in a 50 μl reaction. A quick blunting kit (NEB) was used to convert 30 μl of the digested sample to 5`-phosphorylated, blunt-ended DNA in a 50-μl reaction mixture; the reaction was performed with 30 μl of digested sample, 5 μl 10X blunting buffer, 5 μl 1 mM dNTP mix, 2 μl blunting enzyme mix and 8 μl sterile dH2O at room temperature for 30 min. A 3`-adenine overhang was added to the resulting samples in a 50-μl reaction with 32 μl blunt-ended DNA sample, 5 μl Klenow buffer (10X), 10 μl dATP (1 mM), 3 μl Klenow fragments (3` → 5` exo-, 5U/μl) and sterile dH2O to the final volume at 37°C for 1 h. Then 2 μl of 100 nM P1 and P2 adapter with a 3- to 5-bp plant-specific index (barcode) at the 5` end and a thymine overhang at the 3` end was added to each sample in a 50-μl reaction. The sequence of P1 and P2 adaptors: P1F: 5`-ACACTCTTTCCCTACACGACGCTCTTCCGATCTxxxT-3`(xxx indicated barcode); P1R: 5`phos-yyyAGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT-3` (yyy, reverse complement of xxx); P2F: 5`phos-AGATCGGAAGAGCGGTTCAGCAGGAATGCCGAG-3`; P2R: 5`-CTCGGCATTCCTGCTGAACCGCTCTTCCGATCTT-3`. A ligation reaction was carried out overnight at 16°C with T4 DNA ligase and 16 samples with different plant indices were pooled into one. DNA fragments from 400 to 500 bp (including the ~120-bp adaptor) were separated on a 1.5% agarose gel and purified using a MiniElute gel extraction kit (Qiagen). Finally, all pooled samples were amplified with Phusion High-Fidelity PCR Master Mix (NEB) for 18 cycles in a 100-μl reaction including 20 μl Phusion master mix, 5 μl of 10 μM modified Solexa amplification primer mix (AP1 and AP1; 2006 Illumina, Inc. , all right reserved) and sterile dH2O to the final volume. The AP1 and AP2 primers contained Illumina sequencing primer sites. The sequences are: AP1: 5`-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT-3`; P2: 5`-CAAGCAGAAGACGGCATACGAGATCGGTCTCGGCATTCCTGCTGAACCGCTCTTCCGATCT-3`; the underlined sequences are identical to Illumina sequencing primer sites. PCR products were repurified using the QIAquick PCR purification kit (Qiagen) and sequenced on a genome analyzer II instrument. All of these experiments were performed at Beijing's Biomarker Technologies Co. Ltd. (http://www.biomarker.com.cn/english/).
SNP identification and F1 plant genotyping were performed according to the method of Pfender et al. , with a few modifications. A number of Perl scripts (Biomarker Technologies Co. Ltd.) were programmed to conduct the analysis. In brief, low-quality data were discarded (five bases with Q score < 20) first, and Solexa sequences were assigned to the 102 plants according to their given index. The first 30 bp of each read (designated as RAD tags) were employed for subsequent analysis. For SNP marker identification, a cluster analysis was performed for both parents' data together. RAD tags were compared and nearly identical tags, with one or two mismatches (SNPs or 1- to 2-bp InDels), were assigned to one cluster. Clusters with > 200 or < 5 reads were discarded. More than one mismatch on the 30-bp sequence of the same RAD tag was considered a haplotype and regarded as one potential SNP marker in the subsequent analyses. In one RAD tag cluster, mismatches among different plants in the F1 population were considered putative polymorphisms and the different mismatches were regarded as multiple alleles. The parental genotypes for each RAD tag cluster were also analyzed according to the origin of the 30-bp tag sequence. To genotype all 100 F1 plants, their 30-bp sequences were also clustered and analyzed separately following the strategy applied in the cluster analysis of the two parent RAD tags. The genotypes for each RAD tag cluster of a single F1 plant were then determined by the identity between them and the corresponding clusters in the two parents.
Linkage map construction
Because of the lack of an anchor marker in this study, we first identified a set of SNP markers to assign the 19 grapevine chromosomes to19 LGs. This was performed in two steps: 1) we marked the segregation patterns of all 1,814 SNP markers as ab × cd, ef × eg, hk × hk, lm × ll and nn × np. Three types of markers, ab × cd, ef × eg and hk × hk, which could be mapped to both parental linkage maps, were regarded as candidate anchor markers; 2) the two representative 30-bp sequences (because all alleles of a SNP marker had two nearly identical 30-bp sequences, we could take the sequence of any allele representing the genotype of this SNP marker) of these candidate anchor markers were aligned with the sequence of the 12x genomic assembly of V. vinifera Pinot noir PN40024 using local BLAST software. The positions of each sequence for one SNP marker on the genome were identified by their highest number of hits. Three strict criteria were used to screen the candidate anchor marker: 1) it had to show a normal segregation ratio among the 100 F1 progeny; 2) both 30-bp end sequences had to align with the same chromosome position in the reference PN40024 genome; 3) the distance between the positions for the two end sequences on the reference genome had to fall between 200 and 500 bp (the expected size of the digested fragments was ~300–400 bp). The strategy for alignment of RAD tags with the reference genome was also used for the 1,646 SNP markers with the 19 chromosomes for subsequent comparison.
The double pseudo-test cross strategy of Grattapaglia and Sederoff  was applied, using JoinMap® 4.0 software, during the map construction. After data had been imported, a “CP” model was used for data mining. The ratio of marker segregation was calculated by Chi-square test. Markers showing significantly distorted segregation (P-value < 0.001) were excluded from the map construction. The genotypes of the 1,814 SNP markers were analyzed for linkage and recombination by applying the Kosambi function to estimate genetic map distances. To group all 1,814 markers, logarithm of odds (LOD) score thresholds ≥ 7 were used. After the LGs had been computed, their number was assigned according to the anchor markers mapped on them. The integrated map for both male and female plants was computed using the ‘Combine Group for Map Integration’ function.
Results and discussion
Selection of suitable restriction enzymes for RAD sequencing library construction
In this study, we did not sequence the whole genome of all F1 plants; rather, we sequenced the two ends of the ~300- to 400-bp RAD tags to simplify the grape genome and increase sequencing efficiency. Thus, selection of a suitable restriction enzyme for DNA digestion was key. Theoretically, two characteristics are required for an appropriate restriction enzyme: 1) because the NGS technology can only cover 75 to 100 bp of DNA at each end concurrently, the enzyme must be able to digest the genome of interest to a suitable size (e.g. ~300–400 bp); 2) the number of digested fragments of the expected size should be sufficient for subsequent manipulation (100,000–150,000 RAD tags). The V. vinifera Pinot noir PN40024 genome sequence was taken as the reference to search for an appropriate restriction enzyme.
SNP markers and their characteristics
SNP modulation for the F1 population
No. of reads
As already noted, the main advantages of NGS technology are low cost and high throughput. However, it also has a very serious disadvantage in its high probability of sequence error . To overcome this problem, high coverage of a specific sequence must be obtained. We digested the DNA and only then sequenced the RAD tags, greatly reducing the size of the genome. Jaillon et al. (2007) claimed that grapevine harbors a sequence that is ~470 Mb. During the genotyping of our 102 plants, we only manipulated ~80,709 valid clusters and each contained an ~70-bp sequence. Thus the grape genome was simplified to ~5.65 Mb (80,709 × 70 bp). This amounts to an ~83-fold reduction compared with the original 470 Mb reference genome, resulting in the requirement of very little data to achieve high coverage. According to our data, the average coverage for each tag was 17.0-fold in an individual plant. Moreover, because all sequence tags were from the two parents, Beihong and Z180, the number of alleles for each locus was ≤ 4. The total coverage for each tag at the population level was ~469-fold, leading us to adjust the SNPs in some loci where their coverage in an individual plant was insufficient. In addition, with these and subsequent strict criterions, we found the coverage of clusters corresponding to final SNP markers on the genetic map showed almost larger than 7 in an individual plant; only 24 showed from 5- to 7-fold coverage. Based on the above analyses, we concluded that the applied strategy provides high-throughput and high-quality identification of SNPs.
There were a number of possible patterns for the polymorphic markers in an F1 population (ab × cd, ef × eg, hk × hk, lm × ll, nn × np and aa × bb). However, the last pattern, aa × bb, could not be applied to the genetic map construction due to its lack of segregation in our F1 population, even though it probably constituted the largest proportion of all marker types. Thus, calculation of the segregating patterns for all loci would be necessary before a linkage map could be constructed. In addition, despite a high average coverage for the predicted RAD tag clusters, there were still a number of RAD tag clusters with low coverage in some F1 plants. To increase the accuracy of our data, only the clusters showing three or more fold coverage of > 80% of the F1 plants were used for subsequent development of SNP markers. We screened all 21,599 polymorphic clusters based on the above criteria and obtained 1,814 valid SNP markers with segregating patterns of ab × cd, ef × eg, hk × hk, lm × ll or nn × np (note that if two polymorphic clusters came from the same MseI-digested fragment, they were regarded as one marker). In addition to the coverage of the sequence data, the integrity for each locus among these 100 F1 individuals and their two parents was a key parameter in controlling map quality. We therefore investigated the data on missing rate for these plants, and found full integrity for the two parents, Z180 and Beihong, and 92.3% integrity on average for the 100 F1 plants. For a single SNP marker, the lowest integrity was ~85%, meeting the requirement for LG construction. Of these 1,814 SNP markers, 1,545 were homozygous for one parent and heterozygous for the other (960 for lm × ll and 585 for nn × np), constituting 85.2% of all selected SNP markers. However, the other three types of markers that could be mapped on both female and male linkage maps only amounted to 14.8% (ab × cd: 77, ef × eg: 171 and hk × hk: 21). This indicated that at most, 269 SNP markers could be used as shared markers for the integration of the two parents’ maps into one.
Because all of the SNP markers in this study were uniquely developed and no LG information was available, we identified a set of anchor markers that would indicate their chromosomal location. As described in Materials and Methods, the chromosome location of the 269 markers with ab × cd, ef × eg and hk × hk segregation patterns were detected according to their sequence alignment to the grape reference genome. After a series of strict selections and calculations, 212 markers clearly showed their chromosome location (Additional file 1: Table S1). Of these anchor markers, two were located on random chromosomes because the grape genomic sequence has not been completely assembled. The lowest number of anchor markers was on chromosome 15, with only two being usable for map construction (Additional file 1: Table S1). The average number of anchor markers for each chromosome was ~11.2 and only one chromosome had < 5 markers. This indicated that these anchor markers were sufficient for LG assignment.
Genetic map for 19 linkage groups (LGs)
Number of markers
Genetic sizes (cM)
Chromosome size (Mb)a
A number of future studies can be based on the high-density genetic map developed in this work. First, several excellent traits exist in one of the two parents. Thus, a given trait might be improved by selection of markers which are linked to elite loci or alleles after QTL detection. Moreover, several excellent traits might be combined in one grape plant, thereby producing a new cultivar, through a series of crosses and marker-assisted selection (MAS). Second, compared to other genetic maps for grape, there are two obvious advantages: high density and complete sequence information for all markers (Additional file 3: Table S2). These advantages could greatly benefit comparative mapping and genome assembly. The markers' combined 60-bp sequences mapped to the LGs could be used as anchors for the genome. Although the genome sequence of grapevine was published several years ago, it still has a number of gaps and random sequences [30, 31]. In this study, a set of markers could be aligned to the random chromosomes of V. vinifera Pinot noir PN40024 (data not shown). According to their positions on LGs, it might be easy to put the random chromosomes into the common one. On the other hand, the published grape genome is only for V. vinifera, and the genome structures of different Vitis species are expected to be more or less different due to the long evolutionary history of the Vitaceae . Thus, comparing the genome characteristics of the different species could give us a better understanding of grape. The 1,646 mapped markers' combined 60-bp sequences could be used as shared anchors to compare genetic and physical maps (Additional file 3: Table S2). These studies might facilitate use of the grape genomic resource.
Comparison of genetic and physical maps
To compare the genetic and physical maps, we investigated the locations of all 1,814 SNP markers on the reference genome. The high-quality 30-bp sequences from both ends of each SNP marker were employed for the location search by aligning them to the reference genome. A total of 1,456 SNP markers showed a match between their two ends and the same positions (intervals of 200–500 bp) on the reference genome; 106 markers only showed a match for one end to one position on the reference genome, while the other end had no match; the remaining 252 markers showed no match to the reference genome, showed a conflict in matching positions for the two ends, or were mapped on the random genome. To increase accuracy, only the first type of markers (1,456 SNP markers) was used to compare the genetic and physical maps.
Number of common markers between genetic and physical maps for 19 individual chromosomes
Number of markers
Moreover, using the high-quality, high-density genetic map, we could investigate not only the variation in genome structure among different species but also the variation that occurred during interspecies crosses. There have been a number of studies on interspecies hybridization in the past decade. Chromosome rearrangement, retrotransposon activation and SSR mutations have been seen in interspecies crosses between different types of Brassica species . For our plant material, the parents came from two separate interspecies crosses (V. monticola × V. riparia and V. vinifera × V. amurensis), and thus the population contained four grape pedigrees. An overview comparative analysis of the genetic map and the reference genome (Figure 7) reveals a number of markers in some regions that were not in the same order. Moreover, we achieved a similar result by comparing the Z180, Beihong and integrated genetic maps (Additional file 2: Figure S1). The variations among the different species might be the first reason for this non-uniformity; however, genomic variation occurring due to Vitis interspecies crosses might also exist because we observed variations in the positions of a number of markers as a uniform block between the male genetic map and the V. vinifera physical map; the male parent (Beihong) harboured half of the V. vinifera pedigree. Thus, with our detailed and complete investigation of the genetic map, more knowledge of the variation among different species and interspecies crosses can be obtained in the future.
We constructed a genetic map of a Z180 × Beihong F1 population of high density and quality. According to the analysis of the SNPs and their sequence information, we conclude that next generation RAD sequencing is a powerful strategy for genotyping. With further characterization of the genetic map, variations and conservation between the genetic map and reference genome were clearly detected. This genetic map is expected to be useful for QTL detection, sequence assembly and genome structure comparisons.
Financial support for this work was provided by the National Natural Science Foundation of China (NSFC accession No.: 31171931 and 31130047).
- Wu CD: Grape products and oral health. J Nutr. 2009, 139 (9): 1818S-1823S. 10.3945/jn.109.107854.PubMedPubMed CentralView ArticleGoogle Scholar
- Pezzuto JM, Venkatasubramanian V, Hamad M, Morris KR: Unraveling the relationship between grapes and health. J Nutr. 2009, 139 (9): 1783S-1787S. 10.3945/jn.109.107458.PubMedPubMed CentralView ArticleGoogle Scholar
- Dohadwala MM, Vita JA: Grapes and cardiovascular disease. J Nutr. 2009, 139 (9): 1788S-1793S. 10.3945/jn.109.107474.PubMedPubMed CentralView ArticleGoogle Scholar
- Bertelli AA, Das DK: Grapes, wines, resveratrol, and heart health. J Cardiovasc Pharmacol. 2009, 54 (6): 468-476. 10.1097/FJC.0b013e3181bfaff3.PubMedView ArticleGoogle Scholar
- Lodhi MA, Daly MJ, Ye GN, Weeden NF, Reisch BI: A molecular marker based linkage map of Vitis. Genome. 1995, 38 (4): 786-794. 10.1139/g95-100.PubMedView ArticleGoogle Scholar
- Grando MS, Bellin D, Edwards KJ, Pozzi C, Stefanini M, Velasco R: Molecular linkage maps of Vitis vinifera L. and Vitis riparia Mchx. Theor Appl Genet. 2003, 106 (7): 1213-1224.PubMedGoogle Scholar
- Adam-Blondon AF, Roux C, Claux D, Butterlin G, Merdinoglu D, This P: Mapping 245 SSR markers on the Vitis vinifera genome: a tool for grape genetics. Theor Appl Genet. 2004, 109 (5): 1017-1027. 10.1007/s00122-004-1704-y.PubMedView ArticleGoogle Scholar
- Doucleff M, Jin Y, Gao F, Riaz S, Krivanek AF, Walker MA: A genetic linkage map of grape, utilizing Vitis rupestris and Vitis arizonica. Theor Appl Genet. 2004, 109 (6): 1178-1187. 10.1007/s00122-004-1728-3.PubMedView ArticleGoogle Scholar
- Riaz S, Dangl GS, Edwards KJ, Meredith CP: A microsatellite marker based framework linkage map of Vitis vinifera L. Theor Appl Genet. 2004, 108 (5): 864-872. 10.1007/s00122-003-1488-5.PubMedView ArticleGoogle Scholar
- Doligez A, Adam-Blondon AF, Cipriani G, Di Gaspero G, Laucou V, Merdinoglu D, Meredith CP, Riaz S, Roux C, This P: An integrated SSR map of grapevine based on five mapping populations. Theor Appl Genet. 2006, 113 (3): 369-382. 10.1007/s00122-006-0295-1.PubMedView ArticleGoogle Scholar
- Lamoureux D, Bernole A, Le Clainche I, Tual S, Thareau V, Paillard S, Legeai F, Dossat C, Wincker P, Oswald M, et al: Anchoring of a large set of markers onto a BAC library for the development of a draft physical map of the grapevine genome. Theor Appl Genet. 2006, 113 (2): 344-356. 10.1007/s00122-006-0301-7.PubMedView ArticleGoogle Scholar
- Lowe KM, Walker MA: Genetic linkage map of the interspecific grape rootstock cross Ramsey (Vitis champinii) x Riparia Gloire (Vitis riparia). Theor Appl Genet. 2006, 112 (8): 1582-1592. 10.1007/s00122-006-0264-8.PubMedView ArticleGoogle Scholar
- Di Gaspero G, Cipriani G, Adam-Blondon AF, Testolin R: Linkage maps of grapevine displaying the chromosomal locations of 420 microsatellite markers and 82 markers for R-gene candidates. Theor Appl Genet. 2007, 114 (7): 1249-1263. 10.1007/s00122-007-0516-2.PubMedView ArticleGoogle Scholar
- Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, et al: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring Pinot Noir bacterial artificial chromosome contigs. Genetics. 2007, 176 (4): 2637-2650. 10.1534/genetics.106.067462.PubMedPubMed CentralView ArticleGoogle Scholar
- Marguerit E, Boury C, Manicki A, Donnart M, Butterlin G, Nemorin A, Wiedemann-Merdinoglu S, Merdinoglu D, Ollat N, Decroocq S: Genetic dissection of sex determinism, inflorescence morphology and downy mildew resistance in grapevine. Theor Appl Genet. 2009, 118 (7): 1261-1278. 10.1007/s00122-009-0979-4.PubMedView ArticleGoogle Scholar
- Zhang J, Hausmann L, Eibach R, Welter LJ, Topfer R, Zyprian EM: A framework map from grapevine V3125 (Vitis vinifera 'Schiava grossa' x 'Riesling') x rootstock cultivar 'Borner' (Vitis riparia x Vitis cinerea) to localize genetic determinants of phylloxera root resistance. Theor Appl Genet. 2009, 119 (6): 1039-1051. 10.1007/s00122-009-1107-1.PubMedView ArticleGoogle Scholar
- Gerber S, Mariette S, Streiff R, Bodenes C, Kremer A: Comparison of microsatellites and amplified fragment length polymorphism markers for parentage analysis. Mol Ecol. 2000, 9 (8): 1037-1048. 10.1046/j.1365-294x.2000.00961.x.PubMedView ArticleGoogle Scholar
- Woodhead M, Russell J, Squirrell J, Hollingsworth PM, Mackenzie K, Gibby M, Powell W: Comparative analysis of population genetic structure in Athyrium distentifolium (Pteridophyta) using AFLPs and SSRs from anonymous and transcribed gene regions. Mol Ecol. 2005, 14 (6): 1681-1695. 10.1111/j.1365-294X.2005.02543.x.PubMedView ArticleGoogle Scholar
- Divne AM, Allen M: A DNA microarray system for forensic SNP analysis. Forensic Sci Int. 2005, 154 (2–3): 111-121.PubMedView ArticleGoogle Scholar
- Vossen RH, Aten E, Roos A, den Dunnen JT: High-resolution melting analysis (HRMA): more than just sequence variant screening. Hum Mutat. 2009, 30 (6): 860-866. 10.1002/humu.21019.PubMedView ArticleGoogle Scholar
- Wang N, Wang Y, Tian F, King GJ, Zhang C, Long Y, Shi L, Meng J: A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING. New Phytol. 2008, 180 (4): 751-765. 10.1111/j.1469-8137.2008.02619.x.PubMedView ArticleGoogle Scholar
- Wang N, Shi L, Tian F, Ning H, Wu X, Long Y, Meng J: Assessment of FAE1 polymorphisms in three Brassica species using EcoTILLING and their association with differences in seed erucic acid contents. BMC Plant Biol. 2010, 10: 137-10.1186/1471-2229-10-137.PubMedPubMed CentralView ArticleGoogle Scholar
- Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, Selker EU, Cresko WA, Johnson EA: Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One. 2008, 3 (10): e3376-10.1371/journal.pone.0003376.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfender WF, Saha MC, Johnson EA, Slabaugh MB: Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theor Appl Genet. 2011, 122 (8): 1467-1480. 10.1007/s00122-011-1546-3.PubMedView ArticleGoogle Scholar
- Chutimanitsakun Y, Nipper RW, Cuesta-Marcos A, Cistue L, Corey A, Filichkina T, Johnson EA, Hayes PM: Construction and application for QTL analysis of a Restriction Site Associated DNA (RAD) linkage map in barley. BMC Genomics. 2011, 12: 4-10.1186/1471-2164-12-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Scaglione D, Acquadro A, Portis E, Tirone M, Knapp SJ, Lanteri S: RAD tag sequencing as a source of SNP markers in Cynara cardunculus L. BMC Genomics. 2012, 13: 3-10.1186/1471-2164-13-3.PubMedPubMed CentralView ArticleGoogle Scholar
- Barchi L, Lanteri S, Portis E, Acquadro A, Vale G, Toppino L, Rotino GL: Identification of SNP and SSR markers in eggplant using RAD tag sequencing. BMC Genomics. 2011, 12: 304-10.1186/1471-2164-12-304.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang X, Feng Q, Qian Q, Zhao Q, Wang L, Wang A, Guan J, Fan D, Weng Q, Huang T, et al: High-throughput genotyping by whole-genome resequencing. Genome Res. 2009, 19 (6): 1068-1076. 10.1101/gr.089516.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: mapping strategy and RAPD markers. Genetics. 1994, 137 (4): 1121-1137.PubMedPubMed CentralGoogle Scholar
- Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS One. 2007, 2 (12): e1326-10.1371/journal.pone.0001326.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): 463-467. 10.1038/nature06148.PubMedView ArticleGoogle Scholar
- Mardis ER: The impact of next-generation sequencing technology on genetics. Trends Genet. 2008, 24 (3): 133-141. 10.1016/j.tig.2007.12.007.PubMedView ArticleGoogle Scholar
- Yang S, Zhang X, Yue JX, Tian D, Chen JQ: Recent duplications dominate NBS-encoding gene expansion in two woody species. Mol Genet Genomics. 2008, 280 (3): 187-198. 10.1007/s00438-008-0355-0.PubMedView ArticleGoogle Scholar
- Zou J, Fu D, Gong H, Qian W, Xia W, Pires JC, Li R, Long Y, Mason AS, Yang TJ, et al: De novo genetic variation associated with retrotransposon activation, genomic rearrangements and trait variation in a recombinant inbred line population of Brassica napus derived from interspecific hybridization with Brassica rapa. Plant J. 2011, 68 (2): 212-224. 10.1111/j.1365-313X.2011.04679.x.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.