De novo sequencing of sunflower genome for SNP discovery using RAD (Restriction site Associated DNA) approach
© Pegadaraju et al.; licensee BioMed Central Ltd. 2013
Received: 19 February 2013
Accepted: 9 August 2013
Published: 15 August 2013
Application of Single Nucleotide Polymorphism (SNP) marker technology as a tool in sunflower breeding programs offers enormous potential to improve sunflower genetics, and facilitate faster release of sunflower hybrids to the market place. Through a National Sunflower Association (NSA) funded initiative, we report on the process of SNP discovery through reductive genome sequencing and local assembly of six diverse sunflower inbred lines that represent oil as well as confection types.
A combination of Restriction site Associated DNA Sequencing (RAD-Seq) protocols and Illumina paired-end sequencing chemistry generated high quality 89.4 M paired end reads from the six lines which represent 5.3 GB of the sequencing data. Raw reads from the sunflower line, RHA 464 were assembled de novo to serve as a framework reference genome. About 15.2 Mb of sunflower genome distributed over 42,267 contigs were obtained upon assembly of RHA 464 sequencing data, the contig lengths ranged from 200 to 950 bp with an N50 length of 393 bp. SNP calling was performed by aligning sequencing data from the six sunflower lines to the assembled reference RHA 464. On average, 1 SNP was located every 143 bp of the sunflower genome sequence. Based on several filtering criteria, a final set of 16,467 putative sequence variants with characteristics favorable for Illumina Infinium Genotyping Technology (IGT) were mined from the sequence data generated across six diverse sunflower lines.
Here we report the molecular and computational methodology involved in SNP development for a complex genome like sunflower lacking reference assembly, offering an attractive tool for molecular breeding purposes in sunflower.
Domestic sunflower (Helianthus annuus L., 2n = 2x = 34, haploid genome size ~3.5 Gbp) is native to North America  and widely cultivated as oilseed and confection crop types. Besides being economically important, sunflower also serves as model in ecological and evolutionary studies [2–4]. A major focus in both public and private sunflower breeding programs has been to develop sunflower hybrid varieties with improved yield, oil content and resistance to a wide range of diseases. Breeding new hybrids by conventional practices mostly is slow and uncertain; however, application of molecular markers can improve efficiency of plant selection, saving time and providing accuracy in a breeding program [5–8]. A wide range of molecular markers such as RFLP, AFLP, SSR and TRAP developed in sunflower have successfully enabled construction of high density genetic maps [9–15] and led to identification of molecular markers linked to disease resistance genes [16–19]. However, in general, practical usage of these markers for routine breeding purposes is limited due to high assay cost, low reproducibility, and lack of QTL validation studies .
In recent years, SNP markers have gained popularity in crop breeding programs due to their low cost, high throughput efficiency, and abundance. Particularly in association mapping studies, SNPs are the preferred marker type since they involves scanning whole genomes with extremely high marker densities to identify closely linked markers to causal polymorphisms [21, 22]. It is estimated that due to low linkage disequilibrium and high haplotype diversity, SNPs in the order of several thousand would be needed to successfully conduct genome wide association analysis in sunflower .
Large-scale discovery of genome-wide distributed SNPs can be effectively conducted with the aid of massively parallel, next-generation sequencing (NGS) technologies . Several studies that involve whole genome sequencing (WGS) efforts have led to the successful discovery of SNPs in Arabidopsis , humans , and Medicago . NGS technologies have also been extended for SNP discovery in large and complex genomes that lack an assembled reference genome [28, 29]. A common approach in these situations is the use of a complexity reduction strategy that is designed to selectively interrogate a small percentage of the target genome [30, 31]. By restricting sequencing on a smaller fraction of the genome, overall sequencing costs are reduced compared to WGS strategies, while still identifying a large amount of genetic variation. For instance, by using an RNA sequencing approach (RNA-Seq) on tissues from two diverse maize inbred lines, more than 4900 SNPs associated with 2,400 genes were identified and validated . Similarly, previous work in sunflower produced nearly 10,000 SNPs with RNA-Seq .
In addition, a cadre of methods have been developed that involve the usage of restriction enzymes on genomic DNA for complexity reduction. These strategies can be at the nucleotide level and are viewed as simple and highly efficient methods in plant and animal genome sequencing studies. One such method, CRoPS (complexity reduction of polymorphic sequences) can overcome the problems associated with highly duplicated regions in complex genomes that hamper the process of SNP identification . Restriction site-associated DNA sequencing (RAD-Seq) is an emerging method for SNP detection in genomes and is based on identifying polymorphic variants adjacent to restriction enzyme digestion sites [34, 35]. Application of RAD-Seq for identifying genetic variants has been demonstrated in a variety of species with and without an available reference genome [36, 37]. More recently paired-end RAD-Seq (RAD-PE) has been used in a variety of efforts for both genome assembly and SNP marker development [38, 39].
Furthermore, RAD-Seq approach has also been exploited in wide range of other studies such as association mapping , population genetics inferences [41–44], genetic mapping [35, 45, 46] and in estimation of allele frequencies . RAD-Seq differs from RNA-Seq in that non-transcribed loci are also sequenced, thus affording us an opportunity to broaden the known SNPs in sunflower to include those outside of transcribed regions. Here we demonstrate the use of paired-end RAD sequencing to enable efficient, cost-effective, high throughput marker development in H. annuus, a major oil crop without an assembled genome sequence. Results on the use of this sequence resource for detection of sequence variation and design of SNP marker panels for Illumina Infinium Genotyping Technology (IIGT) are also discussed.
Results and discussion
Paired-end RAD-Seq and de novo assembly
To promote SNP identification in low-copy, gene-rich regions of the 3.5 Gbp sunflower genome, a species expected to contain upwards of 80% retrotransposon content , the 5-methylcytosine (5mC) sensitive type II nuclease PstI (5’CTGCA/G’3) was selected for RAD-Seq in each of the six lines. Numerous studies have documented retroelement-dense regions of plant genome are often subjected to cytosine methylation of CpG, CpNpG and CpNpN nucleotides [49–51]. Restriction enzymes such as PstI, which do not cleave 5mC-modified DNA, have been shown to specifically sample the hypomethylated genomic fraction of plant genomes .
Paired-end RAD-Seq assembly statistics
Number of RHA 464 contigs assembled
Contigs removed due to plastid homology
Number of contigs retained
Total assembly length (bp)
Minimum contig length (bp)
Maximum contig length (bp)
N50 Contig Length (bp)
N90 Contig Length (bp)
Comparison of RAD-Seq assemblies from RHA 464 to preexisting sunflower unigenes at the Dana Farber Cancer Institute (DFCI) confirmed high sequence identity between RAD assemblies and known sunflower genomic sequences. A representative alignment is shown in Additional file 3 illustrating the match of a single paired-end RAD-Seq contig with tentative consensus EST TC57527 from the H. annuus DFCI EST database. The high sequence coverage inherent in paired-end RAD-Seq minimizes sequencing and assembly errors, as each nucleotide in the contig is derived from the consensus of many overlapping Illumina sequence reads. Experimental paired-end RAD-Seq studies using the sequenced B73 maize cultivar have placed the single base accuracy for paired-end RAD-Seq contigs at approximately 99.95% (unpublished data). This accuracy in RAD-Seq contig assembly is a key advantage for applications involving genome assembly and downstream marker development.
Helianthus SNP filtering and statistics
Helianthus samples sequenced:
Total number of SNP variants identified:
Total possible SNP genotypes in population:
SNP genotypes with high confidence call:
SNP genotypes with missing or low quality data:
SNP loci with < 50% genotype data:
SNPs passing initial filters:
Number of fixed genotype calls:
Number of heterozygous genotype calls:
SNP Loci with insufficient flanking sequence for IIGT*:
SNP Loci with nearby polymorphism (< 50 bp):
SNP Loci meeting all defined IIGT assay design criteria:
The application of next-generation DNA sequencing to generate large numbers of genetic markers has revolutionized plant breeding, facilitating both molecular genetic research and marker-assisted selection efforts. We have demonstrated paired-end RAD-Seq is an efficient and cost-effective means of SNP discovery in a species with a complex, highly repetitive genome. From less than a flowcell of Illumina paired-end sequence data we sequenced six diverse sunflower lines, assembled over 40,000 high-quality sequence contigs with an N50 contig length of 393 nucleotides, mined more than 100,000 sequence variations from the sunflower genome and identified 16,467 candidate SNPs suitable for downstream genotyping. The RAD-Seq method is appropriate for studies where many thousands of SNPs need to be rapidly identified at a low cost, in a format suitable for high-throughput genotyping.
Plant material and DNA extraction
Sunflower inbred lines (TX16R, CR29, SEEDS2000 B-Line, HA 467, RHA 468, and RHA 464) were grown under laboratory greenhouse conditions for four weeks, all true leaves were harvested and lyophilized prior to DNA extraction. DNA was extracted from 40 mg of each inbred line with the DNeasy 96 Plant Kit (Qiagen) using a modified protocol. Tissue was pulverized with 3 mm beads in a Harbil shaker. Buffer AP1 with DX and RNaseA was added to the tissue, 500 μL per sample, and incubated at 55°C for 60 min. Buffer AP2 was added, 150 μL per well, and incubated at −20°C for 15 min. AP3/E was combined with supernatant, 600 μL and 400 μL respectively, and then added to the binding plates. The rest of the extraction was carried out according to kit instructions. DNA was eluted in a final volume of 50 μL.
DNA was quantified using the PicoGreen kit (Molecular Probes) according to the kit instructions. A standard curve was made using quantified λ DNA from 100 to 0 ng/μL. A 1/200 dilution of Picogreen reagent in 1x TE (provided in kit) was mixed with 2 μL of isolated DNA, briefly vortexed, and incubated in the dark for 5 min. Assays were performed in black 96-well Fluotrac plates and fluorescence was measured with a Spectramax Gemini XPS (Molecular Devices) using 485 nm excitation and 538 nm emission.
RAD library preparation protocols
Genomic DNA from six selected sunflower inbred lines (TX16R, CR29, SEEDS2000 B-Line, HA 467, RHA 468, and RHA 464) was digested with the restriction endonuclease PstI and processed into RAD libraries similar to the method of Baird, et al., 2008. Briefly, ~300 ng of genomic DNA was digested for 60 min at 37°C in a 50 μL reaction with 20 units (U) of PstI (New England Biolabs [NEB]). After digestion, samples were heat-inactivated for 20 min at 65°C followed by addition of 2.0 μL of 100 nM P1 Adapter(s), a modified Solexa© adapter (Illumina, Inc.). PstI P1 adapters each contained a unique multiplex sequence index (barcode) which is read during the first four nucleotides of the Illumina sequence read. 100 nM P1adaptors were added to each sample along with 1 μL of 10 mM rATP (Promega), 1 μL 10× NEB Buffer 4, 1.0 μL (1000 U) T4 DNA Ligase (high concentration, Enzymatics, Inc), and 5 μL H2O which was then incubated at room temperature (RT) for 20 min. Samples were again heat-inactivated for 20 min at 65°C, pooled and randomly sheared with a Bioruptor (Diagenode) to an average size of 500 bp. Samples were then run out on a 1.5% agarose (Sigma), 0.5X TBE gel, and DNA 300 bp to 800 bp was isolated using a MinElute Gel Extraction Kit (Qiagen). End blunting enzymes (Enzymatics, Inc) were then used to polish the ends of the DNA. Samples were then purified using a MinElute column (Qiagen) and 15 U of Klenow exo− (Enzymatics) was used to add adenine (Fermentas) overhangs on the 3′ end of the DNA at 37°C. After subsequent purification, 1 μL of 10 μM P2 adapter, a divergent modified Solexa© adapter (Illumina, Inc.), was ligated to the obtained DNA fragments at 18°C. Samples were again purified and eluted in 50 μL. The eluate was quantified using a Qubit fluorimeter and 20 ng of this product was used in a PCR amplification with 20 μL Phusion Master Mix (NEB), 5 μL of 10 μM modified Solexa© Amplification primer mix (Illumina, Inc.) and up to 100 μL H2O. Phusion PCR settings followed product guidelines for a total of 18 cycles. Again, samples were gel purified, excising DNA from the 300 to 700 bp size range, and diluted to 1 nM.
A set of RAD libraries generated from lines TX16R, CR29, SEEDS2000 B-Line, HA 467, RHA 468, and RHA 464 was run on an Illumina Genome Analyzer IIx at the University of Oregon High Throughput Sequencing Facility in Eugene, Oregon. Illumina protocols were followed for an asymmetric length paired end sequencing run, with an initial 40 bp read and second 80 bp read.
Bioinformatics – Sequence processing, paired-end RAD-Seq assembly and SNP detection
A combination of open source and proprietary bioinformatics tools was used for processing and sequence analysis. A list of open source programs, versions, and commands used in sequence analysis can be found in a supplemental file (Additional file 4). Initially, raw sequence data produced on two GAIIx sequence lanes were sorted by the appropriate multiplex index (MID) or “barcode” assigned to each sunflower line during RAD-Seq library construction. During de-multiplexing, indexes were trimmed from reads and the remaining sequence segregated to individual sample files. Reads from RHA 464 were then processed to extract low quality sequences. Any sequence with an average phred-scaled quality score below 20 (Q20) over the last 5 base pairs of the read was discarded. Remaining reads were then collapsed into RAD sequence clusters sharing 100% sequence identity across the single end Illumina read. To maximize efficient assembly of sequences we imposed a minimum of 50x and maximum 750x sequence coverage at any RAD sequence cluster. These thresholds were selected for this effort, because single loci with coverage under 50x would be expected to suffer from low sequence coverage ((80 bp × 50)/ 400 bp = 10.0x)) resulting in short and fragmented contig assemblies, while loci with greater than 500 identical SE (Single-end) reads may be composed of high-copy contaminant DNA (plastids) or dosage from multiple genomic loci (e.g. retrotransposon derived sequences). The paired end sequences for each RAD locus were extracted from these selected loci and passed to the Velvet sequence assembler (version 1.0.18) for contig assembly . Contigs not reaching a minimum length of 200 bp were excluded from the assembly.
Sequence reads from TX16R, CR29, Seeds 2000 B-line, HA 467, and RHA 468, were aligned to the reference RHA 464 assembly using the short-read aligner Bowtie (version 0.12.5) . Alignment thresholds were specified which allowed up to three base pair mismatches between the 80 bp Illumina read and the reference (>95% identity). Reads not uniquely mapping (e.g. aligning to more than one contig in the RHA 464 reference) were discarded and not considered in the analysis. Bowtie alignments were piped to SAM tools (version 0.1.14) and reformatted into BAM and pileup files for SNP identification . Sequence variants from pileups were then condensed into a variant call format (VCF) file using custom perl scripts. To be considered for genotyping design, a SNP had to have a minimum sequencing coverage of 4x in at least three lines, with at least 50 bp of flanking genomic sequence surrounding the target SNP. Variants with nearby flanking polymorphisms (within 50 bp of the candidate marker) were also excluded from further consideration for Infinium genotyping design.
This project was supported by the National Sunflower Association. The authors are grateful to Chondra Carlson, Robert Bialozynski, Mark Blackstad and Ruhi Rai for their technical support in the project. Authors would like to thank Nuseed Americas Inc. & USDA for providing CR29, Seeds 2000-B, and TX16R sunflower lines for the sequencing project.
- Harter AV, Gardner KA, Falush D, Lentz DL, Bye RA, Rieseberg LH: Origin of extant domesticated sunflowers in eastern North America. Nature. Nature. 2004, 430: 201-205. 10.1038/nature02710.View ArticlePubMedGoogle Scholar
- Rieseberg LH, Raymond O, Rosenthal DM, Lai Z, Livingstone K, Nakazato T, Durphy JL, Schwarzbach AE, Donovan LA, Lexer C: Major ecological transitions in wild sunflowers facilitated by hybridization. Science. 2003, 301: 1211-1216. 10.1126/science.1086949.View ArticlePubMedGoogle Scholar
- Snow AA, Pilson D, Rieseberg LH, Paulsen M, Pleskac N, Reagon MR, Wolf DE: A Bt transgene reduces herbivory and enhances fecundity in wild sunflowers. Ecol Appl. 2003, 13: 279-286. 10.1890/1051-0761(2003)013[0279:ABTRHA]2.0.CO;2.View ArticleGoogle Scholar
- Lexer C, Lai Z, Rieseberg LH: Candidate gene polymorphisms associated with salt tolerance in wild sunflower hybrids: Implications for the origin of Helianthus paradoxus, a diploid hybrid species. New Phytol. 2004, 161: 225-233.PubMed CentralView ArticlePubMedGoogle Scholar
- Mackill DJ, Nguyen HT, Zhan J: Use of molecular markers in plant improvement programs for rain fed lowland rice. Field Crops Res. 1999, 64: 177-185. 10.1016/S0378-4290(99)00058-1.View ArticleGoogle Scholar
- McCouch SR, Doerge RW: QTL mapping in rice. Trends Genet. 1995, 11: 482-487. 10.1016/S0168-9525(00)89157-X.View ArticlePubMedGoogle Scholar
- Snowdon RJ, Friedt W: Molecular markers in Brassica oilseed breeding: Current status and future possibilities. Plant Breeding. 2004, 123: 1-8. 10.1111/j.1439-0523.2003.00968.x.View ArticleGoogle Scholar
- Collard BCY, Jahufer MZZ, Brouwer JB, Pang ECK: An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. Euphytica. 2005, 142: 169-196. 10.1007/s10681-005-1681-5.View ArticleGoogle Scholar
- Berry ST, Leon AJ, Hanfrey CC, Challis P, Burkholz A: Molecular marker analysis of Helianthus-annuus L. 2. Construction of an RFLP linkage map for cultivated sunflower. Theor Appl Genet. 1995, 91: 195-199.View ArticlePubMedGoogle Scholar
- Gentzbittel LF, Vear YX, Zhang A, Berville , Nicolas P: Development of a consensus linkage RFLP map of cultivated sunflower (Helianthus-annuus L). Theor App Genet. 1995, 90: 1079-1086.View ArticleGoogle Scholar
- Rieseberg LH, Choi HC, Chan R, Spore C: Genomic map of a diploid hybrid species. Heredity. 1993, 70: 285-293. 10.1038/hdy.1993.41.View ArticleGoogle Scholar
- Gedil MA, Wye C, Berry S, Segers B, Peleman J: An integrated restriction fragment length polymorphism - amplified fragment length polymorphism linkage map for cultivated sunflower. Genome. 2001, 44: 213-221. 10.1139/g00-111.View ArticlePubMedGoogle Scholar
- Heesacker A, Kishore VK, Gao W, Tang S, Kolkman JM: SSRs and INDELs mined from the sunflower EST database: Abundance, polymorphisms, and cross-taxa utility. Theor Appl Genet. 2008, 117: 1021-1029. 10.1007/s00122-008-0841-0.View ArticlePubMedGoogle Scholar
- Tang S, Yu JK, Slabaugh MB, Shintani DK, Knapp SJ: Simple sequence repeat map of the sunflower genome. Theor Appl Genet. 2002, 105: 1124-1136. 10.1007/s00122-002-0989-y.View ArticlePubMedGoogle Scholar
- Yu JK, Tang S, Slabaugh MB, Heesacker A, Cole G: Towards a saturated molecular genetic linkage map for cultivated sunflower. Crop Sci. 2003, 43: 367-387. 10.2135/cropsci2003.0367.View ArticleGoogle Scholar
- Tang S, Kishore VK, Knapp SJ: PCR-multiplexes for a genome-wide framework of simple sequence repeat marker loci in cultivated sunflower. Theor Appl Genet. 2003, 107: 6-19.PubMedGoogle Scholar
- Qi LL, Hulke BS, Vick BA, Gulya TJ: Molecular mapping of the rust resistance gene R4 to a large NBS-LRR cluster on linkage group 13 of sunflower. Theor Appl Genet. 2011, 123: 351-358. 10.1007/s00122-011-1588-6.View ArticlePubMedGoogle Scholar
- Lawson WR, Jan CC, Shatte T, Smith L, Kong GA, Kochman JK: DNA markers linked to the R2 rust resistance gene in sunflower (Helianthus annus L.) Facilitate anticipatory breeding for this disease variant. Mol Breeding. 2011, 28: 569-576. 10.1007/s11032-010-9506-1.View ArticleGoogle Scholar
- Bachlava E, Radwan OE, Abratti G, Tang S, Gao W, Heesacker AF, Bazzalo ME, Zambelli A, Leon AJ, Knapp SJ: Downy mildew (Pl8 and Pl14) and rust (RAdv) resistance genes reside in close proximity to tandemly duplicated clusters of non-TIR-like NBS-LRR-encoding genes on sunflower Chromosomes 1 and 13. Theor Appl Genet. 2011, 122: 1211-1221. 10.1007/s00122-010-1525-0.View ArticlePubMedGoogle Scholar
- Bernardo R: Molecular markers and selection for complex traits in plants: Learning from the Last 20 Years. Crop Sci. 2008, 48: 1649-1664. 10.2135/cropsci2008.03.0131.View ArticleGoogle Scholar
- Ganal MW, Altmann T, Röder MS: SNP identification in crop plants. Curr Opin Plant Biol. 2009, 12: 211-217. 10.1016/j.pbi.2008.12.009.View ArticlePubMedGoogle Scholar
- Rafalski A: Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol. 2002, 5: 94-100. 10.1016/S1369-5266(02)00240-6.View ArticlePubMedGoogle Scholar
- Zhu C, Gore M, Buckler ES, Yu J: Status and prospects of association mapping in plants. Plant Gen. 2008, 1: 5-20. 10.3835/plantgenome2008.02.0089.View ArticleGoogle Scholar
- Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML: Genomewide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011, 12: 499-510. 10.1038/nrg3012.View ArticlePubMedGoogle Scholar
- Weber AP, Weber KL, Carr K, Wilkerson C, Ohlrogge JB: Sampling the Arabidopsis transcriptome with massively parallel pyrosequencing. Plant Physiol. 2007, 144: 32-42. 10.1104/pp.107.096677.PubMed CentralView ArticlePubMedGoogle Scholar
- Wheeler DA, Srinivasan M, Egholm M, Shen Y, Chen L, McGuire A, He W, Chen YJ, Makhijani V, Roth GT, Gomes X, Tartaro K, Niazi F, Turcotte CL, Irzyk GP L, Chinault C, Song XZ, Liu Y, Yuan Y, Nazareth L, Qin X, Muzny DM, Margulies M, Weinstock GM, Gibbs RA, Rothberg JM: The complete genome of an individual by massively parallel DNA sequencing. Nature. 2008, 452: 872-876. 10.1038/nature06884.View ArticlePubMedGoogle Scholar
- Cheung F, Haas BJ, Goldberg SMD, May GD, Xiao Y, Town CD: Sequencing Medicago truncatula expressed sequenced tags using 454 Life Sciences technology. BMC Genomics. 2006, 7: 272-10.1186/1471-2164-7-272.PubMed CentralView ArticlePubMedGoogle Scholar
- You FM, Huo N, Deal KR, Gu YQ, Luo M-C, McGuire PE, Dvorak J, Anderson OD: Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence. BMC Genomics. 2011, 12: 59-10.1186/1471-2164-12-59.PubMed CentralView ArticlePubMedGoogle Scholar
- Bachlava E, Taylor CA, Tang S, Bowers JE M, Burke JM, Knapp SJ: SNP discovery and development of a high density genotyping array for sunflower. PLoS ONE. 2012, 7: e29814-10.1371/journal.pone.0029814.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu X, Ren C, Joshi T, Vuong T, Xu D, Nguyen HT: SNP discovery by high throughput sequencing in soybean. BMC Genomics. 2010, 11: 469-10.1186/1471-2164-11-469.PubMed CentralView ArticlePubMedGoogle Scholar
- Varshney RK, Nayak SN, May GD, Jackson SA: Next-generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol. 2009, 27: 522-530. 10.1016/j.tibtech.2009.05.006.View ArticlePubMedGoogle Scholar
- Barbazuk WB, Emrich SJ, Chen HD, Li L, Schnable PS: SNP discovery via 454 transcriptome sequencing. Plant J. 2007, 51: 910-918. 10.1111/j.1365-313X.2007.03193.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Orsouw NJV, Hogers R, Janssen A, Yalcin F, Snoeijers S, Verstege E, Schneiders H, Van Der PH, Van OJ, Verstegen H: Complexity reduction of polymorphic sequences (CRoPS): A novel approach for large-scale polymorphism discovery in complex genomes. PLoS ONE. 2007, 14: e1172-View ArticleGoogle Scholar
- Miller M, Dunham J, Amores A, Cresko W, Johnson E: Rapid and cost effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 2007, 17: 240-248. 10.1101/gr.5681207.PubMed CentralView ArticlePubMedGoogle 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: e3376-10.1371/journal.pone.0003376.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Chutimanitsakun Y, Nipper R, Cuesta-Marcos A, Cistue L, Corey A, Filichkina T, Johnson E, Hayes P: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Etter PD, Preston JL, Bassham S, Cresko WA, Johnson EA: Local de novo assembly of RAD paired-end contigs using short sequencing reads. PLoS ONE. 2011, 6: e18561-10.1371/journal.pone.0018561.PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson JC, Wang S, Wu Y, Li X, Antony G, White FF, Yu J: Single-nucleotide polymorphism discovery by high-throughput sequencing in sorghum. BMC Genomics. 2011, 12: 352-10.1186/1471-2164-12-352.PubMed CentralView ArticlePubMedGoogle Scholar
- Parchman TL, Gompert Z, Mudge J, Schilkey FD, Benkman CW, Buerkle CA: Genome-wide association genetics of an adaptive trait in lodgepole pine. Mol Ecol. 2012, 21: 2991-3005. 10.1111/j.1365-294X.2012.05513.x.View ArticlePubMedGoogle Scholar
- Emerson KJ, Merz CR, Catchen JM: Resolving postglacial phylogeography using high-throughput sequencing. Proc Natl Acad Sci USA. 2010, 107: 16196-16200. 10.1073/pnas.1006538107.PubMed CentralView ArticlePubMedGoogle Scholar
- Gompert Z, Forister ML, Fordyce JA, Nice CC, Williamson RJ, Buerkle CA: Bayesian analysis of molecular variance in pyrosequences quantifies population genetic structure across the genome of Lycaeides butterflies. Mol Ecol. 2010, 19: 2455-2473.View ArticlePubMedGoogle Scholar
- Hohenlohe PA, Bassham S, Etter PD, Stiffler N, Johnson EA, Cresko WA: Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags. PLoS Genetics. 2010, 6: e1000862-10.1371/journal.pgen.1000862.PubMed CentralView ArticlePubMedGoogle Scholar
- Arnold B, Corbett-Detig RB, Hartl D, Bomblies K: RADseq underestimates diversity and introduces genealogical biases due to nonrandom haplotype sampling. Mol Ecol. 2013, 22: 3179-3190. 10.1111/mec.12276.View ArticlePubMedGoogle Scholar
- Andolfatto P, Davison D, Erezyilmaz D: Multiplexed shotgun genotyping for rapid and efficient genetic mapping. Genome Res. 2011, 21: 610-617. 10.1101/gr.115402.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfender W, Saha M, Johnson E, Slabaugh M: Mapping with RAD (restriction site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. TAG Theor Appl Genet. 2011, 122: 1467-1480. 10.1007/s00122-011-1546-3.View ArticlePubMedGoogle Scholar
- Tassell CPV, Smith TPL, Matukumalli LK: SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat Methods. 2008, 5: 247-252. 10.1038/nmeth.1185.View ArticlePubMedGoogle Scholar
- Kane NC, Gill N, King MG, Bowers JE, Berges H, Gouzy J, Bachlava E, Langlade NB, Lai Z, Stewart M, Burke JM, Vincourt P, Knapp SJ, Rieseberg LH: Progress towards a reference genome for sunflower. Botany. 2011, 89: 429-437. 10.1139/b11-032.View ArticleGoogle Scholar
- Palmer LE, Rabinowicz PD, O’Shaughnessy AL, Balija VS, Nascimento LU, Dike S, De La BM, Martienssen RA, McCombie WR: Maize genome sequencing by methylation filtration. Science. 2003, 302: 2115-2117. 10.1126/science.1091265.View ArticlePubMedGoogle Scholar
- Rabinowicz PD, Citek R, Budiman MA, Nunberg A, Bedell JA, Lakey N, O’Shaughnessy AL, Nascimento LU, McCombie WR, Martienssen RA: Differential methylation of genes and repeats in land plants. Genome Res. 2005, 15: 1431-1440. 10.1101/gr.4100405.PubMed CentralView ArticlePubMedGoogle Scholar
- Ma XF, Gustafson JP: Timing and rate of genome variation in triticale following allopolyploidization. Genome. 2006, 49: 950-958. 10.1139/G06-078.View ArticlePubMedGoogle Scholar
- Reyna-Lopez GE, Simpson J, Ruiz-Herrera J: Differences in DNA methylation patterns are detectable during the dimorphic transition of fungi by amplification of restriction polymorphisms. Mol Gen Genet. 1997, 253: 703-710. 10.1007/s004380050374.View ArticlePubMedGoogle Scholar
- Hulke BS, Miller JF, Gulya TJ: Registration of the restorer oilseed sunflower germplasm RHA 464 possessing genes for resistance to downy mildew and sunflower rust. J Plant Reg. 2010, 4: 249-254. 10.3198/jpr2010.02.0114crg.View ArticleGoogle Scholar
- Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei W, Qi X, Wang L, Zhang Y, Hua W, Li D, Lv H, Zhang X: Characterization of the sesame (Sesamum indicum L.) Global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers. BMC Genomics. 2011, 12: 451-10.1186/1471-2164-12-451.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: 1000 Genome Project Data Processing Subgroup: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu A, Burke JM: Patterns of nucleotide diversity in wild and cultivated sunflower. Genetics. 2006, 173: 321-330. 10.1534/genetics.105.051110.PubMed CentralView ArticlePubMedGoogle Scholar
- Kolkman JM, Berry ST, Leon AJ, Slabaugh MB, Tang S, Gao W, Shintani DK, Burke JM, Knapp SJ: Single nucleotide polymorphisms and linkage disequilibrium in sunflower. Genetics. 2007, 177: 457-468. 10.1534/genetics.107.074054.PubMed CentralView ArticlePubMedGoogle 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.